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.

High-yield identification of pathogenic NF1 variants by skin fibroblast transcriptome screening after apparently normal diagnostic DNA testing.

Neurofibromatosis type 1 (NF1) is caused by inactivating mutations in NF1. Due to the size, complexity, and high mutation rate at the NF1 locus, the identification of causative variants can be challenging. To obtain a molecular diagnosis in 15 individuals meeting diagnostic criteria for NF1, we performed transcriptome analysis (RNA-seq) on RNA obtained from cultured skin fibroblasts. In each case, routine molecular DNA diagnostics had failed to identify a disease-causing variant in NF1. A pathogenic variant or abnormal mRNA splicing was identified in 13 cases: 6 deep intronic variants and 2 transposon insertions causing noncanonical splicing, 3 postzygotic changes, 1 branch point mutation and, in 1 case, abnormal splicing for which the responsible DNA change remains to be identified. These findings helped resolve the molecular findings for an additional 17 individuals in multiple families with NF1, demonstrating the utility of skin-fibroblast-based transcriptome analysis for molecular diagnostics. RNA-seq improves mutation detection in NF1 and provides a powerful complementary approach to DNA-based methods. Importantly, our approach is applicable to other genetic disorders, particularly those caused by a wide variety of variants in a limited number of genes and specifically for individuals in whom routine molecular DNA diagnostics did not identify the causative variant.

Defining the phenotypical spectrum associated with variants in TUBB2A.

BACKGROUND: Variants in genes belonging to the tubulin superfamily account for a heterogeneous spectrum of brain malformations referred to as tubulinopathies. Variants in TUBB2A have been reported in 10 patients with a broad spectrum of brain imaging features, ranging from a normal cortex to polymicrogyria, while one patient has been reported with progressive atrophy of the cerebellar vermis. METHODS: In order to further refine the phenotypical spectrum associated with TUBB2A, clinical and imaging features of 12 patients with pathogenic TUBB2A variants, recruited via the international network of the authors, were reviewed. RESULTS: We report 12 patients with eight novel and one recurrent variants spread throughout the TUBB2A gene but encoding for amino acids clustering at the protein surface. Eleven patients (91.7%) developed seizures in early life. All patients suffered from intellectual disability, and 11 patients had severe motor developmental delay, with 4 patients (36.4 %) being non-ambulatory. The cerebral cortex was normal in five individuals and showed dysgyria of variable severity in seven patients. Associated brain malformations were less frequent in TUBB2A patients compared with other tubulinopathies. None of the patients had progressive cerebellar atrophy. CONCLUSION: The imaging phenotype associated with pathogenic variants in TUBB2A is highly variable, ranging from a normal cortex to extensive dysgyria with associated brain malformations. For recurrent variants, no clear genotype-phenotype correlations could be established, suggesting the role of additional modifiers.

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.

Heterogeneous clinical phenotypes and cerebral malformations reflected by rotatin cellular dynamics.

Recessive mutations in RTTN, encoding the protein rotatin, were originally identified as cause of polymicrogyria, a cortical malformation. With time, a wide variety of other brain malformations has been ascribed to RTTN mutations, including primary microcephaly. Rotatin is a centrosomal protein possibly involved in centriolar elongation and ciliogenesis. However, the function of rotatin in brain development is largely unknown and the molecular disease mechanism underlying cortical malformations has not yet been elucidated. We performed both clinical and cell biological studies, aimed at clarifying rotatin function and pathogenesis. Review of the 23 published and five unpublished clinical cases and genomic mutations, including the effect of novel deep intronic pathogenic mutations on RTTN transcripts, allowed us to extrapolate the core phenotype, consisting of intellectual disability, short stature, microcephaly, lissencephaly, periventricular heterotopia, polymicrogyria and other malformations. We show that the severity of the phenotype is related to residual function of the protein, not only the level of mRNA expression. Skin fibroblasts from eight affected individuals were studied by high resolution immunomicroscopy and flow cytometry, in parallel with in vitro expression of RTTN in HEK293T cells. We demonstrate that rotatin regulates different phases of the cell cycle and is mislocalized in affected individuals. Mutant cells showed consistent and severe mitotic failure with centrosome amplification and multipolar spindle formation, leading to aneuploidy and apoptosis, which could relate to depletion of neuronal progenitors often observed in microcephaly. We confirmed the role of rotatin in functional and structural maintenance of primary cilia and determined that the protein localized not only to the basal body, but also to the axoneme, proving the functional interconnectivity between ciliogenesis and cell cycle progression. Proteomics analysis of both native and exogenous rotatin uncovered that rotatin interacts with the neuronal (non-muscle) myosin heavy chain subunits, motors of nucleokinesis during neuronal migration, and in human induced pluripotent stem cell-derived bipolar mature neurons rotatin localizes at the centrosome in the leading edge. This illustrates the role of rotatin in neuronal migration. These different functions of rotatin explain why RTTN mutations can lead to heterogeneous cerebral malformations, both related to proliferation and migration defects.

Basic Fibroblast Growth Factor Induces Adipogenesis in Orbital Fibroblasts: Implications for the Pathogenesis of Graves' Orbitopathy.

BACKGROUND: Basic fibroblast growth factor (bFGF) has been implicated in the pathogenesis of Graves' orbitopathy (GO). It stimulates several processes, including hyaluronan synthesis, involved in orbital tissue volume expansion and may act synergistically with platelet-derived growth factor (PDGF)-BB. PDGF-BB is known to stimulate adipogenesis in orbital fibroblasts, but the effect of bFGF on adipogenesis in orbital fibroblasts is so far unknown. This study was conducted to determine whether (i) bFGF induces adipogenesis in orbital fibroblasts, (ii) bFGF and PDGF-BB together exert an additive or synergistic effect on adipogenesis, and (iii) treatment directed at bFGF- and PDGF-BB signaling may potentially be of interest for the treatment of GO. METHODS: Orbital fibroblasts from GO patients and controls were cultured in adipocyte differentiation medium with or without bFGF and/or PDGF-BB at different concentrations. Adipogenesis was determined by Oil Red O staining and messenger RNA expression of the late adipocyte differentiation markers cell death-inducing DFFA-like effector C (CIDEC) and adiponectin (ADIPOQ). To demonstrate involvement of FGF-receptor and PDGF-receptor signaling, experiments were also conducted in the presence of dasatinib (inhibitor of PDGF-receptor) or nintedanib (inhibitor of PDGF-receptor and FGF-receptor). RESULTS: bFGF significantly stimulated adipogenesis by orbital fibroblasts, as shown by increased Oil Red O staining and CIDEC and ADIPOQ expression after 14 days of differentiation. Furthermore, an additive effect of bFGF/PDGF-BB co-stimulation on adipogenesis was observed at the lowest concentration (12.5 ng/mL) of the growth factors tested. Nintedanib completely inhibited bFGF-, PDGF-BB-, and bFGF/PDGF-BB-induced adipogenesis, while dasatinib only fully abrogated PDGF-BB-induced adipogenesis. CONCLUSION: bFGF induces adipogenesis in orbital fibroblasts and as such may contribute to GO. The additive effect of bFGF and PDGF-BB on adipogenesis, along with the observed inhibitory effects of dasatinib and nintedanib, point at independent receptor-mediated effects. This supports the hypothesis that multi-target directed therapy might be more efficient in the treatment of GO.

MACF1 Mutations Encoding Highly Conserved Zinc-Binding Residues of the GAR Domain Cause Defects in Neuronal Migration and Axon Guidance.

To date, mutations in 15 actin- or microtubule-associated genes have been associated with the cortical malformation lissencephaly and variable brainstem hypoplasia. During a multicenter review, we recognized a rare lissencephaly variant with a complex brainstem malformation in three unrelated children. We searched our large brain-malformation databases and found another five children with this malformation (as well as one with a less severe variant), analyzed available whole-exome or -genome sequencing data, and tested ciliogenesis in two affected individuals. The brain malformation comprised posterior predominant lissencephaly and midline crossing defects consisting of absent anterior commissure and a striking W-shaped brainstem malformation caused by small or absent pontine crossing fibers. We discovered heterozygous de novo missense variants or an in-frame deletion involving highly conserved zinc-binding residues within the GAR domain of MACF1 in the first eight subjects. We studied cilium formation and found a higher proportion of mutant cells with short cilia than of control cells with short cilia. A ninth child had similar lissencephaly but only subtle brainstem dysplasia associated with a heterozygous de novo missense variant in the spectrin repeat domain of MACF1. Thus, we report variants of the microtubule-binding GAR domain of MACF1 as the cause of a distinctive and most likely pathognomonic brain malformation. A gain-of-function or dominant-negative mechanism appears likely given that many heterozygous mutations leading to protein truncation are included in the ExAC Browser. However, three de novo variants in MACF1 have been observed in large schizophrenia cohorts.