PURPOSE:The gold standard for identification of post-zygotic variants (PZVs) is droplet digital polymerase chain reaction 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 (n = 145 trios) and Simons Simplex Collection (n = 405 families). Candidate mosaic variants were validated using deep amplicon sequencing or droplet digital polymerase chain reaction. 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 Australian Cerebral Palsy Biobank, 16.6% of probands and 20.7% of parents had at least 1 true-positive somatic or pGoM variant, respectively. A large proportion of PZVs detected in Simons Simplex Collection parents (79.8%) and child (94.5%) were not previously reported. We reclassified 3.7% to 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.

We implicated the X-chromosome THOC2 gene, which encodes the largest subunit of the highly-conserved TREX ( Tr anscription- Ex port) 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.

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.

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.

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.

Objective To identify the causative gene in a large unsolved family with genetic epilepsy with febrile seizures plus (GEFS+), we sequenced the genomes of family members, and then determined the contribution of the identified gene to the pathogenicity of epilepsies by examining sequencing data from 2,772 additional patients. Methods We performed whole genome sequencing of 3 members of a GEFS+ family. Subsequently, whole exome sequencing data from 1,165 patients with epilepsy from the Epi4K dataset and 1,329 Australian patients with epilepsy from the Epi25 dataset were interrogated. Targeted resequencing was performed on 278 patients with febrile seizures or GEFS+ phenotypes. Variants were validated and familial segregation examined by Sanger sequencing. Results Eight previously unreported missense variants were identified in SLC32A1, coding for the vesicular inhibitory amino acid cotransporter VGAT. Two variants cosegregated with the phenotype in 2 large GEFS+ families containing 8 and 10 affected individuals, respectively. Six further variants were identified in smaller families with GEFS+ or idiopathic generalized epilepsy (IGE). Conclusion Missense variants in SLC32A1 cause GEFS+ and IGE. These variants are predicted to alter gamma-aminobutyric acid (GABA) transport into synaptic vesicles, leading to altered neuronal inhibition. Examination of further epilepsy cohorts will determine the full genotype-phenotype spectrum associated with SLC32A1 variants.

PCDH19 is a nonclustered protocadherin molecule involved in axon bundling, synapse function, and transcriptional coregulation. Pathogenic variants in PCDH19 cause infantile-onset epilepsy known as PCDH19-clustering epilepsy or PCDH19-CE. Recent advances in DNA-sequencing technologies have led to a significant increase in the number of reported PCDH19-CE variants, many of uncertain significance. We aimed to determine the best approaches for assessing the disease relevance of missense variants in PCDH19. The application of the American College of Medical Genetics and Association for Molecular Pathology (ACMG-AMP) guidelines was only 50% accurate. Using a training set of 322 known benign or pathogenic missense variants, we identified MutPred2, MutationAssessor, and GPP as the best performing in silico tools. We generated a protein structural model of the extracellular domain and assessed 24 missense variants. We also assessed 24 variants using an in vitro reporter assay. A combination of these tools was 93% accurate in assessing known pathogenic and benign PCDH19 variants. We increased the accuracy of the ACMG-AMP classification of 45 PCDH19 variants from 50% to 94%, using these tools. In summary, we have developed a robust toolbox for the assessment of PCDH19 variant pathogenicity to improve the accuracy of PCDH19-CE variant classification.

Cerebral palsy (CP) is the most common cause of childhood physical disability, with incidence between 1/500 and 1/700 births in the developed world. Despite increasing evidence for a major contribution of genetics to CP aetiology, genetic testing is currently not performed systematically. We assessed the diagnostic rate of genome sequencing (GS) in a clinically unselected cohort of 150 singleton CP patients, with CP confirmed at >4 years of age. Clinical grade GS was performed on the proband and variants were filtered, and classified according to American College of Medical Genetics and Genomics–Association for Molecular Pathology (ACMG-AMP) guidelines. Variants classified as pathogenic or likely pathogenic (P/LP) were further assessed for their contribution to CP. In total, 24.7% of individuals carried a P/LP variant(s) causing or increasing risk of CP, with 4.7% resolved by copy number variant analysis and 20% carrying single nucleotide or indel variants. A further 34.7% carried one or more rare, high impact variants of uncertain significance (VUS) in variation intolerant genes. Variants were identified in a heterogeneous group of genes, including genes associated with hereditary spastic paraplegia, clotting and thrombophilic disorders, small vessel disease, and other neurodevelopmental disorders. Approximately 1/2 of individuals were classified as likely to benefit from changed clinical management as a result of genetic findings. In addition, no significant association between genetic findings and clinical factors was detectable in this cohort, suggesting that systematic sequencing of CP will be required to avoid missed diagnoses.

The pioneering discovery research of X-linked intellectual disability (XLID) genes has benefitted thousands of individuals worldwide; however, approximately 30% of XLID families still remain unresolved. We postulated that noncoding variants that affect gene regulation or splicing may account for the lack of a genetic diagnosis in some cases. Detecting pathogenic, gene-regulatory variants with the same sensitivity and specificity as structural and coding variants is a major challenge for Mendelian disorders. Here, we describe three pedigrees with suggestive XLID where distinctive phenotypes associated with known genes guided the identification of three different noncoding variants. We used comprehensive structural, single-nucleotide, and repeat expansion analyses of genome sequencing. RNA-Seq from patient-derived cell lines, reverse-transcription polymerase chain reactions, Western blots, and reporter gene assays were used to confirm the functional effect of three fundamentally different classes of pathogenic noncoding variants: a retrotransposon insertion, a novel intronic splice donor, and a canonical splice variant of an untranslated exon. In one family, we excluded a rare coding variant in ARX, a known XLID gene, in favor of a regulatory noncoding variant in OFD1 that correlated with the clinical phenotype. Our results underscore the value of genomic research on unresolved XLID families to aid novel, pathogenic noncoding variant discovery.

USP9X is an X-chromosome gene that escapes X-inactivation. Loss or compromised function of USP9X leads to neurodevelopmental disorders in males and females. While males are impacted primarily by hemizygous partial loss-of-function missense variants, in females de novo heterozygous complete loss-of-function mutations predominate, and give rise to the clinically recognisable USP9X-female syndrome. Here we provide evidence of the contribution of USP9X missense and small in-frame deletion variants in USP9X-female syndrome also. We scrutinise the pathogenicity of eleven such variants, ten of which were novel. Combined application of variant prediction algorithms, protein structure modelling, and assessment under clinically relevant guidelines universally support their pathogenicity. The core phenotype of this cohort overlapped with previous descriptions of USP9X-female syndrome, but exposed heightened variability. Aggregate phenotypic information of 35 currently known females with predicted pathogenic variation in USP9X reaffirms the clinically recognisable USP9X-female syndrome, and highlights major differences when compared to USP9X-male associated neurodevelopmental disorders.

The major and most well-studied genetic cause of Fragile-X syndrome (FXS) is expansion of a CGG repeat in the 5'-UTR of the FMR1 gene. Routine testing for this expansion is performed globally. Overall, there is a paucity of intragenic variants explaining FXS, a fact which is being addressed by a more systematic application of whole exome (WES) and whole genome (WGS) sequencing, even in the diagnostic setting. Here we report two families comprising probands with a clinical suspicion of FXS and no CGG repeat expansions. Using WES/WGS we identified deleterious variants within the coding region of FMR1 in both families. In a family from Finland we identified a complex indel c.1021-1028delinsTATTGG in exon 11 of FMR1 which gives rise to a frameshift and a premature termination codon (PTC), p.Asn341Tyrfs*7. Follow-up mRNA and protein studies on a cell line from the proband revealed that although the mRNA levels of FMR1 were not altered, Fragile X Mental Retardation 1 Protein (FMRP) was undetectable. Additionally, we identified a variant, c.881-1G > T, affecting the canonical acceptor splice site of exon 10 of FMR1 in an Australian family. Our findings reinforce the importance of intragenic FMR1 variant testing, particularly in cases with clinical features of FXS and no CGG repeat expansions identified.

Multiple TREX mRNA export complex subunits (e.g., THOC1, THOC2, THOC5, THOC6, THOC7) have now been implicated in neurodevelopmental disorders (NDDs), neurodegeneration and cancer. We previously implicated missense and splicing-defective THOC2 variants in NDDs and a broad range of other clinical features. Here we report 10 individuals from nine families with rare missense THOC2 variants including the first case of a recurrent variant (p.Arg77Cys), and an additional individual with an intragenic THOC2 microdeletion (Del-Ex37-38). Ex vivo missense variant testing and patient-derived cell line data from current and published studies show 9 of the 14 missense THOC2 variants result in reduced protein stability. The splicing-defective and deletion variants result in a loss of small regions of the C-terminal THOC2 RNA binding domain (RBD). Interestingly, reduced stability of THOC2 variant proteins has a flow-on effect on the stability of the multi-protein TREX complex; specifically on the other NDD-associated THOC subunits. Our current, expanded cohort refines the core phenotype of THOC2 NDDs to language disorder and/or ID, with a variable severity, and disorders of growth. A subset of affected individuals' has severe-profound ID, persistent hypotonia and respiratory abnormalities. Further investigations to elucidate the pathophysiological basis for this severe phenotype are warranted.

The need to interpret the pathogenicity of novel missense variants of unknown significance identified in the homeodomain of X-chromosome aristaless-related homeobox (ARX) gene prompted us to assess the utility of conservation and constraint across these domains in multiple genes compared to conventional in vitro functional analysis. Pathogenic missense variants clustered in the homeodomain of ARX contribute to intellectual disability (ID) and epilepsy, with and without brain malformation in affected males. Here we report novel c.1112G>A, p.Arg371Gln and c.1150C>T, p.Arg384Cys variants in male patients with ID and severe seizures. The third case of a male patient with a c.1109C>T, p.Ala370Val variant is perhaps the first example of ID and autism spectrum disorder (ASD), without seizures or brain malformation. We compiled data sets of pathogenic variants from ClinVar and presumed benign variation from gnomAD and demonstrated that the high levels of sequence conservation and constraint of benign variation within the homeodomain impacts upon the ability of publicly available in silico prediction tools to accurately discern likely benign from likely pathogenic variants in these data sets. Despite this, considering the inheritance patterns of the genes and disease variants with the conservation and constraint of disease variants affecting the homeodomain in conjunction with current clinical assessments may assist in predicting the pathogenicity of missense variants, particularly for genes with autosomal recessive and X-linked patterns of disease inheritance, such as ARX. In vitro functional analysis demonstrates that the transcriptional activity of all three variants was diminished compared to ARX-Wt. We review the associated phenotypes of the published cases of patients with ARX homeodomain variants and propose expansion of the ARX-related phenotype to include severe ID and ASD without brain malformations or seizures. We propose that the use of the constraint and conservation data in conjunction with consideration of the patient phenotype and inheritance pattern may negate the need for the experimental functional validation currently required to achieve a diagnosis.

BACKGROUND: The X-chromosome gene USP9X encodes a deubiquitylating enzyme that has been associated with neurodevelopmental disorders primarily in female subjects. USP9X escapes X inactivation, and in female subjects de novo heterozygous copy number loss or truncating mutations cause haploinsufficiency culminating in a recognizable syndrome with intellectual disability and signature brain and congenital abnormalities. In contrast, the involvement of USP9X in male neurodevelopmental disorders remains tentative. METHODS: We used clinically recommended guidelines to collect and interrogate the pathogenicity of 44 USP9X variants associated with neurodevelopmental disorders in males. Functional studies in patient-derived cell lines and mice were used to determine mechanisms of pathology. RESULTS: Twelve missense variants showed strong evidence of pathogenicity. We define a characteristic phenotype of the central nervous system (white matter disturbances, thin corpus callosum, and widened ventricles); global delay with significant alteration of speech, language, and behavior; hypotonia; joint hypermobility; visual system defects; and other common congenital and dysmorphic features. Comparison of in silico and phenotypical features align additional variants of unknown significance with likely pathogenicity. In support of partial loss-of-function mechanisms, using patient-derived cell lines, we show loss of only specific USP9X substrates that regulate neurodevelopmental signaling pathways and a united defect in transforming growth factor signaling. In addition, we find correlates of the male phenotype in Usp9x brain-specific knockout mice, and further resolve loss of hippocannpal-dependent learning and memory. CONCLUSIONS: Our data demonstrate the involvement of USP9X variants in a distinctive neurodevelopmental and behavioral syndrome in male subjects and identify plausible mechanisms of pathogenesis centered on disrupted transforming growth factor beta signaling and hippocampal function.

We show that the loss of Iqsec2 function in mice recapitulates key aspects of the human phenotype, irrespective of the X-inactivation status of the gene between species. Our understanding of the traditional X-chromosome inheritance with heterozygous female sparing needs revisiting. Clinical presentations of mutations in the IQSEC2 gene on the X-chromosome initially implicated to cause non-syndromic intellectual disability (ID) in males have expanded to include early onset seizures in males as well as in females. The molecular pathogenesis is not well understood, nor the mechanisms driving disease expression in heterozygous females. Using a CRISPR/Cas9–edited Iqsec2 KO mouse model, we confirm the loss of Iqsec2 mRNA expression and lack of Iqsec2 protein within the brain of both founder and progeny mice. Both male (52%) and female (46%) Iqsec2 KO mice present with frequent and recurrent seizures. Focusing on Iqsec2 KO heterozygous female mice, we demonstrate increased hyperactivity, altered anxiety and fear responses, decreased social interactions, delayed learning capacity and decreased memory retention/novel recognition, recapitulating psychiatric issues, autistic-like features, and cognitive deficits present in female patients with loss-of-function IQSEC2 variants. Despite Iqsec2 normally acting to activate Arf6 substrate, we demonstrate that mice modelling the loss of Iqsec2 function present with increased levels of activated Arf6. We contend that loss of Iqsec2 function leads to altered regulation of activated Arf6-mediated responses to synaptic signalling and immature synaptic networks. We highlight the importance of IQSEC2 function for females by reporting a novel nonsense variant c.566C > A, p.(S189*) in an elderly female patient with profound intellectual disability, generalised seizures, and behavioural disturbances. Our human and mouse data reaffirm IQSEC2 as another disease gene with an unexpected X-chromosome heterozygous female phenotype. Our Iqsec2 mouse model recapitulates the phenotypes observed in human patients despite the differences in the IQSEC2/Iqsec2 gene X-chromosome inactivation between the species.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Mutations in ATP6V1B2, which encodes the B2 subunit of the vacuolar H + ATPase have previously been associated with Zimmermann-Laband syndrome 2 (ZLS2) and deafness-onychodystrophy (DDOD) syndrome. Recently epilepsy has also been described as a potentially associated phenotype. Here we further uncover the role of ATP61VB2 in epilepsy and report autosomal dominant inheritance of a novel missense variant in ATP6V1B2 in a large Polish family with relatively mild gingival and nail problems, no phalangeal hypoplasia and with generalized epilepsy. In light of our findings and review of the literature, we propose that the ATP6V1B2 gene should be considered in families with autosomal dominant epilepsy both with or without intellectual disability, and that presence of subtle gingival and nail problems may be another characteristic calling card of affected individuals with ATP6V1B2 mutations.