De novo and inherited monoallelic variants in TUBA4A cause ataxia and spasticity.

Alpha-tubulin 4A encoding gene (TUBA4A) has been associated with familial amyotrophic lateral sclerosis (fALS) and fronto-temporal dementia (FTD), based on identification of likely pathogenic variants in patients from distinct ALS and FTD cohorts. By screening a multicentric French cohort of 448 unrelated probands presenting with cerebellar ataxia, we identified ultra-rare TUBA4A missense variants, all being absent from public databases and predicted pathogenic by multiple in-silico tools. In addition, gene burden analyses in the 100,000 genomes project (100KGP) showed enrichment of TUBA4A rare variants in the inherited ataxia group compared to controls (OR: 57.0847 [10.2- 576.7]; p = 4.02 x10-07). Altogether, we report 12 patients presenting with spasticity and/or cerebellar ataxia and harboring a predicted pathogenic TUBA4A missense mutation, including 5 confirmed de novo cases and a mutation previously reported in a large family presenting with spastic ataxia. Cultured fibroblasts from 3 patients harboring distinct TUBA4A missense showed significant alterations in microtubule organisation and dynamics, providing insight of TUBA4A variants pathogenicity. Our data confirm the identification of a hereditary spastic ataxia disease gene with variable age of onset, expanding the clinical spectrum of TUBA4A associated phenotypes.

SpliceAI-visual: a free online tool to improve SpliceAI splicing variant interpretation.

SpliceAI is an open-source deep learning splicing prediction algorithm that has demonstrated in the past few years its high ability to predict splicing defects caused by DNA variations. However, its outputs present several drawbacks: (1) although the numerical values are very convenient for batch filtering, their precise interpretation can be difficult, (2) the outputs are delta scores which can sometimes mask a severe consequence, and (3) complex delins are most often not handled. We present here SpliceAI-visual, a free online tool based on the SpliceAI algorithm, and show how it complements the traditional SpliceAI analysis. First, SpliceAI-visual manipulates raw scores and not delta scores, as the latter can be misleading in certain circumstances. Second, the outcome of SpliceAI-visual is user-friendly thanks to the graphical presentation. Third, SpliceAI-visual is currently one of the only SpliceAI-derived implementations able to annotate complex variants (e.g., complex delins). We report here the benefits of using SpliceAI-visual and demonstrate its relevance in the assessment/modulation of the PVS1 classification criteria. We also show how SpliceAI-visual can elucidate several complex splicing defects taken from the literature but also from unpublished cases. SpliceAI-visual is available as a Google Colab notebook and has also been fully integrated in a free online variant interpretation tool, MobiDetails ( https://mobidetails.iurc.montp.inserm.fr/MD ).

RFC1 nonsense and frameshift variants cause CANVAS: clues for an unsolved pathophysiology.

Cerebellar ataxia, neuropathy and vestibular areflexia syndrome (CANVAS) is an inherited late-onset neurological disease caused by bi-allelic AAGGG pentanucleotide expansions within intron 2 of RFC1. Despite extensive studies, the pathophysiological mechanism of these intronic expansions remains elusive. We screened by clinical exome sequencing two unrelated patients presenting with late-onset ataxia. A repeat-primer polymerase chain reaction was used for RFC1 AAGGG intronic expansion identification. RFC1 mRNA expression was assessed by quantitative reverse transcription-polymerase chain reaction. We identified the first two CANVAS affected patients who are compound heterozygous for RFC1 truncating variants (p.Arg388* and c.575delA, respectively) and a pathological AAGGG expansion. RFC1 expression studies in whole blood showed a significant reduction of RFC1 mRNA for both patients compared to three patients with bi-allelic RFC1 expansions. In conclusion, this observation provides clues that suggest bi-allelic RFC1 conditional loss-of-function as the cause of the disease.

Evaluating the Transition from Targeted to Exome Sequencing: A Guide for Clinical Laboratories.

The transition from targeted to exome or genome sequencing in clinical contexts requires quality standards, such as targeted sequencing, in order to be fully adopted. However, no clear recommendations or methodology have emerged for evaluating this technological evolution. We developed a structured method based on four run-specific sequencing metrics and seven sample-specific sequencing metrics for evaluating the performance of exome sequencing strategies to replace targeted strategies. The indicators include quality metrics and coverage performance on gene panels and OMIM morbid genes. We applied this general strategy to three different exome kits and compared them with a myopathy-targeted sequencing method. After having achieved 80 million reads, all-tested exome kits generated data suitable for clinical diagnosis. However, significant differences in the coverage and PCR duplicates were observed between the kits. These are two main criteria to consider for the initial implementation with high-quality assurance. This study aims to assist molecular diagnostic laboratories in adopting and evaluating exome sequencing kits in a diagnostic context compared to the strategy used previously. A similar strategy could be used to implement whole-genome sequencing for diagnostic purposes.

High rate of hypomorphic variants as the cause of inherited ataxia and related diseases: study of a cohort of 366 families.

PURPOSE: Diagnosis of inherited ataxia and related diseases represents a real challenge given the tremendous heterogeneity and clinical overlap of the various causes. We evaluated the efficacy of molecular diagnosis of these diseases by sequencing a large cohort of undiagnosed families. METHODS: We analyzed 366 unrelated consecutive patients with undiagnosed ataxia or related disorders by clinical exome-capture sequencing. In silico analysis was performed with an in-house pipeline that combines variant ranking and copy-number variant (CNV) searches. Variants were interpreted according to American College of Medical Genetics and Genomics/Association for Molecular Pathology (ACMG/AMP) guidelines. RESULTS: We established the molecular diagnosis in 46% of the cases. We identified 35 mildly affected patients with causative variants in genes that are classically associated with severe presentations. These cases were explained by the occurrence of hypomorphic variants, but also rarely suspected mechanisms such as C-terminal truncations and translation reinitiation. CONCLUSION: A significant fraction of the clinical heterogeneity and phenotypic overlap is explained by hypomorphic variants that are difficult to identify and not readily predicted. The hypomorphic C-terminal truncation and translation reinitiation mechanisms that we identified may only apply to few genes, as it relies on specific domain organization and alterations. We identified PEX10 and FASTKD2 as candidates for translation reinitiation accounting for mild disease presentation.

SLC12A2 variants cause a neurodevelopmental disorder or cochleovestibular defect.

The SLC12 gene family consists of SLC12A1-SLC12A9, encoding electroneutral cation-coupled chloride co-transporters. SCL12A2 has been shown to play a role in corticogenesis and therefore represents a strong candidate neurodevelopmental disorder gene. Through trio exome sequencing we identified de novo mutations in SLC12A2 in six children with neurodevelopmental disorders. All had developmental delay or intellectual disability ranging from mild to severe. Two had sensorineural deafness. We also identified SLC12A2 variants in three individuals with non-syndromic bilateral sensorineural hearing loss and vestibular areflexia. The SLC12A2 de novo mutation rate was demonstrated to be significantly elevated in the deciphering developmental disorders cohort. All tested variants were shown to reduce co-transporter function in Xenopus laevis oocytes. Analysis of SLC12A2 expression in foetal brain at 16-18 weeks post-conception revealed high expression in radial glial cells, compatible with a role in neurogenesis. Gene co-expression analysis in cells robustly expressing SLC12A2 at 16-18 weeks post-conception identified a transcriptomic programme associated with active neurogenesis. We identify SLC12A2 de novo mutations as the cause of a novel neurodevelopmental disorder and bilateral non-syndromic sensorineural hearing loss and provide further data supporting a role for this gene in human neurodevelopment.

The Study of a 231 French Patient Cohort Significantly Extends the Mutational Spectrum of the Two Major Usher Genes MYO7A and USH2A.

Usher syndrome is an autosomal recessive disorder characterized by congenital hearing loss combined with retinitis pigmentosa, and in some cases, vestibular areflexia. Three clinical subtypes are distinguished, and MYO7A and USH2A represent the two major causal genes involved in Usher type I, the most severe form, and type II, the most frequent form, respectively. Massively parallel sequencing was performed on a cohort of patients in the context of a molecular diagnosis to confirm clinical suspicion of Usher syndrome. We report here 231 pathogenic MYO7A and USH2A genotypes identified in 73 Usher type I and 158 Usher type II patients. Furthermore, we present the ACMG classification of the variants, which comprise all types. Among them, 68 have not been previously reported in the literature, including 12 missense and 16 splice variants. We also report a new deep intronic variant in USH2A. Despite the important number of molecular studies published on these two genes, we show that during the course of routine genetic diagnosis, undescribed variants continue to be identified at a high rate. This is particularly pertinent in the current era, where therapeutic strategies based on DNA or RNA technologies are being developed.

The CYSMA web server: An example of integrative tool for in silico analysis of missense variants identified in Mendelian disorders.

Exome sequencing used for molecular diagnosis of Mendelian disorders considerably increases the number of missense variants of unclear significance, whose pathogenicity can be assessed by a variety of prediction tools. As the performance of algorithms may vary according to the datasets, complementary specific resources are needed to improve variant interpretation. As a model, we were interested in the cystic fibrosis transmembrane conductance regulator gene (CFTR) causing cystic fibrosis, in which at least 40% of missense variants are reported. Cystic fibrosis missense analysis (CYSMA) is a new web server designed for online estimation of the pathological relevance of CFTR missense variants. CYSMA generates a set of computationally derived data, ranging from evolutionary conservation to functional observations from three-dimensional structures, provides all available allelic frequencies, clinical observations, and references for functional studies. Compared to software classically used in analysis pipelines on a dataset of 141 well-characterized missense variants, CYSMA was the most efficient tool to discriminate benign missense variants, with a specificity of 85%, and very good sensitivity of 89%. These results suggest that such integrative tools could be adapted to numbers of genes involved in Mendelian disorders to improve the interpretation of missense variants identified in the context of diagnosis.

Pathogenicity of novel atypical variants leading to choroideremia as determined by functional analyses.

Choroideremia is a monogenic X-linked recessive chorioretinal disease linked to pathogenic variants in the CHM gene. These variants are commonly base-pair changes, frameshifts, or large deletions. However, a few rare or unusual events comprising large duplications, a retrotransposon insertion, a pseudo-exon activation, and two c-98 promoter substitutions have also been described. Following an exhaustive molecular diagnosis, we identified and characterized three novel atypical disease-causing variants in three unrelated male patients. One is a first-ever reported Alu insertion within CHM and the other two are nucleotide substitutions, c.-90C>G and c.-108A>G, affecting highly conserved promoter positions. RNA analysis combined with western blot and functional assays of patient cells established the pathogenicity of the Alu insertion and the c.-90C>G alteration. Furthermore, luciferase reporter assays suggested a CHM transcription defect associated with the c.-90C>G and c.-108A>G variants. These findings broaden our knowledge of the mutational spectrum and the transcriptional regulation of the CHM gene.

CFTR-France, a national relational patient database for sharing genetic and phenotypic data associated with rare CFTR variants.

Most of the 2,000 variants identified in the CFTR (cystic fibrosis transmembrane regulator) gene are rare or private. Their interpretation is hampered by the lack of available data and resources, making patient care and genetic counseling challenging. We developed a patient-based database dedicated to the annotations of rare CFTR variants in the context of their cis- and trans-allelic combinations. Based on almost 30 years of experience of CFTR testing, CFTR-France (https://cftr.iurc.montp.inserm.fr/cftr) currently compiles 16,819 variant records from 4,615 individuals with cystic fibrosis (CF) or CFTR-RD (related disorders), fetuses with ultrasound bowel anomalies, newborns awaiting clinical diagnosis, and asymptomatic compound heterozygotes. For each of the 736 different variants reported in the database, patient characteristics and genetic information (other variations in cis or in trans) have been thoroughly checked by a dedicated curator. Combining updated clinical, epidemiological, in silico, or in vitro functional data helps to the interpretation of unclassified and the reassessment of misclassified variants. This comprehensive CFTR database is now an invaluable tool for diagnostic laboratories gathering information on rare variants, especially in the context of genetic counseling, prenatal and preimplantation genetic diagnosis. CFTR-France is thus highly complementary to the international database CFTR2 focused so far on the most common CF-causing alleles.

Whole USH2A Gene Sequencing Identifies Several New Deep Intronic Mutations.

Deep intronic mutations leading to pseudoexon (PE) insertions are underestimated and most of these splicing alterations have been identified by transcript analysis, for instance, the first deep intronic mutation in USH2A, the gene most frequently involved in Usher syndrome type II (USH2). Unfortunately, analyzing USH2A transcripts is challenging and for 1.8%-19% of USH2 individuals carrying a single USH2A recessive mutation, a second mutation is yet to be identified. We have developed and validated a DNA next-generation sequencing approach to identify deep intronic variants in USH2A and evaluated their consequences on splicing. Three distinct novel deep intronic mutations have been identified. All were predicted to affect splicing and resulted in the insertion of PEs, as shown by minigene assays. We present a new and attractive strategy to identify deep intronic mutations, when RNA analyses are not possible. Moreover, the bioinformatics pipeline developed is independent of the gene size, implying the possible application of this approach to any disease-linked gene. Finally, an antisense morpholino oligonucleotide tested in vitro for its ability to restore splicing caused by the c.9959-4159A>G mutation provided high inhibition rates, which are indicative of its potential for molecular therapy.

Enrichment of LOVD-USHbases with 152 USH2A genotypes defines an extensive mutational spectrum and highlights missense hotspots.

Alterations of USH2A, encoding usherin, are responsible for more than 70% of cases of Usher syndrome type II (USH2), a recessive disorder that combines moderate to severe hearing loss and retinal degeneration. The longest USH2A transcript encodes usherin isoform b, a 5,202-amino-acid transmembrane protein with an exceptionally large extracellular domain consisting notably of a Laminin N-terminal domain and numerous Laminin EGF-like (LE) and Fibronectin type III (FN3) repeats. Mutations of USH2A are scattered throughout the gene and mostly private. Annotating these variants is therefore of major importance to correctly assign pathogenicity. We have extensively genotyped a novel cohort of 152 Usher patients and identified 158 different mutations, of which 93 are newly described. Pooling this new data with the existing pathogenic variants already incorporated in USHbases reveals several previously unappreciated features of the mutational spectrum. We show that parts of the protein are more likely to tolerate single amino acid variations, whereas others constitute pathogenic missense hotspots. We have found, in repeated LE and FN3 domains, a nonequal distribution of the missense mutations that highlights some crucial positions in usherin with possible consequences for the assessment of the pathogenicity of the numerous missense variants identified in USH2A.

Validation of Nanopore long-read sequencing to resolve RPGR ORF15 genotypes in individuals with X-linked retinitis pigmentosa.

X-linked retinitis pigmentosa (XLRP) is characterized by progressive vision loss leading to legal blindness in males and a broad severity spectrum in carrier females. Pathogenic alterations of the retinitis pigmentosa GTPase regulator gene (RPGR) are responsible for over 70% of XLRP cases. In the retina, the RPGRORF15 transcript includes a terminal exon, called ORF15, that is altered in the large majority of RPGR-XLRP cases. Unfortunately, due to its highly repetitive sequence, ORF15 represents a considerable challenge in terms of sequencing for molecular diagnostic laboratories. However, in a recent preliminary work Yahya et al. reported a long-read sequencing approach seeming promising. Here, the aim of the study was to validate and integrate this new sequencing strategy in a routine screening workflow. For that purpose, we performed a masked test on 52 genomic DNA samples from male and female individuals carrying 32 different pathogenic ORF15 variations including 20 located in the highly repetitive region of the exon. For the latter, we have obtained a detection rate of 80-85% in males and 60-80% in females after bioinformatic analyses. These numbers raised to 100% for both status after adding a complementary visual inspection of ORF15 long-reads. In accordance with these results, and considering the frequency of ORF15 pathogenic variations in XLRP, we suggest that a long-read screening of ORF15 should be systematically considered before any other sequencing approach in subjects with a diagnosis compatible with XLRP.

A classification model relative to splicing for variants of unknown clinical significance: application to the CFTR gene.

Molecular diagnosis of cystic fibrosis and cystic fibrosis transmembrane regulator (CFTR)-related disorders led to the worldwide identification of nearly 1,900 sequence variations in the CFTR gene that consist mainly of private point mutations and small insertions/deletions. Establishing their effect on the function of the encoded protein and therefore their involvement in the disease is still challenging and directly impacts genetic counseling. In this context, we built a decision tree following the international guidelines for the classification of variants of unknown clinical significance (VUCS) in the CFTR gene specifically focused on their consequences on splicing. We applied general and specific criteria, including comprehensive review of literature and databases, familial genetics data, and thorough in silico studies. This model was tested on 15 intronic and exonic VUCS identified in our cohort. Six variants were classified as probably nonpathogenic considering their impact on splicing and eight as probably pathogenic, which include two apparent missense mutations. We assessed the validity of our method by performing minigenes studies and confirmed that 93% (14/15) were correctly classified. We provide in this study a high-performance method that can play a full role in interpreting the results of molecular diagnosis in emergency context, when functional studies are not achievable.

The contribution of GPR98 and DFNB31 genes to a Spanish Usher syndrome type 2 cohort.

BACKGROUND: Usher syndrome type 2 (USH2) is an autosomal recessive disease characterized by moderate to severe hearing loss and retinitis pigmentosa. To date, three disease-causing genes have been identified, USH2A, GPR98, and DFNB31, of which USH2A is clearly the major contributor. The aim of this work was to determine the contribution of GPR98 and DFNB31 genes in a Spanish cohort of USH2A negative patients using exhaustive molecular analysis, including sequencing, dosage, and splicing analysis. METHODS: Linkage analysis was performed to prioritize the gene to study, followed by sequencing of exons and intron-exon boundaries of the selected gene, GPR98 (90 exons) or DFNB31 (12 exons). Functional splicing analyses and comparative genomic hybridization array to detect large rearrangements were performed when appropriate. RESULTS: We confirmed that mutations in GPR98 contribute a significant but minor role to Usher syndrome type 2. In a group of patients referred for molecular diagnosis, 43 had been found to be positive for USH2A mutations, the remaining 19 without USH2A alterations were screened, and seven different mutations were identified in the GPR98 gene in seven patients (five in the homozygous state), of which six were novel. All detected mutations result in a truncated protein; deleterious missense mutations were not found. No pathological mutations were identified in the DFNB31 gene. CONCLUSIONS: In Spain, USH2A and GPR98 are responsible for 95.8% and 5.2% of USH2 mutated cases, respectively. DFNB31 plays a minor role in the Spanish population. There was a group of patients in whom no mutation was found. These findings confirm the importance of including at least GPR98 analysis for comprehensive USH2 molecular diagnosis.

Identification and in vivo functional investigation of a HOMER2 nonstop variant causing hearing loss.

DFNA68 is a rare subtype of autosomal dominant nonsyndromic hearing impairment caused by heterozygous alterations in the HOMER2 gene. To date, only 5 pathogenic or likely pathogenic coding variants, including two missense substitutions (c.188 C > T and c.587 G > C), a single base pair duplication (c.840dupC) and two short deletions (c.592_597delACCACA and c.832_836delCCTCA) have been described in 5 families. In this study, we report a novel HOMER2 variation, identified by massively parallel sequencing, in a Sicilian family suffering from progressive dominant hearing loss over 3 generations. This novel alteration is a nonstop substitution (c.1064 A > G) that converts the translational termination codon (TAG) of the gene into a tryptophan codon (TGG) and is predicted to extend the HOMER2 protein by 10 amino acids. RNA analyses from the proband suggested that HOMER2 transcripts carrying the nonstop variant escaped the non-stop decay pathway. Finally, in vivo studies using a zebrafish animal model and behavioral tests clearly established the deleterious impact of this novel HOMER2 alteration on hearing function. This study identifies the fourth causal variation responsible for DFNA68 and describes a simple in vivo approach to assess the pathogenicity of candidate HOMER2 variants.

Usher syndrome type 2 caused by activation of an USH2A pseudoexon: implications for diagnosis and therapy.

USH2A sequencing in three affected members of a large family, referred for the recessive USH2 syndrome, identified a single pathogenic alteration in one of them and a different mutation in the two affected nieces. As the patients carried a common USH2A haplotype, they likely shared a mutation not found by standard sequencing techniques. Analysis of RNA from nasal cells in one affected individual identified an additional pseudoexon (PE) resulting from a deep intronic mutation. This was confirmed by minigene assay. This is the first example in Usher syndrome (USH) with a mutation causing activation of a PE. The finding of this alteration in eight other individuals of mixed European origin emphasizes the importance of including RNA analysis in a comprehensive diagnostic service. Finally, this mutation, which would not have been found by whole-exome sequencing, could offer, for the first time in USH, the possibility of therapeutic correction by antisense oligonucleotides (AONs).

Non-USH2A mutations in USH2 patients.

We have systematically analyzed the two known minor genes involved in Usher syndrome type 2, DFNB31 and GPR98, for mutations in a cohort of 31 patients not linked to USH2A. PDZD7, an Usher syndrome type 2 (USH2) related gene, was analyzed when indicated. We found that mutations in GPR98 contribute significantly to USH2. We report 17 mutations in 10 individuals, doubling the number of GPR98 mutations reported to date. In contrast to mutations in usherin, the mutational spectrum of GPR98 predominantly results in a truncated protein product. This is true even when the mutation affects splicing, and we have incorporated a splicing reporter minigene assay to show this, where appropriate. Only two mutations were found which we believe to be genuine missense changes. Discrepancy in the mutational spectrum between GPR98 and USH2A is discussed. Only two patients were found with mutations in DFNB31, showing that mutations of this gene contribute to only a very small extent to USH2. Close examination of the clinical details, where available, for patients in whom no mutation was found in USH2A, GPR98, or DFNB31, showed that most of them had atypical features. In effect, these three genes account for the vast majority of USH2 patients and their analysis provide a robust pathway for routine molecular diagnosis.

Identification of the First Single GSDME Exon 8 Structural Variants Associated with Autosomal Dominant Hearing Loss.

GSDME, also known as DFNA5, is a gene implicated in autosomal dominant nonsyndromic hearing loss (ADNSHL), affecting, at first, the high frequencies with a subsequent progression over all frequencies. To date, all the GSDME pathogenic variants associated with deafness lead to skipping of exon 8. In two families with apparent ADNSHL, massively parallel sequencing (MPS) integrating a coverage-based method for detection of copy number variations (CNVs) was applied, and it identified the first two causal GSDME structural variants affecting exon 8. The deleterious impact of the c.991-60_1095del variant, which includes the acceptor splice site sequence of exon 8, was confirmed by the study of the proband's transcripts. The second mutational event is a complex rearrangement that deletes almost all of the exon 8 sequence. This study increases the mutational spectrum of the GSDME gene and highlights the crucial importance of MPS data for the detection of GSDME exon 8 deletions, even though the identification of a causal single-exon CNV by MPS analysis is still challenging.

Reclassification of a TMC1 synonymous substitution as a variant disrupting splicing regulatory elements associated with recessive hearing loss.

Alterations of the transmembrane channel-like 1 gene (TMC1) are involved in autosomal recessive and dominant nonsyndromic hearing loss (NSHL). To date, up to 117 causal variants including substitutions, insertions and splice variants have been reported in families from different populations. In a patient suffering from severe prelingual NSHL, we identified, in the homozygous state, the previously considered likely benign synonymous c.627C>T; p.(Leu209=) substitution. We used in silico tools predicting variant-induced alterations of splicing regulatory elements (SREs) and pinpointed this transition as a candidate splice-altering variation. Functional splicing analysis, using a minigene assay, confirmed that the variant altered a critical regulatory sequence which is essential for the exon 11 inclusion in the TMC1 transcripts. This result was reinforced by the analysis of orthologous TMC1 mammalian sequences for which the deleterious effect on the mRNA processing of a native thymidine was always counteracted by the presence of a stronger donor splice site or additional enhancer motifs. This study demonstrates, for the first time, the pathogenicity of the c.627C>T alteration leading to its reclassification as a causal variant impacting SREs and highlights the major importance of exhaustive studies to accurately evaluate the pathogenicity of a variant, regardless of the variation type.