Identification of circular RNAs hosted by the RPGR ORF15 genomic locus.

Mutations in the retina-specific isoform of the gene encoding retinitis pigmentosa GTPase regulator (RPGRorf15) cause X-linked retinitis pigmentosa, a severe and early onset inherited retinal degeneration. The underlying pathogenic mechanisms and variability in disease severity remain to be fully elucidated. The present study examines structural features of the ORF15 exonic region to provide new insights into the disease pathogenesis. Using canine and human RNA samples, we identified several novel RPGR ORF15-like linear RNA transcripts containing cryptic introns (exitrons) within the annotated exon ORF15. Furthermore, using outward-facing primers designed inside exitrons in the ORF15 exonic region, we found many of previously unidentified circular RNAs (circRNAs) that formed via back fusion of linear parts of the RPGRorf15 pre-mRNAs. These circRNAs (resistant to RNAse R treatment) were found in all studied cells and tissues. Notably, some circRNAs were present in cytoplasmic and polysomal RNA fractions. Although certain RPGR circRNAs may be cell type specific, we found some of the same circRNAs expressed in different cell types, suggesting similarities in their biogenesis and functions. Sequence analysis of RPGR circRNAs revealed several remarkable features, including identification of N6-methyladenosine (m6A) consensus sequence motifs and high prevalence of predictive microRNA binding sites pointing to the functional roles of these circRNAs. Our findings also illustrate the presence of non-canonical RPGR circRNA biogenesis pathways independent of the known back splicing mechanism. The obtained data on novel RPGR circRNAs further underline structural complexity of the RPGR ORF15 region and provide a potential molecular basis for the disease phenotypic heterogeneity.

Candidate Genetic Modifiers for RPGR Retinal Degeneration.

Purpose: To define genetic variants associated with variable severity of X-linked progressive retinal atrophy 1 (XLPRA1) caused by a five-nucleotide deletion in canine RPGR exon ORF15. Methods: A genome-wide association study (GWAS) was performed in XLPRA1 phenotype informative pedigree. Whole genome sequencing (WGS) was used for mutational analysis of genes within the candidate genomic region. Retinas of normal and mutant dogs were used for gene expression, gene structure, and RNA duplex analyses. Results: GWAS followed by haplotype phasing identified an approximately 4.6 Mb candidate genomic interval on CFA31 containing seven protein-coding genes expressed in retina (ROBO1, ROBO2, RBM11, NRIP1, HSPA13, SAMSN1, and USP25). Furthermore, we identified and characterized two novel lncRNAs, ROBO1-AS and ROBO2-AS, that display overlapping gene organization with axon guidance pathway genes ROBO1 and ROBO2, respectively, producing sense-antisense gene pairs. Notably, ROBO1-AS and ROBO2-AS act in cis to form lncRNA/mRNA duplexes with ROBO1 and ROBO2, respectively, suggesting important roles for these lncRNAs in the ROBO regulatory network. A subsequent WGS identified candidate genes within the genomic region on CFA31 that might be implicated in modifying severity of XLPRA1. This approach led to discovery of genetic variants in ROBO1, ROBO1-AS, ROBO2-AS, and USP25 that are strongly associated with the XLPRA1 moderate phenotype. Conclusions: The study provides new insights into the genetic basis of phenotypic variation in severity of RPGRorf15-associated retinal degeneration. Our findings suggest an important role for ROBO pathways in disease progression further expanding on our previously reported changes of ROBO1 expression in XLPRA1 retinas.

Critical Decrease in the Level of Axon Guidance Receptor ROBO1 in Rod Synaptic Terminals Is Followed by Axon Retraction.

Purpose: To define remodeling of photoreceptor synaptic terminals and second-order retinal neurons in canine X-linked progressive retinal atrophy 1 caused by a five-nucleotide deletion in the RPGR exon ORF15. Methods: Retinas of normal and mutant dogs were used for gene expression, Western blot, and immunohistochemistry. Cell-specific markers were used to examine disease-dependent retinal remodeling. Results: In mutant retinas, a number of rod axon terminals retract into the outer nuclear layer. This neuritic atrophy preceded significant loss of rods and was evident early in disease. Rod bipolar and horizontal cell processes were found to extend into the outer nuclear layer, where they seemed to form contacts with the spherules of rod photoreceptors. No ectopic rewiring was observed. Because cytoskeletal reorganization was previously shown to underlie photoreceptor axon retraction, we examined normal and mutant retinas for expression of axon guidance receptors ROBO1 and ROBO2, which are known to regulate actin cytoskeleton dynamics. We found that the overall expression of both ROBO1 and ROBO2 is retained at the same level in premature and fully developed normal retinas. However, analysis of predisease and early disease retinas identified markedly decreased levels of ROBO1 in rod spherules compared with controls. In contrast, no differences in ROBO1 signals were noted in cone pedicles in normal and mutant retinas, where ROBO1 levels remained similarly low. Conclusions: Depletion of ROBO1 in rod synaptic terminals correlates with the remodeling of axonal and dendritic processes in the outer retina of dogs with X-linked progressive retinal atrophy 1 and may play a role in the retraction of rod axons.

Mutation-independent rhodopsin gene therapy by knockdown and replacement with a single AAV vector.

Inherited retinal degenerations are caused by mutations in >250 genes that affect photoreceptor cells or the retinal pigment epithelium and result in vision loss. For autosomal recessive and X-linked retinal degenerations, significant progress has been achieved in the field of gene therapy as evidenced by the growing number of clinical trials and the recent commercialization of the first gene therapy for a form of congenital blindness. However, despite significant efforts to develop a treatment for the most common form of autosomal dominant retinitis pigmentosa (adRP) caused by >150 mutations in the rhodopsin (RHO) gene, translation to the clinic has stalled. Here, we identified a highly efficient shRNA that targets human (and canine) RHO in a mutation-independent manner. In a single adeno-associated viral (AAV) vector we combined this shRNA with a human RHO replacement cDNA made resistant to RNA interference and tested this construct in a naturally occurring canine model of RHO-adRP. Subretinal vector injections led to nearly complete suppression of endogenous canine RHO RNA, while the human RHO replacement cDNA resulted in up to 30% of normal RHO protein levels. Noninvasive retinal imaging showed photoreceptors in treated areas were completely protected from retinal degeneration. Histopathology confirmed retention of normal photoreceptor structure and RHO expression in rod outer segments. Long-term (>8 mo) follow-up by retinal imaging and electroretinography indicated stable structural and functional preservation. The efficacy of this gene therapy in a clinically relevant large-animal model paves the way for treating patients with RHO-adRP.

Strong upregulation of inflammatory genes accompanies photoreceptor demise in canine models of retinal degeneration.

We have analyzed the complex pattern of the inflammatory response in early-onset canine models of human retinitis pigmentosa, rcd1, xlpra2 and erd, as well as late-onset xlpra1, in comparative manner. The time course of immune response genes and proteins expression was examined along the timeline of photoreceptors degeneration. Gene expression analysis of the early-onset models prior to and after the peak of photoreceptors death identified the involvement of multiple immune response genes including those encoding constituents of the NLRP3 inflammasome, its substrates, pro-IL1B, pro-IL18, and common components of IL1B, IL18 and TLR4 pathways. Out of two activated caspase-1 cleavage products, IL1B and IL18, only IL1B was detected in rcd1 and xlpra2 while precursor IL18 remained unprocessed in the same protein extract highlighting prominence of IL1B pathway. An overall immune response was most prominent in rcd1 followed by xlpra2 and least prominent in erd. Noticeably, in rcd1 and xlpra2, but not in erd, early induction of the immune response was accompanied by sustained intraretinal migration and activation of retinal microglia. Lastly, delayed activation of the anti-inflammatory factors in all early-onset models was insufficient to counterbalance rapidly progressing inflammation. In contrast to early-onset models, in late-onset xlpra1 retinas a subset of the pro-inflammatory genes was highly upregulated long before any disease-related structural changes occurred, but was counterbalanced by an adequate anti-inflammatory response. Results point out to upregulated immune response accompanying disease progression in animal models of retinal degeneration, and to potential benefits of early anti-inflammatory therapy.

Molecular studies of phenotype variation in canine RPGR-XLPRA1.

PURPOSE: Canine X-linked progressive retinal atrophy 1 (XLPRA1) caused by a mutation in retinitis pigmentosa (RP) GTPase regulator (RPGR) exon ORF15 showed significant variability in disease onset in a colony of dogs that all inherited the same mutant X chromosome. Defective protein trafficking has been detected in XLPRA1 before any discernible degeneration of the photoreceptors. We hypothesized that the severity of the photoreceptor degeneration in affected dogs may be associated with defects in genes involved in ciliary trafficking. To this end, we examined six genes as potential disease modifiers. We also examined the expression levels of 24 genes involved in ciliary trafficking (seven), visual pathway (five), neuronal maintenance genes (six), and cellular stress response (six) to evaluate their possible involvement in early stages of the disease. METHODS: Samples from a pedigree derived from a single XLPRA1-affected male dog outcrossed to unrelated healthy mix-bred or purebred females were used for immunohistochemistry (IHC), western blot, mutational and haplotype analysis, and gene expression (GE). Cell-specific markers were used to examine retinal remodeling in the disease. Single nucleotide polymorphisms (SNPs) spanning the entire RPGR interacting and protein trafficking genes (RAB8A, RPGRIP1L, CEP290, CC2D2A, DFNB31, and RAB11B) were genotyped in the pedigree. Quantitative real-time PCR (qRT-PCR) was used to examine the expression of a total of 24 genes, including the six genes listed. RESULTS: Examination of cryosections from XLPRA1-affected animals of similar age (3-4 years) with different disease severity phenotype revealed mislocalization of opsins and upregulation of the Müller cell gliosis marker GFAP. Four to ten haplotypes per gene were identified in RAB8A, RPGRIP1L, CEP290, CC2D2A, DFNB31, and RAB11B for further assessment as potential genetic modifiers of XLPRA1. No correlation was found between the haplotypes and disease severity. During mutational analysis, several new variants, including a single intronic mutation in RAB8A and three mutations in exon 3 of DFNB31 were described (c.970G>A (V324I), c.978T>C (G326=), and c.985G>A (A329T)). Expression analysis of stress response genes in 16-week-old predisease XLPRA1 retinas revealed upregulation of GFAP but not HSPA5, DDIT3, HSPA4, HSP90B1, or HIF1A. Western blot analysis confirmed GFAP upregulation. In the same predisease group, no significant differences were found in the expression of 18 selected genes (RHO, OPN1LW, OPN1MW, RLBP1, RPGRORF15, RAB8A, RPGRIP1L, CEP290, CC2D2A, DFNB31, RAB11B, CRX, RCVRN, PVALB, CALB1, FGFR1, NTRK2, and NTRK3) involved in neuronal function. CONCLUSIONS: Lack of association between haplotypes of RAB8A, RPGRIP1L, CEP290, CC2D2A, DFNB31, and RAB11B and the disease phenotype suggests that these genes are not genetic modifiers of XLPRA1. Upregulation of GFAP, an established indicator of the Müller cell gliosis, manifests as an important early feature of the disease.