Functional analysis of germline VANGL2 variants using rescue assays of vangl2 knockout zebrafish.

Developmental studies have shown that the evolutionarily conserved Wnt Planar Cell Polarity (PCP) pathway is essential for the development of a diverse range of tissues and organs including the brain, spinal cord, heart and sensory organs, as well as establishment of the left-right body axis. Germline mutations in the highly conserved PCP gene VANGL2 in humans have only been associated with central nervous system malformations, and functional testing to understand variant impact has not been performed. Here we report three new families with missense variants in VANGL2 associated with heterotaxy and congenital heart disease p.(Arg169His), non-syndromic hearing loss p.(Glu465Ala) and congenital heart disease with brain defects p.(Arg135Trp). To test the in vivo impact of these and previously described variants, we have established clinically-relevant assays using mRNA rescue of the vangl2 mutant zebrafish. We show that all variants disrupt Vangl2 function, although to different extents and depending on the developmental process. We also begin to identify that different VANGL2 missense variants may be haploinsufficient and discuss evidence in support of pathogenicity. Together, this study demonstrates that zebrafish present a suitable pipeline to investigate variants of unknown significance and suggests new avenues for investigation of the different developmental contexts of VANGL2 function that are clinically meaningful.

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 ).

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.

CHM mutation spectrum and disease: An update at the time of human therapeutic trials.

Choroideremia is an X-linked inherited retinal disorder (IRD) characterized by the degeneration of retinal pigment epithelium, photoreceptors, choriocapillaris and choroid affecting males with variable phenotypes in female carriers. Unlike other IRD, characterized by a large clinical and genetic heterogeneity, choroideremia shows a specific phenotype with causative mutations in only one gene, CHM. Ongoing gene replacement trials raise further interests in this disorder. We describe here the clinical and genetic data from a French cohort of 45 families, 25 of which carry novel variants, in the context of 822 previously reported choroideremia families. Most of the variants represent loss-of-function mutations with eleven families having large (i.e. ≥6 kb) genomic deletions, 18 small insertions, deletions or insertion deletions, six showing nonsense variants, eight splice site variants and two missense variants likely to affect splicing. Similarly, 822 previously published families carry mostly loss-of-function variants. Recurrent variants are observed worldwide, some of which linked to a common ancestor, others arisen independently in specific CHM regions prone to mutations. Since all exons of CHM may harbor variants, Sanger sequencing combined with quantitative polymerase chain reaction or multiplex ligation-dependent probe amplification experiments are efficient to achieve the molecular diagnosis in patients with typical choroideremia features.

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.

CRB1-Related Retinal Dystrophies in a Cohort of 50 Patients: A Reappraisal in the Light of Specific Müller Cell and Photoreceptor CRB1 Isoforms.

Pathogenic variants in CRB1 lead to diverse recessive retinal disorders from severe Leber congenital amaurosis to isolated macular dystrophy. Until recently, no clear phenotype-genotype correlation and no appropriate mouse models existed. Herein, we reappraise the phenotype-genotype correlation of 50 patients with regards to the recently identified CRB1 isoforms: a canonical long isoform A localized in Müller cells (12 exons) and a short isoform B predominant in photoreceptors (7 exons). Twenty-eight patients with early onset retinal dystrophy (EORD) consistently had a severe Müller impairment, with variable impact on the photoreceptors, regardless of isoform B expression. Among them, two patients expressing wild type isoform B carried one variant in exon 12, which specifically damaged intracellular protein interactions in Müller cells. Thirteen retinitis pigmentosa patients had mainly missense variants in laminin G-like domains and expressed at least 50% of isoform A. Eight patients with the c.498_506del variant had macular dystrophy. In one family homozygous for the c.1562C>T variant, the brother had EORD and the sister macular dystrophy. In contrast with the mouse model, these data highlight the key role of Müller cells in the severity of CRB1-related dystrophies in humans, which should be taken into consideration for future clinical trials.

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.

PHENOTYPIC CHARACTERISTICS OF ROD-CONE DYSTROPHY ASSOCIATED WITH MYO7A MUTATIONS IN A LARGE FRENCH COHORT.

PURPOSE: To document the rod-cone dystrophy phenotype of patients with Usher syndrome type 1 (USH1) harboring MYO7A mutations. METHODS: Retrospective cohort study of 53 patients (42 families) with biallelic MYO7A mutations who underwent comprehensive examination, including functional visual tests and multimodal retinal imaging. Genetic analysis was performed either using a multiplex amplicon panel or through direct sequencing. Data were analyzed with IBM SPSS Statistics software v. 21.0. RESULTS: Fifty different genetic variations including 4 novel were identified. Most patients showed a typical rod-cone dystrophy phenotype, with best-corrected visual acuity and central visual field deteriorating linearly with age. At age 29, binocular visual field demonstrated an average preservation of 50 central degrees, constricting by 50% within 5 years. Structural changes based on spectral domain optical coherence tomography, short wavelength autofluorescence, and near-infrared autofluorescence measurements did not however correlate with age. Our study revealed a higher percentage of epiretinal membranes and cystoid macular edema in patients with MYO7A mutations compared with rod-cone dystrophy patients with other mutations. Subgroup analyses did not reveal substantial genotype-phenotype correlations. CONCLUSION: To the best of our knowledge, this is the largest French cohort of patients with MYO7A mutations reported to date. Functional visual characteristics of this subset of patients followed a linear decline as in other typical rod-cone dystrophy, but structural changes were variable indicating the need for a case-by-case evaluation for prognostic prediction and choice of potential therapies.

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.

Mutations in PLS1, encoding fimbrin, cause autosomal dominant nonsyndromic hearing loss.

Nonsyndromic hearing loss (NSHL), a common sensory disorder, is characterized by high clinical and genetic heterogeneity (i.e., approximately 115 genes and 170 loci so far identified). Nevertheless, almost half of patients submitted for genetic testing fail to receive a conclusive molecular diagnosis. We used next-generation sequencing to identify causal variants in PLS1 (c.805G>A, p.[E269K]; c.713G>T, p.[L238R], and c.383T>C, p.[F128S]) in three unrelated families of European ancestry with autosomal dominant NSHL. PLS1 encodes Plastin 1 (also called fimbrin), one of the most abundant actin-bundling proteins of the stereocilia. In silico protein modeling suggests that all variants destabilize the structure of the actin-binding domain 1, likely reducing the protein's ability to bind F actin. The role of PLS1 gene in hearing function is further supported by the recent demonstration that Pls1-/- mice show a hearing loss phenotype similar to that of our patients. In summary, we report PLS1 as a novel gene for autosomal dominant NSHL, suggesting that this gene is required for normal hearing in humans and mice.

Pathogenicity of a novel missense variant associated with choroideremia and its impact on gene replacement therapy.

Choroideremia (CHM) is an inherited retinal dystrophy characterised by progressive degeneration of photoreceptors, retinal pigment epithelium (RPE) and underlying choroid. It is caused by loss-of-function mutations in CHM, which has an X-linked inheritance, and is thus an ideal candidate for gene replacement strategies. CHM encodes REP1, which plays a key role in the prenylation of Rab GTPases. We recently showed that an induced pluripotent stem cell (iPSc)-derived RPE model for CHM is fully functional and reproduces the underlying prenylation defect. This criterion can thus be used for testing the pathogenic nature of novel variants. Until recently, missense variants were not associated with CHM. Currently, at least nine such variants have been reported but only two have been shown to be pathogenic. We report here the characterisation of the third pathogenic missense CHM variant, p.Leu457Pro. Clinically, the associated phenotype is indistinguishable from that of loss-of-function mutations. By contrast, this missense variant results in wild type CHM expression levels and detectable levels of mutant protein. The prenylation status of patient-specific fibroblasts and iPSc-derived RPE is within the range observed for loss-of-function mutations, consistent with the clinical phenotype. Lastly, considering the current climate of CHM gene therapy, we assayed whether the presence of mutant REP1 could interfere with a gene replacement strategy by testing the prenylation status of patient-specific iPSc-derived RPE following AAV-mediated gene transfer. Our results show that correction of the functional defect is possible and highlight the predictive value of these models for therapy screening prior to inclusion in clinical trials.

Mutation of SLC9A1, encoding the major Na⁺/H⁺ exchanger, causes ataxia-deafness Lichtenstein-Knorr syndrome.

Lichtenstein-Knorr syndrome is an autosomal recessive condition that associates sensorineural hearing loss and cerebellar ataxia. Here, we report the first identification of a gene involved in Lichtenstein-Knorr syndrome. By using a combination of homozygosity mapping and whole-exome sequencing, we identified the homozygous p.Gly305Arg missense mutation in SLC9A1 that segregates with the disease in a large consanguineous family. Mutant glycine 305 is a highly conserved amino acid present in the eighth transmembrane segment of all metazoan orthologues of NHE1, the Na(+)/H(+) exchanger 1, encoded by SLC9A1. We demonstrate that the p.Gly305Arg mutation causes the near complete de-glycosylation, mis-targeting and loss of proton pumping activity of NHE1. The comparison of our family with the phenotypes of spontaneous and knockout Slc9a1 murine models demonstrates that the association between ataxia and hearing loss is caused by complete or near complete loss of function of NHE1 and altered regulation of pHi in the central nervous system.

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.

HGVS Recommendations for the Description of Sequence Variants: 2016 Update.

The consistent and unambiguous description of sequence variants is essential to report and exchange information on the analysis of a genome. In particular, DNA diagnostics critically depends on accurate and standardized description and sharing of the variants detected. The sequence variant nomenclature system proposed in 2000 by the Human Genome Variation Society has been widely adopted and has developed into an internationally accepted standard. The recommendations are currently commissioned through a Sequence Variant Description Working Group (SVD-WG) operating under the auspices of three international organizations: the Human Genome Variation Society (HGVS), the Human Variome Project (HVP), and the Human Genome Organization (HUGO). Requests for modifications and extensions go through the SVD-WG following a standard procedure including a community consultation step. Version numbers are assigned to the nomenclature system to allow users to specify the version used in their variant descriptions. Here, we present the current recommendations, HGVS version 15.11, and briefly summarize the changes that were made since the 2000 publication. Most focus has been on removing inconsistencies and tightening definitions allowing automatic data processing. An extensive version of the recommendations is available online, at http://www.HGVS.org/varnomen.

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.

Novel deletions involving the USH2A gene in patients with Usher syndrome and retinitis pigmentosa.

PURPOSE: The aim of the present work was to identify and characterize large rearrangements involving the USH2A gene in patients with Usher syndrome and nonsyndromic retinitis pigmentosa. METHODS: The multiplex ligation-dependent probe amplification (MLPA) technique combined with a customized array-based comparative genomic hybridization (aCGH) analysis was applied to 40 unrelated patients previously screened for point mutations in the USH2A gene in which none or only one pathologic mutation was identified. RESULTS: We detected six large deletions involving USH2A in six out of the 40 cases studied. Three of the patients were homozygous for the deletion, and the remaining three were compound heterozygous with a previously identified USH2A point mutation. In five of these cases, the patients displayed Usher type 2, and the remaining case displayed nonsyndromic retinitis pigmentosa. The exact breakpoint junctions of the deletions found in USH2A in four of these cases were characterized. CONCLUSIONS: Our study highlights the need to develop improved efficient strategies of mutation screening based upon next generation sequencing (NGS) that reduce cost, time, and complexity and allow simultaneous identification of all types of disease-causing mutations in diagnostic procedures.

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.