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.

Deep intronic variant causes aberrant splicing of ATP7A in a family with a variable occipital horn syndrome phenotype.

Genetic variants in ATP7A are associated with a spectrum of X-linked disorders. In descending order of severity, these are Menkes disease, occipital horn syndrome, and X-linked distal spinal muscular atrophy. After 30 years of diagnostic investigation, we identified a deep intronic ATP7A variant in four males from a family affected to variable degrees by a predominantly skeletal phenotype, featuring bowing of long bones, elbow joints with restricted mobility which dislocate frequently, coarse curly hair, chronic diarrhoea, and motor coordination difficulties. Analysis of whole genome sequencing data from the Genomics England 100,000 Genomes Project following clinical re-evaluation identified a deep intronic ATP7A variant, which was predicted by SpliceAI to have a modest splicing effect. Using a mini-gene splicing assay, we determined that the intronic variant results in aberrant splicing. Sanger sequencing of patient cDNA revealed ATP7A transcripts with exon 5 skipping, or inclusion of a novel intron 4 pseudoexon. In both instances, frameshift leading to premature termination are predicted. Quantification of ATP7A mRNA transcripts using a qPCR assay indicated that the majority of transcripts (86.1 %) have non-canonical splicing, with 68.0 % featuring exon 5 skipping, and 18.1 % featuring the novel pseudoexon. We suggest that the variability of the phenotypes within the affected males results from the stochastic effects of splicing. This deep intronic variant, resulting in aberrant ATP7A splicing, expands the understanding of intronic variation on the ATP7A-related disease spectrum.