Clinical and neurogenetic characterisation of autosomal recessive RBL2-associated progressive neurodevelopmental disorder.

Retinoblastoma (RB) proteins are highly conserved transcriptional regulators that play important roles during development by regulating cell-cycle gene expression. RBL2 dysfunction has been linked to a severe neurodevelopmental disorder. However, to date, clinical features have only been described in six individuals carrying five biallelic predicted loss of function (pLOF) variants. To define the phenotypic effects of RBL2 mutations in detail, we identified and clinically characterized a cohort of 28 patients from 18 families carrying LOF variants in RBL2 , including fourteen new variants that substantially broaden the molecular spectrum. The clinical presentation of affected individuals is characterized by a range of neurological and developmental abnormalities. Global developmental delay and intellectual disability were uniformly observed, ranging from moderate to profound and involving lack of acquisition of key motor and speech milestones in most patients. Frequent features included postnatal microcephaly, infantile hypotonia, aggressive behaviour, stereotypic movements and non-specific dysmorphic features. Common neuroimaging features were cerebral atrophy, white matter volume loss, corpus callosum hypoplasia and cerebellar atrophy. In parallel, we used the fruit fly, Drosophila melanogaster , to investigate how disruption of the conserved RBL2 orthologueue Rbf impacts nervous system function and development. We found that Drosophila Rbf LOF mutants recapitulate several features of patients harboring RBL2 variants, including alterations in the head and brain morphology reminiscent of microcephaly, and perturbed locomotor behaviour. Surprisingly, in addition to its known role in controlling tissue growth during development, we find that continued Rbf expression is also required in fully differentiated post-mitotic neurons for normal locomotion in Drosophila , and that adult-stage neuronal re-expression of Rbf is sufficient to rescue Rbf mutant locomotor defects. Taken together, this study provides a clinical and experimental basis to understand genotype-phenotype correlations in an RBL2 -linked neurodevelopmental disorder and suggests that restoring RBL2 expression through gene therapy approaches may ameliorate aspects of RBL2 LOF patient symptoms.

The impact of inversions across 33,924 families with rare disease from a national genome sequencing project.

Detection of structural variants (SVs) is currently biased toward those that alter copy number. The relative contribution of inversions toward genetic disease is unclear. In this study, we analyzed genome sequencing data for 33,924 families with rare disease from the 100,000 Genomes Project. From a database hosting >500 million SVs, we focused on 351 genes where haploinsufficiency is a confirmed disease mechanism and identified 47 ultra-rare rearrangements that included an inversion (24 bp to 36.4 Mb, 20/47 de novo). Validation utilized a number of orthogonal approaches, including retrospective exome analysis. RNA-seq data supported the respective diagnoses for six participants. Phenotypic blending was apparent in four probands. Diagnostic odysseys were a common theme (>50 years for one individual), and targeted analysis for the specific gene had already been performed for 30% of these individuals but with no findings. We provide formal confirmation of a European founder origin for an intragenic MSH2 inversion. For two individuals with complex SVs involving the MECP2 mutational hotspot, ambiguous SV structures were resolved using long-read sequencing, influencing clinical interpretation. A de novo inversion of HOXD11-13 was uncovered in a family with Kantaputra-type mesomelic dysplasia. Lastly, a complex translocation disrupting APC and involving nine rearranged segments confirmed a clinical diagnosis for three family members and resolved a conundrum for a sibling with a single polyp. Overall, inversions play a small but notable role in rare disease, likely explaining the etiology in around 1/750 families across heterogeneous clinical cohorts.

Reticulon 2 deficiency results in an autosomal recessive distal motor neuropathy with lower limb spasticity.

Heterozygous RTN2 variants have been previously identified in a limited cohort of families affected by autosomal dominant spastic paraplegia (SPG12-OMIM:604805) with a variable age of onset. Nevertheless, the definitive validity of SPG12 remains to be confidently confirmed due to scarcity of supporting evidence. In our study, we identified and validated seven novel or ultra-rare homozygous loss-of-function RTN2 variants in 14 individuals from seven consanguineous families with distal hereditary motor neuropathy (dHMN) using exome, genome and Sanger sequencing coupled with deep-phenotyping. All affected individuals (seven males and seven females, aged 9-50 years) exhibited weakness in the distal upper and lower limbs, lower limb spasticity, hyperreflexia, with an onset in the first decade of life. Nerve conduction studies revealed axonal motor neuropathy with neurogenic changes in the electromyography. Despite a slowly progressive disease course, all patients remained ambulatory over a mean disease duration of 19.71 ± 13.70 years. Characterisation of C. elegans RTN2 homolog loss-of-function variants demonstrated morphological and behavioural differences compared to the parental strain. Treatment of the mutant with an endoplasmic/sarcoplasmic reticulum Ca2+ reuptake inhibitor (2,5-di-tert-butylhydroquinone) rescued key phenotypic differences, suggesting a potential therapeutic benefit for RTN2-disorder. Despite Reticulon-2 being an endoplasmic reticulum (ER)-resident membrane shaping protein, our analysis of patient fibroblast cells did not find significant alterations in ER structure or the response to ER stress. Our findings delineate a distinct form of autosomal recessive dHMN with pyramidal features associated with Reticulon-2 deficiency. This phenotype shares similarities with SIGMAR1-related dHMN, and Silver-like syndromes, providing valuable insights into the clinical spectrum and potential therapeutic strategies for RTN2-related dHMN.

Biallelic NAA60 variants with impaired n-terminal acetylation capacity cause autosomal recessive primary familial brain calcifications.

Primary familial brain calcification (PFBC) is characterized by calcium deposition in the brain, causing progressive movement disorders, psychiatric symptoms, and cognitive decline. PFBC is a heterogeneous disorder currently linked to variants in six different genes, but most patients remain genetically undiagnosed. Here, we identify biallelic NAA60 variants in ten individuals from seven families with autosomal recessive PFBC. The NAA60 variants lead to loss-of-function with lack of protein N-terminal (Nt)-acetylation activity. We show that the phosphate importer SLC20A2 is a substrate of NAA60 in vitro. In cells, loss of NAA60 caused reduced surface levels of SLC20A2 and a reduction in extracellular phosphate uptake. This study establishes NAA60 as a causal gene for PFBC, provides a possible biochemical explanation of its disease-causing mechanisms and underscores NAA60-mediated Nt-acetylation of transmembrane proteins as a fundamental process for healthy neurobiological functioning.

Bi-allelic genetic variants in the translational GTPases GTPBP1 and GTPBP2 cause a distinct identical neurodevelopmental syndrome.

The homologous genes GTPBP1 and GTPBP2 encode GTP-binding proteins 1 and 2, which are involved in ribosomal homeostasis. Pathogenic variants in GTPBP2 were recently shown to be an ultra-rare cause of neurodegenerative or neurodevelopmental disorders (NDDs). Until now, no human phenotype has been linked to GTPBP1. Here, we describe individuals carrying bi-allelic GTPBP1 variants that display an identical phenotype with GTPBP2 and characterize the overall spectrum of GTP-binding protein (1/2)-related disorders. In this study, 20 individuals from 16 families with distinct NDDs and syndromic facial features were investigated by whole-exome (WES) or whole-genome (WGS) sequencing. To assess the functional impact of the identified genetic variants, semi-quantitative PCR, western blot, and ribosome profiling assays were performed in fibroblasts from affected individuals. We also investigated the effect of reducing expression of CG2017, an ortholog of human GTPBP1/2, in the fruit fly Drosophila melanogaster. Individuals with bi-allelic GTPBP1 or GTPBP2 variants presented with microcephaly, profound neurodevelopmental impairment, pathognomonic craniofacial features, and ectodermal defects. Abnormal vision and/or hearing, progressive spasticity, choreoathetoid movements, refractory epilepsy, and brain atrophy were part of the core phenotype of this syndrome. Cell line studies identified a loss-of-function (LoF) impact of the disease-associated variants but no significant abnormalities on ribosome profiling. Reduced expression of CG2017 isoforms was associated with locomotor impairment in Drosophila. In conclusion, bi-allelic GTPBP1 and GTPBP2 LoF variants cause an identical, distinct neurodevelopmental syndrome. Mutant CG2017 knockout flies display motor impairment, highlighting the conserved role for GTP-binding proteins in CNS development across species.

Biallelic NUDT2 variants defective in mRNA decapping cause a neurodevelopmental disease.

Dysfunctional RNA processing caused by genetic defects in RNA processing enzymes has a profound impact on the nervous system, resulting in neurodevelopmental conditions. We characterized a recessive neurological disorder in 18 children and young adults from 10 independent families typified by intellectual disability, motor developmental delay and gait disturbance. In some patients peripheral neuropathy, corpus callosum abnormalities and progressive basal ganglia deposits were present. The disorder is associated with rare variants in NUDT2, a mRNA decapping and Ap4A hydrolysing enzyme, including novel missense and in-frame deletion variants. We show that these NUDT2 variants lead to a marked loss of enzymatic activity, strongly implicating loss of NUDT2 function as the cause of the disorder. NUDT2-deficient patient fibroblasts exhibit a markedly altered transcriptome, accompanied by changes in mRNA half-life and stability. Amongst the most up-regulated mRNAs in NUDT2-deficient cells, we identified host response and interferon-responsive genes. Importantly, add-back experiments using an Ap4A hydrolase defective in mRNA decapping highlighted loss of NUDT2 decapping as the activity implicated in altered mRNA homeostasis. Our results confirm that reduction or loss of NUDT2 hydrolase activity is associated with a neurological disease, highlighting the importance of a physiologically balanced mRNA processing machinery for neuronal development and homeostasis.

Structural and non-coding variants increase the diagnostic yield of clinical whole genome sequencing for rare diseases.

BACKGROUND: Whole genome sequencing is increasingly being used for the diagnosis of patients with rare diseases. However, the diagnostic yields of many studies, particularly those conducted in a healthcare setting, are often disappointingly low, at 25-30%. This is in part because although entire genomes are sequenced, analysis is often confined to in silico gene panels or coding regions of the genome. METHODS: We undertook WGS on a cohort of 122 unrelated rare disease patients and their relatives (300 genomes) who had been pre-screened by gene panels or arrays. Patients were recruited from a broad spectrum of clinical specialties. We applied a bioinformatics pipeline that would allow comprehensive analysis of all variant types. We combined established bioinformatics tools for phenotypic and genomic analysis with our novel algorithms (SVRare, ALTSPLICE and GREEN-DB) to detect and annotate structural, splice site and non-coding variants. RESULTS: Our diagnostic yield was 43/122 cases (35%), although 47/122 cases (39%) were considered solved when considering novel candidate genes with supporting functional data into account. Structural, splice site and deep intronic variants contributed to 20/47 (43%) of our solved cases. Five genes that are novel, or were novel at the time of discovery, were identified, whilst a further three genes are putative novel disease genes with evidence of causality. We identified variants of uncertain significance in a further fourteen candidate genes. The phenotypic spectrum associated with RMND1 was expanded to include polymicrogyria. Two patients with secondary findings in FBN1 and KCNQ1 were confirmed to have previously unidentified Marfan and long QT syndromes, respectively, and were referred for further clinical interventions. Clinical diagnoses were changed in six patients and treatment adjustments made for eight individuals, which for five patients was considered life-saving. CONCLUSIONS: Genome sequencing is increasingly being considered as a first-line genetic test in routine clinical settings and can make a substantial contribution to rapidly identifying a causal aetiology for many patients, shortening their diagnostic odyssey. We have demonstrated that structural, splice site and intronic variants make a significant contribution to diagnostic yield and that comprehensive analysis of the entire genome is essential to maximise the value of clinical genome sequencing.

Lunapark deficiency leads to an autosomal recessive neurodevelopmental phenotype with a degenerative course, epilepsy and distinct brain anomalies.

LNPK encodes a conserved membrane protein that stabilizes the junctions of the tubular endoplasmic reticulum network playing crucial roles in diverse biological functions. Recently, homozygous variants in LNPK were shown to cause a neurodevelopmental disorder (OMIM#618090) in four patients displaying developmental delay, epilepsy and nonspecific brain malformations including corpus callosum hypoplasia and variable impairment of cerebellum. We sought to delineate the molecular and phenotypic spectrum of LNPK-related disorder. Exome or genome sequencing was carried out in 11 families. Thorough clinical and neuroradiological evaluation was performed for all the affected individuals, including review of previously reported patients. We identified 12 distinct homozygous loss-of-function variants in 16 individuals presenting with moderate to profound developmental delay, cognitive impairment, regression, refractory epilepsy and a recognizable neuroimaging pattern consisting of corpus callosum hypoplasia and signal alterations of the forceps minor ('ear-of-the-lynx' sign), variably associated with substantia nigra signal alterations, mild brain atrophy, short midbrain and cerebellar hypoplasia/atrophy. In summary, we define the core phenotype of LNPK-related disorder and expand the list of neurological disorders presenting with the 'ear-of-the-lynx' sign suggesting a possible common underlying mechanism related to endoplasmic reticulum-phagy dysfunction.

Pure cerebellar ataxia due to bi-allelic PRDX3 variants including recurring p.Asp202Asn.

Bi-allelic variants in peroxiredoxin 3 (PRDX3) have only recently been associated with autosomal recessive spinocerebellar ataxia characterized by early onset slowly progressive cerebellar ataxia, variably associated with hyperkinetic and hypokinetic features, accompanied by cerebellar atrophy and occasional olivary and brainstem involvement. Herein, we describe a further simplex case carrying a reported PRDX3 variant as well as two additional cases with novel variants. We report the first Brazilian patient with SCAR32, replicating the pathogenic status of a known variant. All presented cases from the Brazilian and Indian populations expand the phenotypic spectrum of the disease by displaying prominent neuroradiological findings. SCAR32, although rare, should be included in the differential diagnosis of sporadic or recessive childhood and adolescent-onset pure and complex cerebellar ataxia.

Use of genome sequencing to hunt for cryptic second-hit variants: analysis of 31 cases recruited to the 100 000 Genomes Project.

BACKGROUND: Current clinical testing methods used to uncover the genetic basis of rare disease have inherent limitations, which can lead to causative pathogenic variants being missed. Within the rare disease arm of the 100 000 Genomes Project (100kGP), families were recruited under the clinical indication 'single autosomal recessive mutation in rare disease'. These participants presented with strong clinical suspicion for a specific autosomal recessive disorder, but only one suspected pathogenic variant had been identified through standard-of-care testing. Whole genome sequencing (WGS) aimed to identify cryptic 'second-hit' variants. METHODS: To investigate the 31 families with available data that remained unsolved following formal review within the 100kGP, SVRare was used to aggregate structural variants present in <1% of 100kGP participants. Small variants were assessed using population allele frequency data and SpliceAI. Literature searches and publicly available online tools were used for further annotation of pathogenicity. RESULTS: Using these strategies, 8/31 cases were solved, increasing the overall diagnostic yield of this cohort from 10/41 (24.4%) to 18/41 (43.9%). Exemplar cases include a patient with cystic fibrosis harbouring a novel exonic LINE1 insertion in CFTR and a patient with generalised arterial calcification of infancy with complex interlinked duplications involving exons 2-6 of ENPP1. Although ambiguous by short-read WGS, the ENPP1 variant structure was resolved using optical genome mapping and RNA analysis. CONCLUSION: Systematic examination of cryptic variants across a multi-disease cohort successfully identifies additional pathogenic variants. WGS data analysis in autosomal recessive rare disease should consider complex structural and small intronic variants as potentially pathogenic second hits.

CUX1-related neurodevelopmental disorder: deep insights into phenotype-genotype spectrum and underlying pathology.

Heterozygous, pathogenic CUX1 variants are associated with global developmental delay or intellectual disability. This study delineates the clinical presentation in an extended cohort and investigates the molecular mechanism underlying the disorder in a Cux1+/- mouse model. Through international collaboration, we assembled the phenotypic and molecular information for 34 individuals (23 unpublished individuals). We analyze brain CUX1 expression and susceptibility to epilepsy in Cux1+/- mice. We describe 34 individuals, from which 30 were unrelated, with 26 different null and four missense variants. The leading symptoms were mild to moderate delayed speech and motor development and borderline to moderate intellectual disability. Additional symptoms were muscular hypotonia, seizures, joint laxity, and abnormalities of the forehead. In Cux1+/- mice, we found delayed growth, histologically normal brains, and increased susceptibility to seizures. In Cux1+/- brains, the expression of Cux1 transcripts was half of WT animals. Expression of CUX1 proteins was reduced, although in early postnatal animals significantly more than in adults. In summary, disease-causing CUX1 variants result in a non-syndromic phenotype of developmental delay and intellectual disability. In some individuals, this phenotype ameliorates with age, resulting in a clinical catch-up and normal IQ in adulthood. The post-transcriptional balance of CUX1 expression in the heterozygous brain at late developmental stages appears important for this favorable clinical course.

Genome sequencing identifies KMT2E-disrupting cryptic structural variant in a female with O'Donnell-Luria-Rodan syndrome.

We describe a patient from the 100,000 Genomes Project with a complex de novo structural variant within KMT2E leading to O'Donnell-Luria-Rodan syndrome. This case expands the mutational spectrum for this syndrome and highlights the importance of revisiting unsolved cases using better SV prioritisation tools and updated gene panels.

AMFR dysfunction causes autosomal recessive spastic paraplegia in human that is amenable to statin treatment in a preclinical model.

Hereditary spastic paraplegias (HSP) are rare, inherited neurodegenerative or neurodevelopmental disorders that mainly present with lower limb spasticity and muscle weakness due to motor neuron dysfunction. Whole genome sequencing identified bi-allelic truncating variants in AMFR, encoding a RING-H2 finger E3 ubiquitin ligase anchored at the membrane of the endoplasmic reticulum (ER), in two previously genetically unexplained HSP-affected siblings. Subsequently, international collaboration recognized additional HSP-affected individuals with similar bi-allelic truncating AMFR variants, resulting in a cohort of 20 individuals from 8 unrelated, consanguineous families. Variants segregated with a phenotype of mainly pure but also complex HSP consisting of global developmental delay, mild intellectual disability, motor dysfunction, and progressive spasticity. Patient-derived fibroblasts, neural stem cells (NSCs), and in vivo zebrafish modeling were used to investigate pathomechanisms, including initial preclinical therapy assessment. The absence of AMFR disturbs lipid homeostasis, causing lipid droplet accumulation in NSCs and patient-derived fibroblasts which is rescued upon AMFR re-expression. Electron microscopy indicates ER morphology alterations in the absence of AMFR. Similar findings are seen in amfra-/- zebrafish larvae, in addition to altered touch-evoked escape response and defects in motor neuron branching, phenocopying the HSP observed in patients. Interestingly, administration of FDA-approved statins improves touch-evoked escape response and motor neuron branching defects in amfra-/- zebrafish larvae, suggesting potential therapeutic implications. Our genetic and functional studies identify bi-allelic truncating variants in AMFR as a cause of a novel autosomal recessive HSP by altering lipid metabolism, which may potentially be therapeutically modulated using precision medicine with statins.

The prevalence and phenotypic range associated with biallelic PKDCC variants.

PKDCC encodes a component of Hedgehog signalling required for normal chondrogenesis and skeletal development. Although biallelic PKDCC variants have been implicated in rhizomelic shortening of limbs with variable dysmorphic features, this association was based on just two patients. In this study, data from the 100 000 Genomes Project was used in conjunction with exome sequencing and panel-testing results accessed via international collaboration to assemble a cohort of eight individuals from seven independent families with biallelic PKDCC variants. The allelic series included six frameshifts, a previously described splice-donor site variant and a likely pathogenic missense variant observed in two families that was supported by in silico structural modelling. Database queries suggested that the prevalence of this condition is between 1 of 127 and 1 of 721 in clinical cohorts with skeletal dysplasia of unknown aetiology. Clinical assessments, combined with data from previously published cases, indicate a predominantly upper limb involvement. Micrognathia, hypertelorism and hearing loss appear to be commonly co-occurring features. In conclusion, this study strengthens the link between biallelic inactivation of PKDCC and rhizomelic limb-shortening and will enable clinical testing laboratories to better interpret variants in this gene.

Conclusion of diagnostic odysseys due to inversions disrupting GLI3 and FBN1.

Many genetic testing methodologies are biased towards picking up structural variants (SVs) that alter copy number. Copy-neutral rearrangements such as inversions are therefore likely to suffer from underascertainment. In this study, manual review prompted by a virtual multidisciplinary team meeting and subsequent bioinformatic prioritisation of data from the 100K Genomes Project was performed across 43 genes linked to well-characterised skeletal disorders. Ten individuals from three independent families were found to harbour diagnostic inversions. In two families, inverted segments of 1.2/14.8 Mb unequivocally disrupted GLI3 and segregated with skeletal features consistent with Greig cephalopolysyndactyly syndrome. For one family, phenotypic blending was due to the opposing breakpoint lying ~45 kb from HOXA13 In the third family, long suspected to have Marfan syndrome, a 2.0 Mb inversion disrupting FBN1 was identified. These findings resolved lengthy diagnostic odysseys of 9-20 years and highlight the importance of direct interaction between clinicians and data-analysts. These exemplars of a rare mutational class inform future SV prioritisation strategies within the NHS Genomic Medicine Service and similar genome sequencing initiatives. In over 30 years since these two disease-gene associations were identified, large inversions have yet to be described and so our results extend the mutational spectra linked to these conditions.

Biallelic variants in PIGN cause Fryns syndrome, multiple congenital anomalies-hypotonia-seizures syndrome, and neurologic phenotypes: A genotype-phenotype correlation study.

PURPOSE: Biallelic PIGN variants have been described in Fryns syndrome, multiple congenital anomalies-hypotonia-seizure syndrome (MCAHS), and neurologic phenotypes. The full spectrum of clinical manifestations in relation to the genotypes is yet to be reported. METHODS: Genotype and phenotype data were collated and analyzed for 61 biallelic PIGN cases: 21 new and 40 previously published cases. Functional analysis was performed for 2 recurrent variants (c.2679C>G p.Ser893Arg and c.932T>G p.Leu311Trp). RESULTS: Biallelic-truncating variants were detected in 16 patients-10 with Fryns syndrome, 1 with MCAHS1, 2 with Fryns syndrome/MCAHS1, and 3 with neurologic phenotype. There was an increased risk of prenatal or neonatal death within this group (6 deaths were in utero or within 2 months of life; 6 pregnancies were terminated). Incidence of polyhydramnios, congenital anomalies (eg, diaphragmatic hernia), and dysmorphism was significantly increased. Biallelic missense or mixed genotype were reported in the remaining 45 cases-32 showed a neurologic phenotype and 12 had MCAHS1. No cases of diaphragmatic hernia or abdominal wall defects were seen in this group except patient 1 in which we found the missense variant p.Ser893Arg to result in functionally null alleles, suggesting the possibility of an undescribed functionally important region in the final exon. For all genotypes, there was complete penetrance for developmental delay and near-complete penetrance for seizures and hypotonia in patients surviving the neonatal period. CONCLUSION: We have expanded the described spectrum of phenotypes and natural history associated with biallelic PIGN variants. Our study shows that biallelic-truncating variants usually result in the more severe Fryns syndrome phenotype, but neurologic problems, such as developmental delay, seizures, and hypotonia, present across all genotypes. Functional analysis should be considered when the genotypes do not correlate with the predicted phenotype because there may be other functionally important regions in PIGN that are yet to be discovered.

Brain monoamine vesicular transport disease caused by homozygous SLC18A2 variants: A study in 42 affected individuals.

PURPOSE: Brain monoamine vesicular transport disease is an infantile-onset movement disorder that mimics cerebral palsy. In 2013, the homozygous SLC18A2 variant, p.Pro387Leu, was first reported as a cause of this rare disorder, and dopamine agonists were efficient for treating affected individuals from a single large family. To date, only 6 variants have been reported. In this study, we evaluated genotype-phenotype correlations in individuals with biallelic SLC18A2 variants. METHODS: A total of 42 affected individuals with homozygous SLC18A2 variant alleles were identified. We evaluated genotype-phenotype correlations and the missense variants in the affected individuals based on the structural modeling of rat VMAT2 encoded by Slc18a2, with cytoplasm- and lumen-facing conformations. A Caenorhabditis elegans model was created for functional studies. RESULTS: A total of 19 homozygous SLC18A2 variants, including 3 recurrent variants, were identified using exome sequencing. The affected individuals typically showed global developmental delay, hypotonia, dystonia, oculogyric crisis, and autonomic nervous system involvement (temperature dysregulation/sweating, hypersalivation, and gastrointestinal dysmotility). Among the 58 affected individuals described to date, 16 (28%) died before the age of 13 years. Of the 17 patients with p.Pro237His, 9 died, whereas all 14 patients with p.Pro387Leu survived. Although a dopamine agonist mildly improved the disease symptoms in 18 of 21 patients (86%), some affected individuals with p.Ile43Phe and p.Pro387Leu showed milder phenotypes and presented prolonged survival even without treatment. The C. elegans model showed behavioral abnormalities. CONCLUSION: These data expand the phenotypic and genotypic spectra of SLC18A2-related disorders.

Rare disease gene association discovery from burden analysis of the 100,000 Genomes Project data.

To discover rare disease-gene associations, we developed a gene burden analytical framework and applied it to rare, protein-coding variants from whole genome sequencing of 35,008 cases with rare diseases and their family members recruited to the 100,000 Genomes Project (100KGP). Following in silico triaging of the results, 88 novel associations were identified including 38 with existing experimental evidence. We have published the confirmation of one of these associations, hereditary ataxia with UCHL1 , and independent confirmatory evidence has recently been published for four more. We highlight a further seven compelling associations: hypertrophic cardiomyopathy with DYSF and SLC4A3 where both genes show high/specific heart expression and existing associations to skeletal dystrophies or short QT syndrome respectively; monogenic diabetes with UNC13A with a known role in the regulation of ß cells and a mouse model with impaired glucose tolerance; epilepsy with KCNQ1 where a mouse model shows seizures and the existing long QT syndrome association may be linked; early onset Parkinson's disease with RYR1 with existing links to tremor pathophysiology and a mouse model with neurological phenotypes; anterior segment ocular abnormalities associated with POMK showing expression in corneal cells and with a zebrafish model with developmental ocular abnormalities; and cystic kidney disease with COL4A3 showing high renal expression and prior evidence for a digenic or modifying role in renal disease. Confirmation of all 88 associations would lead to potential diagnoses in 456 molecularly undiagnosed cases within the 100KGP, as well as other rare disease patients worldwide, highlighting the clinical impact of a large-scale statistical approach to rare disease gene discovery.

Venus: Elucidating the Impact of Amino Acid Variants on Protein Function Beyond Structure Destabilisation.

Exploring the functional effect of a non-synonymous coding variant at the protein level requires multiple pieces of information to be interpreted appropriately. This is particularly important when embarking on the study of a potentially pathogenic variant linked to a rare or monogenic disease. Whereas accurate protein stability predictions alone are generally informative, other effects, such as disruption of post-translational modifications or weakened ligand binding, may also contribute to the disease phenotype. Furthermore, consideration of nearby variants that are found in the healthy population may strengthen or refute a given mechanistic hypothesis. Whilst there are several bioinformatics tools available that score a genetic variant in terms of deleteriousness, there is no single tool that assembles multiple effects of a variant on the encoded protein, beyond structural stability, and presents them on the structure for inspection. Venus is a web application which, given a protein substitution, rapidly estimates the predicted effect on protein stability of the variant, flags if the variant affects a post-translational modification site, a predicted linear motif or known annotation, and determines the effect on protein stability of variants which affect nearby residues and have been identified in healthy populations. Venus is built upon Michelanglo and the results can be exported to it, allowing them to be annotated and shared with other researchers. Venus is freely accessible at https://venus.cmd.ox.ac.uk and its source code is openly available at https://github.com/CMD-Oxford/Michelanglo-and-Venus.

The Bartter-Gitelman Spectrum: 50-Year Follow-up With Revision of Diagnosis After Whole-Genome Sequencing.

Bartter syndrome (BS) and Gitelman syndrome (GS) are renal tubular disorders affecting sodium, potassium, and chloride reabsorption. Clinical features include muscle cramps and weakness, in association with hypokalemia, hypochloremic metabolic alkalosis, and hyperreninemic hyperaldosteronism. Hypomagnesemia and hypocalciuria are typical of GS, while juxtaglomerular hyperplasia is characteristic of BS. GS is due to SLC12A3 variants, whereas BS is due to variants in SLC12A1, KCNJ1, CLCNKA, CLCNKB, BSND, MAGED2, or CASR. We had the opportunity to follow up one of the first reported cases of a salt-wasting tubulopathy, who based on clinical features was diagnosed with GS. The patient had presented at age 10 years with tetany precipitated by vomiting or diarrhea. She had hypokalemia, a hypochloremic metabolic alkalosis, hyponatremia, mild hypercalcemia, and normomagnesemia, and subsequently developed hypocalciuria and hypomagnesemia. A renal biopsy showed no evidence for juxtaglomerular hyperplasia. She developed chronic kidney failure at age 55 years, and ocular sclerochoroidal calcification, associated with BS and GS, at older than 65 years. Our aim was therefore to establish the genetic diagnosis in this patient using whole-genome sequencing (WGS). Leukocyte DNA was used for WGS analysis, and this revealed a homozygous c.226C > T (p.Arg76Ter) nonsense CLCNKB mutation, thereby establishing a diagnosis of BS type-3. WGS also identified 2 greater than 5-Mb regions of homozygosity that suggested likely mutational heterozygosity in her parents, who originated from a Greek island with fewer than 1500 inhabitants and may therefore have shared a common ancestor. Our results demonstrate the utility of WGS in establishing the correct diagnosis in renal tubular disorders with overlapping phenotypes.