Studying Disease-Associated UBE3A Missense Variants Using Enhanced Sampling Molecular Simulations.

Missense variants in UBE3A underlie neurodevelopmental conditions such as Angelman Syndrome and Autism Spectrum Disorder, but the underlying molecular pathological consequences on protein folding and function are poorly understood. Here, we report a novel, maternally inherited, likely pathogenic missense variant in UBE3A (NM_000462.4(UBE3A_v001):(c.1841T>C) (p.(Leu614Pro))) in a child that presented with myoclonic epilepsy from 14 months, subsequent developmental regression from 16 months, and additional features consistent with Angelman Syndrome. To understand the impact of p.(Leu614Pro) on UBE3A, we used adiabatic biased molecular dynamics and metadynamics simulations to investigate conformational differences from wildtype proteins. Our results suggest that Leu614Pro substitution leads to less efficient binding and substrate processing compared to wildtype. Our results support the use of enhanced sampling molecular simulations to investigate the impact of missense UBE3A variants on protein function that underlies neurodevelopment and human disorders.

Epilepsy surgery in PCDH 19 related developmental and epileptic encephalopathy: A case report.

We report a female child with PCDH19 related developmental and epileptic encephalopathy with drug-resistant seizures, cognitive and language impairment, autism spectrum disorder and sleep dysfunction. Her seizures, which started at 10 months of age, were resistant to multiple anti-seizure medications. Developmental stagnation followed by regression occurred after the onset of recurrent seizures. Her ictal EEGS suggested left temporal lobe origin for her recorded seizures. MRI upon expert re-review showed a subtle abnormality in the left temporal lobe. In view of the severe nature and frequency of her seizures, a left temporal lobectomy was undertaken at the age of 2 years and 3 months. Though her seizure outcome was Engel class 3, her seizure frequency and severity were significantly reduced. She has been seizure-free for 10 months at her last outpatient assessment when she was 4 years and 8 months of age (2 years and 5 months after epilepsy surgery). However she recently had an admission for COVID19 infection, with a breakthrough cluster of seizures. Her developmental trajectory changed, though she is making good progress with her cognitive and language skills.

Pathogenic variants in RNPC3 are associated with hypopituitarism and primary ovarian insufficiency.

PURPOSE: We aimed to investigate the molecular basis underlying a novel phenotype including hypopituitarism associated with primary ovarian insufficiency. METHODS: We used next-generation sequencing to identify variants in all pedigrees. Expression of Rnpc3/RNPC3 was analyzed by in situ hybridization on murine/human embryonic sections. CRISPR/Cas9 was used to generate mice carrying the p.Leu483Phe pathogenic variant in the conserved murine Rnpc3 RRM2 domain. RESULTS: We described 15 patients from 9 pedigrees with biallelic pathogenic variants in RNPC3, encoding a specific protein component of the minor spliceosome, which is associated with a hypopituitary phenotype, including severe growth hormone (GH) deficiency, hypoprolactinemia, variable thyrotropin (also known as thyroid-stimulating hormone) deficiency, and anterior pituitary hypoplasia. Primary ovarian insufficiency was diagnosed in 8 of 9 affected females, whereas males had normal gonadal function. In addition, 2 affected males displayed normal growth when off GH treatment despite severe biochemical GH deficiency. In both mouse and human embryos, Rnpc3/RNPC3 was expressed in the developing forebrain, including the hypothalamus and Rathke's pouch. Female Rnpc3 mutant mice displayed a reduction in pituitary GH content but with no reproductive impairment in young mice. Male mice exhibited no obvious phenotype. CONCLUSION: Our findings suggest novel insights into the role of RNPC3 in female-specific gonadal function and emphasize a critical role for the minor spliceosome in pituitary and ovarian development and function.

Further heterogeneity in Silver-Russell syndrome: PLAG1 deletion in association with a complex chromosomal rearrangement.

Silver-Russell syndrome (SRS) is a rare genetic condition primarily characterized by growth restriction and facial dysmorphisms. While hypomethylation of H19/IGF2:IG-DMR (imprinting control region 1 [IC1]) located at 11p15.5 and maternal uniparental disomy of chromosome 7 (upd[7]mat) are the most common genetic mechanisms responsible for SRS, the expanding body of literature describing alternative causative variants suggests SRS is a highly heterogeneous condition, also involving variation in the HMGA2-PLAG1-IGF2 pathway. We report a familial PLAG1 deletion in association with a complex chromosomal rearrangement. We describe two siblings with differing unbalanced chromosomal rearrangements inherited from a mother with a 5-breakpoint balanced complex rearrangement involving chromosomes 2, 8, and 21. The overlapping but diverse phenotypes in the siblings were characterized by shared SRS-like features, underlined by a PLAG1 whole gene deletion. Genetic analysis and interpretation was further complicated by a meiotic recombination event occurring in one of the siblings. This family adds to the limited literature available on PLAG1-related SRS. We have reviewed all currently known cases aiming to define the associated phenotype and guide future genetic testing strategies. The heterogeneity of SRS is further expanded by the involvement of complex cytogenomic abnormalities, imposing requirements for a comprehensive approach to testing and genetic counseling.

Deletion of ERF and CIC causes abnormal skull morphology and global developmental delay.

The ETS2 repressor factor (ERF) is a transcription factor in the RAS-MEK-ERK signal transduction cascade that regulates cell proliferation and differentiation, and pathogenic sequence variants in the ERF gene cause variable craniosynostosis inherited in an autosomal dominant pattern. The reported ERF variants are largely loss-of-function, implying haploinsufficiency as a primary disease mechanism; however, ERF gene deletions have not been reported previously. Here we describe three probands with macrocephaly, craniofacial dysmorphology, and global developmental delay. Clinical genetic testing for fragile X and other relevant sequencing panels were negative; however, chromosomal microarray identified heterozygous deletions (63.7-583.2 kb) on Chromosome 19q13.2 in each proband that together included five genes associated with Mendelian diseases (ATP1A3, ERF, CIC, MEGF8, and LIPE). Parental testing indicated that the aberrations were apparently de novo in two of the probands and were inherited in the one proband with the smallest deletion. Deletion of ERF is consistent with the reported loss-of-function ERF variants, prompting clinical copy-number-variant classifications of likely pathogenic. Moreover, the recent characterization of heterozygous loss-of-function CIC sequence variants as a cause of intellectual disability and neurodevelopmental disorders inherited in an autosomal dominant pattern is also consistent with the developmental delays and intellectual disabilities identified among the two probands with CIC deletions. Taken together, this case series adds to the previously reported patients with ERF and/or CIC sequence variants and supports haploinsufficiency of both genes as a mechanism for a variable syndromic cranial phenotype with developmental delays and intellectual disability inherited in an autosomal dominant pattern.

Expanding Clinical Presentations Due to Variations in THOC2 mRNA Nuclear Export Factor.

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.

A homozygous UBA5 pathogenic variant causes a fatal congenital neuropathy.

BACKGROUND: UBA5 is the activating enzyme of UFM1 in the ufmylation post-translational modification system. Different neurological phenotypes have been associated with UBA5 pathogenic variants including epilepsy, intellectual disability, movement disorders and ataxia. METHODS AND RESULTS: We describe a large multigenerational consanguineous family presenting with a severe congenital neuropathy causing early death in infancy. Whole exome sequencing and linkage analysis identified a novel homozygous UBA5 NM_024818.3 c.31C>T (p.Arg11Trp) mutation. Protein expression assays in mouse tissue showed similar levels of UBA5 in peripheral nerves to the central nervous system. CRISPR-Cas9 edited HEK (human embrionic kidney) cells homozygous for the UBA5 p.Arg11Trp mutation showed reduced levels of UBA5 protein compared with the wild-type. The mutant p.Arg11Trp UBA5 protein shows reduced ability to activate UFM1. CONCLUSION: This report expands the phenotypical spectrum of UBA5 mutations to include fatal peripheral neuropathy.

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.].

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.

Silver Russel syndrome in an aboriginal patient from Australia.

Silver-Russell syndrome (SRS OMIM 180860) is a rare, albeit well-recognized disorder characterized by severe intrauterine and postnatal growth retardation. It remains a clinical diagnosis with a molecular cause identifiable in approximately 60%-70% of patients. We report a 4-year-old Australian Aboriginal girl who was born at 32 weeks gestation with features strongly suggestive of SRS, after extensive investigation she was referred to our undiagnosed disease program (UDP). Genomic sequencing was performed which identified a heterozygous splice site variant in IGF2 which is predicted to be pathogenic by in-silico studies, paternal allelic origin, de novo status, and RNA studies on fibroblasts. We compare clinical findings with reported patients to add to the knowledge base on IGF2 variants and to promote the engagement of other Australian Aboriginal families in genomic medicine.

Angiopoietin-1 is required for Schlemm's canal development in mice and humans.

Primary congenital glaucoma (PCG) is a leading cause of blindness in children worldwide and is caused by developmental defects in 2 aqueous humor outflow structures, Schlemm's canal (SC) and the trabecular meshwork. We previously identified loss-of-function mutations in the angiopoietin (ANGPT) receptor TEK in families with PCG and showed that ANGPT/TEK signaling is essential for SC development. Here, we describe roles for the major ANGPT ligands in the development of the aqueous outflow pathway. We determined that ANGPT1 is essential for SC development, and that Angpt1-knockout mice form a severely hypomorphic canal with elevated intraocular pressure. By contrast, ANGPT2 was dispensable, although mice deficient in both Angpt1 and Angpt2 completely lacked SC, indicating that ANGPT2 compensates for the loss of ANGPT1. In addition, we identified 3 human subjects with rare ANGPT1 variants within an international cohort of 284 PCG patients. Loss of function in 2 of the 3 patient alleles was observed by functional analysis of ANGPT1 variants in a combined in silico, in vitro, and in vivo approach, supporting a causative role for ANGPT1 in disease. By linking ANGPT1 with PCG, these results highlight the importance of ANGPT/TEK signaling in glaucoma pathogenesis and identify a candidate target for therapeutic development.

Transcriptome-wide effects of a POLR3A gene mutation in patients with an unusual phenotype of striatal involvement.

RNA polymerase III is essential for the transcription of non-coding RNAs, including tRNAs. Mutations in the genes encoding its largest subunits are known to cause hypomyelinating leukodystrophies (HLD7) with pathogenetic mechanisms hypothesised to involve impaired availability of tRNAs. We have identified a founder mutation in the POLR3A gene that leads to aberrant splicing, a premature termination codon and partial deficiency of the canonical full-length transcript. Our clinical and imaging data showed no evidence of the previously reported white matter or cerebellar involvement; instead the affected brain structures included the striatum and red nuclei with the ensuing clinical manifestations. Our transcriptome-wide investigations revealed an overall decrease in the levels of Pol III-transcribed tRNAs and an imbalance in the levels of regulatory ncRNAs such as small nuclear and nucleolar RNAs (snRNAs and snoRNAs). In addition, the Pol III mutation was found to exert complex downstream effects on the Pol II transcriptome, affecting the general regulation of RNA metabolism.

Angiopoietin receptor TEK mutations underlie primary congenital glaucoma with variable expressivity.

Primary congenital glaucoma (PCG) is a devastating eye disease and an important cause of childhood blindness worldwide. In PCG, defects in the anterior chamber aqueous humor outflow structures of the eye result in elevated intraocular pressure (IOP); however, the genes and molecular mechanisms involved in the etiology of these defects have not been fully characterized. Previously, we observed PCG-like phenotypes in transgenic mice that lack functional angiopoietin-TEK signaling. Herein, we identified rare TEK variants in 10 of 189 unrelated PCG families and demonstrated that each mutation results in haploinsufficiency due to protein loss of function. Multiple cellular mechanisms were responsible for the loss of protein function resulting from individual TEK variants, including an absence of normal protein production, protein aggregate formation, enhanced proteasomal degradation, altered subcellular localization, and reduced responsiveness to ligand stimulation. Further, in mice, hemizygosity for Tek led to the formation of severely hypomorphic Schlemm's canal and trabecular meshwork, as well as elevated IOP, demonstrating that anterior chamber vascular development is sensitive to Tek gene dosage and the resulting decrease in angiopoietin-TEK signaling. Collectively, these results identify TEK mutations in patients with PCG that likely underlie disease and are transmitted in an autosomal dominant pattern with variable expressivity.

Mosaic dominant TUBB4A mutation in an inbred family with complicated hereditary spastic paraplegia.

BACKGROUND: Mutations in TUBB4A have been associated with a spectrum of neurological conditions, ranging from the severe hypomyelination with atrophy of the basal ganglia and cerebellum syndrome to the clinically milder dystonia type 4. The presence of movement abnormalities was considered the common hallmark of these disorders. METHODS: Clinical, neurological, and neuroimaging examinations, followed by whole exome sequencing and mutation analysis, were performed in a highly consanguineous pedigree with five affected children. RESULTS: We identified a novel c.568C>T (p.H190Y) TUBB4A mutation that originated de novo in the asymptomatic mother. The affected subjects presented with an early-onset, slowly progressive spastic paraparesis of the lower limbs, ataxia, and brain hypomyelination, in the absence of dystonia or rigidity. CONCLUSIONS: Our study adds complicated hereditary spastic paraplegia to the clinical spectrum of TUBB4A-associated neurological disorders. We establish genotype-phenotype correlations with mutations located in the same region in the tertiary structure of the protein.

Founder p.Arg 446* mutation in the PDHX gene explains over half of cases with congenital lactic acidosis in Roma children.

Investigation of 31 of Roma patients with congenital lactic acidosis (CLA) from Bulgaria identified homozygosity for the R446* mutation in the PDHX gene as the most common cause of the disorder in this ethnic group. It accounted for around 60% of patients in the study and over 25% of all CLA cases referred to the National Genetic Laboratory in Bulgaria. The detection of a homozygous patient from Hungary and carriers among population controls from Romania and Slovakia suggests a wide spread of the mutation in the European Roma population. The clinical phenotype of the twenty R446* homozygotes was relatively homogeneous, with lactic acidosis crisis in the first days or months of life as the most common initial presentation (15/20 patients) and delayed psychomotor development and/or seizures in infancy as the leading manifestations in a smaller group (5/20 patients). The subsequent clinical picture was dominated by impaired physical growth and a very consistent pattern of static cerebral palsy-like encephalopathy with spasticity and severe to profound mental retardation seen in over 80% of cases. Most patients had a positive family history. We propose testing for the R446* mutation in PDHX as a rapid first screening in Roma infants with metabolic acidosis. It will facilitate and accelerate diagnosis in a large proportion of cases, allow early rehabilitation to alleviate the chronic clinical course, and prevent further affected births in high-risk families.

Challenges of diagnostic exome sequencing in an inbred founder population.

Exome sequencing was used as a diagnostic tool in a Roma/Gypsy family with three subjects (one deceased) affected by lissencephaly with cerebellar hypoplasia (LCH), a clinically and genetically heterogeneous diagnostic category. Data analysis identified high levels of unreported inbreeding, with multiple rare/novel "deleterious" variants occurring in the homozygous state in the affected individuals. Step-wise filtering was facilitated by the inclusion of parental samples in the analysis and the availability of ethnically matched control exome data. We identified a novel mutation, p.Asp487Tyr, in the VLDLR gene involved in the Reelin developmental pathway and associated with a rare form of LCH, the Dysequilibrium Syndrome. p.Asp487Tyr is the third reported missense mutation in this gene and the first example of a change affecting directly the functionally crucial ß-propeller domain. An unexpected additional finding was a second unique mutation (p.Asn494His) with high scores of predicted pathogenicity in KCNV2, a gene implicated in a rare eye disorder, retinal cone dystrophy type 3B. This result raised diagnostic and counseling challenges that could be resolved through mutation screening of a large panel of healthy population controls. The strategy and findings of this study may inform the search for new disease mutations in the largest European genetic isolate.

Autosomal-recessive congenital cerebellar ataxia is caused by mutations in metabotropic glutamate receptor 1.

Autosomal-recessive congenital cerebellar ataxia was identified in Roma patients originating from a small subisolate with a known strong founder effect. Patients presented with global developmental delay, moderate to severe stance and gait ataxia, dysarthria, mild dysdiadochokinesia, dysmetria and tremors, intellectual deficit, and mild pyramidal signs. Brain imaging revealed progressive generalized cerebellar atrophy, and inferior vermian hypoplasia and/or a constitutionally small brain were observed in some patients. Exome sequencing, used for linkage analysis on extracted SNP genotypes and for mutation detection, identified two novel (i.e., not found in any database) variants located 7 bp apart within a unique 6q24 linkage region. Both mutations cosegregated with the disease in five affected families, in which all ten patients were homozygous. The mutated gene, GRM1, encodes metabotropic glutamate receptor mGluR1, which is highly expressed in cerebellar Purkinje cells and plays an important role in cerebellar development and synaptic plasticity. The two mutations affect a gene region critical for alternative splicing and the generation of receptor isoforms; they are a 3 bp exon 8 deletion and an intron 8 splicing mutation (c.2652_2654del and c.2660+2T>G, respectively [RefSeq accession number NM_000838.3]). The functional impact of the deletion is unclear and is overshadowed by the splicing defect. Although ataxia lymphoblastoid cell lines expressed GRM1 at levels comparable to those of control cells, the aberrant transcripts skipped exon 8 or ended in intron 8 and encoded various species of nonfunctional receptors either lacking the transmembrane domain and containing abnormal intracellular tails or completely missing the tail. The study implicates mGluR1 in human hereditary ataxia. It also illustrates the potential of the Roma founder populations for mutation identification by exome sequencing.

A novel GEFS+ locus on 12p13.33 in a large Roma family.

We report a clinical and genetic follow-up study of a large consanguineous family from an endogamous Roma/Gypsy sub-isolate, where previous analyses have been inconclusive. Detailed clinical information was collected through extensive field work, repeat interviews and electrophysiological and neuroimaging investigations on 18 affected subjects. The phenotype is compatible with GEFS+, with some unusual features, e.g. GTCS persisting into late adult life and high frequency of focal epilepsy. Updated genealogical information, a dense SNP genome scan and linkage analysis identified a novel GEFS+ locus on 12p13.33, where 13 affected individuals from two branches of the kindred shared an identical haplotype. This haplotype was not found in the 3rd branch or in the remaining 21 Roma epilepsy families in our collection. Genetic heterogeneity and evidence of bilineality were found despite the inbreeding and endogamous nature of the family and population of origin. These data add to the growing evidence of lack of founder effect and significant genetic heterogeneity in epilepsy in the Roma/Gypsy population. Sequencing of the coding regions of three genes linked to neurotransmitter transport and release, SLC6A12, SLC6A13 and ERC1, on 12p did not identify a causative mutation.

Identification of a Gypsy SHOX mutation (p.A170P) in Léri-Weill dyschondrosteosis and Langer mesomelic dysplasia.

We report the clinical and molecular characteristics of 12 Spanish families with multiple members affected with Léri-Weill dyschondrosteosis (LWD) or Langer mesomelic dysplasia (LMD), who present the SHOX (short stature homeobox gene) mutation p.A170P (c.508G>C) in heterozygosity or homozygosity, respectively. In all studied families, the A170P mutation co-segregated with the fully penetrant phenotype of mesomelic limb shortening and Madelung deformity. A shared haplotype around SHOX was observed by microsatellite analysis, confirming the presence of a common ancestor, probably of Gypsy origin, as 11 of the families were of this ethnic group. Mutation screening in 359 Eastern-European Gypsies failed to identify any carriers. For the first time, we have shown SHOX expression in the human growth plate of a 22-week LMD fetus, homozygous for the A170P mutation. Although the mutant SHOX protein was expressed in all zones of the growth plate, the chondrocyte columns in the proliferative zone were disorganized with the chondrocytes occurring in smaller columnal clusters. We have also identified a novel mutation at the same residue, c. 509C>A (p.A170D), in two unrelated Spanish LWD families, which similar to A170P mutation impedes nuclear localization of SHOX. In conclusion, we have identified A170P as the first frequent SHOX mutation in Gypsy LWD and LMD individuals.