Aberrant activation or inhibition of potassium (K+) currents across the plasma membrane of cells has been causally linked to altered neurotransmission, cardiac arrhythmias, endocrine dysfunction, and (more rarely) perturbed developmental processes. The K+ channel subfamily K member 4 (KCNK4), also known as TRAAK (TWIK-related arachidonic acid-stimulated K+ channel), belongs to the mechano-gated ion channels of the TRAAK/TREK subfamily of two-pore-domain (K2P) K+ channels. While K2P channels are well known to contribute to the resting membrane potential and cellular excitability, their involvement in pathophysiological processes remains largely uncharacterized. We report that de novo missense mutations in KCNK4 cause a recognizable syndrome with a distinctive facial gestalt, for which we propose the acronym FHEIG (facial dysmorphism, hypertrichosis, epilepsy, intellectual disability/developmental delay, and gingival overgrowth). Patch-clamp analyses documented a significant gain of function of the identified KCNK4 channel mutants basally and impaired sensitivity to mechanical stimulation and arachidonic acid. Co-expression experiments indicated a dominant behavior of the disease-causing mutations. Molecular dynamics simulations consistently indicated that mutations favor sealing of the lateral intramembrane fenestration that has been proposed to negatively control K+ flow by allowing lipid access to the central cavity of the channel. Overall, our findings illustrate the pleiotropic effect of dysregulated KCNK4 function and provide support to the hypothesis of a gating mechanism based on the lateral fenestrations of K2P channels.

Identification of over 500 epigenetic regulators in humans raises an interesting question regarding how chromatin dysregulation contributes to different diseases. Bromodomain and PHD finger-containing protein 1 (BRPF1) is a multivalent chromatin regulator possessing three histone-binding domains, one non-specific DNA-binding module, and several motifs for interacting with and activating three lysine acetyltransferases. Genetic analyses of fish brpf1 and mouse Brpf1 have uncovered an important role in skeletal, hematopoietic, and brain development, but it remains unclear how BRPF1 is linked to human development and disease. Here, we describe an intellectual disability disorder in ten individuals with inherited or de novo monoallelic BRPF1 mutations. Symptoms include infantile hypotonia, global developmental delay, intellectual disability, expressive language impairment, and facial dysmorphisms. Central nervous system and spinal abnormalities are also seen in some individuals. These clinical features overlap with but are not identical to those reported for persons with KAT6A or KAT6B mutations, suggesting that BRPF1 targets these two acetyltransferases and additional partners in humans. Functional assays showed that the resulting BRPF1 variants are pathogenic and impair acetylation of histone H3 at lysine 23, an abundant but poorly characterized epigenetic mark. We also found a similar deficiency in different lines of Brpf1-knockout mice. These data indicate that aberrations in the chromatin regulator gene BRPF1 cause histone H3 acetylation deficiency and a previously unrecognized intellectual disability syndrome.

We report 15 individuals with de novo pathogenic variants in WDR26. Eleven of the individuals carry loss-of-function mutations, and four harbor missense substitutions. These 15 individuals comprise ten females and five males, and all have intellectual disability with delayed speech, a history of febrile and/or non-febrile seizures, and a wide-based, spastic, and/or stiff-legged gait. These subjects share a set of common facial features that include a prominent maxilla and upper lip that readily reveal the upper gingiva, widely spaced teeth, and a broad nasal tip. Together, these features comprise a recognizable facial phenotype. We compared these features with those of chromosome 1q41q42 microdeletion syndrome, which typically contains WDR26, and noted that clinical features are consistent between the two subsets, suggesting that haploinsufficiency of WDR26 contributes to the pathology of 1q41q42 microdeletion syndrome. Consistent with this, WDR26 loss-of-function single-nucleotide mutations identified in these subjects lead to nonsense-mediated decay with subsequent reduction of RNA expression and protein levels. We derived a structural model of WDR26 and note that missense variants identified in these individuals localize to highly conserved residues of this WD-40-repeat-containing protein. Given that WDR26 mutations have been identified in ∼1 in 2,000 of subjects in our clinical cohorts and that WDR26 might be poorly annotated in exome variant-interpretation pipelines, we would anticipate that this disorder could be more common than currently appreciated.

KCNQ5 is a highly conserved gene encoding an important channel for neuronal function; it is widely expressed in the brain and generates M-type current. Exome sequencing identified de novo heterozygous missense mutations in four probands with intellectual disability, abnormal neurological findings, and treatment-resistant epilepsy (in two of four). Comprehensive analysis of this potassium channel for the four variants expressed in frog oocytes revealed shifts in the voltage dependence of activation, including altered activation and deactivation kinetics. Specifically, both loss-of-function and gain-of-function KCNQ5 mutations, associated with increased excitability and decreased repolarization reserve, lead to pathophysiology.

Large-scale efforts like the ENCODE Project have made tremendous progress in cataloging the genomic binding patterns of DNA-associated proteins (DAPs), such as transcription factors (TFs). However, most chromatin immunoprecipitation-sequencing (ChIP-seq) analyses have focused on a few immortalized cell lines whose activities and physiology differ in important ways from endogenous cells and tissues. Consequently, binding data from primary human tissue are essential to improving our understanding of in vivo gene regulation. Here, we identify and analyze more than 440,000 binding sites using ChIP-seq data for 20 DAPs in two human liver tissue samples. We integrated binding data with transcriptome and phased WGS data to investigate allelic DAP interactions and the impact of heterozygous sequence variation on the expression of neighboring genes. Our tissue-based data set exhibits binding patterns more consistent with liver biology than cell lines, and we describe uses of these data to better prioritize impactful noncoding variation. Collectively, our rich data set offers novel insights into genome function in human liver tissue and provides a valuable resource for assessing disease-related disruptions.

The discovery of expression quantitative trait loci ("eQTLs") can help to unravel genetic contributions to complex traits. We identified genetic determinants of human liver gene expression variation using two independent collections of primary tissue profiled with Agilent (n = 206) and Illumina (n = 60) expression arrays and Illumina SNP genotyping (550K), and we also incorporated data from a published study (n = 266). We found that ∼30% of SNP-expression correlations in one study failed to replicate in either of the others, even at thresholds yielding high reproducibility in simulations, and we quantified numerous factors affecting reproducibility. Our data suggest that drug exposure, clinical descriptors, and unknown factors associated with tissue ascertainment and analysis have substantial effects on gene expression and that controlling for hidden confounding variables significantly increases replication rate. Furthermore, we found that reproducible eQTL SNPs were heavily enriched near gene starts and ends, and subsequently resequenced the promoters and 3'UTRs for 14 genes and tested the identified haplotypes using luciferase assays. For three genes, significant haplotype-specific in vitro functional differences correlated directly with expression levels, suggesting that many bona fide eQTLs result from functional variants that can be mechanistically isolated in a high-throughput fashion. Finally, given our study design, we were able to discover and validate hundreds of liver eQTLs. Many of these relate directly to complex traits for which liver-specific analyses are likely to be relevant, and we identified dozens of potential connections with disease-associated loci. These included previously characterized eQTL contributors to diabetes, drug response, and lipid levels, and they suggest novel candidates such as a role for NOD2 expression in leprosy risk and C2orf43 in prostate cancer. In general, the work presented here will be valuable for future efforts to precisely identify and functionally characterize genetic contributions to a variety of complex traits.

Computational efforts to identify functional elements within genomes leverage comparative sequence information by looking for regions that exhibit evidence of selective constraint. One way of detecting constrained elements is to follow a bottom-up approach by computing constraint scores for individual positions of a multiple alignment and then defining constrained elements as segments of contiguous, highly scoring nucleotide positions. Here we present GERP++, a new tool that uses maximum likelihood evolutionary rate estimation for position-specific scoring and, in contrast to previous bottom-up methods, a novel dynamic programming approach to subsequently define constrained elements. GERP++ evaluates a richer set of candidate element breakpoints and ranks them based on statistical significance, eliminating the need for biased heuristic extension techniques. Using GERP++ we identify over 1.3 million constrained elements spanning over 7% of the human genome. We predict a higher fraction than earlier estimates largely due to the annotation of longer constrained elements, which improves one to one correspondence between predicted elements with known functional sequences. GERP++ is an efficient and effective tool to provide both nucleotide- and element-level constraint scores within deep multiple sequence alignments.

Most studies of genomic disorders have focused on patients with cognitive disability and/or peripheral nervous system defects. In an effort to broaden the phenotypic spectrum of this disease model, we assessed 155 autopsy samples from fetuses with well-defined developmental pathologies in regions predisposed to recurrent rearrangement, by array-based comparative genomic hybridization. We found that 6% of fetal material showed evidence of microdeletion or microduplication, including three independent events that likely resulted from unequal crossing-over between segmental duplications. One of the microdeletions, identified in a fetus with multicystic dysplastic kidneys, encompasses the TCF2 gene on 17q12, previously shown to be mutated in maturity-onset diabetes, as well as in a subset of pediatric renal abnormalities. Fine-scale mapping of the breakpoints in different patient cohorts revealed a recurrent 1.5-Mb de novo deletion in individuals with phenotypes that ranged from congenital renal abnormalities to maturity-onset diabetes of the young type 5. We also identified the reciprocal duplication, which appears to be enriched in samples from patients with epilepsy. We describe the first example of a recurrent genomic disorder associated with diabetes.

We report the generation and analysis of functional data from multiple, diverse experiments performed on a targeted 1% of the human genome as part of the pilot phase of the ENCODE Project. These data have been further integrated and augmented by a number of evolutionary and computational analyses. Together, our results advance the collective knowledge about human genome function in several major areas. First, our studies provide convincing evidence that the genome is pervasively transcribed, such that the majority of its bases can be found in primary transcripts, including non-protein-coding transcripts, and those that extensively overlap one another. Second, systematic examination of transcriptional regulation has yielded new understanding about transcription start sites, including their relationship to specific regulatory sequences and features of chromatin accessibility and histone modification. Third, a more sophisticated view of chromatin structure has emerged, including its inter-relationship with DNA replication and transcriptional regulation. Finally, integration of these new sources of information, in particular with respect to mammalian evolution based on inter- and intra-species sequence comparisons, has yielded new mechanistic and evolutionary insights concerning the functional landscape of the human genome. Together, these studies are defining a path for pursuit of a more comprehensive characterization of human genome function.

Histones mediate dynamic packaging of nuclear DNA in chromatin, a process that is precisely controlled to guarantee efficient compaction of the genome and proper chromosomal segregation during cell division and to accomplish DNA replication, transcription, and repair. Due to the important structural and regulatory roles played by histones, it is not surprising that histone functional dysregulation or aberrant levels of histones can have severe consequences for multiple cellular processes and ultimately might affect development or contribute to cell transformation. Recently, germline frameshift mutations involving the C-terminal tail of HIST1H1E, which is a widely expressed member of the linker histone family and facilitates higher-order chromatin folding, have been causally linked to an as-yet poorly defined syndrome that includes intellectual disability. We report that these mutations result in stable proteins that reside in the nucleus, bind to chromatin, disrupt proper compaction of DNA, and are associated with a specific methylation pattern. Cells expressing these mutant proteins have a dramatically reduced proliferation rate and competence, hardly enter into the S phase, and undergo accelerated senescence. Remarkably, clinical assessment of a relatively large cohort of subjects sharing these mutations revealed a premature aging phenotype as a previously unrecognized feature of the disorder. Our findings identify a direct link between aberrant chromatin remodeling, cellular senescence, and accelerated aging.

SOX6 belongs to a family of 20 SRY-related HMG-box-containing (SOX) genes that encode transcription factors controlling cell fate and differentiation in many developmental and adult processes. For SOX6, these processes include, but are not limited to, neurogenesis and skeletogenesis. Variants in half of the SOX genes have been shown to cause severe developmental and adult syndromes, referred to as SOXopathies. We here provide evidence that SOX6 variants also cause a SOXopathy. Using clinical and genetic data, we identify 19 individuals harboring various types of SOX6 alterations and exhibiting developmental delay and/or intellectual disability; the individuals are from 17 unrelated families. Additional, inconstant features include attention-deficit/hyperactivity disorder (ADHD), autism, mild facial dysmorphism, craniosynostosis, and multiple osteochondromas. All variants are heterozygous. Fourteen are de novo, one is inherited from a mosaic father, and four offspring from two families have a paternally inherited variant. Intragenic microdeletions, balanced structural rearrangements, frameshifts, and nonsense variants are predicted to inactivate the SOX6 variant allele. Four missense variants occur in residues and protein regions highly conserved evolutionarily. These variants are not detected in the gnomAD control cohort, and the amino acid substitutions are predicted to be damaging. Two of these variants are located in the HMG domain and abolish SOX6 transcriptional activity in vitro. No clear genotype-phenotype correlations are found. Taken together, these findings concur that SOX6 haploinsufficiency leads to a neurodevelopmental SOXopathy that often includes ADHD and abnormal skeletal and other features.

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 is predicted to lead to transcript degradation. However, it is not clear how these intron 20 variants alter SCN1A expression 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.

Neurodevelopmental disorders are highly heterogenous conditions resulting from abnormalities of brain architecture and/or function. FBXW7 (F-box and WD-repeat-domain-containing 7), a recognized developmental regulator and tumor suppressor, has been shown to regulate cell-cycle progression and cell growth and survival by targeting substrates including CYCLIN E1/2 and NOTCH for degradation via the ubiquitin proteasome system. We used a genotype-first approach and global data-sharing platforms to identify 35 individuals harboring de novo and inherited FBXW7 germline monoallelic chromosomal deletions and nonsense, frameshift, splice-site, and missense variants associated with a neurodevelopmental syndrome. The FBXW7 neurodevelopmental syndrome is distinguished by global developmental delay, borderline to severe intellectual disability, hypotonia, and gastrointestinal issues. Brain imaging detailed variable underlying structural abnormalities affecting the cerebellum, corpus collosum, and white matter. A crystal-structure model of FBXW7 predicted that missense variants were clustered at the substrate-binding surface of the WD40 domain and that these might reduce FBXW7 substrate binding affinity. Expression of recombinant FBXW7 missense variants in cultured cells demonstrated impaired CYCLIN E1 and CYCLIN E2 turnover. Pan-neuronal knockdown of the Drosophila ortholog, archipelago, impaired learning and neuronal function. Collectively, the data presented herein provide compelling evidence of an F-Box protein-related, phenotypically variable neurodevelopmental disorder associated with monoallelic variants in FBXW7.

Exome 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. Here, 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. We identified variants and created de novo assemblies in each trio, with global metrics indicating these datasets 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. The breadth and quality of variant detection, coupled to finding variants of clinical and research interest in two of six probands with unexplained NDDs, support the hypothesis that long-read genome sequencing can substantially improve rare disease genetic discovery rates.

Located in the critical 1p36 microdeletion region, the chromodomain helicase DNA-binding protein 5 (CHD5) gene encodes a subunit of the nucleosome remodeling and deacetylation (NuRD) complex required for neuronal development. Pathogenic variants in six of nine chromodomain (CHD) genes cause autosomal dominant neurodevelopmental disorders, while CHD5-related disorders are still unknown. Thanks to GeneMatcher and international collaborations, we assembled a cohort of 16 unrelated individuals harboring heterozygous CHD5 variants, all identified by exome sequencing. Twelve patients had de novo CHD5 variants, including ten missense and two splice site variants. Three familial cases had nonsense or missense variants segregating with speech delay, learning disabilities, and/or craniosynostosis. One patient carried a frameshift variant of unknown inheritance due to unavailability of the father. The most common clinical features included language deficits (81%), behavioral symptoms (69%), intellectual disability (64%), epilepsy (62%), and motor delay (56%). Epilepsy types were variable, with West syndrome observed in three patients, generalized tonic-clonic seizures in two, and other subtypes observed in one individual each. Our findings suggest that, in line with other CHD-related disorders, heterozygous CHD5 variants are associated with a variable neurodevelopmental syndrome that includes intellectual disability with speech delay, epilepsy, and behavioral problems as main features.

Calcium (Ca2+) is a universal second messenger involved in synaptogenesis and cell survival; consequently, its regulation is important for neurons. ATPase plasma membrane Ca2+ transporting 1 (ATP2B1) belongs to the family of ATP-driven calmodulin-dependent Ca2+ pumps that participate in the regulation of intracellular free Ca2+. Here, we clinically describe a cohort of 12 unrelated individuals with variants in ATP2B1 and an overlapping phenotype of mild to moderate global development delay. Additional common symptoms include autism, seizures, and distal limb abnormalities. Nine probands harbor missense variants, seven of which were in specific functional domains, and three individuals have nonsense variants. 3D structural protein modeling suggested that the variants have a destabilizing effect on the protein. We performed Ca2+ imaging after introducing all nine missense variants in transfected HEK293 cells and showed that all variants lead to a significant decrease in Ca2+ export capacity compared with the wild-type construct, thus proving their pathogenicity. Furthermore, we observed for the same variant set an incorrect intracellular localization of ATP2B1. The genetic findings and the overlapping phenotype of the probands as well as the functional analyses imply that de novo variants in ATP2B1 lead to a monogenic form of neurodevelopmental disorder.

Massively parallel reporter assays (MPRAs) are useful tools to characterize regulatory elements in human genomes. An aspect of MPRAs that is not typically the focus of analysis is their intrinsic ability to differentiate activity levels for a given sequence element when placed in both of its possible orientations relative to the reporter construct. Here, we describe pervasive strand asymmetry of MPRA signals in data sets from multiple reporter configurations in both published and newly reported data. These effects are reproducible across different cell types and in different treatments within a cell type and are observed both within and outside of annotated regulatory elements. From elements in gene bodies, MPRA strand asymmetry favors the sense strand, suggesting that function related to endogenous transcription is driving the phenomenon. Similarly, we find that within Alu mobile element insertions, strand asymmetry favors the transcribed strand of the ancestral retrotransposon. The effect is consistent across the multiplicity of Alu elements in human genomes and is more pronounced in less diverged Alu elements. We find sequence features driving MPRA strand asymmetry and show its prediction from sequence alone. We see some evidence for RNA stabilization and transcriptional activation mechanisms and hypothesize that the effect is driven by natural selection favoring efficient transcription. Our results indicate that strand asymmetry is a pervasive and reproducible feature in MPRA data. More importantly, the fact that MPRA asymmetry favors naturally transcribed strands suggests that it stems from preserved biological functions that have a substantial, global impact on gene and genome evolution.

Current methods for annotating and interpreting human genetic variation tend to exploit a single information type (for example, conservation) and/or are restricted in scope (for example, to missense changes). Here we describe Combined Annotation-Dependent Depletion (CADD), a method for objectively integrating many diverse annotations into a single measure (C score) for each variant. We implement CADD as a support vector machine trained to differentiate 14.7 million high-frequency human-derived alleles from 14.7 million simulated variants. We precompute C scores for all 8.6 billion possible human single-nucleotide variants and enable scoring of short insertions-deletions. C scores correlate with allelic diversity, annotations of functionality, pathogenicity, disease severity, experimentally measured regulatory effects and complex trait associations, and they highly rank known pathogenic variants within individual genomes. The ability of CADD to prioritize functional, deleterious and pathogenic variants across many functional categories, effect sizes and genetic architectures is unmatched by any current single-annotation method.

The identification of signatures of natural selection has long been used as an approach to understanding the unique features of any given species. Genes within segmental duplications are overlooked in most studies of selection due to the limitations of draft nonhuman genome assemblies and to the methodological reliance on accurate gene trees, which are difficult to obtain for duplicated genes.

To understand the genetic heterogeneity underlying developmental delay, we compared copy number variants (CNVs) in 15,767 children with intellectual disability and various congenital defects (cases) to CNVs in 8,329 unaffected adult controls. We estimate that ∼14.2% of disease in these children is caused by CNVs >400 kb. We observed a greater enrichment of CNVs in individuals with craniofacial anomalies and cardiovascular defects compared to those with epilepsy or autism. We identified 59 pathogenic CNVs, including 14 new or previously weakly supported candidates, refined the critical interval for several genomic disorders, such as the 17q21.31 microdeletion syndrome, and identified 940 candidate dosage-sensitive genes. We also developed methods to opportunistically discover small, disruptive CNVs within the large and growing diagnostic array datasets. This evolving CNV morbidity map, combined with exome and genome sequencing, will be critical for deciphering the genetic basis of developmental delay, intellectual disability and autism spectrum disorders.