DNA was extracted from blood using DNeasy Blood and Tissue Kit (Qiagen Inc., Valencia, CA). Genome-wide SNP screening on the two affected individuals from Family 1 was performed using Affymetrix 6.0 SNP Array (Affymetrix, USA). SNP genotyping data was used to perform homozygosity and copy number analyses. Homozygosity mapping was performed as previously described35 (link). CNV analysis was performed using four algorithms in parallel (Birdsuite, PennCNV, Nexus and Affymetrix Genotyping Console), requiring a minimum of 10 probes and the concordance of at least two algorithms for specificity. To distinguish rare from common variation, results were matched against a catalog of reference CNVs generated via the same pipeline, using 1251 samples from the HapMap project (11 subpopulations, public release #3). For additional stringency, CNVs were further matched with CNVs present in the Database of Genomic Variants (DGV) as well as an internal database of CNVs from over 12000 unrelated research samples (including >1000 from the Middle East with a variety of neurodevelopmental conditions, and >105 healthy Saudi individuals).