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.

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.

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.

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.

When Familial Hearing Loss Means Genetic Heterogeneity: A Model Case Report.

We describe a family with both hearing loss (HL) and thrombocytopenia, caused by pathogenic variants in three genes. The proband was a child with neonatal thrombocytopenia, childhood-onset HL, hyper-laxity and severe myopia. The child's mother (and some of her relatives) presented with moderate thrombocytopenia and adulthood-onset HL. The child's father (and some of his relatives) presented with adult-onset HL. An HL panel analysis, completed by whole exome sequencing, was performed in this complex family. We identified three pathogenic variants in three different genes: MYH9, MYO7A and ACTG1. The thrombocytopenia in the child and her mother is explained by the MYH9 variant. The post-lingual HL in the paternal branch is explained by the MYO7A variant, absent in the proband, while the congenital HL of the child is explained by a de novo ACTG1 variant. This family, in which HL segregates, illustrates that multiple genetic conditions coexist in individuals and make patient care more complex than expected.

A 4.6 Mb Inversion Leading to PCDH15-LINC00844 and BICC1-PCDH15 Fusion Transcripts as a New Pathogenic Mechanism Implicated in Usher Syndrome Type 1.

Usher type 1 syndrome is a rare autosomal recessive disorder involving congenital severe-to-profound hearing loss, development of vision impairment in the first decade, and severe balance difficulties. The PCDH15 gene, one of the five genes implicated in this disease, is involved in 8-20% of cases. In this study, we aimed to identify and characterize the two causal variants in a French patient with typical Usher syndrome clinical features. Massively parallel sequencing-based gene panel and screening for large rearrangements were used, which detected a single multi-exon deletion in the PCDH15 gene. As the second pathogenic event was likely localized in the unscreened regions of the gene, PCDH15 transcripts from cultured nasal cells were analyzed and revealed a loss of junction between exon 13 and exon 14. This aberration could be explained by the identification of two fusion transcripts, PCDH15-LINC00844 and BICC1-PCDH15, originating from a 4.6 Mb inversion. This complex chromosomal rearrangement could not be detected by our diagnostic approach but was instead characterized by long-read sequencing, which offers the possibility of detecting balanced structural variants (SVs). This finding extends our knowledge of the mutational spectrum of the PCDH15 gene with the first ever identification of a large causal paracentric inversion of chromosome 10 and illustrates the utility of screening balanced SVs in an exhaustive molecular diagnostic approach.

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.

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.

Assessment of the latest NGS enrichment capture methods in clinical context.

Enrichment capture methods for NGS are widely used, however, they evolve rapidly and it is necessary to periodically measure their strengths and weaknesses before transfer to diagnostic services. We assessed two recently released custom DNA solution-capture enrichment methods for NGS, namely Illumina NRCCE and Agilent SureSelect(QXT), against a reference method NimbleGen SeqCap EZ Choice on a similar gene panel, sharing 678 kb and 110 genes. Two Illumina MiSeq runs of 12 samples each have been performed, for each of the three methods, using the same 24 patients (affected with sensorineural disorders). Technical outcomes have been computed and compared, including depth and evenness of coverage, enrichment in targeted regions, performance in GC-rich regions and ability to generate consistent variant datasets. While we show that the three methods resulted in suitable datasets for standard DNA variant discovery, we describe significant differences between the results for the above parameters. NimbleGen offered the best depth of coverage and evenness, while NRCCE showed the highest on target levels but high duplicate rates. SureSelect(QXT) showed an overall quality close to that of NimbleGen. The new methods exhibit reduced preparation time but behave differently. These findings will guide laboratories in their choice of library enrichment approach.

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.

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.

Experience of targeted Usher exome sequencing as a clinical test.

We show that massively parallel targeted sequencing of 19 genes provides a new and reliable strategy for molecular diagnosis of Usher syndrome (USH) and nonsyndromic deafness, particularly appropriate for these disorders characterized by a high clinical and genetic heterogeneity and a complex structure of several of the genes involved. A series of 71 patients including Usher patients previously screened by Sanger sequencing plus newly referred patients was studied. Ninety-eight percent of the variants previously identified by Sanger sequencing were found by next-generation sequencing (NGS). NGS proved to be efficient as it offers analysis of all relevant genes which is laborious to reach with Sanger sequencing. Among the 13 newly referred Usher patients, both mutations in the same gene were identified in 77% of cases (10 patients) and one candidate pathogenic variant in two additional patients. This work can be considered as pilot for implementing NGS for genetically heterogeneous diseases in clinical service.

Four-year follow-up of diagnostic service in USH1 patients.

PURPOSE: The purpose of this study was to establish the mutation spectrum of an Usher type I cohort of 61 patients from France and to describe a diagnostic strategy, including a strategy for estimating the pathogenicity of sequence changes. METHODS: To optimize the identification of Usher (USH)-causative mutations, taking into account the genetic heterogeneity, preliminary haplotyping at the five USH1 loci was performed to prioritize the gene to be sequenced, as previously described. Coding exons and flanking intronic sequences were sequenced and, where necessary, semiquantitative PCR and multiplex ligation-dependent probe amplification (MLPA) were performed to detect large genomic rearrangements. RESULTS: Four years ' experience confirms that the chosen approach provides an efficient diagnostic service. Sixty-one patients showed an abnormal genotype in one of the five USH1 genes. Genetic heterogeneity was confirmed, and, although MYO7A remains the major gene, involvement of other genes is considerable. Distribution of missense, splicing, premature termination codons (PTCs; due to point substitution and small deletions/ or insertions), and large genomic alterations was determined among the USH genes and clearly highlights the need to pay special attention to the diagnostic approach and interpretation, depending on the mutated gene. CONCLUSIONS: Over the 4 years of a diagnostic service offering USH1 patient testing, pathogenic genotypes were identified in most cases (>90%). The complexity and heterogeneity of mutations reinforces the need for a comprehensive approach. Because 32% of the mutations are newly described, the results show that a screening strategy based on known mutations would have solved less than 55% of the cases.

Nasal epithelial cells are a reliable source to study splicing variants in Usher syndrome.

We have shown that nasal ciliated epithelium, which can be easily biopsied under local anesthetic, provides a good source of RNA transcripts from eight of the nine known genes that cause Usher syndrome, namely, MYO7A, USH1C, CDH23, PCDH15, USH1G for Usher type 1, and USH2A, GPR98, WHRN for Usher type 2. Furthermore, the known or predicted effect on mRNA splicing of eight variants was faithfully reproduced in the biopsied sample as measured by nested RT-PCR. These included changes at the canonical acceptor site, changes within the noncanonical acceptor site and both synonymous and nonsynonymous amino acid changes. This shows that mRNA analysis by this method will help in assessing the pathogenic effect of variants, which is a major problem in the molecular diagnosis of Usher syndrome.

Ex vivo splicing assays of mutations at noncanonical positions of splice sites in USHER genes.

Molecular diagnosis in Usher syndrome type 1 and 2 patients led to the identification of 21 sequence variations located in noncanonical positions of splice sites in MYO7A, CDH23, USH1C, and USH2A genes. To establish experimentally the splicing pattern of these substitutions, whose impact on splicing is not always predictable by available softwares, ex vivo splicing assays were performed. The branch-point mapping strategy was also used to investigate further a putative branch-point mutation in USH2A intron 43. Aberrant splicing was demonstrated for 16 of the 21 (76.2%) tested sequence variations. The mutations resulted more frequently in activation of a nearby cryptic splice site or use of a de novo splice site than exon skipping (37.5%). This study allowed the reclassification as splicing mutations of one silent (c.7872G>A (p.Glu2624Glu) in CDH23) and four missense mutations (c.2993G>A (p.Arg998Lys) in USH2A, c.592G>A (p.Ala198Thr), c.3503G>C [p.Arg1168Pro], c.5944G>A (p.Gly1982Arg) in MYO7A), whereas it provided clues about a role in structure/function in four other cases: c.802G>A (p.Gly268Arg), c.653T>A (p.Val218Glu) (USH2A), and c.397C>T (p.His133Tyr), c.3502C>T (p.Arg1168Trp) (MYO7A). Our data provide insights into the contribution of splicing mutations in Usher genes and illustrate the need to define accurately their splicing outcome for diagnostic purposes.

Molecular screening of deafness in Algeria: high genetic heterogeneity involving DFNB1 and the Usher loci, DFNB2/USH1B, DFNB12/USH1D and DFNB23/USH1F.

A systematic approach, involving haplotyping and genotyping, to the molecular diagnosis of non-syndromic deafness within 50 families and 9 sporadic cases from Algeria is described. Mutations at the DFNB1 locus (encompassing the GJB2 and GJB6 genes) are responsible for more than half of autosomal recessive prelingual non-syndromic deafness in various populations. A c.35delG mutation can account for up to 85% of GJB2 mutations and two large deletions del(GJB6-D13S1830) and del(GJB6-D13S1854) have also been reported in several population groups. In view of the genetic heterogeneity a strategy was developed which involved direct analysis of DFNB1. In negative familial cases, haplotype analysis was carried out, where possible, to exclude DFNB1 mutations. Following this, haplotype analysis of five Usher syndrome loci, sometimes involved in autosomal non-syndromic hearing loss, was carried out to identify cases in which Usher gene sequencing was indicated. When homozygosity was observed at a locus in a consanguineous family, the corresponding gene was exhaustively sequenced. Pathogenic DFNB1 genotypes were identified in 40% of the cases. Of the 21 cases identified with 2 pathogenic mutations, c.35delG represented 76% of the mutated alleles. The additional mutations were one nonsense, two missense and one splicing mutation. Four additional patients were identified with a single DFNB1 mutation. None carried the large deletions. Three families with non-syndromic deafness carried novel unclassified variants (UVs) in MYO7A (1 family) and CDH23 (2 families) of unknown pathogenic effect. Additionally, molecular diagnosis was carried out on two Usher type I families and pathogenic mutations in MYO7A and PCDH15 were found.

UMD-USHbases: a comprehensive set of databases to record and analyse pathogenic mutations and unclassified variants in seven Usher syndrome causing genes.

Using the Universal Mutation Database (UMD) software, we have constructed "UMD-USHbases", a set of relational databases of nucleotide variations for seven genes involved in Usher syndrome (MYO7A, CDH23, PCDH15, USH1C, USH1G, USH3A and USH2A). Mutations in the Usher syndrome type I causing genes are also recorded in non-syndromic hearing loss cases and mutations in USH2A in non-syndromic retinitis pigmentosa. Usher syndrome provides a particular challenge for molecular diagnostics because of the clinical and molecular heterogeneity. As many mutations are missense changes, and all the genes also contain apparently non-pathogenic polymorphisms, well-curated databases are crucial for accurate interpretation of pathogenicity. Tools are provided to assess the pathogenicity of mutations, including conservation of amino acids and analysis of splice-sites. Reference amino acid alignments are provided. Apparently non-pathogenic variants in patients with Usher syndrome, at both the nucleotide and amino acid level, are included. The UMD-USHbases currently contain more than 2,830 entries including disease causing mutations, unclassified variants or non-pathogenic polymorphisms identified in over 938 patients. In addition to data collected from 89 publications, 15 novel mutations identified in our laboratory are recorded in MYO7A (6), CDH23 (8), or PCDH15 (1) genes. Information is given on the relative involvement of the seven genes, the number and distribution of variants in each gene. UMD-USHbases give access to a software package that provides specific routines and optimized multicriteria research and sorting tools. These databases should assist clinicians and geneticists seeking information about mutations responsible for Usher syndrome.

Large genomic rearrangements within the PCDH15 gene are a significant cause of USH1F syndrome.

PURPOSE: Protocadherin-15 (PCDH15) is one of the five genes currently identified as being mutated in Usher 1 syndrome and defines Usher syndrome type 1F (USH1F). When PCDH15 was systematically analyzed for mutations in a cohort of USH1 patients, a number of deletions were found. Here we characterize these deletions as to extent, position, and breakpoints. METHODS: Microsatellite and single nucleotide polymorphism (SNP) analyses, used in a preliminary survey of an Usher cohort of 31 patients, revealed large deletions in three patients. These deletions were further characterized by semiquantitative PCR assays to narrow down the breakpoints. RESULTS: The analysis of the three large deletions revealed that all six breakpoints are different. The breakpoint junction was identified in one patient and the four other breakpoints were mapped to 4 kb. There were no specific distinguishing features of the isolated breakpoints. CONCLUSIONS: A complete screen of PCDH15 should include a search for large deletions. Failure to screen for gross genomic rearrangements is likely to significantly lower the mutation detection rate. A likely explanation for the high rate of such deletions is the unusual gene structure. PCDH15 gene spans nearly 1 Mb for a corresponding open reading frame (ORF) of 7,021 bp. The intron sizes of PCDH15 are up to 150 kb, and the first three exons of the gene cover 0.42 Mb. The genomic structure of any gene should be taken into consideration when designing a mutation screening strategy.

Molecular epidemiology of DFNB1 deafness in France.

BACKGROUND: Mutations in the GJB2 gene have been established as a major cause of inherited non syndromic deafness in different populations. A high number of sequence variations have been described in the GJB2 gene and the associated pathogenic effects are not always clearly established. The prevalence of a number of mutations is known to be population specific, and therefore population specific testing should be a prerequisite step when molecular diagnosis is offered. Moreover, population studies are needed to determine the contribution of GJB2 variants to deafness. We present our findings from the molecular diagnostic screening of the GJB2 and GJB6 genes over a three year period, together with a population-based study of GJB2 variants. METHODS AND RESULTS: Molecular studies were performed using denaturing High Performance Liquid Chromatograghy (DHPLC) and sequencing of the GJB2 gene. Over the last 3 years we have studied 159 families presenting sensorineural hearing loss, including 84 with non syndromic, stable, bilateral deafness. Thirty families were genotyped with causative mutations. In parallel, we have performed a molecular epidemiology study on more than 3000 dried blood spots and established the frequency of the GJB2 variants in our population. Finally, we have compared the prevalence of the variants in the hearing impaired population with the general population. CONCLUSION: Although a high heterogeneity of sequence variation was observed in patients and controls, the 35delG mutation remains the most common pathogenic mutation in our population. Genetic counseling is dependent on the knowledge of the pathogenicity of the mutations and remains difficult in a number of cases. By comparing the sequence variations observed in hearing impaired patients with those sequence variants observed in general population, from the same ethnic background, we show that the M34T, V37I and R127H variants can not be responsible for profound or severe deafness.