
==== Front
Exp Mol Med
Exp Mol Med
Experimental & Molecular Medicine
1226-3613
2092-6413
Nature Publishing Group UK London

39085355
1284
10.1038/s12276-024-01284-1
Article
Low-level brain somatic mutations in exonic regions are collectively implicated in autism with germline mutations in autism risk genes
Kim Il Bin 12
Kim Myeong-Heui 2
Jung Saehoon 2
Kim Woo Kyeong 2
Lee Junehawk 3
http://orcid.org/0000-0002-5514-4189
Ju Young Seok 2
Webster Maree J. 4
Kim Sanghyeon 4
Kim Ja Hye 5
Kim Hyun Jung 6
http://orcid.org/0000-0003-3513-401X
Kim Junho junho.kim@skku.edu

7
http://orcid.org/0000-0001-5356-0827
Kim Sangwoo swkim@yuhs.ac

8
http://orcid.org/0000-0002-2299-630X
Lee Jeong Ho jhlee4246@kaist.ac.kr

29
1 grid.410886.3 0000 0004 0647 3511 Department of Psychiatry, CHA Gangnam Medical Center, CHA University School of Medicine, Seoul, 06135 Republic of Korea
2 https://ror.org/05apxxy63 grid.37172.30 0000 0001 2292 0500 Graduate School of Medical Science and Engineering, Korea Advanced Institute of Science and Technology, Daejeon, 34141 Republic of Korea
3 grid.249964.4 0000 0001 0523 5253 Center for Supercomputing Applications, Division of National Supercomputing, Korea Institute of Science and Technology Information, Daejeon, 34141 Republic of Korea
4 https://ror.org/01pj5nn22 grid.453353.7 0000 0004 0473 2858 Stanley Medical Research Institute, Laboratory of Brain Research, 9800 Medical Center Drive, Suite C-050, Rockville, MD 20850 USA
5 grid.267370.7 0000 0004 0533 4667 Department of Pediatrics, Asan Medical Center Children’s Hospital, University of Ulsan College of Medicine, Seoul, 05505 Republic of Korea
6 grid.222754.4 0000 0001 0840 2678 Department of Anatomy, Korea University College of Medicine, Seoul, 02841 Republic of Korea
7 https://ror.org/04q78tk20 grid.264381.a 0000 0001 2181 989X Department of Biological Sciences, Sungkyunkwan University, Suwon, 16419 Republic of Korea
8 https://ror.org/01wjejq96 grid.15444.30 0000 0004 0470 5454 Department of Biomedical Systems Informatics and Brain Korea 21 PLUS for Medical Science, Yonsei University College of Medicine, Seoul, 03722 Republic of Korea
9 SoVarGen, SoVarGen, Inc., Daejeon, 34141 Republic of Korea
1 8 2024
1 8 2024
8 2024
56 8 17501762
9 12 2023
15 4 2024
12 5 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Low-level somatic mutations in the human brain are implicated in various neurological disorders. The contribution of low-level brain somatic mutations to autism spectrum disorder (ASD), however, remains poorly understood. Here, we performed high-depth exome sequencing with an average read depth of 559.3x in 181 cortical, cerebellar, and peripheral tissue samples to identify brain somatic single nucleotide variants (SNVs) in 24 ASD subjects and 31 controls. We detected ~2.4 brain somatic SNVs per exome per single brain region, with a variant allele frequency (VAF) as low as 0.3%. The mutational profiles, including the number, signature, and type, were not significantly different between the ASD patients and controls. Intriguingly, when considering genes with low-level brain somatic SNVs and ASD risk genes with damaging germline SNVs together, the merged set of genes carrying either somatic or germline SNVs in ASD patients was significantly involved in ASD-associated pathophysiology, including dendrite spine morphogenesis (p = 0.025), mental retardation (p = 0.012), and intrauterine growth retardation (p = 0.012). Additionally, the merged gene set showed ASD-associated spatiotemporal expression in the early and mid-fetal cortex, striatum, and thalamus (all p < 0.05). Patients with damaging mutations in the merged gene set had a greater ASD risk than did controls (odds ratio = 3.92, p = 0.025, 95% confidence interval = 1.12–14.79). The findings of this study suggest that brain somatic SNVs and germline SNVs may collectively contribute to ASD-associated pathophysiology.

Brain somatic mutations in autism: insights from high-depth sequencing

Autism Spectrum Disorder is a complex condition influenced by various genetic factors, including inherited traits and new changes in genes. This study investigates the role of low-level brain somatic mutations in ASD. The researchers analyzed brain tissues from deceased individuals, both with and without ASD, using high-depth whole-exome sequencing. The results showed that low-level brain somatic mutations, along with inherited genetic variations, contribute to ASD’s genetic makeup. These mutations were found in genes linked to brain development and function. The study emphasizes the need to consider both inherited and somatic mutations to understand ASD’s genetic complexity. Researchers conclude that the interaction between somatic and inherited mutations is crucial in ASD, providing new insights into its genetic basis. This study enhances our understanding of ASD’s genetic diversity and suggests a multifaceted genetic contribution to the disorder. This summary was initially drafted using artificial intelligence, then revised and fact-checked by the author.

Subject terms

Genetics of the nervous system
Genome informatics
Autism spectrum disorders
https://doi.org/10.13039/501100003725 National Research Foundation of Korea (NRF) 2022R1C1C1010430 2021M3A9E4080784 2021M3E5D9021878 RS-2023-00278314 2019R1A3B2066619 Kim Il Bin Lee Junehawk Kim Junho Kim Sangwoo the Suh Kyungbae Foundationissue-copyright-statement© Korean Society for Biochemical and Molecular Biology 2024
==== Body
pmcIntroduction

ASD is characterized by phenotypic diversity ranging from defects in social interaction and communication to restricted and repetitive behaviors, interests, or activities. ASD is a genetically heterogeneous disorder1–3, with contributions from de novo as well as inherited mutations. Although the heritability of ASD is estimated to be up to 50–90%4–6, the genetic accountability for ASD is at most 30%7–9. De novo and inherited germline mutations are considered to contribute to ASD7,10–12. However, the genetic etiology causing a substantial proportion of sporadic ASD cases remains unclear.

Somatic mutations are postzygotic variations in somatic cells present in a single tissue or several tissues in an organism13. Somatic mutations in brain tissues have increasingly been found to be genetic causes of neurological disorders of previously unknown etiology13,14. Recent studies have reported that low-level brain somatic mutations, with a less than 5% variant allele frequency (VAF) that arise from neural stem cell niches during brain development or with aging, contribute to various neurological disorders, such as focal epilepsy15–17, schizophrenia18, and Alzheimer’s disease19. Regarding ASD, several studies using brain tissues have sequenced a limited set of known ASD genes and reported the presence of brain somatic mutations in a few ASD-related genes20,21. Additionally, recent studies using whole-genome sequencing, of which the read depth ranged from 200x to 360x and was much lower than that of high-depth whole-exome sequencing (>500x), showed that high-level somatic mutations with >5% VAF in neural enhancer sequences or large mosaic copy number variations are associated with ASD22,23. These studies support the hypothesis that brain somatic mutations can be implicated in the genetic architecture underlying ASD. However, there is a lack of studies examining whether low-level brain somatic mutations (e.g., less than 5% in VAF) in protein-coding genes are implicated in ASD at the genome-wide level.

Here, we performed high-depth whole-exome sequencing (WES) with an average throughput read-depth of 559.3x to accurately detect brain somatic single nucleotide variants (SNVs) in 181 postmortem tissues. The postmortem tissues comprised 86 specimens of multiple brain regions and 13 specimens of peripheral tissue from 24 ASD subjects, as well as 51 brain and 31 peripheral tissue specimens from 31 neurotypical controls. We performed independent, ultradeep targeted amplicon sequencing (TASeq) to accurately detect low-level brain somatic SNVs. Considering both somatic and germline mutations, we found that low-level brain somatic SNVs are collectively implicated in the pathogenesis of ASD with germline SNVs in ASD risk genes, suggesting that the heterogeneous genetic architecture of ASD may be further explained jointly by somatic as well as germline mutations.

Materials and Methods

Tissue collection

Fresh frozen human brain tissue (cortex and cerebellum) and paired peripheral tissue (heart and liver) were generously provided by the National Institute of Child Health and Human Development (NICHD). Samples were obtained from 24 clinically diagnosed autism spectrum disorder (ASD)-affected individuals and 5 age-matched unaffected controls. To increase the number of controls, we additionally acquired DNA from the human dorsolateral prefrontal cortex and peripheral tissues (spleen) of 26 neurotypical controls, which were generously provided by the Stanley Medical Research Institute (SMRI). The human brain and peripheral tissues were acquired from the NICHD and the SMRI, both of which confirmed official consent from all subjects. The research performed on deidentified postmortem human tissues was approved by the Institutional Review Board of the Korea Advanced Institute of Science and Technology.

DNA extraction and high-depth whole-exome sequencing

Genomic DNA (gDNA) was extracted from the brain and paired peripheral tissues using a QIAamp DNA Mini Kit (Qiagen). We prepared exome libraries following the manufacturer’s protocol (Agilent, Human All Exon V4/V5 + UTR 50 Mb Kit) using up to 1 μg of gDNA as an input. Then, we performed paired-end sequencing on an Illumina HiSeq 2000/2500 instrument (average throughput depth of 559.3×) according to the manufacturer’s instructions using exome libraries that passed quality control (QC-passed). We followed the GATK Best Practices (v3.5) workflow to generate analysis-ready bam files from QC-passed Fastq files. The Fastq files were aligned to a reference genome (GRCh38) using BWA-MEM to generate bam files, and PCR duplicates were marked by Picard. Reads adjacent to indels in the bam fi2000/2500 instrument (RealignerTargetCreator and IndelRealigner from the GATK analysis tools. Finally, we performed recalibration of the base quality scores with BaseRecalibrator from the GATK analysis tools for subsequent variant calling. We publicly deposited our raw fastq files generated from tissues from the National Institute of Child Health and Human Development for 24 ASD subjects and 5 controls (SUB6919131) and also the additional 26 controls from the Stanley Medical Research Institute (SUB6847120).

Somatic SNV calling

Somatic SNVs were independently detected using MuTect224 (v.3.6.0) and RePlow25 for single region-based and multiregion-based candidates, respectively. Both methods were applied to analyze 99 brain and 44 peripheral specimens from 44 individuals (paired cases) and 38 brain-only specimens from 11 individuals (unpaired cases). We ran MuTect2 with default options, including fraction_contamination (default = 0.02). We then checked for cross-contamination between samples using an in-house contamination filter and an independent tool, ContEst26. From the MuTect2 output from the paired cases, we excluded unreliable calls by applying the following criteria: (i) variant allele frequency (VAF) ≥ 20%, (ii) EBscore ≤5, and (iii) variants with all supporting reads located at either end of the reads. We also excluded unreliable calls via manual inspection with IGViewer (v2.3.94), checking for the following: (a) supporting reads with altered alleles that had no other base changes, if they were not heterozygotic/homozygotic SNPs; (b) an average for the second highest BLAT scores for supporting reads of <900; and (c) more than 50% of supporting reads (at least three reads). For brain-only samples (unpaired cases), a more strict depth of <300 and a VAF of ≥10% were applied. In addition to MuTect2 calling, we independently ran RePlow for both paired and unpaired cases (see below for more details) and applied the same filter conditions to the outputs. All putative SNVs in protein-coding regions were annotated with the Ensembl Variant Effect Predictor for the characterization of mutation subtypes using the following simplified categories: (a) LOF_MIS, defined as missense, stop gained, start lost, stop lost, splice donor, or splice acceptor variants, and (b) SYN, defined as synonymous variants. For pathogenicity scoring of variants, the Phred-scaled CADD score (v1.5, GRCh38 model) was used, and a cutoff of >20 was applied for identifying putatively pathogenic variants.

Somatic SNV calling powered by RePlow and barcoded sequencing

The reason for applying RePlow was to rescue SNV candidates filtered out from MuTect2 due to low VAFs. RePlow is a specialized tool designed to detect SNV candidates with low VAFs by utilizing information gleaned from replicates of the same sample. We assumed that if a given SNV candidate was consistently observed from multiple brain regions of the same individual, this would indicate a true variant, despite its low VAF. We considered data from multiple brain regions as pseudoreplicates and applied RePlow to recover such variants. To do so, we first generated initial call sets for each brain region using MuTect2. Then, for each possible pair of brain regions from the same individual (e.g., BA17-BA21), we applied RePlow and selected somatic SNVs identified as true mutations in the initial call sets, which we defined as replicate calls. Among the multiple replicate calls from at least two different region pairs (i.e., three different brain regions), not from the matched peripheral tissue, we confirmed genuine calls using independent ultradeep validation with high-depth sequencing with barcode tags (BCDseq) at an average sequencing read depth of 82,924x. For BCDseq calls, the filter conditions were as follows: (1) >4 supporting barcodes, (2) VAF concordance between WES and BCDseq within a 10-fold change and within a difference of 5%, and (3) conditions where (1) and (2) were satisfied for both regions of an identical variant.

Empirical Bayesian score

We applied EBscore27 to our raw variants from MuTect2 and RePlow calling to distinguish true calls from false-positives based on empirical estimation of the error rates of the variants. The EBscore performance was previously analyzed and described by our group elsewhere19. We utilized WES data from 21 independent healthy controls as a panel of reference samples with which to estimate the error rate of each variant position, and we set the EBscore cutoff value to 5 to be conservative in detecting accurate calls.

Targeted amplicon validation sequencing

Primers were designed with the Primer3 algorithm and synthesized by Macrogen (Seoul, Korea). Amplicons were prepared by two-step PCR using Illumina TruSeq adapters. First, PCRs were carried out using 5 ng of initial template gDNA. The 1st amplicons were analyzed on 2% agarose gels, and bands of the expected sizes were isolated and purified using Mega Quick-Spin Kits (iNtRON, Korea). After purification, 50 ng of the 1st amplicon was used as a template for a second PCR, and the products were subsequently purified with the same purification kit. These 2nd amplicons were quantified using the Bioanalyzer 2100 system (Agilent, USA). QC-passed amplicons were then sequenced on a HiSeq 2500 sequencer (Illumina, USA). The generated Fastq files were aligned to the GRCh38 reference genome by BWA-MEM28, and the reads at the target sites were filtered for MQ20 and BQ30 with the bam-readcount R package. We checked the number and quality of altered alleles via visualization with IGViewer29 (v2.3.94). For the estimation of background error rates, we performed replicate sequencing with previously constructed spike-in samples. If the identified mutations were statistically reliable (p < 0.05) compared to the estimated background error rates, we considered them true calls. The background error rates for the amplicon-based platforms were as follows: T > A (VAF = 0.00312), T > C (VAF = 0.00797), T > G (VAF = 0.000758), C > T (VAF = 0.00407), C > G (VAF = 0.000765), and C > A (VAF = 0.00185)25.

Germline SNV calling

We applied GATK HaplotypeCaller30 (v.3.8–0) for calling germline SNVs and included reliable variants in our analysis based on filtering conditions of a total depth ≥100 and an allele frequency ≥30%. We applied the filter conditions to calls obtained concurrently for all brain and peripheral specimens of each subject. To select ASD risk genes with rare and putatively damaging protein-altering germline SNVs, we applied additional conditions, including an ExAC allele frequency <0.02% and CADD score >20. The ASD risk genes were selected from the Simons Foundation Autism Research Initiative (SFARI) database31 if they were categorized as class 1 to 3.

Random permutation test

We collected all genes with putatively damaging brain somatic mutations and the known ASD risk genes (SFARI class 1 to 3) with damaging germline mutations from 24 ASD subjects and 31 normal controls, which totaled a set of 60 unique genes in our cohort. Then, we used the reference gene lists of axon guidance (GO:0007411), neuron projection guidance (GO:0097), dendrite spine morphogenesis (GO:0060997), cellular component morphogenesis (GO:0032989), calcium ion transport into cytosol (GO:0010524), intellectual disability (DisGeNET), mental retardation (DisGeNET), and intrauterine growth retardation (HP:0001511), all of which were identified in advance from EnrichR32 by analyzing our merged gene set of ASD subjects (n = 18 genes; 7 genes with damaging brain somatic SNVs validated via ultradeep targeted amplicon sequencing and 11 genes with damaging germline SNVs categorized as SFARI class 1 to 3). To simulate the ASD merged gene set overlapping each of the reference genes, we selected 18 random genes from the cohort gene set (n = 60) with 10,000 permutations (random resampling). From the permutation distribution of gene overlaps with the reference gene sets, we estimated significance by comparing the actual overlap count with an inferred overlap count cutoff of 5%. We repeated the same permutation test for the 13 merged genes with damaging brain somatic and germline SNVs from normal controls.

Mutation signature analysis

To determine the contribution of mutation signatures, we pooled all brain somatic SNVs from ASD patients and controls. Then, we formatted the pooled SNVs in VCF files and used them as input files for running Mutalisk33. The input files were compared with reference signatures generated by the tool using multiple likelihood estimation followed by linear regression. The best model of signature combinations for somatic SNVs was suggested from the tool by considering both cosine similarity and Bayesian information.

Protein-protein interactions

Protein‒protein interaction (PPI) datasets were downloaded from the Mentha project34. The dataset integrates PPI information where each protein is given a reliability score for its interaction with another protein. We summed the reliability scores of interactions for each protein to produce a weighted number of interactions per protein according to a previously published method35. We then compared the weighted numbers of interactions between different groups using two-sided Wilcoxon’s rank sum test.

Spatiotemporal gene expression analysis

The RNA-seq data of developmental gene expression for human brains were downloaded from BrainSpan36. The analysis was limited to 10 stages, from early fetal (10–13 weeks postconception) to middle adulthood (up to 45 years), and included expression values (RPKMs) for all 15 brain regions. We generated developmental signatures of gene expression for each of the 150 combination windows from 10 stages and 15 regions according to the overall formulations of a previous method37.

First, we calculated two z scores for each gene i in terms of stage s and region r. Here, eisr, which is the gene’s expression in a specific stage and region, was compared to the expression distributions across all stages at r to obtain zis and across all regions at s to obtain zir, respectively:zis=eisr−mis1.4826⋅MADisandzir=eisr−mir1.4826⋅MADir

where mis and mir are the median expression levels of gene i across all regions at stage s and across all stages at region r, respectively, and MADis and MADir are the median absolute deviations (MADs) of the gene expression across all regions at stage s and across all stages at region r, respectively. These zis and zir were then combined into a meta-z score:zisr=zis+zir2

Finally, the genes with zisr≥ 1.5 were utilized as the expression signatures for each of the 150 combination windows of stages and regions.

We then performed enrichment tests of the 150 spatiotemporal expression signatures for the merged genes (genes with damaging brain somatic SNVs validated via ultradeep targeted amplicon sequencing and genes with damaging germline SNVs categorized as SFARI class 1 to 3) in ASD subjects and normal controls. Fisher’s exact test was applied to calculate the significance of overlaps with all 150 expression signatures. P values were corrected for multiple tests using the Benjamini–Hochberg procedure.

Results

Identification of low-level brain somatic SNVs through deep sequencing

We performed high-depth WES (average throughput read-depth, ASD = 560.6x and control = 552.4x, respectively) on 137 postmortem cortical and cerebellar specimens and 44 peripheral specimens from 24 ASD subjects and 31 neurotypical controls (Supplementary Tables 1 and 2). Among them, 13 of the ASD subjects and all 31 controls had matched peripheral tissues (heart, liver, and spleen) that enabled SNV calling by paired analysis. For the quality control of the raw sequencing data, we confirmed that there was no contamination across the samples using ContEst26 and in-house filters (Supplementary Fig. 1a, b). We then established multiple variants calling pipelines optimized for sample conditions: i) paired or unpaired and ii) single or multiple brain regions for each subject (Fig. 1a). Initially, we detected 270 brain somatic SNVs, with an average VAF of 5.8%, in the deep WES data with the application of the somatic mutation caller Mutect224 and filtration strategies, including read depth, VAF cutoff, empirical Bayesian score27 and manual inspection. To rescue somatic SNVs with low VAFs under the detection limit of the Mutect2 pipeline, we further applied BCDseq and a replication-aware variant caller (RePlow)25 to samples from multiple brain regions to call 62 somatic SNVs with lower VAFs (average 2.7%) (Fig. 1a and Supplementary Fig. 2a). Among them, 54 (87.1%) were undetected in the Mutect2 analysis pipeline. We confirmed that the VAFs of the 54 mutations were concordant among multiple brain regions and also with those from BCDseq, excluding the risk of erroneous calls (Supplementary Fig. 2b). Overall, our comprehensive analysis identified 324 brain somatic SNVs (Supplementary Table 3).Fig. 1 Profiles of brain somatic single nucleotide variants (SNVs).

a Analytic pipeline. For paired and unpaired brain samples, SNVs were called using different callers and subjected to a unified postcall filter. In total, 324 analysis-ready brain somatic SNVs were identified. b Accuracy of the postcall filter. Using 52 random variants, we estimated the precision of the postcall filter to be 80.8%. c Concordance in variant allele frequencies (VAFs) between different sequencing platforms. Strong positive correlations were observed between the VAFs of BCDseq and TASeq and between WES and TASeq. d Similar VAFs among different brain regions. e Similar VAFs between autism spectrum disorder (ASD) patients and controls. Boxplots indicate the median and first and third quartiles; whiskers represent 1.5 times the interquartile range; red diamonds indicate mean VAFs. f Similar SNV counts between ASD subjects and controls across different brain regions. Error bars indicate standard errors.

For validation, we randomly selected and subjected 52 of the 324 somatic SNVs (16%) to ultradeep TASeq, which showed that the precision of detecting low-level brain somatic SNVs (average VAF = 2.4%, standard error = 1.84–3.01) was 80.8% (Fig. 1b and Supplementary Table 4). We then compared VAFs between the WES and TASeq datasets and between the BCDseq and TASeq datasets and found that the VAFs of the validated mutations were concordant between the different sequencing platforms (Fig. 1c). The average VAFs were also similar among different brain regions, including BA9, BA21, BA22, BA17, and the cerebellum (6.0%, 4.9%, 4.8%, 5.3%, and 4.7%, respectively) (Fig. 1d), and between ASD patients and control subjects (5.1% and 5.8%, respectively) (Fig. 1e).

The average numbers of brain somatic SNVs per single brain region were 2.3 for ASD subjects and 2.4 for controls (Fig. 1f). For ASD subjects, there was no sex bias in the number of brain somatic SNVs (average number of males = 1.73, average number of females = 2.50; p = 0.18, Student’s t test). To further investigate sex bias in brain somatic mutations, we collated diverse somatic mutation sources38 and found that there was still no female-enriched somatic mutation burden (somatic mutation numbers per brain per individual, males = 3, females = 2.27; p = 0.48, Student’s t test). We additionally compared the average VAFs of somatic SNVs between paired and unpaired brain samples and found no significant difference between the two (average VAFs, paired samples = 5.3% and unpaired samples = 5.2%, respectively) (Supplementary Fig. 1c). We found no differences among brain regions regarding average throughput read depths (Supplementary Fig. 1d, e). There was no evidence of significant bias according to brain region, paired tissue availability, or disease status. Additionally, regarding VAFs and mutation numbers, there was no significant difference between ASD patients and controls.

For all 55 subjects with an average age of 29.6 years, we estimated the presence of ~2.4 brain somatic SNVs per exome per single brain region, with VAFs as low as 0.3%. By extrapolating this information on a broader genomic scale (~75 million bps to ~3 billion bps), we determined that a total of 96 low-level somatic SNVs per brain may be present at the genome level in ASD patients and controls. Overall, our genetic analyses were able to identify 324 reliable brain somatic SNVs with an average VAF of 5.3% from 24 ASD subjects and 31 neurotypical controls.

Brain somatic SNVs in ASD patients exhibit comparable mutational loads but are enriched in genes associated with high PPI and neurodevelopment

To examine the potential mechanisms of mutagenesis underlying brain somatic SNVs, we performed mutation signature analysis of the 324 brain somatic SNVs detected in the 55 subjects. We utilized Mutalisk33 for the analysis and found that mutation signatures converged on signatures 1 and 5 with cosine similarities of 0.91 and 0.8 in ASD subjects and controls, respectively (Fig. 2a). Mutation signatures 1 and 5 were consistent among subsets of the brain somatic SNVs from different individual mutation callers (Supplementary Fig. 3). Generally, signatures 1 and 5 have been documented in all cancer types and in most cancer samples39,40, indicating that a spontaneous, endogenous mutational process may undergird the mutagenesis of the brain somatic SNVs in both ASD subjects and controls in our exome sequencing dataset.Fig. 2 Comparison of functional profiles of brain somatic SNVs between ASD patients and controls.

a Mutation signature analysis of brain somatic SNVs. Brain somatic SNVs converge on signatures 1 and 5, with cosine similarities of 0.91 and 0.80 in ASD subjects and controls, respectively. b Comparison of mutation burden between MIS_LOF and SYN. There was no biased enrichment of transcript length-normalized mutation counts between genes with MIS_LOF and those with SYN. Boxplots depict the median and first and third quartiles; whiskers represent 1.5 times the interquartile range. c Numbers of annotated SNVs per subject. The results indicate no difference between ASD subjects and controls for MIS_LOF or SYN. Error bars indicate standard errors. d gnomAD exome allele frequency distribution. At least 60% of the subjects were found to carry MIS_LOF, which has not been previously reported in the general population. e Protein‒protein interaction (PPI) levels for MIS_LOF. PPI levels for genes with MIS_LOF and CADD scores >20 were greater in ASD patients than in controls. Error bars indicate standard errors.

Next, we sought to determine whether the types of brain somatic SNVs differed between ASD patients and controls. To do this, we utilized the Ensembl Variant Effect Predictor (database version 97, GRCh38)41. A total of 36.4% of the somatic SNVs (118 out of 324) were annotated within protein-coding regions, among which we specifically chose SNVs presenting as loss-of-function (LOF), missense (MIS), and synonymous (SYN) mutations for further analyses. Between the LOF_MIS and SYN mutations, we examined the mutation burden per gene length42, ensuring no biased enrichment in nonsense or nonsynonymous mutations from our exome sequencing dataset (Fig. 2b). We found no significant differences in the number of LOF_MIS or SYN mutations per subject between the ASD patients and controls (Fig. 2c). Interestingly, more than 60% of the LOF_MIS somatic SNVs that we identified had not previously been reported in the genomAD exome allele frequency database43 (Fig. 2d). These results indicated that both average mutation counts and population-based minor allele frequencies for somatic SNVs did not significantly differ between ASD patients and controls.

Accordingly, we hypothesized that genes with nonsynonymous mutations would be functionally different between ASD patients and controls. In particular, we focused on genes with putatively damaging mutations, for which we applied Phred-scaled CADD scores44 >20 to MIS_LOF (MIS_LOF genes with CADD > 20, ASD = 19 of 51 (37.3%) and control = 12 of 30 (40%)). It was previously reported that highly damaging mutations in ASD-related genes are highly associated with protein‒protein interactions (PPIs)35. Generally, the PPI score of a gene indicates the degree of involvement “per gene” in diverse protein‒protein interactions. To obtain the single PPI score per gene, we calculated a weighted sum of the varying PPI scores assigned to a single gene. We then sought to compare PPI levels between ASD patients and controls for the genes with LOF_MIS mutations after classifying genes by high and low CADD scores. The median level of PPIs was greater in ASD subjects than in controls for LOF_MIS mutations with CADD scores >20 (ASD = 4.03 and control=1.93, p = 0.024, two-sided Wilcoxon rank sum test) (Fig. 2e). These results showed that damaging brain somatic SNVs in ASD patients are more significantly enriched for genes with high PPIs than in controls, perhaps leading to biological network dysfunction related to ASD.

Given the genetic heterogeneity underlying ASD, various mutational impacts may converge on ASD-related pathogenesis45. Thus, we examined whether brain somatic SNVs were found in genes implicated in neurodevelopmental processes or neuropsychiatric disorders. Using ultradeep TASeq, we were able to validate rare damaging brain somatic SNVs in seven different genes with CADD scores >20. Among them, five genes (ADCY5, CENPJ, DVL1, PEAK1, and RGS6) were not considered high-risk ASD candidate genes in the Simons Foundation Autism Research Initiative (SFARI) dataset31 (Table 1 and Supplementary Table 5).Table 1 Subjects carrying damaging brain somatic SNVs.

Subject (n = 6)	Group	Brain somatic SNVs	
Gene	Mutation	VAF (%)	gnomAD_AF (%)	Neurodevelopmental relevance	
4999	ASD	DVL1	NM_001330311.2:c.196 C > A (p.D66Y)	0.3	Not reported	Social recognition of hierarchy and dominance↓ (PMID: 14960015)	
ADCY5	NM_001199642.1:c.2101 G > A (p.H701Y)	1.0	Not reported	Dyskinesia↑ (PMID: 26537056)	
5144	ASD	ERBB3	NM_001982.3:c.1611 T > G (p.N537K)	2.0	Not reported	Social novelty preference↓

(PMID: 21547722)

	
5176	ASD	PEAK1	NM_024776.3:c.890 C > T (p.R297Q)	2.5	0.000402		
5403	ASD	RGS6	NM_001204424.2:c.335 G > A (p.R112H)	3.0	0.00199	Anxiety and Depression↑

(PMID: 24421401)

	
5308	ASD	SLC25A22	NM_001191060.1:c.718 C > T (p.A240T)	3.8	0.000401		
5841	ASD	CENPJ	NM_018451.5:c.3001 G > T (p.Q1001K)	1.7	Not reported	Microcephaly

(PMID: 15793586)

	

However, the newly identified genes were associated with defective neurodevelopmental or psychiatric phenotypes. ADCY5, known to produce cAMP, which regulates neuronal function, was reported to be related to dyskinesia or Parkinsonian-like motor dysfunction in an ADCY5-null mouse model46. De novo mutations in the 5’ donor splice site of ADCY5 were also shown to cause early-onset autosomal dominant chorea and dystonia in a three-generation family study, indicating that ADCY5 haploinsufficiency is involved in these movement disorders47.

A homozygous single-base deletion and missense mutation in CENPJ were found to cause microcephaly via centrosome aberrations, and this gene is normally responsible for protein localization to the spindle poles of mitotic cells during prenatal neurogenesis in mice48.

DVL1 was shown to be associated with deficits in the recognition of social hierarchy and dominance in a DVL1-null mouse model49. Conditional ablation of ERBB3 in the central nervous system was reported to result in a lack of social novelty preference in a mouse model50. RGS6 knockout in mice was found to inhibit anxiety and depression via serotonin-mediated activation of the 5-HT receptor-adenylyl cyclase axis51, suggesting that a gain of function of this gene might be implicated in neuropsychiatric phenotypes52.

Notably, three of the aforementioned genes, ADCY5, PEAK1, and RGS6, showed relatively high constraint levels of 0.14 (90% CI 0.08 –0.25), 0.29 (90% CI 0.2 –0.44), and 0.22 (90% CI 0.13 –0.4), respectively, which were calculated by the expected ratios for loss-of-function mutations based on the GenomAD database43 (Supplementary Table 6). Furthermore, when using a brain cell-type specificity database based on human single-cell RNA sequencing data53, ADCY5 and ERBB3 were found to be enriched in oligodendrocytes, while PEAK1 and RGS6 were enriched in microglia and neurons, respectively (Supplementary Table 6). Taken together, our results suggested that these genes with rare damaging somatic mutations may be associated with the neurodevelopmental pathogenesis of ASD.

Damaging somatic and germline SNVs are collectively implicated in the pathogenesis of ASD

Many ASD genetic studies have reproducibly shown that rare damaging germline SNVs in ASD risk genes contribute to the pathogenesis of ASD54–56. Thus, from the perspective of the heterogeneous genetic architecture of ASD, we sought to scrutinize whether genes with brain somatic SNVs collaboratively affect ASD-related pathogenesis, with known ASD risk genes carrying rare damaging germline SNVs (CADD score >20, ExAC <0.02%). To do this, we extracted data from subjects with putatively damaging germline SNVs in ASD risk genes among ASD patients and controls. We found that 50% (12 of 24) of the ASD subjects and 25.8% (8 of 31) of the controls had putatively damaging germline SNVs in the ASD risk genes categorized as SFARI ASD risk class 1 (strong), 2 (high), or 3 (suggestive) (Supplementary Table 7). We then merged the newly discovered genes with damaging brain somatic SNVs and known ASD risk genes with damaging germline SNVs in ASD patients to determine whether the merged set of genes converged on biological mechanisms related to ASD pathogenesis. The merged gene set consisted of the genes with damaging brain somatic SNVs validated via ultradeep targeted amplicon sequencing (gene count, ASD = 7 and control = 4, respectively) and the SFARI class 1 – 3 genes with damaging germline SNVs (gene count, ASD = 11 and control = 9, respectively). The merged gene sets thus included 18 genes from 24 ASD subjects and 13 genes from 31 controls. Given that pathogenic mutations are likely located in highly expressed genes involved in tissue development57, we hypothesized that the merged gene set with damaging SNVs in ASD patients could overlap with highly expressed genes during fetal brain development, which is associated with the crucial pathobiology of ASD58.

To test our hypothesis, we constructed 150 spatiotemporal signatures of gene expression from 15 brain regions and 10 developmental periods (Fig. 3a). After calculating how many of the merged gene sets (n = 18) overlapped each of the 150 signatures, we found that the genes of ASD patients exhibited significantly high expression levels during the fetal period across ASD-related brain regions, including the cortex59,60, thalamus61, and striatum62 (all p < 0.05, Benjamini–Hochberg correction). In controls, however, the merged gene set (n = 13) with somatic and germline SNVs did not show significantly enriched gene overlaps for signatures of fetal brains.Fig. 3 Damaging somatic and germline SNVs are collectively implicated in ASD-related pathobiology.

a Spatiotemporal gene expression analysis. Spatiotemporal expression analysis was performed for the merged genes (ASD = 18 and control=13) using 150 spatiotemporal gene expression signatures constructed from 15 brain regions and 10 developmental periods. The merged genes consisted of genes with damaging brain somatic SNVs validated via ultradeep targeted amplicon sequencing and genes with damaging germline SNVs categorized as SFARI class 1 – 3. Green and gray boxes indicate significant differences in the number of merged genes overlapping each of the 150 signatures (p < 0.05). P values were corrected for the Benjamini‒Hochberg procedure. b Gene enrichment analysis using random permutation. Gene enrichment analysis was performed for the merged genes (ASD = 18 and control = 13) via random 10,000 permutations. P values were adjusted for the false discovery rate using the Benjamini–Hochberg method. The merged genes consisted of genes with damaging brain somatic SNVs validated via ultradeep targeted amplicon sequencing and genes with damaging germline SNVs categorized as SFARI class 1–3. Reference gene lists of biological terms were obtained from DisGeNET97, Gene Ontology Biological Process98, and Human Phenotype Ontology99. Curating the reference gene lists did not involve any presumptions or the use of specific sets associated with brain function-relevant genes. Genes in the orange box and blue box indicate genes with damaging brain somatic SNVs and genes with damaging germline SNVs, respectively. The red arrow and blue vertical line represent the actual and estimated (cutoff=5%) counts of genes overlapping the reference gene lists, respectively.

To further validate our methodology, we also employed an independent gene expression analysis tool, CSEA63, and found similar results: the merged gene set in ASD subjects was enriched in the mid-fetal cortex (p = 0.021, Bonferroni correction), but the genes in controls were not in any of the spatiotemporal windows (Supplementary Fig. 4). There were no significant findings from the same analysis when only genes carrying germline mutations were used, probably due to the limited cohort size or the partial genetic accountability for the disease association. Taken together, these results indicated that brain somatic and germline SNVs might impact a collection of genes involved in the aberrant development of fetal brain subregions that are particularly crucial for the pathogenesis of ASD.

With the noted accumulation of damaging brain somatic and germline SNVs, we then analyzed biological terms using EnrichR32 and found ASD-related biological terms with adjusted p values < 0.05 (Supplementary Table 8). For all of the ASD-related biological terms, we further validated the significance of each of the gene enrichment results by performing 10,000 random permutations, considering a pool of all of the genes with damaging SNVs and the ASD-risk genes with SFARI class 1 –3 from our exome sequencing dataset: axon guidance (GO:0007411) (p = 0.004), neuron projection guidance (GO:0097485) (p = 0.012), dendrite spine morphogenesis (GO:0060997) (p = 0.025), cellular component morphogenesis (GO:0032989) (p = 0.019), calcium ion transport into cytosol (GO:0010524) (p = 0.012), intellectual disability (p = 0.004), mental retardation (p = 0.012), and intrauterine growth retardation (HP:0001511) (p = 0.012) (Fig. 3b). However, for the merged gene set found in controls, ASD-related or neural development-associated terms were not enriched after the permutation test. The biological process terms for the ASD subjects were reported to be closely related to the pathogenesis of ASD. Briefly, axon pathology, including disruptions in axon growth and projection, has been repeatedly reported in both ASD mouse models64,65 and ASD patients66. ASD patients exhibit reductions in the size and number of dendrites as well as altered dendrite morphology67,68. Aberrations in neuron projections in the human deep cortex during the fetal period have been documented as crucial pathogeneses of ASD69. Calcium ion transport activity was reported to be markedly increased in the brain tissues of autistic patients, implicating altered calcium homeostasis in ASD pathology70. ASD and other neurodevelopmental conditions, including intellectual disability71–73 and mental retardation74,75, are thought to have common genetic etiologies, neural circuit alterations, and brain abnormalities such as synaptic transmission. Finally, previous studies indicate that ASD is frequently accompanied by intrauterine growth retardation76–78, indicating that the disrupted growth of fetal organs is a critical pathogenesis of ASD.

To further validate the biological implications of the merged gene set in the ASD subjects, we additionally collected ASD-related genes of coexpression modules, which included 10,459 genes from the study of transcriptomic coordination in the developing human prefrontal cortex79. Notably, our genes carrying brain somatic and germline mutations still showed enrichment of coexpressed genes (p = 0.05, 10,000 random permutations). Overall, these results suggested that brain somatic SNVs may converge with germline SNVs on the pathogenic features observed in ASD.

Finally, we sought to evaluate the combined contribution of brain somatic and germline mutations in ASD patients compared to that in normal controls80. To do this, we curated novel genes with brain somatic SNVs, as well as ASD risk genes with damaging germline SNVs, and arrayed the subjects according to gene counts (Fig. 4a and Supplementary Fig. 5a). We found that 58.3% (14 of 24) and 25.8% (8 of 31) of the ASD subjects and controls, respectively, carried at least one gene with either a somatic or germline SNV. Then, we examined whether the collection of brain somatic and germline SNVs was associated with an increased risk of ASD. We compared the number of ASD subjects to that of controls according to the presence of genes with damaging somatic or germline mutations: genes with damaging brain somatic SNVs only (6 ASD subjects and 2 controls), genes with damaging germline SNVs only (12 ASD subjects and 8 controls), and the merged gene sets with damaging brain somatic and germline SNVs (14 ASD subjects and 8 controls) (Fig. 4b and Supplementary Fig. 5b). Interestingly, we found that compared with normal controls, ASD patients carrying the merged gene set with damaging brain somatic or germline SNVs had a significantly greater risk of developing ASD, compared to normal controls (p = 0.025, odds ratio = 3.92, 95% confidence interval=1.121 – 14.794, two-sided Fisher’s exact test) (Fig. 4c). However, due to the small size of our cohort (24 ASD individuals and 31 controls), we failed to observe a significant contribution of either genes with damaging brain somatic SNVs only or genes with damaging germline SNVs only to ASD. Taken together, these findings indicated that brain somatic SNVs, in addition to germline SNVs, underlie the heterogeneous genetic architecture of ASD.Fig. 4 Collective contribution of brain somatic SNVs and germline SNVs to ASD.

a Heterogeneous genetic architectures of ASD subjects. A total of 58.3% of the ASD subjects (14 of 24) were found to carry genes with damaging brain somatic SNVs, as validated via targeted amplicon sequencing, and/or genes with damaging germline SNVs, categorized as class 1 to 3. Among the 14 patients, 6 carried 2 or 3 genes with damaging SNVs. b Proportion of ASD subjects with impactful genes with damaging SNVs. Among the 24 ASD subjects, 14 (58.3%) were found to carry genes with damaging somatic and/or germline SNVs. c Assessment of ASD risk was performed by comparing the gene counts of the subsets of ASD subjects (n = 24) and normal controls (n = 31), for the genes with damaging brain somatic SNVs only (6 ASD subjects and 2 controls), the genes with damaging germline SNVs only (12 ASD subjects and 8 controls), and the merged gene set with damaging brain somatic SNVs and damaging germline SNVs (14 ASD subjects and 8 controls). ASD patients carrying the merged gene set with damaging brain somatic and germline SNVs had a significantly greater risk of developing ASD (odds ratio = 3.92, p = 0.025, 95% confidence interval = 1.12 –14.79) than did controls.

Discussion

In the present study, we generated and exploited a unique resource of deep WES data from 181 postmortem brain and peripheral tissue samples from 24 ASD subjects and 31 neurotypical controls to identify low-level somatic mutations and their contribution to ASD. In doing so, we discovered that brain somatic SNVs in novel genes implicated in neurodevelopmental processes, along with germline SNVs in well-known ASD risk genes, and the merged gene set identified from both somatic and germline mutation sources are associated with the pathogenesis of ASD. The merged gene set affected by either the brain somatic or germline SNVs demonstrated significant enrichment for biological process terms, including axon guidance, neuron projection, and dendrite spine morphogenesis, as well as for other neurodevelopmental condition terms, including intellectual disability, mental retardation, and intrauterine growth retardation. The merged gene set also showed higher gene expression levels during fetal brain development and greater mutation burden in ASD subjects than in normal controls.

In ASD patients, we identified deleterious brain somatic SNVs in ADCY5, CENPJ, DVL1, ERBB3, PEAK1, RGS6, and SLC25A22. Interestingly, five of these genes had previously been found to affect neural commitments: ADCY5, which is highly expressed in the striatum, is involved in the modulation of dopaminergic signals and is thus tightly related to motor control81. Indeed, in the present study, a subject (ID4999) with a brain somatic SNV in ADCY5 was reported to suffer from involuntary rhythmic movements. CENPJ regulates neural progenitor division and neuronal migration in the cerebral cortex82. DVL1 is required for normal dendritic development in hippocampal neurons83,84. ERBB3 is known to modulate hippocampal neuroplasticity85,86. RGS6 is known to suppress dopaminergic neurodegeneration and motor dysfunction87. The other two genes, however, have undergone little study in regard to their neural functions and their relationships with neuropsychiatric conditions. The studies that are available have reported that these genes are closely associated with the pathologic features of ASD. PEAK1 is known to be a regulator of cell migration88, and aberrations in neuronal migration have been documented in individuals with ASD89,90. SLC25A22 constitutes a mitochondrial glutamate carrier91, and disruption of mitochondrial carriers by altered calcium signaling has been implicated in ASD pathogenesis92.

Notably, the percentage of VAFs of the brain somatic SNVs found in these genes ranged from 0.3 to 3.8%. These low-level somatic mutation burdens in the brain imply that some cells of specific types that are critical for network activity or oscillation, such as inhibitory neurons93, or cell nonautonomous mechanisms94, may cause defects in the entire brain. How these genes and brain somatic SNVs lead to neural dysfunction or ASD phenotypes remains unclear and warrants further study.

The detection of low-level somatic SNVs generally involves many false-positive calls95. To prevent this, we performed high-depth WES at an average throughput read depth of 559.3x. Previous research using the same brain samples only reached an average read depth of up to 95x in target regions and reported negative findings in detecting significant somatic SNVs in ASD brains21. A few studies performed only targeted sequencing on confined sets of known ASD risk genes, with no further investigation of novel genes20. Additionally, recent studies involving relatively low-depth (~250×) whole-genome sequencing have suggested that somatic mutations in enhancers or mosaic copy number variations may be associated with ASD22,23. Although these studies revealed somatic mutations in noncoding regions in ASD patients, low-level somatic mutations in protein-coding regions are generally difficult to detect because of the low read depth of the current WGS approach. Compared with those studies, the high-depth (559.3x) sequencing approach in the present study allowed us to spot low-level somatic mutations with a VAF as low as 0.3%. Furthermore, low-level somatic mutations were found in neurodevelopmental genes that are possibly associated with the pathobiology of ASD (VAF, DVL1 = 0.3%, ADCY5 = 1.0%, CENPJ = 1.7%, ERBB3 = 2.0%, PEAK1 = 2.5%, RGS6 = 3.0%, and SLC25A22 = 3.8%). Our results thus support that high-depth sequencing is necessary for the detection of brain somatic SNVs with low-level allele frequencies in novel genes, as well as in known ASD risk genes.

By using deep WES followed by strict validation methods, we identified reliable somatic mutations that putatively affect genes key to neural dysfunction related to ASD. To understand the comprehensive genetic architecture of ASD, we also explored rare damaging germline mutations, which are plausibly considered to increase the risk of developing ASD. The approach of combining somatic and germline mutations can in part explain the complex features of the heterogeneous genetic architectures implicated in ASD. Specifically, our results support the idea that the accumulation of brain somatic SNVs during early embryonic or brain development likely contributes to the neural dysregulation and related pathogenesis of ASD in conjunction with predisposed germline SNVs96. Both ASD patients and controls were found to have germline mutations in ASD risk genes at their earliest stage of life, embryos (Fig. 5). As embryos undergo cell division, somatic mutations arise and accumulate in genetic architectures that are already predisposed with germline mutations in ASD risk genes. Consistent with this idea, our unsupervised hierarchical clustering analysis, which considered various features, such as the cell division timing of the somatic mutations, the presence of damaging somatic mutations and damaging germline mutations related to SFARI ASD risk, and neurodevelopmental relevance, showed that the ASD subjects (ID 5144, 4999, and 5841; ID 4917 and 4899) were distinctively clustered, whereas the controls were not (Fig. 5). These findings additionally support the hypothesis that somatic and germline mutations can be collectively implicated in the heterogeneous architecture of ASD, although this needs to be replicated in a larger cohort. Nonetheless, the observed biological convergence between somatic and germline SNVs is consistent with an oligogenic model of ASD96, wherein diverse mutational sources collectively contribute to ASD risk and its genetic heterogeneity. We cautiously address the oligogenic model-based approach we adopted, which does not mean that mutations occur at the same time but rather that they accumulate and contribute to the overall incidence of ASD. Future studies using a larger ASD cohort with postmortem brain samples will be required to further assess the full contribution of brain somatic mutations to ASD.Fig. 5 The accumulation of somatic mutations during embryonic cell division.

Germline mutations are present from the beginning of embryogenesis. Somatic mutations spontaneously arise in multiple brain regions as embryonic cell division progresses, thereby leading to the accumulation of mutations with low-level variant allele frequencies (VAFs) in the genetic architectures of ASD patients. Prepredisposed germline mutations and accumulating somatic mutations can be simultaneously present in the subjects, thus possibly leading to collective damage to ASD-related neurodevelopmental processes, which is in line with the genetic heterogeneity model of ASD. To differentiate between the ASD patients and controls, unsupervised hierarchical clustering analysis was further performed using features such as the cell division timing of the somatic mutations, the presence of damaging somatic mutations, the presence of damaging germline mutations related to the risk of SFARI ASD, and the neurodevelopmental relevance of the damaging somatic mutations. Some of the ASD subjects (ID 5144, 4999, and 5841; ID 4917 and 4899) were found to be distinctively clustered, whereas the controls were not. For a clear depiction, we used a subset of brain somatic mutations identified from an analysis pipeline using Replow and BCDseq.

Supplementary information

Supplementary Information

Supplementary Table 1. Sample information summary

Supplementary Table 2. Whole-exome sequencing information of 181 post-mortem brain and peripheral tissues

Supplementary Table 3. High-confidence brain somatic SNVs

Supplementary Table 4. Somatic SNV validation using targeted amplicon sequencing

Supplementary Table 8. Gene enrichment analysis using EnrichR for the merged gene set with damaging brain somatic and germline SNVs in ASD subjects

Supplementary information

The online version contains supplementary material available at 10.1038/s12276-024-01284-1.

Acknowledgements

We gratefully acknowledge the National Institute of Child Health and Human Development (NICHD) and the Stanley Medical Research Institute (SMRI) for providing brain or peripheral tissue samples. This work was supported by grants from the Suh Kyungbae Foundation (to J.H.L.) and from the National Research Foundation of Korea (NRF) funded by the Korean government, Ministry of Science and ICT (NRF-2019R1A3B2066619 to J.H.L., RS-2023-00278314 to Sw.K., NRF-2022R1C1C1010430 to J.K., NRF-2021M3A9E4080784 to I.B.K., and NRF-2021M3E5D9021878 to I.B.K.).

Author contributions

Study conceptualization, J.H.L., I.B.K.; organizing experimental design, I.B.K.; data generation, I.B.K., M.K., W.K.K., M.J.W., Sh.K. and J.H.K.; data processing. I.B.K., M.K., S.J., J.H.K., H.J.K. and J.K.; annotation of functional regions, I.B.K.; data analysis, I.B.K., M.K., S.J., Y.S.J. and J.L.; statistical analysis, I.B.K., S.J. and J.K.; manuscript preparation, I.B.K., J.K., Sw.K., and J.H.L.; and leading project, J.H.L.

Data availability

The raw high-depth whole-exome sequencing data were deposited in K-BDS (Korea BioData Station) with the accession ID KAP240684.

Competing interests

J.H.L. is a co-founder and CTO of SoVarGen, Inc., which seeks to develop new diagnostics and therapeutics for brain disorders. The remaining authors declare no competing financial interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Il Bin Kim, Myeong-Heui Kim.
==== Refs
References

1. De Rubeis S Buxbaum JD Genetics and genomics of autism spectrum disorder: embracing complexity Hum. Mol .Genet. 2015 24 R24 R31 26188008
De Rubeis, S. & Buxbaum, J. D. Genetics and genomics of autism spectrum disorder: embracing complexity. Hum. Mol .Genet. 24, R24–R31 (2015).26188008
2. Choi L An JY Genetic architecture of autism spectrum disorder: Lessons from large-scale genomic studies Neurosci. Biobehav. Rev. 2021 128 244 257 10.1016/j.neubiorev.2021.06.028 34166716
Choi, L. & An, J. Y. Genetic architecture of autism spectrum disorder: Lessons from large-scale genomic studies. Neurosci. Biobehav. Rev. 128, 244–257, 10.1016/j.neubiorev.2021.06.028 (2021).34166716 10.1016/j.neubiorev.2021.06.028
3. Kim IB Non-coding de novo mutations in chromatin interactions are implicated in autism spectrum disorder Mol. Psychiatry 2022 27 4680 4694 10.1038/s41380-022-01697-2 35840799
Kim, I. B. et al. Non-coding de novo mutations in chromatin interactions are implicated in autism spectrum disorder. Mol. Psychiatry 27, 4680–4694, 10.1038/s41380-022-01697-2 (2022).35840799 10.1038/s41380-022-01697-2
4. Ronald A Genetic heterogeneity between the three components of the autism spectrum: a twin study J. Am. Acad. Child Adolesc. Psychiatry 2006 45 691 699 16721319
Ronald, A. et al. Genetic heterogeneity between the three components of the autism spectrum: a twin study. J. Am. Acad. Child Adolesc. Psychiatry 45, 691–699 (2006).16721319
5. Gaugler T Most genetic risk for autism resides with common variation Nat Genet 2014 46 881 25038753
Gaugler, T. et al. Most genetic risk for autism resides with common variation. Nat Genet 46, 881 (2014).25038753
6. Bai D Association of genetic and environmental factors with autism in a 5-country cohort JAMA Psychiatry 2019 76 1035 1043 31314057
Bai, D. et al. Association of genetic and environmental factors with autism in a 5-country cohort. JAMA Psychiatry 76, 1035–1043 (2019).31314057
7. Iossifov I The contribution of de novo coding mutations to autism spectrum disorder Nature 2014 515 216 25363768
Iossifov, I. et al. The contribution of de novo coding mutations to autism spectrum disorder. Nature 515, 216 (2014).25363768
8. Krumm N Excess of rare, inherited truncating mutations in autism Nat. Genet 2015 47 582 25961944
Krumm, N. et al. Excess of rare, inherited truncating mutations in autism. Nat. Genet 47, 582 (2015).25961944
9. Gratten J Visscher PM Mowry BJ Wray NR Interpreting the role of de novo protein-coding mutations in neuropsychiatric disease Nat. Genet 2013 45 234 23438595
Gratten, J., Visscher, P. M., Mowry, B. J. & Wray, N. R. Interpreting the role of de novo protein-coding mutations in neuropsychiatric disease. Nat. Genet 45, 234 (2013).23438595
10. De Rubeis S Synaptic, transcriptional and chromatin genes disrupted in autism Nature 2014 515 209 25363760
De Rubeis, S. et al. Synaptic, transcriptional and chromatin genes disrupted in autism. Nature 515, 209 (2014).25363760
11. O’Roak BJ Exome sequencing in sporadic autism spectrum disorders identifies severe de novo mutations Nat. Genet 2011 43 585 21572417
O’Roak, B. J. et al. Exome sequencing in sporadic autism spectrum disorders identifies severe de novo mutations. Nat. Genet 43, 585 (2011).21572417
12. Lim ET Rates, distribution and implications of postzygotic mosaic mutations in autism spectrum disorder Nat. Neurosci 2017 20 1217 28714951
Lim, E. T. et al. Rates, distribution and implications of postzygotic mosaic mutations in autism spectrum disorder. Nat. Neurosci 20, 1217 (2017).28714951
13. Poduri A Evrony GD Cai X Walsh CA Somatic mutation, genomic variation, and neurological disease Science 2013 341 1237758 10.1126/science.1237758 23828942
Poduri, A., Evrony, G. D., Cai, X. & Walsh, C. A. Somatic mutation, genomic variation, and neurological disease. Science 341, 1237758, 10.1126/science.1237758 (2013).23828942 10.1126/science.1237758
14. Lee JH Somatic mutations in disorders with disrupted brain connectivity Exp. Mol. Med. 2016 48 e239 10.1038/emm.2016.53 27282107
Lee, J. H. Somatic mutations in disorders with disrupted brain connectivity. Exp. Mol. Med. 48, e239, 10.1038/emm.2016.53 (2016).27282107 10.1038/emm.2016.53
15. Lim JS Brain somatic mutations in MTOR cause focal cortical dysplasia type II leading to intractable epilepsy Nat. Med. 2015 21 395 400 10.1038/nm.3824 25799227
Lim, J. S. et al. Brain somatic mutations in MTOR cause focal cortical dysplasia type II leading to intractable epilepsy. Nat. Med. 21, 395–400, 10.1038/nm.3824 (2015).25799227 10.1038/nm.3824
16. Pelorosso, C. et al. Somatic double-hit in MTOR and RPS6 in hemimegalencephaly with intractable epilepsy. Hum. Mol. Genet (2019). 10.1093/hmg/ddz194
17. Sim NS Precise detection of low-level somatic mutation in resected epilepsy brain tissue Acta Neuropathol. 2019 138 901 912 31377847
Sim, N. S. et al. Precise detection of low-level somatic mutation in resected epilepsy brain tissue. Acta Neuropathol. 138, 901–912 (2019).31377847
18. Kim M-H Low-level brain somatic mutations are implicated in schizophrenia Biol. Psychiatry 2021 90 35 46 33867114
Kim, M.-H. et al. Low-level brain somatic mutations are implicated in schizophrenia. Biol. Psychiatry 90, 35–46 (2021).33867114
19. Park JS Brain somatic mutations observed in Alzheimer’s disease associated with aging and dysregulation of tau phosphorylation Nat. Commun. 2019 10 3090 31300647
Park, J. S. et al. Brain somatic mutations observed in Alzheimer’s disease associated with aging and dysregulation of tau phosphorylation. Nat. Commun. 10, 3090 (2019).31300647
20. D’Gama AM Targeted DNA sequencing from autism spectrum disorder brains implicates multiple genetic mechanisms Neuron 2015 88 910 917 26637798
D’Gama, A. M. et al. Targeted DNA sequencing from autism spectrum disorder brains implicates multiple genetic mechanisms. Neuron 88, 910–917 (2015).26637798
21. Freed D Pevsner J The contribution of mosaic variants to autism spectrum disorder PLoS Genet 2016 12 e1006245 27632392
Freed, D. & Pevsner, J. The contribution of mosaic variants to autism spectrum disorder. PLoS Genet 12, e1006245 (2016).27632392
22. Rodin RE The landscape of somatic mutation in cerebral cortex of autistic and neurotypical individuals revealed by ultra-deep whole-genome sequencing Nat. Neurosci. 2021 24 176 185 33432195
Rodin, R. E. et al. The landscape of somatic mutation in cerebral cortex of autistic and neurotypical individuals revealed by ultra-deep whole-genome sequencing. Nat. Neurosci. 24, 176–185 (2021).33432195
23. Sherman MA Large mosaic copy number variations confer autism risk Nat. Neurosci. 2021 24 197 203 33432194
Sherman, M. A. et al. Large mosaic copy number variations confer autism risk. Nat. Neurosci. 24, 197–203 (2021).33432194
24. Cibulskis K Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples Nat. Biotechnol. 2013 31 213 23396013
Cibulskis, K. et al. Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat. Biotechnol. 31, 213 (2013).23396013
25. Kim J The use of technical replication for detection of low-level somatic mutations in next-generation sequencing Nat. Commun. 2019 10 1047 30837471
Kim, J. et al. The use of technical replication for detection of low-level somatic mutations in next-generation sequencing. Nat. Commun. 10, 1047 (2019).30837471
26. Cibulskis K ContEst: estimating cross-contamination of human samples in next-generation sequencing data Bioinformatics 2011 27 2601 2602 10.1093/bioinformatics/btr446 21803805
Cibulskis, K. et al. ContEst: estimating cross-contamination of human samples in next-generation sequencing data. Bioinformatics 27, 2601–2602, 10.1093/bioinformatics/btr446 (2011).21803805 10.1093/bioinformatics/btr446
27. Shiraishi Y An empirical Bayesian framework for somatic mutation detection from cancer genome sequencing data Nucleic Acids Res. 2013 41 e89 e89 23471004
Shiraishi, Y. et al. An empirical Bayesian framework for somatic mutation detection from cancer genome sequencing data. Nucleic Acids Res. 41, e89–e89 (2013).23471004
28. Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv:1303.3997 (2013).
29. Thorvaldsdóttir H Robinson JT Mesirov JP Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration Brief Bioinforma. 2013 14 178 192
Thorvaldsdóttir, H., Robinson, J. T. & Mesirov, J. P. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinforma. 14, 178–192 (2013).
30. Poplin, R. et al. Scaling accurate genetic variant discovery to tens of thousands of samples. BioRxiv 201178 (2017).
31. Abrahams BS SFARI Gene 2.0: a community-driven knowledgebase for the autism spectrum disorders (ASDs) Mol. Autism 2013 4 36 10.1186/2040-2392-4-36 24090431
Abrahams, B. S. et al. SFARI Gene 2.0: a community-driven knowledgebase for the autism spectrum disorders (ASDs). Mol. Autism 4, 36, 10.1186/2040-2392-4-36 (2013).24090431 10.1186/2040-2392-4-36
32. Kuleshov MV Enrichr: a comprehensive gene set enrichment analysis web server 2016 update Nucleic Acids Res. 2016 44 W90 W97 10.1093/nar/gkw377 27141961
Kuleshov, M. V. et al. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 44, W90–W97, 10.1093/nar/gkw377 (2016).27141961 10.1093/nar/gkw377
33. Lee J Mutalisk: a web-based somatic MUTation AnaLyIS toolKit for genomic, transcriptional and epigenomic signatures Nucleic Acids Res. 2018 46 W102 W108 10.1093/nar/gky406 29790943
Lee, J. et al. Mutalisk: a web-based somatic MUTation AnaLyIS toolKit for genomic, transcriptional and epigenomic signatures. Nucleic Acids Res. 46, W102–W108, 10.1093/nar/gky406 (2018).29790943 10.1093/nar/gky406
34. Calderone A Castagnoli L Cesareni G Mentha: a resource for browsing integrated protein-interaction networks Nat Methods 2013 10 690 23900247
Calderone, A., Castagnoli, L. & Cesareni, G. Mentha: a resource for browsing integrated protein-interaction networks. Nat Methods 10, 690 (2013).23900247
35. Shohat S Ben-David E Shifman S Varying intolerance of gene pathways to mutational classes explain genetic convergence across neuropsychiatric disorders Cell Rep. 2017 18 2217 2227 28249166
Shohat, S., Ben-David, E. & Shifman, S. Varying intolerance of gene pathways to mutational classes explain genetic convergence across neuropsychiatric disorders. Cell Rep. 18, 2217–2227 (2017).28249166
36. Kang HJ Spatio-temporal transcriptome of the human brain Nature 2011 478 483 22031440
Kang, H. J. et al. Spatio-temporal transcriptome of the human brain. Nature 478, 483 (2011).22031440
37. Krishnan A Genome-wide prediction and functional characterization of the genetic basis of autism spectrum disorder Nat. Neurosci. 2016 19 1454 27479844
Krishnan, A. et al. Genome-wide prediction and functional characterization of the genetic basis of autism spectrum disorder. Nat. Neurosci. 19, 1454 (2016).27479844
38. Kim JH Analysis of low-level somatic mosaicism reveals stage and tissue-specific mutational features in human development PLoS Genet 2022 18 e1010404 10.1371/journal.pgen.1010404 36121845
Kim, J. H. et al. Analysis of low-level somatic mosaicism reveals stage and tissue-specific mutational features in human development. PLoS Genet 18, e1010404, 10.1371/journal.pgen.1010404 (2022).36121845 10.1371/journal.pgen.1010404
39. Rahbari R Timing, rates and spectra of human germline mutation Nat. Genet 2016 48 126 133 10.1038/ng.3469 26656846
Rahbari, R. et al. Timing, rates and spectra of human germline mutation. Nat. Genet 48, 126–133, 10.1038/ng.3469 (2016).26656846 10.1038/ng.3469
40. Goldmann JM Parent-of-origin-specific signatures of de novo mutations Nat. Genet 2016 48 935 939 10.1038/ng.3597 27322544
Goldmann, J. M. et al. Parent-of-origin-specific signatures of de novo mutations. Nat. Genet 48, 935–939, 10.1038/ng.3597 (2016).27322544 10.1038/ng.3597
41. McLaren W The Ensembl Variant Effect Predictor Genome Biol. 2016 17 122 10.1186/s13059-016-0974-4 27268795
McLaren, W. et al. The Ensembl Variant Effect Predictor. Genome Biol. 17, 122, 10.1186/s13059-016-0974-4 (2016).27268795 10.1186/s13059-016-0974-4
42. Samocha KE A framework for the interpretation of de novo mutation in human disease Nat. Genet 2014 46 944 25086666
Samocha, K. E. et al. A framework for the interpretation of de novo mutation in human disease. Nat. Genet 46, 944 (2014).25086666
43. Karczewski KJ The mutational constraint spectrum quantified from variation in 141,456 humans Nature 2020 581 434 443 32461654
Karczewski, K. J. et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 581, 434–443 (2020).32461654
44. Rentzsch P Witten D Cooper GM Shendure J Kircher M CADD: predicting the deleteriousness of variants throughout the human genome Nucleic Acids Res. 2019 47 D886 D894 10.1093/nar/gky1016 30371827
Rentzsch, P., Witten, D., Cooper, G. M., Shendure, J. & Kircher, M. CADD: predicting the deleteriousness of variants throughout the human genome. Nucleic Acids Res. 47, D886–D894, 10.1093/nar/gky1016 (2019).30371827 10.1093/nar/gky1016
45. Pinto D Convergence of genes and cellular pathways dysregulated in autism spectrum disorders Am. J Hum. Genet. 2014 94 677 694 24768552
Pinto, D. et al. Convergence of genes and cellular pathways dysregulated in autism spectrum disorders. Am. J Hum. Genet. 94, 677–694 (2014).24768552
46. Iwamoto T Motor dysfunction in type 5 adenylyl cyclase-null mice J Biol Chem 2003 278 16936 16940 10.1074/jbc.C300075200 12665504
Iwamoto, T. et al. Motor dysfunction in type 5 adenylyl cyclase-null mice. J Biol Chem 278, 16936–16940, 10.1074/jbc.C300075200 (2003).12665504 10.1074/jbc.C300075200
47. Carapito R A de novo ADCY5 mutation causes early-onset autosomal dominant chorea and dystonia Mov. Disord. 2015 30 423 427 10.1002/mds.26115 25545163
Carapito, R. et al. A de novo ADCY5 mutation causes early-onset autosomal dominant chorea and dystonia. Mov. Disord. 30, 423–427, 10.1002/mds.26115 (2015).25545163 10.1002/mds.26115
48. Bond J A centrosomal mechanism involving CDK5RAP2 and CENPJ controls brain size Nat Genet 2005 37 353 15793586
Bond, J. et al. A centrosomal mechanism involving CDK5RAP2 and CENPJ controls brain size. Nat Genet 37, 353 (2005).15793586
49. Long J LaPorte P Paylor R Wynshaw‐Boris A Expanded characterization of the social interaction abnormalities in mice lacking Dvl1 Genes, Brain Behav 2004 3 51 62 14960015
Long, J., LaPorte, P., Paylor, R. & Wynshaw‐Boris, A. Expanded characterization of the social interaction abnormalities in mice lacking Dvl1. Genes, Brain Behav 3, 51–62 (2004).14960015
50. Moy SS Deficient NRG1-ERBB signaling alters social approach: relevance to genetic mouse models of schizophrenia J Neurodev. Disord. 2009 1 302 21547722
Moy, S. S. et al. Deficient NRG1-ERBB signaling alters social approach: relevance to genetic mouse models of schizophrenia. J Neurodev. Disord. 1, 302 (2009).21547722
51. Stewart A Regulator of G-protein signaling 6 (RGS6) promotes anxiety and depression by attenuating serotonin-mediated activation of the 5-HT1A receptor-adenylyl cyclase axis FASEB J 2014 28 1735 1744 24421401
Stewart, A. et al. Regulator of G-protein signaling 6 (RGS6) promotes anxiety and depression by attenuating serotonin-mediated activation of the 5-HT1A receptor-adenylyl cyclase axis. FASEB J 28, 1735–1744 (2014).24421401
52. Coban-Akdemir Z Identifying genes whose mutant transcripts cause dominant disease traits by potential gain-of-function alleles Am. J Hum. Genet 2018 103 171 187 10.1016/j.ajhg.2018.06.009 30032986
Coban-Akdemir, Z. et al. Identifying genes whose mutant transcripts cause dominant disease traits by potential gain-of-function alleles. Am. J Hum. Genet 103, 171–187, 10.1016/j.ajhg.2018.06.009 (2018).30032986 10.1016/j.ajhg.2018.06.009
53. McKenzie AT Brain cell type specific gene expression and co-expression network architectures Sci. Rep. 2018 8 1 19 29311619
McKenzie, A. T. et al. Brain cell type specific gene expression and co-expression network architectures. Sci. Rep. 8, 1–19 (2018).29311619
54. Turner TN Genomic patterns of de novo mutation in simplex autism Cell 2017 171 710 722 28965761
Turner, T. N. et al. Genomic patterns of de novo mutation in simplex autism. Cell 171, 710–722 (2017). e712.28965761
55. RK CY Whole genome sequencing resource identifies 18 new candidate genes for autism spectrum disorder Nat. Neurosci. 2017 20 602 611 10.1038/nn.4524 28263302
RK, C. Y. et al. Whole genome sequencing resource identifies 18 new candidate genes for autism spectrum disorder. Nat. Neurosci. 20, 602–611, 10.1038/nn.4524 (2017).28263302 10.1038/nn.4524
56. Yuen RK Genome-wide characteristics of de novo mutations in autism NPJ Genom. Med. 2016 1 16027
Yuen, R. K. et al. Genome-wide characteristics of de novo mutations in autism. NPJ Genom. Med. 1, 16027 (2016).
57. Ding J Systematic analysis of somatic mutations impacting gene expression in 12 tumour types Nat. Commun. 2015 6 8554 10.1038/ncomms9554 26436532
Ding, J. et al. Systematic analysis of somatic mutations impacting gene expression in 12 tumour types. Nat. Commun. 6, 8554, 10.1038/ncomms9554 (2015).26436532 10.1038/ncomms9554
58. Courchesne E Mapping early brain development in autism Neuron 2007 56 399 413 10.1016/j.neuron.2007.10.016 17964254
Courchesne, E. et al. Mapping early brain development in autism. Neuron 56, 399–413, 10.1016/j.neuron.2007.10.016 (2007).17964254 10.1016/j.neuron.2007.10.016
59. Masuda F Motor cortex excitability and inhibitory imbalance in autism spectrum disorder assessed with transcranial magnetic stimulation: a systematic review Transl. Psychiatry 2019 9 110 10.1038/s41398-019-0444-3 30846682
Masuda, F. et al. Motor cortex excitability and inhibitory imbalance in autism spectrum disorder assessed with transcranial magnetic stimulation: a systematic review. Transl. Psychiatry 9, 110, 10.1038/s41398-019-0444-3 (2019).30846682 10.1038/s41398-019-0444-3
60. Khan S Somatosensory cortex functional connectivity abnormalities in autism show opposite trends, depending on direction and spatial scale Brain 2015 138 1394 1409 10.1093/brain/awv043 25765326
Khan, S. et al. Somatosensory cortex functional connectivity abnormalities in autism show opposite trends, depending on direction and spatial scale. Brain 138, 1394–1409, 10.1093/brain/awv043 (2015).25765326 10.1093/brain/awv043
61. Tomasi D Volkow ND Reduced local and increased long-range functional connectivity of the thalamus in autism spectrum disorder Cereb. Cortex 2019 29 573 585 10.1093/cercor/bhx340 29300843
Tomasi, D. & Volkow, N. D. Reduced local and increased long-range functional connectivity of the thalamus in autism spectrum disorder. Cereb. Cortex 29, 573–585, 10.1093/cercor/bhx340 (2019).29300843 10.1093/cercor/bhx340
62. Schuetze M Morphological alterations in the thalamus, striatum, and pallidum in autism spectrum disorder Neuropsychopharmacology 2016 41 2627 2637 10.1038/npp.2016.64 27125303
Schuetze, M. et al. Morphological alterations in the thalamus, striatum, and pallidum in autism spectrum disorder. Neuropsychopharmacology 41, 2627–2637, 10.1038/npp.2016.64 (2016).27125303 10.1038/npp.2016.64
63. Dougherty JD Schmidt EF Nakajima M Heintz N Analytical approaches to RNA profiling data for the identification of genes enriched in specific cells Nucleic Acids Res. 2010 38 4218 4230 10.1093/nar/gkq130 20308160
Dougherty, J. D., Schmidt, E. F., Nakajima, M. & Heintz, N. Analytical approaches to RNA profiling data for the identification of genes enriched in specific cells. Nucleic Acids Res. 38, 4218–4230, 10.1093/nar/gkq130 (2010).20308160 10.1093/nar/gkq130
64. Zikopoulos B Barbas H Changes in prefrontal axons may disrupt the network in autism J. Neurosci 2010 30 14595 14609 21048117
Zikopoulos, B. & Barbas, H. Changes in prefrontal axons may disrupt the network in autism. J. Neurosci. 30, 14595–14609 (2010).21048117
65. Xu Q Autism-associated CHD8 deficiency impairs axon development and migration of cortical neurons Mol. Autism 2018 9 65 10.1186/s13229-018-0244-2 30574290
Xu, Q. et al. Autism-associated CHD8 deficiency impairs axon development and migration of cortical neurons. Mol. Autism 9, 65, 10.1186/s13229-018-0244-2 (2018).30574290 10.1186/s13229-018-0244-2
66. Lazar M Miles LM Babb JS Donaldson JB Axonal deficits in young adults with High Functioning Autism and their impact on processing speed NeuroImage: Clin 2014 4 417 425 24624327
Lazar, M., Miles, L. M., Babb, J. S. & Donaldson, J. B. Axonal deficits in young adults with High Functioning Autism and their impact on processing speed. NeuroImage: Clin 4, 417–425 (2014).24624327
67. Martínez‐Cerdeño V Dendrite and spine modifications in autism and related neurodevelopmental disorders in patients and animal models Dev. Neurobiol. 2017 77 393 404 27390186
Martínez‐Cerdeño, V. Dendrite and spine modifications in autism and related neurodevelopmental disorders in patients and animal models. Dev. Neurobiol. 77, 393–404 (2017).27390186
68. Gouder L Altered spinogenesis in iPSC-derived cortical neurons from patients with autism carrying de novo SHANK3 mutations Sci. Rep. 2019 9 94 10.1038/s41598-018-36993-x 30643170
Gouder, L. et al. Altered spinogenesis in iPSC-derived cortical neurons from patients with autism carrying de novo SHANK3 mutations. Sci. Rep. 9, 94, 10.1038/s41598-018-36993-x (2019).30643170 10.1038/s41598-018-36993-x
69. Willsey AJ Coexpression networks implicate human midfetal deep cortical projection neurons in the pathogenesis of autism Cell 2013 155 997 1007 10.1016/j.cell.2013.10.020 24267886
Willsey, A. J. et al. Coexpression networks implicate human midfetal deep cortical projection neurons in the pathogenesis of autism. Cell 155, 997–1007, 10.1016/j.cell.2013.10.020 (2013).24267886 10.1016/j.cell.2013.10.020
70. Gargus JJ Genetic calcium signaling abnormalities in the central nervous system: seizures, migraine, and autism Ann. N. Y. Acad. Sci 2009 1151 133 156 10.1111/j.1749-6632.2008.03572.x 19154521
Gargus, J. J. Genetic calcium signaling abnormalities in the central nervous system: seizures, migraine, and autism. Ann. N. Y. Acad. Sci. 1151, 133–156, 10.1111/j.1749-6632.2008.03572.x (2009).19154521 10.1111/j.1749-6632.2008.03572.x
71. Abekhoukh S Bardoni B CYFIP family proteins between autism and intellectual disability: links with Fragile X syndrome Front. Cell Neurosci 2014 8 81 10.3389/fncel.2014.00081 24733999
Abekhoukh, S. & Bardoni, B. CYFIP family proteins between autism and intellectual disability: links with Fragile X syndrome. Front. Cell Neurosci. 8, 81, 10.3389/fncel.2014.00081 (2014).24733999 10.3389/fncel.2014.00081
72. van Bon BW Disruptive de novo mutations of DYRK1A lead to a syndromic form of autism and ID Mol. Psychiatry 2016 21 126 132 10.1038/mp.2015.5 25707398
van Bon, B. W. et al. Disruptive de novo mutations of DYRK1A lead to a syndromic form of autism and ID. Mol. Psychiatry 21, 126–132, 10.1038/mp.2015.5 (2016).25707398 10.1038/mp.2015.5
73. Langridge AT Maternal conditions and perinatal characteristics associated with autism spectrum disorder and intellectual disability PLoS One 2013 8 e50963 10.1371/journal.pone.0050963 23308096
Langridge, A. T. et al. Maternal conditions and perinatal characteristics associated with autism spectrum disorder and intellectual disability. PLoS One 8, e50963, 10.1371/journal.pone.0050963 (2013).23308096 10.1371/journal.pone.0050963
74. Guilmatre A Recurrent rearrangements in synaptic and neurodevelopmental genes and shared biologic pathways in schizophrenia, autism, and mental retardation Arch. Gen. Psychiatry 2009 66 947 956 10.1001/archgenpsychiatry.2009.80 19736351
Guilmatre, A. et al. Recurrent rearrangements in synaptic and neurodevelopmental genes and shared biologic pathways in schizophrenia, autism, and mental retardation. Arch. Gen. Psychiatry 66, 947–956, 10.1001/archgenpsychiatry.2009.80 (2009).19736351 10.1001/archgenpsychiatry.2009.80
75. Berkel S Mutations in the SHANK2 synaptic scaffolding gene in autism spectrum disorder and mental retardation Nat. Genet 2010 42 489 491 10.1038/ng.589 20473310
Berkel, S. et al. Mutations in the SHANK2 synaptic scaffolding gene in autism spectrum disorder and mental retardation. Nat. Genet 42, 489–491, 10.1038/ng.589 (2010).20473310 10.1038/ng.589
76. Alesi V 335.4 kb microduplication in chromosome band Xp11.2p11.3 associated with developmental delay, growth retardation, autistic disorder and dysmorphic features Gene 2012 505 384 387 10.1016/j.gene.2012.05.031 22634100
Alesi, V. et al. 335.4 kb microduplication in chromosome band Xp11.2p11.3 associated with developmental delay, growth retardation, autistic disorder and dysmorphic features. Gene 505, 384–387, 10.1016/j.gene.2012.05.031 (2012).22634100 10.1016/j.gene.2012.05.031
77. Pyhälä R Very low birth weight, infant growth, and autism-spectrum traits in adulthood Pediatrics 2014 134 1075 1083 25367538
Pyhälä, R. et al. Very low birth weight, infant growth, and autism-spectrum traits in adulthood. Pediatrics 134, 1075–1083 (2014).25367538
78. Lampi KM Risk of autism spectrum disorders in low birth weight and small for gestational age infants J. Pediatr. 2012 161 830 836 10.1016/j.jpeds.2012.04.058 22677565
Lampi, K. M. et al. Risk of autism spectrum disorders in low birth weight and small for gestational age infants. J. Pediatr. 161, 830–836, 10.1016/j.jpeds.2012.04.058 (2012).22677565 10.1016/j.jpeds.2012.04.058
79. Werling DM Whole-Genome and RNA sequencing reveal variation and transcriptomic coordination in the developing human prefrontal cortex Cell Rep. 2020 31 107489 10.1016/j.celrep.2020.03.053 32268104
Werling, D. M. et al. Whole-Genome and RNA sequencing reveal variation and transcriptomic coordination in the developing human prefrontal cortex. Cell Rep. 31, 107489, 10.1016/j.celrep.2020.03.053 (2020).32268104 10.1016/j.celrep.2020.03.053
80. de la Torre-Ubieta L Won H Stein JL Geschwind DH Advancing the understanding of autism disease mechanisms through genetics Nat. Med. 2016 22 345 27050589
de la Torre-Ubieta, L., Won, H., Stein, J. L. & Geschwind, D. H. Advancing the understanding of autism disease mechanisms through genetics. Nat. Med. 22, 345 (2016).27050589
81. Lee KW Impaired D2 dopamine receptor function in mice lacking type 5 adenylyl cyclase J Neurosci 2002 22 7931 7940 12223546
Lee, K. W. et al. Impaired D2 dopamine receptor function in mice lacking type 5 adenylyl cyclase. J Neurosci. 22, 7931–7940 (2002).12223546
82. Garcez PP Cenpj/CPAP regulates progenitor divisions and neuronal migration in the cerebral cortex downstream of Ascl1 Nat. Commun 2015 6 6474 10.1038/ncomms7474 25753651
Garcez, P. P. et al. Cenpj/CPAP regulates progenitor divisions and neuronal migration in the cerebral cortex downstream of Ascl1. Nat. Commun 6, 6474, 10.1038/ncomms7474 (2015).25753651 10.1038/ncomms7474
83. Rosso SB Sussman D Wynshaw-Boris A Salinas PC Wnt signaling through Dishevelled, Rac and JNK regulates dendritic development Nat. Neurosci. 2005 8 34 15608632
Rosso, S. B., Sussman, D., Wynshaw-Boris, A. & Salinas, P. C. Wnt signaling through Dishevelled, Rac and JNK regulates dendritic development. Nat. Neurosci. 8, 34 (2005).15608632
84. Rosso SB Inestrosa NC WNT signaling in neuronal maturation and synaptogenesis Front. Cell Neurosci 2013 7 103 23847469
Rosso, S. B. & Inestrosa, N. C. WNT signaling in neuronal maturation and synaptogenesis. Front. Cell Neurosci. 7, 103 (2013).23847469
85. Mahar I Subchronic peripheral neuregulin-1 increases ventral hippocampal neurogenesis and induces antidepressant-like effects PLoS One 2011 6 e26610 10.1371/journal.pone.0026610 22028923
Mahar, I. et al. Subchronic peripheral neuregulin-1 increases ventral hippocampal neurogenesis and induces antidepressant-like effects. PLoS One 6, e26610, 10.1371/journal.pone.0026610 (2011).22028923 10.1371/journal.pone.0026610
86. Mahar I Effects of neuregulin-1 administration on neurogenesis in the adult mouse hippocampus, and characterization of immature neurons along the septotemporal axis Sci. Rep. 2016 6 30467 10.1038/srep30467 27469430
Mahar, I. et al. Effects of neuregulin-1 administration on neurogenesis in the adult mouse hippocampus, and characterization of immature neurons along the septotemporal axis. Sci. Rep. 6, 30467, 10.1038/srep30467 (2016).27469430 10.1038/srep30467
87. Luo Z Regulator of G protein Signaling 6 (RGS6) suppresses late-age-onset Nigral Dopaminergic Neurodegeneration, α-Synuclein accumulation and motor dysfunction FASEB J. 2019 33 669.669 669.669
Luo, Z. et al. Regulator of G protein Signaling 6 (RGS6) suppresses late-age-onset Nigral Dopaminergic Neurodegeneration, α-Synuclein accumulation and motor dysfunction. FASEB J. 33, 669.669–669.669 (2019).
88. Bristow JM Reno TA Jo M Gonias SL Klemke RL Dynamic phosphorylation of tyrosine 665 in pseudopodium-enriched atypical kinase 1 (PEAK1) is essential for the regulation of cell migration and focal adhesion turnover J. Biol. Chem. 2013 288 123 131 23105102
Bristow, J. M., Reno, T. A., Jo, M., Gonias, S. L. & Klemke, R. L. Dynamic phosphorylation of tyrosine 665 in pseudopodium-enriched atypical kinase 1 (PEAK1) is essential for the regulation of cell migration and focal adhesion turnover. J. Biol. Chem. 288, 123–131 (2013).23105102
89. Reiner O Karzbrun E Kshirsagar A Kaibuchi K Regulation of neuronal migration, an emerging topic in autism spectrum disorders J. Neurochem. 2016 136 440 456 26485324
Reiner, O., Karzbrun, E., Kshirsagar, A. & Kaibuchi, K. Regulation of neuronal migration, an emerging topic in autism spectrum disorders. J. Neurochem. 136, 440–456 (2016).26485324
90. Wegiel J The neuropathology of autism: defects of neurogenesis and neuronal migration, and dysplastic changes Acta. Neuropathol. 2010 119 755 770 20198484
Wegiel, J. et al. The neuropathology of autism: defects of neurogenesis and neuronal migration, and dysplastic changes. Acta. Neuropathol. 119, 755–770 (2010).20198484
91. Goubert E Inhibition of the mitochondrial glutamate carrier SLC25A22 in astrocytes leads to intracellular glutamate accumulation Front Cell Neurosci 2017 11 149 10.3389/fncel.2017.00149 28620281
Goubert, E. et al. Inhibition of the mitochondrial glutamate carrier SLC25A22 in astrocytes leads to intracellular glutamate accumulation. Front Cell Neurosci. 11, 149, 10.3389/fncel.2017.00149 (2017).28620281 10.3389/fncel.2017.00149
92. Napolioni V Persico AM Porcelli V Palmieri L The mitochondrial aspartate/glutamate carrier AGC1 and calcium homeostasis: physiological links and abnormalities in autism Mol. Neurobiol 2011 44 83 92 10.1007/s12035-011-8192-2 21691713
Napolioni, V., Persico, A. M., Porcelli, V. & Palmieri, L. The mitochondrial aspartate/glutamate carrier AGC1 and calcium homeostasis: physiological links and abnormalities in autism. Mol. Neurobiol. 44, 83–92, 10.1007/s12035-011-8192-2 (2011).21691713 10.1007/s12035-011-8192-2
93. Wang Z Tian C Dhamala M Liu Z A small change in neuronal network topology can induce explosive synchronization transition and activity propagation in the entire network Sci. Rep. 2017 7 561 28373712
Wang, Z., Tian, C., Dhamala, M. & Liu, Z. A small change in neuronal network topology can induce explosive synchronization transition and activity propagation in the entire network. Sci. Rep. 7, 561 (2017).28373712
94. Petrelli F Pucci L Bezzi P Astrocytes and microglia and their potential link with autism spectrum disorders Front Cell Neurosci. 2016 10 21 26903806
Petrelli, F., Pucci, L. & Bezzi, P. Astrocytes and microglia and their potential link with autism spectrum disorders. Front Cell Neurosci. 10, 21 (2016).26903806
95. Shi W Reliability of whole-exome sequencing for assessing intratumor genetic heterogeneity Cell Rep. 2018 25 1446 1457 10.1016/j.celrep.2018.10.046 30404001
Shi, W. et al. Reliability of whole-exome sequencing for assessing intratumor genetic heterogeneity. Cell Rep. 25, 1446–1457, 10.1016/j.celrep.2018.10.046 (2018).30404001 10.1016/j.celrep.2018.10.046
96. Guo H Inherited and multiple de novo mutations in autism/developmental delay risk genes suggest a multifactorial model Mol. Autism 2018 9 64 10.1186/s13229-018-0247-z 30564305
Guo, H. et al. Inherited and multiple de novo mutations in autism/developmental delay risk genes suggest a multifactorial model. Mol. Autism 9, 64, 10.1186/s13229-018-0247-z (2018).30564305 10.1186/s13229-018-0247-z
97. Pinero J DisGeNET: a comprehensive platform integrating information on human disease-associated genes and variants Nucleic Acids Res 2017 45 D833 D839 10.1093/nar/gkw943 27924018
Pinero, J. et al. DisGeNET: a comprehensive platform integrating information on human disease-associated genes and variants. Nucleic Acids Res. 45, D833–D839, 10.1093/nar/gkw943 (2017).27924018 10.1093/nar/gkw943
98. The Gene Ontology, C. The Gene Ontology Resource: 20 years and still GOing strong Nucleic Acids Res. 2019 47 D330 D338 10.1093/nar/gky1055 30395331
The Gene Ontology, C. The Gene Ontology Resource: 20 years and still GOing strong. Nucleic Acids Res. 47, D330–D338, 10.1093/nar/gky1055 (2019).30395331 10.1093/nar/gky1055
99. Kohler S Expansion of the Human Phenotype Ontology (HPO) knowledge base and resources Nucleic Acids Res. 2019 47 D1018 D1027 10.1093/nar/gky1105 30476213
Kohler, S. et al. Expansion of the Human Phenotype Ontology (HPO) knowledge base and resources. Nucleic Acids Res. 47, D1018–D1027, 10.1093/nar/gky1105 (2019).30476213 10.1093/nar/gky1105
