Search

Tuberous sclerosis complex (TSC) is an inherited neurodevelopmental disorder (NDD) with frequent manifestations of epilepsy and autism spectrum disorder (ASD). TSC is caused by inactivating mutations in TSC1 or TSC2 tumor suppressor genes, with encoded proteins hamartin (TSC1) and tuberin (TSC2) forming a functional complex inhibiting mechanistic target of rapamycin complex 1 (mTORC1) signaling. This has led to treatment with allosteric mTORC1 inhibitor rapamycin analogs ("rapalogs") for TSC tumors; however, rapalogs are ineffective for treating neurodevelopmental manifestations. mTORC1 signaling controls protein synthesis by regulating formation of the eIF4F complex, with further modulation by MNK1/2 kinases via phosphorylation of the eIF4F subunit eIF4E. While both these pathways modulate translation, comparing their impact on transcriptome-wide mRNA translation, as well as effects of inhibiting these pathways in TSC has not been explored. Here, employing CRISPR-modified, isogenic TSC2 patient-derived neural progenitor cells (NPCs), we have examined transcriptome-wide changes in mRNA translation upon TSC2 loss. Our results reveal dysregulated translation in TSC2 -Null NPCs, which significantly overlaps with the translatome from TSC1 -Null NPCs. Interestingly, numerous non-monogenic ASD-, NDD-and epilepsy-associated genes identified in patients harboring putative loss-of-function mutations, were translationally suppressed in TSC2 -Null NPCs. Importantly, translation of these ASD- and NDD-associated genes was reversed upon inhibition of either mTORC1 or MNK1/2 signaling using RMC-6272 or eFT-508, respectively. This study establishes the importance of mTORC1-eIF4F- and MNK-eIF4E-sensitive mRNA translation in TSC, ASD and other neurodevelopmental disorders laying the groundwork for evaluating drugs in clinical development that target these pathways as a treatment strategy for these disorders.

Familial dysautonomia (FD) is a rare recessive neurodevelopmental disease caused by a splice mutation in the Elongator acetyltransferase complex subunit 1 (ELP1) gene. This mutation results in a tissue-specific reduction of ELP1 protein, with the lowest levels in the central and peripheral nervous systems (CNS and PNS, respectively). FD patients exhibit complex neurological phenotypes due to the loss of sensory and autonomic neurons. Disease symptoms include decreased pain and temperature perception, impaired or absent myotatic reflexes, proprioceptive ataxia, and progressive retinal degeneration. While the involvement of the PNS in FD pathogenesis has been clearly recognized, the underlying mechanisms responsible for the preferential neuronal loss remain unknown. In this study, we aimed to elucidate the molecular mechanisms underlying FD by conducting a comprehensive transcriptome analysis of neuronal tissues from the phenotypic mouse model TgFD9; Elp1 Δ20/flox . This mouse recapitulates the same tissue-specific ELP1 mis-splicing observed in patients while modeling many of the disease manifestations. Comparison of FD and control transcriptomes from dorsal root ganglion (DRG), trigeminal ganglion (TG), medulla (MED), cortex, and spinal cord (SC) showed significantly more differentially expressed genes (DEGs) in the PNS than the CNS. We then identified genes that were tightly co-expressed and functionally dependent on the level of full-length ELP1 transcript. These genes, defined as ELP1 dose-responsive genes, were combined with the DEGs to generate tissue-specific dysregulated FD signature genes and networks. Within the PNS networks, we observed direct connections between Elp1 and genes involved in tRNA synthesis and genes related to amine metabolism and synaptic signaling. Importantly, transcriptomic dysregulation in PNS tissues exhibited enrichment for neuronal subtype markers associated with peptidergic nociceptors and myelinated sensory neurons, which are known to be affected in FD. In summary, this study has identified critical tissue-specific gene networks underlying the etiology of FD and provides new insights into the molecular basis of the disease., (© 2024. The Author(s).)

Background: Neurofibromatosis 2 (NF2) is an inherited disorder caused by bi-allelic inactivation of the NF2 tumor suppressor gene. NF2-associated tumors, including schwannoma and meningioma, are resistant to chemotherapy, often recurring despite surgery and/or radiation, and have generally shown cytostatic response to signal transduction pathway inhibitors, highlighting the need for improved cytotoxic therapies., Methods: Leveraging data from our previous high-throughput drug screening in NF2 preclinical models, we identified a class of compounds targeting the ubiquitin-proteasome pathway (UPP), and undertook studies using candidate UPP inhibitors, ixazomib/MLN9708, pevonedistat/MLN4924, and TAK-243/MLN7243. Employing human primary and immortalized meningioma (MN) cell lines, CRISPR-modified Schwann cells (SCs), and mouse Nf2-/- SCs, we performed dose response testing, flow cytometry-based Annexin V and cell cycle analyses, and RNA-sequencing to identify potential underlying mechanisms of apoptosis. In vivo efficacy was also assessed in orthotopic NF2-deficient meningioma and schwannoma tumor models., Results: Testing of three UPP inhibitors demonstrated potent reduction in cell viability and induction of apoptosis for ixazomib or TAK-243, but not pevonedistat. In vitro analyses revealed that ixazomib or TAK-243 downregulates expression of c-KIT and PDGFRα, as well as the E3 ubiquitin ligase SKP2 while upregulating genes associated with endoplasmic reticulum stress-mediated activation of the unfolded protein response (UPR). In vivo treatment of mouse models revealed delayed tumor growth, suggesting a therapeutic potential., Conclusions: This study demonstrates the efficacy of proteasomal pathway inhibitors in meningioma and schwannoma preclinical models and lays the groundwork for use of these drugs as a promising novel treatment strategy for NF2 patients., (© The Author(s) 2023. Published by Oxford University Press on behalf of the Society for Neuro-Oncology. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Autism spectrum disorder (ASD) is a heritable neurodevelopmental disorder characterized by deficits in social interactions and communication. Protein-altering variants in many genes have been shown to contribute to ASD; however, understanding the convergence across many genes remains a challenge. We demonstrate that coexpression patterns from 993 human postmortem brains are significantly correlated with the transcriptional consequences of CRISPR perturbations in human neurons. Across 71 ASD risk genes, there was significant tissue-specific convergence implicating synaptic pathways. Tissue-specific convergence was further demonstrated across schizophrenia and atrial fibrillation risk genes. The degree of ASD convergence was significantly correlated with ASD association from rare variation and differential expression in ASD brains. Positively convergent genes showed intolerance to functional mutations and had shorter coding lengths than known risk genes even after removing association with ASD. These results indicate that convergent coexpression can identify potentially novel genes that are unlikely to be discovered by sequencing studies., Competing Interests: The authors declare no competing interests., (© 2023 The Author(s).)

Disruptive mutations in the chromodomain helicase DNA-binding protein 8 gene (CHD8) have been recurrently associated with autism spectrum disorders (ASDs). Here we investigated how chromatin reacts to CHD8 suppression by analyzing a panel of histone modifications in induced pluripotent stem cell-derived neural progenitors. CHD8 suppression led to significant reduction (47.82%) in histone H3K36me3 peaks at gene bodies, particularly impacting on transcriptional elongation chromatin states. H3K36me3 reduction specifically affects highly expressed, CHD8-bound genes and correlates with altered alternative splicing patterns of 462 genes implicated in 'regulation of RNA splicing' and 'mRNA catabolic process'. Mass spectrometry analysis uncovered a novel interaction between CHD8 and the splicing regulator heterogeneous nuclear ribonucleoprotein L (hnRNPL), providing the first mechanistic insights to explain the CHD8 suppression-derived splicing phenotype, partly implicating SETD2, a H3K36me3 methyltransferase. In summary, our results point toward broad molecular consequences of CHD8 suppression, entailing altered histone deposition/maintenance and RNA processing regulation as important regulatory processes in ASD., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Point mutations and structural variants that directly disrupt the coding sequence of MEF2C have been associated with a spectrum of neurodevelopmental disorders (NDDs). However, the impact of MEF2C haploinsufficiency on neurodevelopmental pathways and synaptic processes is not well understood, nor are the complex mechanisms that govern its regulation. To explore the functional changes associated with structural variants that alter MEF2C expression and/or regulation, we generated an allelic series of 204 isogenic human induced pluripotent stem cell (hiPSC)-derived neural stem cells and glutamatergic induced neurons. These neuronal models harbored CRISPR-engineered mutations that involved direct deletion of MEF2C or deletion of the boundary points for topologically associating domains (TADs) and chromatin loops encompassing MEF2C. Systematic profiling of mutation-specific alterations, contrasted to unedited controls that were exposed to the same guide RNAs for each edit, revealed that deletion of MEF2C caused differential expression of genes associated with neurodevelopmental pathways and synaptic function. We also discovered significant reduction in synaptic activity measured by multielectrode arrays (MEAs) in neuronal cells. By contrast, we observed robust buffering against MEF2C regulatory disruption following deletion of a distal 5q14.3 TAD and loop boundary, whereas homozygous loss of a proximal loop boundary resulted in down-regulation of MEF2C expression and reduced electrophysiological activity on MEA that was comparable to direct gene disruption. Collectively, these studies highlight the considerable functional impact of MEF2C deletion in neuronal cells and systematically characterize the complex interactions that challenge a priori predictions of regulatory consequences from structural variants that disrupt three-dimensional genome organization., Competing Interests: Declaration of interests M.E.T. receives research funding and/or reagents from Levo Therapeutics, Microsoft Inc, and Illumina Inc. K.M. initiated employment with Tornado Bio during the writing of this manuscript and is at present an employee of the company. J.F.G. is a founder and member of the scientific advisory board of Triplet Therapeutics, Inc., and has been a paid consultant to Biogen, Inc., Pfizer, Inc., and Wave Biosciences, Inc., (Copyright © 2022 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

The canonical paradigm for converting genetic association to mechanism involves iteratively mapping individual associations to the proximal genes through which they act. In contrast, in the present study we demonstrate the feasibility of extracting biological insights from a very large region of the genome and leverage this strategy to study the genetic influences on autism. Using a new statistical approach, we identified the 33-Mb p-arm of chromosome 16 (16p) as harboring the greatest excess of autism's common polygenic influences. The region also includes the mechanistically cryptic and autism-associated 16p11.2 copy number variant. Analysis of RNA-sequencing data revealed that both the common polygenic influences within 16p and the 16p11.2 deletion were associated with decreased average gene expression across 16p. The transcriptional effects of the rare deletion and diffuse common variation were correlated at the level of individual genes and analysis of Hi-C data revealed patterns of chromatin contact that may explain this transcriptional convergence. These results reflect a new approach for extracting biological insight from genetic association data and suggest convergence of common and rare genetic influences on autism at 16p., (© 2022. The Author(s).)

Chromosome 16p11.2 reciprocal genomic disorder, resulting from recurrent copy-number variants (CNVs), involves intellectual disability, autism spectrum disorder (ASD), and schizophrenia, but the responsible mechanisms are not known. To systemically dissect molecular effects, we performed transcriptome profiling of 350 libraries from six tissues (cortex, cerebellum, striatum, liver, brown fat, and white fat) in mouse models harboring CNVs of the syntenic 7qF3 region, as well as cellular, transcriptional, and single-cell analyses in 54 isogenic neural stem cell, induced neuron, and cerebral organoid models of CRISPR-engineered 16p11.2 CNVs. Transcriptome-wide differentially expressed genes were largely tissue-, cell-type-, and dosage-specific, although more effects were shared between deletion and duplication and across tissue than expected by chance. The broadest effects were observed in the cerebellum (2,163 differentially expressed genes), and the greatest enrichments were associated with synaptic pathways in mouse cerebellum and human induced neurons. Pathway and co-expression analyses identified energy and RNA metabolism as shared processes and enrichment for ASD-associated, loss-of-function constraint, and fragile X messenger ribonucleoprotein target gene sets. Intriguingly, reciprocal 16p11.2 dosage changes resulted in consistent decrements in neurite and electrophysiological features, and single-cell profiling of organoids showed reciprocal alterations to the proportions of excitatory and inhibitory GABAergic neurons. Changes both in neuronal ratios and in gene expression in our organoid analyses point most directly to calretinin GABAergic inhibitory neurons and the excitatory/inhibitory balance as targets of disruption that might contribute to changes in neurodevelopmental and cognitive function in 16p11.2 carriers. Collectively, our data indicate the genomic disorder involves disruption of multiple contributing biological processes and that this disruption has relative impacts that are context specific., Competing Interests: Declaration of interests J.F.G. is a founder and member of the scientific advisory board of Triplet Therapeutics and has been a paid consultant to Biogen, Pfizer, and Wave Biosciences. M.E.T. receives research funding and/or reagents from Illumina, Levo Therapeutics, and Microsoft. E.M. receives research funding and reagents from PTC Therapeutics. The following co-authors are currently employed by for-profit companies or non-profit organizations: P.R. is founder and CEO of Kumuda; K.M. is currently employed by Tornado Bio; T.A. is a founder and CEO of Independent Data Lab and OmicsChart; and A.H. is currently employed by Bionano Genomics., (Copyright © 2022 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Familial dysautonomia (FD), a hereditary sensory and autonomic neuropathy, is caused by a mutation in the Elongator complex protein 1 (ELP1) gene that leads to a tissue-specific reduction of ELP1 protein. Our work to generate a phenotypic mouse model for FD headed to the discovery that homozygous deletion of the mouse Elp1 gene leads to embryonic lethality prior to mid-gestation. Given that FD is caused by a reduction, not loss, of ELP1, we generated two new mouse models by introducing different copy numbers of the human FD ELP1 transgene into the Elp1 knockout mouse (Elp1 -/- ) and observed that human ELP1 expression rescues embryonic development in a dose-dependent manner. We then conducted a comprehensive transcriptome analysis in mouse embryos to identify genes and pathways whose expression correlates with the amount of ELP1. We found that ELP1 is essential for the expression of genes responsible for nervous system development. Further, gene length analysis of the differentially expressed genes showed that the loss of Elp1 mainly impacts the expression of long genes and that by gradually restoring Elongator, their expression is progressively rescued. Finally, through evaluation of co-expression modules, we identified gene sets with unique expression patterns that depended on ELP1 expression., Competing Interests: Conflict of interest The authors declare competing financial interests. Funding: Research support from PTC Therapeutics, Inc. (S.A.S.). Personal financial interests: Susan A. Slaugenhaupt is a paid consultant to PTC Therapeutics and is an inventor on several US and foreign patents and patent applications assigned to the Massachusetts General Hospital, including U.S Patents 8,729,025 and 9,265,766, both entitled ‘Methods for altering mRNA splicing and treating familial dysautonomia by administering benzyladenine’, filed on August 31, 2012, and May 19, 2014, and related to use of kinetin; and U.S. Patent 10,675,475 entitled, ‘Compounds for improving mRNA splicing’ filed on July 14, 2017, and related to use of BPN-15477. Elisabetta Morini, Dadi Gao, Michael E. Talkowski and Susan A. Slaugenhaupt are inventors on an International Patent Application Number PCT/US2021/012103, assigned to Massachusetts General Hospital and entitled “RNA Splicing Modulation” related to the use of BPN-15477 in modulating splicing., (Copyright © 2021 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Cerebral organoids can be used to gain insights into cell type specific processes perturbed by genetic variants associated with neuropsychiatric disorders. However, robust and scalable phenotyping of organoids remains challenging. Here, we perform RNA sequencing on 71 samples comprising 1,420 cerebral organoids from 25 donors, and describe a framework (Orgo-Seq) to integrate bulk RNA and single-cell RNA sequence data. We apply Orgo-Seq to 16p11.2 deletions and 15q11-13 duplications, two loci associated with autism spectrum disorder, to identify immature neurons and intermediate progenitor cells as critical cell types for 16p11.2 deletions. We further applied Orgo-Seq to identify cell type-specific driver genes. Our work presents a quantitative phenotyping framework to integrate multi-transcriptomic datasets for the identification of cell types and cell type-specific co-expressed driver genes associated with neuropsychiatric disorders., (© 2022. The Author(s).)

Dystonia is a neurologic disorder associated with an increasingly large number of genetic variants in many genes, resulting in characteristic disturbances in volitional movement. Dissecting the relationships between these mutations and their functional outcomes is critical in understanding the pathways that drive dystonia pathogenesis. Here we established a pipeline for characterizing an allelic series of dystonia-specific mutations. We used this strategy to investigate the molecular consequences of genetic variation in THAP1, which encodes a transcription factor linked to neural differentiation. Multiple pathogenic mutations associated with dystonia cluster within distinct THAP1 functional domains and are predicted to alter DNA-binding properties and/or protein interactions differently, yet the relative impact of these varied changes on molecular signatures and neural deficits is unclear. To determine the effects of these mutations on THAP1 transcriptional activity, we engineered an allelic series of eight alterations in a common induced pluripotent stem cell background and differentiated these lines into a panel of near-isogenic neural stem cells (n = 94 lines). Transcriptome profiling followed by joint analysis of the most robust signatures across mutations identified a convergent pattern of dysregulated genes functionally related to neurodevelopment, lysosomal lipid metabolism, and myelin. On the basis of these observations, we examined mice bearing Thap1-disruptive alleles and detected significant changes in myelin gene expression and reduction of myelin structural integrity relative to control mice. These results suggest that deficits in neurodevelopment and myelination are common consequences of dystonia-associated THAP1 mutations and highlight the potential role of neuron-glial interactions in the pathogenesis of dystonia., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

Gonadotropin-releasing hormone (GnRH) neurons in the hypothalamus play a key role in the regulation of reproductive function. In this study, we sought an efficient method for generating GnRH neurons from human embryonic and induced pluripotent stem cells (hESC and hiPSC, respectively). First, we found that exposure of primitive neuroepithelial cells, rather than neuroprogenitor cells, to fibroblast growth factor 8 (FGF8), was more effective in generating GnRH neurons. Second, addition of kisspeptin to FGF8 further increased the efficiency rates of GnRH neurogeneration. Third, we generated a fluorescent marker mCherry labeled human embryonic GnRH cell line (mCh-hESC) using a CRISPR-Cas9 targeting approach. Fourth, we examined physiological characteristics of GnRH (mCh-hESC) neurons: similar to GnRH neurons in vivo, they released the GnRH peptide in a pulsatile manner at ~60 min intervals; GnRH release increased in response to high potassium, kisspeptin, estradiol, and neurokinin B challenges; and injection of depolarizing current induced action potentials. Finally, we characterized developmental changes in transcriptomes of GnRH neurons using hESC, hiPSC, and mCh-hESC. The developmental pattern of transcriptomes was remarkably similar among the 3 cell lines. Collectively, human stem cell-derived GnRH neurons will be an important tool for establishing disease models to understand diseases, such as idiopathic hypothalamic hypogonadism, and testing contraceptive drugs., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Endocrine Society.)

Neurofibromatosis Type 2 (NF2) is an autosomal dominant genetic syndrome caused by mutations in the NF2 tumor suppressor gene resulting in multiple schwannomas and meningiomas. There are no FDA approved therapies for these tumors and their relentless progression results in high rates of morbidity and mortality. Through a combination of high throughput screens, preclinical in vivo modeling, and evaluation of the kinome en masse, we identified actionable drug targets and efficacious experimental therapeutics for the treatment of NF2 related schwannomas and meningiomas. These efforts identified brigatinib (ALUNBRIG®), an FDA-approved inhibitor of multiple tyrosine kinases including ALK, to be a potent inhibitor of tumor growth in established NF2 deficient xenograft meningiomas and a genetically engineered murine model of spontaneous NF2 schwannomas. Surprisingly, neither meningioma nor schwannoma cells express ALK. Instead, we demonstrate that brigatinib inhibited multiple tyrosine kinases, including EphA2, Fer and focal adhesion kinase 1 (FAK1). These data demonstrate the power of the de novo unbiased approach for drug discovery and represents a major step forward in the advancement of therapeutics for the treatment of NF2 related malignancies., Competing Interests: No authors have competing interests.

Pre-mRNA splicing is a key controller of human gene expression. Disturbances in splicing due to mutation lead to dysregulated protein expression and contribute to a substantial fraction of human disease. Several classes of splicing modulator compounds (SMCs) have been recently identified and establish that pre-mRNA splicing represents a target for therapy. We describe herein the identification of BPN-15477, a SMC that restores correct splicing of ELP1 exon 20. Using transcriptome sequencing from treated fibroblast cells and a machine learning approach, we identify BPN-15477 responsive sequence signatures. We then leverage this model to discover 155 human disease genes harboring ClinVar mutations predicted to alter pre-mRNA splicing as targets for BPN-15477. Splicing assays confirm successful correction of splicing defects caused by mutations in CFTR, LIPA, MLH1 and MAPT. Subsequent validations in two disease-relevant cellular models demonstrate that BPN-15477 increases functional protein, confirming the clinical potential of our predictions.

Background: Ubiquitination has a central role in numerous biological processes, including cell development, stress responses and ageing. Perturbed ubiquitination has been implicated in human diseases ranging from cancer to neurodegenerative diseases. SIAH1 encodes a RING-type E3 ubiquitin ligase involved in protein ubiquitination. Among numerous other roles, SIAH1 regulates metabotropic glutamate receptor signalling and affects neural cell fate. Moreover, SIAH1 positively regulates Wnt signalling through ubiquitin-mediated degradation of Axin and accumulation of β-catenin., Methods: Trio exome sequencing followed by Sanger validation was undertaken in five individuals with syndromic developmental delay. Three-dimensional structural modelling was used to predict pathogenicity of affected residues. Wnt stimulatory activity was measured by luciferase reporter assays and Axin degradation assays in HEK293 cells transfected with wild-type and mutant SIAH1 expression plasmids., Results: We report five unrelated individuals with shared features of developmental delay, infantile hypotonia, dysmorphic features and laryngomalacia, in whom exome sequencing identified de novo monoallelic variants in SIAH1 . In silico protein modelling suggested alteration of conserved functional sites. In vitro experiments demonstrated loss of Wnt stimulatory activity with the SIAH1 mutants, suggesting variant pathogenicity., Conclusion: Our results lend support to SIAH1 as a candidate Mendelian disease gene for a recognisable syndrome, further strengthening the connection between SIAH1 and neurodevelopmental disorders. Furthermore, the results suggest that dysregulation of the Wnt/β-catenin pathway may be involved in the pathogenesis., Competing Interests: Competing interests: TBP is an employee of GeneDx., (© Author(s) (or their employer(s)) 2021. No commercial re-use. See rights and permissions. Published by BMJ.)

Meningiomas (MNs), arising from the arachnoid/meningeal layer, are nonresponsive to chemotherapies, with ∼50% showing loss of the Neurofibromatosis 2 (NF2) tumor suppressor gene. Previously, we established NF2 loss activates mechanistic target of rapamycin complex 1 (mTORC1) and mechanistic target of rapamycin complex 2 (mTORC2) signaling, leading to clinical trials for NF2 and MN. Recently our omics studies identified activated ephrin (EPH) receptor and Src family kinases upon NF2 loss. Here, we report increased expression of several ligands in NF2-null human arachnoidal cells (ACs) and the MN cell line Ben-Men-1, particularly neuregulin-1/heregulin (NRG1), and confirm increased NRG1 secretion and activation of V-ERB-B avian erythroblastic leukemia viral oncogene homolog 3 (ERBB3) receptor kinase. Conditioned-medium from NF2-null ACs or exogenous NRG1 stimulated ERBB3, EPHA2, and mTORC1/2 signaling, suggesting pathway crosstalk. NF2-null cells treated with an ERBB3-neutralizing antibody partially downregulated mTOR pathway activation but showed no effect on viability. mTORC1/2 inhibitor treatment decreased NRG1 expression and downregulated ERBB3 while re-activating pAkt T308, suggesting a mechanism independent of NRG1-ERBB3 but likely involving activation of another upstream receptor kinase. Transcriptomics after mTORC1/2 inhibition confirmed decreased ERBB3/ERBB4 while revealing increased expression of insulin-like growth factor receptor 1 (IGF1R). Drug treatment co-targeting mTORC1/2 and IGF1R/insulin receptor attenuated pAkt T308 and showed synergistic effects on viability. Our findings indicate potential autocrine signaling where NF2 loss leads to secretion/activation of NRG1-ERBB3 signaling. mTORC1/2 inhibition downregulates NRG1-ERBB3, while upregulating pAkt T308 through an adaptive response involving IGF1R/insulin receptor and co-targeting these pathways may prove effective for treatment of NF2-deficient MN., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

Over the last two years, remarkable gene discovery efforts have implicated disruption of pathways involving gene regulatory functions and neuronal processes in autism spectrum disorder (ASD), and more broadly defined neurodevelopmental disorders (NDDs). Functional studies in the developing brain and across cell types demonstrate that the spatiotemporal expression patterns of many of these genes coalesce on subnetworks with distinct developmental trajectories. Here, we review the convergent biological processes derived from gene discovery and functional genomics in ASD and NDD from 2018-2020. We further probe the mechanistic insights that suggest these frequently perturbed pathways are interconnected and, ultimately, converge on specific functional deficits in human neurodevelopment., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Somatic expansion of the Huntington's disease (HD) CAG repeat drives the rate of a pathogenic process ultimately resulting in neuronal cell death. Although mechanisms of toxicity are poorly delineated, transcriptional dysregulation is a likely contributor. To identify modifiers that act at the level of CAG expansion and/or downstream pathogenic processes, we tested the impact of genetic knockout, in Htt Q111 mice, of Hdac2 or Hdac3 in medium-spiny striatal neurons that exhibit extensive CAG expansion and exquisite disease vulnerability. Both knockouts moderately attenuated CAG expansion, with Hdac2 knockout decreasing nuclear huntingtin pathology. Hdac2 knockout resulted in a substantial transcriptional response that included modification of transcriptional dysregulation elicited by the Htt Q111 allele, likely via mechanisms unrelated to instability suppression. Our results identify novel modifiers of different aspects of HD pathogenesis in medium-spiny neurons and highlight a complex relationship between the expanded Htt allele and Hdac2 with implications for targeting transcriptional dysregulation in HD., Competing Interests: MK, SE, MA, JS, MS, LH, JN, AS, RM, JW, MT No competing interests declared, DF DMF is a member of the scientific advisory board of Psy Therapeutics. SH SJH is a member of the scientific advisory board of Psy Therapeutics, Frequency Therapeutics and Souvien Therapeutics, and former member of the scientific advisory board of Rodin Therapeutics that is focused on HDAC2 inhibitors, none of whom were involved in the present study. SJH has also received speaking or consulting fees from Amgen, AstraZeneca, Biogen, Merck, Regenacy Pharmaceuticals, as well as sponsored research or gift funding from AstraZeneca, JW Pharmaceuticals, and Vesigen unrelated to the content of this manuscript. His financial interests were reviewed and are managed by Massachusetts General Hospital and Partners HealthCare in accordance with their conflict of interest policies. VW VCW is a scientific advisory board member of Triplet Therapeutics, a company developing new therapeutic approaches to address triplet repeat disorders such Huntington's disease and Myotonic Dystrophy and of LoQus23 Therapeutics. Her financial interests in Triplet Therapeutics, were reviewed and are managed by Massachusetts General Hospital and Partners HealthCare in accordance with their conflict of interest policies., (© 2020, Kovalenko et al.)

Objective: To investigate the genotype-phenotype correlation between neurofibromatosis 1 (NF1) germline mutations and imaging features of neurofibromas on whole-body MRI (WBMRI) by using radiomics image analysis techniques., Materials and Methods: Twenty-nine patients with NF1 who had known germline mutations determined by targeted next-generation sequencing were selected from a previous WBMRI study using coronal short tau inversion recovery sequence. Each tumor was segmented in WBMRI and a set of 59 imaging features was calculated using our in-house volumetric image analysis platform, 3DQI. A radiomics heatmap of 59 imaging features was analyzed to investigate the per-tumor and per-patient associations between the imaging features and mutation domains and mutation types. Linear mixed-effect models and one-way analysis of variance tests were performed to assess the similarity of tumor imaging features within mutation groups, between mutation groups, and between randomly selected groups., Results: A total of 218 neurofibromas (97 discrete neurofibromas and 121 plexiform neurofibromas) were identified in 19 of the 29 patients. The unsupervised hierarchical clustering in heatmap analysis revealed 6 major image feature patterns that were significantly correlated with gene mutation domains and types with strong to very strong associations of genotype-phenotype correlations in both per-tumor and per-patient studies ( p < 0.05, Cramer V > 0.5), whereas tumor size and locations showed no correlations with imaging features ( p = 0.79 and p = 0.42, respectively). The statistical analyses revealed that the number of significantly different features (SDFs) within mutation groups were significantly lower than those between mutation groups (mutation domains: 10.9 ± 9.5% vs 31.9 ± 23.8% and mutation types: 31.8 ± 30.7% vs 52.6 ± 29.3%). The first and second quartile p values of within-patient groups were more than 2 times higher than those between-patient groups. However, the numbers of SDFs between randomly selected groups were much lower (approximately 5.2%)., Conclusion: This preliminary study identified the NF1 radiogenomics linkage between NF1 causative mutations and MRI radiomic features, i.e., the correlation between NF1 genotype and imaging phenotype on WBMRI., (© 2020 American Academy of Neurology.)

Background: MBD5, encoding the methyl-CpG-binding domain 5 protein, has been proposed as a necessary and sufficient driver of the 2q23.1 microdeletion syndrome. De novo missense and protein-truncating variants from exome sequencing studies have directly implicated MBD5 in the etiology of autism spectrum disorder (ASD) and related neurodevelopmental disorders (NDDs). However, little is known concerning the specific function(s) of MBD5., Methods: To gain insight into the complex interactions associated with alteration of MBD5 in individuals with ASD and related NDDs, we explored the transcriptional landscape of MBD5 haploinsufficiency across multiple mouse brain regions of a heterozygous hypomorphic Mbd5 +/GT mouse model, and compared these results to CRISPR-mediated mutations of MBD5 in human iPSC-derived neuronal models., Results: Gene expression analyses across three brain regions from Mbd5 +/GT mice showed subtle transcriptional changes, with cortex displaying the most widespread changes following Mbd5 reduction, indicating context-dependent effects. Comparison with MBD5 reduction in human neuronal cells reinforced the context-dependence of gene expression changes due to MBD5 deficiency. Gene co-expression network analyses revealed gene clusters that were associated with reduced MBD5 expression and enriched for terms related to ciliary function., Limitations: These analyses included a limited number of mouse brain regions and neuronal models, and the effects of the gene knockdown are subtle. As such, these results will not reflect the full extent of MBD5 disruption across human brain regions during early neurodevelopment in ASD, or capture the diverse spectrum of cell-type-specific changes associated with MBD5 alterations., Conclusions: Our study points to modest and context-dependent transcriptional consequences of Mbd5 disruption in the brain. It also suggests a possible link between MBD5 and perturbations in ciliary function, which is an established pathogenic mechanism in developmental disorders and syndromes.

Objective: Defects in ion channels and neurotransmitter receptors are implicated in developmental and epileptic encephalopathy (DEE). Metabotropic glutamate receptor 7 (mGluR7), encoded by GRM7, is a presynaptic G-protein-coupled glutamate receptor critical for synaptic transmission. We previously proposed GRM7 as a candidate disease gene in two families with neurodevelopmental disorders (NDDs). One additional family has been published since. Here, we describe three additional families with GRM7 biallelic variants and deeply characterize the associated clinical neurological and electrophysiological phenotype and molecular data in 11 affected individuals from six unrelated families., Methods: Exome sequencing and family-based rare variant analyses on a cohort of 220 consanguineous families with NDDs revealed three families with GRM7 biallelic variants; three additional families were identified through literature search and collaboration with a clinical molecular laboratory., Results: We compared the observed clinical features and variants of 11 affected individuals from the six unrelated families. Identified novel deleterious variants included two homozygous missense variants (c.2671G>A:p.Glu891Lys and c.1973G>A:p.Arg685Gln) and one homozygous stop-gain variant (c.1975C>T:p.Arg659Ter). Developmental delay, neonatal- or infantile-onset epilepsy, and microcephaly were universal. Three individuals had hypothalamic-pituitary-axis dysfunction without pituitary structural abnormality. Neuroimaging showed cerebral atrophy and hypomyelination in a majority of cases. Two siblings demonstrated progressive loss of myelination by 2 years in both and an acquired microcephaly pattern in one. Five individuals died in early or late childhood., Conclusion: Detailed clinical characterization of 11 individuals from six unrelated families demonstrates that rare biallelic GRM7 pathogenic variants can cause DEEs, microcephaly, hypomyelination, and cerebral atrophy., (© 2020 The Authors. Annals of Clinical and Translational Neurology published by Wiley Periodicals, Inc on behalf of American Neurological Association.)

Background: Self-reported consanguinity is associated with risk for schizophrenia (SZ) in several inbred populations, but estimates using DNA-based coefficients of inbreeding are unavailable. Further, it is not known whether recessively inherited risk mutations can be identified through homozygosity by descent (HBD) mapping., Methods: We studied self-reported and DNA-based estimates of inbreeding among Egyptian patients with SZ (n = 421, DSM IV criteria) and adult controls without psychosis (n = 301), who were evaluated using semi-structured diagnostic interview schedules and genotyped using the Illumina Infinium PsychArray. Following quality control checks, coefficients of inbreeding (F) and regions of homozygosity (ROH) were estimated using PLINK software for HBD analysis. Exome sequencing was conducted in selected cases., Results: Inbreeding was associated with schizophrenia based on self-reported consanguinity (χ 2  = 4.506, 1 df, p = 0.034) and DNA-based estimates for inbreeding (F); the latter with a significant F × age interaction (β = 32.34, p = 0.0047). The association was most notable among patients older than age 40 years. Eleven ROH were over-represented in cases on chromosomes 1, 3, 6, 11, and 14; all but one region is novel for schizophrenia risk. Exome sequencing identified six recessively-acting genes in ROH with loss-of-function variants; one of which causes primary hereditary microcephaly., Conclusions: We propose consanguinity as an age-dependent risk factor for SZ in Egypt. HBD mapping is feasible for SZ in adequately powered samples., Competing Interests: Declaration of competing interest The Authors have declared no conflicts of interest in relation to the subject of this study., (Copyright © 2019 Elsevier B.V. All rights reserved.)

Background: Tuberous sclerosis complex (TSC) is a neurodevelopmental disorder with frequent occurrence of epilepsy, autism spectrum disorder (ASD), intellectual disability (ID), and tumors in multiple organs. The aberrant activation of mTORC1 in TSC has led to treatment with mTORC1 inhibitor rapamycin as a lifelong therapy for tumors, but TSC-associated neurocognitive manifestations remain unaffected by rapamycin., Methods: Here, we generated patient-specific, induced pluripotent stem cells (iPSCs) from a TSC patient with a heterozygous, germline, nonsense mutation in exon 15 of TSC1 and established an isogenic set of heterozygous (Het), null and corrected wildtype (Corr-WT) iPSCs using CRISPR/Cas9-mediated gene editing. We differentiated these iPSCs into neural progenitor cells (NPCs) and examined neurodevelopmental phenotypes, signaling and changes in gene expression by RNA-seq., Results: Differentiated NPCs revealed enlarged cell size in TSC1-Het and Null NPCs, consistent with mTORC1 activation. TSC1-Het and Null NPCs also revealed enhanced proliferation and altered neurite outgrowth in a genotype-dependent manner, which was not reversed by rapamycin. Transcriptome analyses of TSC1-NPCs revealed differentially expressed genes that display a genotype-dependent linear response, i.e., genes upregulated/downregulated in Het were further increased/decreased in Null. In particular, genes linked to ASD, epilepsy, and ID were significantly upregulated or downregulated warranting further investigation. In TSC1-Het and Null NPCs, we also observed basal activation of ERK1/2, which was further activated upon rapamycin treatment. Rapamycin also increased MNK1/2-eIF4E signaling in TSC1-deficient NPCs., Conclusion: MEK-ERK and MNK-eIF4E pathways regulate protein translation, and our results suggest that aberrant translation distinct in TSC1/2-deficient NPCs could play a role in neurodevelopmental defects. Our data showing upregulation of these signaling pathways by rapamycin support a strategy to combine a MEK or a MNK inhibitor with rapamycin that may be superior for TSC-associated CNS defects. Importantly, our generation of isogenic sets of NPCs from TSC patients provides a valuable platform for translatome and large-scale drug screening studies. Overall, our studies further support the notion that early developmental events such as NPC proliferation and initial process formation, such as neurite number and length that occur prior to neuronal differentiation, represent primary events in neurogenesis critical to disease pathogenesis of neurodevelopmental disorders such as ASD., Competing Interests: Competing interestsSJH is a member of the scientific advisory boards of Frequency Therapeutics, Psy Therapeutics, and Souvien Therapeutics, none of who were involved in this study. SDS is a scientific advisor for Outermost Therapeutics, Inc., which played no part in the present study. The other authors declare no competing interests., (© The Author(s). 2020.)

The 16p11.2 BP4-BP5 deletion and duplication syndromes are associated with a complex spectrum of neurodevelopmental phenotypes that includes developmental delay and autism spectrum disorder, with a reciprocal effect on head circumference, brain structure and body mass index. Mouse models of the 16p11.2 copy number variant have recapitulated some of the patient phenotypes, while studies in flies and zebrafish have uncovered several candidate contributory genes within the region, as well as complex genetic interactions. We evaluated one of these loci, KCTD13, by modeling haploinsufficiency and complete knockout in mice. In contrast to the zebrafish model, and in agreement with recent data, we found normal brain structure in heterozygous and homozygous mutants. However, recapitulating previously observed genetic interactions, we discovered sex-specific brain volumetric alterations in double heterozygous Kctd13xMvp and Kctd13xLat mice. Behavioral testing revealed a significant deficit in novel object recognition, novel location recognition and social transmission of food preference in Kctd13 mutants. These phenotypes were concomitant with a reduction in density of mature spines in the hippocampus, but potentially independent of RhoA abundance, which was unperturbed postnatally in our mutants. Furthermore, transcriptome analyses from cortex and hippocampus highlighted the dysregulation of pathways important in neurodevelopment, the most significant of which was synaptic formation. Together, these data suggest that KCTD13 contributes to the neurocognitive aspects of patients with the BP4-BP5 deletion, likely through genetic interactions with other loci., (© The Author(s) 2018. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Rare individuals with inactivating mutations in the Huntington's disease gene (HTT) exhibit variable abnormalities that imply essential HTT roles during organ development. Here we report phenotypes produced when increasingly severe hypomorphic mutations in the murine HTT orthologue Htt, (HdhneoQ20, HdhneoQ50, HdhneoQ111), were placed over a null allele (Hdhex4/5). The most severe hypomorphic allele failed to rescue null lethality at gastrulation, while the intermediate, though still severe, alleles yielded recessive perinatal lethality and a variety of fetal abnormalities affecting body size, skin, skeletal and ear formation, and transient defects in hematopoiesis. Comparative molecular analysis of wild-type and Htt-null retinoic acid-differentiated cells revealed gene network dysregulation associated with organ development that nominate polycomb repressive complexes and miRNAs as molecular mediators. Together these findings demonstrate that Htt is required both pre- and post-gastrulation to support normal development., Competing Interests: The authors have declared that no competing interests exist.

Neurofibromatosis 2 (NF2) is a rare tumor suppressor syndrome that manifests with multiple schwannomas and meningiomas. There are no effective drug therapies for these benign tumors and conventional therapies have limited efficacy. Various model systems have been created and several drug targets have been implicated in NF2-driven tumorigenesis based on known effects of the absence of merlin, the product of the NF2 gene. We tested priority compounds based on known biology with traditional dose-concentration studies in meningioma and schwann cell systems. Concurrently, we studied functional kinome and gene expression in these cells pre- and post-treatment to determine merlin deficient molecular phenotypes. Cell viability results showed that three agents (GSK2126458, Panobinostat, CUDC-907) had the greatest activity across schwannoma and meningioma cell systems, but merlin status did not significantly influence response. In vivo, drug effect was tumor specific with meningioma, but not schwannoma, showing response to GSK2126458 and Panobinostat. In culture, changes in both the transcriptome and kinome in response to treatment clustered predominantly based on tumor type. However, there were differences in both gene expression and functional kinome at baseline between meningioma and schwannoma cell systems that may form the basis for future selective therapies. This work has created an openly accessible resource (www.synapse.org/SynodosNF2) of fully characterized isogenic schwannoma and meningioma cell systems as well as a rich data source of kinome and transcriptome data from these assay systems before and after treatment that enables single and combination drug discovery based on molecular phenotype., Competing Interests: The authors have declared that no competing interests exist.

Schwannomatosis has been linked to germline mutations in the SMARCB1 and LZTR1 genes, and is frequently associated with pain.In a cohort study, we assessed the mutation status of 37 patients with clinically diagnosed schwannomatosis and compared to clinical data, whole body MRI (WBMRI), visual analog pain scale, and Short Form 36 (SF-36) bodily pain subscale.We identified a germline mutation in LZTR1 in 5 patients (13.5%) and SMARCB1 in 15 patients (40.5%), but found no germline mutation in 17 patients (45.9%). Peripheral schwannomas were detected in 3 LZTR1-mutant (60%) and 10 SMARCB1-mutant subjects (66.7%). Among those with peripheral tumors, the median tumor number was 4 in the LZTR1 group (median total body tumor volume 30 cc) and 10 in the SMARCB1 group (median volume 85cc), (P=.2915 for tumor number and P = .2289 for volume). mutation was associated with an increased prevalence of spinal schwannomas (100% vs 41%, P = .0197). The median pain score was 3.9/10 in the LZTR1 group and 0.5/10 in the SMARCB1 group (P = .0414), and SF-36 pain-associated quality of life was significantly worse in the LZTR1 group (P = .0106). Pain scores correlated with total body tumor volume (rho = 0.32471, P = .0499), but not with number of tumors (rho = 0.23065, P = .1696).We found no significant difference in quantitative tumor burden between mutational groups, but spinal schwannomas were more common in LZTR1-mutant patients. Pain was significantly higher in LZTR1-mutant than in SMARCB1-mutant patients, though spinal tumor location did not significantly correlate with pain. This suggests a possible genetic association with schwannomatosis-associated pain.

Genetic alterations of ARID1B have been recently recognized as one of the most common mendelian causes of intellectual disability and are associated with both syndromic and non-syndromic phenotypes. The ARID1B protein, a subunit of the chromatin remodeling complex SWI/SNF-A, is involved in the regulation of transcription and multiple downstream cellular processes. We report here the clinical, genetic, and proteomic phenotypes of an individual with a unique apparent de novo mutation of ARID1B due to an intragenic duplication. His neurodevelopmental phenotype includes a severe speech/language disorder with full scale IQ scores 78-98 and scattered academic skill levels, expanding the phenotypic spectrum of ARID1B mutations. Haploinsufficiency of ARID1B was determined both by RNA sequencing and quantitative RT-PCR. Fluorescence in situ hybridization analysis supported an intragenic localization of the ARID1B copy number gain. Principal component analysis revealed marked differentiation of the subject's lymphoblast proteome from that of controls. Of 3426 proteins quantified, 1014 were significantly up- or down-regulated compared to controls (q < 0.01). Pathway analysis revealed highly significant enrichment for canonical pathways of EIF2 and EIF4 signaling, protein ubiquitination, tRNA charging and chromosomal replication, among others. Network analyses revealed down-regulation of: (1) intracellular components involved in organization of membranes, organelles, and vesicles; (2) aspects of cell cycle control, signal transduction, and nuclear protein export; (3) ubiquitination and proteosomal function; and (4) aspects of mRNA synthesis/splicing. Further studies are needed to determine the detailed molecular and cellular mechanisms by which constitutional haploinsufficiency of ARID1B causes syndromic and non-syndromic developmental disabilities., (© 2017 Wiley Periodicals, Inc.)

Genetic variation within the transcription factor TCF4 locus can cause the intellectual disability and developmental disorder Pitt-Hopkins syndrome (PTHS), whereas single-nucleotide polymorphisms within noncoding regions are associated with schizophrenia. These genetic findings position TCF4 as a link between transcription and cognition; however, the neurobiology of TCF4 remains poorly understood. Here, we quantitated multiple distinct TCF4 transcript levels in human induced pluripotent stem cell-derived neural progenitors and differentiated neurons, and PTHS patient fibroblasts. We identify two classes of pharmacological treatments that regulate TCF4 expression: WNT pathway activation and inhibition of class I histone deacetylases. In PTHS fibroblasts, both of these perturbations upregulate a subset of TCF4 transcripts. Finally, using chromatin immunoprecipitation sequencing in conjunction with genome-wide transcriptome analysis, we identified TCF4 target genes that may mediate the effect of TCF4 loss on neuroplasticity. Our studies identify new pharmacological assays, tools, and targets for the development of therapeutics for cognitive disorders.

Arhinia, or absence of the nose, is a rare malformation of unknown etiology that is often accompanied by ocular and reproductive defects. Sequencing of 40 people with arhinia revealed that 84% of probands harbor a missense mutation localized to a constrained region of SMCHD1 encompassing the ATPase domain. SMCHD1 mutations cause facioscapulohumeral muscular dystrophy type 2 (FSHD2) via a trans-acting loss-of-function epigenetic mechanism. We discovered shared mutations and comparable DNA hypomethylation patterning between these distinct disorders. CRISPR/Cas9-mediated alteration of smchd1 in zebrafish yielded arhinia-relevant phenotypes. Transcriptome and protein analyses in arhinia probands and controls showed no differences in SMCHD1 mRNA or protein abundance but revealed regulatory changes in genes and pathways associated with craniofacial patterning. Mutations in SMCHD1 thus contribute to distinct phenotypic spectra, from craniofacial malformation and reproductive disorders to muscular dystrophy, which we speculate to be consistent with oligogenic mechanisms resulting in pleiotropic outcomes.

Integration of exogenous DNA into a host genome represents an important route to generate animal and cellular models for exploration into human disease and therapeutic development. In most models, little is known concerning structural integrity of the transgene, precise site of integration, or its impact on the host genome. We previously used whole-genome and targeted sequencing approaches to reconstruct transgene structure and integration sites in models of Huntington's disease, revealing complex structural rearrangements that can result from transgenesis. Here, we demonstrate in the R6/2 mouse, a widely used Huntington's disease model, that integration of a rearranged transgene with coincident deletion of 5,444 bp of host genome within the gene Gm12695 has striking molecular consequences. Gm12695, the function of which is unknown, is normally expressed at negligible levels in mouse brain, but transgene integration has resulted in cortical expression of a partial fragment (exons 8-11) 3' to the transgene integration site in R6/2. This transcript shows significant expression among the extensive network of differentially expressed genes associated with this model, including synaptic transmission, cell signalling and transcription. These data illustrate the value of sequence-level resolution of transgene insertions and transcription analysis to inform phenotypic characterization of transgenic models utilized in therapeutic research.

ATPase family AAA-domain containing protein 3A (ATAD3A) is a nuclear-encoded mitochondrial membrane protein implicated in mitochondrial dynamics, nucleoid organization, protein translation, cell growth, and cholesterol metabolism. We identified a recurrent de novo ATAD3A c.1582C>T (p.Arg528Trp) variant by whole-exome sequencing (WES) in five unrelated individuals with a core phenotype of global developmental delay, hypotonia, optic atrophy, axonal neuropathy, and hypertrophic cardiomyopathy. We also describe two families with biallelic variants in ATAD3A, including a homozygous variant in two siblings, and biallelic ATAD3A deletions mediated by nonallelic homologous recombination (NAHR) between ATAD3A and gene family members ATAD3B and ATAD3C. Tissue-specific overexpression of bor R534W , the Drosophila mutation homologous to the human c.1582C>T (p.Arg528Trp) variant, resulted in a dramatic decrease in mitochondrial content, aberrant mitochondrial morphology, and increased autophagy. Homozygous null bor larvae showed a significant decrease of mitochondria, while overexpression of bor WT resulted in larger, elongated mitochondria. Finally, fibroblasts of an affected individual exhibited increased mitophagy. We conclude that the p.Arg528Trp variant functions through a dominant-negative mechanism that results in small mitochondria that trigger mitophagy, resulting in a reduction in mitochondrial content. ATAD3A variation represents an additional link between mitochondrial dynamics and recognizable neurological syndromes, as seen with MFN2, OPA1, DNM1L, and STAT2 mutations., (Copyright © 2016 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Recurrent, reciprocal genomic disorders resulting from non-allelic homologous recombination (NAHR) between near-identical segmental duplications (SDs) are a major cause of human disease, often producing phenotypically distinct syndromes. The genomic architecture of flanking SDs presents a challenge for modeling these syndromes; however, the capability to efficiently generate reciprocal copy number variants (CNVs) that mimic NAHR would represent a valuable modeling tool. We describe here a CRISPR/Cas9 genome engineering method, single-guide CRISPR/Cas targeting of repetitive elements (SCORE), to model reciprocal genomic disorders and demonstrate its capabilities by generating reciprocal CNVs of 16p11.2 and 15q13.3, including alteration of one copy-equivalent of the SDs that mediate NAHR in vivo. The method is reproducible, and RNA sequencing reliably clusters transcriptional signatures from human subjects with in vivo CNVs and their corresponding in vitro models. This new approach will provide broad applicability for the study of genomic disorders and, with further development, may also permit efficient correction of these defects.

Development of the human nervous system involves complex interactions among fundamental cellular processes and requires a multitude of genes, many of which remain to be associated with human disease. We applied whole exome sequencing to 128 mostly consanguineous families with neurogenetic disorders that often included brain malformations. Rare variant analyses for both single nucleotide variant (SNV) and copy number variant (CNV) alleles allowed for identification of 45 novel variants in 43 known disease genes, 41 candidate genes, and CNVs in 10 families, with an overall potential molecular cause identified in >85% of families studied. Among the candidate genes identified, we found PRUNE, VARS, and DHX37 in multiple families and homozygous loss-of-function variants in AGBL2, SLC18A2, SMARCA1, UBQLN1, and CPLX1. Neuroimaging and in silico analysis of functional and expression proximity between candidate and known disease genes allowed for further understanding of genetic networks underlying specific types of brain malformations., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Klippel-Feil syndrome is a rare disorder represented by a subgroup of segmentation defects of the vertebrae and characterized by fusion of the cervical vertebrae, low posterior hairline, and short neck with limited motion. Both autosomal dominant and recessive inheritance patterns were reported in families with Klippel-Feil. Mutated genes for both dominant (GDF6 and GDF3) and recessive (MEOX1) forms of Klippel-Feil syndrome have been shown to be involved in somite development via transcription regulation and signaling pathways. Heterotaxy arises from defects in proteins that function in the development of left-right asymmetry of the developing embryo. We describe a consanguineous family with a male proband who presents with classical Klippel-Feil syndrome together with heterotaxy (situs inversus totalis). The present patient also had Sprengel's deformity, deformity of the sternum, and a solitary kidney. Using exome sequencing, we identified a homozygous frameshift mutation (c.299delT; p.L100fs) in RIPPLY2, a gene shown to play a crucial role in somitogenesis and participate in the Notch signaling pathway via negatively regulating Tbx6. Our data confirm RIPPLY2 as a novel gene for autosomal recessive Klippel-Feil syndrome, and in addition-from a mechanistic standpoint-suggest the possibility that mutations in RIPPLY2 could also lead to heterotaxy. © 2015 Wiley Periodicals, Inc., (© 2015 Wiley Periodicals, Inc.)

The CAG repeat expansion in the Huntington's disease gene HTT extends a polyglutamine tract in mutant huntingtin that enhances its ability to facilitate polycomb repressive complex 2 (PRC2). To gain insight into this dominant gain of function, we mapped histone modifications genome-wide across an isogenic panel of mouse embryonic stem cell (ESC) and neuronal progenitor cell (NPC) lines, comparing the effects of Htt null and different size Htt CAG mutations. We found that Htt is required in ESC for the proper deposition of histone H3K27me3 at a subset of 'bivalent' loci but in NPC it is needed at 'bivalent' loci for both the proper maintenance and the appropriate removal of this mark. In contrast, Htt CAG size, though changing histone H3K27me3, is prominently associated with altered histone H3K4me3 at 'active' loci. The sets of ESC and NPC genes with altered histone marks delineated by the lack of huntingtin or the presence of mutant huntingtin, though distinct, are enriched in similar pathways with apoptosis specifically highlighted for the CAG mutation. Thus, the manner by which huntingtin function facilitates PRC2 may afford mutant huntingtin with multiple opportunities to impinge upon the broader machinery that orchestrates developmentally appropriate chromatin status., (© The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Context: Pituitary stalk interruption syndrome (PSIS) is a rare, congenital anomaly of the pituitary gland characterized by pituitary gland insufficiency, thin or discontinuous pituitary stalk, anterior pituitary hypoplasia, and ectopic positioning of the posterior pituitary gland (neurohypophysis). The clinical presentation of patients with PSIS varies from isolated growth hormone (GH) deficiency to combined pituitary insufficiency and accompanying extrapituitary findings. Mutations in HESX1, LHX4, OTX2, SOX3, and PROKR2 have been associated with PSIS in less than 5% of cases; thus, the underlying genetic etiology for the vast majority of cases remains to be determined., Objective: We applied whole-exome sequencing (WES) to a consanguineous family with two affected siblings who have pituitary gland insufficiency and radiographic findings of hypoplastic (thin) pituitary gland, empty sella, ectopic neurohypophysis, and interrupted pitiutary stalk-characteristic clinical diagnostic findings of PSIS., Design and Participants: WES was applied to two affected and one unaffected siblings., Results: WES of two affected and one unaffected sibling revealed a unique homozygous missense mutation in GPR161, which encodes the orphan G protein-coupled receptor 161, a protein responsible for transducing extracellular signals across the plasma membrane into the cell., Conclusion: Mutations of GPR161 may be implicated as a potential novel cause of PSIS.

Truncating mutations of chromodomain helicase DNA-binding protein 8 (CHD8), and of many other genes with diverse functions, are strong-effect risk factors for autism spectrum disorder (ASD), suggesting multiple mechanisms of pathogenesis. We explored the transcriptional networks that CHD8 regulates in neural progenitor cells (NPCs) by reducing its expression and then integrating transcriptome sequencing (RNA sequencing) with genome-wide CHD8 binding (ChIP sequencing). Suppressing CHD8 to levels comparable with the loss of a single allele caused altered expression of 1,756 genes, 64.9% of which were up-regulated. CHD8 showed widespread binding to chromatin, with 7,324 replicated sites that marked 5,658 genes. Integration of these data suggests that a limited array of direct regulatory effects of CHD8 produced a much larger network of secondary expression changes. Genes indirectly down-regulated (i.e., without CHD8-binding sites) reflect pathways involved in brain development, including synapse formation, neuron differentiation, cell adhesion, and axon guidance, whereas CHD8-bound genes are strongly associated with chromatin modification and transcriptional regulation. Genes associated with ASD were strongly enriched among indirectly down-regulated loci (P < 10(-8)) and CHD8-bound genes (P = 0.0043), which align with previously identified coexpression modules during fetal development. We also find an intriguing enrichment of cancer-related gene sets among CHD8-bound genes (P < 10(-10)). In vivo suppression of chd8 in zebrafish produced macrocephaly comparable to that of humans with inactivating mutations. These data indicate that heterozygous disruption of CHD8 precipitates a network of gene-expression changes involved in neurodevelopmental pathways in which many ASD-associated genes may converge on shared mechanisms of pathogenesis.

GAPO syndrome (OMIM#230740) is the acronym for growth retardation, alopecia, pseudoanodontia, and optic atrophy. About 35 cases have been reported, making it among one of the rarest recessive conditions. Distinctive craniofacial features including alopecia, rarefaction of eyebrows and eyelashes, frontal bossing, high forehead, mid-facial hypoplasia, hypertelorism, and thickened eyelids and lips make GAPO syndrome a clinically recognizable phenotype. While this genomic study was in progress mutations in ANTXR1 were reported to cause GAPO syndrome. In our study we performed whole exome sequencing (WES) for five affected individuals from three Turkish kindreds segregating the GAPO trait. Exome sequencing analysis identified three novel homozygous mutations including; one frame-shift (c.1220&#95;1221insT; p.Ala408Cysfs*2), one splice site (c.411A>G; p.Gln137Gln), and one non-synonymous (c.1150G>A; p.Gly384Ser) mutation in the ANTXR1 gene. Our studies expand the allelic spectrum in this rare condition and potentially provide insight into the role of ANTXR1 in the regulation of the extracellular matrix., (© 2014 Wiley Periodicals, Inc.)

Genome editing has attracted wide interest for the generation of cellular models of disease using human pluripotent stem cells and other cell types. CRISPR-Cas systems and TALENs can target desired genomic sites with high efficiency in human cells, but recent publications have led to concern about the extent to which these tools may cause off-target mutagenic effects that could potentially confound disease-modeling studies. Using CRISPR-Cas9 and TALEN targeted human pluripotent stem cell clones, we performed whole-genome sequencing at high coverage in order to assess the degree of mutagenesis across the entire genome. In both types of clones, we found that off-target mutations attributable to the nucleases were very rare. From this analysis, we suggest that, although some cell types may be at risk for off-target mutations, the incidence of such effects in human pluripotent stem cells may be sufficiently low and thus not a significant concern for disease modeling and other applications., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Reciprocal copy-number variation (CNV) of a 593 kb region of 16p11.2 is a common genetic cause of autism spectrum disorder (ASD), yet it is not completely penetrant and can manifest in a wide array of phenotypes. To explore its molecular consequences, we performed RNA sequencing of cerebral cortex from mouse models with CNV of the syntenic 7qF3 region and lymphoblast lines from 34 members of 7 multiplex ASD-affected families harboring the 16p11.2 CNV. Expression of all genes in the CNV region correlated well with their DNA copy number, with no evidence of dosage compensation. We observed effects on gene expression outside the CNV region, including apparent positional effects in cis and in trans at genomic segments with evidence of physical interaction in Hi-C chromosome conformation data. One of the most significant positional effects was telomeric to the 16p11.2 CNV and includes the previously described "distal" 16p11.2 microdeletion. Overall, 16p11.2 CNV was associated with altered expression of genes and networks that converge on multiple hypotheses of ASD pathogenesis, including synaptic function (e.g., NRXN1, NRXN3), chromatin modification (e.g., CHD8, EHMT1, MECP2), transcriptional regulation (e.g., TCF4, SATB2), and intellectual disability (e.g., FMR1, CEP290). However, there were differences between tissues and species, with the strongest effects being consistently within the CNV region itself. Our analyses suggest that through a combination of indirect regulatory effects and direct effects on nuclear architecture, alteration of 16p11.2 genes disrupts expression networks that involve other genes and pathways known to contribute to ASD, suggesting an overlap in mechanisms of pathogenesis., (Copyright © 2014 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

CLP1 is a RNA kinase involved in tRNA splicing. Recently, CLP1 kinase-dead mice were shown to display a neuromuscular disorder with loss of motor neurons and muscle paralysis. Human genome analyses now identified a CLP1 homozygous missense mutation (p.R140H) in five unrelated families, leading to a loss of CLP1 interaction with the tRNA splicing endonuclease (TSEN) complex, largely reduced pre-tRNA cleavage activity, and accumulation of linear tRNA introns. The affected individuals develop severe motor-sensory defects, cortical dysgenesis, and microcephaly. Mice carrying kinase-dead CLP1 also displayed microcephaly and reduced cortical brain volume due to the enhanced cell death of neuronal progenitors that is associated with reduced numbers of cortical neurons. Our data elucidate a neurological syndrome defined by CLP1 mutations that impair tRNA splicing. Reduction of a founder mutation to homozygosity illustrates the importance of rare variations in disease and supports the clan genomics hypothesis., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Structural Genomics aims to elucidate protein structures to identify their functions. Unfortunately, the variation of just a few residues can be enough to alter activity or binding specificity and limit the functional resolution of annotations based on sequence and structure; in enzymes, substrates are especially difficult to predict. Here, large-scale controls and direct experiments show that the local similarity of five or six residues selected because they are evolutionarily important and on the protein surface can suffice to identify an enzyme activity and substrate. A motif of five residues predicted that a previously uncharacterized Silicibacter sp. protein was a carboxylesterase for short fatty acyl chains, similar to hormone-sensitive-lipase-like proteins that share less than 20% sequence identity. Assays and directed mutations confirmed this activity and showed that the motif was essential for catalysis and substrate specificity. We conclude that evolutionary and structural information may be combined on a Structural Genomics scale to create motifs of mixed catalytic and noncatalytic residues that identify enzyme activity and substrate specificity.

Motivation: The constraints under which sequence, structure and function coevolve are not fully understood. Bringing this mutual relationship to light can reveal the molecular basis of binding, catalysis and allostery, thereby identifying function and rationally guiding protein redesign. Underlying these relationships are the epistatic interactions that occur when the consequences of a mutation to a protein are determined by the genetic background in which it occurs. Based on prior data, we hypothesize that epistatic forces operate most strongly between residues nearby in the structure, resulting in smooth evolutionary importance across the structure., Methods and Results: We find that when residue scores of evolutionary importance are distributed smoothly between nearby residues, functional site prediction accuracy improves. Accordingly, we designed a novel measure of evolutionary importance that focuses on the interaction between pairs of structurally neighboring residues. This measure that we term pair-interaction Evolutionary Trace yields greater functional site overlap and better structure-based proteome-wide functional predictions., Conclusions: Our data show that the structural smoothness of evolutionary importance is a fundamental feature of the coevolution of sequence, structure and function. Mutations operate on individual residues, but selective pressure depends in part on the extent to which a mutation perturbs interactions with neighboring residues. In practice, this principle led us to redefine the importance of a residue in terms of the importance of its epistatic interactions with neighbors, yielding better annotation of functional residues, motivating experimental validation of a novel functional site in LexA and refining protein function prediction., Contact: lichtarge@bcm.edu., Supplementary Information: Supplementary data are available at Bioinformatics online.

Automated annotation of protein function is challenging. As the number of sequenced genomes rapidly grows, the overwhelming majority of protein products can only be annotated computationally. If computational predictions are to be relied upon, it is crucial that the accuracy of these methods be high. Here we report the results from the first large-scale community-based critical assessment of protein function annotation (CAFA) experiment. Fifty-four methods representing the state of the art for protein function prediction were evaluated on a target set of 866 proteins from 11 organisms. Two findings stand out: (i) today's best protein function prediction algorithms substantially outperform widely used first-generation methods, with large gains on all types of targets; and (ii) although the top methods perform well enough to guide experiments, there is considerable need for improvement of currently available tools.