Generation of induced pluripotent stem cell lines from two unrelated patients affected by intellectual disability carrying homozygous variants in SGIP1.

Intellectual disability (ID) is a diverse neurodevelopmental condition and almost half of the cases have a genetic etiology. SGIP1 acts as an endocytic protein that influences the signaling of receptors in neuronal systems related to energy homeostasis through its interaction with endophilins. This study focuses on the generation and characterization of induced pluripotent stem cells (iPSC) from two unrelated patients due to a frameshift variant (c.764dupA, NM_032291.4) and a splice donor site variant (c.74 + 1G > A, NM_032291.4) in the SGIP1 gene.

SCN1A-deficient excitatory neuronal networks display mutation-specific phenotypes.

Dravet syndrome is a severe epileptic encephalopathy, characterized by (febrile) seizures, behavioural problems and developmental delay. Eighty per cent of patients with Dravet syndrome have a mutation in SCN1A, encoding Nav1.1. Milder clinical phenotypes, such as GEFS+ (generalized epilepsy with febrile seizures plus), can also arise from SCN1A mutations. Predicting the clinical phenotypic outcome based on the type of mutation remains challenging, even when the same mutation is inherited within one family. This clinical and genetic heterogeneity adds to the difficulties of predicting disease progression and tailoring the prescription of anti-seizure medication. Understanding the neuropathology of different SCN1A mutations may help to predict the expected clinical phenotypes and inform the selection of best-fit treatments. Initially, the loss of Na+-current in inhibitory neurons was recognized specifically to result in disinhibition and consequently seizure generation. However, the extent to which excitatory neurons contribute to the pathophysiology is currently debated and might depend on the patient clinical phenotype or the specific SCN1A mutation. To examine the genotype-phenotype correlations of SCN1A mutations in relation to excitatory neurons, we investigated a panel of patient-derived excitatory neuronal networks differentiated on multi-electrode arrays. We included patients with different clinical phenotypes, harbouring various SCN1A mutations, along with a family in which the same mutation led to febrile seizures, GEFS+ or Dravet syndrome. We hitherto describe a previously unidentified functional excitatory neuronal network phenotype in the context of epilepsy, which corresponds to seizurogenic network prediction patterns elicited by proconvulsive compounds. We found that excitatory neuronal networks were affected differently, depending on the type of SCN1A mutation, but did not segregate according to clinical severity. Specifically, loss-of-function mutations could be distinguished from missense mutations, and mutations in the pore domain could be distinguished from mutations in the voltage sensing domain. Furthermore, all patients showed aggravated neuronal network responses at febrile temperatures compared with controls. Finally, retrospective drug screening revealed that anti-seizure medication affected GEFS+ patient- but not Dravet patient-derived neuronal networks in a patient-specific and clinically relevant manner. In conclusion, our results indicate a mutation-specific excitatory neuronal network phenotype, which recapitulates the foremost clinically relevant features, providing future opportunities for precision therapies.

Genetic diagnostic approach to intellectual disability and multiple congenital anomalies in Indonesia.

Intellectual disability (ID) and multiple congenital anomalies (MCA) are major contributors to infant mortality, childhood morbidity, and long-term disability, with multifactorial aetiology including genetics. We aim to set a diagnostic approach for genetic evaluation of patients with ID and MCA, which can be applied efficiently with a good diagnostic rate in Indonesia or other low resources settings. Out of 131 ID cases, twenty-three individuals with ID/global developmental delay (GDD) and MCA were selected from two-steps of dysmorphology screening and evaluation. Genetic analysis included chromosomal microarray (CMA) analysis, targeted panel gene sequencing, and exome sequencing (ES). CMA revealed conclusive results for seven individuals. Meanwhile, two out of four cases were diagnosed by targeted gene sequencing. Five out of seven individuals were diagnosed using ES testing. Based on the experience, a novel and comprehensive flowchart combining thorough physical and dysmorphology evaluation, followed by suitable genetic tests is proposed as a diagnostic approach to elucidate the genetic factor(s) of ID/GDD and MCA in low resources settings such as Indonesia.

A complex structural variant near SOX3 causes X-linked split-hand/foot malformation.

Split-hand/foot malformation (SHFM) is a congenital limb defect most typically presenting with median clefts in hands and/or feet, that can occur in a syndromic context as well as in isolated form. SHFM is caused by failure to maintain normal apical ectodermal ridge function during limb development. Although several genes and contiguous gene syndromes are implicated in the monogenic etiology of isolated SHFM, the disorder remains genetically unexplained for many families and associated genetic loci. We describe a family with isolated X-linked SHFM, for which the causative variant could be detected after a diagnostic journey of 20 years. We combined well-established approaches including microarray-based copy number variant analysis and fluorescence in situ hybridization coupled with optical genome mapping and whole genome sequencing. This strategy identified a complex structural variant (SV) comprising a 165-kb gain of 15q26.3 material ([GRCh37/hg19] chr15:99795320-99960362dup) inserted in inverted position at the site of a 38-kb deletion on Xq27.1 ([GRCh37/hg19] chrX:139481061-139518989del). In silico analysis suggested that the SV disrupts the regulatory framework on the X chromosome and may lead to SOX3 misexpression. We hypothesize that SOX3 dysregulation in the developing limb disturbed the fine balance between morphogens required for maintaining AER function, resulting in SHFM in this family.

Missense MED12 variants in 22 males with intellectual disability: From nonspecific symptoms to complete syndromes.

We describe the phenotype of 22 male patients (20 probands) carrying a hemizygous missense variant in MED12. The phenotypic spectrum is very broad ranging from nonspecific intellectual disability (ID) to the three well-known syndromes: Opitz-Kaveggia syndrome, Lujan-Fryns syndrome, or Ohdo syndrome. The identified variants were randomly distributed throughout the gene (p = 0.993, χ2 test), but mostly outside the functional domains (p = 0.004; χ2 test). Statistical analyses did not show a correlation between the MED12-related phenotypes and the locations of the variants (p = 0.295; Pearson correlation), nor the protein domain involved (p = 0.422; Pearson correlation). In conclusion, establishing a genotype-phenotype correlation in MED12-related diseases remains challenging. Therefore, we think that patients with a causative MED12 variant are currently underdiagnosed due to the broad patients' clinical presentations.

Large-Scale Functional Assessment of Genes Involved in Rare Diseases with Intellectual Disabilities Unravels Unique Developmental and Behaviour Profiles in Mouse Models.

Major progress has been made over the last decade in identifying novel genes involved in neurodevelopmental disorders, although the task of elucidating their corresponding molecular and pathophysiological mechanisms, which are an essential prerequisite for developing therapies, has fallen far behind. We selected 45 genes for intellectual disabilities to generate and characterize mouse models. Thirty-nine of them were based on the frequency of pathogenic variants in patients and literature reports, with several corresponding to de novo variants, and six other candidate genes. We used an extensive screen covering the development and adult stages, focusing specifically on behaviour and cognition to assess a wide range of functions and their pathologies, ranging from basic neurological reflexes to cognitive abilities. A heatmap of behaviour phenotypes was established, together with the results of selected mutants. Overall, three main classes of mutant lines were identified based on activity phenotypes, with which other motor or cognitive deficits were associated. These data showed the heterogeneity of phenotypes between mutation types, recapitulating several human features, and emphasizing the importance of such systematic approaches for both deciphering genetic etiological causes of ID and autism spectrum disorders, and for building appropriate therapeutic strategies.

Biallelic loss of EMC10 leads to mild to severe intellectual disability.

The endoplasmic reticulum membrane protein complex subunit 10 (EMC10) is a highly conserved protein responsible for the post-translational insertion of tail-anchored membrane proteins into the endoplasmic reticulum in a defined topology. Two biallelic variants in EMC10 have previously been associated with a neurodevelopmental disorder. Utilizing exome sequencing and international data sharing we have identified 10 affected individuals from six independent families with five new biallelic loss-of-function and one previously reported recurrent EMC10 variants. This report expands the molecular and clinical spectrum of EMC10 deficiency, provides a comprehensive dysmorphological assessment and highlights an overlap between the clinical features of EMC10-and EMC1-related disease.

Loss-of-function variants in the schizophrenia risk gene SETD1A alter neuronal network activity in human neurons through the cAMP/PKA pathway.

Heterozygous loss-of-function (LoF) mutations in SETD1A, which encodes a subunit of histone H3 lysine 4 methyltransferase, cause a neurodevelopmental syndrome and increase the risk for schizophrenia. Using CRISPR-Cas9, we generate excitatory/inhibitory neuronal networks from human induced pluripotent stem cells with a SETD1A heterozygous LoF mutation (SETD1A+/-). Our data show that SETD1A haploinsufficiency results in morphologically increased dendritic complexity and functionally increased bursting activity. This network phenotype is primarily driven by SETD1A haploinsufficiency in glutamatergic neurons. In accordance with the functional changes, transcriptomic profiling reveals perturbations in gene sets associated with glutamatergic synaptic function. At the molecular level, we identify specific changes in the cyclic AMP (cAMP)/Protein Kinase A pathway pointing toward a hyperactive cAMP pathway in SETD1A+/- neurons. Finally, by pharmacologically targeting the cAMP pathway, we are able to rescue the network deficits in SETD1A+/- cultures. Our results demonstrate a link between SETD1A and the cAMP-dependent pathway in human neurons.

Generation of induced pluripotent stem cell lines carrying monoallelic (UCSFi001-A-60) or biallelic (UCSFi001-A-61; UCSFi001-A-62) frameshift variants in CACNA1A using CRISPR/Cas9.

CACNA1A encodes a P/Q-type voltage-gated calcium channel. Heterozygous loss-of-function variants in this gene have been associated with episodic ataxia type 2. In this study, we used CRISPR/Cas9 to generate isogenic human induced pluripotent stem cell lines with a gene-dosage dependent deficiency of CACNA1A. We obtained one clone with monoallelic (UCSFi001-A-60) and two clones with biallelic (UCSFi001-A-61; UCSFi001-A-62) frameshift variants in CACNA1A. All three lines showed expression of pluripotency markers and a normal karyotype.

Phenotypic and mutational spectrum of ROR2-related Robinow syndrome.

Robinow syndrome is characterized by a triad of craniofacial dysmorphisms, disproportionate-limb short stature, and genital hypoplasia. A significant degree of phenotypic variability seems to correlate with different genes/loci. Disturbances of the noncanonical WNT-pathway have been identified as the main cause of the syndrome. Biallelic variants in ROR2 cause an autosomal recessive form of the syndrome with distinctive skeletal findings. Twenty-two patients with a clinical diagnosis of autosomal recessive Robinow syndrome were screened for variants in ROR2 using multiple molecular approaches. We identified 25 putatively pathogenic ROR2 variants, 16 novel, including single nucleotide variants and exonic deletions. Detailed phenotypic analyses revealed that all subjects presented with a prominent forehead, hypertelorism, short nose, abnormality of the nasal tip, brachydactyly, mesomelic limb shortening, short stature, and genital hypoplasia in male patients. A total of 19 clinical features were present in more than 75% of the subjects, thus pointing to an overall uniformity of the phenotype. Disease-causing variants in ROR2, contribute to a clinically recognizable autosomal recessive trait phenotype with multiple skeletal defects. A comprehensive quantitative clinical evaluation of this cohort delineated the phenotypic spectrum of ROR2-related Robinow syndrome. The identification of exonic deletion variant alleles further supports the contention of a loss-of-function mechanism in the etiology of the syndrome.

THUMPD1 bi-allelic variants cause loss of tRNA acetylation and a syndromic neurodevelopmental disorder.

Covalent tRNA modifications play multi-faceted roles in tRNA stability, folding, and recognition, as well as the rate and fidelity of translation, and other cellular processes such as growth, development, and stress responses. Mutations in genes that are known to regulate tRNA modifications lead to a wide array of phenotypes and diseases including numerous cognitive and neurodevelopmental disorders, highlighting the critical role of tRNA modification in human disease. One such gene, THUMPD1, is involved in regulating tRNA N4-acetylcytidine modification (ac4C), and recently was proposed as a candidate gene for autosomal-recessive intellectual disability. Here, we present 13 individuals from 8 families who harbor rare loss-of-function variants in THUMPD1. Common phenotypic findings included global developmental delay, speech delay, moderate to severe intellectual deficiency, behavioral abnormalities such as angry outbursts, facial dysmorphism, and ophthalmological abnormalities. We demonstrate that the bi-allelic variants identified cause loss of function of THUMPD1 and that this defect results in a loss of ac4C modification in small RNAs, and of individually purified tRNA-Ser-CGA. We further corroborate this effect by showing a loss of tRNA acetylation in two CRISPR-Cas9-generated THUMPD1 KO cell lines. In addition, we also show the resultant amino acid substitution that occurs in a missense THUMPD1 allele identified in an individual with compound heterozygous variants results in a marked decrease in THUMPD1 stability and RNA-binding capacity. Taken together, these results suggest that the lack of tRNA acetylation due to THUMPD1 loss of function results in a syndromic form of intellectual disability associated with developmental delay, behavioral abnormalities, hearing loss, and facial dysmorphism.

Cadherin-13 is a critical regulator of GABAergic modulation in human stem-cell-derived neuronal networks.

Activity in the healthy brain relies on a concerted interplay of excitation (E) and inhibition (I) via balanced synaptic communication between glutamatergic and GABAergic neurons. A growing number of studies imply that disruption of this E/I balance is a commonality in many brain disorders; however, obtaining mechanistic insight into these disruptions, with translational value for the patient, has typically been hampered by methodological limitations. Cadherin-13 (CDH13) has been associated with autism and attention-deficit/hyperactivity disorder. CDH13 localizes at inhibitory presynapses, specifically of parvalbumin (PV) and somatostatin (SST) expressing GABAergic neurons. However, the mechanism by which CDH13 regulates the function of inhibitory synapses in human neurons remains unknown. Starting from human-induced pluripotent stem cells, we established a robust method to generate a homogenous population of SST and MEF2C (PV-precursor marker protein) expressing GABAergic neurons (iGABA) in vitro, and co-cultured these with glutamatergic neurons at defined E/I ratios on micro-electrode arrays. We identified functional network parameters that are most reliably affected by GABAergic modulation as such, and through alterations of E/I balance by reduced expression of CDH13 in iGABAs. We found that CDH13 deficiency in iGABAs decreased E/I balance by means of increased inhibition. Moreover, CDH13 interacts with Integrin-ß1 and Integrin-ß3, which play opposite roles in the regulation of inhibitory synaptic strength via this interaction. Taken together, this model allows for standardized investigation of the E/I balance in a human neuronal background and can be deployed to dissect the cell-type-specific contribution of disease genes to the E/I balance.

Imbalanced autophagy causes synaptic deficits in a human model for neurodevelopmental disorders.

Macroautophagy (hereafter referred to as autophagy) is a finely tuned process of programmed degradation and recycling of proteins and cellular components, which is crucial in neuronal function and synaptic integrity. Mounting evidence implicates chromatin remodeling in fine-tuning autophagy pathways. However, this epigenetic regulation is poorly understood in neurons. Here, we investigate the role in autophagy of KANSL1, a member of the nonspecific lethal complex, which acetylates histone H4 on lysine 16 (H4K16ac) to facilitate transcriptional activation. Loss-of-function of KANSL1 is strongly associated with the neurodevelopmental disorder Koolen-de Vries Syndrome (KdVS). Starting from KANSL1-deficient human induced-pluripotent stem cells, both from KdVS patients and genome-edited lines, we identified SOD1 (superoxide dismutase 1), an antioxidant enzyme, to be significantly decreased, leading to a subsequent increase in oxidative stress and autophagosome accumulation. In KANSL1-deficient neurons, autophagosome accumulation at excitatory synapses resulted in reduced synaptic density, reduced GRIA/AMPA receptor-mediated transmission and impaired neuronal network activity. Furthermore, we found that increased oxidative stress-mediated autophagosome accumulation leads to increased MTOR activation and decreased lysosome function, further preventing the clearing of autophagosomes. Finally, by pharmacologically reducing oxidative stress, we could rescue the aberrant autophagosome formation as well as synaptic and neuronal network activity in KANSL1-deficient neurons. Our findings thus point toward an important relation between oxidative stress-induced autophagy and synapse function, and demonstrate the importance of H4K16ac-mediated changes in chromatin structure to balance reactive oxygen species- and MTOR-dependent autophagy.Abbreviations: APO: apocynin; ATG: autophagy related; BAF: bafilomycin A1; BSO: buthionine sulfoximine; CV: coefficient of variation; DIV: days in vitro; H4K16ac: histone 4 lysine 16 acetylation; iPSC: induced-pluripotent stem cell; KANSL1: KAT8 regulatory NSL complex subunit 1; KdVS: Koolen-de Vries Syndrome; LAMP1: lysosomal associated membrane protein 1; MAP1LC3/LC3: microtubule associated protein 1 light chain 3; MEA: micro-electrode array; MTOR: mechanistic target of rapamycin kinase; NSL complex: nonspecific lethal complex; 8-oxo-dG: 8-hydroxydesoxyguanosine; RAP: rapamycin; ROS: reactive oxygen species; sEPSCs: spontaneous excitatory postsynaptic currents; SOD1: superoxide dismutase 1; SQSTM1/p62: sequestosome 1; SYN: synapsin; WRT: wortmannin.

Human Induced Pluripotent Stem Cell-Based Modelling of Spinocerebellar Ataxias.

Dominant spinocerebellar ataxias (SCAs) constitute a large group of phenotypically and genetically heterogeneous disorders that mainly present with dysfunction of the cerebellum as their main hallmark. Although animal and cell models have been highly instrumental for our current insight into the underlying disease mechanisms of these neurodegenerative disorders, they do not offer the full human genetic and physiological context. The advent of human induced pluripotent stem cells (hiPSCs) and protocols to differentiate these into essentially every cell type allows us to closely model SCAs in a human context. In this review, we systematically summarize recent findings from studies using hiPSC-based modelling of SCAs, and discuss what knowledge has been gained from these studies. We conclude that hiPSC-based models are a powerful tool for modelling SCAs as they contributed to new mechanistic insights and have the potential to serve the development of genetic therapies. However, the use of standardized methods and multiple clones of isogenic lines are essential to increase validity and reproducibility of the insights gained.

Biallelic in-frame deletion of SOX4 is associated with developmental delay, hypotonia and intellectual disability.

Intellectual disability (ID) represents an extremely heterogeneous group of disorders, characterized by significant limitations in intellectual function and adaptive behavior. Among the monogenic causes, autosomal recessive genes (ARID) are responsible for more than 50% of ID. Here, we report a novel in-frame homozygous deletion variant [c.730_753del; p.(Ala244_Gly251del)] in SOX4 (sex-determining region Y-related high-mobility group box 4), segregating with moderate to severe ID, hypotonia, and developmental delay in a Pakistani family. Our identified variant p.(Ala244_Gly251del) is predicted to remove evolutionarily conserved residues from the interdomain region and may destabilize the protein secondary structure. SOX4 belongs to group C of the SOX transcription regulating family known to be involved in early embryo development. Single-cell RNA data analysis of developing telencephalon revealed highly overlapping expression of SOX4 with SOX11 and DCX, known neurogenesis regulators. Our study expands the mutational landscape of SOX4 and the repertoire of the known genetic causes of ARID.

The complexities of CACNA1A in clinical neurogenetics.

Variants in CACNA1A are classically related to episodic ataxia type 2, familial hemiplegic migraine type 1, and spinocerebellar ataxia type 6. Over the years, CACNA1A has been associated with a broader spectrum of phenotypes. Targeted analysis and unbiased sequencing of CACNA1A result not only in clear molecular diagnoses, but also in large numbers of variants of uncertain significance (VUS), or likely pathogenic variants with a phenotype that does not directly match the CACNA1A spectrum. Over the last years, targeted and clinical exome sequencing in our center has identified 41 CACNA1A variants. Ultimately, variants were considered pathogenic or likely pathogenic in 23 cases, with most phenotypes ranging from episodic or progressive ataxia to more complex ataxia syndromes, as well as intellectual disability and epilepsy. In two cases, the causality of the variant was discarded based on non-segregation or an alternative diagnosis. In the remaining 16 cases, the variant was classified as uncertain, due to lack of opportunities for segregation analysis or uncertain association with a non-classic phenotype. Phenotypic variability and the large number of VUS make CACNA1A a challenging gene for neurogenetic diagnostics. Accessible functional read-outs are clearly needed, especially in cases with a non-classic phenotype.

Analysis of a Set of KDM5C Regulatory Genes Mutated in Neurodevelopmental Disorders Identifies Temporal Coexpression Brain Signatures.

Dysregulation of transcriptional pathways is observed in multiple forms of neurodevelopmental disorders (NDDs), such as intellectual disability (ID), epilepsy and autism spectrum disorder (ASD). We previously demonstrated that the NDD genes encoding lysine-specific demethylase 5C (KDM5C) and its transcriptional regulators Aristaless related-homeobox (ARX), PHD Finger Protein 8 (PHF8) and Zinc Finger Protein 711 (ZNF711) are functionally connected. Here, we show their relation to each other with respect to the expression levels in human and mouse datasets and in vivo mouse analysis indicating that the coexpression of these syntenic X-chromosomal genes is temporally regulated in brain areas and cellular sub-types. In co-immunoprecipitation assays, we found that the homeotic transcription factor ARX interacts with the histone demethylase PHF8, indicating that this transcriptional axis is highly intersected. Furthermore, the functional impact of pathogenic mutations of ARX, KDM5C, PHF8 and ZNF711 was tested in lymphoblastoid cell lines (LCLs) derived from children with varying levels of syndromic ID establishing the direct correlation between defects in the KDM5C-H3K4me3 pathway and ID severity. These findings reveal novel insights into epigenetic processes underpinning NDD pathogenesis and provide new avenues for assessing developmental timing and critical windows for potential treatments.

Human neuronal networks on micro-electrode arrays are a highly robust tool to study disease-specific genotype-phenotype correlations in vitro.

Micro-electrode arrays (MEAs) are increasingly used to characterize neuronal network activity of human induced pluripotent stem cell (hiPSC)-derived neurons. Despite their gain in popularity, MEA recordings from hiPSC-derived neuronal networks are not always used to their full potential in respect to experimental design, execution, and data analysis. Therefore, we benchmarked the robustness of MEA-derived neuronal activity patterns from ten healthy individual control lines, and uncover comparable network phenotypes. To achieve standardization, we provide recommendations on experimental design and analysis. With such standardization, MEAs can be used as a reliable platform to distinguish (disease-specific) network phenotypes. In conclusion, we show that MEAs are a powerful and robust tool to uncover functional neuronal network phenotypes from hiPSC-derived neuronal networks, and provide an important resource to advance the hiPSC field toward the use of MEAs for disease phenotyping and drug discovery.

Bi-allelic variants in the ER quality-control mannosidase gene EDEM3 cause a congenital disorder of glycosylation.

EDEM3 encodes a protein that converts Man8GlcNAc2 isomer B to Man7-5GlcNAc2. It is involved in the endoplasmic reticulum-associated degradation pathway, responsible for the recognition of misfolded proteins that will be targeted and translocated to the cytosol and degraded by the proteasome. In this study, through a combination of exome sequencing and gene matching, we have identified seven independent families with 11 individuals with bi-allelic protein-truncating variants and one individual with a compound heterozygous missense variant in EDEM3. The affected individuals present with an inherited congenital disorder of glycosylation (CDG) consisting of neurodevelopmental delay and variable facial dysmorphisms. Experiments in human fibroblast cell lines, human plasma, and mouse plasma and brain tissue demonstrated decreased trimming of Man8GlcNAc2 isomer B to Man7GlcNAc2, consistent with loss of EDEM3 enzymatic activity. In human cells, Man5GlcNAc2 to Man4GlcNAc2 conversion is also diminished with an increase of Glc1Man5GlcNAc2. Furthermore, analysis of the unfolded protein response showed a reduced increase in EIF2AK3 (PERK) expression upon stimulation with tunicamycin as compared to controls, suggesting an impaired unfolded protein response. The aberrant plasma N-glycan profile provides a quick, clinically available test for validating variants of uncertain significance that may be identified by molecular genetic testing. We propose to call this deficiency EDEM3-CDG.

De novo and bi-allelic variants in AP1G1 cause neurodevelopmental disorder with developmental delay, intellectual disability, and epilepsy.

Adaptor protein (AP) complexes mediate selective intracellular vesicular trafficking and polarized localization of somatodendritic proteins in neurons. Disease-causing alleles of various subunits of AP complexes have been implicated in several heritable human disorders, including intellectual disabilities (IDs). Here, we report two bi-allelic (c.737C>A [p.Pro246His] and c.1105A>G [p.Met369Val]) and eight de novo heterozygous variants (c.44G>A [p.Arg15Gln], c.103C>T [p.Arg35Trp], c.104G>A [p.Arg35Gln], c.229delC [p.Gln77Lys∗11], c.399_400del [p.Glu133Aspfs∗37], c.747G>T [p.Gln249His], c.928-2A>C [p.?], and c.2459C>G [p.Pro820Arg]) in AP1G1, encoding gamma-1 subunit of adaptor-related protein complex 1 (AP1γ1), associated with a neurodevelopmental disorder (NDD) characterized by mild to severe ID, epilepsy, and developmental delay in eleven families from different ethnicities. The AP1γ1-mediated adaptor complex is essential for the formation of clathrin-coated intracellular vesicles. In silico analysis and 3D protein modeling simulation predicted alteration of AP1γ1 protein folding for missense variants, which was consistent with the observed altered AP1γ1 levels in heterologous cells. Functional studies of the recessively inherited missense variants revealed no apparent impact on the interaction of AP1γ1 with other subunits of the AP-1 complex but rather showed to affect the endosome recycling pathway. Knocking out ap1g1 in zebrafish leads to severe morphological defect and lethality, which was significantly rescued by injection of wild-type AP1G1 mRNA and not by transcripts encoding the missense variants. Furthermore, microinjection of mRNAs with de novo missense variants in wild-type zebrafish resulted in severe developmental abnormalities and increased lethality. We conclude that de novo and bi-allelic variants in AP1G1 are associated with neurodevelopmental disorder in diverse populations.