Microcephaly gene links trithorax and REST/NRSF to control neural stem cell proliferation and differentiation.

Microcephaly is a neurodevelopmental disorder causing significantly reduced cerebral cortex size. Many known microcephaly gene products localize to centrosomes, regulating cell fate and proliferation. Here, we identify and characterize a nuclear zinc finger protein, ZNF335/NIF-1, as a causative gene for severe microcephaly, small somatic size, and neonatal death. Znf335 null mice are embryonically lethal, and conditional knockout leads to severely reduced cortical size. RNA-interference and postmortem human studies show that ZNF335 is essential for neural progenitor self-renewal, neurogenesis, and neuronal differentiation. ZNF335 is a component of a vertebrate-specific, trithorax H3K4-methylation complex, directly regulating REST/NRSF, a master regulator of neural gene expression and cell fate, as well as other essential neural-specific genes. Our results reveal ZNF335 as an essential link between H3K4 complexes and REST/NRSF and provide the first direct genetic evidence that this pathway regulates human neurogenesis and neuronal differentiation.

A recurrent missense variant in the E3 ubiquitin ligase substrate recognition subunit FEM1B causes a rare syndromic neurodevelopmental disorder.

PURPOSE: Fem1 homolog B (FEM1B) acts as a substrate recognition subunit for ubiquitin ligase complexes belonging to the CULLIN 2-based E3 family. Several biological functions have been proposed for FEM1B, including a structurally resolved function as a sensor for redox cell status by controlling mitochondrial activity, but its implication in human disease remains elusive. METHODS: To understand the involvement of FEM1B in human disease, we made use of Matchmaker exchange platforms to identify individuals with de novo variants in FEM1B and performed their clinical evaluation. We performed functional validation using primary neuronal cultures and in utero electroporation assays, as well as experiments on patient's cells. RESULTS: Five individuals with a recurrent de novo missense variant in FEM1B were identified: NM_015322.5:c.377G>A NP_056137.1:p.(Arg126Gln) (FEM1BR126Q). Affected individuals shared a severe neurodevelopmental disorder with behavioral phenotypes and a variable set of malformations, including brain anomalies, clubfeet, skeletal abnormalities, and facial dysmorphism. Overexpression of the FEM1BR126Q variant but not FEM1B wild-type protein, during mouse brain development, resulted in delayed neuronal migration of the target cells. In addition, the individuals' cells exhibited signs of oxidative stress and induction of type I interferon signaling. CONCLUSION: Overall, our data indicate that p.(Arg126Gln) induces aberrant FEM1B activation, resulting in a gain-of-function mechanism associated with a severe syndromic developmental disorder in humans.

Expanding the genetics and phenotypes of ocular congenital cranial dysinnervation disorders.

PURPOSE: To identify genetic etiologies and genotype/phenotype associations for unsolved ocular congenital cranial dysinnervation disorders (oCCDDs). METHODS: We coupled phenotyping with exome or genome sequencing of 467 probands (550 affected and 1108 total individuals) with genetically unsolved oCCDDs, integrating analyses of pedigrees, human and animal model phenotypes, and de novo variants to identify rare candidate single nucleotide variants, insertion/deletions, and structural variants disrupting protein-coding regions. Prioritized variants were classified for pathogenicity and evaluated for genotype/phenotype correlations. RESULTS: Analyses elucidated phenotypic subgroups, identified pathogenic/likely pathogenic variant(s) in 43/467 probands (9.2%), and prioritized variants of uncertain significance in 70/467 additional probands (15.0%). These included known and novel variants in established oCCDD genes, genes associated with syndromes that sometimes include oCCDDs (e.g., MYH10, KIF21B, TGFBR2, TUBB6), genes that fit the syndromic component of the phenotype but had no prior oCCDD association (e.g., CDK13, TGFB2), genes with no reported association with oCCDDs or the syndromic phenotypes (e.g., TUBA4A, KIF5C, CTNNA1, KLB, FGF21), and genes associated with oCCDD phenocopies that had resulted in misdiagnoses. CONCLUSION: This study suggests that unsolved oCCDDs are clinically and genetically heterogeneous disorders often overlapping other Mendelian conditions and nominates many candidates for future replication and functional studies.

Recessive variants in COL25A1 gene as novel cause of arthrogryposis multiplex congenita with ocular congenital cranial dysinnervation disorder.

A proper interaction between muscle-derived collagen XXV and its motor neuron-derived receptors protein tyrosine phosphatases σ and δ (PTP σ/δ) is indispensable for intramuscular motor innervation. Despite this, thus far, pathogenic recessive variants in the COL25A1 gene had only been detected in a few patients with isolated ocular congenital cranial dysinnervation disorders. Here we describe five patients from three unrelated families with recessive missense and splice site COL25A1 variants presenting with a recognizable phenotype characterized by arthrogryposis multiplex congenita with or without an ocular congenital cranial dysinnervation disorder phenotype. The clinical features of the older patients remained stable over time, without central nervous system involvement. This study extends the phenotypic and genotypic spectrum of COL25A1 related conditions, and further adds to our knowledge of the complex process of intramuscular motor innervation. Our observations indicate a role for collagen XXV in regulating the appropriate innervation not only of extraocular muscles, but also of bulbar, axial, and limb muscles in the human.

Decreased ACKR3 (CXCR7) function causes oculomotor synkinesis in mice and humans.

Oculomotor synkinesis is the involuntary movement of the eyes or eyelids with a voluntary attempt at a different movement. The chemokine receptor CXCR4 and its ligand CXCL12 regulate oculomotor nerve development; mice with loss of either molecule have oculomotor synkinesis. In a consanguineous family with congenital ptosis and elevation of the ptotic eyelid with ipsilateral abduction, we identified a co-segregating homozygous missense variant (c.772G>A) in ACKR3, which encodes an atypical chemokine receptor that binds CXCL12 and functions as a scavenger receptor, regulating levels of CXCL12 available for CXCR4 signaling. The mutant protein (p.V258M) is expressed and traffics to the cell surface but has a lower binding affinity for CXCL12. Mice with loss of Ackr3 have variable phenotypes that include misrouting of the oculomotor and abducens nerves. All embryos show oculomotor nerve misrouting, ranging from complete misprojection in the midbrain, to aberrant peripheral branching, to a thin nerve, which aberrantly innervates the lateral rectus (as seen in Duane syndrome). The abducens nerve phenotype ranges from complete absence, to aberrant projections within the orbit, to a normal trajectory. Loss of ACKR3 in the midbrain leads to downregulation of CXCR4 protein, consistent with reports that excess CXCL12 causes ligand-induced degradation of CXCR4. Correspondingly, excess CXCL12 applied to ex vivo oculomotor slices causes axon misrouting, similar to inhibition of CXCR4. Thus, ACKR3, through its regulation of CXCL12 levels, is an important regulator of axon guidance in the oculomotor system; complete loss causes oculomotor synkinesis in mice, while reduced function causes oculomotor synkinesis in humans.

TUBB3 Arg262His causes a recognizable syndrome including CFEOM3, facial palsy, joint contractures, and early-onset peripheral neuropathy.

Microtubules are formed from heterodimers of alpha- and beta-tubulin, each of which has multiple isoforms encoded by separate genes. Pathogenic missense variants in multiple different tubulin isoforms cause brain malformations. Missense mutations in TUBB3, which encodes the neuron-specific beta-tubulin isotype, can cause congenital fibrosis of the extraocular muscles type 3 (CFEOM3) and/or malformations of cortical development, with distinct genotype-phenotype correlations. Here, we report fourteen individuals from thirteen unrelated families, each of whom harbors the identical NM_006086.4 (TUBB3):c.785G>A (p.Arg262His) variant resulting in a phenotype we refer to as the TUBB3 R262H syndrome. The affected individuals present at birth with ptosis, ophthalmoplegia, exotropia, facial weakness, facial dysmorphisms, and, in most cases, distal congenital joint contractures, and subsequently develop intellectual disabilities, gait disorders with proximal joint contractures, Kallmann syndrome (hypogonadotropic hypogonadism and anosmia), and a progressive peripheral neuropathy during the first decade of life. Subsets may also have vocal cord paralysis, auditory dysfunction, cyclic vomiting, and/or tachycardia at rest. All fourteen subjects share a recognizable set of brain malformations, including hypoplasia of the corpus callosum and anterior commissure, basal ganglia malformations, absent olfactory bulbs and sulci, and subtle cerebellar malformations. While similar, individuals with the TUBB3 R262H syndrome can be distinguished from individuals with the TUBB3 E410K syndrome by the presence of congenital and acquired joint contractures, an earlier onset peripheral neuropathy, impaired gait, and basal ganglia malformations.

MACF1 Mutations Encoding Highly Conserved Zinc-Binding Residues of the GAR Domain Cause Defects in Neuronal Migration and Axon Guidance.

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.

Phenotype delineation of ZNF462 related syndrome.

Zinc finger protein 462 (ZNF462) is a relatively newly discovered vertebrate specific protein with known critical roles in embryonic development in animal models. Two case reports and a case series study have described the phenotype of 10 individuals with ZNF462 loss of function variants. Herein, we present 14 new individuals with loss of function variants to the previous studies to delineate the syndrome of loss of function in ZNF462. Collectively, these 24 individuals present with recurring phenotypes that define a multiple congenital anomaly syndrome. Most have some form of developmental delay (79%) and a minority has autism spectrum disorder (33%). Characteristic facial features include ptosis (83%), down slanting palpebral fissures (58%), exaggerated Cupid's bow/wide philtrum (54%), and arched eyebrows (50%). Metopic ridging or craniosynostosis was found in a third of study participants and feeding problems in half. Other phenotype characteristics include dysgenesis of the corpus callosum in 25% of individuals, hypotonia in half, and structural heart defects in 21%. Using facial analysis technology, a computer algorithm applying deep learning was able to accurately differentiate individuals with ZNF462 loss of function variants from individuals with Noonan syndrome and healthy controls. In summary, we describe a multiple congenital anomaly syndrome associated with haploinsufficiency of ZNF462 that has distinct clinical characteristics and facial features.

PSMD12 haploinsufficiency in a neurodevelopmental disorder with autistic features.

Protein homeostasis is tightly regulated by the ubiquitin proteasome pathway. Disruption of this pathway gives rise to a host of neurological disorders. Through whole exome sequencing (WES) in families with neurodevelopmental disorders, we identified mutations in PSMD12, a core component of the proteasome, underlying a neurodevelopmental disorder with intellectual disability (ID) and features of autism spectrum disorder (ASD). We performed WES on six affected siblings from a multiplex family with ID and autistic features, the affected father, and two unaffected mothers, and a trio from a simplex family with one affected child with ID and periventricular nodular heterotopia. We identified an inherited heterozygous nonsense mutation in PSMD12 (NM_002816: c.367C>T: p.R123X) in the multiplex family and a de novo nonsense mutation in the same gene (NM_002816: c.601C>T: p.R201X) in the simplex family. PSMD12 encodes a non-ATPase regulatory subunit of the 26S proteasome. We confirm the association of PSMD12 with ID, present the first cases of inherited PSMD12 mutation, and demonstrate the heterogeneity of phenotypes associated with PSMD12 mutations.

Deficient activity of alanyl-tRNA synthetase underlies an autosomal recessive syndrome of progressive microcephaly, hypomyelination, and epileptic encephalopathy.

Aminoacyl-transfer RNA (tRNA) synthetases ligate amino acids to specific tRNAs and are essential for protein synthesis. Although alanyl-tRNA synthetase (AARS) is a synthetase implicated in a wide range of neurological disorders from Charcot-Marie-Tooth disease to infantile epileptic encephalopathy, there have been limited data on their pathogenesis. Here, we report loss-of-function mutations in AARS in two siblings with progressive microcephaly with hypomyelination, intractable epilepsy, and spasticity. Whole-exome sequencing identified that the affected individuals were compound heterozygous for mutations in AARS gene, c.2067dupC (p.Tyr690Leufs*3) and c.2738G>A (p.Gly913Asp). A lymphoblastoid cell line developed from one of the affected individuals showed a strong reduction in AARS abundance. The mutations decrease aminoacylation efficiency by 70%-90%. The p.Tyr690Leufs*3 mutation also abolished editing activity required for hydrolyzing misacylated tRNAs, thereby increasing errors during aminoacylation. Our study has extended potential mechanisms underlying AARS-related disorders to include destabilization of the protein, aminoacylation dysfunction, and defective editing activity.

Mutations in QARS, encoding glutaminyl-tRNA synthetase, cause progressive microcephaly, cerebral-cerebellar atrophy, and intractable seizures.

Progressive microcephaly is a heterogeneous condition with causes including mutations in genes encoding regulators of neuronal survival. Here, we report the identification of mutations in QARS (encoding glutaminyl-tRNA synthetase [QARS]) as the causative variants in two unrelated families affected by progressive microcephaly, severe seizures in infancy, atrophy of the cerebral cortex and cerebellar vermis, and mild atrophy of the cerebellar hemispheres. Whole-exome sequencing of individuals from each family independently identified compound-heterozygous mutations in QARS as the only candidate causative variants. QARS was highly expressed in the developing fetal human cerebral cortex in many cell types. The four QARS mutations altered highly conserved amino acids, and the aminoacylation activity of QARS was significantly impaired in mutant cell lines. Variants p.Gly45Val and p.Tyr57His were located in the N-terminal domain required for QARS interaction with proteins in the multisynthetase complex and potentially with glutamine tRNA, and recombinant QARS proteins bearing either substitution showed an over 10-fold reduction in aminoacylation activity. Conversely, variants p.Arg403Trp and p.Arg515Trp, each occurring in a different family, were located in the catalytic core and completely disrupted QARS aminoacylation activity in vitro. Furthermore, p.Arg403Trp and p.Arg515Trp rendered QARS less soluble, and p.Arg403Trp disrupted QARS-RARS (arginyl-tRNA synthetase 1) interaction. In zebrafish, homozygous qars loss of function caused decreased brain and eye size and extensive cell death in the brain. Our results highlight the importance of QARS during brain development and that epilepsy due to impairment of QARS activity is unusually severe in comparison to other aminoacyl-tRNA synthetase disorders.

Somatic mutations in cerebral cortical malformations.

BACKGROUND: Although there is increasing recognition of the role of somatic mutations in genetic disorders, the prevalence of somatic mutations in neurodevelopmental disease and the optimal techniques to detect somatic mosaicism have not been systematically evaluated. METHODS: Using a customized panel of known and candidate genes associated with brain malformations, we applied targeted high-coverage sequencing (depth, ≥200×) to leukocyte-derived DNA samples from 158 persons with brain malformations, including the double-cortex syndrome (subcortical band heterotopia, 30 persons), polymicrogyria with megalencephaly (20), periventricular nodular heterotopia (61), and pachygyria (47). We validated candidate mutations with the use of Sanger sequencing and, for variants present at unequal read depths, subcloning followed by colony sequencing. RESULTS: Validated, causal mutations were found in 27 persons (17%; range, 10 to 30% for each phenotype). Mutations were somatic in 8 of the 27 (30%), predominantly in persons with the double-cortex syndrome (in whom we found mutations in DCX and LIS1), persons with periventricular nodular heterotopia (FLNA), and persons with pachygyria (TUBB2B). Of the somatic mutations we detected, 5 (63%) were undetectable with the use of traditional Sanger sequencing but were validated through subcloning and subsequent sequencing of the subcloned DNA. We found potentially causal mutations in the candidate genes DYNC1H1, KIF5C, and other kinesin genes in persons with pachygyria. CONCLUSIONS: Targeted sequencing was found to be useful for detecting somatic mutations in patients with brain malformations. High-coverage sequencing panels provide an important complement to whole-exome and whole-genome sequencing in the evaluation of somatic mutations in neuropsychiatric disease. (Funded by the National Institute of Neurological Disorders and Stroke and others.).

Mammalian target of rapamycin pathway mutations cause hemimegalencephaly and focal cortical dysplasia.

Focal malformations of cortical development, including focal cortical dysplasia (FCD) and hemimegalencephaly (HME), are important causes of intractable childhood epilepsy. Using targeted and exome sequencing on DNA from resected brain samples and nonbrain samples from 53 patients with FCD or HME, we identified pathogenic germline and mosaic mutations in multiple PI3K/AKT pathway genes in 9 patients, and a likely pathogenic variant in 1 additional patient. Our data confirm the association of DEPDC5 with sporadic FCD but also implicate this gene for the first time in HME. Our findings suggest that modulation of the mammalian target of rapamycin pathway may hold promise for malformation-associated epilepsy.

Novel loss-of-function variants in DIAPH1 associated with syndromic microcephaly, blindness, and early onset seizures.

Exome sequencing identified homozygous loss-of-function variants in DIAPH1 (c.2769delT; p.F923fs and c.3145C>T; p.R1049X) in four affected individuals from two unrelated consanguineous families. The affected individuals in our report were diagnosed with postnatal microcephaly, early-onset epilepsy, severe vision impairment, and pulmonary symptoms including bronchiectasis and recurrent respiratory infections. A heterozygous DIAPH1 mutation was originally reported in one family with autosomal dominant deafness. Recently, however, a homozygous nonsense DIAPH1 mutation (c.2332C4T; p.Q778X) was reported in five siblings in a single family affected by microcephaly, blindness, early onset seizures, developmental delay, and bronchiectasis. The role of DIAPH1 was supported using parametric linkage analysis, RNA and protein studies in their patients' cell lines and further studies in human neural progenitors cells and a diap1 knockout mouse. In this report, the proband was initially brought to medical attention for profound metopic synostosis. Additional concerns arose when his head circumference did not increase after surgical release at 5 months of age and he was diagnosed with microcephaly and epilepsy at 6 months of age. Clinical exome analysis identified a homozygous DIAPH1 mutation. Another homozygous DIAPH1 mutation was identified in the research exome analysis of a second family with three siblings presenting with a similar phenotype. Importantly, no hearing impairment is reported in the homozygous affected individuals or in the heterozygous carrier parents in any of the families demonstrating the autosomal recessive microcephaly phenotype. These additional families provide further evidence of the likely causal relationship between DIAPH1 mutations and a neurodevelopmental disorder.

Exome sequencing and functional validation in zebrafish identify GTDC2 mutations as a cause of Walker-Warburg syndrome.

Whole-exome sequencing (WES), which analyzes the coding sequence of most annotated genes in the human genome, is an ideal approach to studying fully penetrant autosomal-recessive diseases, and it has been very powerful in identifying disease-causing mutations even when enrollment of affected individuals is limited by reduced survival. In this study, we combined WES with homozygosity analysis of consanguineous pedigrees, which are informative even when a single affected individual is available, to identify genetic mutations responsible for Walker-Warburg syndrome (WWS), a genetically heterogeneous autosomal-recessive disorder that severely affects the development of the brain, eyes, and muscle. Mutations in seven genes are known to cause WWS and explain 50%-60% of cases, but multiple additional genes are expected to be mutated because unexplained cases show suggestive linkage to diverse loci. Using WES in consanguineous WWS-affected families, we found multiple deleterious mutations in GTDC2 (also known as AGO61). GTDC2's predicted role as an uncharacterized glycosyltransferase is consistent with the function of other genes that are known to be mutated in WWS and that are involved in the glycosylation of the transmembrane receptor dystroglycan. Therefore, to explore the role of GTDC2 loss of function during development, we used morpholino-mediated knockdown of its zebrafish ortholog, gtdc2. We found that gtdc2 knockdown in zebrafish replicates all WWS features (hydrocephalus, ocular defects, and muscular dystrophy), strongly suggesting that GTDC2 mutations cause WWS.

Familial recurrences of FOXG1-related disorder: Evidence for mosaicism.

FOXG1-related disorders are caused by heterozygous mutations in FOXG1 and result in a spectrum of neurodevelopmental phenotypes including postnatal microcephaly, intellectual disability with absent speech, epilepsy, chorea, and corpus callosum abnormalities. The recurrence risk for de novo mutations in FOXG1-related disorders is assumed to be low. Here, we describe three unrelated sets of full siblings with mutations in FOXG1 (c.515_577del63, c.460dupG, and c.572T > G), representing familial recurrence of the disorder. In one family, we have documented maternal somatic mosaicism for the FOXG1 mutation, and all of the families presumably represent parental gonadal (or germline) mosaicism. To our knowledge, mosaicism has not been previously reported in FOXG1-related disorders. Therefore, this report provides evidence that germline mosaicism for FOXG1 mutations is a likely explanation for familial recurrence and should be considered during recurrence risk counseling for families of children with FOXG1-related disorders.

Human mutations in NDE1 cause extreme microcephaly with lissencephaly [corrected].

Genes disrupted in human microcephaly (meaning "small brain") define key regulators of neural progenitor proliferation and cell-fate specification. In comparison, genes mutated in human lissencephaly (lissos means smooth and cephalos means brain) highlight critical regulators of neuronal migration. Here, we report two families with extreme microcephaly and grossly simplified cortical gyral structure, a condition referred to as microlissencephaly, and show that they carry homozygous frameshift mutations in NDE1, which encodes a multidomain protein that localizes to the centrosome and mitotic spindle poles. Both human mutations in NDE1 truncate the C-terminal NDE1domains, which are essential for interactions with cytoplasmic dynein and thus for regulation of cytoskeletal dynamics in mitosis and for cell-cycle-dependent phosphorylation of NDE1 by Cdk1. We show that the patient NDE1 proteins are unstable, cannot bind cytoplasmic dynein, and do not localize properly to the centrosome. Additionally, we show that CDK1 phosphorylation at T246, which is within the C-terminal region disrupted by the mutations, is required for cell-cycle progression from the G2 to the M phase. The role of NDE1 in cell-cycle progression probably contributes to the profound neuronal proliferation defects evident in Nde1-null mice and patients with NDE1 mutations, demonstrating the essential role of NDE1 in human cerebral cortical neurogenesis.

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory disorders.

Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN) cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generate single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. We evaluate enhancer activity for 59 elements using an in vivo transgenic assay and validate 44 (75%), demonstrating that single cell accessibility can be a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieve significant reduction in our variant search space and nominate candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 - as well as candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work delivers non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

Gene identification for ocular congenital cranial motor neuron disorders using human sequencing, zebrafish screening, and protein binding microarrays.

Purpose: To functionally evaluate novel human sequence-derived candidate genes and variants for unsolved ocular congenital cranial dysinnervation disorders (oCCDDs). Methods: Through exome and genome sequencing of a genetically unsolved human oCCDD cohort, we previously identified variants in 80 strong candidate genes. Here, we further prioritized a subset of these (43 human genes, 57 zebrafish genes) using a G0 CRISPR/Cas9-based knockout assay in zebrafish and generated F2 germline mutants for seventeen. We tested the functionality of variants of uncertain significance in known and novel candidate transcription factor-encoding genes through protein binding microarrays. Results: We first demonstrated the feasibility of the G0 screen by targeting known oCCDD genes phox2a and mafba. 70-90% of gene-targeted G0 zebrafish embryos recapitulated germline homozygous null-equivalent phenotypes. Using this approach, we then identified three novel candidate oCCDD genes (SEMA3F, OLIG2, and FRMD4B) with putative contributions to human and zebrafish cranial motor development. In addition, protein binding microarrays demonstrated reduced or abolished DNA binding of human variants of uncertain significance in known and novel sequence-derived transcription factors PHOX2A (p.(Trp137Cys)), MAFB (p.(Glu223Lys)), and OLIG2 (p.(Arg156Leu)). Conclusions: This study nominates three strong novel candidate oCCDD genes (SEMA3F, OLIG2, and FRMD4B) and supports the functionality and putative pathogenicity of transcription factor candidate variants PHOX2A p.(Trp137Cys), MAFB p.(Glu223Lys), and OLIG2 p.(Arg156Leu). Our findings support that G0 loss-of-function screening in zebrafish can be coupled with human sequence analysis and protein binding microarrays to aid in prioritizing oCCDD candidate genes/variants.

Expanding the genetics and phenotypes of ocular congenital cranial dysinnervation disorders.

Purpose: To identify genetic etiologies and genotype/phenotype associations for unsolved ocular congenital cranial dysinnervation disorders (oCCDDs). Methods: We coupled phenotyping with exome or genome sequencing of 467 pedigrees with genetically unsolved oCCDDs, integrating analyses of pedigrees, human and animal model phenotypes, and de novo variants to identify rare candidate single nucleotide variants, insertion/deletions, and structural variants disrupting protein-coding regions. Prioritized variants were classified for pathogenicity and evaluated for genotype/phenotype correlations. Results: Analyses elucidated phenotypic subgroups, identified pathogenic/likely pathogenic variant(s) in 43/467 probands (9.2%), and prioritized variants of uncertain significance in 70/467 additional probands (15.0%). These included known and novel variants in established oCCDD genes, genes associated with syndromes that sometimes include oCCDDs (e.g., MYH10, KIF21B, TGFBR2, TUBB6), genes that fit the syndromic component of the phenotype but had no prior oCCDD association (e.g., CDK13, TGFB2), genes with no reported association with oCCDDs or the syndromic phenotypes (e.g., TUBA4A, KIF5C, CTNNA1, KLB, FGF21), and genes associated with oCCDD phenocopies that had resulted in misdiagnoses. Conclusion: This study suggests that unsolved oCCDDs are clinically and genetically heterogeneous disorders often overlapping other Mendelian conditions and nominates many candidates for future replication and functional studies.