Bi-allelic ACBD6 variants lead to a neurodevelopmental syndrome with progressive and complex movement disorders.

The acyl-CoA-binding domain-containing protein 6 (ACBD6) is ubiquitously expressed, plays a role in the acylation of lipids and proteins and regulates the N-myristoylation of proteins via N-myristoyltransferase enzymes (NMTs). However, its precise function in cells is still unclear, as is the consequence of ACBD6 defects on human pathophysiology. Using exome sequencing and extensive international data sharing efforts, we identified 45 affected individuals from 28 unrelated families (consanguinity 93%) with bi-allelic pathogenic, predominantly loss-of-function (18/20) variants in ACBD6. We generated zebrafish and Xenopus tropicalis acbd6 knockouts by CRISPR/Cas9 and characterized the role of ACBD6 on protein N-myristoylation with myristic acid alkyne (YnMyr) chemical proteomics in the model organisms and human cells, with the latter also being subjected further to ACBD6 peroxisomal localization studies. The affected individuals (23 males and 22 females), aged 1-50 years, typically present with a complex and progressive disease involving moderate-to-severe global developmental delay/intellectual disability (100%) with significant expressive language impairment (98%), movement disorders (97%), facial dysmorphism (95%) and mild cerebellar ataxia (85%) associated with gait impairment (94%), limb spasticity/hypertonia (76%), oculomotor (71%) and behavioural abnormalities (65%), overweight (59%), microcephaly (39%) and epilepsy (33%). The most conspicuous and common movement disorder was dystonia (94%), frequently leading to early-onset progressive postural deformities (97%), limb dystonia (55%) and cervical dystonia (31%). A jerky tremor in the upper limbs (63%), a mild head tremor (59%), parkinsonism/hypokinesia developing with advancing age (32%) and simple motor and vocal tics were among other frequent movement disorders. Midline brain malformations including corpus callosum abnormalities (70%), hypoplasia/agenesis of the anterior commissure (66%), short midbrain and small inferior cerebellar vermis (38% each) as well as hypertrophy of the clava (24%) were common neuroimaging findings. Acbd6-deficient zebrafish and Xenopus models effectively recapitulated many clinical phenotypes reported in patients including movement disorders, progressive neuromotor impairment, seizures, microcephaly, craniofacial dysmorphism and midbrain defects accompanied by developmental delay with increased mortality over time. Unlike ACBD5, ACBD6 did not show a peroxisomal localization and ACBD6-deficiency was not associated with altered peroxisomal parameters in patient fibroblasts. Significant differences in YnMyr-labelling were observed for 68 co- and 18 post-translationally N-myristoylated proteins in patient-derived fibroblasts. N-myristoylation was similarly affected in acbd6-deficient zebrafish and X. tropicalis models, including Fus, Marcks and Chchd-related proteins implicated in neurological diseases. The present study provides evidence that bi-allelic pathogenic variants in ACBD6 lead to a distinct neurodevelopmental syndrome accompanied by complex and progressive cognitive and movement disorders.

Evaluating the association of biallelic OGDHL variants with significant phenotypic heterogeneity.

BACKGROUND: Biallelic variants in OGDHL, encoding part of the α-ketoglutarate dehydrogenase complex, have been associated with highly heterogeneous neurological and neurodevelopmental disorders. However, the validity of this association remains to be confirmed. A second OGDHL patient cohort was recruited to carefully assess the gene-disease relationship. METHODS: Using an unbiased genotype-first approach, we screened large, multiethnic aggregated sequencing datasets worldwide for biallelic OGDHL variants. We used CRISPR/Cas9 to generate zebrafish knockouts of ogdhl, ogdh paralogs, and dhtkd1 to investigate functional relationships and impact during development. Functional complementation with patient variant transcripts was conducted to systematically assess protein functionality as a readout for pathogenicity. RESULTS: A cohort of 14 individuals from 12 unrelated families exhibited highly variable clinical phenotypes, with the majority of them presenting at least one additional variant, potentially accounting for a blended phenotype and complicating phenotypic understanding. We also uncovered extreme clinical heterogeneity and high allele frequencies, occasionally incompatible with a fully penetrant recessive disorder. Human cDNA of previously described and new variants were tested in an ogdhl zebrafish knockout model, adding functional evidence for variant reclassification. We disclosed evidence of hypomorphic alleles as well as a loss-of-function variant without deleterious effects in zebrafish variant testing also showing discordant familial segregation, challenging the relationship of OGDHL as a conventional Mendelian gene. Going further, we uncovered evidence for a complex compensatory relationship among OGDH, OGDHL, and DHTKD1 isoenzymes that are associated with neurodevelopmental disorders and exhibit complex transcriptional compensation patterns with partial functional redundancy. CONCLUSIONS: Based on the results of genetic, clinical, and functional studies, we formed three hypotheses in which to frame observations: biallelic OGDHL variants lead to a highly variable monogenic disorder, variants in OGDHL are following a complex pattern of inheritance, or they may not be causative at all. Our study further highlights the continuing challenges of assessing the validity of reported disease-gene associations and effects of variants identified in these genes. This is particularly more complicated in making genetic diagnoses based on identification of variants in genes presenting a highly heterogenous phenotype such as "OGDHL-related disorders".

Clinical, neuroradiological, and molecular characterization of mitochondrial threonyl-tRNA-synthetase (TARS2)-related disorder.

PURPOSE: Biallelic variants in TARS2, encoding the mitochondrial threonyl-tRNA-synthetase, have been reported in a small group of individuals displaying a neurodevelopmental phenotype but with limited neuroradiological data and insufficient evidence for causality of the variants. METHODS: Exome or genome sequencing was carried out in 15 families. Clinical and neuroradiological evaluation was performed for all affected individuals, including review of 10 previously reported individuals. The pathogenicity of TARS2 variants was evaluated using in vitro assays and a zebrafish model. RESULTS: We report 18 new individuals harboring biallelic TARS2 variants. Phenotypically, these individuals show developmental delay/intellectual disability, regression, cerebellar and cerebral atrophy, basal ganglia signal alterations, hypotonia, cerebellar signs, and increased blood lactate. In vitro studies showed that variants within the TARS2301-381 region had decreased binding to Rag GTPases, likely impairing mTORC1 activity. The zebrafish model recapitulated key features of the human phenotype and unraveled dysregulation of downstream targets of mTORC1 signaling. Functional testing of the variants confirmed the pathogenicity in a zebrafish model. CONCLUSION: We define the clinico-radiological spectrum of TARS2-related mitochondrial disease, unveil the likely involvement of the mTORC1 signaling pathway as a distinct molecular mechanism, and establish a TARS2 zebrafish model as an important tool to study variant pathogenicity.

Clinical, genetic, epidemiologic, evolutionary, and functional delineation of TSPEAR-related autosomal recessive ectodermal dysplasia 14.

TSPEAR variants cause autosomal recessive ectodermal dysplasia (ARED) 14. The function of TSPEAR is unknown. The clinical features, the mutation spectrum, and the underlying mechanisms of ARED14 are poorly understood. Combining data from new and previously published individuals established that ARED14 is primarily characterized by dental anomalies such as conical tooth cusps and hypodontia, like those seen in individuals with WNT10A-related odontoonychodermal dysplasia. AlphaFold-predicted structure-based analysis showed that most of the pathogenic TSPEAR missense variants likely destabilize the ß-propeller of the protein. Analysis of 100000 Genomes Project (100KGP) data revealed multiple founder TSPEAR variants across different populations. Mutational and recombination clock analyses demonstrated that non-Finnish European founder variants likely originated around the end of the last ice age, a period of major climatic transition. Analysis of gnomAD data showed that the non-Finnish European population TSPEAR gene-carrier rate is ∼1/140, making it one of the commonest AREDs. Phylogenetic and AlphaFold structural analyses showed that TSPEAR is an ortholog of drosophila Closca, an extracellular matrix-dependent signaling regulator. We, therefore, hypothesized that TSPEAR could have a role in enamel knot, a structure that coordinates patterning of developing tooth cusps. Analysis of mouse single-cell RNA sequencing (scRNA-seq) data revealed highly restricted expression of Tspear in clusters representing enamel knots. A tspeara -/-;tspearb -/- double-knockout zebrafish model recapitulated the clinical features of ARED14 and fin regeneration abnormalities of wnt10a knockout fish, thus suggesting interaction between tspear and wnt10a. In summary, we provide insights into the role of TSPEAR in ectodermal development and the evolutionary history, epidemiology, mechanisms, and consequences of its loss of function variants.

A homozygous MED11 C-terminal variant causes a lethal neurodegenerative disease.

PURPOSE: The mediator (MED) multisubunit-complex modulates the activity of the transcriptional machinery, and genetic defects in different MED subunits (17, 20, 27) have been implicated in neurologic diseases. In this study, we identified a recurrent homozygous variant in MED11 (c.325C>T; p.Arg109Ter) in 7 affected individuals from 5 unrelated families. METHODS: To investigate the genetic cause of the disease, exome or genome sequencing were performed in 5 unrelated families identified via different research networks and Matchmaker Exchange. Deep clinical and brain imaging evaluations were performed by clinical pediatric neurologists and neuroradiologists. The functional effect of the candidate variant on both MED11 RNA and protein was assessed using reverse transcriptase polymerase chain reaction and western blotting using fibroblast cell lines derived from 1 affected individual and controls and through computational approaches. Knockouts in zebrafish were generated using clustered regularly interspaced short palindromic repeats/Cas9. RESULTS: The disease was characterized by microcephaly, profound neurodevelopmental impairment, exaggerated startle response, myoclonic seizures, progressive widespread neurodegeneration, and premature death. Functional studies on patient-derived fibroblasts did not show a loss of protein function but rather disruption of the C-terminal of MED11, likely impairing binding to other MED subunits. A zebrafish knockout model recapitulates key clinical phenotypes. CONCLUSION: Loss of the C-terminal of MED subunit 11 may affect its binding efficiency to other MED subunits, thus implicating the MED-complex stability in brain development and neurodegeneration.

Biallelic variants in WARS1 cause a highly variable neurodevelopmental syndrome and implicate a critical exon for normal auditory function.

Aminoacyl-tRNA synthetases (ARSs) are essential enzymes for faithful assignment of amino acids to their cognate tRNA. Variants in ARS genes are frequently associated with clinically heterogeneous phenotypes in humans and follow both autosomal dominant or recessive inheritance patterns in many instances. Variants in tryptophanyl-tRNA synthetase 1 (WARS1) cause autosomal dominantly inherited distal hereditary motor neuropathy and Charcot-Marie-Tooth disease. Presently, only one family with biallelic WARS1 variants has been described. We present three affected individuals from two families with biallelic variants (p.Met1? and p.(Asp419Asn)) in WARS1, showing varying severities of developmental delay and intellectual disability. Hearing impairment and microcephaly, as well as abnormalities of the brain, skeletal system, movement/gait, and behavior were variable features. Phenotyping of knocked down wars-1 in a Caenorhabditis elegans model showed depletion is associated with defects in germ cell development. A wars1 knockout vertebrate model recapitulates the human clinical phenotypes, confirms variant pathogenicity, and uncovers evidence implicating the p.Met1? variant as potentially impacting an exon critical for normal hearing. Together, our findings provide consolidating evidence for biallelic disruption of WARS1 as causal for an autosomal recessive neurodevelopmental syndrome and present a vertebrate model that recapitulates key phenotypes observed in patients.

Biallelic variants in KARS1 are associated with neurodevelopmental disorders and hearing loss recapitulated by the knockout zebrafish.

PURPOSE: Pathogenic variants in Lysyl-tRNA synthetase 1 (KARS1) have increasingly been recognized as a cause of early-onset complex neurological phenotypes. To advance the timely diagnosis of KARS1-related disorders, we sought to delineate its phenotype and generate a disease model to understand its function in vivo. METHODS: Through international collaboration, we identified 22 affected individuals from 16 unrelated families harboring biallelic likely pathogenic or pathogenic in KARS1 variants. Sequencing approaches ranged from disease-specific panels to genome sequencing. We generated loss-of-function alleles in zebrafish. RESULTS: We identify ten new and four known biallelic missense variants in KARS1 presenting with a moderate-to-severe developmental delay, progressive neurological and neurosensory abnormalities, and variable white matter involvement. We describe novel KARS1-associated signs such as autism, hyperactive behavior, pontine hypoplasia, and cerebellar atrophy with prevalent vermian involvement. Loss of kars1 leads to upregulation of p53, tissue-specific apoptosis, and downregulation of neurodevelopmental related genes, recapitulating key tissue-specific disease phenotypes of patients. Inhibition of p53 rescued several defects of kars1-/- knockouts. CONCLUSION: Our work delineates the clinical spectrum associated with KARS1 defects and provides a novel animal model for KARS1-related human diseases revealing p53 signaling components as potential therapeutic targets.

A biallelic variant in CLRN2 causes non-syndromic hearing loss in humans.

Deafness, the most frequent sensory deficit in humans, is extremely heterogeneous with hundreds of genes involved. Clinical and genetic analyses of an extended consanguineous family with pre-lingual, moderate-to-profound autosomal recessive sensorineural hearing loss, allowed us to identify CLRN2, encoding a tetraspan protein, as a new deafness gene. Homozygosity mapping followed by exome sequencing identified a 14.96 Mb locus on chromosome 4p15.32p15.1 containing a likely pathogenic missense variant in CLRN2 (c.494C > A, NM_001079827.2) segregating with the disease. Using in vitro RNA splicing analysis, we show that the CLRN2 c.494C > A variant leads to two events: (1) the substitution of a highly conserved threonine (uncharged amino acid) to lysine (charged amino acid) at position 165, p.(Thr165Lys), and (2) aberrant splicing, with the retention of intron 2 resulting in a stop codon after 26 additional amino acids, p.(Gly146Lysfs*26). Expression studies and phenotyping of newly produced zebrafish and mouse models deficient for clarin 2 further confirm that clarin 2, expressed in the inner ear hair cells, is essential for normal organization and maintenance of the auditory hair bundles, and for hearing function. Together, our findings identify CLRN2 as a new deafness gene, which will impact future diagnosis and treatment for deaf patients.

MultiFRAGing: Rapid and Simultaneous Genotyping of Multiple Alleles in a Single Reaction.

Powerful and simple, RNA-guided CRISPR/Cas9 technology is a versatile genome editing tool that has revolutionized targeted mutagenesis. CRISPR-based genome editing has enabled large-scale functional genetic studies through the generation of gene knockouts in a variety of model organisms including zebrafish, and can be used to target multiple genes simultaneously. One of the challenges associated with the large scale application of this technique to zebrafish is the lack of a cost-effective method by which to identify mutants. To address this, we optimized the high-throughput, high-resolution fluorescent PCR-based fragment analysis method to develop MultiFRAGing - a robust and cost-effective method to genotype multiple targets in a single reaction. Our approach can identify indels in up to four targets from a single reaction, which represents a four-fold increase in genotyping throughput. This method can be used by any laboratory with access to capillary electrophoresis-based sequencing equipment.

Expanding the CRISPR Toolbox in Zebrafish for Studying Development and Disease.

The study of model organisms has revolutionized our understanding of the mechanisms underlying normal development, adult homeostasis, and human disease. Much of what we know about gene function in model organisms (and its application to humans) has come from gene knockouts: the ability to show analogous phenotypes upon gene inactivation in animal models. The zebrafish (Danio rerio) has become a popular model organism for many reasons, including the fact that it is amenable to various forms of genetic manipulation. The RNA-guided CRISPR/Cas9-mediated targeted mutagenesis approaches have provided powerful tools to manipulate the genome toward developing new disease models and understanding the pathophysiology of human diseases. CRISPR-based approaches are being used for the generation of both knockout and knock-in alleles, and also for applications including transcriptional modulation, epigenome editing, live imaging of the genome, and lineage tracing. Currently, substantial effort is being made to improve the specificity of Cas9, and to expand the target coverage of the Cas9 enzymes. Novel types of naturally occurring CRISPR systems [Cas12a (Cpf1); engineered variants of Cas9, such as xCas9 and SpCas9-NG], are being studied and applied to genome editing. Since the majority of pathogenic mutations are single point mutations, development of base editors to convert C:G to T:A or A:T to G:C has further strengthened the CRISPR toolbox. In this review, we provide an overview of the increasing number of novel CRISPR-based tools and approaches, including lineage tracing and base editing.