Autozygome-guided exome-first study in a consanguineous cohort with early-onset retinal disease uncovers an isolated RIMS2 phenotype and a retina-enriched RIMS2 isoform.

Leber congenital amaurosis (LCA) and early-onset retinal degeneration (EORD) are inherited retinal diseases (IRD) characterized by early-onset vision impairment. Herein, we studied 15 Saudi families by whole exome sequencing (WES) and run-of-homozygosity (ROH) detection via AutoMap in 12/15 consanguineous families. This revealed (likely) pathogenic variants in 11/15 families (73%). A potential founder variant was found in RPGRIP1. Homozygous pathogenic variants were identified in known IRD genes (ATF6, CRB1, CABP4, RDH12, RIMS2, RPGRIP1, SPATA7). We established genotype-driven clinical reclassifications for ATF6, CABP4, and RIMS2. Specifically, we observed isolated IRD in the individual with the novel RIMS2 variant, and we found a retina-enriched RIMS2 isoform conserved but not annotated in mouse. The latter illustrates potential different phenotypic consequences of pathogenic variants depending on the particular tissue/cell-type specific isoforms they affect. Lastly, a compound heterozygous genotype in GUCY2D in one non-consanguineous family was demonstrated, and homozygous variants in novel candidate genes ATG2B and RUFY3 were found in the two remaining consanguineous families. Reporting these genes will allow to validate them in other IRD cohorts. Finally, the missing heritability of the two unsolved IRD cases may be attributed to variants in non-coding regions or structural variants that remained undetected, warranting future WGS studies.

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 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.

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.