- Open Access
Targeted massively parallel sequencing of autism spectrum disorder-associated genes in a case control cohort reveals rare loss-of-function risk variants
Molecular Autism volume 6, Article number: 43 (2015)
Autism spectrum disorder (ASD) is highly heritable, yet genome-wide association studies (GWAS), copy number variation screens, and candidate gene association studies have found no single factor accounting for a large percentage of genetic risk. ASD trio exome sequencing studies have revealed genes with recurrent de novo loss-of-function variants as strong risk factors, but there are relatively few recurrently affected genes while as many as 1000 genes are predicted to play a role. As such, it is critical to identify the remaining rare and low-frequency variants contributing to ASD.
We have utilized an approach of prioritization of genes by GWAS and follow-up with massively parallel sequencing in a case-control cohort. Using a previously reported ASD noise reduction GWAS analyses, we prioritized 837 RefSeq genes for custom targeting and sequencing. We sequenced the coding regions of those genes in 2071 ASD cases and 904 controls of European white ancestry. We applied comprehensive annotation to identify single variants which could confer ASD risk and also gene-based association analysis to identify sets of rare variants associated with ASD.
We identified a significant over-representation of rare loss-of-function variants in genes previously associated with ASD, including a de novo premature stop variant in the well-established ASD candidate gene RBFOX1. Furthermore, ASD cases were more likely to have two damaging missense variants in candidate genes than controls. Finally, gene-based rare variant association implicates genes functioning in excitatory neurotransmission and neurite outgrowth and guidance pathways including CACNAD2, KCNH7, and NRXN1.
We find suggestive evidence that rare variants in synaptic genes are associated with ASD and that loss-of-function mutations in ASD candidate genes are a major risk factor, and we implicate damaging mutations in glutamate signaling receptors and neuronal adhesion and guidance molecules. Furthermore, the role of de novo mutations in ASD remains to be fully investigated as we identified the first reported protein-truncating variant in RBFOX1 in ASD. Overall, this work, combined with others in the field, suggests a convergence of genes and molecular pathways underlying ASD etiology.
Autism spectrum disorder (ASD) is a set of neurodevelopmental conditions diagnosed on the basis of a triad of symptoms: marked qualitative differences in social interaction, delayed or absent communication, and restricted interests or repetitive behaviors . ASD is highly prevalent, occurring with a prevalence of approximately 1 in 68 8-year-olds in the United States , and confers significant costs to individuals and families. While the underlying physiological causes are largely unknown, twin concordance and family-based studies have shown ASD to be highly heritable, suggesting a strong underlying genetic etiology [3–5]. Therefore, candidate gene and common risk allele discovery has been a topic of much research, and more than 100 genes and 50 genomic loci have been implicated in autism by linkage analyses, copy number variation screens, and genome-wide association studies . These candidates include single genes with mutations resulting in syndromes with ASD phenotypes [7–9], large cytogenetic aberrations , microduplications and microdeletions [11–15], and associated common alleles [16–20].
Overall, however, there has been limited replication across studies, demonstrating lack of support for the hypothesis that ASD is explained by common variants with strong to moderate effects. Instead, the common disease-rare variant (CDRV) hypothesis might better describe the underlying genetic architecture of ASD [21–23]. Identification of such rare genetic variants responsible for ASD risk has only recently been made possible with the application of massively parallel sequencing (MPS) technologies such as whole-genome sequencing (WGS), whole-exome sequencing (WES), and targeted region re-sequencing . These techniques have been successful in the identification of the genetic causes of dozens of Mendelian disorders . In ASD, which is phenotypically and genetically heterogeneous, WES studies thus far have focused primarily on the identification of de novo variation in family-based cohorts [26–32]. While this method has identified several potential ASD risk genes and mutations, each variant occurs at a low frequency in the population and the studies do not replicate well.
In contrast to this WES in trios approach, our strategy to detect rare functional variants in ASD is to apply MPS to genes within regions that have been implicated by common variant GWAS analysis in a large case-control cohort. We prioritized regions of association with GWAS-noise reduction (GWAS-NR) analyses of ASD datasets [16, 17] for follow-up sequencing, since disease-causing variants are likely not the common SNPs genotyped in GWAS. Sequencing of all exons of 837 RefSeq genes, in 2071 cases and 904 controls, revealed a significantly increased rate of stop gain/loss and splice altering loss-of-function (LOF) mutations, in subsets of candidate genes. Among the individual, LOF variants is a novel de novo premature stop gain in RBFOX1 in an ASD case. Our unique strategy supports the hypothesis that LOF variants in previously implicated ASD genes, including RBFOX1, and novel genes underlying the complex genetic architecture of ASD.
Individuals participating in this study have been collected over the course of several years, and ethical protocols with appropriate amendments have been approved. Individuals have been ascertained at the John P. Hussman Institute for Human Genomics (HIHG) at the University of Miami, Miller School of Medicine (Miami, FL), the University of South Carolina (Columbia, SC), and the Center for Human Genetics Research at Vanderbilt University (Nashville, TN). All participants were ascertained using protocols approved by the appropriate Institutional Review Boards, and the entire study falls under the University of Miami (UM) Institutional Review Board (IRB). This study was approved by the UM Medical Sciences IRB Committee B members: Ofelia Alvarez, M.D., Abdul Mian, Ph.D., Jose Castro, M.D., Jean Raymond Dauphin, O.D., Jean Jose, D.O., Norman Klein, J.D., Howard Landy, M.D., Stephen Richman, M.D., Eric Zetka, Pharm.D., and Liza Gordillo, M.S. In addition, samples from the Autism Genetic Resource Exchange (AGRE) and Simons Simplex Collection (SSC) were utilized. The use of AGRE and SSC samples is covered under the UM-IRB approved research protocol.
A total of 2399 unrelated individuals affected with ASD were used in this study. Nine hundred fifty-six individuals were recruited using Institutional Review Board approved protocols through the John P. Hussman Institute for Human Genomics (HIHG) (912 individuals) at the University of Miami, Miller School of Medicine (Miami, FL) and the Center for Human Genetics Research at Vanderbilt University (43 individuals) (Nashville, TN). Criteria for inclusion for ASD-affected individuals were as follows: 1—between 3 and 21 years of age, 2—an ASD diagnosis using DSM-III-R or DSM-IV criteria  supported by the Autism Diagnostic Interview (ADI-R) , and 3—an IQ equivalent >35 or developmental level >18 months as determined by the Vineland Adaptive Behavior Scale (VABS) . DNA was isolated from these individuals from whole blood collected via venipuncture.
Other ASD case DNA samples were included from the Autism Genetic Resource Exchange (14 individuals)  and the Simons Simplex Collection (1429 individuals) .
In addition, 1325 non-autistic control samples were acquired across two sites for this study. First, 525 individuals aged 4 to 21 years old were sampled through the HIHG. Participants were screened for eligibility by a questionnaire to determine whether the individual had been diagnosed with or had a parent or sibling with a developmental, behavioral, neurological, or other disability or physical conditions. If none of those conditions were present, parents of young children or the participants were informed and signed the informed consent and completed the Social Communication Questionnaire [38, 39] to screen for potential ASD. Three hundred twenty-six controls were part of a preterm birth study from cord blood collected at the Centennial Medical Center (Nashville, TN) from women aged from 18 to 40 years old with term pregnancies (>37 weeks gestation) and live, singleton births. The remaining 474 controls were obtained through the BioVu DNA repository at Vanderbilt .
Targeted sequencing probe-set design
Regions to be targeted for sequencing in this study were determined using criteria previously described . Briefly, the GWAS-NR algorithm was applied to two autism datasets and the top 5000 markers assorted into 2680 haplotype blocks  and 141 markers which were not in any block. Following statistical analysis using the truncated product method (TPM) of p values and a Monte Carlo simulation, 1535 genetic loci met a threshold of p ≤ 0.05. These regions were selected for follow-up sequencing in this study.
For target enrichment, we designed an Agilent SureSelect probe library using eArray web-based software (https://earray.chem.agilent.com/earray/) by bait-tiling 120-bp probes overlapping 60 bp. The targeted regions consisted of all exons of 837 RefSeq genes overlapping with blocks with p ≤ 0.05 TPM or nearest to significant intergenic blocks (Additional file 1). Nearest genes were chosen since GWAS-NR, by design, captures association from regions that may not be in strict linkage disequilibrium with a given SNP .
Massively parallel sequencing of targeted regions
Samples were prepared following Agilent and Illumina protocols for in-solution SureSelect enrichment and sequencing. Briefly, 3 μg of genomic DNA, quantified using the Broad Range dsDNA Assay (Life Technologies, Grand Island, NY), is fragmented in a 96-microTUBE plate on the Covaris E210 (Covaris, Woburn, MA), followed by the verification of the peak size distribution ranging from 150 to 200 nucleotides using the LabChip GX (PerkinElmer, Waltham, MA).
The sheared DNA is then prepared for Illumina sequencing in 96-well plates using the Sciclone G3 NGS Liquid Handling Workstation (PerkinElmer, Waltham, MA). Capture hybridization begins with 500 ng of prepped sampled library, mixed with blockers, buffer, and biotinylated RNA probes during a 24-h incubation at 65 °C. Streptavidin beads bind the biotinylated RNA probes to extract the desired captured library. The bead-bound hybridized product is amplified and purified in 96-well plates using a Zephyr Compact Liquid Handling Workstation (PerkinElmer, Waltham, MA). The final product is then verified for a single peak size approximately 300 to 325 bp using the LabChip GX prior to sequencing.
All samples were indexed and multiplexed to run 7–11 samples per lane on the Illumina HiSeq 2000 (Illumina, San Diego, CA) for paired-end 2 × 100 bp sequencing. Raw sequence data was processed using the Illumina Run Time Analysis base calling pipeline, initially with v1.7 and later with v1.8.
Subsequent processing was accomplished with an in-house sequencing data pipeline. This process includes alignment to the hg19 reference genome with the Burrows-Wheeler Aligner (BWA) , quality control with PICARD, and genotype calling performed with the Genome Analysis Toolkit (GATK) .
Sequencing quality control consisted of base quality score recalibration across flow cells and instruments as well as removal of duplicate reads resulting from PCR library amplification bias. For high-confidence genotype calling, we selected variants with a read depth greater than 8X, a call rate of at least 90 % across all samples, a log odds ratio under the trained Gaussian mixture model (VQSLOD) ≥ −3 and normalized, Phred-scaled likelihood of a reference genotype (PL) ≥200. The resulting alterations were annotated with SeattleSeq (http://snp.gs.washington.edu/SeattleSeqAnnotation137/index.jsp) and Polymorphism Phenotyping v2 (PolyPhen-2) (http://genetics.bwh.harvard.edu/pph2) . LOF variants were defined as stop gains, stop losses, or splice site alterations .
Sample quality control
From this total group of 3724 individuals, we excluded 6 samples with close relatedness or Mendelian inconsistencies identified by previously existing GWAS family data , 15 with underperforming sequencing metrics (less than 65 % of targeted bases covered at least 8X or call rates less than 90 % for all bases), and 8 with low (≤90 %) concordance with GWAS.
Racial and ethnic differences and population substructure have been demonstrated to have great impact on rare variant association analyses ; therefore, we set out to homogenize our sample cohort of 3695 using available Illumina 1M and 1M-duo genotyping array data . We performed principal component analysis (PCA) using the “smartpca” script in EIGENSTRAT . When available, we ran PCA on parents’ whole-genome genotyping data or, when parents were not available, on the individuals own. For those without genotyping data, we ran PCA with the sequencing data generated in this study. All participants who had both parents, or themselves when parental data was unavailable, fall within 2 standard deviations (SD) of the mean PC1 and PC2 values among our largest group of samples, clustering near the HapMap CEU population, were selected for further analysis. Additionally, all individuals with sequencing PCs within 2 SD of the mean PC1 and PC2 values of this genetically defined group defined above were also included. This analysis left 2975 European white individuals (2071 cases and 904 controls) for further analysis.
Variant burden and association analyses
Comparison of the mean number of total, missense, and nonsense single nucleotide variants (SNVs) between cases and controls was done by unpaired t tests. Differential burden analysis of the number cases versus controls harboring loss-of-function variants was calculated by the Fisher exact test. We performed differential burden testing of loss-of-function genes across all genes and then restricted this to a list of strong ASD candidate genes using the curated SFARI gene database . There were 45 genes with a gene score from SFARI that were included as the ASD candidate list.
We performed gene-based analyses testing association between targeted genes and ASD using the Optimized Sequence Kernel Association Test (SKAT-O) [49, 50] with default weights, adjusting for the first three PCs as covariates and correcting for gene size and linkage disequilibrium of variants within a gene. SNVs with differential case-control missingness p < 1.0 × 10−6 were excluded from analyses. This test was chosen because it is more powerful than collapsing and pooling methods in the presence of both risk and protective variants [50–53]. Analyses were run on several subsets of the variants including all and rare (minor allele frequency (MAF) ≤0.01) exonic, nonsynonymous, missense, predicted damaging, and LOF SNVs.
A subset of variants was analyzed by traditional capillary sequencing. Variant specific oligos were designed using the Primer3 (v. 0.4.0) (http://fokker.wi.mit.edu/primer3/input.htm)  and the UCSC genome browser (GRCh37/hg19). Sequencing reactions were performed with the Big Dye Terminator v3.1, sequenced on the Applied Biosystems 3730xl DNA Analyzer (Life Technologies, Grand Island, NY) and evaluated in the Sequencher software v4.10.1 (Gene Codes Corporation, Ann Arbor, MI).
Sequencing output and quality control
Targeted MPS was performed to identify potential ASD risk alleles in 2071 ASD cases by comparing to 904 non-autistic controls in GWAS-NR prioritized regions. We generated 30.3 ± 9.7 million passing filter Illumina HiSeq 2 × 100 bp reads per individual sample with only 8.6 ± 5.9 % duplicate reads removed to avoid amplification artifacts. Of the remaining reads, 96.8 ± 1.4 % aligned to the human genome and 81.2 ± 10.2 % aligned to the 29.1 Mb “near target” region, consisting of bases covered by targeting probes and their 200-bp flanks. This highly efficient targeting, sequencing, and alignment lead to average depth coverage of 78.1 ± 25.9 times of the near target bases and 87.6 ± 5.6 % of targeted bases covered at least 10X.
Genotype calling and variant discovery
Overall, we identified 545,916 genomic positions at which at least one of the 2911 individuals analyzed had a non-reference genotype call, meeting our genotype calling quality criteria. Of those, 231,945 SNVs (56.2 %) were not previously annotated in dbSNP137. Genotype call quality and sample matching integrity was determined by comparison of SNP genotyping data available for most samples [17, 18]. The SNV calls between targeted sequencing and genotyping arrays at over 30,000 markers had an average concordance of 99.3 ± 0.01 %.
Functional annotation of variants
We focused our analyses on SNVs affecting only coding exons of the targeted genes to best identify functional variation. We identified 25,966 such SNVs in our cohort including 16,330 (62.9 %) SNVs were not previously reported in dbSNP134. We employed publically available databases to annotate each variant, including PolyPhen2 damaging predictions. When multiple isoforms were identified, the most deleterious SNV was included. To allow for comparison between cases and controls, we normalized the number of variants correcting for differences in sequence coverage and missing data. This normalization procedure gave extremely consistent rates of variation in cases and controls. We identified no difference between cases and controls in the overall number of SNVs nor in subsets of missense, predicted damaging, or nonsense and splice site variants (Table 1).
Gene-based rare variant association
Because we did not identify a global difference in the number of SNVs between cases and controls, we tested the hypothesis that sets of protein-coding SNVs in genes would be associated with ASD using gene-based rare variant set tests with SKAT-O. We restricted our analyses to rare SNVs (MAF ≤0.01) and identified nominal association (p ≤ 0.05) with 16, 19, and 25 genes when examining exonic, nonsynonymous, and damaging SNVs, respectively, (Table 2). However, no gene was significant in any SKAT-O analysis after correction for multiple gene testing (corrected p = 9.5 × 10−5).
Individual variant prioritization
Since the autism phenotype is highly variable and the numbers of candidate genes are many, rare SNVs found in only one or a few individuals may be clinically relevant by leading to an increased risk or causation of the disorder. We employed a two-tier strategy to identify individual SNVs or genes contributing to autism etiology. First, we identified genes with homozygous or compound heterozygous rare (MAF ≤0.01 for each variant allele), damaging (as predicted by SIFT or PolyPhen2) SNVs. Among all individuals, a total of 47 genes were affected by at least two rare damaging variants in the same individual (Additional file 2). Among the genes with two hits uniquely in cases were three candidate genes linked to autism by a literature-curated database, DMD, SYNE1, and TBL1X and two genes implicated in schizophrenia (GRIK4) and intellectual disability (PQBP1). Overall, there are significantly more ASD cases (144/2071) than controls (26/904) carrying genes with two rare predicted damaging SNVs (Fisher’s exact p = 0.0001).
The second strategy identified rare SNVs (MAF ≤ 0.01) causing LOF. These alterations cause nonsense premature stops or stop losses and splice site changes. Among our 2975 individuals, we identified 464 rare LOF variants (Additional file 3). Across all these variants, there was no overall enrichment of LOF in ASD cases (Fisher’s exact p = 0.325). However, when restricting the gene list to only those that have been implicated previously in ASD, there is a significant enrichment of loss-of-function variants in cases with 27 cases carrying such variants compared to only four controls (Fisher’s exact p = 0.032).
Inheritance of loss-of-function SNVs
Since de novo variations in protein-coding genes have been recently shown to be a strong risk factor for autism, we determined the inheritance status of the LOF SNVs in seven ASD candidate genes and an additional ten LOF SNVs in genes with putative neuronal function (Table 3). We performed standard capillary sequencing of the individual in which the SNV was identified along with both parents and available siblings, if any. All 17 variants were validated. Of the 17 variants examined, one, a premature stop codon in RBFOX1, was de novo (Fig. 1), 13 were maternally inherited, and three were paternally inherited. Two of the variants, premature stop alterations GPR110 and DOCK1, were identified in an individual from a family with multiple affected individuals. The variant in GPR110 did not segregate with the ASD, but the LOF SNV in DOCK1 is identified in the affected pro-band and ASD-affected mother and not an unaffected sibling.
ASD has few common variants known to impact genetic risk leading to the hypothesis that rare variants with larger effects in candidate genes are contributing to the phenotype. Since the objective of GWAS-NR is to identify regions of association for follow-up analysis , our study set out to identify such rare variants specifically in previously associated candidate genes. Similar approaches have been successful at identifying the contribution of rare variants in other complex diseases such as type II diabetes , hypertriglyceridemia , and inflammatory bowel disease [57, 58].
The majority of sequencing studies in ASD thus far have been whole-exome sequencing in large trio-based cohorts [26–32]. These efforts are designed specifically to identify de novo variants of high effect and have identified several recurrent de novo loss-of-function changes as strong risk factors for ASD. While the current study was not designed for identification of de novo variants, we undertook a two-pronged approach to identify individual loss-of-function variants that could have major effects. Firstly, we identified a statistically significant increase in the number of genes with at least two rare predicted damaging alterations in the same individual in ASD cases compared to controls (p = 0.0001). Secondly, we found an excess of LOF mutations in ASD candidate genes in our cases (p = 0.032). Overall, this suggests that carrying two damaging hits or LOF in previously identified candidate genes is a risk factor for ASD. This is similar to previous WES studies specifically identifying two-hit and LOF mutations as risk factors in ASD [30, 59].
Among these loss-of-functions, we identified a de novo premature stop variation in the splice regulation factor RBFOX1. This is an excellent candidate gene as several genomic lesions including microdeletions [11, 60, 61] and rare missense mutations  have been reported in ASD cases in RBFOX1. Additionally, its role as a splice regulatory factor  was highlighted in a transcriptome study where aberrant splicing of RBFOX1-dependent alternative exons in the brains of ASD patients was identified . Targets of RBFOX1 splice regulation are enriched for genes involved in neuronal excitation  and cytoskeletal reorganization . This study offers the first evidence of a de novo or LOF sequence variant in RBOFX1 present in an ASD patient.
While the presence of a premature stop or splice altering variants does not necessarily indicate a molecular loss-of-function, the overall analysis points to genes involved in several underlying physiological mechanisms that support a role for excitatory neurotransmission and altered neurite outgrowth and guidance in ASD. In the analysis of genes with two damaging mutations unique to ASD cases, we identified two genes involved in the glutamatergic signaling pathway, the metabotropic glutamate receptor 3, GRM3, and glutamate receptor, ionotropic, kainate 4, GRIK4. Common variants in GRM3 have been shown to have association with ASD  while rare CNVs in GRIK4 have been implicated , and both genes have evidence for association in bipolar disorder  and schizophrenia [68, 69]. IL1RAPL2 is a gene related to ASD candidate gene IL1RAPL1, which regulates the formation and stabilization of glutamatergic synapses , and rare mutations have been implicated in ASD and neurodevelopmental phenotypes [71, 72]. Neuronal adhesion and guidance molecules are also implicated. SEMA6A has not been directly implicated in ASD or other neuropsychiatric disorders, but mouse mutants demonstrate anatomical abnormalities in limbic and cortical cellular organization  and axon guidance . Finally, polyglutamine-binding protein 1 (PQBP1) has been implicated as a splice factor necessary for proper neurite outgrowth , and rare mutations in PQBP1 have also been implicated in X-linked mental retardation .
Likewise, a commonality of function can be identified among the loss-of-function SNVs. . There are 199 genes with LOF variants uniquely in cases, including 17 in genes previously linked to ASD. Among the most intriguing candidates are two stop inducing SNVs in the cell polarity regulator FAT1 and two stop and one splice SNV in the nuclear envelope protein SYNE1 which have previously been implicated by de novo missense variants in ASD cases [26, 29], though the variants in this study were maternally inherited. The glutamate receptor GRIK2 has been associated with ASD and neurodevelopment several studies [11, 66, 77, 78], and we identify a premature stop and two splice change SNVs. We detected two stop and one splice SNV in the well-recognized ASD candidate gene NRXN1 [79–82].
Taking into account all variants in this study, no differences in global burden of variation were identified between cases and controls reflecting the findings of whole-exome studies in which the rate of mutations is similar in cases compared to unaffected controls or siblings [26, 28, 29, 83]. Additionally, no individual coding variants were significantly associated with ASD. As such, we applied rare variant association testing using SKAT-O but did not identify any gene with a statistically significant association of rare variants. Despite having one of the largest reported sample sizes for an ASD case-control sequencing study, it has been suggested that even larger sample sizes on the order of tens of thousands of cases will be required to identify associations of rare variants with complex disease . However, investigation of the nominally significant genes can be informative for potential roles of several genes in ASD risk.
For example, when observing genes with nominal association between ASD and sets of rare exonic variants, the most significant gene is CACNA2D1 (Calcium Channel, Voltage-Dependent, Alpha 2/Delta Subunit) (p = 0.00047). CACNA2D1 regulates influx of calcium ions into the cell upon membrane polarization and Mutations in a related gene (CACNA1C) are known to cause Timothy syndrome which includes neurological and developmental deficits and often autism symptoms . In addition, rare exonic and nonsynonymous variants in the potassium channel KCNH7 showed nominal association (p = 0.043). Variants in this gene have been previously linked to bipolar disorder , schizophrenia , and developmental delay . Rare exonic variants in the well-established ASD and schizophrenia candidate gene NRXN1 [79–81] show nominal association (p = 0.041). Such pleiotropic relationships between variants conferring ASD risk and other neuropsychiatric disorders have been documented [89, 90]. Overall, the nominally significant gene sets point toward synaptic function as critical in risk to ASD, and an increase in sample size and collaboration across multiple datasets may increase our power to detect significant associations with these or other potential candidate genes.
In conclusion, previous studies of the genetic causes of ASD have implicated hundreds of potential loci with none fully explaining the extent of genetic risk. Therefore, it is a necessity to identify the many rare variants in numerous genes contributing to the disorder using existing GWAS data to prioritize regions of ASD-specific susceptibility variants and then find underlying rare risk variants using recently developed massively parallel sequencing technologies. This has the distinct advantage of reduced cost, ease of multiplexing, and specific information content as compared to whole-exome sequencing studies.
We have used this approach to arrive at three observations. First, our evidence is suggestive that accumulation of rare variants in synaptic genes, including CACNAD2, KCNH7, and NRXN1, is associated with ASD, but independent support from future sequencing studies will be necessary to statistically support this idea. Second, we identified an over-representation of two damaging hit across all genes and LOF mutations in ASD candidate genes as a risk factor for ASD and implicate damaging mutations in glutamate signaling receptors and neuronal adhesion and guidance molecules. Finally, we provide supporting evidence of de novo coding mutations and the role of RBFOX1 dysfunction as a potential risk factor for ASD. These observations highlight the heterogeneity of the genetic etiology of ASD but point toward a convergence of pathways and mechanisms which underlie the complex phenotype.
autism spectrum disorder
Genome Analysis Toolkit
genome-wide association study-noise reduction
optimized sequence kernel association test
single nucleotide variant
American Psychiatric Association. Diagnostic and statistical manual of mental disorders. 5th ed. Arlington, VA: American Psychiatric Publishing; 2013.
Developmental Disabilities Monitoring Network Surveillance Year 2010 Principal Investigators, Centers for Disease Control and Prevention (CDC). Prevalence of autism spectrum disorder among children aged 8 years—autism and developmental disabilities monitoring network, 11 sites, United States, 2010. MMWR Surveill Summ. 2010;2014(63):1–21.
Ritvo ER, Freeman BJ, Mason-Brothers A, Mo A, Ritvo AM. Concordance for the syndrome of autism in 40 pairs of afflicted twins. Am J Psychiatry. 1985;142:74–7.
Folstein S, Rutter M. Infantile autism: a genetic study of 21 twin pairs. J Child Psychol Psychiatry. 1977;18:297–321.
Bailey A, Le Couteur A, Gottesman I, Bolton P, Simonoff E, Yuzda E, et al. Autism as a strongly genetic disorder: evidence from a British twin study. Psychol Med. 1995;25:63–77.
Betancur C. Etiological heterogeneity in autism spectrum disorders: more than 100 genetic and genomic disorders and still counting. Brain Res. 2011;1380:42–77.
Wong V. Study of the relationship between tuberous sclerosis complex and autistic disorder. J Child Neurol. 2006;21:199–204.
Amir RE, Van den Veyver IB, Wan M, Tran CQ, Francke U, Zoghbi HY. Rett syndrome is caused by mutations in X-linked MECP2, encoding methyl-CpG-binding protein 2. Nat Genet. 1999;23:185–8.
Wilson HL, Wong AC, Shaw SR, Tse WY, Stapleton GA, Phelan MC, et al. Molecular characterisation of the 22q13 deletion syndrome supports the role of haploinsufficiency of SHANK3/PROSAP2 in the major neurological symptoms. J Med Genet. 2003;40:575–84.
Girirajan S, Brkanac Z, Coe BP, Baker C, Vives L, Vu TH, et al. Relative burden of large CNVs on a range of neurodevelopmental phenotypes. PLoS Genet. 2011;7(11), e1002334.
Griswold AJ, Ma D, Cukier HN, Nations LD, Schmidt MA, Chung RH, et al. Evaluation of copy number variations reveals novel candidate genes in autism spectrum disorder-associated pathways. Hum Mol Genet. 2012;21:3513–23.
Glessner JT, Wang K, Cai G, Korvatska O, Kim CE, Wood S, et al. Autism genome-wide copy number variation reveals ubiquitin and neuronal genes. Nature. 2009;459:569–73.
Pinto D, Pagnamenta AT, Klei L, Anney R, Merico D, Regan R, et al. Functional impact of global rare copy number variation in autism spectrum disorders. Nature. 2010;466:368–72.
Levy D, Ronemus M, Yamrom B, Lee YH, Leotta A, Kendall J, et al. Rare de novo and transmitted copy-number variation in autistic spectrum disorders. Neuron. 2011;70:886–97.
Pinto D, Delaby E, Merico D, Barbosa M, Merikangas A, Klei L, et al. Convergence of genes and cellular pathways dysregulated in autism spectrum disorders. Am J Hum Genet. 2014;94:677–94.
Hussman JP, Chung RH, Griswold AJ, Jaworkski JM, Salyakina D, Ma D, et al. A noise-reduction GWAS analysis implicates altered regulation of neurite outgrowth and guidance in autism. Mol Autism. 2011;2(1):1.
Ma DQ, Salyakina D, Jaworski JM, Konidari I, Whitehead PL, Andersen A, et al. A genome-wide association study of autism reveals a common novel risk locus at 5p14.1. Ann Hum Genet. 2009;73:263–73.
Wang K, Zhang H, Ma D, Bucan M, Glessner JT, Abrahams BS, et al. Common genetic variants on 5p14.1 associate with autism spectrum disorders. Nature. 2009;459:528–33.
Anney R, Klei L, Pinto D, Regan R, Conroy J, Magalhaes TR, et al. A genome-wide scan for common alleles affecting risk for autism. Hum Mol Genet. 2010;19:4072–82.
Weiss LA, Arking DE. Gene Discovery Project of Johns Hopkins & the Autism Consortium, Daly MJ, Chakravarti A. A genome-wide linkage and association scan reveals novel loci for autism. Nature. 2009;461:802–8.
Pritchard JK. Are rare variants responsible for susceptibility to complex diseases? Am J Hum Genet. 2001;69:124–37.
Pritchard JK, Cox NJ. The allelic architecture of human disease genes: common disease-common variant…or not? Hum Mol Genet. 2002;11:2417–23.
Reich DE, Lander ES. On the allelic spectrum of human disease. Trends Genet. 2001;17:502–10.
Metzker ML. Sequencing technologies—the next generation. Nat Rev Genet. 2010;11:31–46.
Bamshad MJ, Ng SB, Bigham AW, Tabor HK, Emond MJ, Nickerson DA, et al. Exome sequencing as a tool for Mendelian disease gene discovery. Nat Rev Genet. 2011;12:745–55.
O’Roak BJ, Deriziotis P, Lee C, Vives L, Schwartz JJ, Girirajan S, et al. Exome sequencing in sporadic autism spectrum disorders identifies severe de novo mutations. Nat Genet. 2011;43:585–9.
O’Roak BJ, Vives L, Girirajan S, Karakoc E, Krumm N, Coe BP, et al. Sporadic autism exomes reveal a highly interconnected protein network of de novo mutations. Nature. 2012;485:246–50.
Iossifov I, Ronemus M, Levy D, Wang Z, Hakker I, Rosenbaum J, et al. De novo gene disruptions in children on the autistic spectrum. Neuron. 2012;74:285–99.
Neale BM, Kou Y, Liu L, Ma’ayan A, Samocha KE, Sabo A, et al. Patterns and rates of exonic de novo mutations in autism spectrum disorders. Nature. 2012;485:242–5.
Lim ET, Raychaudhuri S, Sanders SJ, Stevens C, Sabo A, MacArthur DG, et al. Rare complete knockouts in humans: population distribution and significant role in autism spectrum disorders. Neuron. 2013;77:235–42.
De Rubeis S, He X, Goldberg AP, Poultney CS, Samocha K, Ercument Cicek A, et al. Synaptic, transcriptional and chromatin genes disrupted in autism. Nature. 2014;515:209–15.
Iossifov I, O’Roak BJ, Sanders SJ, Ronemus M, Krumm N, Levy D, et al. The contribution of de novo coding mutations to autism spectrum disorder. Nature. 2014;515:216–21.
American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders (DSM-IV), text revision. Washington, DC: American Psychiatric Press, Inc; 2000.
Rutter M, LeCouteur A, Lord C. Autism Diagnostic Interview, Revised (ADI-R). 2003.
Carter AS, Volkmar FR, Sparrow SS, Wang J, Lord C, Dawson G, et al. The Vineland Adaptive Behavior Scales: supplementary norms for individuals with Autism. J Autism Dev Disord. 1998;25:287–302.
Fischbach GD, Lord C. The Simons Simplex Collection: a resource for identification of autism genetic risk factors. Neuron. 2010;68:192–5.
Rutter M, Bailey A, Berument SK, LeCouteur A, Lord C, Pickles A. Social Communication Questionnaire (SCQ). Los Angeles, CA: Western Psychological Services; 2003.
Chandler S, Charman T, Baird G, Simonoff E, Loucas T, Meldrum D, et al. Validation of the social communication questionnaire in a population cohort of children with autism spectrum disorders. J Am Acad Child Adolesc Psychiatry. 2007;46:1324–32.
Denny JC, Ritchie MD, Basford MA, Pulley JM, Bastarache L, Brown-Gentry K, et al. PheWAS: demonstrating the feasibility of a phenome-wide scan to discover gene-disease associations. Bioinformatics. 2010;26:1205–10.
Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, et al. The structure of haplotype blocks in the human genome. Science. 2002;296:2225–9.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–60.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20:1297–303.
Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7:248–9.
MacArthur DG, Balasubramanian S, Frankish A, Huang N, Morris J, Walter K, et al. A systematic survey of loss-of-function variants in human protein-coding genes. Science. 2012;335:823–8.
O’Connor TD, Kiezun A, Bamshad M, Rich SS, Smith JD, Turner E, et al. Fine-scale patterns of population stratification confound rare variant association tests. PLoS One. 2013;8, e65834.
Patterson N, Price AL, Reich D. Population structure and eigenanalysis. PLoS Genet. 2006;2, e190.
Banerjee-Basu S, Packer A. SFARI gene: an evolving database for the autism research community. Dis Model Mech. 2010;3:133–5.
Wu MC, Lee S, Cai T, Li Y, Boehnke M, Lin X. Rare-variant association testing for sequencing data with the sequence kernel association test. Am J Hum Genet. 2011;89:82–93.
Lee S, Emond MJ, Bamshad MJ, Barnes KC, Rieder MJ, Nickerson DA, et al. Optimal unified approach for rare-variant association testing with application to small-sample case-control whole-exome sequencing studies. Am J Hum Genet. 2012;91:224–37.
Ladouceur M, Zheng HF, Greenwood CM, Richards JB. Empirical power of very rare variants for common traits and disease: results from Sanger sequencing 1998 individuals. Eur J Hum Genet. 2013;21(9):1027–30. 13(4).
Lee S, Wu MC, Lin X. Optimal tests for rare variant effects in sequencing association studies. Biostatistics. 2012;13(4):762–75.
Bodenhofer U, Hochreiter S. Associating complex traits with rare variants identified by NGS: improving power by a position-dependent kernel approach. 2012.
Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3—new capabilities and interfaces. Nucleic Acids Res. 2012;40, e115.
Nejentsev S, Walker N, Riches D, Egholm M, Todd JA. Rare variants of IFIH1, a gene implicated in antiviral responses, protect against type 1 diabetes. Science. 2009;324:387–9.
Johansen CT, Wang J, Lanktree MB, Cao H, McIntyre AD, Ban MR, et al. Excess of rare variants in genes identified by genome-wide association study of hypertriglyceridemia. Nat Genet. 2010;42:684–7.
Rivas MA, Beaudoin M, Gardet A, Stevens C, Sharma Y, Zhang CK, et al. Deep resequencing of GWAS loci identifies independent rare variants associated with inflammatory bowel disease. Nat Genet. 2011;43:1066–73.
Momozawa Y, Mni M, Nakamura K, Coppieters W, Almer S, Amininejad L, et al. Resequencing of positional candidates identifies low frequency IL23R coding variants protecting against inflammatory bowel disease. Nat Genet. 2011;43:43–7.
Kenny EM, Cormican P, Furlong S, Heron E, Kenny G, Fahey C, et al. Excess of rare novel loss-of-function variants in synaptic genes in schizophrenia and autism spectrum disorders. Mol Psychiatry. 2014;19(8):872–9.
Sebat J, Lakshmi B, Malhotra D, Troge J, Lese-Martin C, Walsh T, et al. Strong association of de novo copy number mutations with autism. Science. 2007;316:445–9.
Davis LK, Maltman N, Mosconi MW, Macmillan C, Schmitt L, Moore K, et al. Rare inherited A2BP1 deletion in a proband with autism and developmental hemiparesis. Am J Med Genet A. 2012;158A:1654–61.
Martin CL, Duvall JA, Ilkin Y, Simon JS, Arreaza MG, Wilkes K, et al. Cytogenetic and molecular characterization of A2BP1/FOX1 as a candidate gene for autism. Am J Med Genet B Neuropsychiatr Genet. 2007;144B:869–76.
Underwood JG, Boutz PL, Dougherty JD, Stoilov P, Black DL. Homologues of the Caenorhabditis elegans Fox-1 protein are neuronal splicing regulators in mammals. Mol Cell Biol. 2005;25:10005–16.
Voineagu I, Wang X, Johnston P, Lowe JK, Tian Y, Horvath S, et al. Transcriptomic analysis of autistic brain reveals convergent molecular pathology. Nature. 2011;474:380–4.
Gehman LT, Stoilov P, Maguire J, Damianov A, Lin CH, Shiue L, et al. The splicing regulator Rbfox1 (A2BP1) controls neuronal excitation in the mammalian brain. Nat Genet. 2011;43:706–11.
Casey JP, Magalhaes T, Conroy JM, Regan R, Shah N, Anney R, et al. A novel approach of homozygous haplotype sharing identifies candidate genes in autism spectrum disorder. Hum Genet. 2012;131:565–79.
Knight HM, Walker R, James R, Porteous DJ, Muir WJ, Blackwood DH, et al. GRIK4/KA1 protein expression in human brain and correlation with bipolar disorder risk variant status. Am J Med Genet B Neuropsychiatr Genet. 2012;159B:21–9.
Pickard BS, Knight HM, Hamilton RS, Soares DC, Walker R, Boyd JK, et al. A common variant in the 3′UTR of the GRIK4 glutamate receptor gene affects transcript abundance and protects against bipolar disorder. Proc Natl Acad Sci U S A. 2008;105:14940–5.
Harrison PJ, Lyon L, Sartorius LJ, Burnet PW, Lane TA. The group II metabotropic glutamate receptor 3 (mGluR3, mGlu3, GRM3): expression, function and involvement in schizophrenia. J Psychopharmacol. 2008;22:308–22.
Hayashi T, Yoshida T, Ra M, Taguchi R, Mishina M. IL1RAPL1 associated with mental retardation and autism regulates the formation and stabilization of glutamatergic synapses of cortical neurons through RhoA signaling pathway. PLoS One. 2013;8, e66254.
Piton A, Michaud JL, Peng H, Aradhya S, Gauthier J, Mottron L, et al. Mutations in the calcium-related gene IL1RAPL1 are associated with autism. Hum Mol Genet. 2008;17:3965–74.
Bhat SS, Ladd S, Grass F, Spence JE, Brasington CK, Simensen RJ, et al. Disruption of the IL1RAPL1 gene associated with a pericentromeric inversion of the X chromosome in a patient with mental retardation and autism. Clin Genet. 2008;73:94–6.
Runker AE, O’Tuathaigh C, Dunleavy M, Morris DW, Little GE, Corvin AP, et al. Mutation of Semaphorin-6A disrupts limbic and cortical connectivity and models neurodevelopmental psychopathology. PLoS One. 2011;6, e26488.
Runker AE, Little GE, Suto F, Fujisawa H, Mitchell KJ. Semaphorin-6A controls guidance of corticospinal tract axons at multiple choice points. Neural Dev. 2008;3:34–8. 104-3-34.
Wang Q, Moore MJ, Adelmant G, Marto JA, Silver PA. PQBP1, a factor linked to intellectual disability, affects alternative splicing associated with neurite outgrowth. Genes Dev. 2013;27:615–26.
Kalscheuer VM, Freude K, Musante L, Jensen LR, Yntema HG, Gecz J, et al. Mutations in the polyglutamine binding protein 1 gene cause X-linked mental retardation. Nat Genet. 2003;35:313–5.
Jamain S, Betancur C, Quach H, Philippe A, Fellous M, Giros B, et al. Linkage and association of the glutamate receptor 6 gene with autism. Mol Psychiatry. 2002;7:302–10.
Shaltiel G, Maeng S, Malkesman O, Pearson B, Schloesser RJ, Tragon T, et al. Evidence for the involvement of the kainate receptor subunit GluR6 (GRIK2) in mediating behavioral displays related to behavioral symptoms of mania. Mol Psychiatry. 2008;13:858–72.
Kim HG, Kishikawa S, Higgins AW, Seong IS, Donovan DJ, Shen Y, et al. Disruption of neurexin 1 associated with autism spectrum disorder. Am J Hum Genet. 2008;82:199–207.
Gauthier J, Siddiqui TJ, Huashan P, Yokomaku D, Hamdan FF, Champagne N, et al. Truncating mutations in NRXN2 and NRXN1 in autism spectrum disorders and schizophrenia. Hum Genet. 2011;130:563–73.
Duong L, Klitten LL, Moller RS, Ingason A, Jakobsen KD, Skjodt C, et al. Mutations in NRXN1 in a family multiply affected with brain disorders: NRXN1 mutations and brain disorders. Am J Med Genet B Neuropsychiatry Genet. 2012;159B(3):354–8.
Ching MS, Shen Y, Tan WH, Jeste SS, Morrow EM, Chen X, et al. Deletions of NRXN1 (neurexin-1) predispose to a wide spectrum of developmental disorders. Am J Med Genet B Neuropsychiatr Genet. 2010;153B:937–47.
Sanders SJ, Murtha MT, Gupta AR, Murdoch JD, Raubeson MJ, Willsey AJ, et al. De novo mutations revealed by whole-exome sequencing are strongly associated with autism. Nature. 2012;485:237–41.
Zuk O, Schaffner SF, Samocha K, Do R, Hechter E, Kathiresan S, et al. Searching for missing heritability: designing rare variant association studies. Proc Natl Acad Sci U S A. 2014;111:E455–64.
Splawski I, Timothy KW, Sharpe LM, Decher N, Kumar P, Bloise R, et al. Ca(V)1.2 calcium channel dysfunction causes a multisystem disorder including arrhythmia and autism. Cell. 2004;119:19–31.
Strauss KA, Markx S, Georgi B, Paul SM, Jinks RN, Hoshi T, et al. A population-based study of KCNH7 p.Arg394His and bipolar spectrum disorder. Hum Mol Genet. 2014;23:6395–406.
Heide J, Mann SA, Vandenberg JI. The schizophrenia-associated Kv11.1-3.1 isoform results in reduced current accumulation during repetitive brief depolarizations. PLoS One. 2012;7, e45624.
Belengeanu V, Gamage TH, Farcas S, Stoian M, Andreescu N, Belengeanu A, et al. A de novo 2.3 Mb deletion in 2q24.2q24.3 in a 20-month-old developmentally delayed girl. Gene. 2014;539:168–72.
Carroll LS, Owen MJ. Genetic overlap between autism, schizophrenia and bipolar disorder. Genome Med. 2009;1(10):102.
Cukier HN, Dueker ND, Slifer SH, Lee JM, Whitehead PL, Lalanne E, et al. Exome sequencing of extended families with autism reveals genes shared across neurodevelopmental and neuropsychiatric disorders. Mol Autism. 2014;5:1-2392-5-1.
We thank all the participants and their family members in this study and personnel at the John P. Hussman Institute for Human Genomics (HIHG) including staff at the HIHG Biorepository, autism ascertainment team, and the HIHG Center for Genome Technology. We acknowledge the contribution of ASD samples from Harry Wright and Ruth Abramson from the University of South Carolina approved by the University of South Carolina Human Subjects Board and the contribution of preterm birth controls from Scott Williams and Ramkumar Menon approved by the Vanderbilt Human Research Protection Program. A subset of the participants was ascertained while Dr. Pericak-Vance was a faculty member at Duke University. We are grateful to all of the families at the participating Simons Simplex Collection (SSC) sites, as well as the principal investigators (A. Beaudet, R. Bernier, J. Constantino, E. Cook, E. Fombonne, D. Geschwind, R. Goin-Kochel, E. Hanson, D. Grice, A. Klin, D. Ledbetter, C. Lord, C. Martin, D. Martin, R. Maxim, J. Miles, O. Ousley, K. Pelphrey, B. Peterson, J. Piggot, C. Saulnier, M. State, W. Stone, J. Sutcliffe, C. Walsh, Z. Warren, E. Wijsman). We gratefully acknowledge the resources provided by the Autism Genetic Resource Exchange (AGRE) Consortium and the participating AGRE families. The Autism Genetic Resource Exchange is a program of Autism Speaks and is supported, in part, by grant 1U24MH081810 from the National Institute of Mental Health to Clara M. Lajonchere (PI). This work was supported by the National Institutes of Health (R01MH080647, P01NS026630) and a gift from the John P. Hussman Foundation.
The authors declare that they have no competing interests.
AJG designed the targeting panel, analyzed and interpreted the data, and was the primary writer of the manuscript. ND performed the data filtering and SKAT-O analysis. DVB conducted the primary sequencing bioinformatics processing and interpretation. MAS created the scripts to identify dual-hit genes and create summaries of the data. SHS performed the annotation and data filtering. ERM supervised all bioinformatics and statistical analyses. MLC coordinated the clinical evaluation of all subjects. JMJ retrieved the clinical data and managed the sample intake. WH managed all sequencing experiments. IK managed all sample preparations. PW managed the DNA extraction and allocation. JAR performed all molecular experiments including DNA QC and Sanger validation. JLH provided input into the study design and statistical analysis. JRG provided input into the study design as well as the molecular analysis and interpretation. JPH helped to conceive the project and aid in the analysis and interpretation. MPV contributed to the design of the study, development of the methods, coordination of the statistical and molecular analysis, and interpretation of the data. All authors contributed to the writing and reviewing of the manuscript. All authors read and approved the manuscript.
Additional file 1:
Genes overlapping or nearest GWAS-NR association peaks. File includes gene name and hg19 genomic coordinates of the canonical transcript of each targeted gene.
Additional file 2:
Genes with at least two rare (MAF < 0.01), damaging SNVs. File includes the gene name, total count of rare damaging variants in any individual, and the count of cases and controls with at least two.
Additional file 3:
Rare loss-of-function variants. Table includes position, base change, and frequency of each of the 464 loss-of-function variants identified in our cohort.
Rights and permissions
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 https://creativecommons.org/licenses/by/4.0/.
The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Griswold, A.J., Dueker, N.D., Van Booven, D. et al. Targeted massively parallel sequencing of autism spectrum disorder-associated genes in a case control cohort reveals rare loss-of-function risk variants. Molecular Autism 6, 43 (2015). https://doi.org/10.1186/s13229-015-0034-z
- Autism Spectrum Disorder
- Autism Spectrum Disorder
- Massively Parallel Sequencing
- Autism Spectrum Disorder Patient
- Autism Spectrum Disorder Case