ISOLATED MACULOPATHY AND MODERATE ROD-CONE DYSTROPHY REPRESENT THE MILDER END OF THE RDH12-RELATED RETINAL DYSTROPHY SPECTRUM.

PURPOSE: To describe an isolated maculopathy and an intermediate rod-cone dystrophy phenotype as the milder end of the RDH12-related retinal dystrophy spectrum. METHODS: Seven patients (17-34 years of age) underwent an extensive ophthalmic workup including psychophysical and electrophysiological testing and multimodal imaging. RESULTS: Three patients have isolated macular disease. Best-corrected visual acuity (BCVA) ranges from 20/125 to 20/40 with normal visual fields or only limited central, relative scotomata, and normal full-field ERGs. Both optical coherence tomography scans and autofluorescent imaging hint at relatively better-preserved foveal quality initially. An intermediate rod-cone phenotype in four patients is characterized by a central retinal dystrophy extending just beyond the vascular arcades, characteristic peripapillary sparing, and additional scattered atrophic patches. Again, foveal quality is initially better on optical coherence tomography scans. Best-corrected visual acuity ranges from counting fingers to 20/32. Goldmann visual fields vary from central scotomata to severe generalized abnormalities. ERGs range between mild and severe rod-cone dysfunction. Nine distinct RDH12 pathogenic variants, two of which are novel, are identified. CONCLUSION: The classic phenotype of RDH12-related early-onset retinal dystrophy is expanded to include an isolated maculopathy and intermediate dystrophy phenotype, characterized by its later onset and milder course with a fair visual potential until much later in life, emphasizing the phenotypic heterogeneity of RDH12-related retinopathy.

Functional characterization of the first missense variant in CEP78, a founder allele associated with cone-rod dystrophy, hearing loss, and reduced male fertility.

Inactivating variants in the centrosomal CEP78 gene have been found in cone-rod dystrophy with hearing loss (CRDHL), a particular phenotype distinct from Usher syndrome. Here, we identified and functionally characterized the first CEP78 missense variant c.449T>C, p.(Leu150Ser) in three CRDHL families. The variant was found in a biallelic state in two Belgian families and in a compound heterozygous state-in trans with c.1462-1G>T-in a third German family. Haplotype reconstruction showed a founder effect. Homology modeling revealed a detrimental effect of p.(Leu150Ser) on protein stability, which was corroborated in patients' fibroblasts. Elongated primary cilia without clear ultrastructural abnormalities in sperm or nasal brushes suggest impaired cilia assembly. Two affected males from different families displayed sperm abnormalities causing infertility. One of these is a heterozygous carrier of a complex allele in SPAG17, a ciliary gene previously associated with autosomal recessive male infertility. Taken together, our data indicate that a missense founder allele in CEP78 underlies the same sensorineural CRDHL phenotype previously associated with inactivating variants. Interestingly, the CEP78 phenotype has been possibly expanded with male infertility. Finally, CEP78 loss-of-function variants may have an underestimated role in misdiagnosed Usher syndrome, with or without sperm abnormalities.

Functional characterization of novel MFSD8 pathogenic variants anticipates neurological involvement in juvenile isolated maculopathy.

Biallelic MFSD8 variants are an established cause of severe late-infantile subtype of neuronal ceroid lipofuscinosis (v-LINCL), a severe lysosomal storage disorder, but have also been associated with nonsyndromic adult-onset maculopathy. Here, we functionally characterized two novel MFSD8 variants found in a child with juvenile isolated maculopathy, in order to establish a refined prognosis. ABCA4 locus resequencing was followed by the analysis of other inherited retinal disease genes by whole exome sequencing (WES). Minigene assays and cDNA sequencing were used to assess the effect of a novel MFSD8 splice variant. MFSD8 expression was quantified with qPCR and overexpression studies were analyzed by immunoblotting. Transmission electron microscopy (TEM) was performed on a skin biopsy and ophthalmological and neurological re-examinations were conducted. WES revealed two novel MFSD8 variants: c.[590del];[439+3A>C] p.[Gly197Valfs*2];[Ile67Glufs*3]. Characterization of the c.439+3A>C variant via splice assays showed exon-skipping (p.Ile67Glufs*3), while overexpression studies of the corresponding protein indicated expression of a truncated polypeptide. In addition, a significantly reduced MFSD8 RNA expression was noted in patient's lymphocytes. TEM of a skin biopsy revealed typical v-LINCL lipopigment inclusions while neurological imaging of the proband displayed subtle cerebellar atrophy. Functional characterization demonstrated the pathogenicity of two novel MFSD8 variants, found in a child with an initial diagnosis of juvenile isolated maculopathy but likely evolving to v-LINCL with a protracted disease course. Our study allowed a refined neurological prognosis in the proband and expands the natural history of MFSD8-associated disease.

The N-terminal p.(Ser38Cys) TIMP3 mutation underlying Sorsby fundus dystrophy is a founder mutation disrupting an intramolecular disulfide bond.

Sorsby fundus dystrophy (SFD) is a macular degeneration caused by mutations in TIMP3, the majority of which introduce a novel cysteine. However, the exact molecular mechanisms underlying SFD remain unknown. We aimed to provide novel insights into the functional consequences of a distinct N-terminal mutation. Haplotype reconstruction in three SFD families revealed that the identified c.113C>G, p.(Ser38Cys) mutation is a founder in Belgian and northern French families with a late-onset SFD phenotype. Functional consequences of the p.(Ser38Cys) mutation were investigated by high-resolution Western blot analysis of wild type and mutant TIMP3 using patient fibroblasts and in vitro generated proteins, and by molecular modeling of TIMP3 and its interaction partners. We could not confirm a previous hypothesis on dimerization of mutant TIMP3 proteins. However, we identified aberrant intramolecular disulfide bonding. Our data provide evidence for disruption of the established Cys36-Cys143 disulfide bond and formation of a novel Cys36-Cys38 bond, possibly associated with increased glycosylation of the protein. In conclusion, we propose a novel pathogenetic mechanism underlying the p.(Ser38Cys) TIMP3 founder mutation involving intramolecular disulfide bonding. These results provide new insights into the pathogenesis of SFD and other retinopathies linked to mutations in TIMP3, such as age-related macular degeneration.

Isolated and Syndromic Retinal Dystrophy Caused by Biallelic Mutations in RCBTB1, a Gene Implicated in Ubiquitination.

Inherited retinal dystrophies (iRDs) are a group of genetically and clinically heterogeneous conditions resulting from mutations in over 250 genes. Here, homozygosity mapping and whole-exome sequencing (WES) in a consanguineous family revealed a homozygous missense mutation, c.973C>T (p.His325Tyr), in RCBTB1. In affected individuals, it was found to segregate with retinitis pigmentosa (RP), goiter, primary ovarian insufficiency, and mild intellectual disability. Subsequent analysis of WES data in different cohorts uncovered four additional homozygous missense mutations in five unrelated families in whom iRD segregates with or without syndromic features. Ocular phenotypes ranged from typical RP starting in the second decade to chorioretinal dystrophy with a later age of onset. The five missense mutations affect highly conserved residues either in the sixth repeat of the RCC1 domain or in the BTB1 domain. A founder haplotype was identified for mutation c.919G>A (p.Val307Met), occurring in two families of Mediterranean origin. We showed ubiquitous mRNA expression of RCBTB1 and demonstrated predominant RCBTB1 localization in human inner retina. RCBTB1 was very recently shown to be involved in ubiquitination, more specifically as a CUL3 substrate adaptor. Therefore, the effect on different components of the CUL3 and NFE2L2 (NRF2) pathway was assessed in affected individuals' lymphocytes, revealing decreased mRNA expression of NFE2L2 and several NFE2L2 target genes. In conclusion, our study puts forward mutations in RCBTB1 as a cause of autosomal-recessive non-syndromic and syndromic iRD. Finally, our data support a role for impaired ubiquitination in the pathogenetic mechanism of RCBTB1 mutations.

Correction: Mapping the genomic landscape of inherited retinal disease genes prioritizes genes prone to coding and noncoding copy-number variations.

The original version of this Article contained an error in the spelling of the author Anja K. Mayer, which was incorrectly given as Anja Kathrin Mayer. This has now been corrected in both the PDF and HTML versions of the Article.

Correction: Biallelic sequence and structural variants in RAX2 are a novel cause for autosomal recessive inherited retinal disease.

The original version of this Article contained an incorrect version of Fig. 3, which included two variants initially shown in black text in Fig. 3a that the authors removed from the final manuscript. The correct version of Fig. 3 without the two variants now appears in the PDF and HTML versions of the Article.

Biallelic sequence and structural variants in RAX2 are a novel cause for autosomal recessive inherited retinal disease.

PURPOSE: RAX2 encodes a homeobox-containing transcription factor, in which four monoallelic pathogenic variants have been described in autosomal dominant cone-dominated retinal disease. METHODS: Exome sequencing in a European cohort with inherited retinal disease (IRD) (n = 2086) was combined with protein structure modeling of RAX2 missense variants, bioinformatics analysis of deletion breakpoints, haplotyping of RAX2 variant c.335dup, and clinical assessment of biallelic RAX2-positive cases and carrier family members. RESULTS: Biallelic RAX2 sequence and structural variants were found in five unrelated European index cases, displaying nonsyndromic autosomal recessive retinitis pigmentosa (ARRP) with an age of onset ranging from childhood to the mid-40s (average mid-30s). Protein structure modeling points to loss of function of the novel recessive missense variants and to a dominant-negative effect of the reported dominant RAX2 alleles. Structural variants were fine-mapped to disentangle their underlying mechanisms. Haplotyping of c.335dup in two cases suggests a common ancestry. CONCLUSION: This study supports a role for RAX2 as a novel disease gene for recessive IRD, broadening the mutation spectrum from sequence to structural variants and revealing a founder effect. The identification of biallelic RAX2 pathogenic variants in five unrelated families shows that RAX2 loss of function may be a nonnegligible cause of IRD in unsolved ARRP cases.

Mapping the genomic landscape of inherited retinal disease genes prioritizes genes prone to coding and noncoding copy-number variations.

PurposePart of the hidden genetic variation in heterogeneous genetic conditions such as inherited retinal diseases (IRDs) can be explained by copy-number variations (CNVs). Here, we explored the genomic landscape of IRD genes listed in RetNet to identify and prioritize those genes susceptible to CNV formation.MethodsRetNet genes underwent an assessment of genomic features and of CNV occurrence in the Database of Genomic Variants and literature. CNVs identified in an IRD cohort were characterized using targeted locus amplification (TLA) on extracted genomic DNA.ResultsExhaustive literature mining revealed 1,345 reported CNVs in 81 different IRD genes. Correlation analysis between rankings of genomic features and CNV occurrence demonstrated the strongest correlation between gene size and CNV occurrence of IRD genes. Moreover, we identified and delineated 30 new CNVs in IRD cases, 13 of which are novel and three of which affect noncoding, putative cis-regulatory regions. Finally, the breakpoints of six complex CNVs were determined using TLA in a hypothesis-neutral manner.ConclusionWe propose a ranking of CNV-prone IRD genes and demonstrate the efficacy of TLA for the characterization of CNVs on extracted DNA. Finally, this IRD-oriented CNV study can serve as a paradigm for other genetically heterogeneous Mendelian diseases with hidden genetic variation.

High-throughput PCR assay design for targeted resequencing using primerXL.

BACKGROUND: Although the sequencing landscape is rapidly evolving and sequencing costs are continuously decreasing, whole genome sequencing is still too expensive for use on a routine basis. Targeted resequencing of only the regions of interest decreases both costs and the complexity of the downstream data-analysis. Various target enrichment strategies are available, but none of them obtain the degree of coverage uniformity, flexibility and specificity of PCR-based enrichment. On the other hand, the biggest limitation of target enrichment by PCR is the need to design large numbers of partially overlapping assays to cover the target. RESULTS: To overcome the aforementioned hurdles, we have developed primerXL, a state-of-the-art PCR primer design pipeline for targeted resequencing. It uses an optimized design criteria relaxation cascade and a thorough downstream in silico evaluation process to generate high quality singleplex PCR assays, reducing the need for amplicon normalization, and outperforming other target enrichment strategies and similar primer design tools when considering assay quality, coverage uniformity and target coverage. Results of four different sequencing projects with 2348 amplicons in total covering 470 kb are presented. PrimerXL can be accessed at www.primerxl.org . CONCLUSION: PrimerXL is an state-of-the-art, easy to use primer design webtool capable of generating high-quality targeted resequencing assays. The workflow is fully customizable to suit every researchers' needs, while an innovative relaxation cascade ensures maximal target coverage.

NR5A1 is a novel disease gene for 46,XX testicular and ovotesticular disorders of sex development.

PURPOSE: We aimed to identify the genetic cause in a cohort of 11 unrelated cases and two sisters with 46,XX SRY-negative (ovo)testicular disorders of sex development (DSD). METHODS: Whole-exome sequencing (n = 9), targeted resequencing (n = 4), and haplotyping were performed. Immunohistochemistry of sex-specific markers was performed on patients' gonads. The consequences of mutation were investigated using luciferase assays, localization studies, and RNA-seq. RESULTS: We identified a novel heterozygous NR5A1 mutation, c.274C>T p.(Arg92Trp), in three unrelated patients. The Arg92 residue is highly conserved and located in the Ftz-F1 region, probably involved in DNA-binding specificity and stability. There were no consistent changes in transcriptional activation or subcellular localization. Transcriptomics in patient-derived lymphocytes showed upregulation of MAMLD1, a direct NR5A1 target previously associated with 46,XY DSD. In gonads of affected individuals, ovarian FOXL2 and testicular SRY-independent SOX9 expression observed. CONCLUSIONS: We propose NR5A1, previously associated with 46,XY DSD and 46,XX primary ovarian insufficiency, as a novel gene for 46,XX (ovo)testicular DSD. We hypothesize that p.(Arg92Trp) results in decreased inhibition of the male developmental pathway through downregulation of female antitestis genes, thereby tipping the balance toward testicular differentiation in 46,XX individuals. In conclusion, our study supports a role for NR5A1 in testis differentiation in the XX gonad.Genet Med 19 4, 367-376.

arrEYE: a customized platform for high-resolution copy number analysis of coding and noncoding regions of known and candidate retinal dystrophy genes and retinal noncoding RNAs.

PURPOSE: Our goal was to design a customized microarray, arrEYE, for high-resolution copy number variant (CNV) analysis of known and candidate genes for inherited retinal dystrophy (iRD) and retina-expressed noncoding RNAs (ncRNAs). METHODS: arrEYE contains probes for the full genomic region of 106 known iRD genes, including those implicated in retinitis pigmentosa (RP) (the most frequent iRD), cone-rod dystrophies, macular dystrophies, and an additional 60 candidate iRD genes and 196 ncRNAs. Eight CNVs in iRD genes identified by other techniques were used as positive controls. The test cohort consisted of 57 patients with autosomal dominant, X-linked, or simplex RP. RESULTS: In an RP patient, a novel heterozygous deletion of exons 7 and 8 of the HGSNAT gene was identified: c.634-408_820+338delinsAGAATATG, p.(Glu212Glyfs*2). A known variant was found on the second allele: c.1843G>A, p.(Ala615Thr). Furthermore, we expanded the allelic spectrum of USH2A and RCBTB1 with novel CNVs. CONCLUSION: The arrEYE platform revealed subtle single-exon to larger CNVs in iRD genes that could be characterized at the nucleotide level, facilitated by the high resolution of the platform. We report the first CNV in HGSNAT that, combined with another mutation, leads to RP, further supporting its recently identified role in nonsyndromic iRD.Genet Med 19 4, 457-466.

Hidden Genetic Variation in LCA9-Associated Congenital Blindness Explained by 5'UTR Mutations and Copy-Number Variations of NMNAT1.

Leber congenital amaurosis (LCA) is a severe autosomal-recessive retinal dystrophy leading to congenital blindness. A recently identified LCA gene is NMNAT1, located in the LCA9 locus. Although most mutations in blindness genes are coding variations, there is accumulating evidence for hidden noncoding defects or structural variations (SVs). The starting point of this study was an LCA9-associated consanguineous family in which no coding mutations were found in the LCA9 region. Exploring the untranslated regions of NMNAT1 revealed a novel homozygous 5'UTR variant, c.-70A>T. Moreover, an adjacent 5'UTR variant, c.-69C>T, was identified in a second consanguineous family displaying a similar phenotype. Both 5'UTR variants resulted in decreased NMNAT1 mRNA abundance in patients' lymphocytes, and caused decreased luciferase activity in human retinal pigment epithelial RPE-1 cells. Second, we unraveled pseudohomozygosity of a coding NMNAT1 mutation in two unrelated LCA patients by the identification of two distinct heterozygous partial NMNAT1 deletions. Molecular characterization of the breakpoint junctions revealed a complex Alu-rich genomic architecture. Our study uncovered hidden genetic variation in NMNAT1-associated LCA and emphasized a shift from coding to noncoding regulatory mutations and repeat-mediated SVs in the molecular pathogenesis of heterogeneous recessive disorders such as hereditary blindness.

Flexible, scalable, and efficient targeted resequencing on a benchtop sequencer for variant detection in clinical practice.

The release of benchtop next-generation sequencing (NGS) instruments has paved the way to implement the technology in clinical setting. The need for flexible, qualitative, and cost-efficient workflows is high. We used singleplex-PCR for highly efficient target enrichment, allowing us to reach the quality standards set in Sanger sequencing-based diagnostics. For the library preparation, a modified NexteraXT protocol was used, followed by sequencing on a MiSeq instrument. With an innovative pooling strategy, high flexibility, scalability, and cost-efficiency were obtained, independent of the availability of commercial kits. The approach was validated for ∼250 genes associated with monogenic disorders. An overall sensitivity (>99%) similar to Sanger sequencing was observed in combination with a positive predictive value of >98%. The distribution of coverage was highly uniform, guaranteeing a minimal number of gaps to be filled with alternative methods. ISO15189-accreditation was obtained for the workflow. A major asset of the singleplex PCR-based enrichment is that new targets can be easily implemented. Diagnostic laboratories have validated assays available ensuring that the proposed workflow can easily be adopted. Although our platform was optimized for constitutional variant detection of monogenic disease genes, it is now also used as a model for somatic mutation detection in acquired diseases.

An augmented ABCA4 screen targeting noncoding regions reveals a deep intronic founder variant in Belgian Stargardt patients.

Autosomal-recessive Stargardt disease (STGD1) is hallmarked by a large proportion of patients with a single heterozygous causative variant in the disease gene ABCA4. Braun et al. () reported deep intronic variants of ABCA4 in STGD1 patients with one coding variant, prompting us to perform an augmented screen in 131 Belgian STGD1 patients with one or no ABCA4 variant to uncover deep intronic causal ABCA4 variants. This revealed a second variant in 28.6% of cases. Twenty-six percent of these carry the same causal variant c.4539+2001G>A (V4). Haplotyping in V4 carriers showed a common region of 63 kb, suggestive of a founder mutation. Genotype-phenotype correlations suggest a moderate-to-severe impact of V4 on the STGD1 phenotype. In conclusion, V4 occurs in a high fraction of Belgian STGD1 patients and represents the first deep intronic founder mutation in ABCA4. This emphasizes the importance of augmented molecular genetic testing of ABCA4 in Belgian STGD1.

A nonsense mutation in PDE6H causes autosomal-recessive incomplete achromatopsia.

Achromatopsia (ACHM) is an autosomal-recessive retinal dystrophy characterized by color blindness, photophobia, nystagmus, and severely reduced visual acuity. Its prevalence has been estimated to about 1 in 30,000 individuals. Four genes, GNAT2, PDE6C, CNGA3, and CNGB3, have been implicated in ACHM, and all encode functional components of the phototransduction cascade in cone photoreceptors. Applying a functional-candidate-gene approach that focused on screening additional genes involved in this process in a cohort of 611 index cases with ACHM or other cone photoreceptor disorders, we detected a homozygous single base change (c.35C>G) resulting in a nonsense mutation (p.Ser12(∗)) in PDE6H, encoding the inhibitory γ subunit of the cone photoreceptor cyclic guanosine monophosphate phosphodiesterase. The c.35C>G mutation was present in three individuals from two independent families with a clinical diagnosis of incomplete ACHM and preserved short-wavelength-sensitive cone function. Moreover, we show through immunohistochemical colocalization studies in mouse retina that Pde6h is evenly present in all retinal cone photoreceptors, a fact that had been under debate in the past. These findings add PDE6H to the set of genes involved in autosomal-recessive cone disorders and demonstrate the importance of the inhibitory γ subunit in cone phototransduction.

Early-onset autosomal recessive cerebellar ataxia associated with retinal dystrophy: new human hotfoot phenotype caused by homozygous GRID2 deletion.

PURPOSE: The aim of this study was to identify the genetic cause of early-onset autosomal recessive cerebellar ataxia associated with retinal dystrophy in a consanguineous family. METHODS: An affected 6-month-old child underwent neurological and ophthalmological examinations. Genetic analyses included homozygosity mapping, copy number analysis, conventional polymerase chain reaction, Sanger sequencing, quantitative polymerase chain reaction, and whole-exome sequencing. Expression analysis of GRID2 was performed by quantitative polymerase chain reaction and immunohistochemistry. RESULTS: A homozygous deletion of exon 2 of GRID2 (p.Gly30_Glu81del) was identified in the proband. GRID2 encodes an ionotropic glutamate receptor known to be selectively expressed in cerebellar Purkinje cells. Here, we demonstrated GRID2 expression in human adult retina and retinal pigment epithelium. In addition, Grid2 expression was demonstrated in different stages of murine retinal development. GRID2 immunostaining was shown in murine and human retina. Whole-exome sequencing in the proband did not provide arguments for other disease-causing mutations, supporting the idea that the phenotype observed represents a single clinical entity. CONCLUSION: We identified GRID2 as an underlying disease gene of early-onset autosomal recessive cerebellar ataxia with retinal dystrophy, expanding the clinical spectrum of GRID2 deletion mutants. We demonstrated for the first time GRID2 expression and localization in human and murine retina, providing evidence for a novel functional role of GRID2 in the retina.

Identity-by-descent-guided mutation analysis and exome sequencing in consanguineous families reveals unusual clinical and molecular findings in retinal dystrophy.

PURPOSE: Autosomal recessive retinal dystrophies are clinically and genetically heterogeneous, which hampers molecular diagnosis. We evaluated identity-by-descent-guided Sanger sequencing or whole-exome sequencing in 26 families with nonsyndromic (19) or syndromic (7) autosomal recessive retinal dystrophies to identify disease-causing mutations. METHODS: Patients underwent genome-wide identity-by-descent mapping followed by Sanger sequencing (16) or whole-exome sequencing (10). Whole-exome sequencing data were filtered against identity-by-descent regions and known retinal dystrophy genes. The medical history was reviewed in mutation-positive families. RESULTS: We identified mutations in 14 known retinal dystrophy genes in 20/26 (77%) families: ABCA4, CERKL, CLN3, CNNM4, C2orf71, IQCB1, LRAT, MERTK, NMNAT1, PCDH15, PDE6B, RDH12, RPGRIP1, and USH2A. Whole-exome sequencing in single individuals revealed mutations in either the largest or smaller identity-by-descent regions, and a compound heterozygous genotype in NMNAT1. Moreover, a novel deletion was found in PCDH15. In addition, we identified mutations in CLN3, CNNM4, and IQCB1 in patients initially diagnosed with nonsyndromic retinal dystrophies. CONCLUSION: Our study emphasized that identity-by-descent-guided mutation analysis and/or whole-exome sequencing are powerful tools for the molecular diagnosis of retinal dystrophy. Our approach uncovered unusual molecular findings and unmasked syndromic retinal dystrophies, guiding future medical management. Finally, elucidating ABCA4, LRAT, and MERTK mutations offers potential gene-specific therapeutic perspectives.

Multi-omics analysis in human retina uncovers ultraconserved cis-regulatory elements at rare eye disease loci.

Cross-species genome comparisons have revealed a substantial number of ultraconserved non-coding elements (UCNEs). Several of these elements have proved to be essential tissue- and cell type-specific cis-regulators of developmental gene expression. Here, we characterize a set of UCNEs as candidate CREs (cCREs) during retinal development and evaluate the contribution of their genomic variation to rare eye diseases, for which pathogenic non-coding variants are emerging. Integration of bulk and single-cell retinal multi-omics data reveals 594 genes under potential cis-regulatory control of UCNEs, of which 45 are implicated in rare eye disease. Mining of candidate cis-regulatory UCNEs in WGS data derived from the rare eye disease cohort of Genomics England reveals 178 ultrarare variants within 84 UCNEs associated with 29 disease genes. Overall, we provide a comprehensive annotation of ultraconserved non-coding regions acting as cCREs during retinal development which can be targets of non-coding variation underlying rare eye diseases.

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

BACKGROUND: 5' untranslated regions (5'UTRs) are essential modulators of protein translation. Predicting the impact of 5'UTR variants is challenging and rarely performed in routine diagnostics. Here, we present a combined approach of a comprehensive prioritization strategy and functional assays to evaluate 5'UTR variation in two large cohorts of patients with inherited retinal diseases (IRDs). METHODS: We performed an isoform-level re-analysis of retinal RNA-seq data to identify the protein-coding transcripts of 378 IRD genes with highest expression in retina. We evaluated the coverage of their 5'UTRs by different whole exome sequencing (WES) kits. The selected 5'UTRs were analyzed in whole genome sequencing (WGS) and WES data from IRD sub-cohorts from the 100,000 Genomes Project (n = 2397 WGS) and an in-house database (n = 1682 WES), respectively. Identified variants were annotated for 5'UTR-relevant features and classified into seven categories based on their predicted functional consequence. We developed a variant prioritization strategy by integrating population frequency, specific criteria for each category, and family and phenotypic data. A selection of candidate variants underwent functional validation using diverse approaches. RESULTS: Isoform-level re-quantification of retinal gene expression revealed 76 IRD genes with a non-canonical retina-enriched isoform, of which 20 display a fully distinct 5'UTR compared to that of their canonical isoform. Depending on the probe design, 3-20% of IRD genes have 5'UTRs fully captured by WES. After analyzing these regions in both cohorts, we prioritized 11 (likely) pathogenic variants in 10 genes (ARL3, MERTK, NDP, NMNAT1, NPHP4, PAX6, PRPF31, PRPF4, RDH12, RD3), of which 7 were novel. Functional analyses further supported the pathogenicity of three variants. Mis-splicing was demonstrated for the PRPF31:c.-9+1G>T variant. The MERTK:c.-125G>A variant, overlapping a transcriptional start site, was shown to significantly reduce both luciferase mRNA levels and activity. The RDH12:c.-123C>T variant was found in cis with the hypomorphic RDH12:c.701G>A (p.Arg234His) variant in 11 patients. This 5'UTR variant, predicted to introduce an upstream open reading frame, was shown to result in reduced RDH12 protein but unaltered mRNA levels. CONCLUSIONS: This study demonstrates the importance of 5'UTR variants implicated in IRDs and provides a systematic approach for 5'UTR annotation and validation that is applicable to other inherited diseases.