Skip to main content

Evidence of novel fine-scale structural variation at autism spectrum disorder candidate loci

Abstract

Background

Autism spectrum disorders (ASD) represent a group of neurodevelopmental disorders characterized by a core set of social-communicative and behavioral impairments. Gamma-aminobutyric acid (GABA) is the major inhibitory neurotransmitter in the brain, acting primarily via the GABA receptors (GABR). Multiple lines of evidence, including altered GABA and GABA receptor expression in autistic patients, indicate that the GABAergic system may be involved in the etiology of autism.

Methods

As copy number variations (CNVs), particularly rare and de novo CNVs, have now been implicated in ASD risk, we examined the GABA receptors and genes in related pathways for structural variation that may be associated with autism. We further extended our candidate gene set to include 19 genes and regions that had either been directly implicated in the autism literature or were directly related (via function or ancestry) to these primary candidates. For the high resolution CNV screen we employed custom-designed 244 k comparative genomic hybridization (CGH) arrays. Collectively, our probes spanned a total of 11 Mb of GABA-related and additional candidate regions with a density of approximately one probe every 200 nucleotides, allowing a theoretical resolution for detection of CNVs of approximately 1 kb or greater on average. One hundred and sixty-eight autism cases and 149 control individuals were screened for structural variants. Prioritized CNV events were confirmed using quantitative PCR, and confirmed loci were evaluated on an additional set of 170 cases and 170 control individuals that were not included in the original discovery set. Loci that remained interesting were subsequently screened via quantitative PCR on an additional set of 755 cases and 1,809 unaffected family members.

Results

Results include rare deletions in autistic individuals at JAKMIP1, NRXN1, Neuroligin4Y, OXTR, and ABAT. Common insertion/deletion polymorphisms were detected at several loci, including GABBR2 and NRXN3. Overall, statistically significant enrichment in affected vs. unaffected individuals was observed for NRXN1 deletions.

Conclusions

These results provide additional support for the role of rare structural variation in ASD.

Background

Autism spectrum disorders (ASD) represent a group of neurodevelopmental disorders characterized by a core set of social-communicative and behavioral impairments. There is abundant evidence for a strong genetic contribution to ASD, with concordance rates among MZ twins ranging from 60% to 90% [1]. It is now evident that the underlying genetic architecture of autism is highly complex, with numerous genomic loci contributing to overall risk. Efforts to identify specific risk loci have met with some success, and variants in several genes and genomic regions having been implicated to date. These loci include SHANK3, CNTNAP2, NLGN4X, pCDH10, 16p11.2, NRXNa1, SYNGAP, SEMA5A, and AVPR1A, among others [2–6]. Evidence for the role of common variation in ASD risk was reported within the CDH9/CDH10 region of chr5 on risk [7]. Despite these important advances, phenotypic and genetic heterogeneity have made it challenging to unravel the underlying causes of ASD, and a substantial component of this disorder's genetic etiology remains to be discovered.

Copy number variation and ASD

A number of studies have implicated copy number variations (CNVs) in ASD [8–15]. In addition to specific structural mutations, global enrichment in rare and/or de novo CNVs has been reported within affected vs. unaffected individuals (for example [5, 7–15]). To date, studies examining copy number variation in ASD have relied heavily on array technologies that were originally designed for the purpose of conducting SNP-based genome-wide association studies (GWAS). This is primarily due to the fact these data already exist (that is they were obtained for genotyping purposes). These GWAS arrays typically place a limit on the resolution at which CNVs can be reliably detected at 10 to 50 kb or larger (contingent upon the array and analysis parameters), although smaller events can be detected in certain circumstances. The use of custom designed comparative genomic hybridization (CGH) arrays, combined with a targeted, candidate-gene oriented strategy, can extend the ability to reliably detect variation below 10 kb. As it is evident that structural DNA variation with potential relevance to pathological phenotypes can and does exist at sub-10 kb scales (for example [16]), we sought to examine a subset of autism candidate regions at high resolution. The GABA receptor family of genes and genes within biological pathways associated with GABA are attractive candidates for autism and ASD. GABA is the primary inhibitory neurotransmitter in the human brain, and it is known to assume both excitatory and inhibitory roles during early development stages [17, 18]. Previously our group has reported evidence of both linkage and association with ASD at the GABA receptor genes in the 15q11-13 region [19, 20]. In addition, differential expression of GABA receptors has been observed in ASD cases [21–23], and rare coding variants in GABRB3 have been found that segregate with ASD [24]. In light of the above evidence, we chose to examine 24 GABA receptors and 19 additional GABA-related genes for evidence of sub-5 kb structural variation that might contribute to the autism/ASD risk. We further extended our candidate gene set to include an additional 19 genes and regions that had been previously implicated in autism (Table 1). These loci were either recognized ASD candidates or directly related to ASD candidates through biological function or sequence homology. Targeted genes include members of the Neurexin and Neuroligin gene families, the CDH9/CHD10 region on chromosome 5p, the 16p11.2 regions, CNTNAP2, OXT, OXTR, MECP2, and SHANK3. While most of these candidate genes and regions have been previously examined for structural variation in the context of GWAS array analysis, these loci have not been systematically targeted with custom designed high resolution aCGH in a population of ASD and control individuals to survey fine-scale structural variation.

Table 1 Targeted genes and regions (hg18 coordinates)

Methods

Ethical approval

This research was carried out in accordance with the World Medical Association's Declaration of Helsinki. Participants in the study were ascertained under applicable IRB protocols of the John P. Hussman Institute for Human Genetics (HIHG), University of Miami, Miller School of Medicine, Vanderbilt University Medical Center, and University of California, San Francisco (UCSF) Medical School. Following a description of the study, informed consent was obtained for each individual or, where appropriate, their guardian.

Ascertainment and DNA extraction

Individuals were recruited as part of our ongoing family-based recruitment via support groups, advertisements, and clinical and educational settings. Written informed consent was obtained from parents for all minor children and those who were unable to give consent. In addition, we obtained assent from all participants of the appropriate developmental and chronological age. Core inclusion criteria were as follows: (1) between 3 and 21 years of age, (2) a clinical diagnosis of ASD, (3) an expert clinical determination of ASD diagnosis using DSM-IV criteria supported by the Autism Diagnostic Interview (ADI-R) [25, 26], and (4) an IQ equivalent > 35 or developmental level >18 months as determined by the Vineland Adaptive Behavior Scale (VABS) [27]. A best estimate diagnostic determination was based on review by a three-person panel, including experienced clinical psychologists and a pediatric medical geneticist. This determination included ADI-R results. In those instances where an ADI-R was not available (n = 86; six were parents, 67 were lost to follow-up, refused to complete the interview, or withdrew, and 13 are in progress), a best-estimate diagnosis was assigned by this panel using all available clinical information including clinician summaries, caregiver report, and medical records. Following review of case materials and discussion, panel members reached a consensus. We excluded participants with severe sensory problems (for example visual impairment or hearing loss), significant motor impairments (for example failure to sit by 12 months or walk by 24 months), or identified metabolic, genetic, or progressive neurological disorders. Family history and pedigree information was collected in a standard semi-structured interview from a knowledgeable informant, frequently the mother. Additional clinical data were also collected by reviewing available medical and psychiatric records of affected individuals.

The majority of control participants were recruited by the HIHG at the University of Miami and included 565 children between the ages of 3 and 21 years. Participants with developmental, behavioral, or neurological conditions were excluded, as well as those with first-degree relatives with such disorders. We obtained consent from all participants or, in the case of minors, their parents. Participants meeting these criteria provided a saliva sample. A knowledgeable informant, usually the mother, completed the Social Communication Questionnaire to screen for potential ASDs [28]. The above controls were supplemented with additional set of 65 individuals obtained from the UCSF School of Medicine. Blood samples were obtained from self-reported Non-white Hispanics, ranging from 15 to 37 years of age and reporting no history of chronic disease.

For the CNV discovery stage of the project, DNA from whole blood and/or saliva (Oragene kit) was collected from 173 Caucasian non-Hispanic and Caucasian Hispanic autism cases and 185 racially matched control individuals (≤37 years of age) with no known signs of ASD (matched backgrounds were confirmed with Eigenstrat clustering [29]). As described below, data from some samples were excluded from further analysis based on laboratory QC metrics, resulting in a final sample set of 168 cases and 149 controls. CNV loci that were molecularly confirmed and trended towards significance were subsequently tested within a larger sample set of 2,564 individuals, consisting of 755 affected (728 Caucasian/non-Hispanic, 17 Asian, 7 Hispanics, 3 African-Americans) and 1,809 unaffected individuals who were family members of the 755 cases.

DNA from both blood and saliva sources was extracted using the Autopure (Gentra) system according to the manufacturer's standard protocols. For a subset of six control individuals for which DNA extracted from both blood and saliva sources was available, we ran replicate CGH arrays on both blood and saliva DNA in order to look for any systematic differences evident between the two sources. Although batch effects attributable to numerous experimental factors have been previously noted for CNV studies, quality metrics and total CNV events called indicated no significant DNA-source effects within this CGH-based study. Specifically, Nexus 4.1 (Biodiscovery) noise metrics and CNV copy call numbers showed no systematic differences between blood and saliva groups (Additional file 1: Table S1).

Target gene selection

In addition to the set of GABA-related gene targets forming the primary focus of this study, a set of 19 additional ASD candidates were selected based on several criteria, including a recently completed GWAS [7]. The list included both a set of genes previously implicated in the ASD literature, as well as genes known to interact biologically and/or exhibit significant homology to those genes. As not all putative and/or established ASD loci could be accommodated on the array, the final prioritized list was arrived at by consensus among study collaborators.

CGH array design

We designed a custom 244 k Agilent comparative genomic hybridization (CGH) array targeting our regions of interest. Array hybridization probes were chosen using Agilent earray software https://earray.chem.agilent.com/earray/, with a preference for catalogue probes, where available. Probe design was based on coordinates from the hg18 build of the human reference sequence. Collectively, our probes spanned a total of approximately 11 Mb with an average probe density of one unique probe every 200 nucleotides. Given that most standard CNV calling approaches require four to five consecutive significant probes to identify an aberrant region, our coverage theoretically allows for the detection (on average) of CNVs 800 bp or larger. The GABA family of receptors poses particular challenges for hybridization-based copy number assays. Due to the level of sequence homology resulting from the ancestral relationships among GABA genes, increased background probe noise within conserved regions was anticipated. In an effort to mitigate the influence of assay noise on the CNV calling algorithm, we employed triplication of all internal experimental probes as well as running technical array replicates of a subset of 40 individuals to aid in identifying and excluding problematic regions.

Array processing

Genomic DNA was fragmented, enzymatically labeled, hybridized to arrays, and washed according to the manufacturer's (Agilent) standard protocols. Briefly, 0.75 μg of extracted DNA from a single hybridization reference individual (sample NA10851, available from the US National Institute of General Medical Sciences (NIGMS) Human Genetic Cell Repository) was labeled with Cy3 dye (following recommendation of [30]). An equal amount of extracted DNA from experimental samples was enzymatically labeled with Cy5 dye. Labeled samples for each experimental individual was combined with labeled reference DNA, hybridized to the CGH array for 40 h at 65°C, and washed to remove non-hybridized fragments prior to imaging.

Imaging and analysis

All arrays from this study were imaged using a single Molecular Devices GenePix 4000B microarray scanner, processed with Agilent Feature Extraction 9.5 and analyzed with Nexus (Biosdiscovery) software version 4.1. Preliminary analysis was conducted using Rank Segmentation (an extension of the circular binary segmentation algorithm [31] as implemented in Nexus 4.1. Final results presented in this manuscript are based on the Nexus Rank segmentation algorithm with the gain/loss LogR threshold set to 0.25, maximum inter-probe distance at 25 kb, and the (default) minimum of five probes required per called CNV event.

Quality control

As indicated above, each probe was replicated three times (each at a randomly distributed location) on the CGH array and fluorescent signal intensities were averaged prior to subsequent analysis. For detection of common sample handling errors, markers on the X chromosome were examined to ensure observed sample sex was consistent with the original sample manifest. In addition, sex chromosome markers, as a proxy for true-positive controls, a set of five highly polymorphic CNV locations reported in the Database of Genomic Variation http://projects.tcag.ca/variation/ were used to evaluate the ability of our array design (and associated processing technique) to detect structural DNA variation at scales of 5 kb and below (Table 1 bottom). Because the selected loci harbored reported minor allele frequencies of 15% and above, each of the minor structural alleles was expected to be observed, with high probability, one or more times within our sample of 317 individuals. All five of the targeted CNVs were readily detected within our sample population. Although it was observed as variant in our sample set, one of the five targets chosen from the database, hg18 locus chr16:4710525-4715494, was found upon further examination to harbor extreme high repetitive content making the reliability of hybridization results questionable. We therefore do not recommend its use as a positive control locus in future experiments.

To identify potential outlier samples (with respect to sample quality) total CNV calls per individual were rank sorted. A clear demarcation point was observed at approximately 60 total CNV calls per individual, with a relatively small number of samples exhibiting significantly higher numbers of CNV calls. We suspected that the excess number of calls associated with these arrays was attributable to poor sample quality or processing issues; as evidence of this, the Nexus 4.1 quality metric (a measure of probe-to-probe variation across contiguous probes) among these individuals was elevated compared to that of the non-excluded set. Collectively, our quality criteria resulted in a subset of 41 samples being excluded from all subsequent analysis.

Statistical analysis of CNV calls

Initial CNV calls output from Nexus 4.1 rank segmentation algorithm were converted to compatible format for loading into PLINK [32]. For the purpose of analysis, PLINK provides the option of handling heterogeneous CNV calls at the same genomic locus by treating each unique CNV start and stop coordinate as separate markers. Any copy number gain (or loss) that overlaps the position of a marker assumes a value (0, 1, 2, 3, or 4-plus copies) for the marker. After initial processing in PLINK to produce a .ped file with CNV markers, logistical regression analysis was conducted via R scripting. Results of logistic regression analysis are provided in Additional file 2: Table S2. Regression results corresponding to -LOG10 (P value) ≥ 3 were evaluated as candidates for additional analysis. For analysis of the larger family-based dataset, the GENMOD procedure (SAS vers8.1) was run on the dataset and a generalized estimating equation (GEE) was employed to account for the relationships between the parents and siblings (unaffected) and cases. GEE allows for specification of a correlation matrix of the between-subject relatedness of participants belonging to the same pedigree [33] and is commonly used in the context of genetic epidemiology studies to produce robust variance estimates that otherwise would likely be too small/conservative in the presence of relatedness between observations [34]. For each of the candidate genes evaluated by GEE at the validation stage of the study, Bonferroni correction was used to correct for the number CNV positions examined within the candidate. Global CNV burden was assessed using the mperm option in PLINK to perform a global permutation test (10,000 permutations) on segment number, total kilobase difference, and average kilobase difference.

Molecular confirmation of CNV calls

Prioritized loci were validated by real-time quantitative PCR using TaqMan chemistry (a full list of target locations are provided in Additional file 3: Table S3). Assays were performed in quadruplicate, as recommended by the manufacturer (Applied Biosystems {Life Technologies, Carlsbad, CA}), reactions were run on the ABI-7900HT real-time PCR system and analyzed using CopyCaller software v.1.0 (Life Technologies, Carlsbad, CA). Regions were chosen for validation based on several criteria, including their potential for biological relevance and logistic regression scores (Additional file 2: Table S2). CNVs spanning known exons were given the highest priority, followed by those including regions of evolutionary conservation, as evidenced by interspecies sequence homology (per UCSC Genome Browser, build 36.1, March 2006, conservation tracks). Some predicted copy number variable regions that exhibited statistically significant scores were nevertheless excluded from further analysis because data from technical replicates indicated poor replicability for calls at these regions. Examination of these loci indicated that they tended to exhibit extreme GC content, repetitive element content, or extensive homology to other areas of the genome likely due to common ancestry and subsequent evolutionary conservation (data not shown). In addition to the above criteria, CNV regions were excluded from follow-up if they exibited 75% or greater overlap with known copy number variable regions reported in DGV http://projects.tcag.ca/variation/; however, exceptions were made in cases where the variants encompassed exons within candidate genes and instances where gained or loss alleles were discordant in case/controls. All loci meeting these criteria subjected to TaqMan CNV assays (Applied Biosystems) for further interrogation in additional cases and controls.

Results and discussion

A total of 24 GABA receptors, 19 GABA-associated coding genes, and 19 additional autism candidates were targeted for aCGH assessment (Table 1). After exclusion of outlier samples, a total of 168 case and 149 control individuals were examined by CGH arrays in the initial discovery set. In all, 6,143 CNV events calls were identified across our tested population, representing 1,223 potentially unique variations. The overall number of CNV events did not exhibit significant differences between ASD affected and unaffected individuals, although both the average and total kilobases impacted did show statistical significance (perm P < 0.001; 10,000 permutations; one-sided), consistent with previous reports indicating pathogenic CNVs tend to be larger in size and over-represented in affected individuals (for example [35]). The full list of CNV calls from the Nexus Rank Segmentation algorithm is provided in Additional file 4: Table S4. The logistic regression analysis of case and control allele frequencies were used, in conjunction with additional biological annotation criteria, to prioritize a subset of regions for confirmation via a second molecular platform, TaqMan copy number quantitative PCR assays (Applied Biosystems). In addition to the ASD candidate loci prioritized for follow-up, five control genomic regions, representing known polymorphic locations with minor allele frequencies >10%, were also targeted with TaqMan assays to test the sensitivity of our laboratory technique. The 29 targeted regions and five control loci targeted for follow-up are detailed in Additional file 3: Table S3, along with confirmation results. A diagram of the work flow for discovery and validation is provided in Additional file 5: Figure S1. The sample set used for molecular confirmation and additional population screening was comprised of the individuals in which the CNVs were initially observed (for confirmation), along with an additional 170 Caucasian cases and 170 unaffected ethnically-matched controls that were not present in the original discovery set (Table 2). A subset of independently confirmed loci that remained interesting was run on a larger set of 755 case and 1,809 unaffected family members. GEE analysis results for these data, accounting for relationships between individuals, are provided in Table 3.

Table 2 Taqman case/control confirmation set data
Table 3 GEE results from family-based validation dataset (n = 2,564, 755 affected and 1,809 unaffected individuals)

Targeted CNV call events were confirmed with quantitative PCR at 13 of 29 loci, 45%. This suggests a minimum false-positive rate (at the locus level) of 55% (addressed below).

Phenotypic evaluation

Medical records were examined for a subset of individuals harboring rare CNV events at loci of interest (ABAT, NRXN1, JAKMIP1, GABBR2, NLGN2, OXTR). We examined salient phenotypic and medical history characteristics for these individuals to look for commonalities among carriers of CNVs at particular loci. Characteristics examined included developmental abnormalities in speech and motor development, learning disabilities, and evidence of dysmorphologies. Additional co-morbidities (for example irritable bowel syndrome, asthma) and neuropsychiatric features were also noted. No obvious correspondence between aberrant CNV state and specific phenotypic traits was observed. A summary of available phenotypic and medical characteristics for the individuals examined is provided in Additional file 6: Table S5.

In addition to these rare variants, common structural polymorphisms were observed and validated at multiple loci, including NRXN3, CNTNAP2, GABBR2, SLC38A3, and NLGN2. Of these more frequent events, only the GABBR2 locus approached nominal significance (GEE; P = 0.029; OR = 1.16; 95% CI: 0.7-1.93) in our larger validation set, but the result was not significant at α = 0.025 after correction for multiple testing (Table 3).

The role of structural variation, particularly rare and de novo CNVs, in autism has been established in a number of previous studies. Here we examined a set of autism candidate loci using high-resolution aCGH analysis. Our principal objective was to detect structural variants at these loci that may have been potentially missed by other screening platforms, such as SNP-based arrays. Among the limitations impacting our study was the presence of a high false-positive rate (in excess of 50%) in our initial aCGH discovery set. The most likely explanation for this elevated false-positive rate is that the majority of genes in our study, including both GABA and non-GABA-related targets (for example MBD1-MBD4), belong to families of homologous genes which can lead to aberrant probe hybridization. As the majority of loci in our study belonged to homology groups, we could not formally test for this effect. In addition, the modest overall sample numbers used in this study limited our ability to accurately estimate population allele frequencies for rarer events. Even with these complicating factors, several promising events were detected and independently confirmed by quantitative PCR. Rare structural variations were detected at ABAT, NLGN4-Y, NRXN1, and JAKMIP1 (Table 2). Most notable among these results, statistically significant enrichment of deletions (GEE; P = 0.0215; OR = 2.75; 95% CI: 1.161-6.5; α = 0.025) in affected vs. unaffected individuals was observed at NRXN1 within our larger family-based validation set with after correction for multiple testing (Table 3).

We observed and confirmed an approximate 4 kb deletion in a female autistic indvidual within the JAKMIP1 gene that removes two exons. JAKMIP1, a member of a family of microtubule interacting proteins, is known to interact with GABBR1 protein [36] and knockdown of its expression via siRNA results in increased GABBR2 expression, suggesting a possible role for JAKMIP1 in GABBR2 regulation. Misregulation of JAKMIP1 has also been observed in lymphoblastoid cell lines from autistic individuals [37]. We note, however, that no significant enrichment for the deletion was observed in our larger validation set. More detailed characterization of these deletions in the affected and unaffected individual is warranted before further inference can be made, as the precise boundaries of the deletion events may be different in case and control individuals.

A large (approximately 70 kb) deletion encompassing two exons of the NRXN1 genes was detected in an autistic male (Figure 1). NRXN1 is a neuronal adhesion molecule that has been implicated in a number of neurological disorders, including autism (for example [38–41]. Suggestive SNP association results at the NRXN1 locus were reported in a large collaborative study, but genome-wide significance was not reached [7]. Significant enrichment was previously reported for rare structural events [42]. The removal of this exonic region within an autistic male in our dataset, along with the detection of significant enrichment of NRXN1 deletions in cases vs. unaffected individuals in this study, provides further evidence for NRXN1 as an ASD risk locus. Furthermore, the location of our CNV, along with the placement of previously reported NRXN1 events (Figure 1), suggest that disruptions of the alpha form of NRXN1. The beta form of NRXN1, which is transcribed from a secondary downstream promoter, is less frequently observed to be disrupted by CNV in ASD cases. Combined with the aforementioned genetic association results, our data suggests that both rare and common variants may be contributing to ASD risk at this genomic location.

Figure 1
figure 1

A deletion of approximately 70 kb encompassing two exons (3,4) from the NRXN1 gene (alpha-1-isoform). Previously reported NRXN1 CNV events in ASD cases are shown mapped below [38–40, 42–49]. A subset of NRXN1 structural variants reported in individuals with mild to severe intellectual disability is also depicted [49]. Solid lines indicate deletions (majority) and dashed lines indicate duplication events.

A deletion of the entire Neuroligin 4Y gene within an autistic male was also observed and confirmed. In the larger dataset, the deletion was observed in the original individual's unaffected father as well as within an unaffected male sibling. Two additional unaffected carriers (unrelated) were also observed in the larger dataset. Interestingly, however, we observed two duplication events within the larger validation set, both of which were present in unrelated ASD individuals. Neuroligin 4Y, which resides on the non-recombining portion of the Y chromosome in modern humans was, in the evolutionary past, allelic with Neuroligin 4X (NLGN4X) gene. Several studies have previously implicated NLGN4X in autism (for example [50–52]). Nucleotide mutations in NLGN4Y have been reported in an autistic male and his father, who possessed learning disabilities [53]. While there is currently no data regarding the functional status and possible role of NLGN4Y, EST evidence from the NCBI database suggests that NLGN4Y retains a molecularly active promoter and is expressed in the brain. A recent report demonstrated a potential role for some pseudogenes in regulating the expression of their parent or source genes by providing decoy targets for miRNAs [54]. It is possible that, through this or some other mechanism, NLGN4Y is involved in regulating the expression of the NLGN4X locus in males. Our finding of copy number aberrations of this gene in autistic individuals suggests that the locus warrants further investigation.

We observed a deletion at the ABAT gene locus within multiple ASD individuals, one in our discovery dataset and a second case ASD case in the validation set. Further examination of the locus within the larger, family based set revealed 7 of 755 affected individuals carrying the deletion and five (of 1,809) unaffected individuals (unrelated to the seven affected individuals). Although we observed a trend of enrichment in affected vs. unaffected individuals in our dataset, statistical significance was not met at α = 0.05. The observed trend, however, suggests that ABAT locus, and the GABA system as a whole, warrant further scrutiny for their potential role in ASD risk. ABAT is involved in the catabolism of GABA, has previously been implicated in autism by genetic association, although replication of the association was not observed in the small validation dataset (n = 91 autism trios) [55]. In the one deleted individual from our dataset for which CGH array data was available at the locus, the ABAT deletion removes exon 1 of the gene (Additional file 7: Figure S2). We note that the nearby A2BP1 and GRIN2A loci, which have also been implicated in autism [55], were not included in our deleted region (Additional file 8: Figure S3).

We also observed a previously reported deletion in the OXTR receptor in a single male ASD individual [56]. A smaller duplication region in OXTR was also detected in a separate ASD case.

In addition to the rare copy number events described above, we also observed and confirmed common copy number polymorphisms at SLC38A3, GABBR2, NRXN3, SHANK3, and NLGN2 (Table 2). Among these, only duplication events in GABBR2 achieved nominal statistical significance at α = 0.05, but this result failed to hold up after multiple correction for multiple testing. Furthermore, the relatively high level of background noise of assays targeting this locus makes us approach the GABBR2 finding with caution. Additional examination of the locus through alternative platforms (for example long-read high throughput sequence technologies) will be necessary before further conclusions can be drawn.

Conclusions

These data further indicate the importance of structural variation in ASD risk and provide additional evidence that rare variants at multiple genomic loci are potentially contributing to this common neurodevelopmental disorder. In particular, we report statistically significant enrichment of rare exonic deletions in NRXN1 in autistic vs. non-autistic individuals. The detection and confirmation of structural variants below the 10-50 kb typically afforded by SNP-based GWAS arrays indicate that aCGH remains an important complementary method for CNV detection.

Availability of supporting data

The dataset(s) supporting the results of this article are included within additional supplemental files.

Abbreviations

aCGH:

Array-based Comparative Genomic Hybridization

ADI-R:

Autism Diagnostic Interview

ASD:

Autism Spectrum Disorders

BP:

Base Pair(s)

CGH:

Comparative Genomic Hybridization

CNV:

Copy Number Variation

DNA:

Deoxyribonucleic Acid

DSM-IV:

Diagnostic and Statistical Manual of Mental Disorders IV

GABA:

Gamma-aminobutyric acid

GABR:

Gamma-aminobutyric acid receptors

GEE:

Generalized Estimating Equation

GWAS:

Genome Wide Association Study

HIHG:

Hussman Institute of Human Genomics

IQ:

Intelligence Quotient

MZ:

Monozygotic

NIGMS:

National Institute of General Medical Sciences

PCR:

Polymerase Chain Reaction

UCSF:

University of California, San Francisco

VABS:

Vineland Adaptive Behavior Scale.

References

  1. Bailey A, Le Couteur A, Gottesman I, Bolton P, Simonoff E, Yuzda E, Rutter M: Autism as a strongly genetic disorder: evidence from a British twin study. Psychol Med. 1995, 25: 63-77. 10.1017/S0033291700028099.

    Article  CAS  PubMed  Google Scholar 

  2. Meyer-Lindenberg A, Kolachana B, Gold B, Olsh A, Nicodemus KK, Mattay V, Dean M, Weinberger DR: Genetic variants in AVPR1A linked to autism predict amygdala activation and personality traits in healthy humans. Mol Psychiatry. 2008, 14: 968-975.

    Article  PubMed Central  PubMed  Google Scholar 

  3. Morrow EM, Yoo SY, Flavell SW, Kim TK, Lin Y, Hill RS, Mukaddes NM, Balkhy S, Gascon G, Hashmi A, Al-Saad S, Ware J, Joseph RM, Greenblatt R, Gleason D, Ertelt JA, Apse KA, Bodell A, Partlow JN, Barry B, Yao H, Markianos K, Ferland RJ, Greenberg ME, Walsh CA: Identifying autism loci and genes by tracing recent shared ancestry. Science. 2008, 321: 218-223. 10.1126/science.1157657.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Yan J, Noltner K, Feng J, Li W, Schroer R, Skinner C, Zeng W, Schwartz CE, Sommer SS: Neurexin 1alpha structural variants associated with autism. Neurosci Lett. 2008, 438: 368-370. 10.1016/j.neulet.2008.04.074.

    Article  CAS  PubMed  Google Scholar 

  5. Pinto D, Pagnamenta AT, Klei L, Anney R, Merico D, Regan R, Conroy J, Magalhaes TR, Correia C, Abrahams BS, Almeida J, Bacchelli E, Bader GD, Bailey AJ, Baird G, Battaglia A, Berney T, Bolshakova N, Bolte S, Bolton PF, Bourgeron T, Brennan S, Brian J, Bryson SE, Carson AR, Casallo G, Casey J, Chung BH, Cochrane L, Corsello C: Functional impact of global rare copy number variation in autism spectrum disorders. Nature. 2010, 466: 368-372. 10.1038/nature09146.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. 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-808. 10.1038/nature08490.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Wang K, Zhang H, Ma D, Bucan M, Glessner JT, Abrahams BS, Salyakina D, Imielinski M, Bradfield JP, Sleiman PM, Kim CE, Hou C, Frackelton E, Chiavacci R, Takahashi N, Sakurai T, Rappaport E, Lajonchere CM, Munson J, Estes A, Korvatska O, Piven J, Sonnenblick LI, Alvarez Retuerto AI, Herman EI, Dong H, Hutman T, Sigman M, Ozonoff S, Klin A: Common genetic variants on 5p14.1 associate with autism spectrum disorders. Nature. 2009, 459: 528-533. 10.1038/nature07999.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Rudan I: New technologies provide insights into genetic basis of psychiatric disorders and explain their co-morbidity. Psychiatr Danub. 2010, 22: 190-192.

    PubMed  Google Scholar 

  9. Fernandez BA, Roberts W, Chung B, Weksberg R, Meyn S, Szatmari P, Joseph-George AM, Mackay S, Whitten K, Noble B, Vardy C, Crosbie V, Luscombe S, Tucker E, Turner L, Marshall CR, Scherer SW: Phenotypic spectrum associated with de novo and inherited deletions and duplications at 16p11.2 in individuals ascertained for diagnosis of autism spectrum disorder. J Med Genet. 2010, 47: 195-203. 10.1136/jmg.2009.069369.

    Article  PubMed  Google Scholar 

  10. Bucan M, Abrahams BS, Wang K, Glessner JT, Herman EI, Sonnenblick LI, Alvarez Retuerto AI, Imielinski M, Hadley D, Bradfield JP, Kim C, Gidaya NB, Lindquist I, Hutman T, Sigman M, Kustanovich V, Lajonchere CM, Singleton A, Kim J, Wassink TH, McMahon WM, Owley T, Sweeney JA, Coon H, Nurnberger JI, Li M, Cantor RM, Minshew NJ, Sutcliffe JS, Cook EH: Genome-wide analyses of exonic copy number variants in a family-based study point to novel autism susceptibility genes. PLoS Genet. 2009, 5: e1000536-10.1371/journal.pgen.1000536.

    Article  PubMed Central  PubMed  Google Scholar 

  11. Cho SC, Yim SH, Yoo HK, Kim MY, Jung GY, Shin GW, Kim BN, Hwang JW, Kang JJ, Kim TM, Chung YJ: Copy number variations associated with idiopathic autism identified by whole-genome microarray-based comparative genomic hybridization. Psychiatr Genet. 2009, 19: 177-185. 10.1097/YPG.0b013e32832bdafa.

    Article  PubMed  Google Scholar 

  12. Sykes NH, Toma C, Wilson N, Volpi EV, Sousa I, Pagnamenta AT, Tancredi R, Battaglia A, Maestrini E, Bailey AJ, Monaco AP, International Molecular Genetic Study of Autism Consortium (IMGSAC): Copy number variation and association analysis of SHANK3 as a candidate gene for autism in the IMGSAC collection. Eur J Hum Genet. 2009, 17: 1347-1353. 10.1038/ejhg.2009.47.

    Article  PubMed Central  PubMed  Google Scholar 

  13. Christian SL, Brune CW, Sudi J, Kumar RA, Liu S, Karamohamed S, Badner JA, Matsui S, Conroy J, McQuaid D, Gergel J, Hatchwell E, Gilliam TC, Gershon ES, Nowak NJ, Dobyns WB, Cook EH: Novel submicroscopic chromosomal abnormalities detected in autism spectrum disorder. Biol Psychiatry. 2008, 63: 1111-1117. 10.1016/j.biopsych.2008.01.009.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Marshall CR, Noor A, Vincent JB, Lionel AC, Feuk L, Skaug J, Shago M, Moessner R, Pinto D, Ren Y, Thiruvahindrapduram B, Fiebig A, Schreiber S, Friedman J, Ketelaars CE, Vos YJ, Ficicioglu C, Kirkpatrick S, Nicolson R, Sloman L, Summers A, Gibbons CA, Teebi A, Chitayat D, Weksberg R, Thompson A, Vardy C, Crosbie V, Luscombe S, Baatjes R: Structural variation of chromosomes in autism spectrum disorder. Am J Hum Genet. 2008, 82: 477-488. 10.1016/j.ajhg.2007.12.009.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Sebat J, Lakshmi B, Malhotra D, Troge J, Lese-Martin C, Walsh T, Yamrom B, Yoon S, Krasnitz A, Kendall J, Leotta A, Pai D, Zhang R, Lee YH, Hicks J, Spence SJ, Lee AT, Puura K, Lehtimaki T, Ledbetter D, Gregersen PK, Bregman J, Sutcliffe JS, Jobanputra V, Chung W, Warburton D, King MC, Skuse D, Geschwind DH, Gilliam TC: Strong association of de novo copy number mutations with autism. Science. 2007, 316: 445-449. 10.1126/science.1138659.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Perry GH, Ben-Dor A, Tsalenko A, Sampas N, Rodriguez-Revenga L, Tran CW, Scheffer A, Steinfeld I, Tsang P, Yamada NA, Park HS, Kim JI, Seo JS, Yakhini Z, Laderman S, Bruhn L, Lee C: The fine-scale and complex architecture of human copy-number variation. Am J Hum Genet. 2008, 82: 685-695. 10.1016/j.ajhg.2007.12.010.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Ben-Ari Y, Gaiarsa JL, Tyzio R, Khazipov R: GABA: a pioneer transmitter that excites immature neurons and generates primitive oscillations. Physiol Rev. 2007, 87: 1215-1284. 10.1152/physrev.00017.2006.

    Article  CAS  PubMed  Google Scholar 

  18. Takayama C, Inoue Y: GABAergic signaling in the developing cerebellum. Anat Sci Int. 2004, 79: 124-136. 10.1111/j.1447-073x.2004.00081.x.

    Article  CAS  PubMed  Google Scholar 

  19. Menold MM, Shao Y, Wolpert CM, Donnelly SL, Raiford KL, Martin ER, Ravan SA, Abramson RK, Wright HH, DeLong GR, Cuccaro ML, Pericak-Vance MA, Gilbert JR: Association analysis of chromosome 15 gabaa receptor subunit genes in autistic disorder. J Neurogenet. 2001, 15: 245-259. 10.3109/01677060109167380.

    Article  CAS  PubMed  Google Scholar 

  20. Bass MP, Menold MM, Wolpert CM, Donnelly SL, Ravan SA, Hauser ER, Maddox LO, Vance JM, Abramson RK, Wright HH, Gilbert JR, Cuccaro ML, DeLong GR, Pericak-Vance MA: Genetic studies in autistic disorder and chromosome 15. Neurogenetics. 2000, 2: 219-226. 10.1007/s100489900081.

    Article  CAS  PubMed  Google Scholar 

  21. Fatemi SH, Reutiman TJ, Folsom TD, Thuras PD: GABA(A) receptor downregulation in brains of subjects with autism. J Autism Dev Disord. 2009, 39: 223-230. 10.1007/s10803-008-0646-7.

    Article  PubMed Central  PubMed  Google Scholar 

  22. Fatemi SH, Folsom TD, Reutiman TJ, Thuras PD: Expression of GABA(B) receptors is altered in brains of subjects with autism. Cerebellum. 2009, 8: 64-69. 10.1007/s12311-008-0075-3.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Oblak AL, Gibbs TT, Blatt GJ: Decreased GABA(B) Receptors in the Cingulate Cortex and Fusiform Gyrus in Autism. J Neurochem. 2010, 114: 1414-1423.

    PubMed Central  CAS  PubMed  Google Scholar 

  24. Delahanty RJ, Kang JQ, Brune CW, Kistner EO, Courchesne E, Cox NJ, Cook EH, Macdonald RL, Sutcliffe JS: Maternal transmission of a rare GABRB3 signal peptide variant is associated with autism. Mol Psychiatry. 2011, 16: 86-96. 10.1038/mp.2009.118.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Rutter M, LeCouteur A, Lord C: Autism Diagnostic Interview, Revised (ADI-R). 2003, Los Angeles, CA: Western Psychological Services

    Google Scholar 

  26. Lord C, Rutter M, DiLavore P, Risi S: Autism Diagnostic Observation Schedule-WPS. 1999, Los Angeles, CA: Western Psychological Services, WPSth {Version 1AU Query: please confirm the edition number} edition.

    Google Scholar 

  27. Sparrow SS, Balla D, Cicchetti D: Vineland Adaptive Behavior Scales, Interview Edition. 1984, Circle Pines, MN: AGS Publishing

    Google Scholar 

  28. Rutter M, Bailey A, Berument SK, LeCouteur A, Lord C, Pickles A: Social Communication Questionnaire (SCQ). 2003, Los Angeles, CA: Western Psychological Services

    Google Scholar 

  29. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D: Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006, 38: 904-909. 10.1038/ng1847.

    Article  CAS  PubMed  Google Scholar 

  30. Scherer SW, Lee C, Birney E, Altshuler DM, Eichler EE, Carter NP, Hurles ME, Feuk L: Challenges and standards in integrating surveys of structural variation. Nat Genet. 2007, 39: S7-S15. 10.1038/ng2093.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Olshen AB, Venkatraman ES, Lucito R, Wigler M: Circular binary segmentation for the analysis of array-based DNA copy number data. Biostatistics. 2004, 5: 557-572. 10.1093/biostatistics/kxh008.

    Article  PubMed  Google Scholar 

  32. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81: 559-575. 10.1086/519795.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Liang K, Zeger SL: Longitudinal data analysis using generalized linear models. Biometrika. 1986, 73: 13-22. 10.1093/biomet/73.1.13.

    Article  Google Scholar 

  34. Oxford Statistical Science Series, Diggle PJ, Heagerty P, Liang K, Zeger SL: Analysis of Longitudinal Data. 2002, Oxford Statistical Science Series

    Google Scholar 

  35. Cooper GM, Coe BP, Girirajan S, Rosenfeld JA, Vu TH, Baker C, Williams C, Stalker H, Hamid R, Hannig V, Abdel-Hamid H, Bader P, McCracken E, Niyazov D, Leppig K, Thiese H, Hummel M, Alexander N, Gorski J, Kussmann J, Shashi V, Johnson K, Rehder C, Ballif BC, Shaffer LG, Eichler EE: A copy number variation morbidity map of developmental delay. Nat Genet. 2011, 43: 838-846. 10.1038/ng.909.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Couve A, Restituito S, Brandon JM, Charles KJ, Bawagan H, Freeman KB, Pangalos MN, Calver AR, Moss SJ: Marlin-1, a novel RNA-binding protein associates with GABA receptors. J Biol Chem. 2004, 279: 13934-13943. 10.1074/jbc.M311737200.

    Article  CAS  PubMed  Google Scholar 

  37. Nishimura Y, Martin CL, Vazquez-Lopez A, Spence SJ, Alvarez-Retuerto AI, Sigman M, Steindler C, Pellegrini S, Schanen NC, Warren ST, Geschwind DH: Genome-wide expression profiling of lymphoblastoid cell lines distinguishes different forms of autism and reveals shared pathways. Hum Mol Genet. 2007, 16: 1682-1698. 10.1093/hmg/ddm116.

    Article  CAS  PubMed  Google Scholar 

  38. Kim HG, Kishikawa S, Higgins AW, Seong IS, Donovan DJ, Shen Y, Lally E, Weiss LA, Najm J, Kutsche K, Descartes M, Holt L, Braddock S, Troxell R, Kaplan L, Volkmar F, Klin A, Tsatsanis K, Harris DJ, Noens I, Pauls DL, Daly MJ, MacDonald ME, Morton CC, Quade BJ, Gusella JF: Disruption of neurexin 1 associated with autism spectrum disorder. Am J Hum Genet. 2008, 82: 199-207. 10.1016/j.ajhg.2007.09.011.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Szatmari P, Paterson AD, Zwaigenbaum L, Roberts W, Brian J, Liu XQ, Vincent JB, Skaug JL, Thompson AP, Senman L, Feuk L, Qian C, Bryson SE, Jones MB, Marshall CR, Scherer SW, Vieland VJ, Bartlett C, Mangin LV, Goedken R, Segre A, Pericak-Vance MA, Cuccaro ML, Gilbert JR, Wright HH, Abramson RK, Betancur C, Bourgeron T, Gillberg C, Leboyer M: Mapping autism risk loci using genetic linkage and chromosomal rearrangements. Nat Genet. 2007, 39: 319-328. 10.1038/ng1985.

    Article  CAS  PubMed  Google Scholar 

  40. Ching MS, Shen Y, Tan WH, Jeste SS, Morrow EM, Chen X, Mukaddes NM, Yoo SY, Hanson E, Hundley R, Austin C, Becker RE, Berry GT, Driscoll K, Engle EC, Friedman S, Gusella JF, Hisama FM, Irons MB, Lafiosca T, LeClair E, Miller DT, Neessen M, Picker JD, Rappaport L, Rooney CM, Sarco DP, Stoler JM, Walsh CA, Wolff RR: Children's Hospital Boston Genotype Phenotype Study Group: Deletions of NRXN1 (neurexin-1) predispose to a wide spectrum of developmental disorders. Am J Med Genet B Neuropsychiatr Genet. 2010, 153B: 937-947.

    PubMed Central  CAS  PubMed  Google Scholar 

  41. Gai X, Xie HM, Perin JC, Takahashi N, Murphy K, Wenocur AS, D'arcy M, O'Hara RJ, Goldmuntz E, Grice DE, Shaikh TH, Hakonarson H, Buxbaum JD, Elia J, White PS: Rare structural variation of synapse and neurotransmission genes in autism. Mol Psychiatry. 2011, doi: 10.1038/mp.2011.10

    Google Scholar 

  42. Glessner JT, Wang K, Cai G, Korvatska O, Kim CE, Wood S, Zhang H, Estes A, Brune CW, Bradfield JP, Imielinski M, Frackelton EC, Reichert J, Crawford EL, Munson J, Sleiman PM, Chiavacci R, Annaiah K, Thomas K, Hou C, Glaberson W, Flory J, Otieno F, Garris M, Soorya L, Klei L, Piven J, Meyer KJ, Anagnostou E, Sakurai T: Autism genome-wide copy number variation reveals ubiquitin and neuronal genes. Nature. 2009, 459: 569-573. 10.1038/nature07953.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. Zweier C, de Jong EK, Zweier M, Orrico A, Ousager LB, Collins AL, Bijlsma EK, Oortveld MA, Ekici AB, Reis A, Schenck A, Rauch A: CNTNAP2 and NRXN1 are mutated in autosomal-recessive Pitt-Hopkins-like mental retardation and determine the level of a common synaptic protein in Drosophila. Am J Hum Genet. 2009, 85: 655-666. 10.1016/j.ajhg.2009.10.004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Zahir FR, Baross A, Delaney AD, Eydoux P, Fernandes ND, Pugh T, Marra MA, Friedman JM: A patient with vertebral, cognitive and behavioural abnormalities and a de novo deletion of NRXN1alpha. J Med Genet. 2008, 45: 239-243.

    Article  CAS  PubMed  Google Scholar 

  45. Rujescu D, Ingason A, Cichon S, Pietilainen OP, Barnes MR, Toulopoulou T, Picchioni M, Vassos E, Ettinger U, Bramon E, Murray R, Ruggeri M, Tosato S, Bonetto C, Steinberg S, Sigurdsson E, Sigmundsson T, Petursson H, Gylfason A, Olason PI, Hardarsson G, Jonsdottir GA, Gustafsson O, Fossdal R, Giegling I, Moller HJ, Hartmann AM, Hoffmann P, Crombie C, Fraser G: Disruption of the neurexin 1 gene is associated with schizophrenia. Hum Mol Genet. 2009, 18: 988-996.

    PubMed Central  CAS  PubMed  Google Scholar 

  46. Wisniowiecka-Kowalnik B, Nesteruk M, Peters SU, Xia Z, Cooper ML, Savage S, Amato RS, Bader P, Browning MF, Haun CL, Duda AW 3rd, Cheung SW, Stankiewicz P: Intragenic rearrangements in NRXN1 in three families with autism spectrum disorder, developmental delay, and speech delay. Am J Med Genet B Neuropsychiatr Genet. 2010, 153B: 983-993.

    CAS  PubMed  Google Scholar 

  47. Sanders SJ, Ercan-Sencicek AG, Hus V, Luo R, Murtha MT, Moreno-De-Luca D, Chu SH, Moreau MP, Gupta AR, Thomson SA, Mason CE, Bilguvar K, Celestino-Soper PB, Choi M, Crawford EL, Davis L, Davis Wright NR, Dhodapkar RM, Dicola M, Dilullo NM, Fernandez TV, Fielding-Singh V, Fishman DO, Frahm S, Garagaloyan R, Goh GS, Kammela S, Klei L, Lowe JK, Lund SC: Multiple recurrent de novo CNVs, including duplications of the 7q11.23 Williams Syndrome region, are strongly associated with autism. Neuron. 2011, 70: 863-885. 10.1016/j.neuron.2011.05.002.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  48. Duong L, Klitten LL, Moller RS, Ingason A, Jakobsen KD, Skjodt C, Didriksen M, Hjalgrim H, Werge T, Tommerup N: Mutations in NRXN1 in a family multiply affected with brain disorders: NRXN1 mutations and brain disorders. Am J Med Genet B Neuropsychiatr Genet. 2012, 159B: 354-358. 10.1002/ajmg.b.32036.

    Article  PubMed  Google Scholar 

  49. Gregor A, Albrecht B, Bader I, Bijlsma EK, Ekici AB, Engels H, Hackmann K, Horn D, Hoyer J, Klapecki J, Kohlhase J, Maystadt I, Nagl S, Prott E, Tinschert S, Ullmann R, Wohlleber E, Woods G, Reis A, Rauch A, Zweier C: Expanding the clinical spectrum associated with defects in CNTNAP2 and NRXN1. BMC Med Genet. 2011, 12: 106-10.1186/1471-2350-12-106.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Jamain S, Quach H, Betancur C, Rastam M, Colineaux C, Gillberg IC, Soderstrom H, Giros B, Leboyer M, Gillberg C, Bourgeron T: Mutations of the X-linked genes encoding neuroligins NLGN3 and NLGN4 are associated with autism. Nat Genet. 2003, 34: 27-29. 10.1038/ng1136.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Laumonnier F, Bonnet-Brilhault F, Gomot M, Blanc R, David A, Moizard MP, Raynaud M, Ronce N, Lemonnier E, Calvas P, Laudier B, Chelly J, Fryns JP, Ropers HH, Hamel BCJ, Andres C, Barthelemy C, Moraine C, Briault S: X-linked mental retardation and autism are associated with a mutation in the NLGN4 gene, a member of the neuroligin family. Am J Hum Genet. 2004, 74: 552-557. 10.1086/382137.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  52. Ylisaukko-oja T, Rehnstrom K, Auranen M, Vanhala R, Alen R, Kempas E, Ellonen P, Turunen JA, Makkonen I, Riikonen R, Nieminen-von Wendt T, von Wendt L, Peltonen L, Jarvela I: Analysis of four neuroligin genes as candidates for autism. Eur J Hum Genet. 2005, 13: 1285-1292. 10.1038/sj.ejhg.5201474.

    Article  CAS  PubMed  Google Scholar 

  53. Yan J, Feng J, Schroer R, Li W, Skinner C, Schwartz CE, Cook EH, Sommer SS: Analysis of the neuroligin 4Y gene in patients with autism. Psychiatr Genet. 2008, 18: 204-207. 10.1097/YPG.0b013e3282fb7fe6.

    Article  PubMed  Google Scholar 

  54. Poliseno L, Salmena L, Zhang J, Carver B, Haveman WJ, Pandolfi PP: A coding-independent function of gene and pseudogene mRNAs regulates tumour biology. Nature. 2010, 465: 1033-1038. 10.1038/nature09144.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  55. Barnby G, Abbott A, Sykes N, Morris A, Weeks DE, Mott R, Lamb J, Bailey AJ, Monaco AP: Candidate-gene screening and association analysis at the autism-susceptibility locus on chromosome 16p: evidence of association at GRIN2A and ABAT. Am J Hum Genet. 2005, 76: 950-966. 10.1086/430454.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  56. Gregory SG, Connelly JJ, Towers AJ, Johnson J, Biscocho D, Markunas CA, Lintas C, Abramson RK, Wright HH, Ellis P, Langford CF, Worley G, Delong GR, Murphy SK, Cuccaro ML, Persico A, Pericak-Vance MA: Genomic and epigenetic evidence for oxytocin receptor deficiency in autism. BMC Med. 2009, 7: 62-10.1186/1741-7015-7-62.

    Article  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements

We appreciate the contributions made to this work by the families of individuals with autism and the research staff at collaborating universities. This study was supported by funding from the National Institute of Neurological Disorders and Stroke (P01NS026630) and the National Institute for Mental Health (5R01MH080647-14), the National Alliance of Autism Research, and the Hussman Foundation. A subset of the participants was ascertained while Margaret A. Pericak-Vance, Michael A. Cuccaro, and John R. Gilbert were faculty members at Duke University.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Dale J Hedges.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

DJH, MPV, JLH, and JRG conceived study. ZMK, BLB, HNC, PLW, and LN performed experiments. DJH, KLH, GWB, DM, and JMJ analyzed experimental data. SJS, LN, JML, and MLC evaluated clinical data. SLH, JRO, and MPV contributed samples and reviewed manuscript. DJH wrote the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

13229_2011_36_MOESM1_ESM.XLSX

Additional file 1: Table S1. Examination of Noise metric and call rates of blood and spit derived DNA extractions from same individuals. (XLSX 12 KB)

Additional file 2: Table S2. (XLSX 2 MB)

Additional file 3: Table S3. Validation Targets. (XLSX 11 KB)

Additional file 4: Table S4. (XLSX 660 KB)

Additional file 5: Figure S1. Outline of experimental work flow. (PDF 30 KB)

Additional file 6: Table S5. (XLSX 15 KB)

Additional file 7: Figure S2. A 40 kb deletion removing the first exon of the ABAT gene. (PDF 32 KB)

13229_2011_36_MOESM8_ESM.PDF

Additional file 8: Figure S3. Larger genomic region of ABAT exonic deletion. Although the start 5' flank point of the deletion has not been localized, intervening probes between ABAT and A2BP1 genes indicates that the deletion terminates a significant distance away from either A2BP1 or GRIN2A, both which have been implicated in autism. (PDF 26 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Hedges, D.J., Hamilton-Nelson, K.L., Sacharow, S.J. et al. Evidence of novel fine-scale structural variation at autism spectrum disorder candidate loci. Molecular Autism 3, 2 (2012). https://doi.org/10.1186/2040-2392-3-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/2040-2392-3-2

Keywords