Structure-based network analysis predicts pathogenic variants in human proteins associated with inherited retinal disease.

Advances in gene sequencing technologies have accelerated the identification of genetic variants, but better tools are needed to understand which are causal of disease. This would be particularly useful in fields where gene therapy is a potential therapeutic modality for a disease-causing variant such as inherited retinal disease (IRD). Here, we apply structure-based network analysis (SBNA), which has been successfully utilized to identify variant-constrained amino acid residues in viral proteins, to identify residues that may cause IRD if subject to missense mutation. SBNA is based entirely on structural first principles and is not fit to specific outcome data, which makes it distinct from other contemporary missense prediction tools. In 4 well-studied human disease-associated proteins (BRCA1, HRAS, PTEN, and ERK2) with high-quality structural data, we find that SBNA scores correlate strongly with deep mutagenesis data. When applied to 47 IRD genes with available high-quality crystal structure data, SBNA scores reliably identified disease-causing variants according to phenotype definitions from the ClinVar database. Finally, we applied this approach to 63 patients at Massachusetts Eye and Ear (MEE) with IRD but for whom no genetic cause had been identified. Untrained models built using SBNA scores and BLOSUM62 scores for IRD-associated genes successfully predicted the pathogenicity of novel variants (AUC = 0.851), allowing us to identify likely causative disease variants in 40 IRD patients. Model performance was further augmented by incorporating orthogonal data from EVE scores (AUC = 0.927), which are based on evolutionary multiple sequence alignments. In conclusion, SBNA can used to successfully identify variants as causal of disease in human proteins and may help predict variants causative of IRD in an unbiased fashion.

Coding and non-coding variants in the ciliopathy gene CFAP410 cause early-onset non-syndromic retinal degeneration.

Inherited retinal degenerations are blinding genetic disorders characterized by high genetic and phenotypic heterogeneity. The implementation of next-generation sequencing in routine diagnostics, together with advanced clinical phenotyping including multimodal retinal imaging, have contributed to the increase of reports describing novel genotype-phenotype associations and phenotypic expansions. In this study, we describe sixteen families with early-onset non-syndromic retinal degenerations in which affected probands carried rare bi-allelic variants in CFAP410, a ciliary gene previously associated with syndromic recessive Jeune syndrome. The most common retinal phenotypes were cone-rod and rod-cone dystrophies, but the clinical presentations were unified by their early onset as well as the severe impact on central visual function. Twelve variants were detected (three pathogenic, seven likely pathogenic, two of uncertain significance), eight of which were novel. One deep intronic change, c.373+91A>G, led to the creation of a cryptic splice acceptor site in intron four, followed by the inclusion of a 200- base pair pseudoexon and subsequent premature stop codon formation. To our knowledge this is the first likely pathogenic deep-intronic variant identified in this gene. Meta-analysis of all published and novel CFAP410 variants revealed no clear correlation between the severity of the CFAP410-associated phenotypes and the identified causal variants. This is supported by the fact that the frequently encountered missense variant p.(Arg73Pro), often found in syndromic cases, was also associated with non-syndromic retinal degeneration. This study expands the current knowledge of CFAP410-associated ciliopathy by enriching its mutational landscape and supports its association with non-syndromic retinal degeneration.

Generation of a human induced pluripotent stem cell line (OGIi001) from peripheral blood mononuclear cells of a healthy male donor.

We have successfully derived a novel human induced pluripotent stem cell (hiPSC) line using non-integrative Sendai virus. This hiPSC line was generated from a healthy male adult donor, aged 55, and subjected to thorough characterization and extensive quality control. The analysis confirmed the expression of undifferentiated stem cell markers, demonstrated the ability to differentiate into the three germ layers, and revealed the absence of any chromosomal abnormalities.

RP2 X-LINKED RETINITIS PIGMENTOSA CARRIER STATE PRESENTING WITH VASCULAR LEAKAGE AND UNILATERAL MACULAR ATROPHY.

PURPOSE: We describe the unusual clinical presentation of a 33-year-old woman subsequently identified as a carrier of RP2-associated X-linked retinitis pigmentosa. METHODS: Case report. RESULTS: A 33-year-old woman without a known family history of retinal disease presented with unilateral reduced visual acuity and central scotoma in the left eye. Examination showed underlying macular atrophy in the left eye and a bilateral tapetal-like reflex. Full-field electroretinogram was abnormal in the left eye but normal in the right eye. Notable findings on wide-field imaging included bilateral peripheral vascular leakage on fluorescein angiography and a bilaterally symmetric radial pattern of hyperfluorescence on fundus autofluorescence. Genetic testing demonstrated a pathogenic variant in the gene RP2 confirming that she was a carrier of X-linked retinitis pigmentosa. CONCLUSION: We describe clinical features of the carrier state of RP2-XLRP and expand potential findings to include peripheral vascular leakage. This case highlights the importance of awareness of the carrier state, particularly if a family history cannot be provided.

Structure-based network analysis predicts mutations associated with inherited retinal disease.

With continued advances in gene sequencing technologies comes the need to develop better tools to understand which mutations cause disease. Here we validate structure-based network analysis (SBNA)1,2 in well-studied human proteins and report results of using SBNA to identify critical amino acids that may cause retinal disease if subject to missense mutation. We computed SBNA scores for genes with high-quality structural data, starting with validating the method using 4 well-studied human disease-associated proteins. We then analyzed 47 inherited retinal disease (IRD) genes. We compared SBNA scores to phenotype data from the ClinVar database and found a significant difference between benign and pathogenic mutations with respect to network score. Finally, we applied this approach to 65 patients at Massachusetts Eye and Ear (MEE) who were diagnosed with IRD but for whom no genetic cause was found. Multivariable logistic regression models built using SBNA scores for IRD-associated genes successfully predicted pathogenicity of novel mutations, allowing us to identify likely causative disease variants in 37 patients with IRD from our clinic. In conclusion, SBNA can be meaningfully applied to human proteins and may help predict mutations causative of IRD.

Combined X-linked familial exudative vitreoretinopathy and retinopathy of prematurity phenotype in an infant with mosaic turner syndrome with ring X chromosome.

BACKGROUND: Retinopathy of prematurity (ROP) and familial exudative vitreoretinopathy (FEVR) are two distinct pathologies of retinal angiogenesis with overlapping clinical features. METHODS: Examination, multimodal imaging, and genetic testing were used to guide diagnosis and treatment. RESULTS: We report a combined phenotype of X-linked FEVR and ROP in a 4-month-old girl with mosaic Turner syndrome with ring X chromosome born at 26 weeks gestational age. She was initially diagnosed with atypical ROP with a vitreous band causing a localized traction retinal detachment, inferotemporal to the macula in the right eye, vessels to posterior zone 2 with no clear ridge temporally in the left eye, and fluorescein leakage in both eyes. Due to the suspicion of concurrent FEVR, genetic testing using a vitreoretinopathy panel was performed which revealed a mosaic Turner syndrome associated with 45,X/46,X,r(X), subsequently confirmed by chromosome analysis. The deleted region in the ring X chromosome included the NDP and RS1 genes. The patient was treated with laser photocoagulation of the peripheral avascular retina and sub-Tenon's triamcinolone injection in both eyes, intravitreal injection of bevacizumab in the left eye, and pars plicata vitrectomy in the right eye. CONCLUSIONS: In premature neonates with atypical ROP, a clinical suspicion of concurrent FEVR or similar vasculopathy is important and genetic testing may elucidate a genetic etiology, which could influence management and prognosis. Turner syndrome can be connected with co-occurring Mendelian gene disorders, particularly in individuals with mosaicism. The concurrence of FEVR and ROP appears to result in atypical and possibly more severe phenotypes.

Disparities in Inherited Retinal Degenerations.

To review disparities in the field of inherited retinal degenerations to establish foundations for future discussions oriented toward finding possible solutions. A narrative overview of the literature. Despite collective efforts towards democratization of genetic testing and investigation, genetic databases containing primarily European populations are heavily relied upon. Access to specialized care and other resources is also still not available to all. Recognizing and addressing disparities and inequities within the field of inherited retinal degenerations will improve our care of these patients and our knowledge of their conditions.

Unique Capabilities of Genome Sequencing for Rare Disease Diagnosis.

Background: Causal variants underlying rare disorders may remain elusive even after expansive gene panels or exome sequencing (ES). Clinicians and researchers may then turn to genome sequencing (GS), though the added value of this technique and its optimal use remain poorly defined. We therefore investigated the advantages of GS within a phenotypically diverse cohort. Methods: GS was performed for 744 individuals with rare disease who were genetically undiagnosed. Analysis included review of single nucleotide, indel, structural, and mitochondrial variants. Results: We successfully solved 218/744 (29.3%) cases using GS, with most solves involving established disease genes (157/218, 72.0%). Of all solved cases, 148 (67.9%) had previously had non-diagnostic ES. We systematically evaluated the 218 causal variants for features requiring GS to identify and 61/218 (28.0%) met these criteria, representing 8.2% of the entire cohort. These included small structural variants (13), copy neutral inversions and complex rearrangements (8), tandem repeat expansions (6), deep intronic variants (15), and coding variants that may be more easily found using GS related to uniformity of coverage (19). Conclusion: We describe the diagnostic yield of GS in a large and diverse cohort, illustrating several types of pathogenic variation eluding ES or other techniques. Our results reveal a higher diagnostic yield of GS, supporting the utility of a genome-first approach, with consideration of GS as a secondary or tertiary test when higher-resolution structural variant analysis is needed or there is a strong clinical suspicion for a condition and prior targeted genetic testing has been negative.

Identification of a novel large multigene deletion and a frameshift indel in PDE6B as the underlying cause of early-onset recessive rod-cone degeneration.

A family, with two affected identical twins with early-onset recessive inherited retinal degeneration, was analyzed to determine the underlying genetic cause of pathology. Exome sequencing revealed a rare and previously reported causative variant (c.1923_1969delinsTCTGGG; p.Asn643Glyfs*29) in the PDE6B gene in the affected twins and their unaffected father. Further investigation, using genome sequencing, identified a novel ∼7.5-kb deletion (Chr 4:670,405-677,862del) encompassing the ATP5ME gene, part of the 5' UTR of MYL5, and a 378-bp (Chr 4:670,405-670,782) region from the 3' UTR of PDE6B in the affected twins and their unaffected mother. Both variants segregated with disease in the family. Analysis of the relative expression of PDE6B, in peripheral blood cells, also revealed a significantly lower level of PDE6B transcript in affected siblings compared to a normal control. PDE6B is associated with recessive rod-cone degeneration and autosomal dominant congenital stationary night blindness. Ophthalmic evaluation of these patients showed night blindness, fundus abnormalities, and peripheral vision loss, which are consistent with PDE6B-associated recessive retinal degeneration. These findings suggest that the loss of PDE6B transcript resulting from the compound heterozygous pathogenic variants is the underlying cause of recessive rod-cone degeneration in the study family.

Novel RCBTB1 variants causing later-onset non-syndromic retinal dystrophy with macular chorioretinal atrophy.

BACKGROUND: Variants in RCBTB1 were recently described to cause a retinal dystrophy with only eight families described to date and a predominant phenotype of macular atrophy and peripheral reticular degeneration. Here, we further evaluate the genotypic and phenotypic characteristics of biallelic RCBTB1-associated retinal dystrophy in a North American clinic population. METHODS: A retrospective analysis of genetic and clinical features was performed in individuals with biallelic variants in RCBTB1. RESULTS: Three unrelated individuals of French-Canadian descent with rare biallelic RCBTB1 variants were identified. All individuals shared a novel p.(Ser342Leu) missense variant; one patient was homozygous whereas the other two each possessed a second unique novel variant p.(Gln120*) and p.(Pro224Leu). All three had macula-predominant disease with symptom onset in the fifth decade of life. CONCLUSION: This report adds to the genetic diversity of RCBTB1-associated disease. These cases confirm the later-onset, relative to many other retinal dystrophies, and macular focus of disease described in most cases to-date. They are thus a reminder of considering hereditary disease in the differential for later-onset macular atrophy.

A hidden structural variation in a known IRD gene: a cautionary tale of two new disease candidate genes.

Rod-cone dystrophy (RCD), also known as retinitis pigmentosa, is an inherited condition leading to vision loss, affecting 1 in 3500 people. More than 270 genes are known to be implicated in the inherited retinal degenerations (IRDs), yet genetic diagnosis for ∼30% of IRD of patients remains elusive despite advances in sequencing technologies. The goal of this study was to determine the genetic causality in a family with RCD. Family members were given a full ophthalmic exam at the Retinal Service at Massachusetts Eye and Ear and consented to genetic testing. Whole-exome sequencing (WES) was performed and variants of interest were Sanger-validated. Functional assays were conducted in zebrafish along with splicing assays in relevant cell lines to determine the impact on retinal function. WES identified variants in two potential candidate genes that segregated with disease: GNL3 (G Protein Nucleolar 3) c.1187 + 3A > C and c.1568-8C > A; and PDE4DIP (Phosphodiester 4D Interacting Protein) c.3868G > A (p.Glu1290Lys) and c.4603G > A (p.Ala1535Thr). Both genes were promising candidates based on their retinal involvement (development and interactions with IRD-associated proteins); however, the functional assays did not validate either gene. Subsequent WES reanalysis with an updated bioinformatics pipeline and widened search parameters led to the detection of a 94-bp duplication in PRPF31 (pre-mRNA Processing Factor 31) c.73_266dup (p.Asp56GlyfsTer33) as the causal variant. Our study demonstrates the importance of thorough functional characterization of new disease candidate genes and the value of reanalyzing next-generation sequencing sequence data, which in our case led to identification of a hidden pathogenic variant in a known IRD gene.

The importance of automation in genetic diagnosis: Lessons from analyzing an inherited retinal degeneration cohort with the Mendelian Analysis Toolkit (MATK).

PURPOSE: In Mendelian disease diagnosis, variant analysis is a repetitive, error-prone, and time consuming process. To address this, we have developed the Mendelian Analysis Toolkit (MATK), a configurable, automated variant ranking program. METHODS: MATK aggregates variant information from multiple annotation sources and uses expert-designed rules with parameterized weights to produce a ranked list of potentially causal solutions. MATK performance was measured by a comparison between MATK-aided and human-domain expert analyses of 1060 families with inherited retinal degeneration (IRD), analyzed using an IRD-specific gene panel (589 individuals) and exome sequencing (471 families). RESULTS: When comparing MATK-assisted analysis with expert curation in both the IRD-specific gene panel and exome sequencing (1060 subjects), 97.3% of potential solutions found by experts were also identified by the MATK-assisted analysis (541 solutions identified with MATK of 556 solutions found by conventional analysis). Furthermore, MATK-assisted analysis identified 114 additional potential solutions from the 504 cases unsolved by conventional analysis. CONCLUSION: MATK expedites the process of identification of likely solving variants in Mendelian traits, and reduces variability stemming from human error and researcher bias. MATK facilitates data reanalysis to keep up with the constantly improving annotation sources and next-generation sequencing processing pipelines. The software is open source and available at https://gitlab.com/matthew_maher/mendelanalysis.

Beyond Sector Retinitis Pigmentosa: Expanding the Phenotype and Natural History of the Rhodopsin Gene Codon 106 Mutation (Gly-to-Arg) in Autosomal Dominant Retinitis Pigmentosa.

Sector and pericentral are two rare, regional forms of retinitis pigmentosa (RP). While usually defined as stable or only very slowly progressing, the available literature to support this claim is limited. Additionally, few studies have analyzed the spectrum of disease within a particular genotype. We identified all cases (9 patients) with an autosomal dominant Rhodopsin variant previously associated with sector RP (RHO c.316G > A, p.Gly106Arg) at our institution. Clinical histories were reviewed, and testing included visual fields, multimodal imaging, and electroretinography. Patients demonstrated a broad phenotypic spectrum that spanned regional phenotypes from sector-like to pericentral RP, as well as generalized disease. We also present evidence of significant intrafamilial variability in regional phenotypes. Finally, we present the longest-reported follow-up for a patient with RHO-associated sector-like RP, showing progression from sectoral to pericentral disease over three decades. In the absence of comorbid macular disease, the long-term prognosis for central visual acuity is good. However, we found that significant progression of RHO p.Gly106Arg disease can occur over protracted periods, with impact on peripheral vision. Longitudinal widefield imaging and periodic ERG reassessment are likely to aid in monitoring disease progression.

Broadening INPP5E phenotypic spectrum: detection of rare variants in syndromic and non-syndromic IRD.

Pathogenic variants in INPP5E cause Joubert syndrome (JBTS), a ciliopathy with retinal involvement. However, despite sporadic cases in large cohort sequencing studies, a clear association with non-syndromic inherited retinal degenerations (IRDs) has not been made. We validate this association by reporting 16 non-syndromic IRD patients from ten families with bi-allelic mutations in INPP5E. Additional two patients showed early onset IRD with limited JBTS features. Detailed phenotypic description for all probands is presented. We report 14 rare INPP5E variants, 12 of which have not been reported in previous studies. We present tertiary protein modeling and analyze all INPP5E variants for deleteriousness and phenotypic correlation. We observe that the combined impact of INPP5E variants in JBTS and non-syndromic IRD patients does not reveal a clear genotype-phenotype correlation, suggesting the involvement of genetic modifiers. Our study cements the wide phenotypic spectrum of INPP5E disease, adding proof that sequence defects in this gene can lead to early-onset non-syndromic IRD.

Characterization of the Spectrum of Ophthalmic Changes in Patients With Alagille Syndrome.

Purpose: The purpose of this study was to characterize the phenotypic spectrum of ophthalmic findings in patients with Alagille syndrome. Methods: We conducted a retrospective, observational, multicenter, study on 46 eyes of 23 subjects with Alagille syndrome. We reviewed systemic and ophthalmologic data extracted from medical records, color fundus photography, fundus autofluorescence, optical coherence tomography, visual fields, electrophysiological assessments, and molecular genetic findings. Results: Cardiovascular abnormalities were found in 83% of all cases (of those, 74% had cardiac murmur), whereas 61% had a positive history of hepatobiliary issues, and musculoskeletal anomalies were present in 61% of all patients. Dysmorphic facies were present in 16 patients, with a broad forehead being the most frequent feature. Ocular symptoms were found in 91%, with peripheral vision loss being the most frequent complaint. Median (range) Snellen visual acuity of all eyes was 20/25 (20/20 to hand motion [HM]). Anterior segment abnormalities were present in 74% of the patients; of those, posterior embryotoxon was the most frequent finding. Abnormalities of the optic disc were found in 52%, and peripheral retinal abnormalities were the most frequent ocular finding in this series, found in 96% of all patients. Fifteen JAG1 mutations were identified in 16 individuals; of those, 6 were novel. Conclusions: This study reports a cohort of patients with Alagille syndrome in which peripheral chorioretinal changes were more frequent than posterior embryotoxon, the most frequent ocular finding according to a number of previous studies. We propose that these peripheral chorioretinal changes are a new hallmark to help diagnose this syndrome.

Expanding the phenotypic spectrum in RDH12-associated retinal disease.

Retinol dehydrogenase 12, RDH12, plays a pivotal role in the visual cycle to ensure the maintenance of normal vision. Alterations in activity of this protein result in photoreceptor death and decreased vision beginning at an early age and progressing to substantial vision loss later in life. Here we describe 11 patients with retinal degeneration that underwent next-generation sequencing (NGS) with a targeted panel of all currently known inherited retinal degeneration (IRD) genes and whole-exome sequencing to identify the genetic causality of their retinal disease. These patients display a range of phenotypic severity prompting clinical diagnoses of macular dystrophy, cone-rod dystrophy, retinitis pigmentosa, and early-onset severe retinal dystrophy all attributed to biallelic recessive mutations in RDH12 We report 15 causal alleles and expand the repertoire of known RDH12 mutations with four novel variants: c.215A > G (p.Asp72Gly); c.362T > C (p.Ile121Thr); c.440A > C (p.Asn147Thr); and c.697G > A (p.Val233Ille). The broad phenotypic spectrum observed with biallelic RDH12 mutations has been observed in other genetic forms of IRDs, but the diversity is particularly notable here given the prior association of RDH12 primarily with severe early-onset disease. This breadth emphasizes the importance of broad genetic testing for inherited retinal disorders and extends the pool of individuals who may benefit from imminent gene-targeted therapies.

Copy-number variation contributes 9% of pathogenicity in the inherited retinal degenerations.

PURPOSE: Current sequencing strategies can genetically solve 55-60% of inherited retinal degeneration (IRD) cases, despite recent progress in sequencing. This can partially be attributed to elusive pathogenic variants (PVs) in known IRD genes, including copy-number variations (CNVs), which have been shown as major contributors to unsolved IRD cases. METHODS: Five hundred IRD patients were analyzed with targeted next-generation sequencing (NGS). The NGS data were used to detect CNVs with ExomeDepth and gCNV and the results were compared with CNV detection with a single-nucleotide polymorphism (SNP) array. Likely causal CNV predictions were validated by quantitative polymerase chain reaction (qPCR). RESULTS: Likely disease-causing single-nucleotide variants (SNVs) and small indels were found in 55.6% of subjects. PVs in USH2A (11.6%), RPGR (4%), and EYS (4%) were the most common. Likely causal CNVs were found in an additional 8.8% of patients. Of the three CNV detection methods, gCNV showed the highest accuracy. Approximately 30% of unsolved subjects had a single likely PV in a recessive IRD gene. CONCLUSION: CNV detection using NGS-based algorithms is a reliable method that greatly increases the genetic diagnostic rate of IRDs. Experimentally validating CNVs helps estimate the rate at which IRDs might be solved by a CNV plus a more elusive variant.

Predictive value of genetic testing for inherited retinal diseases in patients with suspected atypical autoimmune retinopathy.

PURPOSE: The clinical features of autoimmune retinopathy (AIR) can resemble and be difficult to differentiate from inherited retinal degenerations (IRDs). Misdiagnosis of an IRD as AIR causes unnecessary treatment with immunosuppressive agents. The purpose of this study is to calculate the predictive value of genetic testing for IRDs in patients with suspected AIR and provide clinical examples where genetic testing has been useful. METHODS: We identified patients seen at MEEI between April 2013 and January 2017 for whom the differentiation of AIR vs. IRDs was difficult based on clinical assessment alone. All patients had some atypical features for AIR, but tested positive for anti-retinal antibodies. Within this group, we identified six patients who had genetic testing for IRDs with the Genetic Eye Disease panel for retinal genes (GEDi-R). We calculated the positive predictive value (PPV) and negative predictive value (NPV) of genetic testing in a population with approximately equal numbers of IRD and AIR patients. RESULTS: Six patients had clinical features that made distinguishing between IRDs and AIR on a clinical basis difficult and were sent for genetic testing: four women and two men with a mean age of 59.5 years. In two of these six patients, genetic diagnoses were made based upon the identification of known pathogenic variants in the common IRD genes USH2A and RHO. Two patients had variants of unknown significance within genes associated with IRDs, and the other two had no relevant genetic findings. Given the 60% sensitivity and 3% false positive rate for GEDi-R testing and assuming a 50% pre-test probability of having an IRD, the PPV for GEDi-R for detecting IRD is 95.2% and the NPV is 70.8%. CONCLUSIONS AND IMPORTANCE: In patients for whom the differential diagnosis of AIR and IRDs is unclear based on clinical information, genetic testing can be a valuable tool when it identifies an IRD, sparing the patient unnecessary immunosuppressive treatment. However, the test has a low NPV so a negative genetic testing result does not confidently exclude IRD as the true diagnosis.

Contribution of noncoding pathogenic variants to RPGRIP1-mediated inherited retinal degeneration.

PURPOSE: With the advent of gene therapies for inherited retinal degenerations (IRDs), genetic diagnostics will have an increasing role in clinical decision-making. Yet the genetic cause of disease cannot be identified using exon-based sequencing for a significant portion of patients. We hypothesized that noncoding pathogenic variants contribute significantly to the genetic causality of IRDs and evaluated patients with single coding pathogenic variants in RPGRIP1 to test this hypothesis. METHODS: IRD families underwent targeted panel sequencing. Unsolved cases were explored by exome and genome sequencing looking for additional pathogenic variants. Candidate pathogenic variants were then validated by Sanger sequencing, quantitative polymerase chain reaction, and in vitro splicing assays in two cell lines analyzed through amplicon sequencing. RESULTS: Among 1722 families, 3 had biallelic loss-of-function pathogenic variants in RPGRIP1 while 7 had a single disruptive coding pathogenic variants. Exome and genome sequencing revealed potential noncoding pathogenic variants in these 7 families. In 6, the noncoding pathogenic variants were shown to lead to loss of function in vitro. CONCLUSION: Noncoding pathogenic variants were identified in 6 of 7 families with single coding pathogenic variants in RPGRIP1. The results suggest that noncoding pathogenic variants contribute significantly to the genetic causality of IRDs and RPGRIP1-mediated IRDs are more common than previously thought.