De novo mutations of the voltage-gated sodium channel gene SCN8A have recently been recognized as a cause of epileptic encephalopathy, which is characterized by refractory seizures with developmental delay and cognitive disability. We previously described the heterozygous SCN8A missense mutation p.Asn1768Asp in a child with epileptic encephalopathy that included seizures, ataxia, and sudden unexpected death in epilepsy (SUDEP). The mutation results in increased persistent sodium current and hyperactivity of transfected neurons. We have characterized a knock-in mouse model expressing this dominant gain-of-function mutation to investigate the pathology of the altered channel in vivo. The mutant channel protein is stable in vivo. Heterozygous Scn8a(N1768D/+) mice exhibit seizures and SUDEP, confirming the causality of the de novo mutation in the proband. Using video/EEG analysis, we detect ictal discharges that coincide with convulsive seizures and myoclonic jerks. Prior to seizure onset, heterozygous mutants are not defective in motor learning or fear conditioning, but do exhibit mild impairment of motor coordination and social discrimination. Homozygous mutant mice exhibit earlier seizure onset than heterozygotes and more rapid progression to death. Analysis of the intermediate phenotype of functionally hemizygous Scn8a(N1768D/-) mice indicates that severity is increased by a double dose of mutant protein and reduced by the presence of wild-type protein. Scn8a(N1768D) mutant mice provide a model of epileptic encephalopathy that will be valuable for studying the in vivo effects of hyperactive Nav1.6 and the response to therapeutic interventions.

Comparisons of whole-genome sequences from ancient and contemporary samples have pointed to several instances of archaic admixture through interbreeding between the ancestors of modern non-Africans and now extinct hominids such as Neanderthals and Denisovans. One implication of these findings is that some adaptive features in contemporary humans may have entered the population via gene flow with archaic forms in Eurasia. Within Africa, fossil evidence suggests that anatomically modern humans (AMH) and various archaic forms coexisted for much of the last 200,000 yr; however, the absence of ancient DNA in Africa has limited our ability to make a direct comparison between archaic and modern human genomes. Here, we use statistical inference based on high coverage whole-genome data (greater than 60×) from contemporary African Pygmy hunter-gatherers as an alternative means to study the evolutionary history of the genus Homo. Using whole-genome simulations that consider demographic histories that include both isolation and gene flow with neighboring farming populations, our inference method rejects the hypothesis that the ancestors of AMH were genetically isolated in Africa, thus providing the first whole genome-level evidence of African archaic admixture. Our inferences also suggest a complex human evolutionary history in Africa, which involves at least a single admixture event from an unknown archaic population into the ancestors of AMH, likely within the last 30,000 yr.

Dravet syndrome (DS) is a rare, devastating form of childhood epilepsy that is often associated with mutations in the voltage-gated sodium channel gene, SCN1A. There is considerable variability in expressivity within families, as well as among individuals carrying the same primary mutation, suggesting that clinical outcome is modulated by variants at other genes. To identify modifier gene variants that contribute to clinical outcome, we sequenced the exomes of 22 individuals at both ends of a phenotype distribution (i.e., mild and severe cognitive condition). We controlled for variation associated with different mutation types by limiting inclusion to individuals with a de novo truncation mutation resulting in SCN1A haploinsufficiency. We performed tests aimed at identifying 1) single common variants that are enriched in either phenotypic group, 2) sets of common or rare variants aggregated in and around genes associated with clinical outcome, and 3) rare variants in 237 candidate genes associated with neuronal excitability. While our power to identify enrichment of a common variant in either phenotypic group is limited as a result of the rarity of mild phenotypes in individuals with SCN1A truncation variants, our top candidates did not map to functional regions of genes, or in genes that are known to be associated with neurological pathways. In contrast, we found a statistically-significant excess of rare variants predicted to be damaging and of small effect size in genes associated with neuronal excitability in severely affected individuals. A KCNQ2 variant previously associated with benign neonatal seizures is present in 3 of 12 individuals in the severe category. To compare our results with the healthy population, we performed a similar analysis on whole exome sequencing data from 70 Japanese individuals in the 1000 genomes project. Interestingly, the frequency of rare damaging variants in the same set of neuronal excitability genes in healthy individuals is nearly as high as in severely affected individuals. Rather than a single common gene/variant modifying clinical outcome in SCN1A-related epilepsies, our results point to the cumulative effect of rare variants with little to no measurable phenotypic effect (i.e., typical genetic background) unless present in combination with a disease-causing truncation mutation in SCN1A.

Whole exome sequencing (WES) of matched tumor-normal pairs in rare tumors has the potential to identify genome-wide mutations and copy number alterations (CNAs). We evaluated 27 rare cancer patients with tumor-normal matching by WES and tumor-only next generation sequencing (NGS) as a comparator. Our goal was to: 1) identify known and novel variants and CNAs in rare cancers with comparison to common cancers; 2) examine differences between germline and somatic variants and how that functionally impacts rare tumors; 3) detect and characterize alleles in biologically relevant genes-pathways that may be of clinical importance but not represented in classical cancer genes. We identified 3343 germline single nucleotide variants (SNVs) and small indel variants-1670 in oncogenes and 1673 in tumor suppressor genes-generating an average of 124 germline variants/case. The number of somatic SNVs and small indels detected in all cases was 523:306 in oncogenes and 217 in tumor suppressor genes. Of the germline variants, six were identified to be pathogenic or likely pathogenic. In the 27 analyzed rare cancer cases, CNAs are variable depending on tumor type, germline pathogenic variants are more common. Cell fate pathway mutations (e.g., Hippo, Notch, Wnt) dominate pathogenesis and double hit (mutation + CNV) represent ~18% cases.

Variants implicated in childhood epilepsy have been identified in all four voltage-gated sodium channels that initiate action potentials in the central nervous system. Previous research has focused on the functional effects of particular variants within the most studied of these channels (NaV1.1, NaV1.2 and NaV1.6); however, there have been few comparative studies across channels to infer the impact of mutations in patients with epilepsy. Here we compare patterns of variation in patient and public databases to test the hypothesis that regions of known functional significance within voltage-gated sodium (NaV) channels have an increased burden of deleterious variants. We assessed mutational burden in different regions of the Nav channels by (1) performing Fisher exact tests on odds ratios to infer excess variants in domains, segments, and loops of each channel in patient databases versus public "control" databases, and (2) comparing the cumulative distribution of variant sites along DNA sequences of each gene in patient and public databases (i.e., independent of protein structure). Patient variant density was concordant among channels in regions known to play a role in channel function, with statistically significant higher patient variant density in S4-S6 and DIII-DIV and an excess of public variants in SI-S3, DI-DII, DII-DIII. On the other hand, channel-specific patterns of patient burden were found in the NaV1.6 inactivation gate and NaV1.1 S5-S6 linkers, while NaV1.2 and NaV1.6 S4-S5 linkers and S5 segments shared patient variant patterns that contrasted with those in NaV1.1. These different patterns may reflect different roles played by the NaV1.6 inactivation gate in action potential propagation, and by NaV1.1 S5-S6 linkers in loss of function and haploinsufficiency. Interestingly, NaV1.2 and NaV1.6 both lack amino acid substitutions over significantly long stretches in both the patient and public databases suggesting that new mutations in these regions may cause embryonic lethality or a non-epileptic disease phenotype.

We investigated the genome of a 5-year-old male who presented with global developmental delay (motor, cognitive, and speech), hypotonia, possibly ataxia, and cerebellar hypoplasia of unknown origin. Whole genome sequencing (WGS) and mRNA sequencing (RNA-seq) were performed on a family having an affected proband, his unaffected parents, and maternal grandfather. To explore the molecular and functional consequences of the variant, we performed cell proliferation assays, quantitative real-time PCR (qRT-PCR) array, immunoblotting, calcium imaging, and neurite outgrowth experiments in SH-SY5Y neuroblastoma cells to compare the properties of the wild-type TATA-box-binding protein factor 1 (TAF1), deletion of TAF1, and TAF1 variant p.Ser1600Gly samples. The whole genome data identified several gene variants. However, the genome sequence data failed to implicate a candidate gene as many of the variants were of unknown significance. By combining genome sequence data with transcriptomic data, a probable candidate variant, p.Ser1600Gly, emerged in TAF1. Moreover, the RNA-seq data revealed a 90:10 extremely skewed X-chromosome inactivation (XCI) in the mother. Our results showed that neuronal ion channel genes were differentially expressed between TAF1 deletion and TAF1 variant p.Ser1600Gly cells, when compared with their respective controls, and that the TAF1 variant may impair neuronal differentiation and cell proliferation. Taken together, our data suggest that this novel variant in TAF1 plays a key role in the development of a recently described X-linked syndrome, TAF1 intellectual disability syndrome, and further extends our knowledge of a potential link between TAF1 deficiency and defects in neuronal cell function.