Genomic architecture of autism from comprehensive whole-genome sequence annotation.

Fully understanding autism spectrum disorder (ASD) genetics requires whole-genome sequencing (WGS). We present the latest release of the Autism Speaks MSSNG resource, which includes WGS data from 5,100 individuals with ASD and 6,212 non-ASD parents and siblings (total n = 11,312). Examining a wide variety of genetic variants in MSSNG and the Simons Simplex Collection (SSC; n = 9,205), we identified ASD-associated rare variants in 718/5,100 individuals with ASD from MSSNG (14.1%) and 350/2,419 from SSC (14.5%). Considering genomic architecture, 52% were nuclear sequence-level variants, 46% were nuclear structural variants (including copy-number variants, inversions, large insertions, uniparental isodisomies, and tandem repeat expansions), and 2% were mitochondrial variants. Our study provides a guidebook for exploring genotype-phenotype correlations in families who carry ASD-associated rare variants and serves as an entry point to the expanded studies required to dissect the etiology in the ∼85% of the ASD population that remain idiopathic.

Gene copy number variation and pediatric mental health/neurodevelopment in a general population.

We assessed the relationship of gene copy number variation (CNV) in mental health/neurodevelopmental traits and diagnoses, physical health and cognition in a community sample of 7100 unrelated children and youth of European or East Asian ancestry (Spit for Science). Clinically significant or susceptibility CNVs were present in 3.9% of participants and were associated with elevated scores on a continuous measure of attention-deficit/hyperactivity disorder (ADHD) traits (P = 5.0 × 10-3), longer response inhibition (a cognitive deficit found in several mental health and neurodevelopmental disorders; P = 1.0 × 10-2) and increased prevalence of mental health diagnoses (P = 1.9 × 10-6, odds ratio: 3.09), specifically ADHD, autism spectrum disorder anxiety and learning problems/learning disorder (P's < 0.01). There was an increased burden of rare deletions in gene-sets related to brain function or expression in brain associated with more ADHD traits. With the current mental health crisis, our data established a baseline for delineating genetic contributors in pediatric-onset conditions.

OTUD7A Regulates Neurodevelopmental Phenotypes in the 15q13.3 Microdeletion Syndrome.

Copy-number variations (CNVs) are strong risk factors for neurodevelopmental and psychiatric disorders. The 15q13.3 microdeletion syndrome region contains up to ten genes and is associated with numerous conditions, including autism spectrum disorder (ASD), epilepsy, schizophrenia, and intellectual disability; however, the mechanisms underlying the pathogenesis of 15q13.3 microdeletion syndrome remain unknown. We combined whole-genome sequencing, human brain gene expression (proteome and transcriptome), and a mouse model with a syntenic heterozygous deletion (Df(h15q13)/+ mice) and determined that the microdeletion results in abnormal development of cortical dendritic spines and dendrite outgrowth. Analysis of large-scale genomic, transcriptomic, and proteomic data identified OTUD7A as a critical gene for brain function. OTUD7A was found to localize to dendritic and spine compartments in cortical neurons, and its reduced levels in Df(h15q13)/+ cortical neurons contributed to the dendritic spine and dendrite outgrowth deficits. Our results reveal OTUD7A as a major regulatory gene for 15q13.3 microdeletion syndrome phenotypes that contribute to the disease mechanism through abnormal cortical neuron morphological development.

Genome sequencing for detection of pathogenic deep intronic variation: A clinical case report illustrating opportunities and challenges.

Variants in JAM3 have been reported in four families manifesting a severe autosomal recessive disorder characterized by hemorrhagic destruction of the brain, subependymal calcification, and cataracts. We describe a 7-year-old male with a similar presentation found by research-based quad genome sequencing to have two novel splicing variants in trans in JAM3, including one deep intronic variant (NM_032801.4: c.256+1260G>C) not detectable by standard exome sequencing. Targeted sequencing of RNA isolated from transformed lymphoblastoid cell lines confirmed that each of the two variants has a deleterious effect on JAM3 mRNA splicing. The role for genome sequencing as a clinical diagnostic test extends to those patients with phenotypes strongly suggestive of a specific Mendelian disorder, especially when the causal genetic variant(s) are not found by a more targeted approach. Barriers to diagnosis via identification of pathogenic deep intronic variation include lack of laboratory consensus regarding in silico splicing prediction tools and limited access to clinically validated confirmatory RNA experiments.

A novel intronic variant in UBE3A identified by genome sequencing in a patient with an atypical presentation of Angelman syndrome.

Angelman syndrome (AS) is a genetic neurodevelopmental disorder caused by loss or deficient expression of UBE3A on the maternally inherited allele. In 10-15% of individuals with a clinical diagnosis of AS, a molecular diagnosis cannot be established with conventional testing. We describe a 13-year-old male with an atypical presentation of AS, who was found to have a novel, maternally inherited, intronic variant in UBE3A (c.3-12T>A) using genome sequencing (GS). Targeted sequencing of RNA isolated from blood confirmed the creation of a new acceptor splice site. These GS results ended a six-year diagnostic odyssey and revealed a 50% recurrence risk for the unaffected parents. This case illustrates a previously unreported splicing variant causing AS. Intronic variants identifiable by GS may account for a proportion of individuals who are suspected of having well-known genetic disorders despite negative prior genetic testing.

Copy number variation in fetal alcohol spectrum disorder.

Fetal alcohol spectrum disorder (FASD) is characterized by a combination of neurological, developmental, and congenital defects that may occur as a consequence of prenatal alcohol exposure. Earlier reports showed that large chromosomal anomalies may link to FASD. Here, we examined the prevalence and types of copy number variations (CNVs) in FASD cases previously diagnosed by a multidisciplinary FASD team in sites across Canada. We genotyped 95 children with FASD and 87 age-matched, typically developing controls on the Illumina Human Omni2.5 SNP (single nucleotide polymorphisms) array platform. We compared their CNVs with those of 10 851 population controls to identify rare CNVs (<0.1% frequency), which may include large unbalanced chromosomal abnormalities, that might be relevant to FASD. In 12/95 (13%) of the FASD cases, rare CNVs were found that impact potentially clinically relevant developmental genes, including the CACNA1H involved in epilepsy and autism, the 3q29 deletion disorder, and others. Our results show that a subset of children diagnosed with FASD have chromosomal deletions and duplications that may co-occur or explain the neurodevelopmental impairments in a diagnosed cohort of FASD individuals. Children suspected to have FASD with or without sentinel facial features of fetal alcohol syndrome and neurodevelopmental delays should potentially be evaluated by a clinical geneticist and possibly have genetic investigations as appropriate to exclude other etiologies.

Disruption of the ASTN2/TRIM32 locus at 9q33.1 is a risk factor in males for autism spectrum disorders, ADHD and other neurodevelopmental phenotypes.

Rare copy number variants (CNVs) disrupting ASTN2 or both ASTN2 and TRIM32 have been reported at 9q33.1 by genome-wide studies in a few individuals with neurodevelopmental disorders (NDDs). The vertebrate-specific astrotactins, ASTN2 and its paralog ASTN1, have key roles in glial-guided neuronal migration during brain development. To determine the prevalence of astrotactin mutations and delineate their associated phenotypic spectrum, we screened ASTN2/TRIM32 and ASTN1 (1q25.2) for exonic CNVs in clinical microarray data from 89 985 individuals across 10 sites, including 64 114 NDD subjects. In this clinical dataset, we identified 46 deletions and 12 duplications affecting ASTN2. Deletions of ASTN1 were much rarer. Deletions near the 3' terminus of ASTN2, which would disrupt all transcript isoforms (a subset of these deletions also included TRIM32), were significantly enriched in the NDD subjects (P = 0.002) compared with 44 085 population-based controls. Frequent phenotypes observed in individuals with such deletions include autism spectrum disorder (ASD), attention deficit hyperactivity disorder (ADHD), speech delay, anxiety and obsessive compulsive disorder (OCD). The 3'-terminal ASTN2 deletions were significantly enriched compared with controls in males with NDDs, but not in females. Upon quantifying ASTN2 human brain RNA, we observed shorter isoforms expressed from an alternative transcription start site of recent evolutionary origin near the 3' end. Spatiotemporal expression profiling in the human brain revealed consistently high ASTN1 expression while ASTN2 expression peaked in the early embryonic neocortex and postnatal cerebellar cortex. Our findings shed new light on the role of the astrotactins in psychopathology and their interplay in human neurodevelopment.

A high-resolution copy-number variation resource for clinical and population genetics.

PURPOSE: Chromosomal microarray analysis to assess copy-number variation has become a first-tier genetic diagnostic test for individuals with unexplained neurodevelopmental disorders or multiple congenital anomalies. More than 100 cytogenetic laboratories worldwide use the new ultra-high resolution Affymetrix CytoScan-HD array to genotype hundreds of thousands of samples per year. Our aim was to develop a copy-number variation resource from a new population sample that would enable more accurate interpretation of clinical genetics data on this microarray platform and others. METHODS: Genotyping of 1,000 adult volunteers who are broadly representative of the Ontario population (as obtained from the Ontario Population Genomics Platform) was performed with the CytoScan-HD microarray system, which has 2.7 million probes. Four independent algorithms were applied to detect copy-number variations. Reproducibility and validation metrics were quantified using sample replicates and quantitative-polymerase chain reaction, respectively. RESULTS: DNA from 873 individuals passed quality control and we identified 71,178 copy-number variations (81 copy-number variations/individual); 9.8% (6,984) of these copy-number variations were previously unreported. After applying three layers of filtering criteria, from our highest confidence copy-number variation data set we obtained >95% reproducibility and >90% validation rates (73% of these copy-number variations overlapped at least one gene). CONCLUSION: The genotype data and annotated copy-number variations for this largely Caucasian population will represent a valuable public resource enabling clinical genetics research and diagnostics.

Comprehensive whole-genome sequence analyses provide insights into the genomic architecture of cerebral palsy.

We performed whole-genome sequencing (WGS) in 327 children with cerebral palsy (CP) and their biological parents. We classified 37 of 327 (11.3%) children as having pathogenic/likely pathogenic (P/LP) variants and 58 of 327 (17.7%) as having variants of uncertain significance. Multiple classes of P/LP variants included single-nucleotide variants (SNVs)/indels (6.7%), copy number variations (3.4%) and mitochondrial mutations (1.5%). The COL4A1 gene had the most P/LP SNVs. We also analyzed two pediatric control cohorts (n = 203 trios and n = 89 sib-pair families) to provide a baseline for de novo mutation rates and genetic burden analyses, the latter of which demonstrated associations between de novo deleterious variants and genes related to the nervous system. An enrichment analysis revealed previously undescribed plausible candidate CP genes (SMOC1, KDM5B, BCL11A and CYP51A1). A multifactorial CP risk profile and substantial presence of P/LP variants combine to support WGS in the diagnostic work-up across all CP and related phenotypes.

Disruption of DDX53 coding sequence has limited impact on iPSC-derived human NGN2 neurons.

BACKGROUND: The X-linked PTCHD1 locus is strongly associated with autism spectrum disorder (ASD). Males who carry chromosome microdeletions of PTCHD1 antisense long non-coding RNA (PTCHD1-AS)/DEAD-box helicase 53 (DDX53) have ASD, or a sub-clinical form called Broader Autism Phenotype. If the deletion extends beyond PTCHD1-AS/DDX53 to the next gene, PTCHD1, which is protein-coding, the individuals typically have ASD and intellectual disability (ID). Three male siblings with a 90 kb deletion that affects only PTCHD1-AS (and not including DDX53) have ASD. We performed a functional analysis of DDX53 to examine its role in NGN2 neurons. METHODS: We used the clustered regularly interspaced short palindromic repeats (CRISPR) gene editing strategy to knock out DDX53 protein by inserting 3 termination codons (3TCs) into two different induced pluripotent stem cell (iPSC) lines. DDX53 CRISPR-edited iPSCs were differentiated into cortical excitatory neurons by Neurogenin 2 (NGN-2) directed differentiation. The functional differences of DDX53-3TC neurons compared to isogenic control neurons with molecular and electrophysiological approaches were assessed. RESULTS: Isogenic iPSC-derived control neurons exhibited low levels of DDX53 transcripts. Transcriptional analysis revealed the generation of excitatory cortical neurons and DDX53 protein was not detected in iPSC-derived control neurons by western blot. Control lines and DDX53-3TC neurons were active in the multi-electrode array, but no overt electrophysiological phenotype in either isogenic line was observed. CONCLUSION: DDX53-3TC mutation does not alter NGN2 neuronal function in these experiments, suggesting that synaptic deficits causing ASD are unlikely in this cell type.

Mutational Landscape of Autism Spectrum Disorder Brain Tissue.

Rare post-zygotic mutations in the brain are now known to contribute to several neurodevelopmental disorders, including autism spectrum disorder (ASD). However, due to the limited availability of brain tissue, most studies rely on estimates of mosaicism from peripheral samples. In this study, we undertook whole exome sequencing on brain tissue from 26 ASD brain donors from the Harvard Brain Tissue Resource Center (HBTRC) and ascertained the presence of post-zygotic and germline mutations categorized as pathological, including those impacting known ASD-implicated genes. Although quantification did not reveal enrichment for post-zygotic mutations compared with the controls (n = 15), a small number of pathogenic, potentially ASD-implicated mutations were identified, notably in TRAK1 and CLSTN3. Furthermore, germline mutations were identified in the same tissue samples in several key ASD genes, including PTEN, SC1A, CDH13, and CACNA1C. The establishment of tissue resources that are available to the scientific community will facilitate the discovery of new mutations for ASD and other neurodevelopmental disorders.

Genome-wide rare variant score associates with morphological subtypes of autism spectrum disorder.

Defining different genetic subtypes of autism spectrum disorder (ASD) can enable the prediction of developmental outcomes. Based on minor physical and major congenital anomalies, we categorize 325 Canadian children with ASD into dysmorphic and nondysmorphic subgroups. We develop a method for calculating a patient-level, genome-wide rare variant score (GRVS) from whole-genome sequencing (WGS) data. GRVS is a sum of the number of variants in morphology-associated coding and non-coding regions, weighted by their effect sizes. Probands with dysmorphic ASD have a significantly higher GRVS compared to those with nondysmorphic ASD (P = 0.03). Using the polygenic transmission disequilibrium test, we observe an over-transmission of ASD-associated common variants in nondysmorphic ASD probands (P = 2.9 × 10-3). These findings replicate using WGS data from 442 ASD probands with accompanying morphology data from the Simons Simplex Collection. Our results provide support for an alternative genomic classification of ASD subgroups using morphology data, which may inform intervention protocols.

A large data resource of genomic copy number variation across neurodevelopmental disorders.

Copy number variations (CNVs) are implicated across many neurodevelopmental disorders (NDDs) and contribute to their shared genetic etiology. Multiple studies have attempted to identify shared etiology among NDDs, but this is the first genome-wide CNV analysis across autism spectrum disorder (ASD), attention deficit hyperactivity disorder (ADHD), schizophrenia (SCZ), and obsessive-compulsive disorder (OCD) at once. Using microarray (Affymetrix CytoScan HD), we genotyped 2,691 subjects diagnosed with an NDD (204 SCZ, 1,838 ASD, 427 ADHD and 222 OCD) and 1,769 family members, mainly parents. We identified rare CNVs, defined as those found in <0.1% of 10,851 population control samples. We found clinically relevant CNVs (broadly defined) in 284 (10.5%) of total subjects, including 22 (10.8%) among subjects with SCZ, 209 (11.4%) with ASD, 40 (9.4%) with ADHD, and 13 (5.6%) with OCD. Among all NDD subjects, we identified 17 (0.63%) with aneuploidies and 115 (4.3%) with known genomic disorder variants. We searched further for genes impacted by different CNVs in multiple disorders. Examples of NDD-associated genes linked across more than one disorder (listed in order of occurrence frequency) are NRXN1, SEH1L, LDLRAD4, GNAL, GNG13, MKRN1, DCTN2, KNDC1, PCMTD2, KIF5A, SYNM, and long non-coding RNAs: AK127244 and PTCHD1-AS. We demonstrated that CNVs impacting the same genes could potentially contribute to the etiology of multiple NDDs. The CNVs identified will serve as a useful resource for both research and diagnostic laboratories for prioritization of variants.

Germline and somatic mutations in STXBP1 with diverse neurodevelopmental phenotypes.

OBJECTIVE: To expand the clinical phenotype associated with STXBP1 gene mutations and to understand the effect of STXBP1 mutations in the pathogenesis of focal cortical dysplasia (FCD). METHODS: Patients with STXBP1 mutations were identified in various ways: as part of a retrospective cohort study of epileptic encephalopathy; through clinical referrals of individuals (10,619) with developmental delay (DD) for chromosomal microarray; and from a collection of 5,205 individuals with autism spectrum disorder (ASD) examined by whole-genome sequencing. RESULTS: Seven patients with heterozygous de novo mutations affecting the coding region of STXBP1 were newly identified. Three cases had radiologic evidence suggestive of FCD. One male patient with early infantile epileptic encephalopathy, DD, and ASD achieved complete seizure remission following resection of dysplastic brain tissue. Examination of excised brain tissue identified mosaicism for STXBP1, providing evidence for a somatic mechanism. Cell-type expression analysis suggested neuron-specific expression. A comprehensive analysis of the published data revealed that 3.1% of severe epilepsy cases carry a pathogenic de novo mutation within STXBP1. By contrast, ASD was rarely associated with mutations in this gene in our large cohorts. CONCLUSIONS: STXBP1 mutations are an important cause of epilepsy and are also rarely associated with ASD. In a case with histologically proven FCD, an STXBP1 somatic mutation was identified, suggesting a role in its etiology. Removing such tissue may be curative for STXBP1-related epilepsy.

Whole genome sequencing resource identifies 18 new candidate genes for autism spectrum disorder.

We are performing whole-genome sequencing of families with autism spectrum disorder (ASD) to build a resource (MSSNG) for subcategorizing the phenotypes and underlying genetic factors involved. Here we report sequencing of 5,205 samples from families with ASD, accompanied by clinical information, creating a database accessible on a cloud platform and through a controlled-access internet portal. We found an average of 73.8 de novo single nucleotide variants and 12.6 de novo insertions and deletions or copy number variations per ASD subject. We identified 18 new candidate ASD-risk genes and found that participants bearing mutations in susceptibility genes had significantly lower adaptive ability (P = 6 × 10-4). In 294 of 2,620 (11.2%) of ASD cases, a molecular basis could be determined and 7.2% of these carried copy number variations and/or chromosomal abnormalities, emphasizing the importance of detecting all forms of genetic variation as diagnostic and therapeutic targets in ASD.