Abstract
The exome sequences of approximately 8,000 children with autism spectrum disorder (ASD) and/or attention deficit hyperactivity disorder (ADHD) and 5,000 controls were analyzed, finding that individuals with ASD and individuals with ADHD had a similar burden of rare protein-truncating variants in evolutionarily constrained genes, both significantly higher than controls. This motivated a combined analysis across ASD and ADHD, identifying microtubule-associated protein 1A (MAP1A) as a new exome-wide significant gene conferring risk for childhood psychiatric disorders.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 / 30 days
cancel any time
Subscribe to this journal
Receive 12 print issues and online access
$209.00 per year
only $17.42 per issue
Buy this article
Purchase on Springer Link
Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
Data availability
Supplementary data are available as supplementary files to this manuscript (see Supplementary Tables 1, 3 and 5) or at the iPSYCH download page: http://ipsych.au.dk/downloads. For inquiries about more detailed data, contact iPSYCH lead investigator A.D.B. (anders@biomed.au.dk).
Code availability
Hail (0.1) and R scripts used to handle and analyze these data are available upon reasonable request from F.K.S. (satterst@broadinstitute.org).
References
Hallmayer, J. et al. Genetic heritability and shared environmental factors among twin pairs with autism. Arch. Gen. Psychiatry 68, 1095–1102 (2011).
Gaugler, T. et al. Most genetic risk for autism resides with common variation. Nat. Genet. 46, 881–885 (2014).
Larsson, H., Chang, Z., D’Onofrio, B. M. & Lichtenstein, P. The heritability of clinically diagnosed attention deficit hyperactivity disorder across the lifespan. Psychol. Med. 44, 2223–2229 (2014).
Pedersen, C. B. et al. The iPSYCH2012 case–cohort sample: new directions for unravelling genetic and environmental architectures of severe mental disorders. Mol. Psychiatry 23, 6–14 (2018).
Hollegaard, M. V. et al. Robustness of genome-wide scanning using archived dried blood spot samples as a DNA source. BMC Genet. 12, 58 (2011).
Hollegaard, M. V. et al. Archived neonatal dried blood spot samples can be used for accurate whole genome and exome-targeted next-generation sequencing. Mol. Genet. Metab. 110, 65–72 (2013).
Grove, J. et al. Identification of common genetic risk variants for autism spectrum disorder. Nat. Genet. 51, 431–444 (2019).
Demontis, D. et al. Discovery of the first genome-wide significant risk loci for attention deficit/hyperactivity disorder. Nat. Genet. 51, 63–75 (2019).
Neale, B. M. et al. Patterns and rates of exonic de novo mutations in autism spectrum disorders. Nature 485, 242–245 (2012).
De Rubeis, S. et al. Synaptic, transcriptional and chromatin genes disrupted in autism. Nature 515, 209–215 (2014).
Iossifov, I. et al. The contribution of de novo coding mutations to autism spectrum disorder. Nature 515, 216–221 (2014).
Samocha, K. E. et al. A framework for the interpretation of de novo mutation in human disease. Nat. Genet. 46, 944–950 (2014).
Sanders, S. J. et al. Insights into autism spectrum disorder genomic architecture and biology from 71 risk loci. Neuron 87, 1215–1233 (2015).
Lek, M. et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature 536, 285–291 (2016).
Kosmicki, J. A. et al. Refining the role of de novo protein-truncating variants in neurodevelopmental disorders by using population reference samples. Nat. Genet. 49, 504–510 (2017).
Jacquemont, S. et al. A higher mutational burden in females supports a ‘female protective model’ in neurodevelopmental disorders. Am. J. Hum. Genet. 94, 415–425 (2014).
Samocha, K. E. et al. Regional missense constraint improves variant deleteriousness prediction. Preprint at bioRxiv https://doi.org/10.1101/148353 (2017).
Martin, H. C. et al. Quantifying the contribution of recessive coding variation to developmental disorders. Science 362, 1161–1164 (2018).
Neyman, J. & Scott, E. On the use of c(α) optimal tests of composite hypothesis. Bull. Int. Stat. Inst. 41, 477–497 (1966).
Neale, B. M. et al. Testing for an unusual distribution of rare variants. PLoS Genet. 7, e1001322 (2011).
Myers, R. A. et al. A population genetic approach to mapping neurological disorder genes using deep resequencing. PLoS Genet. 7, e1001318 (2011).
Deciphering Developmental Disorders Study. Prevalence and architecture of de novo mutations in developmental disorders. Nature 542, 433–438 (2017).
Slager, R. E., Newton, T. L., Vlangos, C. N., Finucane, B. & Elsea, S. H. Mutations in RAI1 associated with Smith–Magenis syndrome. Nat. Genet. 33, 466–468 (2003).
Martin, J. et al. Biological overlap of attention-deficit/hyperactivity disorder and autism spectrum disorder: evidence from copy number variants. J. Am. Acad. Child Adolesc. Psychiatry 53, 761–770 (2014).
Ronald, A., Simonoff, E., Kuntsi, J., Asherson, P. & Plomin, R. Evidence for overlapping genetic influences on autistic and ADHD behaviours in a community twin sample. J. Child Psychol. Psychiatry 49, 535–542 (2008).
Chien, Y. L. et al. ADHD-related symptoms and attention profiles in the unaffected siblings of probands with autism spectrum disorder: focus on the subtypes of autism and Asperger’s disorder. Mol. Autism 8, 37 (2017).
World Health Organization. The ICD-10 Classification of Mental and Behavioural Disorders: Diagnostic Criteria for Research. (World Health Organization, 1993).
Hannelius, U. et al. Phenylketonuria screening registry as a resource for population genetic studies. J. Med. Genet. 42, e60 (2005).
Hollegaard, M. V. et al. Whole genome amplification and genetic analysis after extraction of proteins from dried blood spots. Clin. Chem. 53, 1161–1162 (2007).
Poulsen, J. B. et al. High-quality exome sequencing of whole-genome amplified neonatal dried blood spot DNA. PLoS One 11, e0153253 (2016).
McKenna, A. et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303 (2010).
Kalia, S. S. et al. Recommendations for reporting of secondary findings in clinical exome and genome sequencing, 2016 update (ACMG SFv2.0): a policy statement of the American College of Medical Genetics and Genomics. Genet. Med. 19, 249–255 (2017).
McLaren, W. et al. The Ensembl variant effect predictor. Genome Biol. 17, 122 (2016).
Acknowledgements
The iPSYCH project is funded by the Lundbeck Foundation (grant nos. R102-A9118 and R155-2014-1724) and the universities and university hospitals of Aarhus and Copenhagen. The DNSB resource at the Statens Serum Institut was supported by the Novo Nordisk Foundation. Sequencing of iPSYCH samples was supported by grants from the Simons Foundation (grant no. SFARI 311789 to M.J.D.) and the Stanley Foundation. Other support for this study was received from the National Institute of Mental Health (grant nos. 5U01MH094432-02, 5U01MH111660-02 and U01MH100229 to M.J.D.). Computational resources for handling and statistical analysis of iPSYCH data on the GenomeDK and Computerome HPC facilities were provided by, respectively, the Centre for Integrative Sequencing, iSEQ, Aarhus University, Denmark (grant to A.D.B.) and iPSYCH.
Author information
Authors and Affiliations
Consortia
Contributions
F.K.S. performed the analysis. R.K.W., T.S., E.M.W., F.L., D.D., J.A.K., J.G., D.S.P. and J.B.M. contributed to the analysis. F.K.S., R.K.W., C.S., J.B.-G., M.B.-H., M.N., O.M., D.M.H., T.M.W., P.B.M., A.D.B. and the iPSYCH–Broad Consortium were involved in sample selection, handling, processing and quality control. M.N., O.M., E.B.R., D.M.H., T.M.W., P.B.M., B.M.N., A.D.B. and M.J.D. were the project core principal investigator group. M.J.D. directed the project. B.M.N. and A.D.B. contributed to project direction. F.K.S. and M.J.D. wrote the manuscript.
Corresponding authors
Ethics declarations
Competing interests
T.M.W. has acted as advisor and lecturer to the pharmaceutical company H. Lundbeck A/S. B.M.N. is a member of the scientific advisory board at Deep Genomics and is a consultant for Camp4 Therapeutics, Takeda Pharmaceuticals and Biogen Inc.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Integrated supplementary information
Supplementary Figure 1 Variant groups for constrained rare PTV rate comparison.
a) Mean rare PTVs/person in constrained genes only (pLI ≥ 0.9). b) Mean rare PTVs/person in unconstrained genes only (pLI < 0.9). c) Mean rare synonymous variants/person in constrained genes only (pLI ≥ 0.9). ‘Rare’ as described in the main text. *** denotes p (compared to controls) < 0.001; all p values by logistic regression (Methods). For exact p values, see Supplementary Table 1. Sample numbers are as follows: ASD, no ID: 3,091; ASD+ADHD, no ID: 684; ADHD, no ID: 3,206; ASD, ID: 871; ASD+ADHD, ID: 217; ADHD, ID: 271; Control: 5,002. ID = intellectual disability. Error bars are Poisson standard error.
Supplementary Figure 2 Proportion of individuals within each phenotype category with 0, 1, 2, 3, or 4 crPTVs.
Simulated p value is chi-square test against Poisson expectation using observed mean, based on 10,000 replicates. Note that only one crPTV is counted per person per gene (removes approx. 0.2% of variants). crPTV = constrained (pLI ≥ 0.9) rare protein-truncating variant; ‘rare’ as described in the main text. ID = intellectual disability.
Supplementary Figure 3 Per-person rates of rare missense variants across phenotypes, considering only variants with MPC ≥ 2.
MPC score is a measure of the deleteriousness of a missense variant based on a regional model of constraint15, with higher values indicating a greater degree of deleteriousness. P values are for comparison to controls by logistic regression. ‘Rare’ as described in the main text. ID = intellectual disability.
Supplementary Figure 4 Variant groups for rare missense variant rate comparison.
a) Mean rare missense variants/person, MPC ≥ 2 only. b) Mean rare missense variants/person, MPC < 2 only. c) Mean rare synonymous variants/person, all genes. ‘Rare’ as described in the main text. * denotes p (compared to controls) < 0.05; ** denotes p < 0.01; *** denotes p < 0.001; all p values by logistic regression. For exact p values, see Supplementary Table 1. Sample numbers are as follows: ASD, no ID: 3,091; ASD+ADHD, no ID: 684; ADHD, no ID: 3,206; ASD, ID: 871; ASD+ADHD, ID: 217; ADHD, ID: 271; Control: 5,002. ID = intellectual disability. Error bars are Poisson standard error.
Supplementary Figure 5 Rates of constrained rare synonymous variants, considering ASD cases (n = 2,430) and ADHD cases (n = 2,360) with only a single diagnosis.
a) Rates in all constrained (pLI ≥ 0.9) genes; b) rates in the 212 constrained genes with a published rare de novo PTV in ASD (‘ASD de novo genes’)14. ‘Single’ diagnosis refers to ASD and ADHD cases without comorbid ASD+ADHD or intellectual disability, and without diagnoses of schizophrenia, bipolar disorder, affective disorder, or anorexia. Error bars are Poisson standard error. In the accompanying box, OR and p values are for comparison to controls (n = 5,002) by logistic regression. OR = odds ratio. Range = OR +/- standard error. ‘Rare’ as described in the main text.
Supplementary Figure 6 Quantile-quantile plots for gene association analysis combining Danish data with gnomAD.
a) Protein-truncating variants, b) synonymous variants, and c) missense variants with MPC ≥ 2. Analyses were limited to genes with a case-control odds ratio greater than 1 (based on a two-tailed Fisher’s exact test of rare variant counts, with ‘rare’ as described in the main text). To avoid bias, genes were excluded from the PTV and missense analyses if they had a case-control synonymous odds ratio greater than 1 (excluded 1,615/17,903, or 9.0% of genes). Sample numbers are 8,340 cases (combining the case categories shown in Table 1) and 49,781 controls (combining the 5,002 Danish controls with 44,779 non-Finnish Europeans from the non-psychiatric exome subset of gnomAD). The numbers of genes plotted are a) 3,182, b) 1,615, and c) 957.
Supplementary information
Supplementary Information
Supplementary Figs. 1–6, Supplementary Tables 2, 4 and 6 and Supplementary Note.
Supplementary Table 1
Rates, odds ratios and P values displayed in manuscript.
Supplementary Table 3
Lists of genes used in c-alpha tests, as well as test P values.
Supplementary Table 5
Full results of gene discovery analysis.
Rights and permissions
About this article
Cite this article
Satterstrom, F.K., Walters, R.K., Singh, T. et al. Autism spectrum disorder and attention deficit hyperactivity disorder have a similar burden of rare protein-truncating variants. Nat Neurosci 22, 1961–1965 (2019). https://doi.org/10.1038/s41593-019-0527-8
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41593-019-0527-8
This article is cited by
-
Impairments of cerebellar structure and function in a zebrafish KO of neuropsychiatric risk gene znf536
Translational Psychiatry (2024)
-
Connecting clinical and genetic heterogeneity in ADHD
Nature Genetics (2024)
-
Attention-deficit/hyperactivity disorder
Nature Reviews Disease Primers (2024)
-
Polygenic profiles define aspects of clinical heterogeneity in attention deficit hyperactivity disorder
Nature Genetics (2024)
-
Analysis of structural variation among inbred mouse strains
BMC Genomics (2023)