Web-accessible application for identifying pathogenic transcripts with RNA-seq: Increased sensitivity in diagnosis of neurodevelopmental disorders.

For neurodevelopmental disorders (NDDs), a molecular diagnosis is key for management, predicting outcome, and counseling. Often, routine DNA-based tests fail to establish a genetic diagnosis in NDDs. Transcriptome analysis (RNA sequencing [RNA-seq]) promises to improve the diagnostic yield but has not been applied to NDDs in routine diagnostics. Here, we explored the diagnostic potential of RNA-seq in 96 individuals including 67 undiagnosed subjects with NDDs. We performed RNA-seq on single individuals' cultured skin fibroblasts, with and without cycloheximide treatment, and used modified OUTRIDER Z scores to detect gene expression outliers and mis-splicing by exonic and intronic outliers. Analysis was performed by a user-friendly web application, and candidate pathogenic transcriptional events were confirmed by secondary assays. We identified intragenic deletions, monoallelic expression, and pseudoexonic insertions but also synonymous and non-synonymous variants with deleterious effects on transcription, increasing the diagnostic yield for NDDs by 13%. We found that cycloheximide treatment and exonic/intronic Z score analysis increased detection and resolution of aberrant splicing. Importantly, in one individual mis-splicing was found in a candidate gene nearly matching the individual's specific phenotype. However, pathogenic splicing occurred in another neuronal-expressed gene and provided a molecular diagnosis, stressing the need to customize RNA-seq. Lastly, our web browser application allowed custom analysis settings that facilitate diagnostic application and ranked pathogenic transcripts as top candidates. Our results demonstrate that RNA-seq is a complementary method in the genomic diagnosis of NDDs and, by providing accessible analysis with improved sensitivity, our transcriptome analysis approach facilitates wider implementation of RNA-seq in routine genome diagnostics.

Rare germline heterozygous missense variants in BRCA1-associated protein 1, BAP1, cause a syndromic neurodevelopmental disorder.

Nuclear deubiquitinase BAP1 (BRCA1-associated protein 1) is a core component of multiprotein complexes that promote transcription by reversing the ubiquitination of histone 2A (H2A). BAP1 is a tumor suppressor whose germline loss-of-function variants predispose to cancer. To our knowledge, there are very rare examples of different germline variants in the same gene causing either a neurodevelopmental disorder (NDD) or a tumor predisposition syndrome. Here, we report a series of 11 de novo germline heterozygous missense BAP1 variants associated with a rare syndromic NDD. Functional analysis showed that most of the variants cannot rescue the consequences of BAP1 inactivation, suggesting a loss-of-function mechanism. In T cells isolated from two affected children, H2A deubiquitination was impaired. In matching peripheral blood mononuclear cells, histone H3 K27 acetylation ChIP-seq indicated that these BAP1 variants induced genome-wide chromatin state alterations, with enrichment for regulatory regions surrounding genes of the ubiquitin-proteasome system (UPS). Altogether, these results define a clinical syndrome caused by rare germline missense BAP1 variants that alter chromatin remodeling through abnormal histone ubiquitination and lead to transcriptional dysregulation of developmental genes.

SMPD4 regulates mitotic nuclear envelope dynamics and its loss causes microcephaly and diabetes.

Biallelic loss-of-function variants in SMPD4 cause a rare and severe neurodevelopmental disorder with progressive congenital microcephaly and early death. SMPD4 encodes a sphingomyelinase that hydrolyses sphingomyelin into ceramide at neutral pH and can thereby affect membrane lipid homeostasis. SMPD4 localizes to the membranes of the endoplasmic reticulum and nuclear envelope and interacts with nuclear pore complexes (NPC). We refine the clinical phenotype of loss-of-function SMPD4 variants by describing five individuals from three unrelated families with longitudinal data due to prolonged survival. All individuals surviving beyond infancy developed insulin-dependent diabetes, besides presenting with a severe neurodevelopmental disorder and microcephaly, making diabetes one of the most frequent age-dependent non-cerebral abnormalities. We studied the function of SMPD4 at the cellular and organ levels. Knock-down of SMPD4 in human neural stem cells causes reduced proliferation rates and prolonged mitosis. Moreover, SMPD4 depletion results in abnormal nuclear envelope breakdown and reassembly during mitosis and decreased post-mitotic NPC insertion. Fibroblasts from affected individuals show deficient SMPD4-specific neutral sphingomyelinase activity, without changing (sub)cellular lipidome fractions, which suggests a local function of SMPD4 on the nuclear envelope. In embryonic mouse brain, knockdown of Smpd4 impairs cortical progenitor proliferation and induces premature differentiation by altering the balance between neurogenic and proliferative progenitor cell divisions. We hypothesize that, in individuals with SMPD4-related disease, nuclear envelope bending, which is needed to insert NPCs in the nuclear envelope, is impaired in the absence of SMPD4 and interferes with cerebral corticogenesis and survival of pancreatic beta cells.

CLEC16A interacts with retromer and TRIM27, and its loss impairs endosomal trafficking and neurodevelopment.

CLEC16A is a membrane-associated C-type lectin protein that functions as a E3-ubiquitin ligase. CLEC16A regulates autophagy and mitophagy, and reportedly localizes to late endosomes. GWAS studies have associated CLEC16A SNPs to various auto-immune and neurological disorders, including multiple sclerosis and Parkinson disease. Studies in mouse models imply a role for CLEC16A in neurodegeneration. We identified bi-allelic CLEC16A truncating variants in siblings from unrelated families presenting with a severe neurodevelopmental disorder including microcephaly, brain atrophy, corpus callosum dysgenesis, and growth retardation. To understand the function of CLEC16A in neurodevelopment we used in vitro models and zebrafish embryos. We observed CLEC16A localization to early endosomes in HEK293T cells. Mass spectrometry of human CLEC16A showed interaction with endosomal retromer complex subunits and the endosomal ubiquitin ligase TRIM27. Expression of the human variant leading to C-terminal truncated CLEC16A, abolishes both its endosomal localization and interaction with TRIM27, suggesting a loss-of-function effect. CLEC16A knockdown increased TRIM27 adhesion to early endosomes and abnormal accumulation of endosomal F-actin, a sign of disrupted vesicle sorting. Mutagenesis of clec16a by CRISPR-Cas9 in zebrafish embryos resulted in accumulated acidic/phagolysosome compartments, in neurons and microglia, and dysregulated mitophagy. The autophagocytic phenotype was rescued by wild-type human CLEC16A but not the C-terminal truncated CLEC16A. Our results demonstrate that CLEC16A closely interacts with retromer components and regulates endosomal fate by fine-tuning levels of TRIM27 and polymerized F-actin on the endosome surface. Dysregulation of CLEC16A-mediated endosomal sorting is associated with neurodegeneration, but it also causes accumulation of autophagosomes and unhealthy mitochondria during brain development.

De novo MCM6 variants in neurodevelopmental disorders: a recognizable phenotype related to zinc binding residues.

The minichromosome maintenance (MCM) complex acts as a DNA helicase during DNA replication, and thereby regulates cell cycle progression and proliferation. In addition, MCM-complex components localize to centrosomes and play an independent role in ciliogenesis. Pathogenic variants in genes coding for MCM components and other DNA replication factors have been linked to growth and developmental disorders as Meier-Gorlin syndrome and Seckel syndrome. Trio exome/genome sequencing identified the same de novo MCM6 missense variant p.(Cys158Tyr) in two unrelated individuals that presented with overlapping phenotypes consisting of intra-uterine growth retardation, short stature, congenital microcephaly, endocrine features, developmental delay and urogenital anomalies. The identified variant affects a zinc binding cysteine in the MCM6 zinc finger signature. This domain, and specifically cysteine residues, are essential for MCM-complex dimerization and the induction of helicase activity, suggesting a deleterious effect of this variant on DNA replication. Fibroblasts derived from the two affected individuals showed defects both in ciliogenesis and cell proliferation. We additionally traced three unrelated individuals with de novo MCM6 variants in the oligonucleotide binding (OB)-fold domain, presenting with variable (neuro)developmental features including autism spectrum disorder, developmental delay, and epilepsy. Taken together, our findings implicate de novo MCM6 variants in neurodevelopmental disorders. The clinical features and functional defects related to the zinc binding residue resemble those observed in syndromes related to other MCM components and DNA replication factors, while de novo OB-fold domain missense variants may be associated with more variable neurodevelopmental phenotypes. These data encourage consideration of MCM6 variants in the diagnostic arsenal of NDD.

TMX2 Is a Crucial Regulator of Cellular Redox State, and Its Dysfunction Causes Severe Brain Developmental Abnormalities.

The redox state of the neural progenitors regulates physiological processes such as neuronal differentiation and dendritic and axonal growth. The relevance of endoplasmic reticulum (ER)-associated oxidoreductases in these processes is largely unexplored. We describe a severe neurological disorder caused by bi-allelic loss-of-function variants in thioredoxin (TRX)-related transmembrane-2 (TMX2); these variants were detected by exome sequencing in 14 affected individuals from ten unrelated families presenting with congenital microcephaly, cortical polymicrogyria, and other migration disorders. TMX2 encodes one of the five TMX proteins of the protein disulfide isomerase family, hitherto not linked to human developmental brain disease. Our mechanistic studies on protein function show that TMX2 localizes to the ER mitochondria-associated membranes (MAMs), is involved in posttranslational modification and protein folding, and undergoes physical interaction with the MAM-associated and ER folding chaperone calnexin and ER calcium pump SERCA2. These interactions are functionally relevant because TMX2-deficient fibroblasts show decreased mitochondrial respiratory reserve capacity and compensatory increased glycolytic activity. Intriguingly, under basal conditions TMX2 occurs in both reduced and oxidized monomeric form, while it forms a stable dimer under treatment with hydrogen peroxide, recently recognized as a signaling molecule in neural morphogenesis and axonal pathfinding. Exogenous expression of the pathogenic TMX2 variants or of variants with an in vitro mutagenized TRX domain induces a constitutive TMX2 polymerization, mimicking an increased oxidative state. Altogether these data uncover TMX2 as a sensor in the MAM-regulated redox signaling pathway and identify it as a key adaptive regulator of neuronal proliferation, migration, and organization in the developing brain.

Loss of SMPD4 Causes a Developmental Disorder Characterized by Microcephaly and Congenital Arthrogryposis.

Sphingomyelinases generate ceramide from sphingomyelin as a second messenger in intracellular signaling pathways involved in cell proliferation, differentiation, or apoptosis. Children from 12 unrelated families presented with microcephaly, simplified gyral pattern of the cortex, hypomyelination, cerebellar hypoplasia, congenital arthrogryposis, and early fetal/postnatal demise. Genomic analysis revealed bi-allelic loss-of-function variants in SMPD4, coding for the neutral sphingomyelinase-3 (nSMase-3/SMPD4). Overexpression of human Myc-tagged SMPD4 showed localization both to the outer nuclear envelope and the ER and additionally revealed interactions with several nuclear pore complex proteins by proteomics analysis. Fibroblasts from affected individuals showed ER cisternae abnormalities, suspected for increased autophagy, and were more susceptible to apoptosis under stress conditions, while treatment with siSMPD4 caused delayed cell cycle progression. Our data show that SMPD4 links homeostasis of membrane sphingolipids to cell fate by regulating the cross-talk between the ER and the outer nuclear envelope, while its loss reveals a pathogenic mechanism in microcephaly.

Biallelic NDC1 variants that interfere with ALADIN binding are associated with neuropathy and triple A-like syndrome.

Nuclear pore complexes (NPCs) regulate nucleocytoplasmic transport and are anchored in the nuclear envelope by the transmembrane nucleoporin NDC1. NDC1 is essential for post-mitotic NPC assembly and the recruitment of ALADIN to the nuclear envelope. While no human disorder has been associated to one of the three transmembrane nucleoporins, biallelic variants in AAAS, encoding ALADIN, cause triple A syndrome (Allgrove syndrome). Triple A syndrome, characterized by alacrima, achalasia, and adrenal insufficiency, often includes progressive demyelinating polyneuropathy and other neurological complaints. In this report, diagnostic exome and/or RNA sequencing was performed in seven individuals from four unrelated consanguineous families with AAAS-negative triple A syndrome. Molecular and clinical studies followed to elucidate the pathogenic mechanism. The affected individuals presented with intellectual disability, motor impairment, severe demyelinating with secondary axonal polyneuropathy, alacrima, and achalasia. None of the affected individuals has adrenal insufficiency. All individuals presented with biallelic NDC1 in-frame deletions or missense variants that affect amino acids and protein domains required for ALADIN binding. No other significant variants associated with the phenotypic features were reported. Skin fibroblasts derived from affected individuals show decreased recruitment of ALADIN to the NE and decreased post-mitotic NPC insertion, confirming pathogenicity of the variants. Taken together, our results implicate biallelic NDC1 variants in the pathogenesis of polyneuropathy and a triple A-like disorder without adrenal insufficiency, by interfering with physiological NDC1 functions, including the recruitment of ALADIN to the NPC.

Biallelic DAB1 Variants Are Associated With Mild Lissencephaly and Cerebellar Hypoplasia.

OBJECTIVE: We aimed to identify pathogenic variants in a girl with epilepsy, developmental delay, cerebellar ataxia, oral motor difficulty, and structural brain abnormalities with the use of whole-exome sequencing. METHODS: Whole-exome trio analysis and molecular functional studies were performed in addition to the clinical findings and neuroimaging studies. RESULTS: Brain MRI showed mild pachygyria, hypoplasia of the cerebellar vermis, and abnormal foliation of the cerebellar vermis, suspected for a variant in one of the genes of the Reelin pathway. Trio whole-exome sequencing and additional functional studies were performed to identify the pathogenic variants. Trio whole-exome sequencing revealed compound heterozygous splice variants in DAB1, both affecting the highly conserved functional phosphotyrosine-binding domain. Expression studies in patient-derived cells showed loss of normal transcripts, confirming pathogenicity. CONCLUSIONS: We conclude that these variants are very likely causally related to the cerebral phenotype and propose to consider loss-of-function DAB1 variants in patients with RELN-like cortical malformations.

Reduction of oxidative stress suppresses poly-GR-mediated toxicity in zebrafish embryos.

The hexanucleotide (G4C2)-repeat expansion in the C9ORF72 gene is the most common pathogenic cause of frontotemporal dementia (FTD) and amyotrophic lateral sclerosis (ALS). This repeat expansion can be translated into dipeptide repeat proteins (DPRs), and distribution of the poly-GR DPR correlates with neurodegeneration in postmortem C9FTD/ALS brains. Here, we assessed poly-GR toxicity in zebrafish embryos, using an annexin A5-based fluorescent transgenic line (secA5) that allows for detection and quantification of apoptosis in vivo. Microinjection of RNA encoding poly-GR into fertilized oocytes evoked apoptosis in the brain and abnormal motor neuron morphology in the trunk of 1-4-days postfertilization embryos. Poly-GR can be specifically detected in protein homogenates from injected zebrafish and in the frontal cortexes of C9FTD/ALS cases. Poly-GR expression further elevated MitoSOX levels in zebrafish embryos, indicating oxidative stress. Inhibition of reactive oxygen species using Trolox showed full suppression of poly-GR toxicity. Our study indicates that poly-GR can exert its toxicity via oxidative stress. This zebrafish model can be used to find suppressors of poly-GR toxicity and identify its molecular targets underlying neurodegeneration observed in C9FTD/ALS.

Multidisciplinary interaction and MCD gene discovery. The perspective of the clinical geneticist.

The increasing pace of gene discovery in the last decade has brought a major change in the way the genetic causes of brain malformations are being diagnosed. Unbiased genomic screening has gained the first place in the diagnostic protocol of a child with congenital (brain) anomalies and the detected variants are matched with the phenotypic presentation afterwards. This process is defined as "reverse phenotyping". Screening of DNA, through copy number variant analysis of microarrays and analysis of exome data on different platforms, obtained from the index patient and both parents has become a routine approach in many centers worldwide. Clinicians are used to multidisciplinary team interaction in patient care and disease management and this explains why the majority of research that has led to the discovery of new genetic disorders nowadays proceeds from clinical observations to genomic analysis and to data exchange facilitated by open access sharing databases. However, the relevance of multidisciplinary team interaction has not been object of systematic research in the field of brain malformations. This review will illustrate some examples of how diagnostically driven questions through multidisciplinary interaction, among clinical and preclinical disciplines, can be successful in the discovery of new genes related to brain malformations. The first example illustrates the setting of interaction among neurologists, geneticists and neuro-radiologists. The second illustrates the importance of interaction among clinical dysmorphologists for pattern recognition of syndromes with multiple congenital anomalies. The third example shows how fruitful it can be to step out of the "clinical comfort zone", and interact with basic scientists in applying emerging technologies to solve the diagnostic puzzles.