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

Direct comparison of distinct naive pluripotent states in human embryonic stem cells

Stem cells & developmental biolog

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

New noncoding base pair mutation at the identical locus as the original NCMD/MCDR1 in a Mexican family, suggesting a mutational hotspot

Purpose:To clinically and molecularly study a newly found family with North Carolina macular dystrophy (NCMD/MCDR1) from Mexico. Methods:This retrospective study comprised 6 members of a 3-generation Mexican family with NCMD. Clinical ophthalmic examinations, including fundus imaging, spectral-domain optical coherence tomography, electroretinography, and electrooculography, were performed. Genotyping with polymorphic markers in the MCDR1 region was performed to determine haplotypes. Whole-genome sequencing (WGS) was performed followed by variant filtering and copy number variant analysis. Results:Four subjects from 3 generations were found to have macular abnormalities. The proband presented with lifelong bilateral vision impairment with bilaterally symmetric vitelliform Best disease-like appearing macular lesions. Her 2 children had bilateral large macular coloboma-like malformations, consistent with autosomal dominant NCMD. The 80-year-old mother of the proband had drusen-like lesions consistent with grade 1 NCMD. WGS and subsequent Sanger sequencing found a point mutation at chr6:99593030G>C (hg38) in the noncoding region of the DNase I site thought to be a regulatory element of the retinal transcription factor gene PRDM13. This mutation is the identical site/nucleotide as in the original NCMD family (#765) but is a guanine to cytosine change rather than a guanine to thymine mutation, as found in the original NCMD family. Conclusions:We report a new noncoding mutation at the same locus (chr6:99593030G>C) involving the same DNase I site regulating the retinal transcription factor gene PRDM13. This suggests that this site, chr6:99593030, is a mutational hotspot

Comparative analysis of different culture conditions inducing naive pluripotency in human embryonic stem cells

Stem cells & developmental biolog

A novel duplication involving PRDM13 in a Turkish family supports its role in North Carolina macular dystrophy (NCMD/MCDR1)

Purpose: To clinically and molecularly investigate a new family with North Carolina macular dystrophy (NCMD) from Turkey, a previously unreported geographic origin for this phenotype. Methods: Clinical ophthalmic examinations, including fundus imaging and spectral domain-optical coherence tomography (SD-OCT), were performed on eight members of a two-generation non-consanguineous family from southern Turkey. Whole genome sequencing (WGS) was performed on two affected subjects, followed by variant filtering and copy number variant (CNV) analysis. Junction PCR and Sanger sequencing were used to confirm and characterize the duplication involving PRDM13 at the nucleotide level. The underlying mechanism was assessed with in silico analyses. Results: The proband presented with lifelong bilateral vision impairment and displayed large grade 3 coloboma-like central macular lesions. Five of her six children showed similar macular malformations, consistent with autosomal dominant NCMD. The severity grades in the six affected individuals from two generations are not evenly distributed. CNV analysis of WGS data of the two affected family members, followed by junction PCR and Sanger sequencing, revealed a novel 56.2 kb tandem duplication involving PRDM13 (chr6:99560265-99616492dup, hg38) at the MCDR1 locus. This duplication cosegregates with the NCMD phenotype in the five affected children. No other (likely) pathogenic variants in known inherited retinal disease genes were found in the WGS data. Bioinformatics analyses of the breakpoints suggest a replicative-based repair mechanism underlying the duplication. Conclusions: We report a novel tandem duplication involving the PRDM13 gene in a family with NCMD from a previously unreported geographic region. The duplication size is the smallest that has been reported thus far and may correlate with the particular phenotype

Correction to: 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

Multi-omics approach dissects cis-regulatory mechanisms underlying North Carolina macular dystrophy, a retinal enhanceropathy.

North Carolina macular dystrophy (NCMD) is a rare autosomal-dominant disease affecting macular development. The disease is caused by non-coding single-nucleotide variants (SNVs) in two hotspot regions near PRDM13 and by duplications in two distinct chromosomal loci, overlapping DNase I hypersensitive sites near either PRDM13 or IRX1. To unravel the mechanisms by which these variants cause disease, we first established a genome-wide multi-omics retinal database, RegRet. Integration of UMI-4C profiles we generated on adult human retina then allowed fine-mapping of the interactions of the PRDM13 and IRX1 promoters and the identification of eighteen candidate cis-regulatory elements (cCREs), the activity of which was investigated by luciferase and Xenopus enhancer assays. Next, luciferase assays showed that the non-coding SNVs located in the two hotspot regions of PRDM13 affect cCRE activity, including two NCMD-associated non-coding SNVs that we identified herein. Interestingly, the cCRE containing one of these SNVs was shown to interact with the PRDM13 promoter, demonstrated in vivo activity in Xenopus, and is active at the developmental stage when progenitor cells of the central retina exit mitosis, suggesting that this region is a PRDM13 enhancer. Finally, mining of single-cell transcriptional data of embryonic and adult retina revealed the highest expression of PRDM13 and IRX1 when amacrine cells start to synapse with retinal ganglion cells, supporting the hypothesis that altered PRDM13 or IRX1 expression impairs interactions between these cells during retinogenesis. Overall, this study provides insight into the cis-regulatory mechanisms of NCMD and supports that this condition is a retinal enhanceropathy

Multi-omics approach dissects cis-regulatory mechanisms underlying North Carolina macular dystrophy, a retinal enhanceropathy

North Carolina macular dystrophy (NCMD) is a rare autosomal-dominant disease affecting macular development. The disease is caused by non-coding single-nucleotide variants (SNVs) in two hotspot regions near PRDM13 and by duplications in two distinct chromosomal loci, overlapping DNase I hypersensitive sites near either PRDM13 or IRX1. To unravel the mechanisms by which these variants cause disease, we first established a genome-wide multi-omics retinal database, RegRet. Integration of UMI-4C profiles we generated on adult human retina then allowed fine-mapping of the interactions of the PRDM13 and IRX1 promoters and the identification of eighteen candidate cis-regulatory elements (cCREs), the activity of which was investigated by luciferase and Xenopus enhancer assays. Next, luciferase assays showed that the non-coding SNVs located in the two hotspot regions of PRDM13 affect cCRE activity, including two NCMD-associated non-coding SNVs that we identified herein. Interestingly, the cCRE containing one of these SNVs was shown to interact with the PRDM13 promoter, demonstrated in vivo activity in Xenopus, and is active at the developmental stage when progenitor cells of the central retina exit mitosis, suggesting that this region is a PRDM13 enhancer. Finally, mining of single-cell transcriptional data of embryonic and adult retina revealed the highest expression of PRDM13 and IRX1 when amacrine cells start to synapse with retinal ganglion cells, supporting the hypothesis that altered PRDM13 or IRX1 expression impairs interactions between these cells during retinogenesis. Overall, this study provides insight into the cis-regulatory mechanisms of NCMD and supports that this condition is a retinal enhanceropathy