Mutations in retinitis pigmentosa GTPase regulator (RPGR) cause severe retinal ciliopathy, X-linked retinitis pigmentosa. Although two major alternatively spliced isoforms, RPGRex1-19 and RPGRORF15, are expressed, the relative importance of these isoforms in disease pathogenesis is unclear. Here, we analyzed fibroblast samples from eight patients and found that all of them form longer cilia than normal controls, albeit to different degrees. Although all mutant RPGRORF15 messenger RNAs (mRNAs) are unstable, their steady-state levels were similar or higher than those in the control cells, suggesting there may be increased transcription. Three of the fibroblasts that had higher levels of mutant RPGRORF15 mRNA also exhibited significantly higher levels of RPGRex1-19 mRNA. Four samples with unaltered RPGRex1-19 levels carried mutations in RPGRORF15 that resulted in this isoform being relatively less stable. Thus, in all cases, the RPGRex1-19/RPGRORF15 isoform ratio was increased, and this was highly correlative to the cilia extension defect. Moreover, overexpression of RPGRex1-19 (mimicking the increase in RPGRex1-19 to RPGRORF15 isoform ratio) or RPGRORF15 (mimicking reduction of the ratio) resulted in significantly longer or shorter cilia, respectively. Notably, the cilia length defect appears to be attributable to both the loss of the wild-type RPGRORF15 protein and to the higher levels of the RPGRex1-19 isoform, indicating that the observed defect is due to the altered isoform ratios. These results suggest that maintaining the optimal RPGRex1-9 to RPGRORF15 ratio is critical for cilia growth and that designing strategies that focus on the best ways to restore the RPGRex1-19/RPGRORF15 ratio may lead to better therapeutic outcomes.
Mutations in retinitis pigmentosa GTPase regulator (RPGR) cause severe retinal ciliopathy, X-linked retinitis pigmentosa. Although two major alternatively spliced isoforms, RPGRex1-19 and RPGRORF15, are expressed, the relative importance of these isoforms in disease pathogenesis is unclear. Here, we analyzed fibroblast samples from eight patients and found that all of them form longer cilia than normal controls, albeit to different degrees. Although all mutant RPGRORF15 messenger RNAs (mRNAs) are unstable, their steady-state levels were similar or higher than those in the control cells, suggesting there may be increased transcription. Three of the fibroblasts that had higher levels of mutant RPGRORF15 mRNA also exhibited significantly higher levels of RPGRex1-19 mRNA. Four samples with unaltered RPGRex1-19 levels carried mutations in RPGRORF15 that resulted in this isoform being relatively less stable. Thus, in all cases, the RPGRex1-19/RPGRORF15 isoform ratio was increased, and this was highly correlative to the cilia extension defect. Moreover, overexpression of RPGRex1-19 (mimicking the increase in RPGRex1-19 to RPGRORF15 isoform ratio) or RPGRORF15 (mimicking reduction of the ratio) resulted in significantly longer or shorter cilia, respectively. Notably, the cilia length defect appears to be attributable to both the loss of the wild-type RPGRORF15 protein and to the higher levels of the RPGRex1-19 isoform, indicating that the observed defect is due to the altered isoform ratios. These results suggest that maintaining the optimal RPGRex1-9 to RPGRORF15 ratio is critical for cilia growth and that designing strategies that focus on the best ways to restore the RPGRex1-19/RPGRORF15 ratio may lead to better therapeutic outcomes.
Primary cilia are evolutionarily conserved microtubule-based non-motile protrusions of the plasma membrane of most vertebrate cell types (1). They are nucleated from the basal bodies (mother centrioles) that dock at the apical plasma membrane in quiescent cells. The proximal region of the cilia, also called transition zone, acts as a gate. The microtubules further extend to form a distal axoneme surrounded by the ciliary membrane. The primary cilium regulates critical developmental and homeostatic signaling cascades by targeting several integral and peripheral membrane proteins as well as soluble moieties along the axoneme via a conserved mechanism termed intraflagellar transport (2,3). The trafficking is also regulated by macromolecular complexes of several proteins including retinitis pigmentosa GTPase regulator (RPGR), centrosomal protein 290 kDa (CEP290), centrosomal protein 110 kDa (CP110), Usher Syndrome proteins and Bardet Biedl Syndrome proteins (4,5). Mutations in these proteins are associated with severe developmental and homeostatic disorders collectively termed ciliopathies (6). These include Meckel–Gruber Syndrome, Joubert Syndrome, Senior–Loken Syndrome, Usher Syndrome and some forms of non-syndromic photoreceptor degenerations.Photoreceptors (rods and cones) have a modified sensory cilium in the form of the light-sensing outer segment. The photoreceptors are one of the highly metabolically active cell types in the body. They depend upon stringently controlled proteins and lipid trafficking from their site of synthesis in the inner segment to the outer segment for light detection (7). Perturbations in the trafficking machinery due to functional deficits in ciliary proteins lead to severe photoreceptor degeneration (2,3,8). One such ciliary trafficking regulator is RPGR.Mutations in RPGR are associated with X-linked forms of retinitis pigmentosa (RP), a debilitating blindness characterized by severe and progressive photoreceptor degeneration (9–12). RPGR mutations account for the majority (~70%) of X-linked retinitis pigmentosa (XLRP) cases and 12–15% of simplex RP, with no family history (12–14). In addition, RPGR mutations are associated with cone–rod degeneration, macular dystrophy and extraocular presentations, such as sperm abnormalities, hearing defects and primary cilia dyskinesia (15–19). Overall, RPGR mutations are a leading cause of blindness accounting for almost 20% of all RP cases.The RPGR gene encodes two major alternatively spliced isoforms: RPGR (19 exons) and RPGR (terminating in intron 15; exon ORF15) (20–23). The terminal exon 19 of the RPGRex1-19 isoform carries an isoprenylation signal involved in its localization to cilia (24,25). Exon ORF15 of the RPGR isoform contains purine-rich repeats that encode for a glutamic acid and glycine-rich domain (26). As majority of RPGR mutations are located in exon ORF15, the RPGR isoform has been considered the major disease-associated isoform in humans (26). The great majority of exon ORF15 mutations are predicted to prematurely truncate the RPGRORF15 protein.Both RPGR isoforms are expressed throughout the body, but are particularly abundant in the retina. Previous studies have identified mouse and canine models of RPGR mutations and have suggested the involvement of RPGR in regulating ciliary trafficking by acting as a gatekeeper at the transition zone (27–31). Loss of RPGR alters protein trafficking to the photoreceptor outer segment, but permits its initial generation. RPGR depletion alters cilia length in cultured cell lines and fibroblasts (32). However, the effect on the ciliary defects of different RPGR mutations is unclear.To investigate the effect of the human disease-causing mutations in RPGR on its expression and on cilia formation, we analyzed RPGR isoform expression, cilia formation and ciliary protein trafficking in patient-derived cells. Our studies unravel a new potential pathogenic mechanism of RPGRORF15 disease, which may assist in developing new therapies for these conditions.
Results
Ciliary defects in RPGRORF15-mutant fibroblasts
Given that RPGR regulates cilia length (32), we examined cilia formation in the patient and control fibroblasts. We generated primary fibroblasts from skin biopsy specimens from eight male RPGRpatients (RPGR-1, RPGR-2, RPGR-3, RPGR-4, RPGR-5, RPGR-6, RPGR-7 and RPGR-8) with frameshift and truncation mutations in exon ORF15 (Fig. 1A and B). The patients were all from families previously confirmed to have X-linked RP due to RPGR mutations (33,34); their ages at skin biopsy ranged from 21 to 66 years. The patients had ophthalmoscopic features of pigmentary retinopathy, abnormal visual fields and reduced or non-detectable electroretinograms (33,34). As the control samples from the same family members could not be obtained, we included six random asymptomatic control samples without RPGR mutations (CTL-1, CTL-2, CTL-3, CTL-4, CTL-5 and CTL-6).
Figure 1
RPGR mutations. Schematic representation (A) and predicted effect of RPGR mutations (B) in XLRP patients. *: optic cups of RPGR-3 were generated.
RPGR mutations. Schematic representation (A) and predicted effect of RPGR mutations (B) in XLRP patients. *: optic cups of RPGR-3 were generated.Although the control and RPGR-2 fibroblasts grew cilia in the range of 1–1.5 μm, we detected two groups of mutations with longer cilia: cilia were elongated to a range of 2.5–3.5 μm range in the RPGR-1, RPGR-3, RPGR-4 and RPGR-5 fibroblasts and to a lesser extent (1.8–2.5 μm) in RPGR-6, RPGR-7 and RPGR-8 (Fig. 2A and B). As predicted, we did not detect an effect on the number of ciliated cells.
Figure 2
Ciliary defects in fibroblasts. (A) Control (CTL-1 and CTL-2) and mutant fibroblasts were stained for acetylated α-tubulin (red, ciliary marker) and 4′,6-diamidino-2-phenylindole (DAPI) (blue, nucleus). Scale: 5 μm. (B) The cilia length was quantified using the ImageJ software (n > 100). The CTL data represents the average cilia length of all control fibroblasts. ****P < 0.0001; ***P < 0.002; **P < 0.01; *P < 0.05.
Ciliary defects in fibroblasts. (A) Control (CTL-1 and CTL-2) and mutant fibroblasts were stained for acetylated α-tubulin (red, ciliary marker) and 4′,6-diamidino-2-phenylindole (DAPI) (blue, nucleus). Scale: 5 μm. (B) The cilia length was quantified using the ImageJ software (n > 100). The CTL data represents the average cilia length of all control fibroblasts. ****P < 0.0001; ***P < 0.002; **P < 0.01; *P < 0.05.
Correlation of cilia length to RPGR isoform levels in fibroblasts
We next asked whether the observed ciliary phenotype was associated with an effect of exon ORF15 mutations on RPGR expression. To this end, we performed quantitative reverse transcriptase-polymerase chain reaction (RT-PCR) analysis of the two major isoforms. We found significantly higher levels (8–10-fold) of the RPGR isoform in RPGR-3, RPGR-4 and RPGR-5 fibroblasts samples when compared with the other RPGR mutants or the control fibroblasts (Fig. 3A). We also detected higher levels (2–3-fold) of the mutant RPGR isoform in RPGR-3, RPGR-4 and RPGR-5. Interestingly, the level of the mutant RPGR in RPGR-1 was significantly reduced when compared with the controls (Fig. 3B).
Figure 3
RPGR isoform levels in patient-derived fibroblasts. (A, B) RPGR (A) or RPGR levels in the mutant fibroblasts relative to controls (CTL) were analyzed. (C) The ratio of RPGR to RPGR mRNA levels was determined. RPGR-1, RPGR-3, RPGR-4 and RPGR-5 showed significantly higher ratio when compared with control (CTL) cells. ns, not significant. Data are mean ± SD from five independent experiments. (D) Positive correlation of the cilium length with the RPGR/RPGR ratio was observed. Each dot represents average cilium length of >100 cilia compared with the RPGR/RPGR ratio from three independent experiments. Black line shows a linear fit through the data. Inset displays the coefficient of determination (R2) and significance (P-value).
RPGR isoform levels in patient-derived fibroblasts. (A, B) RPGR (A) or RPGR levels in the mutant fibroblasts relative to controls (CTL) were analyzed. (C) The ratio of RPGR to RPGR mRNA levels was determined. RPGR-1, RPGR-3, RPGR-4 and RPGR-5 showed significantly higher ratio when compared with control (CTL) cells. ns, not significant. Data are mean ± SD from five independent experiments. (D) Positive correlation of the cilium length with the RPGR/RPGR ratio was observed. Each dot represents average cilium length of >100 cilia compared with the RPGR/RPGR ratio from three independent experiments. Black line shows a linear fit through the data. Inset displays the coefficient of determination (R2) and significance (P-value).Given that RPGR-1, RPGR-3, RPGR-4 and RPGR-5 also formed relatively longer cilia, we tested whether there was a correlation between the RPGR/RPGR ratio and the cilia length. Although the RPGR-2, RPGR-6, RPGR-7, RPGR-8 and control fibroblasts showed an RPGR/RPGR ratio of ~10-fold, in RPGR-1, RPGR-3, RPGR-4 and RPGR-5 fibroblasts, there was an increased ratio of between 18- and 35-fold (Fig. 3C). Furthermore, we found that the length of the cilia was highly correlative to RPGR ratio, with the coefficient of determination (R2) of ~83% (Fig. 3D). Using RPGR isoform-specific antibodies (Supplementary Fig. S1A), we also detected a significant increase (18–20-fold) in the ratio of the RPGRex1-19/RPGRORF15 protein levels in RPGR-1, RPGR-3, RPGR-4 and RPGR-5 when compared with other RPGR mutations or controls (Supplementary Fig. S1B and C). We then checked whether such phenomenon of increased RPGR/RPGR ratio is also conserved in retinal microenvironment. Using retinal organoids derived from reprogrammed control and RPGR-3 human-induced pluripotent stem cells (iPSCs; Supplementary Fig. S2A and B) (West et al. 2020, manuscript submitted), we found significant (~2.8-fold) increase in the RPGR ratio in the mutant sample as compared with the controls (Supplementary Fig. S2C).These results suggest that RPGRORF15 mutations tested in this study affect not only mutant RPGR expression levels but also the RPGR isoform levels, albeit to different extents. Although in fibroblast samples from all RPGRpatients there was a positive correlation between the RPGR ratio and cilia length, in RPGR-1, RPGR-3, RPGR-4 and RPGR-5 samples, there were distinct effects on the RPGR messenger RNA (mRNA) levels. The RPGR-1 mutation results in significantly reduced mutant RPGR expression, but does not affect the RPGR levels. On the other hand, although the mutant RPGR is expressed at relatively higher levels in RPGR-3, RPGR-4 and RPGR-5 fibroblast samples, they exhibit significantly higher levels of RPGR. To further understand the ciliary defects, we focused on RPGR-1, RPGR-3, RPGR-4 and RPGR-5 samples. We also used RPGR-2 as it showed similar RPGR expression as the control fibroblasts.
Localization of cilia length regulators in RPGR-mutant fibroblasts
To investigate the ciliary defect in the fibroblasts with elongated cilia, we assessed the distribution of ciliary length regulators CEP290 and distribution of glutamylated tubulin. CEP290 mutations are a frequent cause of Leber congenital amaurosis and syndromic ciliopathies (35–37). Immunofluorescence staining of RPGR-1, RPGR-3, RPGR-4 and RPGR-5 with a previously reported CEP290 antibody demonstrated ~50% reduction in the levels of CEP290 at the basal body (determined by the co-localization with basal body marker γ-tubulin) and concomitant increase in the cytosolic CEP290 distribution when compared with RPGR-2 or control fibroblasts (Fig. 4A and B). However, the number of CEP290-positive cells was not affected in the mutant as compared with the control fibroblasts (Fig. 4C). These results are consistent with previous studies that demonstrated longer cilia in CEP290-mutant patient fibroblasts and in Cep290-mutant mice (38,39) and suggest an association between the observed cilia length differences and CEP290 levels in the RPGR-mutant fibroblasts.
Figure 4
Distribution of ciliary markers. CEP290 (A–C) and GT335 (D–F) staining in the indicated fibroblasts was assessed. (A) The cells were stained with CEP290 (green) and γ-tubulin (basal body marker, red) antibodies and the ratio of the CEP290-positive signal in the cytoplasm to that at the basal body (B) and the number of CEP290-positive cells (C) were quantified using ImageJ (n > 100). (D) The cells were stained with GT335 (green, glutamylated tubulin) and ARL13B (red, cilia marker). (E) The distribution of GT335 (green signal intensity) along the length of the cilium when compared with ARL13B (red intensity) (distance, μm) was quantified using IMARIS software (Oxford Instruments). a.u., arbitrary units.
Distribution of ciliary markers. CEP290 (A–C) and GT335 (D–F) staining in the indicated fibroblasts was assessed. (A) The cells were stained with CEP290 (green) and γ-tubulin (basal body marker, red) antibodies and the ratio of the CEP290-positive signal in the cytoplasm to that at the basal body (B) and the number of CEP290-positive cells (C) were quantified using ImageJ (n > 100). (D) The cells were stained with GT335 (green, glutamylated tubulin) and ARL13B (red, cilia marker). (E) The distribution of GT335 (green signal intensity) along the length of the cilium when compared with ARL13B (red intensity) (distance, μm) was quantified using IMARIS software (Oxford Instruments). a.u., arbitrary units.Tubulin polyglutamylation is a modification of the stable ciliary microtubules, which can be detected using the GT335 antibody (40,41). We found that the distribution of GT335 was altered in the RPGR-mutant fibroblasts. Immunofluorescence analysis revealed that GT335 immunoreactivity was restricted to the proximal region of the cilia in RPGR-1, RPGR-3, RPGR-4 and RPGR-5 (within ~1.0 μm) (Fig. 4D and E) when compared with RPGR-2 and the controls, in which GT335 staining is distributed across the whole cilium (Fig. 4F).
Expression of wild-type RPGRORF15 in mutant fibroblasts
We then sought to determine whether the observed cilia length defects are due to the lack of the wild-type RPGRORF15 protein. Transfection of the wild-type human green fluorescent protein (GFP)-RPGRORF15-encoding complementary DNA (cDNA) into RPGR-3, RPGR-4 and RPGR-5 did not completely reverse the elongated cilia phenotype when compared with cells transfected with only GFP. Although the RPGR-3, RPGR-4 and RPGR-5 fibroblasts transfected with GFP alone formed longer cilia than the CTL (control) fibroblasts, the cilia length reduced from 2.2–2.4 μm to 1.3–1.5 μm after transfection with GFP-RPGRORF15 (Fig. 5A and B). RPGR-1, which had relatively modest (~1.8 μm) increase in the cilia length and RPGR/RPGR ratio showed almost complete reversal (~1 μm) of the elongated cilia phenotype upon GFP-RPGRORF15 expression. These data strongly suggest that the normal RPGRORF15 protein function partly contributes to the observed cilia length changes. Samples with relatively higher RPGR/RPGR ratio did not show complete reversal of the cilia length phenotype. Further support of this hypothesis comes from the observation that expressing RPGRORF15 in the control or RPGR-2 (~20% longer cilia) cells significantly reduced the cilia length 10–20% less than the average lengths of the GFP-expressing control cells. All cells transfected with RPGRORF15-encoding cDNA overexpressed RPGR to similar levels (4–5-fold higher) when compared with GFP-expressing control cells (Supplementary Fig. S3A) and did not alter the RPGR levels (Supplementary Fig. S3B).
Figure 5
Effect of wild-type RPGRORF15 overexpression on the cilia defect. (A) The fibroblasts were transiently transfected with plasmids-encoding GFP or GFP-RPGRORF15 and processed for staining with the GFP (green) or acetylated α-tubulin (red, ciliary marker). Arrows point to the ciliated cells marked in the inset in green channel. The nuclei are stained with DAPI (blue). (B) The cilia length was quantified using ImageJ (n > 100). The graph shows the cilia length relative to the length of the control (CTL) fibroblasts transfected with only GFP-encoding cDNA. The statistical significance of the cilia length in the GFP or GFP-RPGRORF15-transfected mutant fibroblasts cells was determined relative to the GFP-transfected control fibroblasts. **P < 0.01; ***P < 0.0002; ****P < 0.0001; ns, not significant. Note that the GFP-RPGRORF15-transfected RPGR-2 and CTL cells have significantly shorter cilia when compared with the GFP-transfected counterparts.
Effect of wild-type RPGRORF15 overexpression on the cilia defect. (A) The fibroblasts were transiently transfected with plasmids-encoding GFP or GFP-RPGRORF15 and processed for staining with the GFP (green) or acetylated α-tubulin (red, ciliary marker). Arrows point to the ciliated cells marked in the inset in green channel. The nuclei are stained with DAPI (blue). (B) The cilia length was quantified using ImageJ (n > 100). The graph shows the cilia length relative to the length of the control (CTL) fibroblasts transfected with only GFP-encoding cDNA. The statistical significance of the cilia length in the GFP or GFP-RPGRORF15-transfected mutant fibroblasts cells was determined relative to the GFP-transfected control fibroblasts. **P < 0.01; ***P < 0.0002; ****P < 0.0001; ns, not significant. Note that the GFP-RPGRORF15-transfected RPGR-2 and CTL cells have significantly shorter cilia when compared with the GFP-transfected counterparts.
Some RPGR-mutant isoforms are less stable
Varying mRNA levels of the RPGR isoforms could be due to the mutations affecting the purine-rich repeat region of the genomic DNA or RPGR mRNA stability. As some RPGR mutations (RPGR-2, RPGR-6, RPGR-7 and RPGR-8) did not affect the isoform ratio, we concluded that the DNA architecture in the region of those mutations is not sufficient to cause the observed changes. We then tested whether the effect on RPGR isoform ratio is due to changes in its stability. We performed transcription inhibition assays using actinomycin D to assess the half-lives of the RPGR isoforms. Our studies revealed that the mutant RPGR levels rapidly decrease to ~10% of the total levels within the first hour of actinomycin D treatment and are at minimal levels by 4 h in RPGR-3, RPGR-4 and RPGR-5 samples. Their levels also remained significantly lower than the control at all time points. However, the rate of decline of the mutant RPGR levels in RPGR-1 and RPGR-2 samples was overall similar to that of the control RPGR mRNA. There was, however, a slightly faster decline of the mutant RPGR levels in RPGR-1 and RPGR-2 at 6 and 8 h of actinomycin D treatment when compared with the control (Fig. 6A). As predicted, we did not detect a significant difference in the decay of the RPGR mRNA in the patient fibroblasts relative to the controls (Fig. 6B). Based on these studies, we found that the half-life of the mutant RPGR in RPGR-1 and RPGR-2 after actinomycin D treatment relative to vehicle treatment was similar to that in control fibroblasts (~2.3 and 2.2 h for RPGR-1 and RPGR-2, respectively, and 1.6 h for control) (Fig. 4C), whereas the half-life of the mutant RPGR in RPGR-3, RPGR-4 and RPGR-5 was 0.67, 0.51 and 0.28 h, respectively. Linear regression analysis further validated that the decay kinetics of mutant RPGR in RPGR-3, RPGR-4 and RPGR-5 does not follow a linear correlation (Fig. 6C).
Figure 6
Mutant mRNA decay in RPGR fibroblasts. Quantitative PCR analysis of RPGR (A) or RPGR (B) isoform was performed in fibroblasts treated with actinomycin D for indicated time. The levels of the mutant isoforms were calculated relative to the levels in the control (CTL) fibroblasts at the same time points. Statistical significance was determined by assessing the mRNA levels at each time point relative to that of the control (CTL) at the same time point. ****P < 0.0001; **P < 0.01; *P < 0.05. No significant changes in the degradation of RPGR mRNA were detected at the same time points (B). (C) The half-life and linear regression of the mutant RPGR and RPGR mRNAs was calculated by determining the time (in hours) taken by the isoforms to reduce to half the original levels.
Mutant mRNA decay in RPGR fibroblasts. Quantitative PCR analysis of RPGR (A) or RPGR (B) isoform was performed in fibroblasts treated with actinomycin D for indicated time. The levels of the mutant isoforms were calculated relative to the levels in the control (CTL) fibroblasts at the same time points. Statistical significance was determined by assessing the mRNA levels at each time point relative to that of the control (CTL) at the same time point. ****P < 0.0001; **P < 0.01; *P < 0.05. No significant changes in the degradation of RPGR mRNA were detected at the same time points (B). (C) The half-life and linear regression of the mutant RPGR and RPGR mRNAs was calculated by determining the time (in hours) taken by the isoforms to reduce to half the original levels.
RPGR knockdown in mutant fibroblasts
If RPGRORF15 only partly reversed the cilia length, the increase in the RPGR levels could also be relevant to the observed cilia elongation in the RPGR-mutant fibroblasts. To test this, we overexpressed RPGR in human telomerase reverse transcriptase-retinal pigmented epithelium (hTERT-RPE1) cells by transient transfection. Our analysis revealed 3–4-fold longer cilia in the RPGR overexpressing cells as compared with the control cells (Supplementary Fig. S4A). Thus, we hypothesized that reducing the RPGR levels will alleviate the elongated cilia phenotype. To this end, we performed small interfering ribonucleic acid (siRNA)-mediated depletion of RPGRex1-19. We designed three siRNAs (siRNA-1, siRNA-2 and siRNA-3) against RPGR. All siRNAs reduced the expression of RPGRex1-19, with maximal inhibition observed with siRNA-3 (Supplementary Fig. S4B), when compared with the levels expressed in cells transfected with a scrambled sequence control. Moreover, the siRNAs did not affect RPGRORF15 levels in the fibroblasts (Supplementary Fig. S4C). Transfection of RPGR-siRNAs followed by immunofluorescence analysis in RPGR-3 and RPGR-5 revealed an almost complete reversal of the cilia elongation phenotype with the cilia length comparable with that of the scrambled siRNA-treated control fibroblasts. RPGR-4 showed ~90% reduction in cilia length compared with scrambled siRNA-treated cells (Fig. 7A and B).
Figure 7
Effect of RPGR knockdown on the cilia defect. (A) The fibroblasts were transiently transfected with RPGR-siRNA-3 or scrambled siRNA followed by GFP (green, inset) and acetylated α-tubulin staining. (B) The cilia length was quantified represented as the length relative to the scrambled siRNA-treated control (CTL) cells. Scale: 5 μm; ns, not significant.
Effect of RPGR knockdown on the cilia defect. (A) The fibroblasts were transiently transfected with RPGR-siRNA-3 or scrambled siRNA followed by GFP (green, inset) and acetylated α-tubulin staining. (B) The cilia length was quantified represented as the length relative to the scrambled siRNA-treated control (CTL) cells. Scale: 5 μm; ns, not significant.
Discussion
Although RPGR gene augmentation therapy has been successful in naturally occurring canine mutations and murine models (42,43) and is currently in human clinical trials, efficacy of the treatment across the large spectrum of human exon ORF15 mutations is currently not known. We show that all eight RPGR mutations tested here cause cilia length defects of varying extents, lead to mutant mRNA instability and four mutations lead to a significant increase in the RPGR/RPGR ratio. Our results indicate that three out of four exon ORF15 mutations that result in an increase in the RPGR mRNA ratio have shorter half-life due to relatively faster mutant mRNA decay, yet have relatively higher steady-state levels of the mutant RPGR mRNA and protein. The precise mechanism of the transcriptional complexity observed in this study is not clear. Nonetheless, these results point to transcriptional adaptation in which the unstable mRNA transcripts can result in a feedback response to increase the levels of the protein. RPGR transcriptional adaptation could increase the RPGR protein levels; however, at the same time also increase RPGR levels, as observed in this study. Such transcriptional adaptation may result from the need to compensate for the decaying mutant RPGR transcripts (44–46); however, its relevance to retinal disease remains to be determined. We also found that the expression of wild-type RPGRORF15 did not affect the elevated RPGR levels in fibroblasts. This suggests that the mechanism of increase in the RPGR levels is not due to the absence of the wild-type RPGRORF15 and implicates the expression of the mutant ORF15 isoforms. The observed effects on cilia growth, on the other hand, seem to be due in part to the lack of the RPGRORF15 protein and to the relatively increased expression of RPGRex1-19. This is because RPGRORF15 overexpression or RPGRex1-19 knockdown were able to partly reverse the cilia defect, with RPGR knockdown showing better reversal of the cilia length. These results suggest that the relative levels of both RPGR isoforms are critical for optimal cilia growth. Further support of this hypothesis comes from our observation that the overexpression of RPGR leads to longer cilia while that of RPGR results in shorter cilia. In addition, previous studies showed that overexpression of Rpgr or RPGR in mice results in severe retinal degeneration (47,48).What is the mechanism of cilia elongation in RPGR-mutant samples? We found that the imbalanced RPGR ratio subsequently results in higher overall levels of the RPGR proteins. This could lead to increased context-dependent function of RPGR-containing protein complexes at the cilia and alteration of the gatekeeping function at the transition zone. We previously showed that RPGR regulates the composition of the ciliary outer segments of photoreceptors by affecting its ability to maintain the barrier at the transition zone (30). RPGRex1-19 is also involved in regulating ciliary cargo trafficking in cultured cells as well as mouse photoreceptors (24). Moreover, both RPGR isoforms interact genetically and physically with CEP290 and this interaction is critical for photoreceptor outer segment viability and function (24). Hence, RPGR–CEP290 interaction could take place in overlapping as well as distinct complexes consisting of the different isoforms. We, therefore, propose that RPGR mutations that perturb ratio of the two isoforms and subsequently the integrity or recruitment of the multiprotein complexes result in altered ciliary trafficking. Further studies are underway to delineate the mechanism of such effects.Mutations in retinal ciliary proteins have been mostly associated with reduced cilia length and photoreceptor outer segment degeneration. Knockdown of rpgr in zebrafish or in cultured human RPE cell lines results in shorter cilia (49–52). Although the fibroblasts derived from RPGRpatients exhibited longer cilia, longer photoreceptor outer segments were not observed in iPSC-derived optic cups or in RPGRpatients. These results could reflect cell-type specific differences in the utilization of distinct ciliary regulation pathways. Previous clinical findings and animal model studies also showed that RPGR mutations do not prevent the initial generation of the photoreceptor outer segments, but that they undergo progressive degeneration (28,29,31). Of note, context-dependent cilia length regulation and outer segment defects have been reported with CEP290-patient fibroblasts and CP110-mouse fibroblasts (39,53).The RPGR-2, RPGR-3, RPGR-4 and RPGR-5 are similar in nature (small frameshifts), they produce profoundly different effects. These could partly be due to the location of these mutations: while the RPGR-2 mutation is located towards the 3′-region of the non-repetitive exon 15, the RPGR-3, RPGR-4 and RPGR-5 mutations lie in the 5′-region of the repetitive purine-rich region. The repetitive region of these mutations may be involved in the recognition or recruitment of mRNA decay factors or splicing effectors. On the other hand, the RPGR-1 mutation, which lies within exon 15 does not show increased RPGR levels, yet the steady-state levels of the mutant RPGR mRNA are significantly reduced. As the mRNA stability profile of this mutant is similar to that of the control, there may be a different mechanism underlying reduced mutant RPGR levels. However, the resultant increase in the RPGR ratio indicates the similar cilia extension defect as RPGR-3, RPGR-4 and RPGR-5.The partial effect of overexpression of RPGR in the mutant fibroblasts on the rescue of the ciliary phenotype potentially also raises questions that they may be relevant to the ongoing clinical trials that utilize RPGR gene augmentation strategy. These studies are based upon the rationale that the RPGR mutations result in a complete or partial loss of function. Although our study was in progress, first results of a phase 1/2 study reported encouraging results with 7 out of 18 patients showing early improvement in visual function after RPGR gene augmentation therapy (54). However, at least one of the patients with the mutation c.2405_2406delAG (RPGR-3) did not reveal reliable clinical improvement. Our results predict that this mutation would cause a significant disruption in the RPGR isoform ratio and abnormally longer cilia that might not be completely reversed by RPGR overexpression. This study underscores the complexity of RPGR gene defects and suggests that gene therapy strategies may need to focus on ways to restore the appropriate RPGR ratio, taking into account the consequences of different ORF15 mutations; ORF15 overexpression may be a less effective therapeutic strategy for patients with certain RPGR mutations.
Materials and Methods
Study approval
The skin biopsies were obtained from previously published XLRP patients and unaffected individuals after written informed consent. All procedures followed the Declaration of Helsinki and approved by the Institutional Review Board of the University of Pennsylvania and UMass Medical School.
Plasmid DNA constructs and antibodies
The plasmids-encoding humanRPGRORF15 or RPGRex1-19 have been described (50,51). Antibody details are provided in Supplementary Table S1.
Fibroblast isolation
The biopsies were dissociated as previously described to obtain primary fibroblasts (39). Briefly, skin biopsies were washed in iodine followed by three washes in phosphate-buffered saline (Invitrogen). The sample was minced with a sharp scalpel in a culture plate and incubated with Dulbecco’s-modified Eagle’s medium (DMEM)/F12 (Invitrogen) supplemented with 20% fetal bovine serum (FBS) (Sigma-Aldrich, Saint-Louis, MO, USA), penicillin and streptomycin at 37°C in a humidified 5% CO2 incubator. Ten days after harvesting, the skin was removed, and fibroblasts were passaged for storage and analysis. After dissociation, cells were grown in DMEM/F12 media containing 10% FBS at 37°C and 5% CO2 incubator. All analyses were done in the fibroblasts at the same passage stages.
Fibroblast transfection assays
Plasmid DNA transfections was performed using FuGENE 6 (Promega) according to the manufacturer’s protocol. All Dicer substrate siRNAs (RPGRex1-19, positive and negative controls, see Supplementary Table S2) were designed and purchased from Integrated DNA Technology and transfected at 20 nM using Lipofectamine RNAiMAX Reagent (Thermo Fisher Scientific), as per manufacturer’s instructions. Cells were plated and transfected 24 h later at 50–70% confluency and induced to form primary cilia 24 h post-transfection by serum starvation for 24 h.
RNA extraction and quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) analysis
Total RNAs were isolated with QIAzol Lysis Reagent (Qiagen) and precipitated according to the manufacturer’s instructions. Reverse transcription was performed with 1 μg of total RNAs with the SuperScript™ First-Strand Synthesis System for RT-PCR (Invitrogen). The resulting cDNA was used to perform gene expression analysis using the BioRad CFX96 qPCR instrument and Power SYBR Green PCR Master Mix (Thermo Fisher Scientific). The isoform-specific quantitative PCR (qPCR) primers are described in Supplementary Table S2.
RNA stability assay
Messenger RNA decay after transcriptional inhibition with actinomycin D treatment was performed as previously described (55). Briefly, the fibroblasts were plated at 50% confluency in six-well plates; 24 h later, cells were treated with 10 μg/ml actinomycin D (Sigma-Aldrich) in DMEM/F12 media for 1, 2, 4, 6 and 8 h. Total mRNA was extracted and subjected to RT-qPCR. Experiments were conducted in triplicates and mRNA decay rate was determined by non-linear regression curve fitting (one-phase decay) using GraphPad Prism Software.
Immunoblotting
The cells were lyzed on ice in Pierce RIPA Buffer (Thermo Scientific) containing Halt™ Protease Inhibitor Cocktail (100X, Thermo Fisher). Lysates were sonicated and centrifugated at 12 000g for 15 min. Supernatants were quantified and 20 μg of total protein lysate was incubated for 5 min at 95°C with Laemmli sample buffer. Proteins were separated on a 4–20% Mini-PROTEAN® TGX™ Precast Protein Gel (BioRad) and transferred to a nitrocellulose membrane (BioRad). The blots were processed using LI-COR western blotting instructions and blocking, and primary and secondary antibodies solutions were prepared in Intercept Blocking buffer (LI-COR). Primary antibodies were incubated overnight at 4°C and secondary antibodies were incubated at room temperature for 2 h (see supplementary table for antibodies details and dilution information). Protein expression was detected using LI-COR Odyssey Fc detection system and quantified with Image Studio Lite quantification software (LI-COR).
Immunofluorescence, microscopy and co-localization analysis
The primary fibroblasts were plated on coverslips (AmScope CS-R18-100) and fixed in pre-chilled 100% methanol at −20°C for 5 min or in 4% paraformaldehyde/1X Dulbecco’s phosphate-buffered saline (DPBS; potassium chloride 0.2 g/L, potassium dihydrogen phosphate 0.2 g/L, sodium chloride 8 g/L and disodium hydrogen phosphate 1.15 g/L) for 10 min, washed three times in DPBS and blocked in DPBS containing 5% normal goat serum and 0.5% Triton X-100 for 1 h. Primary antibody solution was prepared in blocking solution and incubated overnight at 4°C. After washing in DPBS, cells were incubated in blocking solution containing Alexa-488-conjugated and Alexa-546-conjugated (Invitrogen) for 1 h. Finally, cells were washed in DPBS and incubated for 5 min with DAPI for nuclear staining. Coverslips were mounted with Antifade Mounting Media (Vectashield) and all images were visualized over 63× objective using a Leica DM6 Thunder microscope with a 16 bit monochrome camera.Cilia length was measured on single z-stacks using the ImageJ software and line and measure tool. An average of 15 cilia were measured in at least five different coverslip area using acetylated tubulin or ADP-ribosylation factor like 13B (ARL13B) ciliary markers. The IMARIS platform was used to process the two channel intensities simultaneously and to measure the degree of overlap of the two channels. After setting intensity threshold, the basal body marker γ-tubulin was used to define a region of interest and a co-localization channel was built on different z-stacks to obtain the percentage of co-localized signal and intensity. When needed, adjacent optical sections were merged in the z-axis to form a projected image (ImageJ). Co-localization value was plotted using the GraphPad Prism software.
Statistics
All data are presented as means ± standard derivation from the mean. Two groups comparison was analyzed by Student’s t-tests using the Prism GraphPad software. Multiple groups comparison was performed with one-way analysis of variance. Differences between groups were considered statistically significant if P < 0.05. The statistical significance is denoted with asterisks (*P < 0.05; **P < 0.01; ***P < 0.005; ****P < 0.001).
Funding
National Institutes of Health (EY022372 and EY029050 to H.K., EY028273 and P30EY01583 to V.R.M.C. and EY017549 to S.G.J. and A.V.C.); Medical Research Council UK (MR/T002735/1, MR/J004553/1); European Research Council (323147-STEMRD); Retina UK (GR576 to R.R.A.); Research to Prevent Blindness Unrestricted Grant Funds to Scheie Eye Institute.Click here for additional data file.
Authors: Wei Chieh Huang; Alan F Wright; Alejandro J Roman; Artur V Cideciyan; Forbes D Manson; Dina Y Gewaily; Sharon B Schwartz; Sam Sadigh; Maria P Limberis; Peter Bell; James M Wilson; Anand Swaroop; Samuel G Jacobson Journal: Invest Ophthalmol Vis Sci Date: 2012-08-15 Impact factor: 4.799
Authors: A J Mears; S Hiriyanna; R Vervoort; B Yashar; L Gieser; S Fahrner; S P Daiger; J R Heckenlively; P A Sieving; A F Wright; A Swaroop Journal: Am J Hum Genet Date: 2000-09-01 Impact factor: 11.025
Authors: Anneke I den Hollander; Robert K Koenekoop; Suzanne Yzer; Irma Lopez; Maarten L Arends; Krysta E J Voesenek; Marijke N Zonneveld; Tim M Strom; Thomas Meitinger; Han G Brunner; Carel B Hoyng; L Ingeborgh van den Born; Klaus Rohrschneider; Frans P M Cremers Journal: Am J Hum Genet Date: 2006-07-11 Impact factor: 11.025
Authors: Zuzanna Bukowy-Bieryłło; Ewa Ziętkiewicz; Niki Tomas Loges; Mariana Wittmer; Maciej Geremek; Heike Olbrich; Manfred Fliegauf; Katarzyna Voelkel; Ewa Rutkiewicz; Jonathan Rutland; Lucy Morgan; Andrzej Pogorzelski; James Martin; Eric Haan; Wolfgang Berger; Heymut Omran; Michał Witt Journal: Pediatr Pulmonol Date: 2012-08-06
Authors: Jason Charng; Artur V Cideciyan; Samuel G Jacobson; Alexander Sumaroka; Sharon B Schwartz; Malgorzata Swider; Alejandro J Roman; Rebecca Sheplock; Manisha Anand; Marc C Peden; Hemant Khanna; Elise Heon; Alan F Wright; Anand Swaroop Journal: Hum Mol Genet Date: 2016-12-15 Impact factor: 6.150
Authors: J Schlegel; J Hoffmann; D Röll; B Müller; S Günther; W Zhang; A Janise; C Vössing; B Fühler; J Neidhardt; H Khanna; B Lorenz; K Stieger Journal: Transl Res Date: 2018-08-24 Impact factor: 7.012
Authors: B S Pawlyk; O V Bulgakov; X Sun; M Adamian; X Shu; A J Smith; E L Berson; R R Ali; S Khani; A F Wright; M A Sandberg; T Li Journal: Gene Ther Date: 2015-09-08 Impact factor: 5.250
Authors: Laura Moreno-Leon; Marco A Quezada-Ramirez; Evan Bilsbury; Courtney Kiss; Andrea Guerin; Hemant Khanna Journal: Front Genet Date: 2022-08-19 Impact factor: 4.772
Authors: Vlasta Hadalin; Maja Šuštar; Marija Volk; Aleš Maver; Jana Sajovic; Martina Jarc-Vidmar; Borut Peterlin; Marko Hawlina; Ana Fakin Journal: Genes (Basel) Date: 2021-03-29 Impact factor: 4.096