While numerous studies have implicated copy number variants (CNVs) in a range of neurological phenotypes, the impact relative to disease severity has been difficult to ascertain due to small sample sizes, lack of phenotypic details, and heterogeneity in platforms used for discovery. Using a customized microarray enriched for genomic hotspots, we assayed for large CNVs among 1,227 individuals with various neurological deficits including dyslexia (376), sporadic autism (350), and intellectual disability (ID) (501), as well as 337 controls. We show that the frequency of large CNVs (>1 Mbp) is significantly greater for ID–associated phenotypes compared to autism (p = 9.58×10−11, odds ratio = 4.59), dyslexia (p = 3.81×10−18, odds ratio = 14.45), or controls (p = 2.75×10−17, odds ratio = 13.71). There is a striking difference in the frequency of rare CNVs (>50 kbp) in autism (10%, p = 2.4×10−6, odds ratio = 6) or ID (16%, p = 3.55×10−12, odds ratio = 10) compared to dyslexia (2%) with essentially no difference in large CNV burden among dyslexia patients compared to controls. Rare CNVs were more likely to arise de novo (64%) in ID when compared to autism (40%) or dyslexia (0%). We observed a significantly increased large CNV burden in individuals with ID and multiple congenital anomalies (MCA) compared to ID alone (p = 0.001, odds ratio = 2.54). Our data suggest that large CNV burden positively correlates with the severity of childhood disability: ID with MCA being most severely affected and dyslexics being indistinguishable from controls. When autism without ID was considered separately, the increase in CNV burden was modest compared to controls (p = 0.07, odds ratio = 2.33).
Deletions and duplications, termed copy number variants (CNVs), have been implicated in a variety of neurodevelopmental disorders including intellectual disability (ID), autism, and schizophrenia. Our understanding of the relevance of large, rare CNVs in a range of neurodevelopmental phenotypes, varying in severity and prevalence, has been difficult because these studies were restricted to the analysis of one disorder at a time using different CNV detection platforms, insufficient sample sizes, and a lack of detailed clinical information. We tested 1,227 individuals with different neurological diseases including dyslexia, autism, and ID using the same CNV detection platform. We observed striking differences in CNV burden and inheritance characteristics among these cohorts and show that ID is the primary correlate of large CNV burden. This correlation is well illustrated by a comparison of autism patients with and without ID—where the latter show only modest increases in large CNV burden compared to controls. We also find significant depletion in the frequency of large CNVs in dyslexia compared to the other cohorts. Further studies on larger sets of individuals using high-resolution arrays and next-generation sequencing are warranted for a detailed understanding of the relative contribution of genetic variants to neurodevelopmental disorders.
Citation: Girirajan S, Brkanac Z, Coe BP, Baker C, Vives L, Vu TH, et al. (2011) Relative Burden of Large CNVs on a Range of Neurodevelopmental Phenotypes. PLoS Genet7(11): e1002334. https://doi.org/10.1371/journal.pgen.1002334
Editor: Harry T. Orr, University of Minnesota, United States of America
Received: April 11, 2011; Accepted: August 24, 2011; Published: November 10, 2011
Copyright: © 2011 Girirajan et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by grants from the National Institutes of Health 5R01 HD054562 to WHR and 1R01 HD065285 to EEE, a Simons Foundation grant to EEE, and in part by the Cancer Genomics Shared Resource at the Winship Cancer Institute, Emory University, to CSM. This study also makes use of data generated by the Wellcome Trust Case-Control Consortium. Funding for the project was provided by the Wellcome Trust under awards 076113 and 085475. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: EEE is on the scientific advisory board for Pacific Biosciences.
Recent studies have implicated large, rare CNVs in a range of neurodevelopmental disorders including intellectual disability (ID) , , autism , , schizophrenia , , bipolar disorder , , epilepsy , , and attention deficit hyperactivity disorder (ADHD) , . Several themes have emerged from these studies: first, a significant enrichment for rare CNVs in individuals with the disease compared to unaffected controls was observed, independently, for each of these disorders; second, the same recurrent CNVs are associated with different neuropsychiatric phenotypes; and third, locus heterogeneity is substantial as many distinct variants can lead to similar phenotypes.
Our understanding of the relevance of rare CNVs across a broad spectrum of neurodevelopmental disorders, varying in severity and prevalence, is limited as previous studies were restricted to the analysis of one phenotype at a time and each of such studies was performed using different CNV genotyping methodologies with distinct platform-specific biases, making comparisons difficult. We undertook a systematic analysis of 1,227 cases and 337 controls to assess the relative contribution of CNVs in three phenotypically distinct neurodevelopmental disorders. We designed a whole-genome custom microarray targeted to genomic hotspots for comparative genomic hybridization (CGH) to identify potentially pathogenic CNVs that contribute to ID, autism, and dyslexia.
We analyzed 1,227 individuals ascertained for three neurodevelopmental disorders: 376 dyslexic children with a verbal IQ (VIQ) ≥90 on the Wechsler Intelligence Scale for Children  and dyslexia defined as poor performance and IQ-performance discrepancy in one or more of a set of standardized reading measures, 350 cases with sporadic autism from the Simons Simplex Collection (SSC), and 501 cases with ID. We used 337 NIMH control individuals for comparison. Further, based on the presence or absence of ID (full-scale IQ score cutoff of 70), autism cases were divided into those with ID (n = 97) or without ID (n = 253) (see Materials and Methods). Based on the presence of multiple congenital anomalies (MCA), individuals with ID were divided into those with ID only—i.e. idiopathic ID (n = 428)—and those with ID and MCA (n = 73).
All copy number variation analyses were performed using a custom microarray with a high probe density (∼2.6 kbp) targeted to 107 genomic hotspot regions  (∼251 Mbp) and a median probe spacing of ∼36 kbp in the genomic backbone (see Materials and Methods, Table S1). We used a Hidden Markov Model (HMM)-based algorithm to identify deletions and duplications. We restricted our analysis to CNVs >50 kbp in size to reduce false positive calls and validated all relevant CNVs using a second custom designed high-density array. To empirically determine the validation rate of the array at different genomic regions, we examined 118 CNVs detected in 24 samples and confirmed 117 events (>99% accuracy, see Table S2). While we were easily able to detect smaller events in the hotspot regions, the specificity of the array restricted our CNV discovery to >50 kbp in hotspot-associated regions and to >300 kbp in regions not associated with genomic hotspots (Figure S1).
Analysis for large CNV burden in neurodevelopmental phenotypes
After quality control (QC) filtering and manual curation, we obtained 5,086 CNVs in 1,395 out of 1,564 individuals (89.2%) with high-quality array CGH data (Table 1; Datasets S1, S2, S3, S4). Using these data, we compared the CNV enrichment between the multiple cohorts tested. We found a significant excess of large CNVs (>1 Mbp) in individuals with ID (p = 2.75×10−17, odds ratio = 13.71) or autism (p = 0.012, odds ratio = 2.99) when compared to controls analyzed on the same microarray platform (Figure 1). The frequency of large CNVs among children with dyslexia was similar to controls (p = 0.64, odds ratio = 0.94), although this might indicate a lack of statistical power in our study to detect any subtle enrichment (power >0.8 to detect 4.2% increase in burden) for large CNVs in dyslexia.
(A) The figure shows the population frequency of the largest CNV (as a survivor function) in individuals with ID, autism, dyslexia, and controls. (B) Population frequency of the largest CNV is shown for ID, ID with MCA, autism with ID, autism without ID, dyslexia, and NIMH control individuals. (C) Histograms depicting deletions per individual at each size range are shown. Note that 35 NIMH control samples carried an approximately 560 kbp deletion involving PRAME on distal 22q11.2. (D) Duplications per individual at each size range are shown. The hotspot chip has higher coverage over segmental duplication regions and therefore there is an expected abundance of duplications per individual compared to deletions.
Within the neurodevelopmental disorder cohorts, a comparison showed a significantly greater large CNV burden in individuals with ID compared to autism (p = 9.58×10−11, odds ratio = 4.58) or dyslexia (p = 3.81×10−18, odds ratio = 14.45). When we partitioned the ID cohort into subsets with and without MCA, we observed a significantly increased large CNV burden in individuals with ID/MCA compared to ID alone (p = 0.001, odds ratio = 2.54). This trend was also observed when individuals with autism were separated into those with ID and without ID (Figure 1), although not statistically significant (p = 0.102, odds ratio = 2.1). When compared to controls, we noted a trend for increase in large CNV burden for autism without ID (p = 0.07, odds ratio = 2.33) as well as autism with ID (p = 0.0048, odds ratio = 4.85). In addition, a gene-based analysis showed an incremental increase in the proportion of disrupted genes and average gene density per CNV with higher estimates for the ID/MCA cohort as compared to ID alone or autism (Table 1). We also note that within the cohorts no bias towards deletions or duplications was observed in relation to phenotypic severity or variability (Tables S3, S4, S5, S6, S7). Overall, our results suggest a positive correlation of the severity of the phenotype to the size and gene density of CNVs.
Rare CNVs in dyslexia, autism, and intellectual disability phenotypes
To identify rare CNVs of likely pathogenic significance, we compared the pattern of CNVs from dyslexia, autism, ID, and NIMH control cohorts to a map developed from an expanded set of 8,329 normal individuals genotyped with Illumina microarrays and to the publicly available Database of Genomic Variants  (see Materials and Methods). We eliminated common copy number polymorphisms and CNVs from our cases if they had a reciprocal overlap of 50% or more of their length with CNVs found in these 8,329 controls. After filtering, we compared the groups. We found a significant increase of rare CNVs in individuals with autism (35/336, 10%; p = 2.4×10−6, odds ratio = 6) or ID (69/431, 16%; p = 3.55×10−12, odds ratio = 10) compared to individuals with dyslexia (6/322, 2%) (Figure 2A, Table 2). In fact, when analyzed separately, the frequency of rare CNVs in NIMH controls (6/306, 2%) was not different compared to dyslexia (p = 0.57, odds ratio = 0.94) (Table S8).
(A) The proportion of rare CNVs as a function of size is shown for NIMH controls and dyslexia, autism, and ID cohorts. To identify rare CNVs, we compared the pattern of CNVs from each of these cohorts to the CNV frequency map from 8,329 controls genotyped on Illumina arrays. (B) The proportion of de novo occurrence of CNVs among the three cohorts is shown for each size range. Note that the CNVs from the dyslexia cohort are all inherited. DNA from parents of NIMH controls was not available and hence not tested for de novo CNV frequency.
Given the high population prevalence of dyslexia , we then relaxed our selection to include events present at an allele frequency of <0.1% in controls (8/8,635) and identified four additional CNVs—i.e., a total of 10 CNV events (Table 3). The analysis of hotspot regions identified only one individual with dyslexia who carried a 15q11.2 BP1–BP2 deletion, which has previously been associated with ID , schizophrenia , , and epilepsy ; however, this deletion was also observed in 25/8,635 of our total control individuals. None of the seven deletions and three duplications detected in our dyslexia cohort mapped to candidate loci known to be associated with dyslexia .
Analysis of 336 individuals from the SSC autism cohort showed that 35 individuals (10%) carried 36 rare CNVs (680 RefSeq genes, median size = 662 kbp) and about 58% (21/36) of these CNVs mapped to genomic hotspots (Table 2). Only eight of the events (all hotspot sites) associated with genomic disorders, including 22q11.2 deletion (TBX1, DiGeorge syndrome), 17p12 duplication (PMP22, Charcot-Marie-Tooth disease), and 15q11.2q13.1 duplication (UBE3A and SNRPN) (Table S9). In addition, as reported previously , , , the autism-associated proximal 16p11.2 deletion (TBX6) was observed in approximately 1% (3/336) of all autism cases analyzed. Interestingly, one case with a de novo 16p11.2 deletion also inherited a 2 Mbp duplication 22q11.2 (TBX1) from the mother.
Among 431 cases with ID (358 cases with ID only and 73 cases with ID plus MCA), 69 individuals carried 77 rare CNVs (2,215 RefSeq genes, median size = 1.5 Mbp) that were either of known pathogenic significance or not observed in a total set of 8,635 controls, and 32% (25/77, median size = 1.42 Mbp) of these variants localized to genomic hotspot regions (Table 2). This is a significant enrichment for rare CNVs in the ID cohort compared to autism (p = 0.019, odds ratio = 1.6) or dyslexia cohorts (p = 3.55×10−12, odds ratio = 10). Interestingly, 20/77 CNVs (16 hotspot and four non-hotspot sites) mapped to a known genomic disorder site, including those associated with variable phenotypes such as 15q13.1q13.3 (CHRNA7), 16p11.2 proximal (TBX6; two cases) and distal (SH2B1) hotspots, 16p13.11 (MYH11; three cases), 17q12 (TCF2), and 3q29 (DLG1) as well as syndromic regions such as 7q11.23 (Williams syndrome), 17q21.31 (MAPT), 5q35 (Sotos syndrome), 8p23.1, 22q13 (Phelan-McDermid syndrome)  and 1p36 .
We next sought to determine whether these rare CNVs were inherited or if they arose de novo in the probands. Parental DNA samples were available to investigate inheritance for 90 out of 123 rare CNVs detected in all three disease cohorts (Table S8). In four cases, only maternal DNA was available. We find that 44/90 CNVs arose de novo and a majority (77%, 34/44) of these de novo CNVs were large (>1 Mbp). Overall, we find a greater proportion of de novo events in ID (64%, 30/47) compared to autism (40%, 14/35; p = 0.027, odds ratio = 2.6) or dyslexia (0/8; p = 0.0009, odds ratio = infinity) (Figure 2B). These data are suggestive of a general trend of increased de novo rates and CNV size with increased severity of the disorder.
Novel, rare CNVs reveal potential candidate genes
We then focused on rare CNVs involving single genes or regions of potential interest. In the dyslexia cohort, two unrelated families carried CNVs on chromosome 7q11.23 that involved the autism susceptibility candidate 2 (AUTS2, MIM# 607270). A 669 kbp duplication that included AUTS2 and WBSCR17 was transmitted from an affected father to the daughter and an approximately 84 kbp deletion was transmitted from the affected paternal grandmother through the unaffected father to the proband (Figure 3). In addition, we also identified a 354 kbp deletion encompassing AUTS2 in one individual with idiopathic ID, pervasive developmental delay, partial epilepsy, and left hemihypertrophy. An approximately 1.2 Mbp deletion encompassing the eyes shut drosophila homolog gene (EYS, MIM# 612424) on chromosome 6q12 was detected in an affected proband and several unaffected family members. Although autosomal recessive single-nucleotide mutations in EYS have been reported in patients with retinitis pigmentosa , , the role of heterozygous microdeletions involving this gene is unknown.
(A) The father has features of dyslexia with a verbal IQ (VIQ) of 122, WATTa 93, WRAT3spb 83 and WIATspc 94. The daughter's scores are VIQ 122, WIDd 86, WATT 91, WRAT3sp 90, and WIATsp 94. (B) In this pedigree the 84 kbp deletion within AUTS2 is transmitted to the proband (VIQ 122, WATT 97, WRAT3sp 96, WIATsp 94) from the affected paternal grandmother (VIQ 118, WATT 88) through the unaffected father. (C) A deletion within IMMP2L is shown for this family. The deletion is transmitted to the proband (VIQ 111, WID 83, WATT 82, WRAT3sp 85, WIATsp 81) from his affected father (VIQ 84, WRAT3sp 68, WIAT-2sp 66). Interestingly, IMMP2L variants have been associated with Tourette syndrome, ADHD, and autism. (D) A 1.2 Mbp deletion within EYS is shown in several family members of this large pedigree. While the proband (VIQ 107, WID 56, WATT 81, WRAT3sp 87, WIATsp 87) and his affected brother (VIQ 101, WID 66, WATT 82, WRAT3sp 78, WIATsp 85) carried the deletion, so did many other unaffected relatives, including the father. Although no inference can be drawn for its role in dyslexia, as the deletion does not segregate with the phenotype, recessive mutations in EYS have been associated with retinitis pigmentosa. aWATT- WRMT-R Woodcock Reading Mastery Test – Revised; Word Attack subtest . A measure of untimed reading of single non-words. bWRAT3sp - Wide Range Achievement Tests – Third Addition; Spelling subtest . Spelling of single words from dictation in writing. cWIAT(2)sp - Wechsler Individual Achievement Test (2nd edition); Spelling subtest . Spelling of single words from dictation in writing. dWID - WRMT-R Woodcock Reading Mastery Test – Revised; Word Identification subtest . A measure of untimed reading of single words.
We also identified a 471 kbp deletion involving IMMP2L inherited by the proband from the affected mother (Figure 3). Deletions involving IMMP2L have been associated with ADHD , autism , and Tourette syndrome . Recently, Pagnamenta and colleagues also reported a 594 kbp IMMP2L-DOCK4 deletion resulting in a fusion transcript and an intragenic DOCK4 deletion segregating with dyslexia . Our results are best interpreted within the context of candidate gene identification in dyslexia. Although at least nine chromosomal loci are associated with dyslexia, for two of these loci the candidate genes were identified on the basis of a rare balanced chromosomal translocations disrupting ROBO1  and DYX1C1/EKN1 , . More recently, a Danish Cytogenetic Registry study of all cases with chromosomal translocations identified additional novel dyslexia candidate genes affirming the value of rare structural variants in understanding the genetics of dyslexia . Our study is the first to systematically characterize rare CNVs in dyslexia and thus evaluate the contribution of rare deletions and duplications to this common genetic disorder.
Within the autism cohort, several novel deletions and duplications involving neurologically-relevant genes were identified. A 5 Mbp de novo deletion involving FOXP1 on chromosome 3p14.1 was identified in an individual with features of idiopathic autism (full-scale IQ = 75). An additional 6.6 Mbp de novo deletion overlapping FOXP1 was also identified in an individual with idiopathic ID (Figure 4). A review of the DECIPHER database revealed a similar-sized deletion disrupting FOXP1 in an individual with developmental delay, sensorineural deafness, hypotonia, club foot, and dislocation of hip. Recently, FOXP1 was implicated in autism, ID, and language impairment , , . It is believed that FOXP1interacts with FOXP2 and CNTNAP2, both implicated in speech disorders and autism , . The overlapping 1.16 Mbp region of the deletion common to both autism and ID indicates a potential involvement of FOXP1 in pathways related to both of these disorders. Other variants involving functionally relevant genes include 7q36.2 deletion and duplication (DPP6), 17q23.3 duplication (SCN4A), and 17q21.32 duplication (WNT3 and WNT9B).
Two deletions (5 Mbp and 6.6 Mbp) are shown intersecting at a common region of 1.16 Mbp containing FOXP1. Note that the deletion in the autism individual also covers ROBO2 and CNTN3.
Our analysis of the ID cohort was enriched for singleton events often involving genes related to developmental or neurological functions including SYNPR, GABRA, AUTS2, FOXP1, FKBP6, COBL, and FMR1. However, within the same cohort we were also able to detect novel overlapping deletions (3.5 Mbp and 9 Mbp) on chromosome 9p24 in two unrelated individuals (Figure 5A). Both cases exhibited clinical features of ID and Pervasive Developmental Delay-Not Otherwise Specified. The distal breakpoints of these deletions map to segmental duplications while the proximal end maps within a high density of repeat elements. A survey of this region in the DECIPHER database  revealed about 15 cases with overlapping deletions. Variable clinical presentations and heterogeneity of deletion breakpoints preclude further genotype-phenotype correlation studies for this region (Figure S2). We also identified a nonrecurrent 6q16 deletion (chr6: 100,383,567-103,310,184) that potentially narrows the critical region for this recently described Prader-Willi-like syndrome  to approximately 2.9 Mbp. The refined critical region contains only five genes including the obesity-associated SIM1  and the autism-associated GRIK2  (Figure 5B). About 70% of children with 6q16 deletion manifest obesity ; however, our case with the smaller deletion, encompassing SIM1, showed no evidence of obesity at 10 years of age (Table S10).
(A) Overlapping deletions on chromosome 9p24 are shown. Deletions of this region containing DMRT1 and DMRT3 have also been associated with urogenital abnormalites and sex reversal. (B) A ∼3 Mbp deletion on 6q16 is shown encompassing SIM1. Larger deletions of this region were previously reported and have been associated with obesity. Orange bars denote deletions and gray regions are non-deleted regions from previous studies . Taken together with other published studies, this deletion narrows down the critical region (dotted box) to about 2.9 Mbp.
Multiple large CNVs associate with phenotypic severity
We find that 8/69 (11.6%) cases in the ID cohort carried more than one large, rare CNV and all of these individuals presented with severe clinical features (Tables S6, S7, S8). A striking difference (p = 0.008, odds ratio = 11.2) in multiple CNV rates was also observed when the ID cohort was divided into those with severe MCA (44%) and those with idiopathic ID (6.7%). Notable examples are co-occurrences of a 3.4 Mbp 16p13.11 duplication and a 3.3 Mbp deletion on chromosome 4q25 involving PITX2 in a case with features of Rieger syndrome  (Figure 6) and a 17p13.3 deletion (YWHAE, Miller-Dieker syndrome) and 3q29 duplication (DLG1) in a child with cryptorchidism, ventricular septal defect, and seizures. The 3q29 duplication is a recurrent interstitial rearrangement ,  potentially mediated by flanking segmental duplications of high sequence identity (27 kbp size, 96% identity). The 17p13.3 deletion is a previously reported nonrecurrent rearrangement associated with Miller-Dieker syndrome , . In contrast, only 1/35 (2.9%) autism cases and none of dyslexia individuals carried another large CNV. This observation suggests that the severity of the phenotypes can be influenced by more than one large, rare CNV co-occurring in the same individual. During this analysis we considered the possibility of a derivative chromosome representing an unbalanced translocation possibly creating the impression of multiple CNVs in our cases. We carefully reviewed available chromosomal analysis data (G-banded karyotyping or FISH) for each of the individuals with two hits reported in our study. We did find one case with two hits where apparent CNVs represent a derivative chromosome inherited from a balanced translocation carrier parent (Table S8D).
A 3.3 Mbp deletion containing PITX2 as well as the 3.4 Mbp paternally inherited 16p13.33 duplication is shown for an individual with ID plus MCA. This individual has features of Rieger syndrome including visual defects, mild hypotonia, right congenital glaucoma, left microophalmia, and anterior segment dysgenesis. Other features include cleft uvula, hypodontia and conical teeth, hyperplasia of frenulum of tongue, midface hypoplasia, strabismus, and deafness.
Initial discoveries of significant enrichment of rare CNVs for ID and autism led to testing the CNV basis for other behavioral and neurodevelopmental disorders of varying population frequency and severity, such as schizophrenia, ADHD, epilepsy, bipolar disorder, and Tourette syndrome. However, comparisons between these studies have been difficult due to differences in study design, insufficient sample sizes, and lack of detailed phenotype information. In this study, we compared 1,564 individuals (cases and controls) on a single platform of relatively modest density with the same type of detection bias. We utilized the duplication architecture of the human genome to custom design a DNA oligonucleotide microarray enriched for genomic hotspots, i.e., regions flanked by high-identity segmental duplications. This array has an advantage over several other commercial arrays in that there is a 25-fold enrichment for recurrent events in the genomic hotspots compared to the rest of the genome . Therefore, fewer samples are required to identify several unrelated individuals with the same pathogenic mutation. We find that our array has a comparable diagnostic yield of 16% for the ID cohort compared to other clinical chromosomal microarray studies reported in the literature (Figure S3).
In strong agreement with previous studies, our data suggest that multiple, rare CNVs contribute to the etiology of autism and ID. In contrast, we find no increase in large pathogenic CNVs in individuals with dyslexia compared to controls. Notwithstanding, our analysis revealed novel regions of potential relevance to the etiology of dyslexia. Two unrelated children (2/322, 0.6%) with dyslexia carried CNVs encompassing AUTS2, both inherited from a parent. While the phenotype of dyslexia segregated with the AUTS2 duplication in the first family (Figure 3A), in the second family the deletion was inherited from affected grandmother through unaffected father (Figure 3B). This could be due to a phenomenon described as “compensation”, where some adults that reported difficulties with reading in childhood no longer evidence signs of dyslexia . Previous studies of de novo chromosomal translocations and inversions identified breakpoints within AUTS2 in individuals with autism and/or ID phenotypes , , , , . More recently, unique AUTS2 deletions and duplications were observed in Juvenile Myoclonic Epilepsy  and ADHD , . It is interesting to note that ADHD and dyslexia are frequently comorbid and may have shared genetic risk factors , , . When our study is taken together with recent CNV studies of ADHD , , AUTS2 CNVs were observed in 5/2,306 combined cases and 3/46,947 unscreened controls (p = 1.12×10−5, odds ratio = 33.9), indicating that AUTS2 might have an important role in pathways related to cognition. While the function of AUTS2 is still unclear, it is strongly expressed in fetal and adult brains, particularly in the frontal, parietal, and temporal lobes . Interestingly, AUTS2 and the 7q11.2 region were identified as having the strongest statistical signal for positive selection in early modern humans as compared to the Neanderthal genome , indicating that AUTS2 might be important for a specialized human function such as cognition.
The CNV profile we observed in individuals with dyslexia was essentially the same as that in control individuals. This is not surprising if we take into consideration that all the subjects in our dyslexia sample had a VIQ above the 25%ile and the mean VIQ of the cohort was 110 (2/3 standard deviations above the general population mean), and given that we have shown that the CNV profile correlates with the severity of ID. The genes involved in dyslexia are likely to affect more specialized cognitive functions, may not adversely affect general intelligence, and may be more amenable to discovery with high-density arrays capable of detecting single gene or single exon CNVs or SNP microarrays that can leverage SNP allele frequency information in addition to signal intensity. In addition, all hybridization-based platforms fail to detect copy number neutral changes, such as balanced chromosomal rearrangements and inversions. This is particularly germane to dyslexia where a large number of candidate genes have been identified through mapping of translocation breakpoints , .
A comparison of rare de novo CNV rates for autism shows that our estimates (4%, 14/336) fall within a range of 4–10% reported previously by other large-scale, high-density array studies , , , . This suggests that no platform-specific bias exists for large variants and also that the contribution of large CNVs is consistent across all studies for autism. We find a significantly greater enrichment for large CNVs, higher de novo rates, and a higher frequency of two rare CNV hits in individuals with ID-associated phenotypes compared to autism or dyslexia. This observation is exemplified by the fact that individuals with autism with ID have more large CNVs than those with autism only. We also find a significant difference between individuals with autism versus those with dyslexia. Sanders and colleagues recently analyzed 1,124 SSC families affected with autism spectrum disorder (ASD). Using stepwise linear models, they evaluated the relationship between intellectual functioning, sex, and the number of genes within rare, de novo CNVs. While the number of genes affected correlated with the size of the de novo CNV, the authors did not find a strong correlation of the Autism Diagnostic Observation Schedule (ADOS) combined severity score (p = 0.25, R2 = 0.005) or full-scale IQ (p = 0.02, R2 = 0.08) with the size of the CNV. In contrast, we considered all large CNVs (common and rare, de novo and transmitted) identified in a relatively smaller sample size and essentially bifurcated the autism cohort using a full-scale IQ score cutoff of 70. There was also a greater enrichment of two hits in the ID cohort (11.6%) compared to the autism cohort (2.8%). In fact, one individual carrying a 16p11.2 deletion with autism and features of ID also has a maternally inherited 22q11.2 duplication (TBX1) providing further evidence for the two-hit hypothesis we previously proposed for severe developmental delay . Further, the frequency of two hits was even more striking when only individuals with ID/MCA were considered (44%), albeit the number of cases is few. We believe these data provide support for an incremental effect of CNV size and number on the severity of phenotypic outcome.
Our experimental design is biased towards interrogating hotspot regions in the human genome. A comparison to recently reported studies ,  suggests that the majority of false-negative calls will reside within non-hotspot regions due to a lack of probe coverage (<10 probes). While the detection power of our array increases with the size of the variant, we would certainly miss smaller and intragenic CNVs, for example in autism candidate genes such as NRXN1 , , CACNA1C, SLC4A10, MAGI1 , SYNGAP1, DLGAP2 , NLGN1, ASTN2 , and exonic copy number variants in ASPM, DPP10, CNTNAP2, A2BP1, PCDH9 , and PTCHD1 . While we find no excess of large CNVs in dyslexia, there is still the possibility that large CNVs are relevant in some familial cases of the disease as well as occasional sporadic cases. Further studies are warranted for a more detailed analysis of all the three neurodevelopmental cohorts using high-resolution arrays and next-generation exome and/or whole-genome sequencing. While it can be difficult to compare data derived from different microarrays, there is value to multiple array platforms and cross-platform validation. The depositing of the resulting data into publicly available databases will facilitate the continued elucidation of recurring clinically significant CNV and genotype-phenotype correlations.
Materials and Methods
Patients from each of study cohort were recruited after appropriate human subjects approval and informed consent. Informed consent was also obtained to publish photographs.
Patient ascertainment based on severity of phenotypes
DNA samples were obtained from cases ascertained for three neurodevelopmental disorders of varying severity: (1) ID/developmental delay and MCA, (2) dyslexia or reading impairment, and (3) idiopathic autism. We defined severity of clinical features based on presence or absence of ID (IQ<70) for the autism group and congenital malformation for the ID group. Our dyslexia cohort had no ID or congenital malformation cases; as an IQ≤90 and the presence of congenital malformations were exclusion criteria. Individuals with idiopathic autism were partitioned into those with autism and ID (IQ<70) and those without ID (IQ>70). For the ID cohort, those individuals with brain malformations, gross craniofacial dysmorphology, cardiac defects, and neurological deficits were separated into an ID plus MCA (ID/MCA) group. Thus, in the order of severity, the ID/MCA cohort is considered the most severe, followed by ID only, autism with ID, autism without ID, dyslexia, and normal controls. However, we note that although the individuals with dyslexia do not have ID, they have severe impairments in core phonological measures leading to significantly reduced reading abilities despite normal IQ (IQ≥90). Detailed descriptions of each of the cohorts are given below.
Ascertainment of individuals with dyslexia or reading disability
For the dyslexia subject set, children were considered eligible for the study if they met researcher-defined criteria based on test scores from a standardized battery of tests. DNA samples were obtained from two cohorts. The first cohort included probands aged 6 to 16 from 198 families who were initially ascertained at the University of Washington (UW) multidisciplinary Learning Disability Center (UWLDC) under protocols approved by the UW Institutional Review Board. For the UWLDC cohort, probands were required to have a prorated VIQ at or above 90 (≥25%ile) on the Wechsler Intelligence Scale for Children – 3rd edition , with performance below the age-specific population mean and at least one standard deviation below the VIQ on one or more out of 10 research measures of reading, writing, or spelling. As a group, on average, probands met the impairment criteria between 6 and 7 measures. As expected by ascertainment requirements, the average VIQ of probands in this cohort was 110 (≥75%ile) , . Siblings older than 6.5 years were invited to participate, and additional family members were added using a sequential sampling strategy to extend pedigrees through family members with the most extreme impairment values on the same 10 research measures. Detailed recruitment and evaluation procedures for the UWLDC cohort were described earlier ,  (see Table S4A).
For the second cohort, 178 children aged 5 to 12 were recruited from a special K-6 school for students with dyslexia or via their direct relatives in the Atlanta area (The Schenck School, Atlanta, GA). For this cohort, children were required to have a psychological battery of tests completed by a licensed psychologist and usually have a diagnosis of a reading disability. Based on strong verbal comprehension score, perceptual reasoning score, Peabody picture vocabulary test, or other cognitive tests that measure intelligence, these children have average to above-average intelligence. Both cohorts were composed of individuals with >90% Caucasian ethnicity with an approximately equal number of males and females. Except for ADHD, children with other psychiatric and neurological disorders, moderate to severe receptive language disorders, developmental disabilities, or other conditions known to affect cognition were excluded based on parental questionnaire. Clinical details are shown in Table S4B.
Ascertainment of individuals with features of autism with or without ID
For the autism cohort, families were identified through the SSC (www.sfari.org) . The Simons Foundation-funded SSC includes families with no more than one child with autism ascertained through 12 data collection sites across North America. Of the 350 individuals included in this study 297 (85%) are of Caucasian ethnicity. Inclusion criteria in the collection requires that the child with autism meet ASD criteria on the ADOS , on the Autism Diagnostic Interview, Revised (ADI-R) , and meet expert clinical judgment. Nonverbal IQ estimate must also be greater than 35. Children with significant hearing, vision, or motor problems, significant birth complications (e.g. extended NICU stay), or with a diagnosis of ASD-related disorders, such as Fragile X, were excluded. Children with a relative (up to third degree) with ASD or sibling who showed ASD-related symptoms were also excluded. Diagnostic evaluations, cognitive assessment, and phenotypic characterization were conducted at each site with data collection, data entry, and data validation methods standardized across sites to ensure reliability of sample collection. We further partitioned the autism cohort into those associated with ID (average full scale IQ = 49) consisting of 97 cases (73 males, 24 females; median age, 12 years) and those without ID (average full scale IQ = 98.9) comprising 253 cases (228 males, 25 females; median age, 11 years and 11 months). Clinical details are shown in Table S5.
Ascertainment of individuals with intellectual disability with or without congenital malformation
The idiopathic ID cohort was selected from individuals admitted to the IRCCS Associazione Oasi Maria Santissima and screened for ID according to the Diagnostic and Statistical Manual of Mental Disorders-IV-Text Revision (DSM-IV-TR) criteria. This cohort consists of 428 cases (153 females, 275 males; median age, 15 years) of Caucasian ethnicity with idiopathic ID and previously excluded for common causes of ID, including Fragile X syndrome, trisomies 21 and 13. In addition, classical genomic syndromes such as Smith-Magenis, DiGeorge, Prader-Willi/Angelman, and Williams syndromes, if recognized by clinical evaluation, were followed up for confirmation using targeted multiplex ligation-dependent probe amplification and excluded. We note that cases with phenotypic variability that escape clinical detection might not have been excluded. Typically, idiopathic cases of ID with no classical constellation of clinical features suggestive of a known disorder or those with mild to moderate ID without significant congenital malformation were included in this cohort. Clinical details are shown in Table S6.
Individuals with features of ID with MCA not necessarily assigned to a specific syndrome were evaluated and recruited at the University of Torino. This cohort consists of 73 individuals (32 females and 41 males) of Caucasian ethnicity with a median age of 2 years at diagnosis. Clinical features of these individuals included brain malformations, craniofacial dysmorphology, and neurological deficits along with variable ID (Table S7). Informed consent was obtained from all the subjects included in both the studies.
Ascertainment of normal controls
The control cohort consisted of 337 DNA samples obtained from the Rutgers University Cell and DNA Repository (www.rucdr.org). These individuals were ascertained by the NIMH Genetics Initiative  through an online self-report based on the Composite International Diagnostic Instrument Short-Form (CIDI-SF)  and screened specifically for eight mental health disorders, including major depression, bipolar disorder, and psychosis, but were not screened for dyslexia and therefore not ideal for such comparisons. Those who did not meet DSM-IV criteria for major depression, denied a history of bipolar disorder or psychosis, and reported exclusively European origins were included , .
Additionally, CNV data from 8,329 additional cell line and blood-derived controls were used to assess the frequency of our putative pathogenic CNVs in a larger population of neurologically normal individuals. These data were derived primarily from genome-wide association studies of non-neurological phenotypes. Although these data were not ascertained specifically for neurological disorders, they consist of adult individuals providing informed consent. Specifically, datasets from the following sources were included in our analysis: Human Genome Diversity Project , ; National Institute of Neurological Disorders and Stroke (NINDS) (dbGaP accession no. phs000089) , ; Pharmacogenomics and Risk of Cardiovascular Disease (PARC/PARC2) , ; parents of asthmatic children courtesy of Stephanie London ; Fred Hutchinson Cancer Research Center (prerelease data provided courtesy of Aaron Aragaki, Charles Kooperberg, and Rebecca Jackson as part of an ongoing genome-wide association study to identify genetic components of hip fracture in the Women's Health Initiative); InCHIANTI (data provided by InCHIANTI study of aging, www.inchiantistudy.net) , ; and the Wellcome Trust Case Control Consortium phase 2 (National Blood Service) . All samples were genotyped on Illumina arrays using methodology described previously   and either natively processed in hg18 or re-mapped after CNV calling (NINDS and PARC) to hg18 using the UCSC LiftOver tool (http://genome.ucsc.edu).
Array CGH and analysis
We designed custom targeted hotspot v1.0 arrays comprised of 135,000 probes (by Roche NimbleGen) with higher density probe coverage (median probe spacing 2.6 kbp) in the genomic hotspots (regions flanked by segmental duplications) and a lower probe density in the genomic backbone (median probe spacing 36 kbp). All microarray hybridization experiments were performed as described previously , using a single unaffected male (GM15724 from Coriell) as reference. All validation experiments were performed using two custom array designs: (1) a custom targeted 4×180 K Agilent chip with median probe spacing of 2 kbp in the genomic hotspots and whole-genome backbone coverage of one probe every 36 kbp (Agilent Technologies) and (2) a custom targeted 3×720 K NimbleGen or 2×400 K Agilent chip with median probe spacing of 500 bp in the genomic hotspots and probe spacing of 14 kbp in the genomic backbone.
All arrays were analyzed by mapping probe coordinates to the human genome assembly Build 36 (hg18). Using chromosome-specific means and standard deviations, normalized log intensity ratios for each sample were transformed into z-scores. These z-scores were then classified as “increased”, “normal”, or “decreased” in copy number using a three-state HMM. The HMM was applied using HMMSeg . For each sample, HMM state assignments of probes were merged into segments if consecutive probes of the same state less than 50 kbp apart. If two segments of the same state were separated by an intervening sequence of ≤5 probes and ≤10 kbp, both segments and intervening sequence were called as a single variant. Further, we employed stringent QC measures and empirically estimated post-HMM filtering thresholds (absolute z-score >1.5 and >10 probes) to increase the specificity of our experiments. With these filtering criteria, we were able to thoroughly scan HMM outputs for CNV events and manually check the validity of each call by examining the normalized log intensity ratios across a chromosome. For the Agilent arrays, data analysis was performed following feature extraction using DNA analytics with ADM-2 setting according to the manufacturer's instructions. All CNVs calls were visually inspected in the UCSC genome browser.
First, we carried out validation on 24 samples from the developmental delay cohort and confirmed 117/118 HMM-inferred calls with a validation rate of 99.15%. Next, we validated 84 CNVs from an independent set of cases both validated using fluorescence in situ hybridization (FISH) and locus-specific custom high-density arrays  (Girirajan and Eichler, unpublished). We also validated all 44 calls from the autism cohort, 22 calls from the dyslexia cohort, and 78 calls from the developmental delay cohort. In addition, in an analysis of 517 individuals with epilepsy using this array design, 61/63 CNVs were validated on a different array platform .
Probe densities and CNV detection threshold of hotspot v1 chip. Although the median density of the chip was designed to be approximately 2.6 kbp in the genomic hotspots and 36 kbp in genomic backbone, the limitations of the chip design (probe assignment restricted to only up to five mismatches) precluded uniform distribution of the probes throughout the genome. Therefore, the actual probe density varied across regions of the human genome. (a) The plot shows the size-wise distribution of CNVs and the density of array probes targeted to the genomic hotspots and non-hotspot regions. Note that non-hotspot regions contain two different probe densities (probe spacing of 20,000 bp and about 30,000 bp labeled a & b) and the hotspots (shaded) are covered every 10,000 bp. (b) The histogram shows the CNV detection threshold for the hotspot v1 chip. The data represent all the CNV calls obtained from analyzing cases with intellectual disability (ID). The number of CNVs detected in the aggregate at different size thresholds is shown on the Y-axis. We utilized a threshold of >10 probes, >1.5 z-score, and >50 kbp for CNV detection analysis. While we were easily able to detect events >50 kbp in the hotspot regions, we were only able to call variants ranging from 150 kbp (hotspot-associated CNVs) and 300 kbp onwards (mostly non-hotspot CNVs) with confidence (i.e., able to validate) for the non-hotspot regions.
Comparison of 9p24 deletions identified in our study to DECIPHER database. The figure shows genome browser snapshot of the 9p24 region with red and blue bars denoting deletions and duplications respectively in the DECIPHER database. Major clinical features observed in these patients are also shown. Please note that the black bars at the top denote CNVs identified in our study. PDD-NOS, Pervasive developmental delay-not otherwise specified; ID/DD, intellectual disability/developmental delay.
Diagnostic yield of different microarray reports from literature. Histograms show the number of rare CNVs (usually disease-associated) observed under different diagnostic centers. Data is shown for sample sizes >50. Data obtained from Table 2 of Miller et al., AJHG.
Chromosomal regions and probes targeted in the hotspot chip. List of all regions targeted in the hotspot chip. The list is derived from the original curation of segmental duplication regions and genomic hotspots from Bailey et al, 2002 .
Confirmation of CNVs arrays using custom high-density arrays. Validation of CNVs identified using NimbleGen hotspotv1 arrays (12×135 K) using a higher density 3×720K NimbleGen or 2×400K Agilent arrays.
Global analysis for CNVs in neurodevelopmental disorders. Summary and characteristics of deletions and duplications identified in the four cohorts studied.
(A) Characteristics of dyslexia cases from UW. Dyslexia measurements and scores of children tested at the UW. WATT- WRMT-R Woodcock Reading Mastery Test – Revised; Word Attack subtest. A measure of untimed reading of single non-words. WRAT3sp - Wide Range Achievement Tests – Third Addition; Spelling subtest. Spelling of single words from dictation in writing. WIAT(2)sp - Wechsler Individual Achievement Test (2nd edition); Spelling subtest. Spelling of single words from dictation in writing. WID - WRMT-R Woodcock Reading Mastery Test – Revised; Word Identification subtest. A measure of untimed reading of single words. (B) Characteristics of dyslexia cases recruited from Atlanta. Gender, phenotype and age of children recruited from the Atlanta collection is shown.
Clinical details of cases with autism (both with and without ID). Clinical features of individuals recruited through the Simons Simplex Collection (n = 350) are shown.
Clinical details of cases with idiopathic ID. Clinical features of individuals with ID are shown. These individuals were recruited and evaluated at the IRCCS Associazione Oasi Maria Santissima, Troina.
Clinical features of cases with ID plus MCA. Presenting clinical features, additional malformations, and preliminary cytogenetic and genetic evaluations performed for these cases are shown.
(A) Comparison of rare CNV rates in the cohorts studied. (B) Rare CNVs in dyslexia, autism, and ID. (C) Inheritance of rare CNVs in the disease cohorts. (D) Individuals with two rare copy number variants (two hits).
Genomic disorders identified in 1,227 cases with ID, autism, and ID/MCA. Genomic disorders are defined as copy number variants that are previously identified to be associated significantly in individuals with disease compared to controls. These CNVs can either map within genomic hotspot (HS) or non-hotspot regions (non-HS).
We thank Catarina Campbell, Mark Matsushita, John Wolff, Arthur Ko, Kenneth Mark, Tonia Brown, Nik Krumm, Noelani Laycock, and David A. Hughes for valuable discussions and technical assistance. A full list of the investigators who contributed to the generation of the data is available from www.wtccc.org.uk. EEE is an investigator with the Howard Hughes Medical Institute.
Conceived and designed the experiments: SG EEE. Performed the experiments: SG CB LV THV NS. Analyzed the data: SG. Contributed reagents/materials/analysis tools: SG EEE BPC ZB STW WHR CSM CR MF GBF MS RB. Wrote the paper: SG EEE.
- 1. Sharp AJ, Hansen S, Selzer RR, Cheng Z, Regan R, et al. (2006) Discovery of previously unidentified genomic disorders from the duplication architecture of the human genome. Nat Genet 38: 1038–1042.AJ SharpS. HansenRR SelzerZ. ChengR. Regan2006Discovery of previously unidentified genomic disorders from the duplication architecture of the human genome.Nat Genet3810381042
- 2. de Vries BB, Pfundt R, Leisink M, Koolen DA, Vissers LE, et al. (2005) Diagnostic genome profiling in mental retardation. Am J Hum Genet 77: 606–616.BB de VriesR. PfundtM. LeisinkDA KoolenLE Vissers2005Diagnostic genome profiling in mental retardation.Am J Hum Genet77606616
- 3. Marshall CR, Noor A, Vincent JB, Lionel AC, Feuk L, et al. (2008) Structural variation of chromosomes in autism spectrum disorder. Am J Hum Genet 82: 477–488.CR MarshallA. NoorJB VincentAC LionelL. Feuk2008Structural variation of chromosomes in autism spectrum disorder.Am J Hum Genet82477488
- 4. Sebat J, Lakshmi B, Malhotra D, Troge J, Lese-Martin C, et al. (2007) Strong association of de novo copy number mutations with autism. Science 316: 445–449.J. SebatB. LakshmiD. MalhotraJ. TrogeC. Lese-Martin2007Strong association of de novo copy number mutations with autism.Science316445449
- 5. International Schizophrenia Consortium (2008) Rare chromosomal deletions and duplications increase risk of schizophrenia. Nature 455: 237–241.International Schizophrenia Consortium2008Rare chromosomal deletions and duplications increase risk of schizophrenia.Nature455237241
- 6. Walsh T, McClellan JM, McCarthy SE, Addington AM, Pierce SB, et al. (2008) Rare structural variants disrupt multiple genes in neurodevelopmental pathways in schizophrenia. Science 320: 539–543.T. WalshJM McClellanSE McCarthyAM AddingtonSB Pierce2008Rare structural variants disrupt multiple genes in neurodevelopmental pathways in schizophrenia.Science320539543
- 7. Craddock N, Hurles ME, Cardin N, Pearson RD, Plagnol V, et al. (2010) Genome-wide association study of CNVs in 16,000 cases of eight common diseases and 3,000 shared controls. Nature 464: 713–720.N. CraddockME HurlesN. CardinRD PearsonV. Plagnol2010Genome-wide association study of CNVs in 16,000 cases of eight common diseases and 3,000 shared controls.Nature464713720
- 8. Grozeva D, Kirov G, Ivanov D, Jones IR, Jones L, et al. (2010) Rare copy number variants: a point of rarity in genetic risk for bipolar disorder and schizophrenia. Arch Gen Psychiatry 67: 318–327.D. GrozevaG. KirovD. IvanovIR JonesL. Jones2010Rare copy number variants: a point of rarity in genetic risk for bipolar disorder and schizophrenia.Arch Gen Psychiatry67318327
- 9. Helbig I, Mefford HC, Sharp AJ, Guipponi M, Fichera M, et al. (2009) 15q13.3 microdeletions increase risk of idiopathic generalized epilepsy. Nat Genet 41: 160–162.I. HelbigHC MeffordAJ SharpM. GuipponiM. Fichera200915q13.3 microdeletions increase risk of idiopathic generalized epilepsy.Nat Genet41160162
- 10. Mefford HC, Muhle H, Ostertag P, von Spiczak S, Buysse K, et al. (2010) Genome-wide copy number variation in epilepsy: novel susceptibility loci in idiopathic generalized and focal epilepsies. PLoS Genet 6: e1000962.HC MeffordH. MuhleP. OstertagS. von SpiczakK. Buysse2010Genome-wide copy number variation in epilepsy: novel susceptibility loci in idiopathic generalized and focal epilepsies.PLoS Genet6e1000962
- 11. Williams NM, Zaharieva I, Martin A, Langley K, Mantripragada K, et al. (2010) Rare chromosomal deletions and duplications in attention-deficit hyperactivity disorder: a genome-wide analysis. Lancet 376: 1401–1408.NM WilliamsI. ZaharievaA. MartinK. LangleyK. Mantripragada2010Rare chromosomal deletions and duplications in attention-deficit hyperactivity disorder: a genome-wide analysis.Lancet37614011408
- 12. Elia J, Gai X, Xie HM, Perin JC, Geiger E, et al. (2010) Rare structural variants found in attention-deficit hyperactivity disorder are preferentially associated with neurodevelopmental genes. Mol Psychiatry 15: 637–646.J. EliaX. GaiHM XieJC PerinE. Geiger2010Rare structural variants found in attention-deficit hyperactivity disorder are preferentially associated with neurodevelopmental genes.Mol Psychiatry15637646
- 13. Wechsler D (1991) Wechsler intelligence scale for children - third edition (WISC-III). San Antonio: The Psychological Corporation. D. Wechsler1991Wechsler intelligence scale for children - third edition (WISC-III)San AntonioThe Psychological Corporation
- 14. Bailey JA, Gu Z, Clark RA, Reinert K, Samonte RV, et al. (2002) Recent segmental duplications in the human genome. Science 297: 1003–1007.JA BaileyZ. GuRA ClarkK. ReinertRV Samonte2002Recent segmental duplications in the human genome.Science29710031007
- 15. Zhang J, Feuk L, Duggan GE, Khaja R, Scherer SW (2006) Development of bioinformatics resources for display and analysis of copy number and other structural variants in the human genome. Cytogenet Genome Res 115: 205–214.J. ZhangL. FeukGE DugganR. KhajaSW Scherer2006Development of bioinformatics resources for display and analysis of copy number and other structural variants in the human genome.Cytogenet Genome Res115205214
- 16. Pennington BF (1990) The genetics of dyslexia. J Child Psychol Psychiatry 31: 193–201.BF Pennington1990The genetics of dyslexia.J Child Psychol Psychiatry31193201
- 17. Mefford HC, Cooper GM, Zerr T, Smith JD, Baker C, et al. (2009) A method for rapid, targeted CNV genotyping identifies rare variants associated with neurocognitive disease. Genome Res. HC MeffordGM CooperT. ZerrJD SmithC. Baker2009A method for rapid, targeted CNV genotyping identifies rare variants associated with neurocognitive disease.Genome Res
- 18. Kirov G, Grozeva D, Norton N, Ivanov D, Mantripragada KK, et al. (2009) Support for the involvement of large copy number variants in the pathogenesis of schizophrenia. Hum Mol Genet 18: 1497–1503.G. KirovD. GrozevaN. NortonD. IvanovKK Mantripragada2009Support for the involvement of large copy number variants in the pathogenesis of schizophrenia.Hum Mol Genet1814971503
- 19. Stefansson H, Rujescu D, Cichon S, Pietilainen OP, Ingason A, et al. (2008) Large recurrent microdeletions associated with schizophrenia. Nature 455: 232–236.H. StefanssonD. RujescuS. CichonOP PietilainenA. Ingason2008Large recurrent microdeletions associated with schizophrenia.Nature455232236
- 20. de Kovel CG, Trucks H, Helbig I, Mefford HC, Baker C, et al. (2010) Recurrent microdeletions at 15q11.2 and 16p13.11 predispose to idiopathic generalized epilepsies. Brain 133: 23–32.CG de KovelH. TrucksI. HelbigHC MeffordC. Baker2010Recurrent microdeletions at 15q11.2 and 16p13.11 predispose to idiopathic generalized epilepsies.Brain1332332
- 21. Paracchini S, Scerri T, Monaco AP (2007) The genetic lexicon of dyslexia. Annu Rev Genomics Hum Genet 8: 57–79.S. ParacchiniT. ScerriAP Monaco2007The genetic lexicon of dyslexia.Annu Rev Genomics Hum Genet85779
- 22. Kumar RA, KaraMohamed S, Sudi J, Conrad DF, Brune C, et al. (2008) Recurrent 16p11.2 microdeletions in autism. Hum Mol Genet 17: 628–638.RA KumarS. KaraMohamedJ. SudiDF ConradC. Brune2008Recurrent 16p11.2 microdeletions in autism.Hum Mol Genet17628638
- 23. Weiss LA, Shen Y, Korn JM, Arking DE, Miller DT, et al. (2008) Association between microdeletion and microduplication at 16p11.2 and autism. N Engl J Med 358: 667–675.LA WeissY. ShenJM KornDE ArkingDT Miller2008Association between microdeletion and microduplication at 16p11.2 and autism.N Engl J Med358667675
- 24. Phelan MC, Rogers RC, Saul RA, Stapleton GA, Sweet K, et al. (2001) 22q13 deletion syndrome. Am J Med Genet 101: 91–99.MC PhelanRC RogersRA SaulGA StapletonK. Sweet200122q13 deletion syndrome.Am J Med Genet1019199
- 25. Gajecka M, Mackay KL, Shaffer LG (2007) Monosomy 1p36 deletion syndrome. Am J Med Genet C Semin Med Genet 145C: 346–356.M. GajeckaKL MackayLG Shaffer2007Monosomy 1p36 deletion syndrome.Am J Med Genet C Semin Med Genet145C346356
- 26. Abd El-Aziz MM, Barragan I, O'Driscoll CA, Goodstadt L, Prigmore E, et al. (2008) EYS, encoding an ortholog of Drosophila spacemaker, is mutated in autosomal recessive retinitis pigmentosa. Nat Genet 40: 1285–1287.MM Abd El-AzizI. BarraganCA O'DriscollL. GoodstadtE. Prigmore2008EYS, encoding an ortholog of Drosophila spacemaker, is mutated in autosomal recessive retinitis pigmentosa.Nat Genet4012851287
- 27. Collin RW, Littink KW, Klevering BJ, van den Born LI, Koenekoop RK, et al. (2008) Identification of a 2 Mb human ortholog of Drosophila eyes shut/spacemaker that is mutated in patients with retinitis pigmentosa. Am J Hum Genet 83: 594–603.RW CollinKW LittinkBJ KleveringLI van den BornRK Koenekoop2008Identification of a 2 Mb human ortholog of Drosophila eyes shut/spacemaker that is mutated in patients with retinitis pigmentosa.Am J Hum Genet83594603
- 28. Maestrini E, Pagnamenta AT, Lamb JA, Bacchelli E, Sykes NH, et al. (2010) High-density SNP association study and copy number variation analysis of the AUTS1 and AUTS5 loci implicate the IMMP2L-DOCK4 gene region in autism susceptibility. Mol Psychiatry 15: 954–968.E. MaestriniAT PagnamentaJA LambE. BacchelliNH Sykes2010High-density SNP association study and copy number variation analysis of the AUTS1 and AUTS5 loci implicate the IMMP2L-DOCK4 gene region in autism susceptibility.Mol Psychiatry15954968
- 29. Petek E, Windpassinger C, Vincent JB, Cheung J, Boright AP, et al. (2001) Disruption of a novel gene (IMMP2L) by a breakpoint in 7q31 associated with Tourette syndrome. Am J Hum Genet 68: 848–858.E. PetekC. WindpassingerJB VincentJ. CheungAP Boright2001Disruption of a novel gene (IMMP2L) by a breakpoint in 7q31 associated with Tourette syndrome.Am J Hum Genet68848858
- 30. Pagnamenta AT, Bacchelli E, de Jonge MV, Mirza G, Scerri TS, et al. (2010) Characterization of a family with rare deletions in CNTNAP5 and DOCK4 suggests novel risk loci for autism and dyslexia. Biol Psychiatry 68: 320–328.AT PagnamentaE. BacchelliMV de JongeG. MirzaTS Scerri2010Characterization of a family with rare deletions in CNTNAP5 and DOCK4 suggests novel risk loci for autism and dyslexia.Biol Psychiatry68320328
- 31. Hannula-Jouppi K, Kaminen-Ahola N, Taipale M, Eklund R, Nopola-Hemmi J, et al. (2005) The axon guidance receptor gene ROBO1 is a candidate gene for developmental dyslexia. PLoS Genet 1: e50.K. Hannula-JouppiN. Kaminen-AholaM. TaipaleR. EklundJ. Nopola-Hemmi2005The axon guidance receptor gene ROBO1 is a candidate gene for developmental dyslexia.PLoS Genet1e50
- 32. Nopola-Hemmi J, Taipale M, Haltia T, Lehesjoki AE, Voutilainen A, et al. (2000) Two translocations of chromosome 15q associated with dyslexia. J Med Genet 37: 771–775.J. Nopola-HemmiM. TaipaleT. HaltiaAE LehesjokiA. Voutilainen2000Two translocations of chromosome 15q associated with dyslexia.J Med Genet37771775
- 33. Taipale M, Kaminen N, Nopola-Hemmi J, Haltia T, Myllyluoma B, et al. (2003) A candidate gene for developmental dyslexia encodes a nuclear tetratricopeptide repeat domain protein dynamically regulated in brain. Proc Natl Acad Sci U S A 100: 11553–11558.M. TaipaleN. KaminenJ. Nopola-HemmiT. HaltiaB. Myllyluoma2003A candidate gene for developmental dyslexia encodes a nuclear tetratricopeptide repeat domain protein dynamically regulated in brain.Proc Natl Acad Sci U S A1001155311558
- 34. Buonincontri R, Bache I, Silahtaroglu A, Elbro C, Nielsen AM, et al. (2011) A cohort of balanced reciprocal translocations associated with dyslexia: identification of two putative candidate genes at DYX1. Behav Genet 41: 125–133.R. BuonincontriI. BacheA. SilahtarogluC. ElbroAM Nielsen2011A cohort of balanced reciprocal translocations associated with dyslexia: identification of two putative candidate genes at DYX1.Behav Genet41125133
- 35. O'Roak BJ, Deriziotis P, Lee C, Vives L, Schwartz JJ, et al. (2011) Exome sequencing in sporadic autism spectrum disorders identifies severe de novo mutations. Nat Genet 43: 585–589.BJ O'RoakP. DeriziotisC. LeeL. VivesJJ Schwartz2011Exome sequencing in sporadic autism spectrum disorders identifies severe de novo mutations.Nat Genet43585589
- 36. Hamdan FF, Daoud H, Rochefort D, Piton A, Gauthier J, et al. (2010) De novo mutations in FOXP1 in cases with intellectual disability, autism, and language impairment. Am J Hum Genet 87: 671–678.FF HamdanH. DaoudD. RochefortA. PitonJ. Gauthier2010De novo mutations in FOXP1 in cases with intellectual disability, autism, and language impairment.Am J Hum Genet87671678
- 37. Horn D, Kapeller J, Rivera-Brugues N, Moog U, Lorenz-Depiereux B, et al. (2010) Identification of FOXP1 deletions in three unrelated patients with mental retardation and significant speech and language deficits. Hum Mutat 31: E1851–1860.D. HornJ. KapellerN. Rivera-BruguesU. MoogB. Lorenz-Depiereux2010Identification of FOXP1 deletions in three unrelated patients with mental retardation and significant speech and language deficits.Hum Mutat31E18511860
- 38. Teramitsu I, Kudo LC, London SE, Geschwind DH, White SA (2004) Parallel FoxP1 and FoxP2 expression in songbird and human brain predicts functional interaction. J Neurosci 24: 3152–3163.I. TeramitsuLC KudoSE LondonDH GeschwindSA White2004Parallel FoxP1 and FoxP2 expression in songbird and human brain predicts functional interaction.J Neurosci2431523163
- 39. Vernes SC, Newbury DF, Abrahams BS, Winchester L, Nicod J, et al. (2008) A functional genetic link between distinct developmental language disorders. N Engl J Med 359: 2337–2345.SC VernesDF NewburyBS AbrahamsL. WinchesterJ. Nicod2008A functional genetic link between distinct developmental language disorders.N Engl J Med35923372345
- 40. Firth HV, Richards SM, Bevan AP, Clayton S, Corpas M, et al. (2009) DECIPHER: Database of Chromosomal Imbalance and Phenotype in Humans Using Ensembl Resources. Am J Hum Genet 84: 524–533.HV FirthSM RichardsAP BevanS. ClaytonM. Corpas2009DECIPHER: Database of Chromosomal Imbalance and Phenotype in Humans Using Ensembl Resources.Am J Hum Genet84524533
- 41. Bonaglia MC, Ciccone R, Gimelli G, Gimelli S, Marelli S, et al. (2008) Detailed phenotype-genotype study in five patients with chromosome 6q16 deletion: narrowing the critical region for Prader-Willi-like phenotype. Eur J Hum Genet 16: 1443–1449.MC BonagliaR. CicconeG. GimelliS. GimelliS. Marelli2008Detailed phenotype-genotype study in five patients with chromosome 6q16 deletion: narrowing the critical region for Prader-Willi-like phenotype.Eur J Hum Genet1614431449
- 42. Holder JL Jr, Butte NF, Zinn AR (2000) Profound obesity associated with a balanced translocation that disrupts the SIM1 gene. Hum Mol Genet 9: 101–108.JL Holder JrNF ButteAR Zinn2000Profound obesity associated with a balanced translocation that disrupts the SIM1 gene.Hum Mol Genet9101108
- 43. Jamain S, Betancur C, Quach H, Philippe A, Fellous M, et al. (2002) Linkage and association of the glutamate receptor 6 gene with autism. Mol Psychiatry 7: 302–310.S. JamainC. BetancurH. QuachA. PhilippeM. Fellous2002Linkage and association of the glutamate receptor 6 gene with autism.Mol Psychiatry7302310
- 44. Strungaru MH, Dinu I, Walter MA (2007) Genotype-phenotype correlations in Axenfeld-Rieger malformation and glaucoma patients with FOXC1 and PITX2 mutations. Invest Ophthalmol Vis Sci 48: 228–237.MH StrungaruI. DinuMA Walter2007Genotype-phenotype correlations in Axenfeld-Rieger malformation and glaucoma patients with FOXC1 and PITX2 mutations.Invest Ophthalmol Vis Sci48228237
- 45. Ballif BC, Theisen A, Coppinger J, Gowans GC, Hersh JH, et al. (2008) Expanding the clinical phenotype of the 3q29 microdeletion syndrome and characterization of the reciprocal microduplication. Mol Cytogenet 1: 8.BC BallifA. TheisenJ. CoppingerGC GowansJH Hersh2008Expanding the clinical phenotype of the 3q29 microdeletion syndrome and characterization of the reciprocal microduplication.Mol Cytogenet18
- 46. Lisi EC, Hamosh A, Doheny KF, Squibb E, Jackson B, et al. (2008) 3q29 interstitial microduplication: a new syndrome in a three-generation family. Am J Med Genet A 146A: 601–609.EC LisiA. HamoshKF DohenyE. SquibbB. Jackson20083q29 interstitial microduplication: a new syndrome in a three-generation family.Am J Med Genet A146A601609
- 47. Cardoso C, Leventer RJ, Ward HL, Toyo-Oka K, Chung J, et al. (2003) Refinement of a 400-kb critical region allows genotypic differentiation between isolated lissencephaly, Miller-Dieker syndrome, and other phenotypes secondary to deletions of 17p13.3. Am J Hum Genet 72: 918–930.C. CardosoRJ LeventerHL WardK. Toyo-OkaJ. Chung2003Refinement of a 400-kb critical region allows genotypic differentiation between isolated lissencephaly, Miller-Dieker syndrome, and other phenotypes secondary to deletions of 17p13.3.Am J Hum Genet72918930
- 48. Nagamani SC, Zhang F, Shchelochkov OA, Bi W, Ou Z, et al. (2009) Microdeletions including YWHAE in the Miller-Dieker syndrome region on chromosome 17p13.3 result in facial dysmorphisms, growth restriction, and cognitive impairment. J Med Genet 46: 825–833.SC NagamaniF. ZhangOA ShchelochkovW. BiZ. Ou2009Microdeletions including YWHAE in the Miller-Dieker syndrome region on chromosome 17p13.3 result in facial dysmorphisms, growth restriction, and cognitive impairment.J Med Genet46825833
- 49. Itsara A, Cooper GM, Baker C, Girirajan S, Li J, et al. (2009) Population analysis of large copy number variants and hotspots of human genetic disease. Am J Hum Genet 84: 148–161.A. ItsaraGM CooperC. BakerS. GirirajanJ. Li2009Population analysis of large copy number variants and hotspots of human genetic disease.Am J Hum Genet84148161
- 50. Raskind W (2001) Current understanding of the genetic basis of reading and spelling disability. Learning Disability Quarterly 24: 141–157.W. Raskind2001Current understanding of the genetic basis of reading and spelling disability.Learning Disability Quarterly24141157
- 51. Bakkaloglu B, O'Roak BJ, Louvi A, Gupta AR, Abelson JF, et al. (2008) Molecular cytogenetic analysis and resequencing of contactin associated protein-like 2 in autism spectrum disorders. Am J Hum Genet 82: 165–173.B. BakkalogluBJ O'RoakA. LouviAR GuptaJF Abelson2008Molecular cytogenetic analysis and resequencing of contactin associated protein-like 2 in autism spectrum disorders.Am J Hum Genet82165173
- 52. de la Barra F, Skoknic V, Alliende A, Raimann E, Cortes F, et al. (1986) [Twins with autism and mental retardation associated with balanced (7;20) chromosomal translocation]. Rev Chil Pediatr 57: 549–554.F. de la BarraV. SkoknicA. AlliendeE. RaimannF. Cortes1986[Twins with autism and mental retardation associated with balanced (7;20) chromosomal translocation].Rev Chil Pediatr57549554
- 53. Huang XL, Zou YS, Maher TA, Newton S, Milunsky JM (2010) A de novo balanced translocation breakpoint truncating the autism susceptibility candidate 2 (AUTS2) gene in a patient with autism. Am J Med Genet A 152A: 2112–2114.XL HuangYS ZouTA MaherS. NewtonJM Milunsky2010A de novo balanced translocation breakpoint truncating the autism susceptibility candidate 2 (AUTS2) gene in a patient with autism.Am J Med Genet A152A21122114
- 54. Kalscheuer VM, FitzPatrick D, Tommerup N, Bugge M, Niebuhr E, et al. (2007) Mutations in autism susceptibility candidate 2 (AUTS2) in patients with mental retardation. Hum Genet 121: 501–509.VM KalscheuerD. FitzPatrickN. TommerupM. BuggeE. Niebuhr2007Mutations in autism susceptibility candidate 2 (AUTS2) in patients with mental retardation.Hum Genet121501509
- 55. Sultana R, Yu CE, Yu J, Munson J, Chen D, et al. (2002) Identification of a novel gene on chromosome 7q11.2 interrupted by a translocation breakpoint in a pair of autistic twins. Genomics 80: 129–134.R. SultanaCE YuJ. YuJ. MunsonD. Chen2002Identification of a novel gene on chromosome 7q11.2 interrupted by a translocation breakpoint in a pair of autistic twins.Genomics80129134
- 56. Fletcher J, Shaywitz S, Shaywitz B (1999) Comorbidity of learning and attention disorders. Separate but equal. Pediatr Clin North Am 46: 885–897.J. FletcherS. ShaywitzB. Shaywitz1999Comorbidity of learning and attention disorders. Separate but equal.Pediatr Clin North Am46885897
- 57. Willcutt EG, Pennington BF, Olson RK, DeFries JC (2007) Understanding comorbidity: a twin study of reading disability and attention-deficit/hyperactivity disorder. Am J Med Genet B Neuropsychiatr Genet 144B: 709–714.EG WillcuttBF PenningtonRK OlsonJC DeFries2007Understanding comorbidity: a twin study of reading disability and attention-deficit/hyperactivity disorder.Am J Med Genet B Neuropsychiatr Genet144B709714
- 58. Willcutt EG, Betjemann RS, McGrath LM, Chhabildas NA, Olson RK, et al. (2010) Etiology and neuropsychology of comorbidity between RD and ADHD: the case for multiple-deficit models. Cortex 46: 1345–1361.EG WillcuttRS BetjemannLM McGrathNA ChhabildasRK Olson2010Etiology and neuropsychology of comorbidity between RD and ADHD: the case for multiple-deficit models.Cortex4613451361
- 59. Bedogni F, Hodge RD, Nelson BR, Frederick EA, Shiba N, et al. (2010) Autism susceptibility candidate 2 (Auts2) encodes a nuclear protein expressed in developing brain regions implicated in autism neuropathology. Gene Expr Patterns 10: 9–15.F. BedogniRD HodgeBR NelsonEA FrederickN. Shiba2010Autism susceptibility candidate 2 (Auts2) encodes a nuclear protein expressed in developing brain regions implicated in autism neuropathology.Gene Expr Patterns10915
- 60. Green RE, Krause J, Briggs AW, Maricic T, Stenzel U, et al. (2010) A draft sequence of the Neandertal genome. Science 328: 710–722.RE GreenJ. KrauseAW BriggsT. MaricicU. Stenzel2010A draft sequence of the Neandertal genome.Science328710722
- 61. Levy D, Ronemus M, Yamrom B, Lee YH, Leotta A, et al. (2011) Rare de novo and transmitted copy-number variation in autistic spectrum disorders. Neuron 70: 886–897.D. LevyM. RonemusB. YamromYH LeeA. Leotta2011Rare de novo and transmitted copy-number variation in autistic spectrum disorders.Neuron70886897
- 62. Pinto D, Pagnamenta AT, Klei L, Anney R, Merico D, et al. (2010) Functional impact of global rare copy number variation in autism spectrum disorders. Nature 466: 368–372.D. PintoAT PagnamentaL. KleiR. AnneyD. Merico2010Functional impact of global rare copy number variation in autism spectrum disorders.Nature466368372
- 63. Sanders SJ, Ercan-Sencicek AG, Hus V, Luo R, Murtha MT, et al. (2011) Multiple Recurrent De Novo CNVs, Including Duplications of the 7q11.23 Williams Syndrome Region, Are Strongly Associated with Autism. Neuron 70: 863–885.SJ SandersAG Ercan-SencicekV. HusR. LuoMT Murtha2011Multiple Recurrent De Novo CNVs, Including Duplications of the 7q11.23 Williams Syndrome Region, Are Strongly Associated with Autism.Neuron70863885
- 64. Girirajan S, Rosenfeld JA, Cooper GM, Antonacci F, Siswara P, et al. (2010) A recurrent 16p12.1 microdeletion supports a two-hit model for severe developmental delay. Nat Genet 42: 203–209.S. GirirajanJA RosenfeldGM CooperF. AntonacciP. Siswara2010A recurrent 16p12.1 microdeletion supports a two-hit model for severe developmental delay.Nat Genet42203209
- 65. Kim HG, Kishikawa S, Higgins AW, Seong IS, Donovan DJ, et al. (2008) Disruption of neurexin 1 associated with autism spectrum disorder. Am J Hum Genet 82: 199–207.HG KimS. KishikawaAW HigginsIS SeongDJ Donovan2008Disruption of neurexin 1 associated with autism spectrum disorder.Am J Hum Genet82199207
- 66. Szatmari P, Paterson AD, Zwaigenbaum L, Roberts W, Brian J, et al. (2007) Mapping autism risk loci using genetic linkage and chromosomal rearrangements. Nat Genet 39: 319–328.P. SzatmariAD PatersonL. ZwaigenbaumW. RobertsJ. Brian2007Mapping autism risk loci using genetic linkage and chromosomal rearrangements.Nat Genet39319328
- 67. Glessner JT, Wang K, Cai G, Korvatska O, Kim CE, et al. (2009) Autism genome-wide copy number variation reveals ubiquitin and neuronal genes. Nature 459: 569–573.JT GlessnerK. WangG. CaiO. KorvatskaCE Kim2009Autism genome-wide copy number variation reveals ubiquitin and neuronal genes.Nature459569573
- 68. Bucan M, Abrahams BS, Wang K, Glessner JT, Herman EI, et al. (2009) Genome-wide analyses of exonic copy number variants in a family-based study point to novel autism susceptibility genes. PLoS Genet 5: e1000536.M. BucanBS AbrahamsK. WangJT GlessnerEI Herman2009Genome-wide analyses of exonic copy number variants in a family-based study point to novel autism susceptibility genes.PLoS Genet5e1000536
- 69. Noor A, Whibley A, Marshall CR, Gianakopoulos PJ, Piton A, et al. (2010) Disruption at the PTCHD1 Locus on Xp22.11 in Autism spectrum disorder and intellectual disability. Sci Transl Med 2: 49ra68.A. NoorA. WhibleyCR MarshallPJ GianakopoulosA. Piton2010Disruption at the PTCHD1 Locus on Xp22.11 in Autism spectrum disorder and intellectual disability.Sci Transl Med249ra68
- 70. Brkanac Z, Chapman NH, Igo RP Jr, Matsushita MM, Nielsen K, et al. (2008) Genome Scan of a Nonword Repetition Phenotype in Families with Dyslexia: Evidence for Multiple Loci. Behav Genet. Z. BrkanacNH ChapmanRP Igo JrMM MatsushitaK. Nielsen2008Genome Scan of a Nonword Repetition Phenotype in Families with Dyslexia: Evidence for Multiple Loci.Behav Genet
- 71. Rubenstein K, Matsushita M, Berninger VW, Raskind WH, Wijsman EM (2011) Genome scan for spelling deficits: effects of verbal IQ on models of transmission and trait gene localization. Behav Genet 41: 31–42.K. RubensteinM. MatsushitaVW BerningerWH RaskindEM Wijsman2011Genome scan for spelling deficits: effects of verbal IQ on models of transmission and trait gene localization.Behav Genet413142
- 72. Berninger V (2006) A developmental approach to learning disabilities. In: Siegel I, Renninger A, editors. Handbook of Child Psychology, Vol IV, Child Psychology and Practice. New York: John Wiley and Sons. pp. 420–452.V. Berninger2006A developmental approach to learning disabilities.I. SiegelA. RenningerHandbook of Child Psychology, Vol IV, Child Psychology and PracticeNew YorkJohn Wiley and Sons420452
- 73. Fischbach GD, Lord C (2010) The Simons Simplex Collection: a resource for identification of autism genetic risk factors. Neuron 68: 192–195.GD FischbachC. Lord2010The Simons Simplex Collection: a resource for identification of autism genetic risk factors.Neuron68192195
- 74. Lord C, Risi S, Lambrecht L, Cook EH Jr, Leventhal BL, et al. (2000) The autism diagnostic observation schedule-generic: a standard measure of social and communication deficits associated with the spectrum of autism. J Autism Dev Disord 30: 205–223.C. LordS. RisiL. LambrechtEH Cook JrBL Leventhal2000The autism diagnostic observation schedule-generic: a standard measure of social and communication deficits associated with the spectrum of autism.J Autism Dev Disord30205223
- 75. Lord C, Rutter M, Le Couteur A (1994) Autism Diagnostic Interview-Revised: a revised version of a diagnostic interview for caregivers of individuals with possible pervasive developmental disorders. J Autism Dev Disord 24: 659–685.C. LordM. RutterA. Le Couteur1994Autism Diagnostic Interview-Revised: a revised version of a diagnostic interview for caregivers of individuals with possible pervasive developmental disorders.J Autism Dev Disord24659685
- 76. Moldin SO (2003) NIMH Human Genetics Initiative: 2003 update. Am J Psychiatry 160: 621–622.SO Moldin2003NIMH Human Genetics Initiative: 2003 update.Am J Psychiatry160621622
- 77. Kessler RC, Ustun TB (2004) The World Mental Health (WMH) Survey Initiative Version of the World Health Organization (WHO) Composite International Diagnostic Interview (CIDI). Int J Methods Psychiatr Res 13: 93–121.RC KesslerTB Ustun2004The World Mental Health (WMH) Survey Initiative Version of the World Health Organization (WHO) Composite International Diagnostic Interview (CIDI).Int J Methods Psychiatr Res1393121
- 78. Talati A, Fyer AJ, Weissman MM (2008) A comparison between screened NIMH and clinically interviewed control samples on neuroticism and extraversion. Mol Psychiatry 13: 122–130.A. TalatiAJ FyerMM Weissman2008A comparison between screened NIMH and clinically interviewed control samples on neuroticism and extraversion.Mol Psychiatry13122130
- 79. Baum AE, Akula N, Cabanero M, Cardona I, Corona W, et al. (2008) A genome-wide association study implicates diacylglycerol kinase eta (DGKH) and several other genes in the etiology of bipolar disorder. Mol Psychiatry 13: 197–207.AE BaumN. AkulaM. CabaneroI. CardonaW. Corona2008A genome-wide association study implicates diacylglycerol kinase eta (DGKH) and several other genes in the etiology of bipolar disorder.Mol Psychiatry13197207
- 80. Li B, Leal SM (2008) Methods for detecting associations with rare variants for common diseases: application to analysis of sequence data. Am J Hum Genet 83: 311–321.B. LiSM Leal2008Methods for detecting associations with rare variants for common diseases: application to analysis of sequence data.Am J Hum Genet83311321
- 81. Simon-Sanchez J, Scholz S, Fung HC, Matarin M, Hernandez D, et al. (2007) Genome-wide SNP assay reveals structural genomic variation, extended homozygosity and cell-line induced alterations in normal individuals. Hum Mol Genet 16: 1–14.J. Simon-SanchezS. ScholzHC FungM. MatarinD. Hernandez2007Genome-wide SNP assay reveals structural genomic variation, extended homozygosity and cell-line induced alterations in normal individuals.Hum Mol Genet16114
- 82. Albert MA, Danielson E, Rifai N, Ridker PM (2001) Effect of statin therapy on C-reactive protein levels: the pravastatin inflammation/CRP evaluation (PRINCE): a randomized trial and cohort study. JAMA 286: 64–70.MA AlbertE. DanielsonN. RifaiPM Ridker2001Effect of statin therapy on C-reactive protein levels: the pravastatin inflammation/CRP evaluation (PRINCE): a randomized trial and cohort study.JAMA2866470
- 83. Simon JA, Lin F, Hulley SB, Blanche PJ, Waters D, et al. (2006) Phenotypic predictors of response to simvastatin therapy among African-Americans and Caucasians: the Cholesterol and Pharmacogenetics (CAP) Study. Am J Cardiol 97: 843–850.JA SimonF. LinSB HulleyPJ BlancheD. Waters2006Phenotypic predictors of response to simvastatin therapy among African-Americans and Caucasians: the Cholesterol and Pharmacogenetics (CAP) Study.Am J Cardiol97843850
- 84. Melzer D, Perry JR, Hernandez D, Corsi AM, Stevens K, et al. (2008) A genome-wide association study identifies protein quantitative trait loci (pQTLs). PLoS Genet 4: e1000072.D. MelzerJR PerryD. HernandezAM CorsiK. Stevens2008A genome-wide association study identifies protein quantitative trait loci (pQTLs).PLoS Genet4e1000072
- 85. Cooper GM, Coe BP, Girirajan S, Rosenfeld JA, Vu TH, et al. (2011) A copy number variation morbidity map of developmental delay. Nat Genet 43: 838–846.GM CooperBP CoeS. GirirajanJA RosenfeldTH Vu2011A copy number variation morbidity map of developmental delay.Nat Genet43838846
- 86. Selzer RR, Richmond TA, Pofahl NJ, Green RD, Eis PS, et al. (2005) Analysis of chromosome breakpoints in neuroblastoma at sub-kilobase resolution using fine-tiling oligonucleotide array CGH. Genes Chromosomes Cancer 44: 305–319.RR SelzerTA RichmondNJ PofahlRD GreenPS Eis2005Analysis of chromosome breakpoints in neuroblastoma at sub-kilobase resolution using fine-tiling oligonucleotide array CGH.Genes Chromosomes Cancer44305319
- 87. Day N, Hemmaplardh A, Thurman RE, Stamatoyannopoulos JA, Noble WS (2007) Unsupervised segmentation of continuous genomic data. Bioinformatics 23: 1424–1426.N. DayA. HemmaplardhRE ThurmanJA StamatoyannopoulosWS Noble2007Unsupervised segmentation of continuous genomic data.Bioinformatics2314241426