De novo missense variants in exon 9 of SEPHS1 cause a neurodevelopmental condition with developmental delay, poor growth, hypotonia, and dysmorphic features.

Selenophosphate synthetase (SEPHS) plays an essential role in selenium metabolism. Two mammalian SEPHS paralogues, SEPHS1 and SEPHS2, share high sequence identity and structural homology with SEPHS. Here, we report nine individuals from eight families with developmental delay, growth and feeding problems, hypotonia, and dysmorphic features, all with heterozygous missense variants in SEPHS1. Eight of these individuals had a recurrent variant at amino acid position 371 of SEPHS1 (p.Arg371Trp, p.Arg371Gln, and p.Arg371Gly); seven of these variants were known to be de novo. Structural modeling and biochemical assays were used to understand the effect of these variants on SEPHS1 function. We found that a variant at residue Trp352 results in local structural changes of the C-terminal region of SEPHS1 that decrease the overall thermal stability of the enzyme. In contrast, variants of a solvent-exposed residue Arg371 do not impact enzyme stability and folding but could modulate direct protein-protein interactions of SEPSH1 with cellular factors in promoting cell proliferation and development. In neuronal SH-SY5Y cells, we assessed the impact of SEPHS1 variants on cell proliferation and ROS production and investigated the mRNA expression levels of genes encoding stress-related selenoproteins. Our findings provided evidence that the identified SEPHS1 variants enhance cell proliferation by modulating ROS homeostasis. Our study supports the hypothesis that SEPHS1 plays a critical role during human development and provides a basis for further investigation into the molecular mechanisms employed by SEPHS1. Furthermore, our data suggest that variants in SEPHS1 are associated with a neurodevelopmental disorder.

Assortative mating and parental genetic relatedness contribute to the pathogenicity of variably expressive variants.

We examined more than 97,000 families from four neurodevelopmental disease cohorts and the UK Biobank to identify phenotypic and genetic patterns in parents contributing to neurodevelopmental disease risk in children. We identified within- and cross-disorder correlations between six phenotypes in parents and children, such as obsessive-compulsive disorder (R = 0.32-0.38, p < 10-126). We also found that measures of sub-clinical autism features in parents are associated with several autism severity measures in children, including biparental mean Social Responsiveness Scale scores and proband Repetitive Behaviors Scale scores (regression coefficient = 0.14, p = 3.38 × 10-4). We further describe patterns of phenotypic similarity between spouses, where spouses show correlations for six neurological and psychiatric phenotypes, including a within-disorder correlation for depression (R = 0.24-0.68, p < 0.001) and a cross-disorder correlation between anxiety and bipolar disorder (R = 0.09-0.22, p < 10-92). Using a simulated population, we also found that assortative mating can lead to increases in disease liability over generations and the appearance of "genetic anticipation" in families carrying rare variants. We identified several families in a neurodevelopmental disease cohort where the proband inherited multiple rare variants in disease-associated genes from each of their affected parents. We further identified parental relatedness as a risk factor for neurodevelopmental disorders through its inverse relationship with variant pathogenicity and propose that parental relatedness modulates disease risk by increasing genome-wide homozygosity in children (R = 0.05-0.26, p < 0.05). Our results highlight the utility of assessing parent phenotypes and genotypes toward predicting features in children who carry rare variably expressive variants and implicate assortative mating as a risk factor for increased disease severity in these families.

Evidence for 28 genetic disorders discovered by combining healthcare and research data.

De novo mutations in protein-coding genes are a well-established cause of developmental disorders1. However, genes known to be associated with developmental disorders account for only a minority of the observed excess of such de novo mutations1,2. Here, to identify previously undescribed genes associated with developmental disorders, we integrate healthcare and research exome-sequence data from 31,058 parent-offspring trios of individuals with developmental disorders, and develop a simulation-based statistical test to identify gene-specific enrichment of de novo mutations. We identified 285 genes that were significantly associated with developmental disorders, including 28 that had not previously been robustly associated with developmental disorders. Although we detected more genes associated with developmental disorders, much of the excess of de novo mutations in protein-coding genes remains unaccounted for. Modelling suggests that more than 1,000 genes associated with developmental disorders have not yet been described, many of which are likely to be less penetrant than the currently known genes. Research access to clinical diagnostic datasets will be critical for completing the map of genes associated with developmental disorders.

TRIM71 mutations cause a neurodevelopmental syndrome featuring ventriculomegaly and hydrocephalus.

Congenital hydrocephalus (CH), characterized by cerebral ventriculomegaly, is one of the most common reasons for pediatric brain surgery. Recent studies have implicated lin-41 (lineage variant 41)/TRIM71 (tripartite motif 71) as a candidate CH risk gene, however, TRIM71 variants have not been systematically examined in a large patient cohort or conclusively linked with an OMIM syndrome. Through cross-sectional analysis of the largest assembled cohort of patients with cerebral ventriculomegaly, including neurosurgically-treated CH (totaling 2,697 parent-proband trios and 8,091 total exomes), we identified 13 protein-altering de novo variants (DNVs) in TRIM71 in unrelated children exhibiting variable ventriculomegaly, CH, developmental delay, dysmorphic features, and other structural brain defects including corpus callosum dysgenesis and white matter hypoplasia. Eight unrelated patients were found to harbor arginine variants, including two recurrent missense DNVs, at homologous positions in RPXGV motifs of different NHL domains. Seven additional patients with rare, damaging, unphased or transmitted variants of uncertain significance were also identified. NHL-domain variants of TRIM71 exhibited impaired binding to the canonical TRIM71 target CDKN1A; other variants failed to direct the subcellular localization of TRIM71 to processing bodies. Single-cell transcriptomic analysis of human embryos revealed expression of TRIM71 in early first-trimester neural stem cells of the brain. These data show TRIM71 is essential for human brain morphogenesis and that TRIM71 mutations cause a novel neurodevelopmental syndrome featuring ventriculomegaly and CH.

Non-coding region variants upstream of MEF2C cause severe developmental disorder through three distinct loss-of-function mechanisms.

Clinical genetic testing of protein-coding regions identifies a likely causative variant in only around half of developmental disorder (DD) cases. The contribution of regulatory variation in non-coding regions to rare disease, including DD, remains very poorly understood. We screened 9,858 probands from the Deciphering Developmental Disorders (DDD) study for de novo mutations in the 5' untranslated regions (5' UTRs) of genes within which variants have previously been shown to cause DD through a dominant haploinsufficient mechanism. We identified four single-nucleotide variants and two copy-number variants upstream of MEF2C in a total of ten individual probands. We developed multiple bespoke and orthogonal experimental approaches to demonstrate that these variants cause DD through three distinct loss-of-function mechanisms, disrupting transcription, translation, and/or protein function. These non-coding region variants represent 23% of likely diagnoses identified in MEF2C in the DDD cohort, but these would all be missed in standard clinical genetics approaches. Nonetheless, these variants are readily detectable in exome sequence data, with 30.7% of 5' UTR bases across all genes well covered in the DDD dataset. Our analyses show that non-coding variants upstream of genes within which coding variants are known to cause DD are an important cause of severe disease and demonstrate that analyzing 5' UTRs can increase diagnostic yield. We also show how non-coding variants can help inform both the disease-causing mechanism underlying protein-coding variants and dosage tolerance of the gene.

Variants in the degron of AFF3 are associated with intellectual disability, mesomelic dysplasia, horseshoe kidney, and epileptic encephalopathy.

The ALF transcription factor paralogs, AFF1, AFF2, AFF3, and AFF4, are components of the transcriptional super elongation complex that regulates expression of genes involved in neurogenesis and development. We describe an autosomal dominant disorder associated with de novo missense variants in the degron of AFF3, a nine amino acid sequence important for its binding to ubiquitin ligase, or with de novo deletions of this region. The sixteen affected individuals we identified, along with two previously reported individuals, present with a recognizable pattern of anomalies, which we named KINSSHIP syndrome (KI for horseshoe kidney, NS for Nievergelt/Savarirayan type of mesomelic dysplasia, S for seizures, H for hypertrichosis, I for intellectual disability, and P for pulmonary involvement), partially overlapping the AFF4-associated CHOPS syndrome. Whereas homozygous Aff3 knockout mice display skeletal anomalies, kidney defects, brain malformations, and neurological anomalies, knockin animals modeling one of the microdeletions and the most common of the missense variants identified in affected individuals presented with lower mesomelic limb deformities like KINSSHIP-affected individuals and early lethality, respectively. Overexpression of AFF3 in zebrafish resulted in body axis anomalies, providing some support for the pathological effect of increased amount of AFF3. The only partial phenotypic overlap of AFF3- and AFF4-associated syndromes and the previously published transcriptome analyses of ALF transcription factors suggest that these factors are not redundant and each contributes uniquely to proper development.

Heterozygous ANKRD17 loss-of-function variants cause a syndrome with intellectual disability, speech delay, and dysmorphism.

ANKRD17 is an ankyrin repeat-containing protein thought to play a role in cell cycle progression, whose ortholog in Drosophila functions in the Hippo pathway as a co-factor of Yorkie. Here, we delineate a neurodevelopmental disorder caused by de novo heterozygous ANKRD17 variants. The mutational spectrum of this cohort of 34 individuals from 32 families is highly suggestive of haploinsufficiency as the underlying mechanism of disease, with 21 truncating or essential splice site variants, 9 missense variants, 1 in-frame insertion-deletion, and 1 microdeletion (1.16 Mb). Consequently, our data indicate that loss of ANKRD17 is likely the main cause of phenotypes previously associated with large multi-gene chromosomal aberrations of the 4q13.3 region. Protein modeling suggests that most of the missense variants disrupt the stability of the ankyrin repeats through alteration of core structural residues. The major phenotypic characteristic of our cohort is a variable degree of developmental delay/intellectual disability, particularly affecting speech, while additional features include growth failure, feeding difficulties, non-specific MRI abnormalities, epilepsy and/or abnormal EEG, predisposition to recurrent infections (mostly bacterial), ophthalmological abnormalities, gait/balance disturbance, and joint hypermobility. Moreover, many individuals shared similar dysmorphic facial features. Analysis of single-cell RNA-seq data from the developing human telencephalon indicated ANKRD17 expression at multiple stages of neurogenesis, adding further evidence to the assertion that damaging ANKRD17 variants cause a neurodevelopmental disorder.

Non-immune hydrops fetalis is associated with bi-allelic pathogenic variants in the MYB Binding Protein 1a (MYBBP1A) gene.

Non-immune hydrops fetalis (NIHF) is a rare entity characterized by excessive accumulation of fluid within the fetal extravascular compartments and body cavities. Here we present two intrauterine fetal demises with NIHF presenting with oligohydramnios, cystic hygroma, pleural effusion, and generalized hydrops with predominance of subcutaneous edema. The fetuses also presented with ascites, severe and precocious IUGR and skeletal anomalies. Whole exome sequencing was applied in order to screen for a possible genetic cause. The results identified biallelic variants in MYBBP1A in both fetuses. A previous report described another case with a similar phenotype having compound heterozygous variants in the same gene. The protein encoded by MYBBP1A is involved in several cellular processes including the synthesis of ribosomal DNA, the response to nucleolar stress, and tumor suppression. Our functional protein analysis through immunohistochemistry indicates that MYBBP1A is a gene expressed during fetal stages. Altogether, we concluded that MYBBP1A is associated with the development of hydrops fetalis. More cases and further studies are necessary to understand the role of this gene and the mechanism associated with NIHF.

Bi-allelic ADARB1 Variants Associated with Microcephaly, Intellectual Disability, and Seizures.

The RNA editing enzyme ADAR2 is essential for the recoding of brain transcripts. Impaired ADAR2 editing leads to early-onset epilepsy and premature death in a mouse model. Here, we report bi-allelic variants in ADARB1, the gene encoding ADAR2, in four unrelated individuals with microcephaly, intellectual disability, and epilepsy. In one individual, a homozygous variant in one of the double-stranded RNA-binding domains (dsRBDs) was identified. In the others, variants were situated in or around the deaminase domain. To evaluate the effects of these variants on ADAR2 enzymatic activity, we performed in vitro assays with recombinant proteins in HEK293T cells and ex vivo assays with fibroblasts derived from one of the individuals. We demonstrate that these ADAR2 variants lead to reduced editing activity on a known ADAR2 substrate. We also demonstrate that one variant leads to changes in splicing of ADARB1 transcript isoforms. These findings reinforce the importance of RNA editing in brain development and introduce ADARB1 as a genetic etiology in individuals with intellectual disability, microcephaly, and epilepsy.

Variants in SCAF4 Cause a Neurodevelopmental Disorder and Are Associated with Impaired mRNA Processing.

RNA polymerase II interacts with various other complexes and factors to ensure correct initiation, elongation, and termination of mRNA transcription. One of these proteins is SR-related CTD-associated factor 4 (SCAF4), which is important for correct usage of polyA sites for mRNA termination. Using exome sequencing and international matchmaking, we identified nine likely pathogenic germline variants in SCAF4 including two splice-site and seven truncating variants, all residing in the N-terminal two thirds of the protein. Eight of these variants occurred de novo, and one was inherited. Affected individuals demonstrated a variable neurodevelopmental disorder characterized by mild intellectual disability, seizures, behavioral abnormalities, and various skeletal and structural anomalies. Paired-end RNA sequencing on blood lymphocytes of SCAF4-deficient individuals revealed a broad deregulation of more than 9,000 genes and significant differential splicing of more than 2,900 genes, indicating an important role of SCAF4 in mRNA processing. Knockdown of the SCAF4 ortholog CG4266 in the model organism Drosophila melanogaster resulted in impaired locomotor function, learning, and short-term memory. Furthermore, we observed an increased number of active zones in larval neuromuscular junctions, representing large glutamatergic synapses. These observations indicate a role of CG4266 in nervous system development and function and support the implication of SCAF4 in neurodevelopmental phenotypes. In summary, our data show that heterozygous, likely gene-disrupting variants in SCAF4 are causative for a variable neurodevelopmental disorder associated with impaired mRNA processing.

De Novo Variants in the ATPase Module of MORC2 Cause a Neurodevelopmental Disorder with Growth Retardation and Variable Craniofacial Dysmorphism.

MORC2 encodes an ATPase that plays a role in chromatin remodeling, DNA repair, and transcriptional regulation. Heterozygous variants in MORC2 have been reported in individuals with autosomal-dominant Charcot-Marie-Tooth disease type 2Z and spinal muscular atrophy, and the onset of symptoms ranges from infancy to the second decade of life. Here, we present a cohort of 20 individuals referred for exome sequencing who harbor pathogenic variants in the ATPase module of MORC2. Individuals presented with a similar phenotype consisting of developmental delay, intellectual disability, growth retardation, microcephaly, and variable craniofacial dysmorphism. Weakness, hyporeflexia, and electrophysiologic abnormalities suggestive of neuropathy were frequently observed but were not the predominant feature. Five of 18 individuals for whom brain imaging was available had lesions reminiscent of those observed in Leigh syndrome, and five of six individuals who had dilated eye exams had retinal pigmentary abnormalities. Functional assays revealed that these MORC2 variants result in hyperactivation of epigenetic silencing by the HUSH complex, supporting their pathogenicity. The described set of morphological, growth, developmental, and neurological findings and medical concerns expands the spectrum of genetic disorders resulting from pathogenic variants in MORC2.

Functionally impaired RPL8 variants associated with Diamond-Blackfan anemia and a Diamond-Blackfan anemia-like phenotype.

Diamond-Blackfan anemia is a rare genetic disease characterized by erythroblastopenia and a large spectrum of developmental anomalies. The vast majority of the cases genetically described are linked to heterozygous pathogenic variants in more than 20 ribosomal protein genes. Here we report an atypical clinical case of DBA associated with a missense variant in RPL8, which encodes RPL8/uL2, a protein of the 60S large ribosomal subunit. RPL8 has been previously implicated as a candidate disease gene in one patient with DBA bearing another type of missense variant; however, evidence for pathogenicity was limited to computational tools. Using functional studies in lymphoblastoid cells as well as yeast models, we show that the RPL8 variants detected in these two patients encode functionally deficient proteins that affect ribosome production and are therefore likely pathogenic. We propose to include RPL8 in the list of DBA-associated genes.

De novo variants in MPP5 cause global developmental delay and behavioral changes.

Membrane Protein Palmitoylated 5 (MPP5) is a highly conserved apical complex protein essential for cell polarity, fate and survival. Defects in cell polarity are associated with neurologic disorders including autism and microcephaly. MPP5 is essential for neurogenesis in animal models, but human variants leading to neurologic impairment have not been described. We identified three patients with heterozygous MPP5 de novo variants (DNV) and global developmental delay (GDD) and compared their phenotypes and magnetic resonance imaging (MRI) to ascertain how MPP5 DNV leads to GDD. All three patients with MPP5 DNV experienced GDD with language delay/regression and behavioral changes. MRI ranged from normal to decreased gyral folding and microcephaly. The effects of MPP5 depletion on the developing brain were assessed by creating a heterozygous conditional knock out (het CKO) murine model with central nervous system (CNS)-specific Nestin-Cre drivers. In the het CKO model, Mpp5 depletion led to microcephaly, decreased cerebellar volume and cortical thickness. Het CKO mice had decreased ependymal cells and Mpp5 at the apical surface of cortical ventricular zone compared with wild type. Het CKO mice also failed to maintain progenitor pools essential for neurogenesis. The proportion of cortical cells undergoing apoptotic cell death increased, suggesting that cell death reduces progenitor population and neuron number. Het CKO mice also showed behavioral changes, similar to our patients. To our knowledge, this is the first report to show that variants in MPP5 are associated with GDD, behavioral abnormalities and language regression/delay. Murine modeling shows that neurogenesis is likely altered in these individuals, with cell death and skewed cellular composition playing significant roles.

Loss of Oxidation Resistance 1, OXR1, Is Associated with an Autosomal-Recessive Neurological Disease with Cerebellar Atrophy and Lysosomal Dysfunction.

We report an early-onset autosomal-recessive neurological disease with cerebellar atrophy and lysosomal dysfunction. We identified bi-allelic loss-of-function (LoF) variants in Oxidative Resistance 1 (OXR1) in five individuals from three families; these individuals presented with a history of severe global developmental delay, current intellectual disability, language delay, cerebellar atrophy, and seizures. While OXR1 is known to play a role in oxidative stress resistance, its molecular functions are not well established. OXR1 contains three conserved domains: LysM, GRAM, and TLDc. The gene encodes at least six transcripts, including some that only consist of the C-terminal TLDc domain. We utilized Drosophila to assess the phenotypes associated with loss of mustard (mtd), the fly homolog of OXR1. Strong LoF mutants exhibit late pupal lethality or pupal eclosion defects. Interestingly, although mtd encodes 26 transcripts, severe LoF and null mutations can be rescued by a single short human OXR1 cDNA that only contains the TLDc domain. Similar rescue is observed with the TLDc domain of NCOA7, another human homolog of mtd. Loss of mtd in neurons leads to massive cell loss, early death, and an accumulation of aberrant lysosomal structures, similar to what we observe in fibroblasts of affected individuals. Our data indicate that mtd and OXR1 are required for proper lysosomal function; this is consistent with observations that NCOA7 is required for lysosomal acidification.

Heterozygous variants in MYH10 associated with neurodevelopmental disorders and congenital anomalies with evidence for primary cilia-dependent defects in Hedgehog signaling.

PURPOSE: Nonmuscle myosin II complexes are master regulators of actin dynamics that play essential roles during embryogenesis with vertebrates possessing 3 nonmuscle myosin II heavy chain genes, MYH9, MYH10, and MYH14. As opposed to MYH9 and MYH14, no recognizable disorder has been associated with MYH10. We sought to define the clinical characteristics and molecular mechanism of a novel autosomal dominant disorder related to MYH10. METHODS: An international collaboration identified the patient cohort. CAS9-mediated knockout cell models were used to explore the mechanism of disease pathogenesis. RESULTS: We identified a cohort of 16 individuals with heterozygous MYH10 variants presenting with a broad spectrum of neurodevelopmental disorders and variable congenital anomalies that affect most organ systems and were recapitulated in animal models of altered MYH10 activity. Variants were typically de novo missense changes with clustering observed in the motor domain. MYH10 knockout cells showed defects in primary ciliogenesis and reduced ciliary length with impaired Hedgehog signaling. MYH10 variant overexpression produced a dominant-negative effect on ciliary length. CONCLUSION: These data presented a novel genetic cause of isolated and syndromic neurodevelopmental disorders related to heterozygous variants in the MYH10 gene with implications for disrupted primary cilia length control and altered Hedgehog signaling in disease pathogenesis.

Phenotypic continuum between POLE-related recessive disorders: A case report and literature review.

POLE is a pleiotropic gene with phenotypic expression of pathogenic variants depending on the type of variant, impact on the protein, and mode of inheritance. Heterozygous missense variants located within the exonuclease domain have been shown to result in polymerase proofreading-associated polyposis (PPAP) which is characterized by an increased risk for colon polyps and colorectal cancer. Biallelic variants resulting in markedly reduced amounts of normal protein have been reported in two separate recessive pediatric syndromes: facial dysmorphism, immunodeficiency, livedo, and short stature as well as intrauterine growth restriction, metaphyseal dysplasia, adrenal hypoplasia congenital, and genital anomalies. Here we report two siblings identified to have POLE c.1686 + 32C > G in trans with POLE p.(Glu709*) via exome sequencing. A detailed review of the reported phenotypes in these two siblings and from available literature revealed that individuals with biallelic POLE pathogenic variants resulting in partial loss-of-function present with a similar phenotype: short stature and facial dysmorphism with or without immunodeficiency. These data suggest a phenotypic continuum between the previously reported POLE-related recessive disorders.

First case of pan-suture craniosynostosis due to de novo mosaic KAT6A mutation.

A nonverbal 3-year-old male with a complex past medical history was referred to pediatric neurosurgery for evaluation of Chiari I malformation. A full clinical evaluation suggested that the "Chiari" was a secondary change caused by craniocerebral disproportion that was the result of delayed pan-sutural craniosynostosis. Given his unknown cause of craniosynostosis, whole-exome sequencing (WES) was performed. WES revealed a de novo, somatic mosaic variant in the KAT6A gene. This report discusses importance of keeping a broad differential in the setting of referral for Chiari I malformation and presents a unique case of craniosynostosis. Additionally, it emphasizes the value of utilizing genetic testing for complex craniofacial cases with unknown causes to provide clinical answers and guide clinical management.

The ubiquitin ligase UBE3B, disrupted in intellectual disability and absent speech, regulates metabolic pathways by targeting BCKDK.

Kaufman oculocerebrofacial syndrome (KOS) is a recessive neurodevelopmental disorder characterized by intellectual disability and lack of speech. KOS is caused by inactivating mutations in UBE3B, but the underlying biological mechanisms are completely unknown. We found that loss of Ube3b in mice resulted in growth retardation, decreased grip strength, and loss of vocalization. The brains of Ube3b-/- mice had hypoplasia of the corpus callosum, enlarged ventricles, and decreased thickness of the somatosensory cortex. Ube3b-/- cortical neurons had abnormal dendritic morphology and synapses. We identified 22 UBE3B interactors and found that branched-chain α-ketoacid dehydrogenase kinase (BCKDK) is an in vivo UBE3B substrate. Since BCKDK targets several metabolic pathways, we profiled plasma and cortical metabolomes from Ube3b-/- mice. Nucleotide metabolism and the tricarboxylic acid cycle were among the pathways perturbed. Substrate-induced mitochondrial respiration was reduced in skeletal muscle but not in liver of Ube3b-/- mice. To assess the relevance of these findings to humans, we identified three KOS patients who had compound heterozygous UBE3B mutations. We discovered changes in metabolites from similar pathways in plasma from these patients. Collectively, our results implicate a disease mechanism in KOS, suggest that it is a metabolic encephalomyopathy, and provide an entry to targeted therapies.

TAOK1 is associated with neurodevelopmental disorder and essential for neuronal maturation and cortical development.

Thousand and one amino-acid kinase 1 (TAOK1) is a MAP3K protein kinase, regulating different mitogen-activated protein kinase pathways, thereby modulating a multitude of processes in the cell. Given the recent finding of TAOK1 involvement in neurodevelopmental disorders (NDDs), we investigated the role of TAOK1 in neuronal function and collected a cohort of 23 individuals with mostly de novo variants in TAOK1 to further define the associated NDD. Here, we provide evidence for an important role for TAOK1 in neuronal function, showing that altered TAOK1 expression levels in the embryonic mouse brain affect neural migration in vivo, as well as neuronal maturation in vitro. The molecular spectrum of the identified TAOK1 variants comprises largely truncating and nonsense variants, but also missense variants, for which we provide evidence that they can have a loss of function or dominant-negative effect on TAOK1, expanding the potential underlying causative mechanisms resulting in NDD. Taken together, our data indicate that TAOK1 activity needs to be properly controlled for normal neuronal function and that TAOK1 dysregulation leads to a neurodevelopmental disorder mainly comprising similar facial features, developmental delay/intellectual disability and/or variable learning or behavioral problems, muscular hypotonia, infant feeding difficulties, and growth problems.