Application of multiple mosaic callers improves post-zygotic mutation detection from exome sequencing data.

PURPOSE: The gold standard for identification of post-zygotic variants (PZVs) is droplet digital PCR (ddPCR) or high-depth sequencing across multiple tissues types. These approaches are yet to be systematically implemented for monogenic disorders. We developed PZV detection pipelines for correct classification of de novo variants. METHOD: Our pipelines detect PZV in parents (gonosomal mosaicism "pGoM") and children (somatic mosaicism, "M3"). We applied them to research exome sequencing (ES) data from The Australian Cerebral Palsy Biobank (ACPB, n=145 trios) and Simons Simplex Collection (SSC, n=405 families). Candidate mosaic variants were validated using deep amplicon sequencing or ddPCR. RESULTS: 69.2% (M3trio), 63.9% (M3single) and 92.7% (pGoM) of detected variants were validated, with 48.6%, 56.7% and 26.2% of variants respectively meeting strict criteria for mosaicism. In the ACPB, 16.6% of probands and 20.7% of parents had at least one true positive somatic or pGoM variant respectively. A large proportion of PZVs detected in SSC parents (79.8%) and child (94.5%) were not previously reported. We reclassified 3.7-8.0% of germline de novo variants as mosaic. CONCLUSION: Many PZVs were incorrectly classified as germline variants or missed by previous approaches. Systematic application of our pipelines could increase genetic diagnostic rate, improve estimates of recurrence risk in families, and benefit novel disease gene identification.

RNA variant assessment using transactivation and transdifferentiation.

Understanding the impact of splicing and nonsense variants on RNA is crucial for the resolution of variant classification as well as their suitability for precision medicine interventions. This is primarily enabled through RNA studies involving transcriptomics followed by targeted assays using RNA isolated from clinically accessible tissues (CATs) such as blood or skin of affected individuals. Insufficient disease gene expression in CATs does however pose a major barrier to RNA based investigations, which we show is relevant to 1,436 Mendelian disease genes. We term these "silent" Mendelian genes (SMGs), the largest portion (36%) of which are associated with neurological disorders. We developed two approaches to induce SMG expression in human dermal fibroblasts (HDFs) to overcome this limitation, including CRISPR-activation-based gene transactivation and fibroblast-to-neuron transdifferentiation. Initial transactivation screens involving 40 SMGs stimulated our development of a highly multiplexed transactivation system culminating in the 6- to 90,000-fold induction of expression of 20/20 (100%) SMGs tested in HDFs. Transdifferentiation of HDFs directly to neurons led to expression of 193/516 (37.4%) of SMGs implicated in neurological disease. The magnitude and isoform diversity of SMG expression following either transactivation or transdifferentiation was comparable to clinically relevant tissues. We apply transdifferentiation and/or gene transactivation combined with short- and long-read RNA sequencing to investigate the impact that variants in USH2A, SCN1A, DMD, and PAK3 have on RNA using HDFs derived from affected individuals. Transactivation and transdifferentiation represent rapid, scalable functional genomic solutions to investigate variants impacting SMGs in the patient cell and genomic context.

Heterozygous loss-of-function variants in DOCK4 cause neurodevelopmental delay and microcephaly.

Neurons form the basic anatomical and functional structure of the nervous system, and defects in neuronal differentiation or formation of neurites are associated with various psychiatric and neurodevelopmental disorders. Dynamic changes in the cytoskeleton are essential for this process, which is, inter alia, controlled by the dedicator of cytokinesis 4 (DOCK4) through the activation of RAC1. Here, we clinically describe 7 individuals (6 males and one female) with variants in DOCK4 and overlapping phenotype of mild to severe global developmental delay. Additional symptoms include coordination or gait abnormalities, microcephaly, nonspecific brain malformations, hypotonia and seizures. Four individuals carry missense variants (three of them detected de novo) and three individuals carry null variants (two of them maternally inherited). Molecular modeling of the heterozygous missense variants suggests that the majority of them affect the globular structure of DOCK4. In vitro functional expression studies in transfected Neuro-2A cells showed that all missense variants impaired neurite outgrowth. Furthermore, Dock4 knockout Neuro-2A cells also exhibited defects in promoting neurite outgrowth. Our results, including clinical, molecular and functional data, suggest that loss-of-function variants in DOCK4 probable cause a variable spectrum of a novel neurodevelopmental disorder with microcephaly.

Proteomic analysis of the developing mammalian brain links PCDH19 to the Wnt/ß-catenin signalling pathway.

Clustering Epilepsy (CE) is a neurological disorder caused by pathogenic variants of the Protocadherin 19 (PCDH19) gene. PCDH19 encodes a protein involved in cell adhesion and Estrogen Receptor α mediated-gene regulation. To gain further insights into the molecular role of PCDH19 in the brain, we investigated the PCDH19 interactome in the developing mouse hippocampus and cortex. Combined with a meta-analysis of all reported PCDH19 interacting proteins, our results show that PCDH19 interacts with proteins involved in actin, microtubule, and gene regulation. We report CAPZA1, αN-catenin and, importantly, ß-catenin as novel PCDH19 interacting proteins. Furthermore, we show that PCDH19 is a regulator of ß-catenin transcriptional activity, and that this pathway is disrupted in CE individuals. Overall, our results support the involvement of PCDH19 in the cytoskeletal network and point to signalling pathways where PCDH19 plays critical roles.

Compromised transcription-mRNA export factor THOC2 causes R-loop accumulation, DNA damage and adverse neurodevelopment.

We implicated the X-chromosome THOC2 gene, which encodes the largest subunit of the highly-conserved TREX (Transcription-Export) complex, in a clinically complex neurodevelopmental disorder with intellectual disability as the core phenotype. To study the molecular pathology of this essential eukaryotic gene, we generated a mouse model based on a hypomorphic Thoc2 exon 37-38 deletion variant of a patient with ID, speech delay, hypotonia, and microcephaly. The Thoc2 exon 37-38 deletion male (Thoc2Δ/Y) mice recapitulate the core phenotypes of THOC2 syndrome including smaller size and weight, and significant deficits in spatial learning, working memory and sensorimotor functions. The Thoc2Δ/Y mouse brain development is significantly impacted by compromised THOC2/TREX function resulting in R-loop accumulation, DNA damage and consequent cell death. Overall, we suggest that perturbed R-loop homeostasis, in stem cells and/or differentiated cells in mice and the patient, and DNA damage-associated functional alterations are at the root of THOC2 syndrome.

Multiomic analysis implicates nuclear hormone receptor signalling in clustering epilepsy.

Clustering Epilepsy (CE) is an epileptic disorder with neurological comorbidities caused by heterozygous variants of the X chromosome gene Protocadherin 19 (PCDH19). Recent studies have implicated dysregulation of the Nuclear Hormone Receptor (NHR) pathway in CE pathogenesis. To obtain a comprehensive overview of the impact and mechanisms of loss of PCDH19 function in CE pathogenesis, we have performed epigenomic, transcriptomic and proteomic analysis of CE relevant models. Our studies identified differential regulation and expression of Androgen Receptor (AR) and its targets in CE patient skin fibroblasts. Furthermore, our cell culture assays revealed the repression of PCDH19 expression mediated through ERα and the co-regulator FOXA1. We also identified a protein-protein interaction between PCDH19 and AR, expanding upon the intrinsic link between PCDH19 and the NHR pathway. Together, these results point to a novel mechanism of NHR signaling in the pathogenesis of CE that can be explored for potential therapeutic options.

Macrocephaly and developmental delay caused by missense variants in RAB5C.

Rab GTPases are important regulators of intracellular vesicular trafficking. RAB5C is a member of the Rab GTPase family that plays an important role in the endocytic pathway, membrane protein recycling and signaling. Here we report on 12 individuals with nine different heterozygous de novo variants in RAB5C. All but one patient with missense variants (n = 9) exhibited macrocephaly, combined with mild-to-moderate developmental delay. Patients with loss of function variants (n = 2) had an apparently more severe clinical phenotype with refractory epilepsy and intellectual disability but a normal head circumference. Four missense variants were investigated experimentally. In vitro biochemical studies revealed that all four variants were damaging, resulting in increased nucleotide exchange rate, attenuated responsivity to guanine exchange factors and heterogeneous effects on interactions with effector proteins. Studies in C. elegans confirmed that all four variants were damaging in vivo and showed defects in endocytic pathway function. The variant heterozygotes displayed phenotypes that were not observed in null heterozygotes, with two shown to be through a dominant negative mechanism. Expression of the human RAB5C variants in zebrafish embryos resulted in defective development, further underscoring the damaging effects of the RAB5C variants. Our combined bioinformatic, in vitro and in vivo experimental studies and clinical data support the association of RAB5C missense variants with a neurodevelopmental disorder characterized by macrocephaly and mild-to-moderate developmental delay through disruption of the endocytic pathway.

Aicardi Syndrome Is a Genetically Heterogeneous Disorder.

Aicardi Syndrome (AIC) is a rare neurodevelopmental disorder recognized by the classical triad of agenesis of the corpus callosum, chorioretinal lacunae and infantile epileptic spasms syndrome. The diagnostic criteria of AIC were revised in 2005 to include additional phenotypes that are frequently observed in this patient group. AIC has been traditionally considered as X-linked and male lethal because it almost exclusively affects females. Despite numerous genetic and genomic investigations on AIC, a unifying X-linked cause has not been identified. Here, we performed exome and genome sequencing of 10 females with AIC or suspected AIC based on current criteria. We identified a unique de novo variant, each in different genes: KMT2B, SLF1, SMARCB1, SZT2 and WNT8B, in five of these females. Notably, genomic analyses of coding and non-coding single nucleotide variants, short tandem repeats and structural variation highlighted a distinct lack of X-linked candidate genes. We assessed the likely pathogenicity of our candidate autosomal variants using the TOPflash assay for WNT8B and morpholino knockdown in zebrafish (Danio rerio) embryos for other candidates. We show expression of Wnt8b and Slf1 are restricted to clinically relevant cortical tissues during mouse development. Our findings suggest that AIC is genetically heterogeneous with implicated genes converging on molecular pathways central to cortical development.

Redefining cerebral palsies as a diverse group of neurodevelopmental disorders with genetic aetiology.

Cerebral palsy is a clinical descriptor covering a diverse group of permanent, non-degenerative disorders of motor function. Around one-third of cases have now been shown to have an underlying genetic aetiology, with the genetic landscape overlapping with those of neurodevelopmental disorders including intellectual disability, epilepsy, speech and language disorders and autism. Here we review the current state of genomic testing in cerebral palsy, highlighting the benefits for personalized medicine and the imperative to consider aetiology during clinical diagnosis. With earlier clinical diagnosis now possible, we emphasize the opportunity for comprehensive and early genomic testing as a crucial component of the routine diagnostic work-up in people with cerebral palsy.

Genomic and phenotypic characterization of 404 individuals with neurodevelopmental disorders caused by CTNNB1 variants.

PURPOSE: Germline loss-of-function variants in CTNNB1 cause neurodevelopmental disorder with spastic diplegia and visual defects (NEDSDV; OMIM 615075) and are the most frequent, recurrent monogenic cause of cerebral palsy (CP). We investigated the range of clinical phenotypes owing to disruptions of CTNNB1 to determine the association between NEDSDV and CP. METHODS: Genetic information from 404 individuals with collectively 392 pathogenic CTNNB1 variants were ascertained for the study. From these, detailed phenotypes for 52 previously unpublished individuals were collected and combined with 68 previously published individuals with comparable clinical information. The functional effects of selected CTNNB1 missense variants were assessed using TOPFlash assay. RESULTS: The phenotypes associated with pathogenic CTNNB1 variants were similar. A diagnosis of CP was not significantly associated with any set of traits that defined a specific phenotypic subgroup, indicating that CP is not additional to NEDSDV. Two CTNNB1 missense variants were dominant negative regulators of WNT signaling, highlighting the utility of the TOPFlash assay to functionally assess variants. CONCLUSION: NEDSDV is a clinically homogeneous disorder irrespective of initial clinical diagnoses, including CP, or entry points for genetic testing.

Bi-allelic variants in neuronal cell adhesion molecule cause a neurodevelopmental disorder characterized by developmental delay, hypotonia, neuropathy/spasticity.

Cell adhesion molecules are membrane-bound proteins predominantly expressed in the central nervous system along principal axonal pathways with key roles in nervous system development, neural cell differentiation and migration, axonal growth and guidance, myelination, and synapse formation. Here, we describe ten affected individuals with bi-allelic variants in the neuronal cell adhesion molecule NRCAM that lead to a neurodevelopmental syndrome of varying severity; the individuals are from eight families. This syndrome is characterized by developmental delay/intellectual disability, hypotonia, peripheral neuropathy, and/or spasticity. Computational analyses of NRCAM variants, many of which cluster in the third fibronectin type III (Fn-III) domain, strongly suggest a deleterious effect on NRCAM structure and function, including possible disruption of its interactions with other proteins. These findings are corroborated by previous in vitro studies of murine Nrcam-deficient cells, revealing abnormal neurite outgrowth, synaptogenesis, and formation of nodes of Ranvier on myelinated axons. Our studies on zebrafish nrcamaΔ mutants lacking the third Fn-III domain revealed that mutant larvae displayed significantly altered swimming behavior compared to wild-type larvae (p < 0.03). Moreover, nrcamaΔ mutants displayed a trend toward increased amounts of α-tubulin fibers in the dorsal telencephalon, demonstrating an alteration in white matter tracts and projections. Taken together, our study provides evidence that NRCAM disruption causes a variable form of a neurodevelopmental disorder and broadens the knowledge on the growing role of the cell adhesion molecule family in the nervous system.

Mutations disrupting neuritogenesis genes confer risk for cerebral palsy.

In addition to commonly associated environmental factors, genomic factors may cause cerebral palsy. We performed whole-exome sequencing of 250 parent-offspring trios, and observed enrichment of damaging de novo mutations in cerebral palsy cases. Eight genes had multiple damaging de novo mutations; of these, two (TUBA1A and CTNNB1) met genome-wide significance. We identified two novel monogenic etiologies, FBXO31 and RHOB, and showed that the RHOB mutation enhances active-state Rho effector binding while the FBXO31 mutation diminishes cyclin D levels. Candidate cerebral palsy risk genes overlapped with neurodevelopmental disorder genes. Network analyses identified enrichment of Rho GTPase, extracellular matrix, focal adhesion and cytoskeleton pathways. Cerebral palsy risk genes in enriched pathways were shown to regulate neuromotor function in a Drosophila reverse genetics screen. We estimate that 14% of cases could be attributed to an excess of damaging de novo or recessive variants. These findings provide evidence for genetically mediated dysregulation of early neuronal connectivity in cerebral palsy.

Definition and diagnosis of cerebral palsy in genetic studies: a systematic review.

AIM: To conduct a systematic review of phenotypic definition and case ascertainment in published genetic studies of cerebral palsy (CP) to inform guidelines for the reporting of such studies. METHOD: Inclusion criteria comprised genetic studies of candidate genes, with CP as the outcome, published between 1990 and 2019 in the PubMed, Embase, and BIOSIS Citation Index databases. RESULTS: Fifty-seven studies met the inclusion criteria. We appraised how CP was defined, the quality of information on case ascertainment, and compliance with international consensus guidelines. Seven studies (12%) were poorly described, 33 studies (58%) gave incomplete information, and 17 studies (30%) were well described. Missing key information precluded determining how many studies complied with the definition by Rosenbaum et al. Only 18 out of 57 studies (32%) were compliant with the Surveillance of Cerebral Palsy in Europe (SCPE) international guidelines on defining CP. INTERPRETATION: Limited compliance with international consensus guidelines on phenotypic definition and mediocre reporting of CP case ascertainment hinders the comparison of results among genetic studies of CP (including meta-analyses), thereby limiting the quality, interpretability, and generalizability of study findings. Compliance with the SCPE guidelines is important for ongoing gene discovery efforts in CP, given the potential for misclassification of unrelated neurological conditions as CP.

Erratum: Author Correction: Pathogenic copy number variants that affect gene expression contribute to genomic burden in cerebral palsy.

[This corrects the article DOI: 10.1038/s41525-018-0073-4.].

Non-self mutation: double-stranded RNA elicits antiviral pathogenic response in a Drosophila model of expanded CAG repeat neurodegenerative diseases.

Inflammation is activated prior to symptoms in neurodegenerative diseases, providing a plausible pathogenic mechanism. Indeed, genetic and pharmacological ablation studies in animal models of several neurodegenerative diseases demonstrate that inflammation is required for pathology. However, while there is growing evidence that inflammation-mediated pathology may be the common mechanism underlying neurodegenerative diseases, including those due to dominantly inherited expanded repeats, the proximal causal agent is unknown. Expanded CAG.CUG repeat double-stranded RNA causes inflammation-mediated pathology when expressed in Drosophila. Repeat dsRNA is recognized by Dicer-2 as a foreign or 'non-self' molecule triggering both antiviral RNA and RNAi pathways. Neither of the RNAi pathway cofactors R2D2 nor loquacious are necessary, indicating antiviral RNA activation. RNA modification enables avoidance of recognition as 'non-self' by the innate inflammatory surveillance system. Human ADAR1 edits RNA conferring 'self' status and when co-expressed with expanded CAG.CUG dsRNA in Drosophila the pathology is lost. Cricket Paralysis Virus protein CrPV-1A is a known antagonist of Argonaute-2 in Drosophila antiviral defense. CrPV-1A co-expression also rescues pathogenesis, confirming anti-viral-RNA response. Repeat expansion mutation therefore confers 'non-self' recognition of endogenous RNA, thereby providing a proximal, autoinflammatory trigger for expanded repeat neurodegenerative diseases.

Genetic or Other Causation Should Not Change the Clinical Diagnosis of Cerebral Palsy.

High throughput sequencing is discovering many likely causative genetic variants in individuals with cerebral palsy. Some investigators have suggested that this changes the clinical diagnosis of cerebral palsy and that these individuals should be removed from this diagnostic category. Cerebral palsy is a neurodevelopmental disorder diagnosed on clinical signs, not etiology. All nonprogressive permanent disorders of movement and posture attributed to disturbances that occurred in the developing fetal and infant brain can be described as "cerebral palsy." This definition of cerebral palsy should not be changed, whatever the cause. Reasons include stability, utility and accuracy of cerebral palsy registers, direct access to services, financial and social support specifically offered to families with cerebral palsy, and community understanding of the clinical diagnosis. Other neurodevelopmental disorders, for example, epilepsy, have not changed the diagnosis when genomic causes are found. The clinical diagnosis of cerebral palsy should remain, should prompt appropriate genetic studies and can subsequently be subclassified by etiology.

Pathogenic copy number variants that affect gene expression contribute to genomic burden in cerebral palsy.

Cerebral palsy (CP) is the most frequent movement disorder of childhood affecting 1 in 500 live births in developed countries. We previously identified likely pathogenic de novo or inherited single nucleotide variants (SNV) in 14% (14/98) of trios by exome sequencing and a further 5% (9/182) from evidence of outlier gene expression using RNA sequencing. Here, we detected copy number variants (CNV) from exomes of 186 unrelated individuals with CP (including our original 98 trios) using the CoNIFER algorithm. CNV were validated with Illumina 850 K SNP arrays and compared with RNA-Seq outlier gene expression analysis from lymphoblastoid cell lines (LCL). Gene expression was highly correlated with gene dosage effect. We resolved an additional 3.7% (7/186) of this cohort with pathogenic or likely pathogenic CNV while a further 7.7% (14/186) had CNV of uncertain significance. We identified recurrent genomic rearrangements previously associated with CP due to 2p25.3 deletion, 22q11.2 deletions and duplications and Xp monosomy. We also discovered a deletion of a single gene, PDCD6IP, and performed additional zebrafish model studies to support its single allele loss in CP aetiology. Combined SNV and CNV analysis revealed pathogenic and likely pathogenic variants in 22.7% of unselected individuals with CP.

Analysis of 182 cerebral palsy transcriptomes points to dysregulation of trophic signalling pathways and overlap with autism.

Cerebral palsy (CP) is the most common motor disability of childhood. It is characterised by permanent, non-progressive but not unchanging problems with movement, posture and motor function, with a highly heterogeneous clinical spectrum and frequent neurodevelopmental comorbidities. The aetiology of CP is poorly understood, despite recent reports of a genetic contribution in some cases. Here we demonstrate transcriptional dysregulation of trophic signalling pathways in patient-derived cell lines from an unselected cohort of 182 CP-affected individuals using both differential expression analysis and weighted gene co-expression network analysis (WGCNA). We also show that genes differentially expressed in CP, as well as network modules significantly correlated with CP status, are enriched for genes associated with ASD. Combining transcriptome and whole exome sequencing (WES) data for this CP cohort likely resolves an additional 5% of cases separated to the 14% we have previously reported as resolved by WES. Collectively, these results support a convergent molecular abnormality in CP and ASD.