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.

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.

RTN2 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 the scarcity of supporting evidence. In this 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 and hyperreflexia, with 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. Characterization of Caenorhabditis elegans RTN2 homologous loss-of-function variants demonstrated morphological and behavioural differences compared with 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 RTN2 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 RTN2 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 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.

Bi-allelic loss-of-function variants in BCAS3 cause a syndromic neurodevelopmental disorder.

BCAS3 microtubule-associated cell migration factor (BCAS3) is a large, highly conserved cytoskeletal protein previously proposed to be critical in angiogenesis and implicated in human embryogenesis and tumorigenesis. Here, we established BCAS3 loss-of-function variants as causative for a neurodevelopmental disorder. We report 15 individuals from eight unrelated families with germline bi-allelic loss-of-function variants in BCAS3. All probands share a global developmental delay accompanied by pyramidal tract involvement, microcephaly, short stature, strabismus, dysmorphic facial features, and seizures. The human phenotype is less severe compared with the Bcas3 knockout mouse model and cannot be explained by angiogenic defects alone. Consistent with being loss-of-function alleles, we observed absence of BCAS3 in probands' primary fibroblasts. By comparing the transcriptomic and proteomic data based on probands' fibroblasts with those of the knockout mouse model, we identified similar dysregulated pathways resulting from over-representation analysis, while the dysregulation of some proposed key interactors could not be confirmed. Together with the results from a tissue-specific Drosophila loss-of-function model, we demonstrate a vital role for BCAS3 in neural tissue development.

A human importin-ß-related disorder: Syndromic thoracic aortic aneurysm caused by bi-allelic loss-of-function variants in IPO8.

Importin 8, encoded by IPO8, is a ubiquitously expressed member of the importin-ß protein family that translocates cargo molecules such as proteins, RNAs, and ribonucleoprotein complexes into the nucleus in a RanGTP-dependent manner. Current knowledge of the cargoes of importin 8 is limited, but TGF-ß signaling components such as SMAD1-4 have been suggested to be among them. Here, we report that bi-allelic loss-of-function variants in IPO8 cause a syndromic form of thoracic aortic aneurysm (TAA) with clinical overlap with Loeys-Dietz and Shprintzen-Goldberg syndromes. Seven individuals from six unrelated families showed a consistent phenotype with early-onset TAA, motor developmental delay, connective tissue findings, and craniofacial dysmorphic features. A C57BL/6N Ipo8 knockout mouse model recapitulates TAA development from 8-12 weeks onward in both sexes but most prominently shows ascending aorta dilatation with a propensity for dissection in males. Compliance assays suggest augmented passive stiffness of the ascending aorta in male Ipo8-/- mice throughout life. Immunohistological investigation of mutant aortic walls reveals elastic fiber disorganization and fragmentation along with a signature of increased TGF-ß signaling, as evidenced by nuclear pSmad2 accumulation. RT-qPCR assays of the aortic wall in male Ipo8-/- mice demonstrate decreased Smad6/7 and increased Mmp2 and Ccn2 (Ctgf) expression, reinforcing a role for dysregulation of the TGF-ß signaling pathway in TAA development. Because importin 8 is the most downstream TGF-ß-related effector implicated in TAA pathogenesis so far, it offers opportunities for future mechanistic studies and represents a candidate drug target for TAA.

SCUBE3 loss-of-function causes a recognizable recessive developmental disorder due to defective bone morphogenetic protein signaling.

Signal peptide-CUB-EGF domain-containing protein 3 (SCUBE3) is a member of a small family of multifunctional cell surface-anchored glycoproteins functioning as co-receptors for a variety of growth factors. Here we report that bi-allelic inactivating variants in SCUBE3 have pleiotropic consequences on development and cause a previously unrecognized syndromic disorder. Eighteen affected individuals from nine unrelated families showed a consistent phenotype characterized by reduced growth, skeletal features, distinctive craniofacial appearance, and dental anomalies. In vitro functional validation studies demonstrated a variable impact of disease-causing variants on transcript processing, protein secretion and function, and their dysregulating effect on bone morphogenetic protein (BMP) signaling. We show that SCUBE3 acts as a BMP2/BMP4 co-receptor, recruits the BMP receptor complexes into raft microdomains, and positively modulates signaling possibly by augmenting the specific interactions between BMPs and BMP type I receptors. Scube3-/- mice showed craniofacial and dental defects, reduced body size, and defective endochondral bone growth due to impaired BMP-mediated chondrogenesis and osteogenesis, recapitulating the human disorder. Our findings identify a human disease caused by defective function of a member of the SCUBE family, and link SCUBE3 to processes controlling growth, morphogenesis, and bone and teeth development through modulation of BMP signaling.

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.

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.

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.

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.

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.

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.

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.

Recurrent de novo missense variants in GNB2 can cause syndromic intellectual disability.

PURPOSE: Binding proteins (G-proteins) mediate signalling pathways involved in diverse cellular functions and comprise Gα and Gßγ units. Human diseases have been reported for all five Gß proteins. A de novo missense variant in GNB2 was recently reported in one individual with developmental delay/intellectual disability (DD/ID) and dysmorphism. We aim to confirm GNB2 as a neurodevelopmental disease gene, and elucidate the GNB2-associated neurodevelopmental phenotype in a patient cohort. METHODS: We discovered a GNB2 variant in the index case via exome sequencing and sought individuals with GNB2 variants via international data-sharing initiatives. In silico modelling of the variants was assessed, along with multiple lines of evidence in keeping with American College of Medical Genetics and Genomics guidelines for interpretation of sequence variants. RESULTS: We identified 12 unrelated individuals with five de novo missense variants in GNB2, four of which are recurrent: p.(Ala73Thr), p.(Gly77Arg), p.(Lys89Glu) and p.(Lys89Thr). All individuals have DD/ID with variable dysmorphism and extraneurologic features. The variants are located at the universally conserved shared interface with the Gα subunit, which modelling suggests weaken this interaction. CONCLUSION: Missense variants in GNB2 cause a congenital neurodevelopmental disorder with variable syndromic features, broadening the spectrum of multisystem phenotypes associated with variants in genes encoding G-proteins.

De Novo Heterozygous POLR2A Variants Cause a Neurodevelopmental Syndrome with Profound Infantile-Onset Hypotonia.

The RNA polymerase II complex (pol II) is responsible for transcription of all ∼21,000 human protein-encoding genes. Here, we describe sixteen individuals harboring de novo heterozygous variants in POLR2A, encoding RPB1, the largest subunit of pol II. An iterative approach combining structural evaluation and mass spectrometry analyses, the use of S. cerevisiae as a model system, and the assessment of cell viability in HeLa cells allowed us to classify eleven variants as probably disease-causing and four variants as possibly disease-causing. The significance of one variant remains unresolved. By quantification of phenotypic severity, we could distinguish mild and severe phenotypic consequences of the disease-causing variants. Missense variants expected to exert only mild structural effects led to a malfunctioning pol II enzyme, thereby inducing a dominant-negative effect on gene transcription. Intriguingly, individuals carrying these variants presented with a severe phenotype dominated by profound infantile-onset hypotonia and developmental delay. Conversely, individuals carrying variants expected to result in complete loss of function, thus reduced levels of functional pol II from the normal allele, exhibited the mildest phenotypes. We conclude that subtle variants that are central in functionally important domains of POLR2A cause a neurodevelopmental syndrome characterized by profound infantile-onset hypotonia and developmental delay through a dominant-negative effect on pol-II-mediated transcription of DNA.

Biallelic TMEM260 variants cause truncus arteriosus, with or without renal defects.

Only two families have been reported with biallelic TMEM260 variants segregating with structural heart defects and renal anomalies syndrome (SHDRA). With a combination of genome, exome sequencing and RNA studies, we identified eight individuals from five families with biallelic TMEM260 variants. Variants included one multi-exon deletion, four nonsense/frameshifts, two splicing changes and one missense change. Together with the published cases, analysis of clinical data revealed ventricular septal defects (12/12), mostly secondary to truncus arteriosus (10/12), elevated creatinine levels (6/12), horse-shoe kidneys (1/12) and renal cysts (1/12) in patients. Three pregnancies were terminated on detection of severe congenital anomalies. Six patients died between the ages of 6 weeks and 5 years. Using a range of stringencies, carrier frequency for SHDRA was estimated at 0.0007-0.007 across ancestries. In conclusion, this study confirms the genetic basis of SHDRA, expands its known mutational spectrum and clarifies its clinical features. We demonstrate that SHDRA is a severe condition associated with substantial mortality in early childhood and characterised by congenital cardiac malformations with a variable renal phenotype.

Biallelic variants in HPDL cause pure and complicated hereditary spastic paraplegia.

Human 4-hydroxyphenylpyruvate dioxygenase-like (HPDL) is a putative iron-containing non-heme oxygenase of unknown specificity and biological significance. We report 25 families containing 34 individuals with neurological disease associated with biallelic HPDL variants. Phenotypes ranged from juvenile-onset pure hereditary spastic paraplegia to infantile-onset spasticity and global developmental delays, sometimes complicated by episodes of neurological and respiratory decompensation. Variants included bona fide pathogenic truncating changes, although most were missense substitutions. Functionality of variants could not be determined directly as the enzymatic specificity of HPDL is unknown; however, when HPDL missense substitutions were introduced into 4-hydroxyphenylpyruvate dioxygenase (HPPD, an HPDL orthologue), they impaired the ability of HPPD to convert 4-hydroxyphenylpyruvate into homogentisate. Moreover, three additional sets of experiments provided evidence for a role of HPDL in the nervous system and further supported its link to neurological disease: (i) HPDL was expressed in the nervous system and expression increased during neural differentiation; (ii) knockdown of zebrafish hpdl led to abnormal motor behaviour, replicating aspects of the human disease; and (iii) HPDL localized to mitochondria, consistent with mitochondrial disease that is often associated with neurological manifestations. Our findings suggest that biallelic HPDL variants cause a syndrome varying from juvenile-onset pure hereditary spastic paraplegia to infantile-onset spastic tetraplegia associated with global developmental delays.

An ancestral 10-bp repeat expansion in VWA1 causes recessive hereditary motor neuropathy.

The extracellular matrix comprises a network of macromolecules such as collagens, proteoglycans and glycoproteins. VWA1 (von Willebrand factor A domain containing 1) encodes a component of the extracellular matrix that interacts with perlecan/collagen VI, appears to be involved in stabilizing extracellular matrix structures, and demonstrates high expression levels in tibial nerve. Vwa1-deficient mice manifest with abnormal peripheral nerve structure/function; however, VWA1 variants have not previously been associated with human disease. By interrogating the genome sequences of 74 180 individuals from the 100K Genomes Project in combination with international gene-matching efforts and targeted sequencing, we identified 17 individuals from 15 families with an autosomal-recessive, non-length dependent, hereditary motor neuropathy and rare biallelic variants in VWA1. A single disease-associated allele p.(G25Rfs*74), a 10-bp repeat expansion, was observed in 14/15 families and was homozygous in 10/15. Given an allele frequency in European populations approaching 1/1000, the seven unrelated homozygote individuals ascertained from the 100K Genomes Project represents a substantial enrichment above expected. Haplotype analysis identified a shared 220 kb region suggesting that this founder mutation arose >7000 years ago. A wide age-range of patients (6-83 years) helped delineate the clinical phenotype over time. The commonest disease presentation in the cohort was an early-onset (mean 2.0 ± 1.4 years) non-length-dependent axonal hereditary motor neuropathy, confirmed on electrophysiology, which will have to be differentiated from other predominantly or pure motor neuropathies and neuronopathies. Because of slow disease progression, ambulation was largely preserved. Neurophysiology, muscle histopathology, and muscle MRI findings typically revealed clear neurogenic changes with single isolated cases displaying additional myopathic process. We speculate that a few findings of myopathic changes might be secondary to chronic denervation rather than indicating an additional myopathic disease process. Duplex reverse transcription polymerase chain reaction and immunoblotting using patient fibroblasts revealed that the founder allele results in partial nonsense mediated decay and an absence of detectable protein. CRISPR and morpholino vwa1 modelling in zebrafish demonstrated reductions in motor neuron axonal growth, synaptic formation in the skeletal muscles and locomotive behaviour. In summary, we estimate that biallelic variants in VWA1 may be responsible for up to 1% of unexplained hereditary motor neuropathy cases in Europeans. The detailed clinical characterization provided here will facilitate targeted testing on suitable patient cohorts. This novel disease gene may have previously evaded detection because of high GC content, consequential low coverage and computational difficulties associated with robustly detecting repeat-expansions. Reviewing previously unsolved exomes using lower QC filters may generate further diagnoses.