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.

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.

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.

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.

Focal epilepsy of probable temporal lobe origin in a Gypsy family showing linkage to a novel locus on 7p21.3.

We aimed to characterise the phenotype and perform genetic studies in a family of Roma/Gypsy ethnicity, affected by epilepsy. The mean age at onset of epilepsy was 9 years and seizures persisted into adulthood. Antecedent febrile convulsions were rare. Seizure semiology and EEG findings suggested mesial temporal lobe origin with no evidence of hippocampal sclerosis. Seizures frequently generalised. Family structure suggested autosomal-dominant inheritance with incomplete penetrance. Linkage analysis identified a single novel locus on 7p21.3, corresponding to the expected maximum in the family. Previously reported temporal lobe epilepsy (TLE) loci were definitely excluded. The minimal shared haplotype of 2.4cM (1.3Mb) was not observed in other affected families or controls from the same population. Three brain-expressed validated genes in the critical region represent potential candidates. We have identified an epilepsy syndrome with temporal lobe seizures commonly evolving to generalised convulsions. Linkage to 7p21.3 adds up to a total of five currently known FTLE loci.

LTBP2 and CYP1B1 mutations and associated ocular phenotypes in the Roma/Gypsy founder population.

Primary congenital glaucoma (PCG) is a genetically heterogeneous autosomal recessive disorder, which is an important cause of blindness in childhood. The first known gene, CYP1B1, accounts for a variable proportion of cases in most populations. A second gene, LTBP2, was recently reported in association with a syndrome, in which glaucoma is secondary to lens dislocation. We report on the molecular and clinical profile of 34 families diagnosed as PCG, all originating from the Roma/Gypsy founder population. Comprehensive sequencing analysis revealed a level of heterogeneity unusual for this population, with five CYP1B1 and one ancestral LTBP2 mutation accounting for ∼70% of patients (25 out of 37) and the remainder still unexplained. Homozygosity for the founder LTBP2 p.R299X mutation resulted in a more severe clinical phenotype and poorer outcome despite a markedly higher number of surgical interventions. The genetically homogeneous group of p.R299X homozygotes showed variable phenotypes (presumably also underlying pathogenetic mechanisms), wherein PCG proper with primary dysgenesis of the trabecular meshwork, and Marfan syndrome-like zonular disease with ectopia lentis and later onset secondary glaucoma are two extremes. The spectrum manifestations may occur in different combinations and have a different evolution even within the same sibship or a single patient. Preliminary observations on compounds with mutations in both CYP1B1-LTBP2 suggest that the observed combinations are of no clinical significance and digenic inheritance is unlikely. We provide a population genetics perspective to explain the allelic heterogeneity, comparing the history and geographic distribution of the two major founder mutations--p.R299X/LTBP2 and p.E387K/CYP1B1.

Mosaicism of a missense SCN1A mutation and Dravet syndrome in a Roma/Gypsy family.

SCN1A mutations account for a large proportion of Dravet syndrome patients, and are reported in other cases of epilepsy, such as some families with genetic epilepsy with febrile seizures plus (GEFS+). While most Dravet syndrome cases are caused by de novo mutations, 5% inherit a mutation from a mildly affected or symptom-free parent. Parental mosaicism has been identified, with documented cases involving truncating mutations or gene rearrangements. We describe a Roma/Gypsy family, where a missense mutation in SCN1A, p.D194N, is transmitted from a mosaic GEFS+ father to a child with Dravet syndrome. Mosaicism may be more common than assumed and should be considered regardless of the nature of the mutation.

Partial epilepsy syndrome in a Gypsy family linked to 5q31.3-q32.

PURPOSE: The restricted genetic diversity and homogeneous molecular basis of Mendelian disorders in isolated founder populations have rarely been explored in epilepsy research. Our long-term goal is to explore the genetic basis of epilepsies in one such population, the Gypsies. The aim of this report is the clinical and genetic characterization of a Gypsy family with a partial epilepsy syndrome. METHODS: Clinical information was collected using semistructured interviews with affected subjects and informants. At least one interictal electroencephalography (EEG) recording was performed for each patient and previous data obtained from records. Neuroimaging included structural magnetic resonance imaging (MRI). Linkage and haplotype analysis was performed using the Illumina IVb Linkage Panel, supplemented with highly informative microsatellites in linked regions and Affymetrix SNP 5.0 array data. RESULTS: We observed an early-onset partial epilepsy syndrome with seizure semiology strongly suggestive of temporal lobe epilepsy (TLE), with mild intellectual deficit co-occurring in a large proportion of the patients. Psychiatric morbidity was common in the extended pedigree but did not cosegregate with epilepsy. Linkage analysis definitively excluded previously reported loci, and identified a novel locus on 5q31.3-q32 with an logarithm of the odds (LOD) score of 3 corresponding to the expected maximum in this family. DISCUSSION: The syndrome can be classified as familial temporal lobe epilepsy (FTLE) or possibly a new syndrome with mild intellectual deficit. The linked 5q region does not contain any ion channel-encoding genes and is thus likely to contribute new knowledge about epilepsy pathogenesis. Identification of the mutation in this family and in additional patients will define the full phenotypic spectrum.