The 17p13.1 microdeletion syndrome is a recently described genomic disorder with a core clinical phenotype of intellectual disability, poor to absent speech, dysmorphic features, and a constellation of more variable clinical features, most prominently microcephaly. We identified five subjects with copy-number variants (CNVs) on 17p13.1 for whom we performed detailed clinical and molecular studies. Breakpoint mapping and retrospective analysis of published cases refined the smallest region of overlap (SRO) for microcephaly to a genomic interval containing nine genes. Dissection of this phenotype in zebrafish embryos revealed a complex genetic architecture: dosage perturbation of four genes (ASGR1, ACADVL, DVL2, and GABARAP) impeded neurodevelopment and decreased dosage of the same loci caused a reduced mitotic index in vitro. Moreover, epistatic analyses in vivo showed that dosage perturbations of discrete gene pairings induce microcephaly. Taken together, these studies support a model in which concomitant dosage perturbation of multiple genes within the CNV drive the microcephaly and possibly other neurodevelopmental phenotypes associated with rearrangements in the 17p13.1 SRO.

Chromodomain helicase DNA-binding protein 4 (CHD4) is an ATP-dependent chromatin remodeler involved in epigenetic regulation of gene transcription, DNA repair, and cell cycle progression. Also known as Mi2β, CHD4 is an integral subunit of a well-characterized histone deacetylase complex. Here we report five individuals with de novo missense substitutions in CHD4 identified through whole-exome sequencing and web-based gene matching. These individuals have overlapping phenotypes including developmental delay, intellectual disability, hearing loss, macrocephaly, distinct facial dysmorphisms, palatal abnormalities, ventriculomegaly, and hypogonadism as well as additional findings such as bone fusions. The variants, c.3380G>A (p.Arg1127Gln), c.3443G>T (p.Trp1148Leu), c.3518G>T (p.Arg1173Leu), and c.3008G>A, (p.Gly1003Asp) (GenBank: NM_001273.3), affect evolutionarily highly conserved residues and are predicted to be deleterious. Previous studies in yeast showed the equivalent Arg1127 and Trp1148 residues to be crucial for SNF2 function. Furthermore, mutations in the same positions were reported in malignant tumors, and a de novo missense substitution in an equivalent arginine residue in the C-terminal helicase domain of SMARCA4 is associated with Coffin Siris syndrome. Cell-based studies of the p.Arg1127Gln and p.Arg1173Leu mutants demonstrate normal localization to the nucleus and HDAC1 interaction. Based on these findings, the mutations potentially alter the complex activity but not its formation. This report provides evidence for the role of CHD4 in human development and expands an increasingly recognized group of Mendelian disorders involving chromatin remodeling and modification.

5q31.3 microdeletion syndrome is characterized by neonatal hypotonia, encephalopathy with or without epilepsy, and severe developmental delay, and the minimal critical deletion interval harbors three genes. We describe 11 individuals with clinical features of 5q31.3 microdeletion syndrome and de novo mutations in PURA, encoding transcriptional activator protein Pur-α, within the critical region. These data implicate causative PURA mutations responsible for the severe neurological phenotypes observed in this syndrome.

We investigated complex genomic rearrangements (CGRs) consisting of triplication copy-number variants (CNVs) that were accompanied by extended regions of copy-number-neutral absence of heterozygosity (AOH) in subjects with multiple congenital abnormalities. Molecular analyses provided observational evidence that in humans, post-zygotically generated CGRs can lead to regional uniparental disomy (UPD) due to template switches between homologs versus sister chromatids by using microhomology to prime DNA replication-a prediction of the replicative repair model, MMBIR. Our findings suggest that replication-based mechanisms might underlie the formation of diverse types of genomic alterations (CGRs and AOH) implicated in constitutional disorders.

Genitopatellar syndrome (GPS) is a skeletal dysplasia with cerebral and genital anomalies for which the molecular basis has not yet been determined. By exome sequencing, we found de novo heterozygous truncating mutations in KAT6B (lysine acetyltransferase 6B, formerly known as MYST4 and MORF) in three subjects; then by Sanger sequencing of KAT6B, we found similar mutations in three additional subjects. The mutant transcripts do not undergo nonsense-mediated decay in cells from subjects with GPS. In addition, human pathological analyses and mouse expression studies point to systemic roles of KAT6B in controlling organismal growth and development. Myst4 (the mouse orthologous gene) is expressed in mouse tissues corresponding to those affected by GPS. Phenotypic differences and similarities between GPS, the Say-Barber-Biesecker variant of Ohdo syndrome (caused by different mutations of KAT6B), and Rubinstein-Taybi syndrome (caused by mutations in other histone acetyltransferases) are discussed. Together, the data support an epigenetic dysregulation of the limb, brain, and genital developmental programs.

Nuclear deubiquitinase BAP1 (BRCA1-associated protein 1) is a core component of multiprotein complexes that promote transcription by reversing the ubiquitination of histone 2A (H2A). BAP1 is a tumor suppressor whose germline loss-of-function variants predispose to cancer. To our knowledge, there are very rare examples of different germline variants in the same gene causing either a neurodevelopmental disorder (NDD) or a tumor predisposition syndrome. Here, we report a series of 11 de novo germline heterozygous missense BAP1 variants associated with a rare syndromic NDD. Functional analysis showed that most of the variants cannot rescue the consequences of BAP1 inactivation, suggesting a loss-of-function mechanism. In T cells isolated from two affected children, H2A deubiquitination was impaired. In matching peripheral blood mononuclear cells, histone H3 K27 acetylation ChIP-seq indicated that these BAP1 variants induced genome-wide chromatin state alterations, with enrichment for regulatory regions surrounding genes of the ubiquitin-proteasome system (UPS). Altogether, these results define a clinical syndrome caused by rare germline missense BAP1 variants that alter chromatin remodeling through abnormal histone ubiquitination and lead to transcriptional dysregulation of developmental genes.

Calcium/calmodulin-dependent protein kinase II (CAMK2) is one of the first proteins shown to be essential for normal learning and synaptic plasticity in mice, but its requirement for human brain development has not yet been established. Through a multi-center collaborative study based on a whole-exome sequencing approach, we identified 19 exceedingly rare de novo CAMK2A or CAMK2B variants in 24 unrelated individuals with intellectual disability. Variants were assessed for their effect on CAMK2 function and on neuronal migration. For both CAMK2A and CAMK2B, we identified mutations that decreased or increased CAMK2 auto-phosphorylation at Thr286/Thr287. We further found that all mutations affecting auto-phosphorylation also affected neuronal migration, highlighting the importance of tightly regulated CAMK2 auto-phosphorylation in neuronal function and neurodevelopment. Our data establish the importance of CAMK2A and CAMK2B and their auto-phosphorylation in human brain function and expand the phenotypic spectrum of the disorders caused by variants in key players of the glutamatergic signaling pathway.

To date, mutations in 15 actin- or microtubule-associated genes have been associated with the cortical malformation lissencephaly and variable brainstem hypoplasia. During a multicenter review, we recognized a rare lissencephaly variant with a complex brainstem malformation in three unrelated children. We searched our large brain-malformation databases and found another five children with this malformation (as well as one with a less severe variant), analyzed available whole-exome or -genome sequencing data, and tested ciliogenesis in two affected individuals. The brain malformation comprised posterior predominant lissencephaly and midline crossing defects consisting of absent anterior commissure and a striking W-shaped brainstem malformation caused by small or absent pontine crossing fibers. We discovered heterozygous de novo missense variants or an in-frame deletion involving highly conserved zinc-binding residues within the GAR domain of MACF1 in the first eight subjects. We studied cilium formation and found a higher proportion of mutant cells with short cilia than of control cells with short cilia. A ninth child had similar lissencephaly but only subtle brainstem dysplasia associated with a heterozygous de novo missense variant in the spectrin repeat domain of MACF1. Thus, we report variants of the microtubule-binding GAR domain of MACF1 as the cause of a distinctive and most likely pathognomonic brain malformation. A gain-of-function or dominant-negative mechanism appears likely given that many heterozygous mutations leading to protein truncation are included in the ExAC Browser. However, three de novo variants in MACF1 have been observed in large schizophrenia cohorts.

TRAF7 is a multi-functional protein involved in diverse signaling pathways and cellular processes. The phenotypic consequence of germline TRAF7 variants remains unclear. Here we report missense variants in TRAF7 in seven unrelated individuals referred for clinical exome sequencing. The seven individuals share substantial phenotypic overlap, with developmental delay, congenital heart defects, limb and digital anomalies, and dysmorphic features emerging as key unifying features. The identified variants are de novo in six individuals and comprise four distinct missense changes, including a c.1964G>A (p.Arg655Gln) variant that is recurrent in four individuals. These variants affect evolutionarily conserved amino acids and are located in key functional domains. Gene-specific mutation rate analysis showed that the occurrence of the de novo variants in TRAF7 (p = 2.6 × 10-3) and the recurrent de novo c.1964G>A (p.Arg655Gln) variant (p = 1.9 × 10-8) in our exome cohort was unlikely to have occurred by chance. In vitro analyses of the observed TRAF7 mutations showed reduced ERK1/2 phosphorylation. Our findings suggest that missense mutations in TRAF7 are associated with a multisystem disorder and provide evidence of a role for TRAF7 in human development.

Mutations in several known or putative glycosyltransferases cause glycosylation defects in α-dystroglycan (α-DG), an integral component of the dystrophin glycoprotein complex. The hypoglycosylation reduces the ability of α-DG to bind laminin and other extracellular matrix ligands and is responsible for the pathogenesis of an inherited subset of muscular dystrophies known as the dystroglycanopathies. By exome and Sanger sequencing we identified two individuals affected by a dystroglycanopathy with mutations in β-1,3-N-acetylgalactosaminyltransferase 2 (B3GALNT2). B3GALNT2 transfers N-acetyl galactosamine (GalNAc) in a β-1,3 linkage to N-acetyl glucosamine (GlcNAc). A subsequent study of a separate cohort of individuals identified recessive mutations in four additional cases that were all affected by dystroglycanopathy with structural brain involvement. We show that functional dystroglycan glycosylation was reduced in the fibroblasts and muscle (when available) of these individuals via flow cytometry, immunoblotting, and immunocytochemistry. B3GALNT2 localized to the endoplasmic reticulum, and this localization was perturbed by some of the missense mutations identified. Moreover, knockdown of b3galnt2 in zebrafish recapitulated the human congenital muscular dystrophy phenotype with reduced motility, brain abnormalities, and disordered muscle fibers with evidence of damage to both the myosepta and the sarcolemma. Functional dystroglycan glycosylation was also reduced in the b3galnt2 knockdown zebrafish embryos. Together these results demonstrate a role for B3GALNT2 in the glycosylation of α-DG and show that B3GALNT2 mutations can cause dystroglycanopathy with muscle and brain involvement.

Bromodomain PHD finger transcription factor (BPTF) is the largest subunit of nucleosome remodeling factor (NURF), a member of the ISWI chromatin-remodeling complex. However, the clinical consequences of disruption of this complex remain largely uncharacterized. BPTF is required for anterior-posterior axis formation of the mouse embryo and was shown to promote posterior neuroectodermal fate by enhancing Smad2-activated wnt8 expression in zebrafish. Here, we report eight loss-of-function and two missense variants (eight de novo and two of unknown origin) in BPTF on 17q24.2. The BPTF variants were found in unrelated individuals aged between 2.1 and 13 years, who manifest variable degrees of developmental delay/intellectual disability (10/10), speech delay (10/10), postnatal microcephaly (7/9), and dysmorphic features (9/10). Using CRISPR-Cas9 genome editing of bptf in zebrafish to induce a loss of gene function, we observed a significant reduction in head size of F0 mutants compared to control larvae. Terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) and phospho-histone H3 (PH3) staining to assess apoptosis and cell proliferation, respectively, showed a significant increase in cell death in F0 mutants compared to controls. Additionally, we observed a substantial increase of the ceratohyal angle of the craniofacial skeleton in bptf F0 mutants, indicating abnormal craniofacial patterning. Taken together, our data demonstrate the pathogenic role of BPTF haploinsufficiency in syndromic neurodevelopmental anomalies and extend the clinical spectrum of human disorders caused by ablation of chromatin remodeling complexes.

Diagnosis at the edges of our knowledge calls upon clinicians to be data driven, cross-disciplinary, and collaborative in unprecedented ways. Exact disease recognition, an element of the concept of precision in medicine, requires new infrastructure that spans geography, institutional boundaries, and the divide between clinical care and research. The National Institutes of Health (NIH) Common Fund supports the Undiagnosed Diseases Network (UDN) as an exemplar of this model of precise diagnosis. Its goals are to forge a strategy to accelerate the diagnosis of rare or previously unrecognized diseases, to improve recommendations for clinical management, and to advance research, especially into disease mechanisms. The network will achieve these objectives by evaluating patients with undiagnosed diseases, fostering a breadth of expert collaborations, determining best practices for translating the strategy into medical centers nationwide, and sharing findings, data, specimens, and approaches with the scientific and medical communities. Building the UDN has already brought insights to human and medical geneticists. The initial focus has been on data sharing, establishing common protocols for institutional review boards and data sharing, creating protocols for referring and evaluating patients, and providing DNA sequencing, metabolomic analysis, and functional studies in model organisms. By extending this precision diagnostic model nationally, we strive to meld clinical and research objectives, improve patient outcomes, and contribute to medical science.

Whole-exome sequencing of 13 individuals with developmental delay commonly accompanied by abnormal muscle tone and seizures identified de novo missense mutations enriched within a sub-region of GNB1, a gene encoding the guanine nucleotide-binding protein subunit beta-1, Gβ. These 13 individuals were identified among a base of 5,855 individuals recruited for various undiagnosed genetic disorders. The probability of observing 13 or more de novo mutations by chance among 5,855 individuals is very low (p = 7.1 × 10(-21)), implicating GNB1 as a genome-wide-significant disease-associated gene. The majority of these 13 mutations affect known Gβ binding sites, which suggests that a likely disease mechanism is through the disruption of the protein interface required for Gα-Gβγ interaction (resulting in a constitutively active Gβγ) or through the disruption of residues relevant for interaction between Gβγ and certain downstream effectors (resulting in reduced interaction with the effectors). Strikingly, 8 of the 13 individuals recruited here for a neurodevelopmental disorder have a germline de novo GNB1 mutation that overlaps a set of five recurrent somatic tumor mutations for which recent functional studies demonstrated a gain-of-function effect due to constitutive activation of G protein downstream signaling cascades for some of the affected residues.

ATP synthase, H+ transporting, mitochondrial F1 complex, δ subunit (ATP5F1D; formerly ATP5D) is a subunit of mitochondrial ATP synthase and plays an important role in coupling proton translocation and ATP production. Here, we describe two individuals, each with homozygous missense variants in ATP5F1D, who presented with episodic lethargy, metabolic acidosis, 3-methylglutaconic aciduria, and hyperammonemia. Subject 1, homozygous for c.245C>T (p.Pro82Leu), presented with recurrent metabolic decompensation starting in the neonatal period, and subject 2, homozygous for c.317T>G (p.Val106Gly), presented with acute encephalopathy in childhood. Cultured skin fibroblasts from these individuals exhibited impaired assembly of F1FO ATP synthase and subsequent reduced complex V activity. Cells from subject 1 also exhibited a significant decrease in mitochondrial cristae. Knockdown of Drosophila ATPsynδ, the ATP5F1D homolog, in developing eyes and brains caused a near complete loss of the fly head, a phenotype that was fully rescued by wild-type human ATP5F1D. In contrast, expression of the ATP5F1D c.245C>T and c.317T>G variants rescued the head-size phenotype but recapitulated the eye and antennae defects seen in other genetic models of mitochondrial oxidative phosphorylation deficiency. Our data establish c.245C>T (p.Pro82Leu) and c.317T>G (p.Val106Gly) in ATP5F1D as pathogenic variants leading to a Mendelian mitochondrial disease featuring episodic metabolic decompensation.

Large-scale, population-based genomic studies have provided a context for modern medical genetics. Among such studies, however, African populations have remained relatively underrepresented. The breadth of genetic diversity across the African continent argues for an exploration of local genomic context to facilitate burgeoning disease mapping studies in Africa. We sought to characterize genetic variation and to assess population substructure within a cohort of HIV-positive children from Botswana-a Southern African country that is regionally underrepresented in genomic databases. Using whole-exome sequencing data from 164 Batswana and comparisons with 150 similarly sequenced HIV-positive Ugandan children, we found that 13%-25% of variation observed among Batswana was not captured by public databases. Uncaptured variants were significantly enriched (p = 2.2 × 10-16) for coding variants with minor allele frequencies between 1% and 5% and included predicted-damaging non-synonymous variants. Among variants found in public databases, corresponding allele frequencies varied widely, with Botswana having significantly higher allele frequencies among rare (<1%) pathogenic and damaging variants. Batswana clustered with other Southern African populations, but distinctly from 1000 Genomes African populations, and had limited evidence for admixture with extra-continental ancestries. We also observed a surprising lack of genetic substructure in Botswana, despite multiple tribal ethnicities and language groups, alongside a higher degree of relatedness than purported founder populations from the 1000 Genomes project. Our observations reveal a complex, but distinct, ancestral history and genomic architecture among Batswana and suggest that disease mapping within similar Southern African populations will require a deeper repository of genetic variation and allelic dependencies than presently exists.

Primary defects in lung branching morphogenesis, resulting in neonatal lethal pulmonary hypoplasias, are incompletely understood. To elucidate the pathogenetics of human lung development, we studied a unique collection of samples obtained from deceased individuals with clinically and histopathologically diagnosed interstitial neonatal lung disorders: acinar dysplasia (n = 14), congenital alveolar dysplasia (n = 2), and other lethal lung hypoplasias (n = 10). We identified rare heterozygous copy-number variant deletions or single-nucleotide variants (SNVs) involving TBX4 (n = 8 and n = 2, respectively) or FGF10 (n = 2 and n = 2, respectively) in 16/26 (61%) individuals. In addition to TBX4, the overlapping ∼2 Mb recurrent and nonrecurrent deletions at 17q23.1q23.2 identified in seven individuals with lung hypoplasia also remove a lung-specific enhancer region. Individuals with coding variants involving either TBX4 or FGF10 also harbored at least one non-coding SNV in the predicted lung-specific enhancer region, which was absent in 13 control individuals with the overlapping deletions but without any structural lung anomalies. The occurrence of rare coding variants involving TBX4 or FGF10 with the putative hypomorphic non-coding SNVs implies a complex compound inheritance of these pulmonary hypoplasias. Moreover, they support the importance of TBX4-FGF10-FGFR2 epithelial-mesenchymal signaling in human lung organogenesis and help to explain the histopathological continuum observed in these rare lethal developmental disorders of the lung.

EIF2AK1 and EIF2AK2 encode members of the eukaryotic translation initiation factor 2 alpha kinase (EIF2AK) family that inhibits protein synthesis in response to physiologic stress conditions. EIF2AK2 is also involved in innate immune response and the regulation of signal transduction, apoptosis, cell proliferation, and differentiation. Despite these findings, human disorders associated with deleterious variants in EIF2AK1 and EIF2AK2 have not been reported. Here, we describe the identification of nine unrelated individuals with heterozygous de novo missense variants in EIF2AK1 (1/9) or EIF2AK2 (8/9). Features seen in these nine individuals include white matter alterations (9/9), developmental delay (9/9), impaired language (9/9), cognitive impairment (8/9), ataxia (6/9), dysarthria in probands with verbal ability (6/9), hypotonia (7/9), hypertonia (6/9), and involuntary movements (3/9). Individuals with EIF2AK2 variants also exhibit neurological regression in the setting of febrile illness or infection. We use mammalian cell lines and proband-derived fibroblasts to further confirm the pathogenicity of variants in these genes and found reduced kinase activity. EIF2AKs phosphorylate eukaryotic translation initiation factor 2 subunit 1 (EIF2S1, also known as EIF2α), which then inhibits EIF2B activity. Deleterious variants in genes encoding EIF2B proteins cause childhood ataxia with central nervous system hypomyelination/vanishing white matter (CACH/VWM), a leukodystrophy characterized by neurologic regression in the setting of febrile illness and other stressors. Our findings indicate that EIF2AK2 missense variants cause a neurodevelopmental syndrome that may share phenotypic and pathogenic mechanisms with CACH/VWM.