Search

Abstract We report the results of a comprehensive copy number variant (CNV) reanalysis of 9171 exome sequencing datasets from 5757 families affected by a rare disease (RD). The data reanalysed was extremely heterogeneous, having been generated using 28 different enrichment kits by 42 different research groups across Europe partnering in the Solve-RD project. Each research group had previously undertaken their own analysis of the data but failed to identify disease-causing variants. We applied three CNV calling algorithms to maximise sensitivity, and rare CNVs overlapping genes of interest, provided by four partner European Reference Networks, were taken forward for interpretation by clinical experts. This reanalysis has resulted in a molecular diagnosis being provided to 51 families in this sample, with ClinCNV performing the best of the three algorithms. We also identified partially explanatory pathogenic CNVs in a further 34 individuals. This work illustrates the value of reanalysing ES cold cases for CNVs.

Abstract Exploring the molecular basis of disease severity in rare disease scenarios is a challenging task provided the limitations on data availability. Causative genes have been described for Congenital Myasthenic Syndromes (CMS), a group of diverse minority neuromuscular junction (NMJ) disorders; yet a molecular explanation for the phenotypic severity differences remains unclear. Here, we present a workflow to explore the functional relationships between CMS causal genes and altered genes from each patient, based on multilayer network community detection analysis of complementary biomedical information provided by relevant data sources, namely protein-protein interactions, pathways and metabolomics. Our results show that CMS severity can be ascribed to the personalized impairment of extracellular matrix components and postsynaptic modulators of acetylcholine receptor (AChR) clustering. This work showcases how coupling multilayer network analysis with personalized -omics information provides molecular explanations to the varying severity of rare diseases; paving the way for sorting out similar cases in other rare diseases.

Abstract Objective The Global FKRP Registry is a database for individuals with conditions caused by mutations in the Fukutin‐Related Protein (FKRP) gene: limb girdle muscular dystrophy R9 (LGMDR9, formerly LGMD2I) and congenital muscular dystrophies MDC1C, Muscle–Eye–Brain Disease and Walker–Warburg Syndrome. The registry seeks to further understand the natural history and prevalence of FKRP‐related conditions; aid the rapid identification of eligible patients for clinical studies; and provide a source of information to clinical and academic communities. Methods Registration is patient‐initiated through a secure online portal. Data, reported by both patients and their clinicians, include: age of onset, presenting symptoms, family history, motor function and muscle strength, respiratory and cardiac function, medication, quality of life and pain. Results Of 663 registered participants, 305 were genetically confirmed LGMDR9 patients from 23 countries. A majority of LGMDR9 patients carried the common mutation c.826C > A on one or both alleles; 67.9% were homozygous and 28.5% were compound heterozygous for this mutation. The mean ages of symptom onset and disease diagnosis were higher in individuals homozygous for c.826C > A compared with individuals heterozygous for c.826C > A. This divergence was replicated in ages of loss of running ability, wheelchair‐dependence and ventilation assistance; consistent with the milder phenotype associated with individuals homozygous for c.826C > A. In LGMDR9 patients, 75.1% were currently ambulant and 24.6%, nonambulant (unreported in 0.3%). Cardiac impairment was reported in 23.2% (30/129). Interpretation The Global FKRP Registry enables the collection of patient natural history data, which informs academics, healthcare professionals and industry. It represents a trial‐ready cohort of individuals and is centrally placed to facilitate recruitment to clinical studies.

Abstract Background Limb girdle muscular dystrophies are a group of rare and genetically heterogeneous diseases that share proximal weakness as a common feature; however they are often lacking very specific phenotypic features to allow an accurate differential diagnosis based on the clinical signs only, limiting the diagnostic rate using phenotype driven genetic testing. Next generation sequencing provides an opportunity to obtain molecular diagnoses for undiagnosed patients, as well as identifying novel genetic causes of muscle diseases. We performed whole exome sequencing (WES) on 104 affected individuals from 75 families in who standard gene by gene testing had not yielded a diagnosis. For comparison we also evaluated the diagnostic rate using sequential gene by gene testing for 91 affected individuals from 84 families over a 2 year period. Results Patients selected for WES had undergone more extensive prior testing than those undergoing standard genetic testing and on average had had 8 genes screened already. In this extensively investigated cohort WES identified the genetic diagnosis in 28 families (28/75, 37%), including the identification of the novel gene ZAK and two unpublished genes. WES of a single affected individual with sporadic disease yielded a diagnosis in 13/38 (34%) of cases. In comparison, conventional gene by gene testing provided a genetic diagnosis in 28/84 (33%) families. Titinopathies and collagen VI related dystrophy were the most frequent diagnoses made by WES. Reasons why mutations in known genes were not identified previously included atypical phenotypes, reassignment of pathogenicity of variants, and in one individual mosaicism for a COL6A1 mutation which was undetected by prior direct sequencing. Conclusion WES was able to overcome many limitations of standard testing and achieved a higher rate of diagnosis than standard testing even in this cohort of extensively investigated patients. Earlier application of WES is therefore likely to yield an even higher diagnostic rate. We obtained a high diagnosis rate in simplex cases and therefore such individuals should be included in exome or genome sequencing projects. Disease due to somatic mosaicism may be increasingly recognised due to the increased sensitivity of next generation sequencing techniques to detect low level mosaicism.

Ectodermal dysplasias (EDs) are a large and heterogeneous group of hereditary disorders characterized by abnormalities in structures of ectodermal origin. Incontinentia pigmenti (IP) is an ED characterized by skin lesions evolving over time, as well as dental, nail, and ocular abnormalities. Due to X-linked dominant inheritance IP symptoms can only be seen in female individuals while affected males die during development in utero. We observed a family of horses, in which several mares developed signs of a skin disorder reminiscent of human IP. Cutaneous manifestations in affected horses included the development of pruritic, exudative lesions soon after birth. These developed into wart-like lesions and areas of alopecia with occasional wooly hair re-growth. Affected horses also had streaks of darker and lighter coat coloration from birth. The observation that only females were affected together with a high number of spontaneous abortions suggested an X-linked dominant mechanism of transmission. Using next generation sequencing we sequenced the whole genome of one affected mare. We analyzed the sequence data for non-synonymous variants in candidate genes and found a heterozygous nonsense variant in the X-chromosomal IKBKG gene (c.184C>T; p.Arg62*). Mutations in IKBKG were previously reported to cause IP in humans and the homologous p.Arg62* variant has already been observed in a human IP patient. The comparative data thus strongly suggest that this is also the causative variant for the observed IP in horses. To our knowledge this is the first large animal model for IP.

Several previous studies have investigated the role of common promoter variants in the vascular endothelial growth factor (VEGF) gene in causing congenital cardiovascular malformation (CVM). However, results have been discrepant between studies and no study to date has comprehensively characterised variation throughout the gene. We genotyped 771 CVM cases, of whom 595 had the outflow tract malformation Tetralogy of Fallot (TOF), and carried out TDT and case-control analyses using haplotype-tagging SNPs in VEGF. We carried out a meta-analysis of previous case-control or family-based studies that had typed VEGF promoter SNPs, which included an additional 570 CVM cases. To identify rare variants potentially causative of CVM, we carried out mutation screening in all VEGF exons and splice sites in 93 TOF cases. There was no significant effect of any VEGF haplotype-tagging SNP on the risk of CVM in our analyses of 771 probands. When the results of this and all previous studies were combined, there was no significant effect of the VEGF promoter SNPs rs699947 (OR 1.05 [95% CI 0.95-1.17]); rs1570360 (OR 1.17 [95% CI 0.99-1.26]); and rs2010963 (OR 1.04 [95% CI 0.93-1.16]) on the risk of CVM in 1341 cases. Mutation screening of 93 TOF cases revealed no VEGF coding sequence variants and no changes at splice consensus sequences. Genetic variation in VEGF appears to play a small role, if any, in outflow tract CVM susceptibility.

Duchenne and Becker muscular dystrophies, caused by pathogenic variants in DMD, are the most common inherited neuromuscular conditions in childhood. These diseases follow an X-linked recessive inheritance pattern, and mainly males are affected. The most prevalent pathogenic variants in the DMD gene are copy number variants (CNVs), and most patients achieve their genetic diagnosis through Multiplex Ligation-dependent Probe Amplification (MLPA) or exome sequencing. Here, we investigated a female patient presenting with muscular dystrophy who remained genetically undiagnosed after MLPA and exome sequencing. RNA sequencing (RNAseq) from the patient’s muscle biopsy identified an 85% reduction in DMD expression compared to 116 muscle samples included in the cohort. A de novo balanced translocation between chromosome 17 and the X chromosome (t(X;17)(p21.1;q23.2)) disrupting the DMD and BCAS3 genes was identified through trio whole genome sequencing (WGS). The combined analysis of RNAseq and WGS played a crucial role in the detection and characterisation of the disease-causing variant in this patient, who had been undiagnosed for over two decades. This case illustrates the diagnostic odyssey of female DMD patients with complex structural variants that are not detected by current panel or exome sequencing analysis.

Summary: TSPEAR variants cause autosomal recessive ectodermal dysplasia (ARED) 14. The function of TSPEAR is unknown. The clinical features, the mutation spectrum, and the underlying mechanisms of ARED14 are poorly understood. Combining data from new and previously published individuals established that ARED14 is primarily characterized by dental anomalies such as conical tooth cusps and hypodontia, like those seen in individuals with WNT10A-related odontoonychodermal dysplasia. AlphaFold-predicted structure-based analysis showed that most of the pathogenic TSPEAR missense variants likely destabilize the β-propeller of the protein. Analysis of 100000 Genomes Project (100KGP) data revealed multiple founder TSPEAR variants across different populations. Mutational and recombination clock analyses demonstrated that non-Finnish European founder variants likely originated around the end of the last ice age, a period of major climatic transition. Analysis of gnomAD data showed that the non-Finnish European population TSPEAR gene-carrier rate is ∼1/140, making it one of the commonest AREDs. Phylogenetic and AlphaFold structural analyses showed that TSPEAR is an ortholog of drosophila Closca, an extracellular matrix-dependent signaling regulator. We, therefore, hypothesized that TSPEAR could have a role in enamel knot, a structure that coordinates patterning of developing tooth cusps. Analysis of mouse single-cell RNA sequencing (scRNA-seq) data revealed highly restricted expression of Tspear in clusters representing enamel knots. A tspeara−/−;tspearb−/− double-knockout zebrafish model recapitulated the clinical features of ARED14 and fin regeneration abnormalities of wnt10a knockout fish, thus suggesting interaction between tspear and wnt10a. In summary, we provide insights into the role of TSPEAR in ectodermal development and the evolutionary history, epidemiology, mechanisms, and consequences of its loss of function variants.

Missense variants in RNA-binding proteins underlie many diseases. Here the authors report an oculopharyngeal muscular dystrophy caused by heterozygous frameshift mutations in HNRNPA2B1 that alter its nucleocytoplasmic transport dynamics and result in cytoplasmic accumulation of hnRNPA2 protein.

Abstract Background Goltz syndrome (GS) is a X-linked disorder defined by defects of mesodermal- and ectodermal-derived structures and caused by PORCN mutations. Features include striated skin-pigmentation, ocular and skeletal malformations and supernumerary or hypoplastic nipples. Generally, GS is associated with in utero lethality in males and most of the reported male patients show mosaicism (only three non-mosaic surviving males have been described so far). Also, precise descriptions of neurological deficits in GS are rare and less severe phenotypes might not only be caused by mosaicism but also by less pathogenic mutations suggesting the need of a molecular genetics and functional work-up of these rare variants. Results We report two cases: one girl suffering from typical skin and skeletal abnormalities, developmental delay, microcephaly, thin corpus callosum, periventricular gliosis and drug-resistant epilepsy caused by a PORCN nonsense-mutation (c.283C > T, p.Arg95Ter). Presence of these combined neurological features indicates that CNS-vulnerability might be a guiding symptom in the diagnosis of GS patients. The other patient is a boy with a supernumerary nipple and skeletal anomalies but also, developmental delay, microcephaly, cerebral atrophy with delayed myelination and drug-resistant epilepsy as predominant features. Skin abnormalities were not observed. Genotyping revealed a novel PORCN missense-mutation (c.847G > C, p.Asp283His) absent in the Genome Aggregation Database (gnomAD) but also identified in his asymptomatic mother. Given that non-random X-chromosome inactivation was excluded in the mother, fibroblasts of the index had been analyzed for PORCN protein-abundance and -distribution, vulnerability against additional ER-stress burden as well as for protein secretion revealing changes. Conclusions Our combined findings may suggest incomplete penetrance for the p.Asp283His variant and provide novel insights into the molecular etiology of GS by adding impaired ER-function and altered protein secretion to the list of pathophysiological processes resulting in the clinical manifestation of GS.

Summary: The Solve-RD project objectives include solving undiagnosed rare diseases (RD) through collaborative research on shared genome-phenome datasets. The RD-Connect Genome-Phenome Analysis Platform (GPAP), for data collation and analysis, and the European Genome-Phenome Archive (EGA), for file storage, are two key components of the Solve-RD infrastructure. Clinical researchers can identify candidate genetic variants within the RD-Connect GPAP and, thanks to the developments presented here as part of joint ELIXIR activities, are able to remotely visualize the corresponding alignments stored at the EGA. The Global Alliance for Genomics and Health (GA4GH) htsget streaming application programming interface (API) is used to retrieve alignment slices, which are rendered by an integrated genome viewer (IGV) instance embedded in the GPAP. As a result, it is no longer necessary for over 11,000 datasets to download large alignment files to visualize them locally. This work highlights the advantages, from both the user and infrastructure perspectives, of implementing interoperability standards for establishing federated genomics data networks.

ABSTRACT Background: Congenital myasthenic syndromes (CMS) have some phenotypic overlap with seronegative myasthenia gravis (SNMG). Objective: The aim of this single center study was to assess the minimum occurrence of CMS misdiagnosed as double SNMG in a Brazilian cohort. Methods: The genetic analysis of the most common mutations in CHRNE, RAPSN, and DOK7 genes was used as the main screening tool. Results: We performed genetic analysis in 22 patients with a previous diagnosis of ‘double’ SNMG. In this study, one CMS patient was confirmed due to the presence of compound heterozygous variants in the CHRNE gene (c.130insG/p.Cys210Phe). Conclusions: This study confirmed that CMS due to CHNRE mutations can be mistaken for SNMG. In addition, our study estimated the prevalence of misdiagnosed CMS to be 4.5% in ‘double’ SNMG patients of our center. Based on our findings, genetic screening could be helpful in the diagnostic workup of patients with ‘double’ SNMG in whom differential diagnosis is recommended.

Summary: Microtubules help building the cytoskeleton of neurons and other cells. Several components of the gamma-tubulin (γ-tubulin) complex have been previously reported in human neurodevelopmental diseases. We describe two siblings from a consanguineous Turkish family with dysmorphic features, developmental delay, brain malformation, and epilepsy carrying a homozygous mutation (p.Glu311Lys) in TUBGCP2 encoding the γ-tubulin complex 2 (GCP2) protein. This variant is predicted to disrupt the electrostatic interaction of GCP2 with GCP3. In primary fibroblasts carrying the variant, we observed a faint delocalization of γ-tubulin during the cell cycle but normal GCP2 protein levels. Through mass spectrometry, we observed dysregulation of multiple proteins involved in the assembly and organization of the cytoskeleton and the extracellular matrix, controlling cellular adhesion and of proteins crucial for neuronal homeostasis including axon guidance. In summary, our functional and proteomic studies link TUBGCP2 and the γ-tubulin complex to the development of the central nervous system in humans.

Abstract Background Dystroglycanopathies are a clinically and genetically heterogeneous group of disorders that are typically characterised by limb-girdle muscle weakness. Mutations in 18 different genes have been associated with dystroglycanopathies, the encoded proteins of which typically modulate the binding of α-dystroglycan to extracellular matrix ligands by altering its glycosylation. This results in a disruption of the structural integrity of the myocyte, ultimately leading to muscle degeneration. Methods Deep phenotypic information was gathered using the PhenoTips online software for 1001 patients with unexplained limb-girdle muscle weakness from 43 different centres across 21 European and Middle Eastern countries. Whole-exome sequencing with at least 250 ng DNA was completed using an Illumina exome capture and a 38 Mb baited target. Genes known to be associated with dystroglycanopathies were analysed for disease-causing variants. Results Suspected pathogenic variants were detected in DPM3, ISPD, POMT1 and FKTN in one patient each, in POMK in two patients, in GMPPB in three patients, in FKRP in eight patients and in POMT2 in ten patients. This indicated a frequency of 2.7% for the disease group within the cohort of 1001 patients with unexplained limb-girdle muscle weakness. The phenotypes of the 27 patients were highly variable, yet with a fundamental presentation of proximal muscle weakness and elevated serum creatine kinase. Conclusions Overall, we have identified 27 patients with suspected pathogenic variants in dystroglycanopathy-associated genes. We present evidence for the genetic and phenotypic diversity of the dystroglycanopathies as a disease group, while also highlighting the advantage of incorporating next-generation sequencing into the diagnostic pathway of rare diseases.

Abstract Background Late-onset Pompe disease is a rare genetic neuromuscular disorder caused by a primary deficiency of α-glucosidase and the associated accumulation of glycogen in lysosomal vacuoles. The deficiency of α-glucosidase can often be detected using an inexpensive and readily accessible dried blood spot test when Pompe disease is suspected. Like several neuromuscular disorders, Pompe disease typically presents with progressive weakness of limb-girdle muscles and respiratory insufficiency. Due to the phenotypic heterogeneity of these disorders, however, it is often difficult for clinicians to reach a diagnosis for patients with Pompe disease. Six hundred and six patients from a European population were recruited onto our study. Inclusion criteria stipulated that index cases must present with limb-girdle weakness or elevated serum creatine kinase activity. Whole exome sequencing with at least 250 ng DNA was completed using an Illumina exome capture and a 38 Mb baited target. A panel of 169 candidate genes for limb-girdle weakness was analysed for disease-causing variants. Results A total of 35 variants within GAA were detected. Ten distinct variants in eight unrelated index cases (and four siblings not sequenced in our study) were considered disease-causing, with the patients presenting with heterogeneous phenotypes. The eight unrelated individuals were compound heterozygotes for two variants. Six patients carried the intronic splice site c.-13 T > G transversion and two of the six patients also carried the exonic p.Glu176ArgfsTer45 frameshift. Four of the ten variants were novel in their association with Pompe disease. Conclusions Here, we highlight the advantage of using whole exome sequencing as a tool for detecting, diagnosing and treating patients with rare, clinically variable genetic disorders.

Background and objective: Recessive mutations in the SIL1 gene cause Marinesco-Sjögren syndrome (MSS), a rare neuropediatric disorder. MSS-patients typically present with congenital cataracts, intellectual disability, cerebellar ataxia and progressive vacuolar myopathy. However, atypical clinical presentations associated with SIL1 mutations have been described over the last years; compound heterozygosity of SIL1 missense mutations even resulted in a phenotype not fulfilling the clinical diagnostic criteria of MSS. Thus, a read-out system to evaluate reliably the pathogenicity of amino acid changes in SIL1 is needed. Here, we aim to provide suitable cellular biomarkers enabling the robust evaluation of pathogenicity of SIL1 mutations.Methods: Five SIL1 variants including one polymorphism (p.K132Q), three known pathogenic mutations (p.V231_I232del, p.G312R, and p.L457P) and one ambiguous missense variant (p.R92W) were studied along with the wild-type proteins in Hek293 in vitro models by cell biological assays, immunoprecipitation, immunoblotting, and immunofluorescence as well as electron microscopy. Moreover, the SIL1-interactomes were interrogated by tandem-affinity-purification and subsequent mass spectrometry.Results: Our combined studies confirmed the pathogenicity of p.V231_I232del, p.G312R, and p.L457P by showing instability of the proteins as well as tendency to form aggregates. This observation is in line with altered structure of the ER-Golgi system and vacuole formation upon expression of these pathogenic SIL1-mutants as well as the presence of oxidative or ER-stress. Reduced cellular fitness along with abnormal mitochondrial architecture could also be observed. Notably, both the polymorphic p.K132Q and the ambiguous p.R92W variants did not elicit such alterations. Study of the SIL1-interactome identified POC1A as a novel binding partner of wild-type SIL1; the interaction is disrupted upon the presence of pathogenic mutants but not influenced by the presence of benign variants. Disrupted SIL1-POC1A interaction is associated with centrosome disintegration.Conclusions: We developed a combination of cellular outcome measures to evaluate the pathogenicity of SIL1 variants in suitable in vitro models and demonstrated that the p. R92W missense variant is a polymorphism rather than a pathogenic mutation leading to MSS.

OBJECTIVE:Rare variants in certain transcription factors involved in cardiac development cause Mendelian forms of congenital heart disease. The purpose of this study was to systematically assess the frequency of rare transcription factor variants in sporadic patients with the cardiac outflow tract malformation tetralogy of Fallot (TOF). METHODS AND RESULTS:We sequenced the coding, 5'UTR, and 3'UTR regions of twelve transcription factor genes implicated in cardiac outflow tract development (NKX2.5, GATA4, ISL1, TBX20, MEF2C, BOP/SMYD1, HAND2, FOXC1, FOXC2, FOXH, FOXA2 and TBX1) in 93 non-syndromic, non-Mendelian TOF cases. We also analysed Illumina Human 660W-Quad SNP Array data for copy number variants in these genes; none were detected. Four of the rare variants detected have previously been shown to affect transactivation in in vitro reporter assays: FOXC1 p.P297S, FOXC2 p.Q444R, FOXH1 p.S113T and TBX1 p.P43_G61del PPPPRYDPCAAAAPGAPGP. Two further rare variants, HAND2 p.A25_A26insAA and FOXC1 p.G378_G380delGGG, A488_491delAAAA, affected transactivation in in vitro reporter assays. Each of these six functionally significant variants was present in a single patient in the heterozygous state; each of the four for which parental samples were available were maternally inherited. Thus in the 93 TOF cases we identified six functionally significant mutations in the secondary heart field transcriptional network. SIGNIFICANCE:This study indicates that rare genetic variants in the secondary heart field transcriptional network with functional effects on protein function occur in 3-13% of patients with TOF. This is the first report of a functionally significant HAND2 mutation in a patient with congenital heart disease.