Multidisciplinary team directed analysis of whole genome sequencing reveals pathogenic non-coding variants in molecularly undiagnosed inherited retinal dystrophies

Clinician-Driven Reanalysis of Exome Sequencing Data From Patients With Inherited Retinal Diseases

Importance

Despite advances in next-generation sequencing (NGS), a significant proportion of patients with inherited retinal disease (IRD) remain undiagnosed after initial genetic testing. Exome sequencing (ES) reanalysis in the clinical setting has been suggested as one method for improving diagnosis of IRD.

Objective

To investigate the association of clinician-led reanalysis of ES data, which incorporates updated clinical information and comprehensive bioinformatic analysis, with the diagnostic yield in a cohort of patients with IRDs in Korea.

Design, Setting, and Participants

This was a multicenter prospective cohort study involving 264 unrelated patients with IRDs, conducted in Korea between March 2018 and February 2020. Comprehensive ophthalmologic examinations and ES analyses were performed, and ES data were reanalyzed by an IRD specialist for single nucleotide variants, copy number variants, mobile element insertions, and mitochondrial variants. Data were analyzed from March to July 2023.

Main Outcomes and Measures

Diagnostic rate of conventional bioinformatic analysis and clinician-driven ES reanalysis.

Results

A total of 264 participants (151 [57.2%] male; mean [SD] age at genetic testing, 33.6 [18.9] years) were enrolled, including 129 patients (48.9%) with retinitis pigmentosa and 26 patients (9.8%) with Stargardt disease or macular dystrophy. Initial bioinformatic analysis diagnosed 166 patients (62.9%). Clinician-driven reanalysis identified the molecular cause of diseases in an additional 22 patients, corresponding to an 8.3-percentage point increase in diagnostic rate. Key factors associated with new molecular diagnoses included clinical phenotype updates (4 patients) and detection of previously overlooked variation, such as structural variants (9 patients), mitochondrial variants (3 patients), filtered or not captured variants (4 patients), and noncanonical splicing variants (2 patients). Among the 22 patients, variants in 7 patients (31.8%) were observed in the initial analysis but not reported to patients, while those in the remaining 15 patients (68.2%) were newly detected by the ES reanalysis.

Conclusions and Relevance

In this cohort study, clinician-centered reanalysis of ES data was associated with improved molecular diagnostic yields in patients with IRD. This approach is important for uncovering missed genetic causes of retinal disease.

PHYH c.678+5G>T Leads to In-Frame Exon Skipping and Is Associated With Attenuated Refsum Disease

Purpose

To investigate the molecular effect of the variant PHYH:c.678+5G>T. This variant has conflicting interpretations in the ClinVar database and a maximum allele frequency of 0.0045 in the South Asian population in gnomAD.

Methods

We recruited patients from Moorfields Eye Hospital (London, UK) and Buenos Aires, Argentina, who were diagnosed with retinitis pigmentosa and found to have biallelic variants in PHYH, with at least one being c.678+5G>T. Total RNA was purified from PaxGene RNA-stabilized whole-blood samples, followed by reverse transcription to cDNA, PCR amplification of the canonical PHYH transcript, Oxford Nanopore Technologies library preparation, and single-molecule amplicon sequencing.

Results

Four patients provided a blood sample. One patient had isolated retinitis pigmentosa and three had mild extraocular findings. Blood phytanic acid levels were normal in two patients, mildly elevated in one, and markedly high in the fourth. Retinal evaluation showed an intact ellipsoid zone as well as preserved autofluorescence in the macular region in three of the four patients. In all patients, we observed in-frame skipping of exons 5 and 6 in 31.1% to 88.4% of the amplicons and a smaller proportion (0% to 11.3% of amplicons) skipping exon 6 only.

Conclusions

We demonstrate a significant effect of PHYH:c.678+5G>T on splicing of the canonical transcript. The in-frame nature of this may be in keeping with a mild presentation and higher prevalence in the general population. These data support the classification of the variant as pathogenic, and patients harboring a biallelic genotype should undergo phytanic acid testing.

Spectrum of Genetic Variants in the Most Common Genes Causing Inherited Retinal Disease in a Large Molecularly Characterized United Kingdom Cohort

PURPOSE

Inherited retinal disease (IRD) is a leading cause of blindness. Recent advances in gene-directed therapies highlight the importance of understanding the genetic basis of these disorders. This study details the molecular spectrum in a large United Kingdom (UK) IRD patient cohort.

DESIGN

Retrospective study of electronic patient records.

PARTICIPANTS

Patients with IRD who attended the Genetics Service at Moorfields Eye Hospital between 2003 and July 2020, in whom a molecular diagnosis was identified.

METHODS

Genetic testing was undertaken via a combination of single-gene testing, gene panel testing, whole exome sequencing, and more recently, whole genome sequencing. Likely disease-causing variants were identified from entries within the genetics module of the hospital electronic patient record (OpenEyes Electronic Medical Record). Analysis was restricted to only genes listed in the Genomics England PanelApp R32 Retinal Disorders panel (version 3.24), which includes 412 genes associated with IRD. Manual curation ensured consistent variant annotation and included only plausible disease-associated variants.

MAIN OUTCOME MEASURES

Detailed analysis was performed for variants in the 5 most frequent genes (ABCA4, USH2A, RPGR, PRPH2, and BEST1), as well as for the most common variants encountered in the IRD study cohort.

RESULTS

We identified 4415 individuals from 3953 families with molecularly diagnosed IRD (variants in 166 genes). Of the families, 42.7% had variants in 1 of the 5 most common IRD genes. Complex disease alleles contributed to disease in 16.9% of affected families with ABCA4-associated retinopathy. USH2A exon 13 variants were identified in 43% of affected individuals with USH2A-associated IRD. Of the RPGR variants, 71% were clustered in the ORF15 region. PRPH2 and BEST1 variants were associated with a range of dominant and recessive IRD phenotypes. Of the 20 most prevalent variants identified, 5 were not in the most common genes; these included founder variants in CNGB3, BBS1, TIMP3, EFEMP1, and RP1.

CONCLUSIONS

We describe the most common pathogenic IRD alleles in a large single-center multiethnic UK cohort and the burden of disease, in terms of families affected, attributable to these variants. Our findings will inform IRD diagnoses in future patients and help delineate the cohort of patients eligible for gene-directed therapies under development.

FINANCIAL DISCLOSURE(S)

Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Combining a prioritization strategy and functional studies nominates 5'UTR variants underlying inherited retinal disease

BACKGROUND

5' untranslated regions (5'UTRs) are essential modulators of protein translation. Predicting the impact of 5'UTR variants is challenging and rarely performed in routine diagnostics. Here, we present a combined approach of a comprehensive prioritization strategy and functional assays to evaluate 5'UTR variation in two large cohorts of patients with inherited retinal diseases (IRDs).

METHODS

We performed an isoform-level re-analysis of retinal RNA-seq data to identify the protein-coding transcripts of 378 IRD genes with highest expression in retina. We evaluated the coverage of their 5'UTRs by different whole exome sequencing (WES) kits. The selected 5'UTRs were analyzed in whole genome sequencing (WGS) and WES data from IRD sub-cohorts from the 100,000 Genomes Project (n = 2397 WGS) and an in-house database (n = 1682 WES), respectively. Identified variants were annotated for 5'UTR-relevant features and classified into seven categories based on their predicted functional consequence. We developed a variant prioritization strategy by integrating population frequency, specific criteria for each category, and family and phenotypic data. A selection of candidate variants underwent functional validation using diverse approaches.

RESULTS

Isoform-level re-quantification of retinal gene expression revealed 76 IRD genes with a non-canonical retina-enriched isoform, of which 20 display a fully distinct 5'UTR compared to that of their canonical isoform. Depending on the probe design, 3-20% of IRD genes have 5'UTRs fully captured by WES. After analyzing these regions in both cohorts, we prioritized 11 (likely) pathogenic variants in 10 genes (ARL3, MERTK, NDP, NMNAT1, NPHP4, PAX6, PRPF31, PRPF4, RDH12, RD3), of which 7 were novel. Functional analyses further supported the pathogenicity of three variants. Mis-splicing was demonstrated for the PRPF31:c.-9+1G>T variant. The MERTK:c.-125G>A variant, overlapping a transcriptional start site, was shown to significantly reduce both luciferase mRNA levels and activity. The RDH12:c.-123C>T variant was found in cis with the hypomorphic RDH12:c.701G>A (p.Arg234His) variant in 11 patients. This 5'UTR variant, predicted to introduce an upstream open reading frame, was shown to result in reduced RDH12 protein but unaltered mRNA levels.

CONCLUSIONS

This study demonstrates the importance of 5'UTR variants implicated in IRDs and provides a systematic approach for 5'UTR annotation and validation that is applicable to other inherited diseases.