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.

Functional characterization of a novel non-coding mutation "Ghent +49A > G" in the iron-responsive element of L-ferritin causing hereditary hyperferritinaemia-cataract syndrome.

Hereditary hyperferritinaemia-cataract syndrome (HHCS) is a rare disorder usually caused by heterozygous mutations in the iron-responsive element (IRE) in the 5' untranslated region (5'UTR) of the L-ferritin gene (FTL), disturbing the binding of iron regulatory proteins (IRPs) and the post-transcriptional regulation of ferritin expression. Here, the proband of a consanguineous family displayed moderate bilateral cataracts and elevated serum ferritin in the absence of iron overload. The parents and siblings showed variable degrees of mild bilateral cataracts combined with elevated levels of circulating ferritin. Sequencing of FTL identified a novel 5'UTR mutation c.-151A > G, also named "Ghent +49A > G". The zygosity of the mutation, occurring in homozygous and heterozygous state in the proband and other affected family members respectively, correlated well with severity of ophthalmological and hematological manifestations. The substitution is expected to impair the secondary structure of the upper IRE stem. Functional characterization of +49A > G by electrophoretic mobility shift assays demonstrated a reduced binding affinity for IRP1 compared to the wild-type IRE of FTL. Overall, we have expanded the repertoire of deleterious biallelic FTL IRE mutations in HHCS with this novel +49A > G mutation, the zygosity of which correlated well with the disease expression.