Photoreceptor Progenitor mRNA Analysis Reveals Exon Skipping Resulting from the ABCA4 c.5461-10T→C Mutation in Stargardt Disease

Understanding and Rescuing the Splicing Defect Caused by the Frequent ABCA4 Variant c.4253+43G>A Underlying Stargardt Disease

Pathogenic variants in ABCA4 are the underlying molecular cause of Stargardt disease (STGD1), an autosomal recessive macular dystrophy characterized by a progressive loss of central vision. Among intronic ABCA4 variants, c.4253+43G>A is frequently detected in STGD1 cases and is classified as a hypomorphic allele, generally associated with late-onset cases. This variant was previously reported to alter splicing regulatory sequences, but the splicing outcome is not fully understood yet. In this study, we attempted to better understand its effect on splicing and to rescue the aberrant splicing via antisense oligonucleotides (AONs). Wild-type and c.4253+43G>A variant-harboring maxigene vectors revealed additional skipping events, which were not previously detected upon transfection in HEK293T cells. To restore exon inclusion, we designed a set of 27 AONs targeting either splicing silencer motifs or the variant region and screened these in maxigene-transfected HEK293T cells. Candidate AONs able to promote exon inclusion were selected for further testing in patient-derived photoreceptor precursor cells. Surprisingly, no robust splicing modulation was observed in this model system. Overall, this research helped to adequately characterize the splicing alteration caused by the c.4253+43G>A variant, although future development of AON-mediated exon inclusion therapy for ABCA4 is needed.

Therapeutic strategies for aberrant splicing in cancer and genetic disorders

Accurate pre-mRNA splicing is essential for proper protein translation; however, aberrant splicing is commonly observed in the context of cancer and genetic disorders. Notably, in genetic diseases, these splicing abnormalities often play a pivotal role. Substantial challenges persist in accurately identifying and classifying disease-induced aberrant splicing, as well as in development of targeted therapeutic strategies. In this review, we examine prevalent forms of aberrant splicing and explore potential therapeutic approaches aimed at addressing these splicing-related diseases. This summary contributes to a deeper understanding of the complexities about aberrant splicing and provide a foundation for the development of effective therapeutic interventions in the field of genetic disorders and cancer.

Towards Uncovering the Role of Incomplete Penetrance in Maculopathies through Sequencing of 105 Disease-Associated Genes

Inherited macular dystrophies (iMDs) are a group of genetic disorders, which affect the central region of the retina. To investigate the genetic basis of iMDs, we used single-molecule Molecular Inversion Probes to sequence 105 maculopathy-associated genes in 1352 patients diagnosed with iMDs. Within this cohort, 39.8% of patients were considered genetically explained by 460 different variants in 49 distinct genes of which 73 were novel variants, with some affecting splicing. The top five most frequent causative genes were ABCA4 (37.2%), PRPH2 (6.7%), CDHR1 (6.1%), PROM1 (4.3%) and RP1L1 (3.1%). Interestingly, variants with incomplete penetrance were revealed in almost one-third of patients considered solved (28.1%), and therefore, a proportion of patients may not be explained solely by the variants reported. This includes eight previously reported variants with incomplete penetrance in addition to CDHR1:c.783G>A and CNGB3:c.1208G>A. Notably, segregation analysis was not routinely performed for variant phasing-a limitation, which may also impact the overall diagnostic yield. The relatively high proportion of probands without any putative causal variant (60.2%) highlights the need to explore variants with incomplete penetrance, the potential modifiers of disease and the genetic overlap between iMDs and age-related macular degeneration. Our results provide valuable insights into the genetic landscape of iMDs and warrant future exploration to determine the involvement of other maculopathy genes.

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.

QR-1011 restores defective ABCA4 splicing caused by multiple severe ABCA4 variants underlying Stargardt disease

Stargardt disease type 1 (STGD1), the most common form of hereditary macular dystrophy, can be caused by biallelic combinations of over 2200 variants in the ABCA4 gene. This leads to reduced or absent ABCA4 protein activity, resulting in toxic metabolite accumulation in the retina and damage of the retinal pigment epithelium and photoreceptors. Approximately 21% of all ABCA4 variants that contribute to disease influence ABCA4 pre-mRNA splicing. This emphasizes the need for therapies to restore disrupted ABCA4 splicing and halt STGD1 progression. Previously, QR-1011, an antisense oligonucleotide (AON), successfully corrected splicing abnormalities and restored normal ABCA4 protein translation in human retinal organoids carrying the prevalent disease-causing variant c.5461-10T>C in ABCA4. Here, we investigated whether QR-1011 could also correct splicing in four less common non-canonical splice site (NCSS) variants flanking ABCA4 exon 39: c.5461-8T>G, c.5461-6T>C, c.5584+5G>A and c.5584+6T>C. We administered QR-1011 and three other AONs to midigene-transfected cells and demonstrate that QR-1011 had the most pronounced effect on splicing compared to the others. Moreover, QR-1011 significantly increased full-length ABCA4 transcript levels for c.5461-8T>G and c.5584+6T>C. Splicing restoration could not be achieved in the other two variants, suggesting their more severe effect on splicing. Overall, QR-1011, initially developed for a single ABCA4 variant, exhibited potent splice correction capabilities for two additional severe NCSS variants nearby. This suggests the possibility of a broader therapeutic impact of QR-1011 extending beyond its original target and highlights the potential for treating a larger population of STGD1 patients affected by multiple severe ABCA4 variants with a single AON.

Stargardt disease-associated missense and synonymous ABCA4 variants result in aberrant splicing

Missense variants in ABCA4 constitute ~50% of causal variants in Stargardt disease (STGD1). Their pathogenicity is attributed to their direct effect on protein function, whilst their potential impact on pre-mRNA splicing disruption remains poorly understood. Interestingly, synonymous ABCA4 variants have previously been classified as 'severe' variants based on in silico analyses. Here, we systemically investigated the role of synonymous and missense variants in ABCA4 splicing by combining computational predictions and experimental assays. To identify variants of interest, we used SpliceAI to ascribe defective splice predictions on a dataset of 5579 biallelic STGD1 probands. We selected those variants with predicted delta scores for acceptor/donor gain > 0.20, and no previous reports on their effect on splicing. Fifteen ABCA4 variants were selected, 4 of which were predicted to create a new splice acceptor site and 11 to create a new splice donor site. In addition, three variants of interest with delta scores < 0.20 were included. The variants were introduced in wild-type midigenes that contained 4-12 kb of ABCA4 genomic sequence, which were subsequently expressed in HEK293T cells. By using RT-PCR and Sanger sequencing, we identified splice aberrations for 16 of 18 analyzed variants. SpliceAI correctly predicted the outcomes for 15 out of 18 variants, illustrating its reliability in predicting the impact of coding ABCA4 variants on splicing. Our findings highlight a causal role for coding ABCA4 variants in splicing aberrations, improving the severity assessment of missense and synonymous ABCA4 variants, and guiding to new treatment strategies for STGD1.

Targeted sequencing and in vitro splice assays shed light on ABCA4-associated retinopathies missing heritability

The ABCA4 gene is the most frequently mutated Mendelian retinopathy-associated gene. Biallelic variants lead to a variety of phenotypes, however, for thousands of cases the underlying variants remain unknown. Here, we aim to shed further light on the missing heritability of ABCA4-associated retinopathy by analyzing a large cohort of macular dystrophy probands. A total of 858 probands were collected from 26 centers, of whom 722 carried no or one pathogenic ABCA4 variant, while 136 cases carried two ABCA4 alleles, one of which was a frequent mild variant, suggesting that deep-intronic variants (DIVs) or other cis-modifiers might have been missed. After single molecule molecular inversion probes (smMIPs)-based sequencing of the complete 128-kb ABCA4 locus, the effect of putative splice variants was assessed in vitro by midigene splice assays in HEK293T cells. The breakpoints of copy number variants (CNVs) were determined by junction PCR and Sanger sequencing. ABCA4 sequence analysis solved 207 of 520 (39.8%) naive or unsolved cases and 70 of 202 (34.7%) monoallelic cases, while additional causal variants were identified in 54 of 136 (39.7%) probands carrying two variants. Seven novel DIVs and six novel non-canonical splice site variants were detected in a total of 35 alleles and characterized, including the c.6283-321C>G variant leading to a complex splicing defect. Additionally, four novel CNVs were identified and characterized in five alleles. These results confirm that smMIPs-based sequencing of the complete ABCA4 gene provides a cost-effective method to genetically solve retinopathy cases and that several rare structural and splice altering defects remain undiscovered in Stargardt disease cases.

ABCA4 Variant c.5714+5G>A in Trans With Null Alleles Results in Primary RPE Damage

Purpose

To determine the disease pathogenesis associated with the frequent ABCA4 variant c.5714+5G>A (p.[=,Glu1863Leufs*33]).

Methods

Patient-derived photoreceptor precursor cells were generated to analyze the effect of c.5714+5G>A on splicing and perform a quantitative analysis of c.5714+5G>A products. Patients with c.5714+5G>A in trans with a null allele (i.e., c.5714+5G>A patients; n = 7) were compared with patients with two null alleles (i.e., double null patients; n = 11); with a special attention to the degree of RPE atrophy (area of definitely decreased autofluorescence and the degree of photoreceptor impairment (outer nuclear layer thickness and pattern electroretinography amplitude).

Results

RT-PCR of mRNA from patient-derived photoreceptor precursor cells showed exon 40 and exon 39/40 deletion products, as well as the normal transcript. Quantification of products showed 52.4% normal and 47.6% mutant ABCA4 mRNA. Clinically, c.5714+5G>A patients displayed significantly better structural and functional preservation of photoreceptors (thicker outer nuclear layer, presence of tubulations, higher pattern electroretinography amplitude) than double null patients with similar degrees of RPE loss, whereas double null patients exhibited signs of extensive photoreceptor ,damage even in the areas with preserved RPE.

Conclusions

The prototypical STGD1 sequence of events of primary RPE and secondary photoreceptor damage is congruous with c.5714+5G>A, but not the double null genotype, which implies different and genotype-dependent disease mechanisms. We hypothesize that the relative photoreceptor sparing in c.5714+5G>A patients results from the remaining function of the ABCA4 transporter originating from the normally spliced product, possibly by decreasing the direct bisretinoid toxicity on photoreceptor membranes.

Stargardt disease-associated in-frame ABCA4 exon 17 skipping results in significant ABCA4 function

BACKGROUND

ABCA4, the gene implicated in Stargardt disease (STGD1), contains 50 exons, of which 17 contain multiples of three nucleotides. The impact of in-frame exon skipping is yet to be determined. Antisense oligonucleotides (AONs) have been investigated in Usher syndrome-associated genes to induce skipping of in-frame exons carrying severe variants and mitigate their disease-linked effect. Upon the identification of a STGD1 proband carrying a novel exon 17 canonical splice site variant, the activity of ABCA4 lacking 22 amino acids encoded by exon 17 was examined, followed by design of AONs able to induce exon 17 skipping.

METHODS

A STGD1 proband was compound heterozygous for the splice variant c.2653+1G>A, that was predicted to result in in-frame skipping of exon 17, and a null variant [c.735T>G, p.(Tyr245*)]. Clinical characteristics of this proband were studied using multi-modal imaging and complete ophthalmological examination. The aberrant splicing of c.2653+1G>A was investigated in vitro in HEK293T cells with wild-type and mutant midigenes. The residual activity of the mutant ABCA4 protein lacking Asp864-Gly885 encoded by exon 17 was analyzed with all-trans-retinal-activated ATPase activity assay, along with its subcellular localization. To induce exon 17 skipping, the effect of 40 AONs was examined in vitro in WT WERI-Rb-1 cells and 3D human retinal organoids.

RESULTS

Late onset STGD1 in the proband suggests that c.2653+1G>A does not have a fully deleterious effect. The in vitro splice assay confirmed that this variant leads to ABCA4 transcripts without exon 17. ABCA4 Asp864_Gly863del was stable and retained 58% all-trans-retinal-activated ATPase activity compared to WT ABCA4. This sequence is located in an unstructured linker region between transmembrane domain 6 and nucleotide-binding domain-1 of ABCA4. AONs were designed to possibly reduce pathogenicity of severe variants harbored in exon 17. The best AON achieved 59% of exon 17 skipping in retinal organoids.

CONCLUSIONS

Exon 17 deletion in ABCA4 does not result in the absence of protein activity and does not cause a severe STGD1 phenotype when in trans with a null allele. By applying AONs, the effect of severe variants in exon 17 can potentially be ameliorated by exon skipping, thus generating partial ABCA4 activity in STGD1 patients.

Retinal organoids provide unique insights into molecular signatures of inherited retinal disease throughout retinogenesis

The demand for induced pluripotent stem cells (iPSC)-derived retinal organoid and retinal pigment epithelium (RPE) models for the modelling of inherited retinopathies has increased significantly in the last decade. These models are comparable with foetal retinas up until the later stages of retinogenesis, expressing all of the key neuronal markers necessary for retinal function. These models have proven to be invaluable in the understanding of retinogenesis, particular in the context of patient-specific diseases. Inherited retinopathies are infamously described as clinically and phenotypically heterogeneous, such that developing gene/mutation-specific animal models in each instance of retinal disease is not financially or ethically feasible. Further to this, many animal models are insufficient in the study of disease pathogenesis due to anatomical differences and failure to recapitulate human disease phenotypes. In contrast, iPSC-derived retinal models provide a high throughput platform which is physiologically relevant for studying human health and disease. They also serve as a platform for drug screening, gene therapy approaches and in vitro toxicology of novel therapeutics in pre-clinical studies. One unique characteristic of stem cell-derived retinal models is the ability to mimic in vivo retinogenesis, providing unparalleled insights into the effects of pathogenic mutations in cells of the developing retina, in a highly accessible way. This review aims to give the reader an overview of iPSC-derived retinal organoids and/or RPE in the context of disease modelling of several inherited retinopathies including Retinitis Pigmentosa, Stargardt disease and Retinoblastoma. We describe the ability of each model to recapitulate in vivo disease phenotypes, validate previous findings from animal models and identify novel pathomechanisms that underpin individual IRDs.

Combined Single Gene Testing and Genome Sequencing as an Effective Diagnostic Approach for Anophthalmia and Microphthalmia Patients

Anophthalmia and microphthalmia (A/M) are among the most severe congenital developmental eye disorders. Despite the advancements in genome screening technologies, more than half of A/M patients do not receive a molecular diagnosis. We included seven consanguineous families affected with A/M from Pakistani cohort and an unknown molecular basis. Single gene testing of FOXE3 was performed, followed by genome sequencing for unsolved probands in order to establish a genetic diagnosis for these families. All seven families were provided with a genetic diagnosis. The identified variants were all homozygous, classified as (likely) pathogenic and present in an A/M-associated gene. Targeted FOXE3 sequencing revealed two previously reported pathogenic FOXE3 variants in four families. In the remaining families, genome sequencing revealed a known pathogenic PXDN variant, a novel 13bp deletion in VSX2, and one novel deep intronic splice variant in PXDN. An in vitro splice assay was performed for the PXDN splice variant which revealed a severe splicing defect. Our study confirmed the utility of genome sequencing as a diagnostic tool for A/M-affected individuals. Furthermore, the identification of a novel deep intronic pathogenic variant in PXDN highlights the role of non-coding variants in A/M-disorders and the value of genome sequencing for the identification of this type of variants.

Detailed analysis of an enriched deep intronic ABCA4 variant in Irish Stargardt disease patients

Over 15% of probands in a large cohort of more than 1500 inherited retinal degeneration patients present with a clinical diagnosis of Stargardt disease (STGD1), a recessive form of macular dystrophy caused by biallelic variants in the ABCA4 gene. Participants were clinically examined and underwent either target capture sequencing of the exons and some pathogenic intronic regions of ABCA4, sequencing of the entire ABCA4 gene or whole genome sequencing. ABCA4 c.4539 + 2028C > T, p.[= ,Arg1514Leufs*36] is a pathogenic deep intronic variant that results in a retina-specific 345-nucleotide pseudoexon inclusion. Through analysis of the Irish STGD1 cohort, 25 individuals across 18 pedigrees harbour ABCA4 c.4539 + 2028C > T and another pathogenic variant. This includes, to the best of our knowledge, the only two homozygous patients identified to date. This provides important evidence of variant pathogenicity for this deep intronic variant, highlighting the value of homozygotes for variant interpretation. 15 other heterozygous incidents of this variant in patients have been reported globally, indicating significant enrichment in the Irish population. We provide detailed genetic and clinical characterization of these patients, illustrating that ABCA4 c.4539 + 2028C > T is a variant of mild to intermediate severity. These results have important implications for unresolved STGD1 patients globally with approximately 10% of the population in some western countries claiming Irish heritage. This study exemplifies that detection and characterization of founder variants is a diagnostic imperative.

A Novel Intronic Deletion in PDE6B Causes Autosomal Recessive Retinitis Pigmentosa by Interfering with RNA Splicing

INTRODUCTION

Retinitis pigmentosa (RP) is a rare degenerative retinal disease caused by mutations in approximately seventy genes. Currently, despite the availability of large-scale DNA sequencing technologies, ∼30-40% of patients still cannot be diagnosed at the molecular level. In this study, we investigated a novel intronic deletion of PDE6B, encoding the beta subunit of phosphodiesterase 6 in association with recessive RP.

METHODS

Three unrelated consanguineous families were recruited from the northwestern part of Pakistan. Whole exome sequencing was performed for the proband of each family, and the data were analyzed according to an in-house computer pipeline. Relevant DNA variants in all available members of these families were assessed through Sanger sequencing. A minigene-based splicing assay was also performed.

RESULTS

The clinical phenotype for all patients was compatible with rod cone degeneration, with the onset during childhood. Whole exome sequencing revealed a homozygous 18 bp intronic deletion (NM_000283.3:c.1921-20_1921-3del) in PDE6B, which co-segregated with disease in 10 affected individuals. In vitro splicing tests showed that this deletion causes aberrant RNA splicing of the gene, leading to the in-frame deletion of 6 codons and, likely, to disease.

CONCLUSION

Our findings further expand the mutational spectrum of the PDE6B gene.

Antisense oligonucleotide therapy corrects splicing in the common Stargardt disease type 1-causing variant ABCA4 c.5461-10T>C

Stargardt disease type 1 (STGD1) is the most common hereditary form of maculopathy and remains untreatable. STGD1 is caused by biallelic variants in the ABCA4 gene, which encodes the ATP-binding cassette (type 4) protein (ABCA4) that clears toxic byproducts of the visual cycle. The c.5461-10T>C p.[Thr1821Aspfs∗6,Thr1821Valfs∗13] variant is the most common severe disease-associated variant, and leads to exon skipping and out-of-frame ABCA4 transcripts that prevent translation of functional ABCA4 protein. Homozygous individuals typically display early onset STGD1 and are legally blind by early adulthood. Here, we applied antisense oligonucleotides (AONs) to promote exon inclusion and restore wild-type RNA splicing of ABCA4 c.5461-10T>C. The effect of AONs was first investigated in vitro using an ABCA4 midigene model. Subsequently, the best performing AONs were administered to homozygous c.5461-10T>C 3D human retinal organoids. Isoform-specific digital polymerase chain reaction revealed a significant increase in correctly spliced transcripts after treatment with the lead AON, QR-1011, up to 53% correct transcripts at a 3 μM dose. Furthermore, western blot and immunohistochemistry analyses identified restoration of ABCA4 protein after treatment. Collectively, we identified QR-1011 as a potent splice-correcting AON and a possible therapeutic intervention for patients harboring the severe ABCA4 c.5461-10T>C variant.

Long-Read Sequencing Identifies Novel Pathogenic Intronic Variants in Gitelman Syndrome

BACKGROUND

Gitelman syndrome is a salt-losing tubulopathy characterized by hypokalemic alkalosis and hypomagnesemia. It is caused by homozygous recessive or compound heterozygous pathogenic variants in SLC12A3 , which encodes the Na + -Cl - cotransporter (NCC). In up to 10% of patients with Gitelman syndrome, current genetic techniques detect only one specific pathogenic variant. This study aimed to identify a second pathogenic variant in introns, splice sites, or promoters to increase the diagnostic yield.

METHODS

Long-read sequencing of SLC12A3 was performed in 67 DNA samples from individuals with suspected Gitelman syndrome in whom a single likely pathogenic or pathogenic variant was previously detected. In addition, we sequenced DNA samples from 28 individuals with one variant of uncertain significance or no candidate variant. Midigene splice assays assessed the pathogenicity of novel intronic variants.

RESULTS

A second likely pathogenic/pathogenic variant was identified in 45 (67%) patients. Those with two likely pathogenic/pathogenic variants had a more severe electrolyte phenotype than other patients. Of the 45 patients, 16 had intronic variants outside of canonic splice sites (nine variants, mostly deep intronic, six novel), whereas 29 patients had an exonic variant or canonic splice site variant. Midigene splice assays of the previously known c.1670-191C>T variant and intronic candidate variants demonstrated aberrant splicing patterns.

CONCLUSION

Intronic pathogenic variants explain an important part of the missing heritability in Gitelman syndrome. Long-read sequencing should be considered in diagnostic workflows for Gitelman syndrome.

Whole genome sequencing for USH2A-associated disease reveals several pathogenic deep-intronic variants that are amenable to splice correction

A significant number of individuals with a rare disorder such as Usher syndrome (USH) and (non-)syndromic autosomal recessive retinitis pigmentosa (arRP) remain genetically unexplained. Therefore, we assessed subjects suspected of USH2A-associated disease and no or mono-allelic USH2A variants using whole genome sequencing (WGS) followed by an improved pipeline for variant interpretation to provide a conclusive diagnosis. One hundred subjects were screened using WGS to identify causative variants in USH2A or other USH/arRP-associated genes. In addition to the existing variant interpretation pipeline, a particular focus was put on assessing splice-affecting properties of variants, both in silico and in vitro. Also structural variants were extensively addressed. For variants resulting in pseudoexon inclusion, we designed and evaluated antisense oligonucleotides (AONs) using minigene splice assays and patient-derived photoreceptor precursor cells. Biallelic variants were identified in 49 of 100 subjects, including novel splice-affecting variants and structural variants, in USH2A or arRP/USH-associated genes. Thirteen variants were shown to affect USH2A pre-mRNA splicing, including four deep-intronic USH2A variants resulting in pseudoexon inclusion, which could be corrected upon AON treatment. We have shown that WGS, combined with a thorough variant interpretation pipeline focused on assessing pre-mRNA splicing defects and structural variants, is a powerful method to provide subjects with a rare genetic condition, a (likely) conclusive genetic diagnosis. This is essential for the development of future personalized treatments and for patients to be eligible for such treatments.

Minigene Splicing Assays and Long-Read Sequencing to Unravel Pathogenic Deep-Intronic Variants in PAX6 in Congenital Aniridia

PAX6 haploinsufficiency causes aniridia, a congenital eye disorder that involves the iris, and foveal hypoplasia. Comprehensive screening of the PAX6 locus, including the non-coding regions, by next-generation sequencing revealed four deep-intronic variants with potential effects on pre-RNA splicing. Nevertheless, without a functional analysis, their pathogenicity could not be established. We aimed to decipher their impact on the canonical PAX6 splicing using in vitro minigene splicing assays and nanopore-based long-read sequencing. Two multi-exonic PAX6 constructs were generated, and minigene assays were carried out. An aberrant splicing pattern was observed for two variants in intron 6, c.357+136G>A and c.357+334G>A. In both cases, several exonization events, such as pseudoexon inclusions and partial intronic retention, were observed due to the creation or activation of new/cryptic non-canonical splicing sites, including a shared intronic donor site. In contrast, two variants identified in intron 11, c.1032+170A>T and c.1033-275A>C, seemed not to affect splicing processes. We confirmed the high complexity of alternative splicing of PAX6 exon 6, which also involves unreported cryptic intronic sites. Our study highlights the importance of integrating functional studies into diagnostic algorithms to decipher the potential implication of non-coding variants, usually classified as variants of unknown significance, thus allowing variant reclassification to achieve a conclusive genetic diagnosis.

Correction of the Splicing Defect Caused by a Recurrent Variant in ABCA4 (c.769-784C>T) That Underlies Stargardt Disease

Stargardt disease is an inherited retinal disease caused by biallelic mutations in the ABCA4 gene, many of which affect ABCA4 splicing. In this study, nine antisense oligonucleotides (AONs) were designed to correct pseudoexon (PE) inclusion caused by a recurrent deep-intronic variant in ABCA4 (c.769-784C>T). First, the ability of AONs to skip the PE from the final ABCA4 mRNA transcript was assessed in two cellular models carrying the c.769-784C>T variant: a midigene assay using HEK293T cells and patient-derived fibroblasts. Based on the splicing-correcting ability of each individual AON, the three most efficacious AONs targeting independent regions of the PE were selected for a final assessment in photoreceptor precursor cells (PPCs). The final analysis in the PPC model confirmed high efficacy of AON2, -5, and -7 in promoting PE exclusion. Among the three AONs, AON2 is chosen as the lead candidate for further optimization, hereby showcasing the high potential of AONs to correct aberrant splicing events driven by deep-intronic variants.

Characterising splicing defects of ABCA4 variants within exons 13-50 in patient-derived fibroblasts

The ATP-binding cassette subfamily A member 4 gene (ABCA4)-associated retinopathy, Stargardt disease, is the most common monogenic inherited retinal disease. Given the pathogenicity of numerous ABCA4 variants is yet to be examined and a significant proportion (more than 15%) of ABCA4 variants are categorized as splice variants in silico, we therefore established a fibroblast-based splice assay to analyze ABCA4 variants in an Australian Stargardt disease cohort and characterize the pathogenic mechanisms of ABCA4 variants. A cohort of 67 patients clinically diagnosed with Stargardt disease was recruited. Genomic DNA was analysed using a commercial panel for ABCA4 variant detection and the consequences of ABCA4 variants were predicted in silico. Dermal fibroblasts were propagated from skin biopsies, total RNA was extracted and the ABCA4 transcript was amplified by RT-PCR. Our analysis identified a total of 67 unique alleles carrying 74 unique variants. The most prevalent splice-affecting complex allele c.[5461-10T>C; 5603A>T] was carried by 10% of patients in a compound heterozygous state. ABCA4 transcripts from exon 13 to exon 50 were readily detected in fibroblasts. In this region, aberrant splicing was evident in 10 out of 57 variant transcripts (18%), carried by 19 patients (28%). Patient-derived fibroblasts provide a feasible platform for identification of ABCA4 splice variants located within exons 13-50. Experimental evidence of aberrant splicing contributes to the pathogenic classification for ABCA4 variants. Moreover, identification of variants that affect splicing processes provides opportunities for intervention, in particular antisense oligonucleotide-mediated splice correction.

Phenotyping of ABCA4 Retinopathy by Machine Learning Analysis of Full-Field Electroretinography

Purpose

Biallelic pathogenic variants in ABCA4 are the commonest cause of monogenic retinal disease. The full-field electroretinogram (ERG) quantifies severity of retinal dysfunction. We explored application of machine learning in ERG interpretation and in genotype-phenotype correlations.

Methods

International standard ERGs in 597 cases of ABCA4 retinopathy were classified into three functional phenotypes by human experts: macular dysfunction alone (group 1), or with additional generalized cone dysfunction (group 2), or both cone and rod dysfunction (group 3). Algorithms were developed for automatic selection and measurement of ERG components and for classification of ERG phenotype. Elastic-net regression was used to quantify severity of specific ABCA4 variants based on effect on retinal function.

Results

Of the cohort, 57.6%, 7.4%, and 35.0% fell into groups 1, 2, and 3 respectively. Compared with human experts, automated classification showed overall accuracy of 91.8% (SE, 0.169), and 96.7%, 39.3%, and 93.8% for groups 1, 2, and 3. When groups 2 and 3 were combined, the average holdout group accuracy was 93.6% (SE, 0.142). A regression model yielded phenotypic severity scores for the 47 commonest ABCA4 variants.

Conclusions

This study quantifies prevalence of phenotypic groups based on retinal function in a uniquely large single-center cohort of patients with electrophysiologically characterized ABCA4 retinopathy and shows applicability of machine learning. Novel regression-based analyses of ABCA4 variant severity could identify individuals predisposed to severe disease.

Translational Relevance

Machine learning can yield meaningful classifications of ERG data, and data-driven scoring of genetic variants can identify patients likely to benefit most from future therapies.

Scrutinizing pathogenicity of the USH2A c.2276 G > T; p.(Cys759Phe) variant

The USH2A variant c.2276 G > T (p.(Cys759Phe)) has been described by many authors as a frequent cause of autosomal recessive retinitis pigmentosa (arRP). However, this is in contrast with the description of two asymptomatic individuals homozygous for this variant. We therefore assessed pathogenicity of the USH2A c.2276 G > T variant using extensive genetic and functional analyses. Whole genome sequencing and optical genome mapping were performed for three arRP cases homozygous for USH2A c.2276 G > T to exclude alternative genetic causes. A minigene splice assay was designed to investigate the effect of c.2276 G > T on pre-mRNA splicing, in presence or absence of the nearby c.2256 T > C variant. Moreover, an ush2a^{p.(Cys771Phe)} zebrafish knock-in model mimicking human p.(Cys759Phe) was generated and characterized using functional and immunohistochemical analyses. Besides the homozygous c.2276 G > T USH2A variant, no alternative genetic causes were identified. Evaluation of the ush2a^{p.(Cys771Phe)} zebrafish model revealed strongly reduced levels of usherin expression at the photoreceptor periciliary membrane, increased levels of rhodopsin localization in the photoreceptor cell body and decreased electroretinogram (ERG) b-wave amplitudes compared to wildtype controls. In conclusion, we confirmed pathogenicity of USH2A c.2276 G > T (p.(Cys759Phe)). Consequently, cases homozygous for c.2276 G > T can now receive a definite genetic diagnosis and can be considered eligible for receiving future QR-421a-mediated exon 13 skipping therapy.

Noncanonical Splice Site and Deep Intronic FRMD7 Variants Activate Cryptic Exons in X-linked Infantile Nystagmus

Purpose

We aim to report noncoding pathogenic variants in patients with FRMD7-related infantile nystagmus (FIN).

Methods

Genome sequencing (n = 2 families) and reanalysis of targeted panel next generation sequencing (n = 2 families) was performed in genetically unsolved cases of suspected FIN. Previous sequence analysis showed no pathogenic coding variants in genes associated with infantile nystagmus. SpliceAI, SpliceRover, and Alamut consensus programs were used to annotate noncoding variants. Minigene splicing assay was performed to confirm aberrant splicing. In silico analysis of exonic splicing enhancer and silencer was also performed.

Results

FRMD7 intronic variants were identified based on genome sequencing and targeted next-generation sequencing analysis. These included c.285-12A>G (pedigree 1), c.284+63T>A (pedigrees 2 and 3), and c. 383-1368A>G (pedigree 4). All variants were absent in gnomAD, and the both c.285-12A>G and c.284+63T>A variants were predicted to enhance new splicing acceptor gains with SpliceAI, SpliceRover, and Alamut consensus approaches. However, the c.383-1368 A>G variant only had a significant impact score on the SpliceRover program. The c.383-1368A>G variant was predicted to promote pseudoexon inclusion by binding of exonic splicing enhancer. Aberrant exonizations were validated through minigene constructs, and all variants were segregated in the families.

Conclusions

Deep learning–based annotation of noncoding variants facilitates the discovery of hidden genetic variations in patients with FIN. This study provides evidence of effectiveness of combined deep learning–based splicing tools to identify hidden pathogenic variants in previously unsolved patients with infantile nystagmus.

Translational Relevance

These results demonstrate robust analysis using two deep learning splicing predictions and in vitro functional study can lead to finding hidden genetic variations in unsolved patients.

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.

Usher syndrome type IV: clinically and molecularly confirmed by novel ARSG variants

Usher syndrome (USH) is an autosomal recessively inherited disease characterized by sensorineural hearing loss (SNHL) and retinitis pigmentosa (RP) with or without vestibular dysfunction. It is highly heterogeneous both clinically and genetically. Recently, variants in the arylsulfatase G (ARSG) gene have been reported to underlie USH type IV. This distinct type of USH is characterized by late-onset RP with predominantly pericentral and macular changes, and late onset SNHL without vestibular dysfunction. In this study, we describe the USH type IV phenotype in three unrelated subjects. We identified three novel pathogenic variants, two novel likely pathogenic variants, and one previously described pathogenic variant in ARSG. Functional experiments indicated a loss of sulfatase activity of the mutant proteins. Our findings confirm that ARSG variants cause the newly defined USH type IV and support the proposed extension of the phenotypic USH classification.

Photoreceptor degeneration in ABCA4-associated retinopathy and its genetic correlates

BACKGROUND

Outcome measures sensitive to disease progression are needed for ATP-binding cassette, sub-family A, member 4–associated (ABCA4-associated) retinopathy. We aimed to quantify ellipsoid zone (EZ) loss and photoreceptor degeneration beyond EZ-loss in ABCA4-associated retinopathy and investigate associations between photoreceptor degeneration, genotype, and age.

METHODS

We analyzed 132 eyes from 66 patients (of 67 enrolled) with molecularly confirmed ABCA4-associated retinopathy from a prospective natural history study with a median [IQR] follow-up of 4.2 years [3.1, 5.1]. Longitudinal spectral-domain optical coherence tomography volume scans (37 B-scans, 30° × 15°) were segmented using a deep learning (DL) approach. For genotype-phenotype analysis, a model of ABCA4 variants was applied with the age of criterion EZ-loss (6.25 mm²) as the dependent variable.

RESULTS

Patients exhibited an average (square-root-transformed) EZ-loss progression rate of [95% CI] 0.09 mm/y [0.06, 0.11]. Outer nuclear layer (ONL) thinning extended beyond the area of EZ-loss. The average distance from the EZ-loss boundary to normalization of ONL thickness (to ±2 z score units) was 3.20° [2.53, 3.87]. Inner segment (IS) and outer segment (OS) thinning was less pronounced, with an average distance from the EZ-loss boundary to layer thickness normalization of 1.20° [0.91, 1.48] for the IS and 0.60° [0.49, 0.72] for the OS. An additive model of allele severity explained 52.7% of variability in the age of criterion EZ-loss.

CONCLUSION

Patients with ABCA4-associated retinopathy exhibited significant alterations of photoreceptors outside of EZ-loss. DL-based analysis of photoreceptor laminae may help monitor disease progression and estimate the severity of ABCA4 variants.

TRIAL REGISTRATION

ClinicalTrials.gov identifier: NCT01736293.

FUNDING

National Eye Institute Intramural Research Program and German Research Foundation grant PF950/1-1.

Structure and function of ABCA4 and its role in the visual cycle and Stargardt macular degeneration

ABCA4 is a member of the superfamily of ATP-binding cassette (ABC) transporters that is preferentially localized along the rim region of rod and cone photoreceptor outer segment disc membranes. It uses the energy from ATP binding and hydrolysis to transport N-retinylidene-phosphatidylethanolamine (N-Ret-PE), the Schiff base adduct of retinal and phosphatidylethanolamine, from the lumen to the cytoplasmic leaflet of disc membranes. This ensures that all-trans-retinal and excess 11-cis-retinal are efficiently cleared from photoreceptor cells thereby preventing the accumulation of toxic retinoid compounds. Loss-of-function mutations in the gene encoding ABCA4 cause autosomal recessive Stargardt macular degeneration, also known as Stargardt disease (STGD1), and related autosomal recessive retinopathies characterized by impaired central vision and an accumulation of lipofuscin and bis-retinoid compounds. High resolution structures of ABCA4 in its substrate and nucleotide free state and containing bound N-Ret-PE or ATP have been determined by cryo-electron microscopy providing insight into the molecular architecture of ABCA4 and mechanisms underlying substrate recognition and conformational changes induced by ATP binding. The expression and functional characterization of a large number of disease-causing missense ABCA4 variants have been determined. These studies have shed light into the molecular mechanisms underlying Stargardt disease and a classification that reliably predicts the effect of a specific missense mutation on the severity of the disease. They also provide a framework for developing rational therapeutic treatments for ABCA4-associated diseases.

Stargardt disease and progress in therapeutic strategies

Background: Stargardt disease (STGD1) is an autosomal recessive retinal dystrophy due to mutations in ABCA4, characterized by subretinal deposition of lipofuscin-like substances and bilateral centrifugal vision loss. Despite the tremendous progress made in the understanding of STGD1, there are no approved treatments to date. This review examines the challenges in the development of an effective STGD1 therapy.Materials and Methods: A literature review was performed through to June 2021 summarizing the spectrum of retinal phenotypes in STGD1, the molecular biology of ABCA4 protein, the in vivo and in vitro models used to investigate the mechanisms of ABCA4 mutations and current clinical trials.Results: STGD1 phenotypic variability remains an challenge for clinical trial design and patient selection. Pre-clinical development of therapeutic options has been limited by the lack of animal models reflecting the diverse phenotypic spectrum of STDG1. Patient-derived cell lines have facilitated the characterization of splice mutations but the clinical presentation is not always predicted by the effect of specific mutations on retinoid metabolism in cellular models. Current therapies primarily aim to delay vision loss whilst strategies to restore vision are less well developed.Conclusions: STGD1 therapy development can be accelerated by a deeper understanding of genotype-phenotype correlations.

An Overview of the Genetics of ABCA4 Retinopathies, an Evolving Story

Stargardt disease (STGD1) and ABCA4 retinopathies (ABCA4R) are caused by pathogenic variants in the ABCA4 gene inherited in an autosomal recessive manner. The gene encodes an importer flippase protein that prevents the build-up of vitamin A derivatives that are toxic to the RPE. Diagnosing ABCA4R is complex due to its phenotypic variability and the presence of other inherited retinal dystrophy phenocopies. ABCA4 is a large gene, comprising 50 exons; to date > 2000 variants have been described. These include missense, nonsense, splicing, structural, and deep intronic variants. Missense variants account for the majority of variants in ABCA4. However, in a significant proportion of patients with an ABCA4R phenotype, a second variant in ABCA4 is not identified. This could be due to the presence of yet unknown variants, or hypomorphic alleles being incorrectly classified as benign, or the possibility that the disease is caused by a variant in another gene. This underlines the importance of accurate genetic testing. The pathogenicity of novel variants can be predicted using in silico programs, but these rely on databases that are not ethnically diverse, thus highlighting the need for studies in differing populations. Functional studies in vitro are useful towards assessing protein function but do not directly measure the flippase activity. Obtaining an accurate molecular diagnosis is becoming increasingly more important as targeted therapeutic options become available; these include pharmacological, gene-based, and cell replacement-based therapies. The aim of this review is to provide an update on the current status of genotyping in ABCA4 and the status of the therapeutic approaches being investigated.

Analysis of Pathogenic Variants Correctable With CRISPR Base Editing Among Patients With Recessive Inherited Retinal Degeneration

Importance

Many common inherited retinal diseases are not easily treated with gene therapy. Gene editing with base editors may allow the targeted repair of single-nucleotide transition variants in DNA and RNA. It is unknown how many patients have pathogenic variants that are correctable with a base editing strategy.

Objective

To assess the prevalence and spectrum of pathogenic single-nucleotide variants amenable to base editing in common large recessively inherited genes that are associated with inherited retinal degeneration.

Design, Setting, and Participants

In this retrospective cross-sectional study, nonidentifiable records of patients with biallelic pathogenic variants of genes associated with inherited retinal degeneration between July 2013 and December 2019 were analyzed using data from the Oxford University Hospitals Medical Genetics Laboratories, the Leiden Open Variation Database, and previously published studies. Six candidate genes (ABCA4, CDH23, CEP290, EYS, MYO7A, and USH2A), which were determined to be the most common recessive genes with coding sequences not deliverable in a single adeno-associated viral vector, were examined. Data were analyzed from April 16 to May 11, 2020.

Main Outcomes and Measures

Proportion of alleles with a pathogenic transition variant that is potentially correctable with a base editing strategy and proportion of patients with a base-editable allele.

Results

A total of 12 369 alleles from the Leiden Open Variation Database and 179 patients who received diagnoses through the genetic service of the Oxford University Hospitals Medical Genetics Laboratories were analyzed. Editable variants accounted for 53% of all pathogenic variants in the candidate genes contained in the Leiden Open Variation Database. The proportion of pathogenic alleles that were editable varied by gene; 63.1% of alleles in ABCA4, 62.7% of alleles in CDH23, 53.8% of alleles in MYO7A, 41.6% of alleles in CEP290, 37.3% of alleles in USH2A, and 22.2% of alleles in EYS were editable. The 5 most common editable pathogenic variants of each gene accounted for a mean (SD) of 19.1% (9.5%) of all pathogenic alleles within each gene. In the Oxford cohort, 136 of 179 patients (76.0%) had at least 1 editable allele. A total of 53 of 107 patients (49.5%) with biallelic pathogenic variants in the gene ABCA4 and 16 of 56 patients (28.6%) with biallelic pathogenic variants in the gene USH2A had 1 of the 5 most common editable alleles.

Conclusions and Relevance

This study found that pathogenic variants amenable to base editing commonly occur in inherited retinal degeneration. These findings, if generalized to other cohorts, provide an approach for developing base editing therapies to treat retinal degeneration not amenable to gene therapy.

Molecular Inversion Probe-Based Sequencing of USH2A Exons and Splice Sites as a Cost-Effective Screening Tool in USH2 and arRP Cases

A substantial proportion of subjects with autosomal recessive retinitis pigmentosa (arRP) or Usher syndrome type II (USH2) lacks a genetic diagnosis due to incomplete USH2A screening in the early days of genetic testing. These cases lack eligibility for optimal genetic counseling and future therapy. USH2A defects are the most frequent cause of USH2 and are also causative in individuals with arRP. Therefore, USH2A is an important target for genetic screening. The aim of this study was to assess unscreened or incompletely screened and unexplained USH2 and arRP cases for (likely) pathogenic USH2A variants. Molecular inversion probe (MIP)-based sequencing was performed for the USH2A exons and their flanking regions, as well as published deep-intronic variants. This was done to identify single nucleotide variants (SNVs) and copy number variants (CNVs) in 29 unscreened or partially pre-screened USH2 and 11 partially pre-screened arRP subjects. In 29 out of these 40 cases, two (likely) pathogenic variants were successfully identified. Four of the identified SNVs and one CNV were novel. One previously identified synonymous variant was demonstrated to affect pre-mRNA splicing. In conclusion, genetic diagnoses were obtained for a majority of cases, which confirms that MIP-based sequencing is an effective screening tool for USH2A. Seven unexplained cases were selected for future analysis with whole genome sequencing.

Cis-acting modifiers in the ABCA4 locus contribute to the penetrance of the major disease-causing variant in Stargardt disease

Over 1200 variants in the ABCA4 gene cause a wide variety of retinal disease phenotypes, the best known of which is autosomal recessive Stargardt disease (STGD1). Disease-causing variation encompasses all mutation categories, from large copy number variants to very mild, hypomorphic missense variants. The most prevalent disease-causing ABCA4 variant, present in ~ 20% of cases of European descent, c.5882G > A p.(Gly1961Glu), has been a subject of controversy since its minor allele frequency (MAF) is as high as ~ 0.1 in certain populations, questioning its pathogenicity, especially in homozygous individuals. We sequenced the entire ~140Kb ABCA4 genomic locus in an extensive cohort of 644 bi-allelic, i.e. genetically confirmed, patients with ABCA4 disease and analyzed all variants in 140 compound heterozygous and 10 homozygous cases for the p.(Gly1961Glu) variant. A total of 23 patients in this cohort additionally harbored the deep intronic c.769-784C > T variant on the p.(Gly1961Glu) allele, which appears on a specific haplotype in ~ 15% of p.(Gly1961Glu) alleles. This haplotype was present in 5/7 of homozygous cases, where the p.(Gly1961Glu) was the only known pathogenic variant. Three cases had an exonic variant on the same allele with the p.(Gly1961Glu). Patients with the c.[769-784C > T;5882G > A] complex allele exhibit a more severe clinical phenotype, as seen in compound heterozygotes with some more frequent ABCA4 mutations, e.g. p.(Pro1380Leu). Our findings indicate that the c.769-784C > T variant is major cis-acting modifier of the p.(Gly1961Glu) allele. The absence of such additional allelic variation on most p.(Gly1961Glu) alleles largely explains the observed paucity of affected homozygotes in the population.

Antisense oligonucleotide-based treatment of retinitis pigmentosa caused by USH2A exon 13 mutations

Mutations in USH2A are among the most common causes of syndromic and non-syndromic retinitis pigmentosa (RP). The two most recurrent mutations in USH2A, c.2299delG and c.2276G > T, both reside in exon 13. Skipping exon 13 from the USH2A transcript presents a potential treatment modality in which the resulting transcript is predicted to encode a slightly shortened usherin protein. Morpholino-induced skipping of ush2a exon 13 in zebrafish ush2a^rmc1 mutants resulted in the production of usherinΔexon 13 protein and a completely restored retinal function. Antisense oligonucleotides were investigated for their potential to selectively induce human USH2A exon 13 skipping. Lead candidate QR-421a induced a concentration-dependent exon 13 skipping in induced pluripotent stem cell (iPSC)-derived photoreceptor precursors from an Usher syndrome patient homozygous for the c.2299delG mutation. Mouse surrogate mQR-421a reached the retinal outer nuclear layer after a single intravitreal injection and induced a detectable level of exon skipping until at least 6 months post-injection. In conclusion, QR-421a-induced exon skipping proves to be a highly promising treatment option for RP caused by mutations in USH2A exon 13.

The Alter Retina: Alternative Splicing of Retinal Genes in Health and Disease

Alternative splicing of mRNA is an essential mechanism to regulate and increase the diversity of the transcriptome and proteome. Alternative splicing frequently occurs in a tissue- or time-specific manner, contributing to differential gene expression between cell types during development. Neural tissues present extremely complex splicing programs and display the highest number of alternative splicing events. As an extension of the central nervous system, the retina constitutes an excellent system to illustrate the high diversity of neural transcripts. The retina expresses retinal specific splicing factors and produces a large number of alternative transcripts, including exclusive tissue-specific exons, which require an exquisite regulation. In fact, a current challenge in the genetic diagnosis of inherited retinal diseases stems from the lack of information regarding alternative splicing of retinal genes, as a considerable percentage of mutations alter splicing or the relative production of alternative transcripts. Modulation of alternative splicing in the retina is also instrumental in the design of novel therapeutic approaches for retinal dystrophies, since it enables precision medicine for specific mutations.

The role of multimodal imaging and vision function testing in ABCA4-related retinopathies and their relevance to future therapeutic interventions

The aim of this review article is to describe the specific features of Stargardt disease and ABCA4 retinopathies (ABCA4R) using multimodal imaging and functional testing and to highlight their relevance to potential therapeutic interventions. Standardised measures of tissue loss, tissue function and rate of change over time using formal structured deep phenotyping in Stargardt disease and ABCA4R are key in diagnosis, and prognosis as well as when selecting cohorts for therapeutic intervention. In addition, a meticulous documentation of natural history will be invaluable in the future to compare treated with untreated retinas. Despite the familiarity with the term Stargardt disease, this eponymous classification alone is unhelpful when evaluating ABCA4R, as the ABCA4 gene is associated with a number of phenotypes, and a range of severity. Multimodal imaging, psychophysical and electrophysiologic measurements are necessary in diagnosing and characterising these differing retinopathies. A wide range of retinal dystrophy phenotypes are seen in association with ABCA4 mutations. In this article, these will be referred to as ABCA4R. These different phenotypes and the existence of phenocopies present a significant challenge to the clinician. Careful phenotypic characterisation coupled with the genotype enables the clinician to provide an accurate diagnosis, associated inheritance pattern and information regarding prognosis and management. This is particularly relevant now for recruiting to therapeutic trials, and in the future when therapies become available. The importance of accurate genotype-phenotype correlation studies cannot be overemphasised. This approach together with segregation studies can be vital in the identification of causal mutations when variants in more than one gene are being considered as possible. In this article, we give an overview of the current imaging, psychophysical and electrophysiological investigations, as well as current therapeutic research trials for retinopathies associated with the ABCA4 gene.

Mutation Screening in the miR-183/96/182 Cluster in Patients With Inherited Retinal Dystrophy

Inherited retinal dystrophy (IRD) is a heterogenous blinding eye disease and affects more than 200,000 Americans and millions worldwide. By far, 270 protein-coding genes have been identified to cause IRD when defective. However, only one microRNA (miRNA), miR-204, has been reported to be responsible for IRD when a point-mutation occurs in its seed sequence. Previously, we identified that a conserved, polycistronic, paralogous miRNA cluster, the miR-183/96/182 cluster, is highly specifically expressed in all photoreceptors and other sensory organs; inactivation of this cluster in mice resulted in syndromic IRD with multi-sensory defects. We hypothesized that mutations in the miR-183/96/182 cluster in human cause IRD. To test this hypothesis, we perform mutation screening in the pre-miR-183, -96, -182 in >1000 peripheral blood DNA samples of patients with various forms of IRD. We identified six sequence variants, three in pre-miR-182 and three in pre-miR-96. These variants are in the pre-miRNA-182 or -96, but not in the mature miRNAs, and are unlikely to be the cause of the IRD in these patients. In spite of this, the nature and location of these sequence variants in the pre-miRNAs suggest that some may have impact on the biogenesis and maturation of miR-182 or miR-96 and potential roles in the susceptibility to diseases. Although reporting on negative results so far, our study established a system for mutation screening in the miR-183/96/182 cluster in human for a continued effort to unravel and provides deeper insight into the potential roles of miR-183/96/182 cluster in human diseases.

Structural Variants Create New Topological-Associated Domains and Ectopic Retinal Enhancer-Gene Contact in Dominant Retinitis Pigmentosa

The cause of autosomal-dominant retinitis pigmentosa (adRP), which leads to loss of vision and blindness, was investigated in families lacking a molecular diagnosis. A refined locus for adRP on Chr17q22 (RP17) was delineated through genotyping and genome sequencing, leading to the identification of structural variants (SVs) that segregate with disease. Eight different complex SVs were characterized in 22 adRP-affected families with >300 affected individuals. All RP17 SVs had breakpoints within a genomic region spanning YPEL2 to LINC01476. To investigate the mechanism of disease, we reprogrammed fibroblasts from affected individuals and controls into induced pluripotent stem cells (iPSCs) and differentiated them into photoreceptor precursor cells (PPCs) or retinal organoids (ROs). Hi-C was performed on ROs, and differential expression of regional genes and a retinal enhancer RNA at this locus was assessed by qPCR. The epigenetic landscape of the region, and Hi-C RO data, showed that YPEL2 sits within its own topologically associating domain (TAD), rich in enhancers with binding sites for retinal transcription factors. The Hi-C map of RP17 ROs revealed creation of a neo-TAD with ectopic contacts between GDPD1 and retinal enhancers, and modeling of all RP17 SVs was consistent with neo-TADs leading to ectopic retinal-specific enhancer-GDPD1 accessibility. qPCR confirmed increased expression of GDPD1 and increased expression of the retinal enhancer that enters the neo-TAD. Altered TAD structure resulting in increased retinal expression of GDPD1 is the likely convergent mechanism of disease, consistent with a dominant gain of function. Our study highlights the importance of SVs as a genomic mechanism in unsolved Mendelian diseases.

Detailed Phenotyping and Therapeutic Strategies for Intronic ABCA4 Variants in Stargardt Disease

Stargardt disease is a progressive retinal disorder caused by bi-allelic mutations in the ABCA4 gene that encodes the ATP-binding cassette, subfamily A, member 4 transporter protein. Over the past few years, we and others have identified several pathogenic variants that reside within the introns of ABCA4, including a recurrent variant in intron 36 (c.5196+1137G>A) of which the pathogenicity so far remained controversial. Detailed clinical characterization of this variant confirmed its pathogenic nature, and classified it as an allele of intermediate severity. Moreover, we discovered several additional ABCA4 variants clustering in intron 36. Several of these variants resulted in aberrant splicing of ABCA4, i.e., the inclusion of pseudoexons, while the splicing defects caused by the recurrent c.5196+1137G>A variant strongly increased upon differentiation of patient-derived induced pluripotent stem cells into retina-like cells. Finally, all splicing defects could be rescued by the administration of antisense oligonucleotides that were designed to specifically block the pseudoexon insertion, including rescue in 3D retinal organoids harboring the c.5196+1137G>A variant. Our data illustrate the importance of intronic variants in ABCA4 and expand the therapeutic possibilities for overcoming splicing defects in Stargardt disease.

Splicing mutations in inherited retinal diseases

Mutations which induce aberrant transcript splicing represent a distinct class of disease-causing genetic variants in retinal disease genes. Such mutations may either weaken or erase regular splice sites or create novel splice sites which alter exon recognition. While mutations affecting the canonical GU-AG dinucleotides at the splice donor and splice acceptor site are highly predictive to cause a splicing defect, other variants in the vicinity of the canonical splice sites or those affecting additional cis-acting regulatory sequences within exons or introns are much more difficult to assess or even to recognize and require additional experimental validation. Splicing mutations are unique in that the actual outcome for the transcript (e.g. exon skipping, pseudoexon inclusion, intron retention) and the encoded protein can be quite different depending on the individual mutation. In this article, we present an overview on the current knowledge about and impact of splicing mutations in inherited retinal diseases. We introduce the most common sub-classes of splicing mutations including examples from our own work and others and discuss current strategies for the identification and validation of splicing mutations, as well as therapeutic approaches, open questions, and future perspectives in this field of research.

The genetic architecture of Stargardt macular dystrophy (STGD1): a longitudinal 40-year study in a genetic isolate

Stargardt disease (STGD1) is a form of inherited retinal dystrophy attributed to variants affecting function of the large ABCA4 gene and is arguably the most complex monogenic disease. Therapeutic trials in patients depend on identifying causal ABCA4 variants in trans, which is complicated by extreme allelic and clinical heterogeneity. We report the genetic architecture of STGD1 in the young genetically isolated population of Newfoundland, Canada. Population-based clinical recruitment over several decades yielded 29 STGD1 and STGD1-like families (15 multiplex, 14 singleton). Family interviews and public archival records reveal the vast majority of pedigree founders to be of English extraction. Full gene sequencing and haplotype analysis yielded a high solve rate (38/41 cases; 92.7%) for STGD1 and identified 16 causative STGD1 alleles, including a novel deletion (NM_000350.3: ABCA4 c.67-1delG). Several STGD1 alleles of European origin (including NM_000350.3: ABCA4 c.5714 + 5G>A and NM_000350.3: ABCA4 c.5461-10T>C) have drifted to a relatively high population frequency due to founder effect. We report on retinal disease progression in homozygous patients, providing valuable allele-specific insights. The least involved retinal disease is seen in patients homozygous for c.5714 + 5G>A variant, a so-called "mild" variant which is sufficient to precipitate a STGD1 phenotype in the absence of other pathogenic variants in the coding region and intron/exon boundaries of ABCA4. The most severe retinal disease is observed in cases with ABCA4 c.[5461-10T>C;5603A>T] complex allele. We discuss the advantages of determining genetic architecture in genetic isolates in order to begin to meet the grand challenge of human genetics.

Molecular Analysis of the ABCA4 Gene Mutations in Patients with Stargardt Disease Using Human Hair Follicles

ABCA4 gene mutations are the cause of a spectrum of ABCA4 retinopathies, and the most common juvenile macular degeneration is called Stargardt disease. ABCA4 has previously been observed almost exclusively in the retina. Therefore, studying the functional consequences of ABCA4 variants has required advanced molecular analysis techniques. The aim of the present study was to evaluate whether human hair follicles may be used for molecular analysis of the ABCA4 gene splice-site variants in patients with ABCA4 retinopathies. We assessed ABCA4 expression in hair follicles and skin at mRNA and protein levels by means of real-time PCR and Western blot analyses, respectively. We performed cDNA sequencing to reveal the presence of full-length ABCA4 transcripts and analyzed ABCA4 transcripts from three patients with Stargardt disease carrying different splice-site ABCA4 variants: c.5312+1G>A, c.5312+2T>G and c.5836-3C>A. cDNA analysis revealed that c.5312+1G>A, c.5312+2T>G variants led to the skipping of exon 37, and the c.5836-3C>A variant resulted in the insertion of 30 nucleotides into the transcript. Our results strongly argue for the use of hair follicles as a model for the molecular analysis of the pathogenicity of ABCA4 variants in patients with ABCA4 retinopathies.

Clinical spectrum, genetic complexity and therapeutic approaches for retinal disease caused by ABCA4 mutations

The ABCA4 protein (then called a “rim protein”) was first identified in 1978 in the rims and incisures of rod photoreceptors. The corresponding gene, ABCA4, was cloned in 1997, and variants were identified as the cause of autosomal recessive Stargardt disease (STGD1). Over the next two decades, variation in ABCA4 has been attributed to phenotypes other than the classically defined STGD1 or fundus flavimaculatus, ranging from early onset and fast progressing cone-rod dystrophy and retinitis pigmentosa-like phenotypes to very late onset cases of mostly mild disease sometimes resembling, and confused with, age-related macular degeneration. Similarly, analysis of the ABCA4 locus uncovered a trove of genetic information, including >1200 disease-causing mutations of varying severity, and of all types – missense, nonsense, small deletions/insertions, and splicing affecting variants, of which many are located deep-intronic. Altogether, this has greatly expanded our understanding of complexity not only of the diseases caused by ABCA4 mutations, but of all Mendelian diseases in general. This review provides an in depth assessment of the cumulative knowledge of ABCA4-associated retinopathy – clinical manifestations, genetic complexity, pathophysiology as well as current and proposed therapeutic approaches.

Tritan color vision deficiency may be associated with an OPN1SW splicing defect and haploinsufficiency

Here we present evidence implicating disrupted RNA splicing as a potential cause of inherited tritan color vision. Initially we tested 51 subjects for color vision deficiencies. One made significant tritan errors; the others were classified as normal trichromats. The putative tritan subject was the only one of the 51 subjects found to be heterozygous for an OPN1SW gene mutation that disrupts RNA splicing in an in vitro assay. In order to gather further support for the role of the splicing mutation in tritan color vision, the putative tritan subject's mother and sister were examined. They also made tritan errors and had the same OPN1SW gene mutation.

In or Out? New Insights on Exon Recognition through Splice-Site Interdependency

Noncanonical splice-site mutations are an important cause of inherited diseases. Based on in vitro and stem-cell-based studies, some splice-site variants show a stronger splice defect than expected based on their predicted effects, suggesting that other sequence motifs influence the outcome. We investigated whether splice defects due to human-inherited-disease-associated variants in noncanonical splice-site sequences in ABCA4, DMD, and TMC1 could be rescued by strengthening the splice site on the other side of the exon. Noncanonical 5′- and 3′-splice-site variants were selected. Rescue variants were introduced based on an increase in predicted splice-site strength, and the effects of these variants were analyzed using in vitro splice assays in HEK293T cells. Exon skipping due to five variants in noncanonical splice sites of exons in ABCA4, DMD, and TMC1 could be partially or completely rescued by increasing the predicted strengths of the other splice site of the same exon. We named this mechanism “splicing interdependency”, and it is likely based on exon recognition by splicing machinery. Awareness of this interdependency is of importance in the classification of noncanonical splice-site variants associated with disease and may open new opportunities for treatments.

Identification of splice defects due to noncanonical splice site or deep-intronic variants in ABCA4

Pathogenic variants in the ATP-binding cassette transporter A4 (ABCA4) gene cause a continuum of retinal disease phenotypes, including Stargardt disease. Noncanonical splice site (NCSS) and deep-intronic variants constitute a large fraction of disease-causing alleles, defining the functional consequences of which remains a challenge. We aimed to determine the effect on splicing of nine previously reported or unpublished NCSS variants, one near exon splice variant and nine deep-intronic variants in ABCA4, using in vitro splice assays in human embryonic kidney 293T cells. Reverse transcription-polymerase chain reaction and Sanger sequence analysis revealed splicing defects for 12 out of 19 variants. Four deep-intronic variants create pseudoexons or elongate the upstream exon. Furthermore, eight NCSS variants cause a partial deletion or skipping of one or more exons in messenger RNAs. Among the 12 variants, nine lead to premature stop codons and predicted truncated ABCA4 proteins. At least two deep-intronic variants affect splice enhancer and silencer motifs and, therefore, these conserved sequences should be carefully evaluated when predicting the outcome of NCSS and deep-intronic variants.

Machine Learning Approaches for the Prioritization of Genomic Variants Impacting Pre-mRNA Splicing

Defects in pre-mRNA splicing are frequently a cause of Mendelian disease. Despite the advent of next-generation sequencing, allowing a deeper insight into a patient's variant landscape, the ability to characterize variants causing splicing defects has not progressed with the same speed. To address this, recent years have seen a sharp spike in the number of splice prediction tools leveraging machine learning approaches, leaving clinical geneticists with a plethora of choices for in silico analysis. In this review, some basic principles of machine learning are introduced in the context of genomics and splicing analysis. A critical comparative approach is then used to describe seven recent machine learning-based splice prediction tools, revealing highly diverse approaches and common caveats. We find that, although great progress has been made in producing specific and sensitive tools, there is still much scope for personalized approaches to prediction of variant impact on splicing. Such approaches may increase diagnostic yields and underpin improvements to patient care.

Prevalence of ABCA4 Deep-Intronic Variants and Related Phenotype in An Unsolved "One-Hit" Cohort with Stargardt Disease

We investigated the prevalence of reported deep-intronic variants in a French cohort of 70 patients with Stargardt disease harboring a monoallelic pathogenic variant on the exonic regions of ABCA4. Direct Sanger sequencing of selected intronic regions of ABCA4 was conducted. Complete phenotypic analysis and correlation with the genotype was performed in case a known intronic pathogenic variant was identified. All other variants found on the analyzed sequences were queried for minor allele frequency and possible pathogenicity by in silico predictions. The second mutated allele was found in 14 (20%) subjects. The three known deep-intronic variants found were c.5196+1137G>A in intron 36 (6 subjects), c.4539+2064C>T in intron 30 (4 subjects) and c.4253+43G>A in intron 28 (4 subjects). Even though the phenotype depends on the compound effect of the biallelic variants, a genotype-phenotype correlation suggests that the c.5196+1137G>A was mostly associated with a mild phenotype and the c.4539+2064C>T with a more severe one. A variable effect was instead associated with the variant c.4253+43G>A. In addition, two novel variants, c.768+508A>G and c.859-245_859-243delinsTGA never associated with Stargardt disease before, were identified and a possible splice defect was predicted in silico. Our study calls for a larger cohort analysis including targeted locus sequencing and 3D protein modeling to better understand phenotype-genotype correlations associated with deep-intronic changes and patients’ selection for clinical trials.

HYPERREFLECTIVE DEPOSITION IN THE BACKGROUND OF ADVANCED STARGARDT DISEASE

PURPOSE

To describe an unusual manifestation of hyperreflective deposits in the subretinal space in a group of patients with clinically and genetically confirmed Stargardt disease.

METHODS

Retrospective review of color fundus, autofluorescence, infrared reflectance, red-free images, and spectral domain optical coherence tomography in 296 clinically diagnosed and genetically confirmed (2 expected disease-causing mutations in ABCA4) patients with Stargardt disease. Full-field electroretinogram (ffERG), medical history, and genotype data (in silico predictions) were further analyzed from the selected cohort.

RESULTS

Eight of 296 patients (2.7%) were found to exhibit small crystalline deposits that were detectable on certain imaging modalities, such as color, infrared reflectance and red-free images, but not autofluorescence. The deposits were most prevalent in the superior region of the macula, and spectral domain optical coherence tomography revealed their presence in the subretinal space. All patients presented with these findings at a notably advanced disease stage with abnormal ffERG and a high proportion of highly deleterious ABCA4 alleles.

CONCLUSION

Hyperreflective subretinal deposits may be a manifestation of advanced ABCA4 disease, particularly in regions susceptible to disease-related changes, such as lipofuscin accumulation.