Splicing mutations in inherited retinal diseases

Variants in UBAP1L lead to autosomal recessive rod-cone and cone-rod dystrophy

PURPOSE

Progressive inherited retinal degenerations (IRDs) affecting rods and cones are clinically and genetically heterogeneous and can lead to blindness with limited therapeutic options. The major gene defects have been identified in subjects of European and Asian descent with only few reports of North African descent.

METHODS

Genome, targeted next-generation, and Sanger sequencing was applied to cohort of ∼4000 IRDs cases. Expression analyses were performed including Chip-seq database analyses, on human-derived retinal organoids (ROs), retinal pigment epithelium cells, and zebrafish. Variants' pathogenicity was accessed using 3D-modeling and/or ROs.

RESULTS

Here, we identified a novel gene defect with three distinct pathogenic variants in UBAP1L in 4 independent autosomal recessive IRD cases from Tunisia. UBAP1L is expressed in the retinal pigment epithelium and retina, specifically in rods and cones, in line with the phenotype. It encodes Ubiquitin-associated protein 1-like, containing a solenoid of overlapping ubiquitin-associated domain, predicted to interact with ubiquitin. In silico and in vitro studies, including 3D-modeling and ROs revealed that the solenoid of overlapping ubiquitin-associated domain is truncated and thus ubiquitin binding most likely abolished secondary to all variants identified herein.

CONCLUSION

Biallelic UBAP1L variants are a novel cause of IRDs, most likely enriched in the North African population.

Splicing defects and CRISPR-Cas9 correction in isogenic homozygous photoreceptor precursors harboring clustered deep-intronic ABCA4 variants

Splicing defects from deep-intronic variants significantly contribute to the mutational spectrum in ABCA4-associated inherited retinal diseases, necessitating functional validation for their pathological classification. Typically, minigene assays in HEK293(T) can qualitatively assess splicing defects, yet they often fail to quantitatively reproduce the resulting mis-splicing patterns, leaving uncertainty on severity and pathogenicity. As a potential cellular model derived from patient cells, photoreceptor precursor cells (PPCs) play a pivotal role in assessing the severity of specific splicing mutations. Nevertheless, the accessibility of biosamples is commonly constrained, and their establishment is costly and laborious. In this study, we combined and investigated the use of a minigene assay and isogenic PPCs, as superior qualitative and more accessible cellular models for the assessment of splicing defects. Specifically, we focused on the clustered c.5196+1013A>G, c.5196+1056A>G, and c.5196+1216C>A deep-intronic variants in intron 36 of ABCA4, comparing their resulting (mis)splicing patterns in minigene-transfected cells and isogenic CRISPR-Cas9-knocked-in PPCs harboring these pathogenic variants in homozygous state. Moreover, we demonstrate the successful correction of these three splicing defects in homozygous mutant PPCs using a single pair of guide RNAs to target Cas9 cleavage, thereby identifying an efficient gene editing strategy for therapeutic applications.

Whole genome sequencing enables new genetic diagnosis for inherited retinal diseases by identifying pathogenic variants

Inherited retinal diseases (IRDs) are a group of common primary retinal degenerative disorders. Conventional genetic testing strategies, such as panel-based sequencing and whole exome sequencing (WES), can only elucidate the genetic etiology in approximately 60% of IRD patients. Studies have suggested that unsolved IRD cases could be attributed to previously undetected structural variants (SVs) and intronic variants in IRD-related genes. The aim of our study was to obtain a definitive genetic diagnosis by employing whole genome sequencing (WGS) in IRD cases where the causative genes were inconclusive following an initial screening by panel sequencing. A total of 271 unresolved IRD patients and their available family members (n = 646) were screened using WGS to identify pathogenic SVs and intronic variants in 792 known ocular disease genes. Overall, 13% (34/271) of IRD patients received a confirmed genetic diagnosis, among which 7% were exclusively attributed to SVs, 4% to a combination of single nucleotide variants (SNVs) and SVs while another 2% were linked to intronic variants. 22 SVs, 3 deep-intronic variants, and 2 non-canonical splice-site variants across 14 IRD genes were identified in the entire cohort. Notably, all of these detected SVs and intronic variants were novel pathogenic variants. Among those, 74% (20/27) of variants were found in genes causally linked to Retinitis Pigmentosa (RP), with the gene EYS being the most frequently affected by SVs. The identification of SVs and intronic variants through WGS enhances the genetic diagnostic yield of IRDs and broadens the mutational spectrum of known IRD-associated genes.

Comprehensive genetic analysis reveals the mutational landscape of ABCA4-associated retinal dystrophy in a Chinese cohort

PURPOSE

To depict the variant profiles of the ABCA4 gene in a large Chinese cohort of patients with ABCA4-associated retinal dystrophy (ABCA4-RD).

METHODS

We recruited 290 unrelated Chinese patients with ABCA4-RD and did ABCA4 mutational screening by a combination of Sanger sequencing, targeted exome sequencing, entire ABCA4 locus sequencing, and whole genome sequencing (WGS). The pathogenicity of variants was assessed using in silico tools or in vitro splicing assays following the American College of Medical Genetics and Genomics guidelines.

RESULTS

Two hundred sixty-eight distinct pathogenic variants were identified, and 57 were novel. In 580 alleles, 22 noncoding region variants outside canonical splice sites and 4 structural variations were found in 44 alleles accounting for 7.6% of all alleles. Bioinformatics analysis showed the complex mechanism of aberrant splicing productsnatural splice site disruption, branch point destruction, and cryptic splice site activation. Correspondingly, minigene assays validated the various abnormal splicing products, including exon skipping, exon elongation, partial exon deletion, and pseudoexon insertion. WGS identified the first inversion variation in ABCA4.

CONCLUSIONS

This study systematically depicted the variant profiles of ABCA4 and revealed the missing alleles of patients with ABCA4-RD in a large Chinese cohort. Our findings demonstrated the complexity of molecular diagnosis of Mendelian diseases and the efficiency of WGS for detecting structural variants.

Association Study Between Kynurenine 3-Monooxygenase (KMO) Gene and Parkinson's Disease Patients

The influence of various risk factors such as aging, intricate cellular molecular processes, and lifestyle factors like smoking, alcohol consumption, caffeine intake, and occupational factors has received increased focus in relation to the risk and development of Parkinson's disease (PD). Limited research has been conducted on the assessment of lifestyle impact on kynurenine 3-monooxygenase (KMO) gene in PD. A total of 164 subjects, including 82 PD cases and 82 healthy individuals, were recruited based on specific inclusion and exclusion criteria. The severity of PD and clinical assessment were evaluated using the Unified Parkinson's Disease Rating Scale (UPDRS) and Hoehn and Yahr (HY) scaling. Sanger sequencing was performed to analyse the KMO gene in the recruited subjects, and case-control studies were conducted. The UPDRS assessment revealed significant impairments in smell, tremors, walking, and posture instability in the late-onset PD cohorts. The HY scaling indicated a higher proportion of late-onset cohorts in stage 2. Moreover, both alcoholic and non-alcoholic groups showed significantly increased levels of 3-HK in late-onset PD. Gene analysis identified missense variants at position g.241593373 T > A (rs752312199) and intronic variants at positions g.241592623A > G (rs640718), g.241592800C > A (rs990388262), g.241592802A > C (rs1350160268), g.241592808 T > C (rs1478255936), and g.241592812G > T (rs948928931). The alterations in the KMO gene were found to influence the levels of kynurenic acid (KYNA) and 3-hydroxykynurenine (3-HK). Genomic analysis revealed a high prevalence of missense mutations in the late-onset PD groups, leading to a decline in 3-HK levels in patients. This leads to the reduction of the progression of disease in late-onset groups which shows that this mutation may lead to the protective effect on the PD subjects. This study suggests the use of KYNA and 3-HK as potential biomarkers in analysing the progression of disease. This study is limited by its small sample size. To overcome this limitation, a larger study involving in greater number of participants is needed to thoroughly investigate the KMO gene and KP metabolites, to enhance our understanding of Parkinson's disease progression, and to enhance diagnostic capabilities.

EyeG2P: an automated variant filtering approach improves efficiency of diagnostic genomic testing for inherited ophthalmic disorders

BACKGROUND

Genomic variant prioritisation is one of the most significant bottlenecks to mainstream genomic testing in healthcare. Tools to improve precision while ensuring high recall are critical to successful mainstream clinical genomic testing, in particular for whole genome sequencing where millions of variants must be considered for each patient.

METHODS

We developed EyeG2P, a publicly available database and web application using the Ensembl Variant Effect Predictor. EyeG2P is tailored for efficient variant prioritisation for individuals with inherited ophthalmic conditions. We assessed the sensitivity of EyeG2P in 1234 individuals with a broad range of eye conditions who had previously received a confirmed molecular diagnosis through routine genomic diagnostic approaches. For a prospective cohort of 83 individuals, we assessed the precision of EyeG2P in comparison with routine diagnostic approaches. For 10 additional individuals, we assessed the utility of EyeG2P for whole genome analysis.

RESULTS

EyeG2P had 99.5% sensitivity for genomic variants previously identified as clinically relevant through routine diagnostic analysis (n=1234 individuals). Prospectively, EyeG2P enabled a significant increase in precision (35% on average) in comparison with routine testing strategies (p<0.001). We demonstrate that incorporation of EyeG2P into whole genome sequencing analysis strategies can reduce the number of variants for analysis to six variants, on average, while maintaining high diagnostic yield.

CONCLUSION

Automated filtering of genomic variants through EyeG2P can increase the efficiency of diagnostic testing for individuals with a broad range of inherited ophthalmic disorders.

Definition of the transcriptional units of inherited retinal disease genes by meta-analysis of human retinal transcriptome data

BACKGROUND

Inherited retinal diseases (IRD) are genetically heterogeneous disorders that cause the dysfunction or loss of photoreceptor cells and ultimately lead to blindness. To date, next-generation sequencing procedures fail to detect pathogenic sequence variants in coding regions of known IRD disease genes in about 30-40% of patients. One of the possible explanations for this missing heritability is the presence of yet unidentified transcripts of known IRD genes. Here, we aimed to define the transcript composition of IRD genes in the human retina by a meta-analysis of publicly available RNA-seq datasets using an ad-hoc designed pipeline.

RESULTS

We analysed 218 IRD genes and identified 5,054 transcripts, 3,367 of which were not previously reported. We assessed their putative expression levels and focused our attention on 435 transcripts predicted to account for at least 5% of the expression of the corresponding gene. We looked at the possible impact of the newly identified transcripts at the protein level and experimentally validated a subset of them.

CONCLUSIONS

This study provides an unprecedented, detailed overview of the complexity of the human retinal transcriptome that can be instrumental in contributing to the resolution of some cases of missing heritability in IRD patients.

A presumed missense variant in the U2AF2 gene causes exon skipping in neurodevelopmental diseases

U2 small nuclear RNA auxiliary factor 2 (U2AF2) is an indispensable pre-mRNA splicing factor in the early process of splicing. Recently, U2AF2 was reported as a novel candidate gene associated with neurodevelopmental disorders. Herein, we report a patient with a novel presumed heterozygous missense variant in the U2AF2 gene (c.603G>T), who has a similar clinical phenotype as the patient reported before, including epilepsy, intellectual disability, language delay, microcephaly, and hypoplastic corpus callosum. We reviewed the phenotypic and genetic spectrum of patients with U2AF2-related neurological diseases, both newly diagnosed and previously reported. To investigate the possible pathogenesis, EBV-immortalized lymphoblastoid cells were derived from the peripheral blood obtained from the patient and control groups. Furthermore, according to the results of WB, RT-PCR, Q-PCR, and cDNA sequencing of RT-PCR products, the presumed missense variant c.603G>T caused exon 6 skipping in the U2AF2 mRNA transcript and led to a truncated protein (p.E163_E201del). Cell Counting Kit-8 (CCK-8) and cell cycle detection demonstrated that the variant c.603G>T inhibited the proliferation of patient lymphocyte cells compared with the control group. This study is aimed at expanding the phenotypic and genetic spectrum of U2AF2-related neurodevelopmental diseases and investigating the potential effects. This is the first report of the possible pathogenesis of a U2AF2 gene pathogenic variant in a patient with neurodevelopmental diseases and shows that a novel presumed missense variant in the U2AF2 gene causes exon skipping.

Recommendations for clinical interpretation of variants found in non-coding regions of the genome

BACKGROUND

The majority of clinical genetic testing focuses almost exclusively on regions of the genome that directly encode proteins. The important role of variants in non-coding regions in penetrant disease is, however, increasingly being demonstrated, and the use of whole genome sequencing in clinical diagnostic settings is rising across a large range of genetic disorders. Despite this, there is no existing guidance on how current guidelines designed primarily for variants in protein-coding regions should be adapted for variants identified in other genomic contexts.

METHODS

We convened a panel of nine clinical and research scientists with wide-ranging expertise in clinical variant interpretation, with specific experience in variants within non-coding regions. This panel discussed and refined an initial draft of the guidelines which were then extensively tested and reviewed by external groups.

RESULTS

We discuss considerations specifically for variants in non-coding regions of the genome. We outline how to define candidate regulatory elements, highlight examples of mechanisms through which non-coding region variants can lead to penetrant monogenic disease, and outline how existing guidelines can be adapted for the interpretation of these variants.

CONCLUSIONS

These recommendations aim to increase the number and range of non-coding region variants that can be clinically interpreted, which, together with a compatible phenotype, can lead to new diagnoses and catalyse the discovery of novel disease mechanisms.

Phenotype-Based Genetic Analysis Reveals Missing Heritability of ABCA4-Related Retinopathy: Deep Intronic Variants and Copy Number Variations

Purpose

To identify the missing heritability of ABCA4-related retinopathy in a Chinese cohort.

Methods

We recruited 33 unrelated patients with ABCA4-related retinopathy carrying a monoallelic variant in ABCA4. All patients underwent ophthalmic examinations. Next-generation sequencing of the whole ABCA4 sequence, including coding and noncoding regions, was performed to detect deep intronic variants (DIVs) and copy number variations (CNVs).

Results

We identified eight missing pathogenic ABCA4 variants in 60.6% of the patients (20/33), which comprised five DIVs and three CNVs. The five DIVs, including four novel (c.1555-816T>G, c.2919-169T>G, c.2919-884G>T, and c.5461-1321A>G) and one reported (c.4539+1100A>G), accounted for the missing alleles in 51.5% of the patients. Minigene assays showed that four novel DIVs activated cryptic splice sites leading to the insertions of pseudoexons. The three novel CNVs consisted of one gross deletion of 1273 bp (exon 2) and two gross duplications covering 25.2 kb (exons 28-43) and 9.4 kb (exons 38-44). The microhomology domains were identified at the breakpoints and revealed the potential mechanisms of CNV formation.

Conclusions

DIVs and CNVs explained approximately two-thirds of the unresolved Chinese cases with ABCA4-related retinopathy. Combining results from phenotypic-directed screening, targeting the whole ABCA4 sequencing and in silico tools can help to identify the missing heritability.

Whole genome sequencing and in vitro splice assays reveal genetic causes for inherited retinal diseases

Inherited retinal diseases (IRDs) are a major cause of visual impairment. These clinically heterogeneous disorders are caused by pathogenic variants in more than 270 genes. As 30-40% of cases remain genetically unexplained following conventional genetic testing, we aimed to obtain a genetic diagnosis in an IRD cohort in which the genetic cause was not found using whole-exome sequencing or targeted capture sequencing. We performed whole-genome sequencing (WGS) to identify causative variants in 100 unresolved cases. After initial prioritization, we performed an in-depth interrogation of all noncoding and structural variants in genes when one candidate variant was detected. In addition, functional analysis of putative splice-altering variants was performed using in vitro splice assays. We identified the genetic cause of the disease in 24 patients. Causative coding variants were observed in genes such as ATXN7, CEP78, EYS, FAM161A, and HGSNAT. Gene disrupting structural variants were also detected in ATXN7, PRPF31, and RPGRIP1. In 14 monoallelic cases, we prioritized candidate noncanonical splice sites or deep-intronic variants that were predicted to disrupt the splicing process based on in silico analyses. Of these, seven cases were resolved as they carried pathogenic splice defects. WGS is a powerful tool to identify causative variants residing outside coding regions or heterozygous structural variants. This approach was most efficient in cases with a distinct clinical diagnosis. In addition, in vitro splice assays provide important evidence of the pathogenicity of rare variants.

Comparison of in silico strategies to prioritize rare genomic variants impacting RNA splicing for the diagnosis of genomic disorders

The development of computational methods to assess pathogenicity of pre-messenger RNA splicing variants is critical for diagnosis of human disease. We assessed the capability of eight algorithms, and a consensus approach, to prioritize 249 variants of uncertain significance (VUSs) that underwent splicing functional analyses. The capability of algorithms to differentiate VUSs away from the immediate splice site as being 'pathogenic' or 'benign' is likely to have substantial impact on diagnostic testing. We show that SpliceAI is the best single strategy in this regard, but that combined usage of tools using a weighted approach can increase accuracy further. We incorporated prioritization strategies alongside diagnostic testing for rare disorders. We show that 15% of 2783 referred individuals carry rare variants expected to impact splicing that were not initially identified as 'pathogenic' or 'likely pathogenic'; one in five of these cases could lead to new or refined diagnoses.

Functional Evaluation of Splicing for Variants of Uncertain Significance in Patients with Inherited Retinal Diseases

Inherited retinal diseases (IRD) comprise a heterogeneous set of clinical and genetic disorders that lead to blindness. Given the emerging opportunities in precision medicine and gene therapy, it has become increasingly important to determine whether DNA variants with uncertain significance (VUS) are responsible for patients’ IRD. This research was performed to assess the functional consequence of six VUS identified in patients with IRD. Clinical assessments included an ophthalmic examination, best-corrected visual acuity, and kinetic perimetry. Imaging was acquired with the Optos ultra-widefield camera and spectral domain optical coherence tomography (SD-OCT). Genetic testing was performed by Molecular Vision Laboratories. VUS that were predicted to alter splicing were analyzed with a minigene assay, which revealed that VUS in the genes OPA1, CNGB1, and CLUAP1 altered spicing mechanisms. Due to emerging gene and cell therapies, these results expand the genotype-phenotype correlations for patients diagnosed with an IRD.

Spectrum of splicing variants in disease genes and the ability of RNA analysis to reduce uncertainty in clinical interpretation

The complexities of gene expression pose challenges for the clinical interpretation of splicing variants. To better understand splicing variants and their contribution to hereditary disease, we evaluated their prevalence, clinical classifications, and associations with diseases, inheritance, and functional characteristics in a 689,321-person clinical cohort and two large public datasets. In the clinical cohort, splicing variants represented 13% of all variants classified as pathogenic (P), likely pathogenic (LP), or variants of uncertain significance (VUSs). Most splicing variants were outside essential splice sites and were classified as VUSs. Among all individuals tested, 5.4% had a splicing VUS. If RNA analysis were to contribute supporting evidence to variant interpretation, we estimated that splicing VUSs would be reclassified in 1.7% of individuals in our cohort. This would result in a clinically significant result (i.e., P/LP) in 0.1% of individuals overall because most reclassifications would change VUSs to likely benign. In ClinVar, splicing VUSs were 4.8% of reported variants and could benefit from RNA analysis. In the Genome Aggregation Database (gnomAD), splicing variants comprised 9.4% of variants in protein-coding genes; most were rare, precluding unambiguous classification as benign. Splicing variants were depleted in genes associated with dominant inheritance and haploinsufficiency, although some genes had rare variants at essential splice sites or had common splicing variants that were most likely compatible with normal gene function. Overall, we describe the contribution of splicing variants to hereditary disease, the potential utility of RNA analysis for reclassifying splicing VUSs, and how natural variation may confound clinical interpretation of splicing variants.

Activation of cryptic donor splice sites by non-coding and coding PAX6 variants contributes to congenital aniridia

BACKGROUND

The paired-domain transcription factor paired box gene 6 (PAX6) causes a wide spectrum of ocular developmental anomalies, including congenital aniridia, Peters anomaly and microphthalmia. Here, we aimed to functionally assess the involvement of seven potentially non-canonical splicing variants on missplicing of exon 6, which represents the main hotspot region for loss-of-function PAX6 variants.

METHODS

By locus-specific analysis of PAX6 using Sanger and/or targeted next-generation sequencing, we screened a Spanish cohort of 106 patients with PAX6-related diseases. Functional splicing assays were performed by in vitro minigene approaches or directly in RNA from patient-derived lymphocytes cell line, when available.

RESULTS

Five out seven variants, including three synonymous changes, one small exonic deletion and one non-canonical splice variant, showed anomalous splicing patterns yielding partial exon skipping and/or elongation.

CONCLUSION

We describe new spliceogenic mechanisms for PAX6 variants mediated by creating or strengthening five different cryptic donor sites at exon 6. Our work revealed that the activation of cryptic PAX6 splicing sites seems to be a recurrent and underestimated cause of aniridia. Our findings pointed out the importance of functional assessment of apparently silent PAX6 variants to uncover hidden genetic alterations and to improve variant interpretation for genetic counselling in aniridia.

Phenotype Driven Analysis of Whole Genome Sequencing Identifies Deep Intronic Variants that Cause Retinal Dystrophies by Aberrant Exonization

Purpose

To demonstrate the effectiveness of combining retinal phenotyping and focused variant filtering from genome sequencing (GS) in identifying deep intronic disease causing variants in inherited retinal dystrophies.

Methods

Affected members from three pedigrees with classical enhanced S-cone syndrome (ESCS; Pedigree 1), congenital stationary night blindness (CSNB; Pedigree 2), and achromatopsia (ACHM; Pedigree 3), respectively, underwent detailed ophthalmologic evaluation, optical coherence tomography, and electroretinography. The probands underwent panel-based genetic testing followed by GS analysis. Minigene constructs (NR2E3, GPR179 and CNGB3) and patient-derived cDNA experiments (NR2E3 and GPR179) were performed to assess the functional effect of the deep intronic variants.

Results

The electrophysiological findings confirmed the clinical diagnosis of ESCS, CSNB, and ACHM in the respective pedigrees. Panel-based testing revealed heterozygous pathogenic variants in NR2E3 (NM_014249.3; c.119-2A>C; Pedigree 1) and CNGB3 (NM_019098.4; c.1148delC/p.Thr383Ilefs*13; Pedigree 3). The GS revealed heterozygous deep intronic variants in Pedigrees 1 (NR2E3; c.1100+1124G>A) and 3 (CNGB3; c.852+4751A>T), and a homozygous GPR179 variant in Pedigree 2 (NM_001004334.3; c.903+343G>A). The identified variants segregated with the phenotype in all pedigrees. All deep intronic variants were predicted to generate a splice acceptor gain causing aberrant exonization in NR2E3 [89 base pairs (bp)], GPR179 (197 bp), and CNGB3 (73 bp); splicing defects were validated through patient-derived cDNA experiments and/or minigene constructs and rescued by antisense oligonucleotide treatment.

Conclusions

Deep intronic mutations contribute to missing heritability in retinal dystrophies. Combining results from phenotype-directed gene panel testing, GS, and in silico splice prediction tools can help identify these difficult-to-detect pathogenic deep intronic variants.