Search

PurposeExome and genome sequencing have proven to be effective tools for the diagnosis of neurodevelopmental disorders (NDDs), but large fractions of NDDs cannot be attributed to currently detectable genetic variation. This is likely, at least in part, a result of the fact that many genetic variants are difficult or impossible to detect through typical short-read sequencing approaches.MethodsHere, we describe a genomic analysis using Pacific Biosciences circular consensus sequencing (CCS) reads, which are both long (>10 kb) and accurate (>99% bp accuracy). We used CCS on six proband-parent trios with NDDs that were unexplained despite extensive testing, including genome sequencing with short reads.ResultsWe identified variants and created de novo assemblies in each trio, with global metrics indicating these data sets are more accurate and comprehensive than those provided by short-read data. In one proband, we identified a likely pathogenic (LP), de novo L1-mediated insertion in CDKL5 that results in duplication of exon 3, leading to a frameshift. In a second proband, we identified multiple large de novo structural variants, including insertion-translocations affecting DGKB and MLLT3, which we show disrupt MLLT3 transcript levels. We consider this extensive structural variation likely pathogenic.ConclusionThe breadth and quality of variant detection, coupled to finding variants of clinical and research interest in two of six probands with unexplained NDDs strongly support the value of long-read genome sequencing for understanding rare disease.

Dravet syndrome (DS) is a developmental and epileptic encephalopathy that results from mutations in the Nav1.1 sodium channel encoded by SCN1A. Most known DS-causing mutations are in coding regions of SCN1A, but we recently identified several disease-associated SCN1A mutations in intron 20 that are within or near to a cryptic and evolutionarily conserved “poison” exon, 20N, whose inclusion leads to transcript degradation. However, it is not clear how these intron 20 variants alter SCN1A transcript processing or DS pathophysiology in an organismal context, nor is it clear how exon 20N is regulated in a tissue-specific and developmental context. We address those questions here by generating an animal model of our index case, NM_006920.4(SCN1A):c.3969+2451G>C, using gene editing to create the orthologous mutation in laboratory mice. Scn1a heterozygous knock-in (+/KI) mice exhibited an ~50% reduction in brain Scn1a mRNA and Nav1.1 protein levels, together with characteristics observed in other DS mouse models, including premature mortality, seizures, and hyperactivity. In brain tissue from adult Scn1a +/+ animals, quantitative RT-PCR assays indicated that ~1% of Scn1a mRNA included exon 20N, while brain tissue from Scn1a +/KI mice exhibited an ~5-fold increase in the extent of exon 20N inclusion. We investigated the extent of exon 20N inclusion in brain during normal fetal development in RNA-seq data and discovered that levels of inclusion were ~70% at E14.5, declining progressively to ~10% postnatally. A similar pattern exists for the homologous sodium channel Nav1.6, encoded by Scn8a. For both genes, there is an inverse relationship between the level of functional transcript and the extent of poison exon inclusion. Taken together, our findings suggest that poison exon usage by Scn1a and Scn8a is a strategy to regulate channel expression during normal brain development, and that mutations recapitulating a fetal-like pattern of splicing cause reduced channel expression and epileptic encephalopathy.Author SummaryDravet syndrome (DS) is a neurological disorder affecting approximately 1:15,700 Americans[1]. While most patients have a mutation in the SCN1A gene encoding Nav1.1 sodium channels, about 20% do not have a mutation identified by exome sequencing. Recently, we identified variants in intron 20N, a noncoding region of SCN1A, in some DS patients [2]. We predicted that these variants alter SCN1A transcript processing, decrease Nav1.1 function, and lead to DS pathophysiology via inclusion of exon 20N, a “poison” exon that leads to a premature stop codon. In this study, we generated a knock-in mouse model, Scn1a+/KI, of one of these variants, NM_006920.4(SCN1A):c.3969+2451G>C, which resides in a genomic region that is extremely conserved across vertebrate species. We found that Scn1a+/KI mice have reduced levels of Scn1a transcript and Nav1.1 protein and develop DS-related phenotypes. Consistent with the poison exon hypothesis, transcripts from brains of Scn1a+/KI mice showed elevated rates of Scn1a exon 20N inclusion. Since Scn1a expression in the brain is regulated developmentally, we next explored the developmental relationship between exon 20N inclusion and Scn1a expression. During normal embryogenesis, when Scn1a expression was low, exon 20N inclusion was high; postnatally, as Scn1a expression increased, there was a corresponding decrease in exon 20N usage. Expression of another voltage-gated sodium channel transcript, Scn8a (Nav1.6), was similarly regulated, with inclusion of a poison exon termed as 18N early in development when Scn8a expression was low, followed by a postnatal decrease in exon 18N inclusion and corresponding increase in Scn8a expression. Together, these data demonstrate that poison exon inclusion is a conserved mechanism to control sodium channel expression in the brain, and that an intronic mutation that disrupts the normal developmental regulation of poison exon inclusion leads to reduced Nav1.1 and DS pathophysiology.