Comprehensive reanalysis for CNVs in ES data from unsolved rare disease cases results in new diagnoses.

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.

Increased Diagnostic Yield by Reanalysis of Whole Exome Sequencing Data in Mitochondrial Disease.

Background: The genetic diagnosis of mitochondrial disorders is complicated by its genetic and phenotypic complexity. Next generation sequencing techniques have much improved the diagnostic yield for these conditions. A cohort of individuals with multiple respiratory chain deficiencies, reported in the literature 10 years ago, had a diagnostic rate of 60% by whole exome sequencing (WES) but 40% remained undiagnosed. Objective: We aimed to identify a genetic diagnosis by reanalysis of the WES data for the undiagnosed arm of this 10-year-old cohort of patients with suspected mitochondrial disorders. Methods: The WES data was transferred and processed by the RD-Connect Genome-Phenome Analysis Platform (GPAP) using their standardized pipeline. Variant prioritisation was carried out on the RD-Connect GPAP. Results: Singleton WES data from 14 individuals was reanalysed. We identified a possible or likely genetic diagnosis in 8 patients (8/14, 57%). The variants identified were in a combination of mitochondrial DNA (n = 1, MT-TN), nuclear encoded mitochondrial genes (n = 2, PDHA1, and SUCLA2) and nuclear genes associated with nonmitochondrial disorders (n = 5, PNPLA2, CDC40, NBAS and SLC7A7). Variants in both the NBAS and CDC40 genes were established as disease causing after the original cohort was published. We increased the diagnostic yield for the original cohort by 15% without generating any further genomic data. Conclusions: In the era of multiomics we highlight that reanalysis of existing WES data is a valid tool for generating additional diagnosis in patients with suspected mitochondrial disease, particularly when more time has passed to allow for new bioinformatic pipelines to emerge, for the development of new tools in variant interpretation aiding in reclassification of variants and the expansion of scientific knowledge on additional genes.

Variants in mitochondrial disease genes are common causes of inherited peripheral neuropathies.

BACKGROUND: Peripheral neuropathies in mitochondrial disease are caused by mutations in nuclear genes encoding mitochondrial proteins, or in the mitochondrial genome. Whole exome or genome sequencing enable parallel testing of nuclear and mtDNA genes, and it has significantly advanced the genetic diagnosis of inherited diseases. Despite this, approximately 40% of all Charcot-Marie-Tooth (CMT) cases remain undiagnosed. METHODS: The genome-phenome analysis platform (GPAP) in RD-Connect was utilised to create a cohort of 2087 patients with at least one Human Phenotype Ontology (HPO) term suggestive of a peripheral neuropathy, from a total of 10,935 patients. These patients' genetic data were then analysed and searched for variants in known mitochondrial disease genes. RESULTS: A total of 1,379 rare variants were identified, 44 of which were included in this study as either reported pathogenic or likely causative in 42 patients from 36 families. The most common genes found to be likely causative for an autosomal dominant neuropathy were GDAP1 and GARS1. We also detected heterozygous likely pathogenic variants in DNA2, MFN2, DNM2, PDHA1, SDHA, and UCHL1. Biallelic variants in SACS, SPG7, GDAP1, C12orf65, UCHL1, NDUFS6, ETFDH and DARS2 and variants in the mitochondrial DNA (mtDNA)-encoded MT-ATP6 and MT-TK were also causative for mitochondrial CMT. Only 50% of these variants were already reported as solved in GPAP. CONCLUSION: Variants in mitochondrial disease genes are frequent in patients with inherited peripheral neuropathies. Due to the clinical overlap between mitochondrial disease and CMT, agnostic exome or genome sequencing have better diagnostic yields than targeted gene panels.

Phenotypic similarity-based approach for variant prioritization for unsolved rare disease: a preliminary methodological report.

Rare diseases (RD) have a prevalence of not more than 1/2000 persons in the European population, and are characterised by the difficulty experienced in obtaining a correct and timely diagnosis. According to Orphanet, 72.5% of RD have a genetic origin although 35% of them do not yet have an identified causative gene. A significant proportion of patients suspected to have a genetic RD receive an inconclusive exome/genome sequencing. Working towards the International Rare Diseases Research Consortium (IRDiRC)'s goal for 2027 to ensure that all people living with a RD receive a diagnosis within one year of coming to medical attention, the Solve-RD project aims to identify the molecular causes underlying undiagnosed RD. As part of this strategy, we developed a phenotypic similarity-based variant prioritization methodology comparing submitted cases with other submitted cases and with known RD in Orphanet. Three complementary approaches based on phenotypic similarity calculations using the Human Phenotype Ontology (HPO), the Orphanet Rare Diseases Ontology (ORDO) and the HPO-ORDO Ontological Module (HOOM) were developed; genomic data reanalysis was performed by the RD-Connect Genome-Phenome Analysis Platform (GPAP). The methodology was tested in 4 exemplary cases discussed with experts from European Reference Networks. Variants of interest (pathogenic or likely pathogenic) were detected in 8.8% of the 725 cases clustered by similarity calculations. Diagnostic hypotheses were validated in 42.1% of them and needed further exploration in another 10.9%. Based on the promising results, we are devising an automated standardized phenotypic-based re-analysis pipeline to be applied to the entire unsolved cases cohort.

Remote visualization of large-scale genomic alignments for collaborative clinical research and diagnosis of rare diseases.

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.

A method for multiplexed full-length single-molecule sequencing of the human mitochondrial genome.

Methods to reconstruct the mitochondrial DNA (mtDNA) sequence using short-read sequencing come with an inherent bias due to amplification and mapping. They can fail to determine the phase of variants, to capture multiple deletions and to cover the mitochondrial genome evenly. Here we describe a method to target, multiplex and sequence at high coverage full-length human mitochondrial genomes as native single-molecules, utilizing the RNA-guided DNA endonuclease Cas9. Combining Cas9 induced breaks, that define the mtDNA beginning and end of the sequencing reads, as barcodes, we achieve high demultiplexing specificity and delineation of the full-length of the mtDNA, regardless of the structural variant pattern. The long-read sequencing data is analysed with a pipeline where our custom-developed software, baldur, efficiently detects single nucleotide heteroplasmy to below 1%, physically determines phase and can accurately disentangle complex deletions. Our workflow is a tool for studying mtDNA variation and will accelerate mitochondrial research.

The RD-Connect Genome-Phenome Analysis Platform: Accelerating diagnosis, research, and gene discovery for rare diseases.

Rare disease patients are more likely to receive a rapid molecular diagnosis nowadays thanks to the wide adoption of next-generation sequencing. However, many cases remain undiagnosed even after exome or genome analysis, because the methods used missed the molecular cause in a known gene, or a novel causative gene could not be identified and/or confirmed. To address these challenges, the RD-Connect Genome-Phenome Analysis Platform (GPAP) facilitates the collation, discovery, sharing, and analysis of standardized genome-phenome data within a collaborative environment. Authorized clinicians and researchers submit pseudonymised phenotypic profiles encoded using the Human Phenotype Ontology, and raw genomic data which is processed through a standardized pipeline. After an optional embargo period, the data are shared with other platform users, with the objective that similar cases in the system and queries from peers may help diagnose the case. Additionally, the platform enables bidirectional discovery of similar cases in other databases from the Matchmaker Exchange network. To facilitate genome-phenome analysis and interpretation by clinical researchers, the RD-Connect GPAP provides a powerful user-friendly interface and leverages tens of information sources. As a result, the resource has already helped diagnose hundreds of rare disease patients and discover new disease causing genes.

High diagnostic rate of trio exome sequencing in consanguineous families with neurogenetic diseases.

Consanguineous marriages have a prevalence rate of 24% in Turkey. These carry an increased risk of autosomal recessive genetic conditions, leading to severe disability or premature death, with a significant health and economic burden. A definitive molecular diagnosis could not be achieved in these children previously, as infrastructures and access to sophisticated diagnostic options were limited. We studied the cause of neurogenetic disease in 246 children from 190 consanguineous families recruited in three Turkish hospitals between 2016 and 2020. All patients underwent deep phenotyping and trio whole exome sequencing, and data were integrated in advanced international bioinformatics platforms. We detected causative variants in 119 known disease genes in 72% of families. Due to overlapping phenotypes 52% of the confirmed genetic diagnoses would have been missed on targeted diagnostic gene panels. Likely pathogenic variants in 27 novel genes in 14% of the families increased the diagnostic yield to 86%. Eighty-two per cent of causative variants (141/172) were homozygous, 11 of which were detected in genes previously only associated with autosomal dominant inheritance. Eight families carried two pathogenic variants in different disease genes. De novo (9.3%), X-linked recessive (5.2%) and compound heterozygous (3.5%) variants were less frequent compared to non-consanguineous populations. This cohort provided a unique opportunity to better understand the genetic characteristics of neurogenetic diseases in a consanguineous population. Contrary to what may be expected, causative variants were often not on the longest run of homozygosity and the diagnostic yield was lower in families with the highest degree of consanguinity, due to the high number of homozygous variants in these patients. Pathway analysis highlighted that protein synthesis/degradation defects and metabolic diseases are the most common pathways underlying paediatric neurogenetic disease. In our cohort 164 families (86%) received a diagnosis, enabling prevention of transmission and targeted treatments in 24 patients (10%). We generated an important body of genomic data with lasting impacts on the health and wellbeing of consanguineous families and economic benefit for the healthcare system in Turkey and elsewhere. We demonstrate that an untargeted next generation sequencing approach is far superior to a more targeted gene panel approach, and can be performed without specialized bioinformatics knowledge by clinicians using established pipelines in populations with high rates of consanguinity.

A systematic study and literature review of parental somatic mosaicism of FBN1 pathogenic variants in Marfan syndrome.

BACKGROUND: A proportion of de novo variants in patients affected by genetic disorders, particularly those with autosomal dominant (AD) inheritance, could be the consequence of somatic mosaicism in one of the progenitors. There is growing evidence that germline and somatic mosaicism are more common and play a greater role in genetic disorders than previously acknowledged. In Marfan syndrome (MFS), caused by pathogenic variants in the fibrillin-1 gene (FBN1) gene, approximately 25% of the disease-causing variants are reported as de novo. Only a few cases of parental mosaicism have been reported in MFS. METHODS: Employing an amplicon-based deep sequencing (ADS) method, we carried out a systematic analysis of 60 parents of 30 FBN1 positive, consecutive patients with MFS with an apparently de novo pathogenic variant. RESULTS: Out of the 60 parents studied (30 families), the majority (n=51, 85%) had a systemic score of 0, seven had a score of 1 and two a score of 2, all due to minor criteria common in the normal population. We detected two families with somatic mosaicism in one of the progenitors, with a rate of 6.6% (2/30) of apparently de novo cases. CONCLUSIONS: The search for parental somatic mosaicism should be routinely implemented in de novo cases of MFS, to offer appropriate genetic and reproductive counselling as well as to reveal masked, isolated clinical signs of MFS in progenitors that may require specific follow-up.

A MT-TL1 variant identified by whole exome sequencing in an individual with intellectual disability, epilepsy, and spastic tetraparesis.

The genetic etiology of intellectual disability remains elusive in almost half of all affected individuals. Within the Solve-RD consortium, systematic re-analysis of whole exome sequencing (WES) data from unresolved cases with (syndromic) intellectual disability (n = 1,472 probands) was performed. This re-analysis included variant calling of mitochondrial DNA (mtDNA) variants, although mtDNA is not specifically targeted in WES. We identified a functionally relevant mtDNA variant in MT-TL1 (NC_012920.1:m.3291T > C; NC_012920.1:n.62T > C), at a heteroplasmy level of 22% in whole blood, in a 23-year-old male with severe intellectual disability, epilepsy, episodic headaches with emesis, spastic tetraparesis, brain abnormalities, and feeding difficulties. Targeted validation in blood and urine supported pathogenicity, with heteroplasmy levels of 23% and 58% in index, and 4% and 17% in mother, respectively. Interestingly, not all phenotypic features observed in the index have been previously linked to this MT-TL1 variant, suggesting either broadening of the m.3291T > C-associated phenotype, or presence of a co-occurring disorder. Hence, our case highlights the importance of underappreciated mtDNA variants identifiable from WES data, especially for cases with atypical mitochondrial phenotypes and their relatives in the maternal line.

Beyond copy number: A new, rapid, and versatile method for sequencing the entire SMN2 gene in SMA patients.

Spinal muscular atrophy (SMA) is caused by bi-allelic loss or pathogenic variants in the SMN1 gene. SMN2, the highly homologous copy of SMN1, is considered the major phenotypic modifier of the disease. Determination of SMN2 copy number is essential to establish robust genotype-phenotype correlations and predict disease evolution, to stratify patients for clinical trials, as well as to define those eligible for treatment. Discordant genotype-phenotype correlations are not uncommon in SMA, some of which are due to intragenic SMN2 variants that may influence the amount of complete SMN transcripts and, therefore, of full-length SMN protein. Detection of these variants is crucial to predict SMA phenotypes in the present scenario of therapeutic advances and with the perspective of SMA neonatal screening and early diagnosis to start treatments. Here, we present a novel, affordable, and versatile method for complete sequencing of the SMN2 gene based on long-range polymerase chain reaction and next-generation sequencing. The method was validated by analyzing samples from 53 SMA patients who lack SMN1, allowing to characterize paralogous, rare variants, and single-nucleotide polymorphisms of SMN2 as well as SMN2-SMN1 hybrid genes. The method identifies partial deletions and can be adapted to determine rare pathogenic variants in patients with at least one SMN1 copy.

Implementation of a Gene Panel for Genetic Diagnosis of Primary Ciliary Dyskinesia. / Implementación de un panel de genes para el diagnóstico genético de la discinesia ciliar primaria.

INTRODUCTION: Primary ciliary dyskinesia (PCD) is characterized by an alteration in the ciliary structure causing difficulty in the clearance of respiratory secretions. Diagnosis is complex and based on a combination of techniques. The objective of this study was to design a gene panel including all known causative genes, and to corroborate their diagnostic utility in a cohort of Spanish patients. METHODS: This was a multicenter cross-sectional study of patients with a high suspicion of PCD, according to European Respiratory Society criteria, designed around a gene panel for massive sequencing using SeqCap EZ capture technology that included 44 genes associated with PCD. RESULTS: We included 79 patients, 53 of whom had a diagnosis of confirmed or highly probable PCD. The sensitivity of the gene panel was 81.1%, with a specificity of 100%. Candidate variants were found in some of the genes of the panel in 43 patients with PCD, 51.2% (22/43) of whom were homozygotes and 48.8% (21/43) compound heterozygotes. The most common causative genes were DNAH5 and CCDC39. We found 52 different variants, 36 of which were not previously described in the literature. CONCLUSIONS: The design and implementation of a tailored gene panel produces a high yield in the genetic diagnosis of PCD. This panel provides a better understanding of the causative factors involved in these patients and lays down the groundwork for future therapeutic approaches.

Phenotypic Variability of Patients With PAX8 Variants Presenting With Congenital Hypothyroidism and Eutopic Thyroid.

PURPOSE: Thyroid dyshormonogenesis is a heterogeneous group of hereditary diseases produced by a total/partial blockage of the biochemical processes of thyroid-hormone synthesis and secretion. Paired box 8 (PAX8) is essential for thyroid morphogenesis and thyroid hormone synthesis. We aimed to identify PAX8 variants in patients with thyroid dyshormonogenesis and to analyze them with in vitro functional studies. PATIENTS AND METHODS: Nine pediatric patients with a eutopic thyroid gland were analyzed by the Catalan screening program for congenital hypothyroidism. Scintigraphies showed absent, low, or normal uptake. Only one patient had a hypoplastic gland. On reevaluation, perchlorate discharge test was negative or compatible with partial iodine-organization deficit. After evaluation, 8 patients showed permanent mild or severe hypothyroidism. Massive-sequencing techniques were used to detect variants in congenital hypothyroidism-related genes. In vitro functional studies were based on transactivating activity of mutant PAX8 on a TG-gene promoter and analyzed by a dual-luciferase assays. RESULTS: We identified 7 heterozygous PAX8 exonic variants and 1 homozygous PAX8 splicing variant in 9 patients with variable phenotypes of thyroid dyshormonogenesis. Five were novel and 5 variants showed a statistically significant impaired transcriptional activity of TG promoter: 51% to 78% vs the wild type. CONCLUSIONS: Nine patients presented with PAX8 candidate variants. All presented with a eutopic thyroid gland and 7 had deleterious variants. The phenotype of affected patients varies considerably, even within the same family; but, all except the homozygous patient presented with a normal eutopic thyroid gland and thyroid dyshormonogenesis. PAX8 functional studies have shown that 6 PAX8 variants are deleterious. Our studies have proven effective in evaluating these variants.

Uncovering Low-Level Maternal Gonosomal Mosaicism in X-Linked Agammaglobulinemia: Implications for Genetic Counseling.

X-linked agammaglobulinemia (XLA) is a clinically and genetically well-defined immunodeficiency and the most common form of agammaglobulinemia. It is characterized by susceptibility to recurrent bacterial infections, profound hypogammaglobulinemia, and few or no circulating B cells. XLA is caused by mutations in the BTK gene, which encodes Bruton's tyrosine kinase (BTK). Because of its X-linked recessive inheritance pattern, XLA virtually only affects males, and the mother is the carrier of the mutation in 80-85% of the males with this condition. In the remaining 15-20% of the cases, the affected male is considered to have a de novo mutation. Here, we present the case of a child with a diagnosis of XLA caused by a missense mutation in the BTK gene (c.494G>A/p.C165Y). Apparently, his mother was wild type for this gene, which implied that the mutation was de novo, but careful analysis of Sanger electropherograms and the use of high-coverage massive parallel sequencing revealed low-level maternal gonosomal mosaicism. The mutation was detected in various samples from the mother (blood, urine, buccal swab, and vaginal swab) at a low frequency of 2-5%, and the status of the patient's mutation changed from de novo to inherited. This study underscores the importance of accurately establishing the parents' status on detection of an apparently de novo mutation in a patient, as inadvertent low-level mosaicism may lead to misinterpretation of the risk of recurrence, vital for genetic counseling.

Severe brain involvement in 5q spinal muscular atrophy type 0.

Spinal muscular atrophy (SMA) type 0 is the most severe form of SMA, associated with the SMN1 gene and manifesting at birth. Most patients die in the first weeks of life. In this work, we present 3 patients with SMA type 0 who survived >1 year and presented diffuse and progressive brain abnormalities on magnetic resonance imaging, which are not usually seen in patients with SMA. Thus, severe brain involvement may likely be the full end manifestation of an already extreme SMA phenotype caused by substantial reduction of the SMN protein in the brain. ANN NEUROL 2019;86:458-462.

Th1-skewed profile and excessive production of proinflammatory cytokines in a NFKB1-deficient patient with CVID and severe gastrointestinal manifestations.

Monoallelic loss-of-function mutations in NFKB1 were recently recognized as the most common monogenic cause of common variable immunodeficiency (CVID). The prototypic clinical phenotype of NFKB1-deficient patients includes common CVID features, such as hypogammaglobulinaemia and sinopulmonary infections, plus other highly variable individual manifestations. Here, we describe a patient with a profound CVID phenotype and severe gastrointestinal manifestations, including chronic and recurrent diarrhoea. Using an NGS customized panel of 323 genes related to primary immunodeficiencies, we identified a novel monoallelic loss-of-function mutation in NFKB1 leading to a truncated protein (c.1149delT/p.Gly384Glu ∗ 48). Interestingly, we also found a rare variant in NOD2 previously associated with Crohn's disease (p.His352Arg). Our patient had hypogammaglobulinaemia with a small number of B cells, most of which were naïve. The most noteworthy findings included marked skewing towards a Th1 phenotype in peripheral blood T cells and excessive production of proinflammatory cytokines (IL-1ß, TNFα). The patient's 6-year-old daughter, a carrier of the NFKB1 mutation, is clinically asymptomatic, but has started to show cellular and molecular changes. This case of NFKB1 deficiency appears to be a combination of immunodeficiency and a hyperinflammatory state. The current situation of the patient's daughter provides a glimpse of the preclinical phase of the condition.

A common structural scaffold in CTD phosphatases that supports distinct catalytic mechanisms.

The phosphorylation and dephosphorylation of the carboxyl-terminal domain (CTD) of the largest RNA polymerase II (RNAPII) subunit is a critical regulatory checkpoint for transcription and mRNA processing. This CTD is unique to eukaryotic organisms and it contains multiple tandem-repeats with the consensus sequence Tyr(1) -Ser(2) -Pro(3) -Thr(4) -Ser(5) -Pro(6) -Ser(7) . Traditionally, CTD phosphatases that use metal-ion-independent (cysteine-based) and metal-ion-assisted (aspartate-based) catalytic mechanisms have been considered to belong to two independent groups. However, using structural comparisons we have identified a common structural scaffold in these two groups of CTD phosphatases. This common scaffold accommodates different catalytic processes with the same substrate specificity, in this case phospho-serine/threonine residues flanked by prolines. Furthermore, this scaffold provides a structural connection between two groups of protein tyrosine phosphatases (PTPs): Cys-based (classes I, II, and III) and Asp-based (class IV) PTPs. Redundancy in catalytic mechanisms is not infrequent and may arise in specific biological settings. To better understand the activity of the CTD phosphatases, we combined our structural analyses with data on CTD phosphatase expression in different human and mouse tissues. The results suggest that aspartate- and cysteine-based CTD-dephosphorylation acts in concert during cellular stress, when high levels of reactive oxygen species can inhibit the nucleophilic function of the catalytic cysteine, as occurs in mental and neurodegenerative disorders like schizophrenia, Alzheimer's and Parkinson's diseases. Moreover, these findings have significant implications for the study of the RNAPII-CTD dephosphorylation in eukaryotes.

Cord blood-derived neuronal cells by ectopic expression of Sox2 and c-Myc.

The finding that certain somatic cells can be directly converted into cells of other lineages by the delivery of specific sets of transcription factors paves the way to novel therapeutic applications. Here we show that human cord blood (CB) CD133(+) cells lose their hematopoietic signature and are converted into CB-induced neuronal-like cells (CB-iNCs) by the ectopic expression of the transcription factor Sox2, a process that is further augmented by the combination of Sox2 and c-Myc. Gene-expression analysis, immunophenotyping, and electrophysiological analysis show that CB-iNCs acquire a distinct neuronal phenotype characterized by the expression of multiple neuronal markers. CB-iNCs show the ability to fire action potentials after in vitro maturation as well as after in vivo transplantation into the mouse hippocampus. This system highlights the potential of CB cells and offers an alternative means to the study of cellular plasticity, possibly in the context of drug screening research and of future cell-replacement therapies.