KIF26A is mutated in the syndrome of congenital hydrocephalus with megacolon.

Human disorders of the enteric nervous system (ENS), e.g., Hirschsprung's disease, are rarely associated with major central nervous system involvement. We describe two families each segregating a different homozygous truncating variant in KIF26A with a unique constellation of severe megacolon that resembles Hirschsprung's disease but lacks aganglionosis as well as brain malformations that range from severe to mild. The intestinal phenotype bears a striking resemblance to that observed in Kif26a-/- mice where KIF26A deficiency was found to cause abnormal GDNF-Ret signaling resulting in failure to establish normal neuronal networks despite myenteric neuronal hyperplasia. Very recently, a range of brain developmental phenotypes were described in patients and mice with KIF26A deficiency and were found to result from abnormal radial migration and increased apoptosis. Our report, therefore, reveals a recognizable autosomal-recessive human KIF26A deficiency phenotype characterized by severe ENS dysfunction and a range of brain malformations.

An HNRNPK-specific DNA methylation signature makes sense of missense variants and expands the phenotypic spectrum of Au-Kline syndrome.

Au-Kline syndrome (AKS) is a neurodevelopmental disorder associated with multiple malformations and a characteristic facial gestalt. The first individuals ascertained carried de novo loss-of-function (LoF) variants in HNRNPK. Here, we report 32 individuals with AKS (26 previously unpublished), including 13 with de novo missense variants. We propose new clinical diagnostic criteria for AKS that differentiate it from the clinically overlapping Kabuki syndrome and describe a significant phenotypic expansion to include individuals with missense variants who present with subtle facial features and few or no malformations. Many gene-specific DNA methylation (DNAm) signatures have been identified for neurodevelopmental syndromes. Because HNRNPK has roles in chromatin and epigenetic regulation, we hypothesized that pathogenic variants in HNRNPK may be associated with a specific DNAm signature. Here, we report a unique DNAm signature for AKS due to LoF HNRNPK variants, distinct from controls and Kabuki syndrome. This DNAm signature is also identified in some individuals with de novo HNRNPK missense variants, confirming their pathogenicity and the phenotypic expansion of AKS to include more subtle phenotypes. Furthermore, we report that some individuals with missense variants have an "intermediate" DNAm signature that parallels their milder clinical presentation, suggesting the presence of an epi-genotype phenotype correlation. In summary, the AKS DNAm signature may help elucidate the underlying pathophysiology of AKS. This DNAm signature also effectively supported clinical syndrome delineation and is a valuable aid for variant interpretation in individuals where a clinical diagnosis of AKS is unclear, particularly for mild presentations.

Biallelic variants in SLC38A3 encoding a glutamine transporter cause epileptic encephalopathy.

The solute carrier (SLC) superfamily encompasses >400 transmembrane transporters involved in the exchange of amino acids, nutrients, ions, metals, neurotransmitters and metabolites across biological membranes. SLCs are highly expressed in the mammalian brain; defects in nearly 100 unique SLC-encoding genes (OMIM: https://www.omim.org) are associated with rare Mendelian disorders including developmental and epileptic encephalopathy and severe neurodevelopmental disorders. Exome sequencing and family-based rare variant analyses on a cohort with neurodevelopmental disorders identified two siblings with developmental and epileptic encephalopathy and a shared deleterious homozygous splicing variant in SLC38A3. The gene encodes SNAT3, a sodium-coupled neutral amino acid transporter and a principal transporter of the amino acids asparagine, histidine, and glutamine, the latter being the precursor for the neurotransmitters GABA and glutamate. Additional subjects with a similar developmental and epileptic encephalopathy phenotype and biallelic predicted-damaging SLC38A3 variants were ascertained through GeneMatcher and collaborations with research and clinical molecular diagnostic laboratories. Untargeted metabolomic analysis was performed to identify novel metabolic biomarkers. Ten individuals from seven unrelated families from six different countries with deleterious biallelic variants in SLC38A3 were identified. Global developmental delay, intellectual disability, hypotonia, and absent speech were common features while microcephaly, epilepsy, and visual impairment were present in the majority. Epilepsy was drug-resistant in half. Metabolomic analysis revealed perturbations of glutamate, histidine, and nitrogen metabolism in plasma, urine, and CSF of selected subjects, potentially representing biomarkers of disease. Our data support the contention that SLC38A3 is a novel disease gene for developmental and epileptic encephalopathy and illuminate the likely pathophysiology of the disease as perturbations in glutamine homeostasis.

Lethal variants in humans: lessons learned from a large molecular autopsy cohort.

BACKGROUND: Molecular autopsy refers to DNA-based identification of the cause of death. Despite recent attempts to broaden its scope, the term remains typically reserved to sudden unexplained death in young adults. In this study, we aim to showcase the utility of molecular autopsy in defining lethal variants in humans. METHODS: We describe our experience with a cohort of 481 cases in whom the cause of premature death was investigated using DNA from the index or relatives (molecular autopsy by proxy). Molecular autopsy tool was typically exome sequencing although some were investigated using targeted approaches in the earlier stages of the study; these include positional mapping, targeted gene sequencing, chromosomal microarray, and gene panels. RESULTS: The study includes 449 cases from consanguineous families and 141 lacked family history (simplex). The age range was embryos to 18 years. A likely causal variant (pathogenic/likely pathogenic) was identified in 63.8% (307/481), a much higher yield compared to the general diagnostic yield (43%) from the same population. The predominance of recessive lethal alleles allowed us to implement molecular autopsy by proxy in 55 couples, and the yield was similarly high (63.6%). We also note the occurrence of biallelic lethal forms of typically non-lethal dominant disorders, sometimes representing a novel bona fide biallelic recessive disease trait. Forty-six disease genes with no OMIM phenotype were identified in the course of this study. The presented data support the candidacy of two other previously reported novel disease genes (FAAH2 and MSN). The focus on lethal phenotypes revealed many examples of interesting phenotypic expansion as well as remarkable variability in clinical presentation. Furthermore, important insights into population genetics and variant interpretation are highlighted based on the results. CONCLUSIONS: Molecular autopsy, broadly defined, proved to be a helpful clinical approach that provides unique insights into lethal variants and the clinical annotation of the human genome.

The recurrent missense mutation p.(Arg367Trp) in YARS1 causes a distinct neurodevelopmental phenotype.

Pathogenic variants in aminoacyl-tRNA synthetases (ARS1) cause a diverse spectrum of autosomal recessive disorders. Tyrosyl tRNA synthetase (TyrRS) is encoded by YARS1 (cytosolic, OMIM*603,623) and is responsible of coupling tyrosine to its specific tRNA. Next to the enzymatic domain, TyrRS has two additional functional domains (N-Terminal TyrRSMini and C-terminal EMAP-II-like domain) which confer cytokine-like functions. Mutations in YARS1 have been associated with autosomal-dominant Charcot-Marie-Tooth (CMT) neuropathy type C and a heterogenous group of autosomal recessive, multisystem diseases. We identified 12 individuals from 6 families with the recurrent homozygous missense variant c.1099C > T;p.(Arg367Trp) (NM_003680.3) in YARS1. This variant causes a multisystem disorder with developmental delay, microcephaly, failure to thrive, short stature, muscular hypotonia, ataxia, brain anomalies, microcytic anemia, hepatomegaly, and hypothyroidism. In silico analyses show that the p.(Arg367Trp) does not affect the catalytic domain responsible of enzymatic coupling, but destabilizes the cytokine-like C-terminal domain. The phenotype associated with p.(Arg367Trp) is distinct from the other biallelic pathogenic variants that reside in different functional domains of TyrRS which all show some common, but also divergent clinical signs [(e.g., p.(Phe269Ser)-retinal anomalies, p.(Pro213Leu)/p.(Gly525Arg)-mild ID, p.(Pro167Thr)-high fatality)]. The diverse clinical spectrum of ARS1-associated disorders is related to mutations affecting the various non-canonical domains of ARS1, and impaired protein translation is likely not the exclusive disease-causing mechanism of YARS1- and ARS1-associated neurodevelopmental disorders. KEY MESSAGES: The missense variant p.(Arg367Trp) in YARS1 causes a distinct multisystem disorder. p.(Arg367Trp) affects a non-canonical domain with cytokine-like functions. Phenotypic heterogeneity associates with the different affected YARS1 domains. Impaired protein translation is likely not the exclusive mechanism of ARS1-associated disorders.

Missense NAA20 variants impairing the NatB protein N-terminal acetyltransferase cause autosomal recessive developmental delay, intellectual disability, and microcephaly.

PURPOSE: N-terminal acetyltransferases modify proteins by adding an acetyl moiety to the first amino acid and are vital for protein and cell function. The NatB complex acetylates 20% of the human proteome and is composed of the catalytic subunit NAA20 and the auxiliary subunit NAA25. In five individuals with overlapping phenotypes, we identified recessive homozygous missense variants in NAA20. METHODS: Two different NAA20 variants were identified in affected individuals in two consanguineous families by exome and genome sequencing. Biochemical studies were employed to assess the impact of the NAA20 variants on NatB complex formation and catalytic activity. RESULTS: Two homozygous variants, NAA20 p.Met54Val and p.Ala80Val (GenBank: NM_016100.4, c.160A>G and c.239C>T), segregated with affected individuals in two unrelated families presenting with developmental delay, intellectual disability, and microcephaly. Both NAA20-M54V and NAA20-A80V were impaired in their capacity to form a NatB complex with NAA25, and in vitro acetylation assays revealed reduced catalytic activities toward different NatB substrates. Thus, both NAA20 variants are impaired in their ability to perform cellular NatB-mediated N-terminal acetylation. CONCLUSION: We present here a report of pathogenic NAA20 variants causing human disease and data supporting an essential role for NatB-mediated N-terminal acetylation in human development and physiology.

ZNF668 deficiency causes a recognizable disorder of DNA damage repair.

The purpose of this study is to describe a Mendelian disorder of DNA damage repair. Phenotypic delineation of two families, one new and one previously published, with overlapping dysmorphic and neurodevelopmental features was undertaken. Functional characterization of DNA damage repair in fibroblasts obtained from the index individuals in each of the two families was pursued. We present new evidence of a distinct disorder caused by biallelic truncating variants in ZNF668 comprising microcephaly, growth deficiency, severe global developmental delay, brain malformation, and distinct facial dysmorphism. DNA damage repair defect was observed in fibroblasts of affected individuals. ZNF668 deficiency in humans results in a recognizable autosomal recessive disorder, which we propose to name ZNF668-related ZMAND (ZNF668-related brain malformation, microcephaly, abnormal growth, neurodevelopmental delay, and dysmorphism). Our results add to the growing list of Mendelian disorders of the DNA damage repair machinery.

Developmental Consequences of Defective ATG7-Mediated Autophagy in Humans.

BACKGROUND: Autophagy is the major intracellular degradation route in mammalian cells. Systemic ablation of core autophagy-related (ATG) genes in mice leads to embryonic or perinatal lethality, and conditional models show neurodegeneration. Impaired autophagy has been associated with a range of complex human diseases, yet congenital autophagy disorders are rare. METHODS: We performed a genetic, clinical, and neuroimaging analysis involving five families. Mechanistic investigations were conducted with the use of patient-derived fibroblasts, skeletal muscle-biopsy specimens, mouse embryonic fibroblasts, and yeast. RESULTS: We found deleterious, recessive variants in human ATG7, a core autophagy-related gene encoding a protein that is indispensable to classical degradative autophagy. Twelve patients from five families with distinct ATG7 variants had complex neurodevelopmental disorders with brain, muscle, and endocrine involvement. Patients had abnormalities of the cerebellum and corpus callosum and various degrees of facial dysmorphism. These patients have survived with impaired autophagic flux arising from a diminishment or absence of ATG7 protein. Although autophagic sequestration was markedly reduced, evidence of basal autophagy was readily identified in fibroblasts and skeletal muscle with loss of ATG7. Complementation of different model systems by deleterious ATG7 variants resulted in poor or absent autophagic function as compared with the reintroduction of wild-type ATG7. CONCLUSIONS: We identified several patients with a neurodevelopmental disorder who have survived with a severe loss or complete absence of ATG7, an essential effector enzyme for autophagy without a known functional paralogue. (Funded by the Wellcome Centre for Mitochondrial Research and others.).

Mutations in TP73 cause impaired mucociliary clearance and lissencephaly.

TP73 belongs to the TP53 family of transcription factors and has therefore been well studied in cancer research. Studies in mice, however, have revealed non-oncogenic activities related to multiciliogenesis. Utilizing whole-exome sequencing analysis in a cohort of individuals with a mucociliary clearance disorder and cortical malformation, we identified homozygous loss-of-function variants in TP73 in seven individuals from five unrelated families. All affected individuals exhibit a chronic airway disease as well as a brain malformation consistent with lissencephaly. We performed high-speed video microscopy, immunofluorescence analyses, and transmission electron microscopy in respiratory epithelial cells after spheroid or air liquid interface culture to analyze ciliary function, ciliary length, and number of multiciliated cells (MCCs). The respiratory epithelial cells studied display reduced ciliary length and basal bodies mislocalized within the cytoplasm. The number of MCCs is severely reduced, consistent with a reduced number of cells expressing the transcription factors crucial for multiciliogenesis (FOXJ1, RFX2). Our data demonstrate that autosomal-recessive deleterious variants in the TP53 family member TP73 cause a mucociliary clearance disorder due to a defect in MCC differentiation.

Biallelic and monoallelic variants in PLXNA1 are implicated in a novel neurodevelopmental disorder with variable cerebral and eye anomalies.

PURPOSE: To investigate the effect of PLXNA1 variants on the phenotype of patients with autosomal dominant and recessive inheritance patterns and to functionally characterize the zebrafish homologs plxna1a and plxna1b during development. METHODS: We assembled ten patients from seven families with biallelic or de novo PLXNA1 variants. We describe genotype-phenotype correlations, investigated the variants by structural modeling, and used Morpholino knockdown experiments in zebrafish to characterize the embryonic role of plxna1a and plxna1b. RESULTS: Shared phenotypic features among patients include global developmental delay (9/10), brain anomalies (6/10), and eye anomalies (7/10). Notably, seizures were predominantly reported in patients with monoallelic variants. Structural modeling of missense variants in PLXNA1 suggests distortion in the native protein. Our zebrafish studies enforce an embryonic role of plxna1a and plxna1b in the development of the central nervous system and the eye. CONCLUSION: We propose that different biallelic and monoallelic variants in PLXNA1 result in a novel neurodevelopmental syndrome mainly comprising developmental delay, brain, and eye anomalies. We hypothesize that biallelic variants in the extracellular Plexin-A1 domains lead to impaired dimerization or lack of receptor molecules, whereas monoallelic variants in the intracellular Plexin-A1 domains might impair downstream signaling through a dominant-negative effect.

Bi-allelic loss-of-function variants in BCAS3 cause a syndromic neurodevelopmental disorder.

BCAS3 microtubule-associated cell migration factor (BCAS3) is a large, highly conserved cytoskeletal protein previously proposed to be critical in angiogenesis and implicated in human embryogenesis and tumorigenesis. Here, we established BCAS3 loss-of-function variants as causative for a neurodevelopmental disorder. We report 15 individuals from eight unrelated families with germline bi-allelic loss-of-function variants in BCAS3. All probands share a global developmental delay accompanied by pyramidal tract involvement, microcephaly, short stature, strabismus, dysmorphic facial features, and seizures. The human phenotype is less severe compared with the Bcas3 knockout mouse model and cannot be explained by angiogenic defects alone. Consistent with being loss-of-function alleles, we observed absence of BCAS3 in probands' primary fibroblasts. By comparing the transcriptomic and proteomic data based on probands' fibroblasts with those of the knockout mouse model, we identified similar dysregulated pathways resulting from over-representation analysis, while the dysregulation of some proposed key interactors could not be confirmed. Together with the results from a tissue-specific Drosophila loss-of-function model, we demonstrate a vital role for BCAS3 in neural tissue development.

MYH1 is a candidate gene for recurrent rhabdomyolysis in humans.

Rhabdomyolysis is a serious medical condition characterized by muscle injury, and there are recognized genetic causes especially in recurrent forms. The majority of these cases, however, remain unexplained. Here, we describe a patient with recurrent rhabdomyolysis in whom extensive clinical testing failed to identify a likely etiology. Whole-exome sequencing revealed a novel missense variant in MYH1, which encodes a major adult muscle fiber protein. Structural biology analysis revealed that the mutated residue is extremely well conserved and is located in the actin binding cleft. Furthermore, immediately adjacent mutations in that cleft in other myosins are pathogenic in humans. Our results are consistent with the finding that MYH1 is mutated in rhabdomyolysis in horses and suggest that this gene should be investigated in cases with recurrent rhabdomyolysis.

Biallelic UBE4A loss-of-function variants cause intellectual disability and global developmental delay.

PURPOSE: To identify novel genes associated with intellectual disability (ID) in four unrelated families. METHODS: Here, through exome sequencing and international collaboration, we report eight individuals from four unrelated families of diverse geographic origin with biallelic loss-of-function variants in UBE4A. RESULTS: Eight evaluated individuals presented with syndromic intellectual disability and global developmental delay. Other clinical features included hypotonia, short stature, seizures, and behavior disorder. Characteristic features were appreciated in some individuals but not all; in some cases, features became more apparent with age. We demonstrated that UBE4A loss-of-function variants reduced RNA expression and protein levels in clinical samples. Mice generated to mimic patient-specific Ube4a loss-of-function variant exhibited muscular and neurological/behavioral abnormalities, some of which are suggestive of the clinical abnormalities seen in the affected individuals. CONCLUSION: These data indicate that biallelic loss-of-function variants in UBE4A cause a novel intellectual disability syndrome, suggesting that UBE4A enzyme activity is required for normal development and neurological function.

KDM5A mutations identified in autism spectrum disorder using forward genetics.

Autism spectrum disorder (ASD) is a constellation of neurodevelopmental disorders with high phenotypic and genetic heterogeneity, complicating the discovery of causative genes. Through a forward genetics approach selecting for defective vocalization in mice, we identified Kdm5a as a candidate ASD gene. To validate our discovery, we generated a Kdm5a knockout mouse model (Kdm5a-/-) and confirmed that inactivating Kdm5a disrupts vocalization. In addition, Kdm5a-/- mice displayed repetitive behaviors, sociability deficits, cognitive dysfunction, and abnormal dendritic morphogenesis. Loss of KDM5A also resulted in dysregulation of the hippocampal transcriptome. To determine if KDM5A mutations cause ASD in humans, we screened whole exome sequencing and microarray data from a clinical cohort. We identified pathogenic KDM5A variants in nine patients with ASD and lack of speech. Our findings illustrate the power and efficacy of forward genetics in identifying ASD genes and highlight the importance of KDM5A in normal brain development and function.

Further delineation of HIDEA syndrome.

Recently, the genetic cause of HIDEA syndrome (hypotonia, hypoventilation, intellectual disability, dysautonomia, epilepsy, and eye abnormalities) was identified as biallelic pathogenic variants in P4HTM, which encodes an atypical member of the prolyl 4-hydroxylases (P4Hs) family of enzymes. We report seven patients from four new families in whom HIDEA was only diagnosed after whole-exome sequencing (WES) revealed novel disease-causing variants in P4HTM. We note the variable phenotypic expressivity of the syndrome except for cognitive impairment/developmental delay, and hypotonia, which seem to be consistent findings. One patient only presented with hypotonia, developmental delay, and abnormal eye movements, which highlights the challenge in diagnosing milder cases with this new syndrome. Other notable features include mild facial dysmorphism, obesity, and brain dysmyelination and atrophy. We conclude that HIDEA is a highly variable syndrome and suspect that a large fraction of patients will be diagnosed via reverse phenotyping after recessive P4HTM variants are identified by agnostic genomic sequencing assays.

Early-infantile onset epilepsy and developmental delay caused by bi-allelic GAD1 variants.

Gamma-aminobutyric acid (GABA) and glutamate are the most abundant amino acid neurotransmitters in the brain. GABA, an inhibitory neurotransmitter, is synthesized by glutamic acid decarboxylase (GAD). Its predominant isoform GAD67, contributes up to ∼90% of base-level GABA in the CNS, and is encoded by the GAD1 gene. Disruption of GAD1 results in an imbalance of inhibitory and excitatory neurotransmitters, and as Gad1-/- mice die neonatally of severe cleft palate, it has not been possible to determine any potential neurological dysfunction. Furthermore, little is known about the consequence of GAD1 disruption in humans. Here we present six affected individuals from six unrelated families, carrying bi-allelic GAD1 variants, presenting with developmental and epileptic encephalopathy, characterized by early-infantile onset epilepsy and hypotonia with additional variable non-CNS manifestations such as skeletal abnormalities, dysmorphic features and cleft palate. Our findings highlight an important role for GAD1 in seizure induction, neuronal and extraneuronal development, and introduce GAD1 as a new gene associated with developmental and epileptic encephalopathy.

Biallelic MFSD2A variants associated with congenital microcephaly, developmental delay, and recognizable neuroimaging features.

Major Facilitator Superfamily Domain containing 2a (MFSD2A) is an essential endothelial lipid transporter at the blood-brain barrier. Biallelic variants affecting function in MFSD2A cause autosomal recessive primary microcephaly 15 (MCPH15, OMIM# 616486). We sought to expand our knowledge of the phenotypic spectrum of MCPH15 and demonstrate the underlying mechanism of inactivation of the MFSD2A transporter. We carried out detailed analysis of the clinical and neuroradiological features of a series of 27 MCPH15 cases, including eight new individuals from seven unrelated families. Genetic investigation was performed through exome sequencing (ES). Structural insights on the human Mfsd2a model and in-vitro biochemical assays were used to investigate the functional impact of the identified variants. All patients had primary microcephaly and severe developmental delay. Brain MRI showed variable degrees of white matter reduction, ventricular enlargement, callosal hypodysgenesis, and pontine and vermian hypoplasia. ES led to the identification of six novel biallelic MFSD2A variants (NG_053084.1, NM_032793.5: c.556+1G>A, c.748G>T; p.(Val250Phe), c.750_753del; p.(Cys251SerfsTer3), c.977G>A; p.(Arg326His), c.1386_1435del; p.(Gln462HisfsTer17), and c.1478C>T; p.(Pro493Leu)) and two recurrent variants (NM_032793.5: c.593C>T; p.(Thr198Met) and c.476C>T; p.(Thr159Met)). All these variants and the previously reported NM_032793.5: c.490C>A; p.(Pro164Thr) resulted in either reduced MFSD2A expression and/or transport activity. Our study further delineates the phenotypic spectrum of MCPH15, refining its clinical and neuroradiological characterization and supporting that MFSD2A deficiency causes early prenatal brain developmental disruption. We also show that poor MFSD2A expression despite normal transporter activity is a relevant pathomechanism in MCPH15.

Analysis of transcript-deleterious variants in Mendelian disorders: implications for RNA-based diagnostics.

BACKGROUND: At least 50% of patients with suspected Mendelian disorders remain undiagnosed after whole-exome sequencing (WES), and the extent to which non-coding variants that are not captured by WES contribute to this fraction is unclear. Whole transcriptome sequencing is a promising supplement to WES, although empirical data on the contribution of RNA analysis to the diagnosis of Mendelian diseases on a large scale are scarce. RESULTS: Here, we describe our experience with transcript-deleterious variants (TDVs) based on a cohort of 5647 families with suspected Mendelian diseases. We first interrogate all families for which the respective Mendelian phenotype could be mapped to a single locus to obtain an unbiased estimate of the contribution of TDVs at 18.9%. We examine the entire cohort and find that TDVs account for 15% of all "solved" cases. We compare the results of RT-PCR to in silico prediction. Definitive results from RT-PCR are obtained from blood-derived RNA for the overwhelming majority of variants (84.1%), and only a small minority (2.6%) fail analysis on all available RNA sources (blood-, skin fibroblast-, and urine renal epithelial cells-derived), which has important implications for the clinical application of RNA-seq. We also show that RNA analysis can establish the diagnosis in 13.5% of 155 patients who had received "negative" clinical WES reports. Finally, our data suggest a role for TDVs in modulating penetrance even in otherwise highly penetrant Mendelian disorders. CONCLUSIONS: Our results provide much needed empirical data for the impending implementation of diagnostic RNA-seq in conjunction with genome sequencing.