Diagnostic genome sequencing improves diagnostic yield: a prospective single-centre study in 1000 patients with inherited eye diseases.

PURPOSE: Genome sequencing (GS) is expected to reduce the diagnostic gap in rare disease genetics. We aimed to evaluate a scalable framework for genome-based analyses 'beyond the exome' in regular care of patients with inherited retinal degeneration (IRD) or inherited optic neuropathy (ION). METHODS: PCR-free short-read GS was performed on 1000 consecutive probands with IRD/ION in routine diagnostics. Complementary whole-blood RNA-sequencing (RNA-seq) was done in a subset of 74 patients. An open-source bioinformatics analysis pipeline was optimised for structural variant (SV) calling and combined RNA/DNA variation interpretation. RESULTS: A definite genetic diagnosis was established in 57.4% of cases. For another 16.7%, variants of uncertain significance were identified in known IRD/ION genes, while the underlying genetic cause remained unresolved in 25.9%. SVs or alterations in non-coding genomic regions made up for 12.7% of the observed variants. The RNA-seq studies supported the classification of two unclear variants. CONCLUSION: GS is feasible in clinical practice and reliably identifies causal variants in a substantial proportion of individuals. GS extends the diagnostic yield to rare non-coding variants and enables precise determination of SVs. The added diagnostic value of RNA-seq is limited by low expression levels of the major IRD disease genes in blood.

Mutations at a split codon in the GTPase-encoding domain of OPA1 cause dominant optic atrophy through different molecular mechanisms.

Exonic (i.e. coding) variants in genes associated with disease can exert pathogenic effects both at the protein and mRNA level, either by altering the amino acid sequence or by affecting pre-mRNA splicing. The latter is often neglected due to the lack of RNA analyses in genetic diagnostic testing. In this study we considered both pathomechanisms and performed a comprehensive analysis of nine exonic nucleotide changes in OPA1, which is the major gene underlying autosomal dominant optic atrophy (DOA) and is characterized by pronounced allelic heterogeneity. We focused on the GTPase-encoding domain of OPA1, which harbors most of the missense variants associated with DOA. Given that the consensus splice sites extend into the exons, we chose a split codon, namely codon 438, for our analyses. Variants at this codon are the second most common cause of disease in our large cohort of DOA patients harboring disease-causing variants in OPA1. In silico splice predictions, heterologous splice assays, analysis of patient's RNA when available, and protein modeling revealed different molecular outcomes for variants at codon 438. The wildtype aspartate residue at amino acid position 438 is directly involved in the dimerization of OPA1 monomers. We found that six amino acid substitutions at codon 438 (i.e. all substitutions of the first and second nucleotide of the codon) destabilized dimerization while only substitutions of the first nucleotide of the codon caused exon skipping. Our study highlights the value of combining RNA analysis and protein modeling approaches to accurately assign patients to future precision therapies.

A duplication on chromosome 16q12 affecting the IRXB gene cluster is associated with autosomal dominant cone dystrophy with early tritanopic color vision defect.

Cone dystrophies are a rare subgroup of inherited retinal dystrophies and hallmarked by color vision defects, low or decreasing visual acuity and central vision loss, nystagmus and photophobia. Applying genome-wide linkage analysis and array comparative genome hybridization, we identified a locus for autosomal dominant cone dystrophy on chromosome 16q12 in four independent multigeneration families. The locus is defined by duplications of variable size with a smallest region of overlap of 608 kb affecting the IRXB gene cluster and encompasses the genes IRX5 and IRX6. IRX5 and IRX6 belong to the Iroquois (Iro) protein family of homeodomain-containing transcription factors involved in patterning and regionalization of embryonic tissue in vertebrates, including the eye and the retina. All patients presented with a unique progressive cone dystrophy phenotype hallmarked by early tritanopic color vision defects. We propose that the disease underlies a misregulation of the IRXB gene cluster on chromosome 16q12 and demonstrate that overexpression of Irx5a and Irx6a, the two orthologous genes in zebrafish, results in visual impairment in 5-day-old zebrafish larvae.

Optical genome mapping and revisiting short-read genome sequencing data reveal previously overlooked structural variants disrupting retinal disease-associated genes.

PURPOSE: Structural variants (SVs) play an important role in inherited retinal diseases (IRD). Although the identification of SVs significantly improved upon the availability of genome sequencing, it is expected that involvement of SVs in IRDs is higher than anticipated. We revisited short-read genome sequencing data to enhance the identification of gene-disruptive SVs. METHODS: Optical genome mapping was performed to improve SV detection in short-read genome sequencing-negative cases. In addition, reanalysis of short-read genome sequencing data was performed to improve the interpretation of SVs and to re-establish SV prioritization criteria. RESULTS: In a monoallelic USH2A case, optical genome mapping identified a pericentric inversion (173 megabase), with 1 breakpoint disrupting USH2A. Retrospectively, the variant could be observed in genome sequencing data but was previously deemed false positive. Reanalysis of short-read genome sequencing data (427 IRD cases) was performed which yielded 30 pathogenic SVs affecting, among other genes, USH2A (n = 15), PRPF31 (n = 3), and EYS (n = 2). Eight of these (>25%) were overlooked during previous analyses. CONCLUSION: Critical evaluation of our findings allowed us to re-establish and improve our SV prioritization and interpretation guidelines, which will prevent missing pathogenic events in future analyses. Our data suggest that more attention should be paid to SV interpretation and the current contribution of SVs in IRDs is still underestimated.

DNAJC30 disease-causing gene variants in a large Central European cohort of patients with suspected Leber's hereditary optic neuropathy and optic atrophy.

BACKGROUND: Leber's hereditary optic neuropathy (LHON) has been considered a prototypical mitochondriopathy and a textbook example for maternal inheritance linked to certain disease-causing variants in the mitochondrial genome. Recently, an autosomal recessive form of LHON (arLHON) has been described, caused by disease-causing variants in the nuclear encoded gene DNAJC30. METHODS AND RESULTS: In this study, we screened the DNAJC30 gene in a large Central European cohort of patients with a clinical diagnosis of LHON or other autosomal inherited optic atrophies (OA). We identified likely pathogenic variants in 35/1202 patients, corresponding to a detection rate of 2.9%. The previously described missense variant c.152A>G;p.(Tyr51Cys) accounts for 90% of disease-associated alleles in our cohort and we confirmed a strong founder effect. Furthermore, we identified two novel pathogenic variants in DNAJC30: the nonsense variant c.610G>T;p.(Glu204*) and the in-frame deletion c.230_232del;p.(His77del). Clinical investigation of the patients with arLHON revealed a younger age of onset, a more frequent bilateral onset and an increased clinically relevant recovery compared with LHON associated with disease-causing variants in the mitochondrial DNA. CONCLUSION: This study expands previous findings on arLHON and emphasises the importance of DNAJC30 in the genetic diagnostics of LHON and OA in European patients.

Genetic and Clinical Profile of Retinopathies Due to Disease-Causing Variants in Leber Congenital Amaurosis (LCA)-Associated Genes in a Large German Cohort.

To report the spectrum of Leber congenital amaurosis (LCA) associated genes in a large German cohort and to delineate their associated phenotype. Local databases were screened for patients with a clinical diagnosis of LCA and for patients with disease-causing variants in known LCA-associated genes independent of their clinical diagnosis. Patients with a mere clinical diagnosis were invited for genetic testing. Genomic DNA was either analyzed in a diagnostic-genetic or research setup using various capture panels for syndromic and non-syndromic IRD (inherited retinal dystrophy) genes. Clinical data was obtained mainly retrospectively. Patients with genetic and phenotypic information were eventually included. Descriptive statistical data analysis was performed. A total of 105 patients (53 female, 52 male, age 3-76 years at the time of data collection) with disease-causing variants in 16 LCA-associated genes were included. The genetic spectrum displayed variants in the following genes: CEP290 (21%), CRB1 (21%), RPE65 (14%), RDH12 (13%), AIPL1 (6%), TULP1 (6%), and IQCB1 (5%), and few cases harbored pathogenic variants in LRAT, CABP4, NMNAT1, RPGRIP1, SPATA7, CRX, IFT140, LCA5, and RD3 (altogether accounting for 14%). The most common clinical diagnosis was LCA (53%, 56/105) followed by retinitis pigmentosa (RP, 40%, 42/105), but also other IRDs were seen (cone-rod dystrophy, 5%; congenital stationary night blindness, 2%). Among LCA patients, 50% were caused by variants in CEP290 (29%) and RPE65 (21%), whereas variants in other genes were much less frequent (CRB1 11%, AIPL1 11%, IQCB1 9%, and RDH12 7%, and sporadically LRAT, NMNAT1, CRX, RD3, and RPGRIP1). In general, the patients showed a severe phenotype hallmarked by severely reduced visual acuity, concentric narrowing of the visual field, and extinguished electroretinograms. However, there were also exceptional cases with best corrected visual acuity as high as 0.8 (Snellen), well-preserved visual fields, and preserved photoreceptors in spectral domain optical coherence tomography. Phenotypic variability was seen between and within genetic subgroups. The study we are presenting pertains to a considerable LCA group, furnishing valuable comprehension of the genetic and phenotypic spectrum. This knowledge holds significance for impending gene therapeutic trials. In this German cohort, CEP290 and CRB1 are the most frequently mutated genes. However, LCA is genetically highly heterogeneous and exhibits clinical variability, showing overlap with other IRDs. For any therapeutic gene intervention, the disease-causing genotype is the primary criterion for treatment access, but the clinical diagnosis, state of the retina, number of to be treated target cells, and the time point of treatment will be crucial.

Biallelic Variants in TULP1 Are Associated with Heterogeneous Phenotypes of Retinal Dystrophy.

Biallelic pathogenic variants in TULP1 are mostly associated with severe rod-driven inherited retinal degeneration. In this study, we analyzed clinical heterogeneity in 17 patients and characterized the underlying biallelic variants in TULP1. All patients underwent thorough ophthalmological examinations. Minigene assays and structural analyses were performed to assess the consequences of splice variants and missense variants. Three patients were diagnosed with Leber congenital amaurosis, nine with early onset retinitis pigmentosa, two with retinitis pigmentosa with an onset in adulthood, one with cone dystrophy, and two with cone-rod dystrophy. Seventeen different alleles were identified, namely eight missense variants, six nonsense variants, one in-frame deletion variant, and two splice site variants. For the latter two, minigene assays revealed aberrant transcripts containing frameshifts and premature termination codons. Structural analysis and molecular modeling suggested different degrees of structural destabilization for the missense variants. In conclusion, we report the largest cohort of patients with TULP1-associated IRD published to date. Most of the patients exhibited rod-driven disease, yet a fraction of the patients exhibited cone-driven disease. Our data support the hypothesis that TULP1 variants do not fold properly and thus trigger unfolded protein response, resulting in photoreceptor death.

Genome-Wide Association Study Identifies Two Common Loci Associated with Pigment Dispersion Syndrome/Pigmentary Glaucoma and Implicates Myopia in its Development.

PURPOSE: To identify genetic variants associated with pigment dispersion syndrome (PDS) and pigmentary glaucoma (PG) in unrelated patients and to further understand the genetic and potentially causal relationships between PDS and associated risk factors. DESIGN: A 2-stage genome-wide association meta-analysis with replication and subsequent in silico analyses including Mendelian randomization. PARTICIPANTS: A total of 574 cases with PG or PDS and 52 627 controls of European descent. METHODS: Genome-wide association analyses were performed in 4 cohorts and meta-analyzed in 3 stages: (1) a discovery meta-analysis was performed in 3 cohorts, (2) replication was performed in the fourth cohort, and (3) all 4 cohorts were meta-analyzed to increase statistical power. Two-sample Mendelian randomization was used to determine whether refractive error and intraocular pressure exert causal effects over PDS. MAIN OUTCOME MEASURES: The association of genetic variants with PDS and whether myopia exerts causal effects over PDS. RESULTS: Significant association was present at 2 novel loci for PDS/PG. These loci and follow-up analyses implicate the genes gamma secretase activator protein (GSAP) (lead single nucleotide polymorphism [SNP]: rs9641220, P = 6.0×10-10) and glutamate metabotropic receptor 5 (GRM5)/TYR (lead SNP: rs661177, P = 3.9×10-9) as important factors in disease risk. Mendelian randomization showed significant evidence that negative refractive error (myopia) exerts a direct causal effect over PDS (P = 8.86×10-7). CONCLUSIONS: Common SNPs relating to the GSAP and GRM5/TYR genes are associated risk factors for the development of PDS and PG. Although myopia is a known risk factor, this study uses genetic data to demonstrate that myopia is, in part, a cause of PDS and PG.

Molecular Properties of Human Guanylate Cyclase-Activating Protein 3 (GCAP3) and Its Possible Association with Retinitis Pigmentosa.

The cone-specific guanylate cyclase-activating protein 3 (GCAP3), encoded by the GUCA1C gene, has been shown to regulate the enzymatic activity of membrane-bound guanylate cyclases (GCs) in bovine and teleost fish photoreceptors, to an extent comparable to that of the paralog protein GCAP1. To date, the molecular mechanisms underlying GCAP3 function remain largely unexplored. In this work, we report a thorough characterization of the biochemical and biophysical properties of human GCAP3, moreover, we identified an isolated case of retinitis pigmentosa, in which a patient carried the c.301G>C mutation in GUCA1C, resulting in the substitution of a highly conserved aspartate residue by a histidine (p.(D101H)). We found that myristoylated GCAP3 can activate GC1 with a similar Ca2+-dependent profile, but significantly less efficiently than GCAP1. The non-myristoylated form did not induce appreciable regulation of GC1, nor did the p.D101H variant. GCAP3 forms dimers under physiological conditions, but at odds with its paralogs, it tends to form temperature-dependent aggregates driven by hydrophobic interactions. The peculiar properties of GCAP3 were confirmed by 2 ms molecular dynamics simulations, which for the p.D101H variant highlighted a very high structural flexibility and a clear tendency to lose the binding of a Ca2+ ion to EF3. Overall, our data show that GCAP3 has unusual biochemical properties, which make the protein significantly different from GCAP1 and GCAP2. Moreover, the newly identified point mutation resulting in a substantially unfunctional protein could trigger retinitis pigmentosa through a currently unknown mechanism.

A Novel, Apparently Silent Variant in MFSD8 Causes Neuronal Ceroid Lipofuscinosis with Marked Intrafamilial Variability.

Variants in MFSD8 can cause neuronal ceroid lipofuscinoses (NCLs) as well as nonsyndromic retinopathy. The mutation spectrum includes mainly missense and stop variants, but splice sites and frameshift variants have also been reported. To date, apparently synonymous substitutions have not been shown to cause MFSD8-associated diseases. We report two closely related subjects from a consanguineous Turkish family who presented classical features of NCLs but demonstrated marked intrafamilial variability in age at the onset and severity of symptoms. In fact, the difference in the onset of first neurologic symptoms was 15 years and that of ophthalmologic symptoms was 12 years. One subject presented an intellectual disability and a considerable cerebellar ataxia syndrome, while the other subject showed no intellectual disability and only a mild atactic syndrome. The diagnostic genetic testing of both subjects based on genome sequencing prioritized a novel, apparently synonymous variant in MFSD8, which was found in homozygosity in both subjects. The variant was not located within an integral part of the splice site consensus sequences. However, the bioinformatic analyses suggested that the mutant allele is more likely to cause exon skipping due to an altered ratio of exonic splice enhancer and silencer motifs. Exon skipping was confirmed in vitro by minigene assays and in vivo by RNA analysis from patient lymphocytes. The mutant transcript is predicted to result in a frameshift and, if translated, in a truncated protein. Synonymous variants are often given a low priority in genetic diagnostics because of their expected lack of functional impact. This study highlights the importance of investigating the impact of synonymous variants on splicing.