Comparison of 5-year progression of retinitis pigmentosa involving the posterior pole among siblings by means of SD-OCT: a retrospective study
BMC Ophthalmology volume 18, Article number: 153 (2018)
The aim of this study is to analyze and compare the progression of photoreceptor atrophy among siblings affected by retinitis pigmentosa by means of spectral SD-OCT.
Fifty three eyes of 27 patients belonging to 12 family clusters were analyzed. To assess the annual progression rate of photoreceptor atrophy, the ellipsoid zone (EZ) line was measured in OCT sections through the fovea. We used multivariate generalized mixed effects to model the rate of progression and its relation to the initial ellipsoid zone line width.
During our 4.84 years (± 1.44) mean follow up time (range 3–7) 53 eyes were examined. The ellipsoid zone line width declined with a yearly average rate of 76.4 μm (4.16% / year) (p-value < 0.0001). Progression rates were poorly correlated within family clusters (p-value = 0.23) and showed statistical difference between affected siblings (p-value = 0.007). There was no correlation between inter-familiar progression rate and mode of inheritance (p-value = 0.98) as well as between age and ellipsoid zone line width among siblings (p-value = 0.91).
RP could be extremely heterogeneous even among siblings: an accurate and sensitive method to follow the progression of the disease is fundamental for future development of clinical trials and therapy strategies.
Retinitis pigmentosa (RP) is a group of inherited retinal disorders leading to vision loss and blindness. It affects approximately 1:4000 individuals with variable modes of inheritance (autosomal dominant, autosomal recessive or X-linked) [1,2,3,4,5].
The peculiarity of RP is its substantial heterogeneity: more than 60 genes are involved, with many possible disease-causing mutations on the same gene, and different clinical outcomes may be linked to the same mutation [5,6,7,8,9,10,11].
Despite this heterogeneity, RP patients have some common clinical features: progressive loss of photoreceptors, typically involving the rod system. The characteristic phenotype includes retinal bone-spicule pigmentation, pallor of the optic disk and attenuation of retinal vessel .
No universally accepted therapies are currently available for RP but many trials are ongoing evaluating different therapeutic approaches. Neuroprotection, stem cells, gene therapy, optogenetics, electrical stimulation and retinal prosthesis represent possible potential future approaches to slow down the progression of the disease or to restore visual function in patients affected by retinal dystrophies [13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28]. In this light, having deeper knowledge of the mechanisms of disease and more sensitive methods to study its progression is becoming increasingly important.
Most of the studies evaluating RP natural course are based on electroretinography (ERG) and visual field (VF) data. Both methods have an important limitation in evaluating short-term progression due to their high test-retest variability [29,30,31,32,33].
The objective of our study was to evaluate disease progression among siblings affected by RP involving the posterior pole and to test whether the progression rate was more similar among subjects within family clusters. For the purpose of our study we measured EZ line changes in a 5-year mean follow up.
We retrospectively reviewed data obtained from outpatients of the Retinal Dystrophies department at the University Eye Clinic of San Paolo Hospital in Milan. Among the whole dataset, siblings affected by RP were selected and, upon informed consent, recruited for the study. Within each family siblings were followed for the same period of time and the same number of visits. Diagnosis of RP, clinical and OCT follow up of at least 3 years, stage of disease involving the posterior pole and clearly measurable EZ band on OCT in all past visits were considered as inclusion criteria.
Exclusion criteria included poor OCT scan quality (media opacities, nystagmus), prior vitreo-retinal surgery and not identifiable EZ band on OCT scans.
Diagnosis of RP was based on clinical signs (characteristic bone spicule pigmentation, optic disc pallor, retinal vessel attenuation, visual field constriction, non detectable scotopic electoretinographic waves) and, if available, confirmed by results of genetic analysis. Genetic tests were made in collaboration with MAGI Human Medical Genetics Institute (Rovereto, Italy).
Demographic data and medical history (including the age of diagnosis, the first symptoms onset, inheritance model) were collected from the database. Extent of disease was considered from the date of first diagnosis.
The study was accomplished in compliance with the Declaration of Helsinki and international guidelines.
Spectral domain optical coherence tomography (SD-OCT)
Retinal imaging was obtained using Spectralis HRA and OCT (Heidelberg Engineering, Heidelberg, Germany). For the study we considered single line scans of 30° and composed of 100 averaged images using the automatic eye tracking software horizontally across the fovea.
All patients meeting the inclusion/exclusion criteria were re-examined. In order to compare OCT scans obtained at the moment of recruitment with previous examinations, since the follow up module was not consistently used, we needed to account for the fact that small displacements of the foveal scan could occur at each visit. Thus, for each follow up visit, we selected the highest quality horizontal OCT scan across the fovea and used each of the selected scans as a reference for a new acquisition at the time of the study, using the progression tool of the Heidelberg Eye Explorer (HEYEX) software.
EZ line measures were taken manually in a masked fashion using the calliper of the HEYEX software by two experienced OCT-readers as shown in Fig. 1 and the values from the two graders were averaged. EZ limit was considered where the hyperreflective band decline to zero. For the purpose of the study horizontal scans were analyzed: limits were nasal and temporal to the fovea.
Also, the presence/absence of cystoid macular edema (CME) and epiretinal membrane (ERM) was recorded. CME was defined as the presence of hypo-reflective spaces visible on at least two consecutive scans. ERM was defined as the presence of a hyper-reflective line adherent to the inner retina often together with underlying waves in the retinal surface layer due to tractional forces.
All scans were performed with dilated pupils (using 1% tropicamide).
Only OCT scans of good quality (higher than 25 dB) were used for the measurements.
Scans were retrieved in anonymized form from our Spectralis database. As explained in the previous section, in the final dataset, each patient had many progression rate measurements derived from the difference in the EZ line width measured from the study visit and from each of the past sections taken as reference, divided by the elapsed time (in years). Under the assumption that the progression rate (but not the EZ line width) could be considered homogeneous for slight displacements, these measurements allowed estimation of the progression rate across multiple follow up visits.
Calculations were performed using generalized linear mixed models. For the progression rate of the EZ line width, we studied the correlation of the progression rates and the EZ line width at each baseline scan. Nested random effects were used to obtain a multilevel model of the error distribution to account for clustered observations within the same subject and the same family. Specifically, the innermost grouping factor was the single eye: this was used to account for repeated measures performed on the same eye, allowing the calculation of a separate mean rate for each eye considered. The second grouping factor was the single subject: this factor accounted for the correlation between the two eyes. Finally, family was the outermost grouping factor. The last two random effects were used to calculate intraclass correlation coefficients (ICCs) within the same family or subject, in order to assess how family clustering could affect the rate of progression. The ICCs were calculated using the following formula (for the Family random effect):
where the numerator is the variance attributed to the grouping factor analyzed, while the denominator is the total variance. Group variances are denoted as σ2 while the residual errors are denoted as ε2. The same formula can be applied to calculate the Subject ICC.
A gamma error distribution with a log link function was used to account for the non-constant variance of the strictly positive value of the progression rate. Such a modelling approach also describes the variation in the rate of change as a proportion of the EZ width in a non-linear fashion.
Similar models were used for the other quantities analyzed, changing the error distribution according to the different variable in study.
Twenty seven patients were recruited for the study (53 eyes): 9 males (33.3%) and 18 females (66.7%).
Twelve patients (44.5%) were affected by autosomal recessive RP (arRP), 9 (33.3%) by autosomal dominant RP (adRP) and 6 (22.2%) by Usher syndrome type II.
Table 1 summarizes demographic data: mean age was 42.40 years (SD ±13.70) and mean follow up time was 4.84 years (range 3–7 years).
Table 2 reports genetic results of siblings included in the study: seven out of 12 families (58%) harboured a genetic variation probably involved in the phenotype. For two of them (family n°5 and 9), the genetic result was not fully informative considering that only one genetic variant in heterozygous state was found in the autosomal recessive gene USH2A that is therefore not sufficient to explain the RP phenotype.
Compound heterozygous state was confirmed by segregation analysis in AR families n°6, 11 and 12.
Mean baseline EZ width was 2345.7 μm, mean EZ width at last follow up visit was 1945.4 μm.
All patients showed an EZ line width decrease, the estimated mean progression rate was 76.4 μm/year (4.16% per year, p-value < 0.0001).
Figure 2 shows the relationship between the initial EZ line width in each B-scan and the progression rate over 1 year. The correlation between repeated observations within the same family, the same subject and the same eye was modeled as nested random effects (see Methods). Although only 4 families out of 12 included siblings with different sexes, we tested the effect of sex on the progression rate. Specifically, we used a log-likelihood ratio test to assess whether the inclusion of a random slope in the mixed effect model allowing a change in rate based on subject’s sex for each family cluster provided a significant increase in the goodness of fit. Since the contribution of the sex was not significant (p = 0.79) we excluded this factor from the model. The curved relationship derives from the Gamma distribution chosen to model the variance and the log function used to link the mean rate to the predictor. The estimated mean progression rate from the model was 76.4 ± 1.16 μm/year (Mean ± SE) with a 38% ± 0.08% rate reduction for every 10 μm EZ line width.
From the same model we calculated the ICCs of the Family and Subject random effects to assess the contribution of each grouping factor in explaining the variability of the progression rate (Table 3). The highest ICC (0.175) was attributable to the Subject grouping, while the ICC was 0.089 for the Family grouping factor and 0.068 for the Eye grouping factor, suggesting that most of the variance can be ascribed to subject differences, with a lesser contribution from the family cluster and the specific eye, although the values of all ICCs calculated are small in magnitude. Figure 3 depicts the variability at the family and subject level (on the log scale of the link function). Notice the high variability of the single subject around the estimated mean log-rate intercept for each family.
No statistically significant correlation was found between age and EZ line width at baseline (p-value = 0.75), even among siblings (p-value = 0.91). Similarly, no relationship was found between the progression rate of EZ line width and the mode of inheritance of the disease (adRP versus arRP, p-value = 0.98). Of note, Usher syndrome was considered in the group of arRP.
Lastly, the percentage of eyes affected by CME did not change during the whole follow up time (33.3%). The percentage of eyes presenting an epiretinal membrane slightly increased (from 49 to 54.7%).
Our study shows that in our cohort of patients the rate of disease progression among siblings affected by RP, evaluated by means of SD-OCT, is not homogenous.
Clinical and genetic heterogeneity of RP is well known. Natural course of RP could be extremely different among patients, with period of stability that can last years followed by a rapid decrease of visual functions within weeks. The landscape of identified genetic modifications responsible for the disease is becoming increasingly varied, with more than 60 genes involved, each possibly carrying many different mutations [1, 5,6,7,8,9,10,11,12].
The peculiarity of our study is the intra-familiar evaluation of the progression of RP.
Different methods can be used to study the progression of RP. Full field electroretinogram (ffERG) showed its usefulness in quantitating the natural course of disease: limits of ffERG are the lower accuracy in detecting progression in late stages of disease and the high test-retest variability. Because of its features, ffERG is considered to be more reliable when evaluating long-term rather than short-term changes.
High test-retest variability also affects subjective functional evaluations with perimetric tests. On one hand, visual field test has a large dynamic range allowing progression detection even in late stages of disease. On the other hand, both static and kinetic perimetry depend on patient’s collaboration. Moreover, while the former is computer-based, the latter, which is the gold standard in routine follow up of RP patients, might be hampered by operator’s ability and experience in performing the test.
For the purpose of our study we evaluated the progression by means of EZ changes at SD-OCT. Recently, EZ line width has been the subject of several studies evaluating the rate of RP progression [34,35,36,37,38]. The advantages of using EZ band as marker of progression are the low test-retest variability, the accuracy in detecting small changes and the straightforward scan acquisition and analysis. The most important limit is that it cannot be used in early and late stages of the disease. In late RP, when the outer retina is completely atrophic, EZ band is not easily detectable. Conversely, in initial RP cases, when outer retinal atrophy does not involve the posterior pole, the margins of the EZ band cannot be identified within classical OCT scans. However, the boundaries of this latter limitation are being progressively blurred by the introduction of wide field imaging .
In general, the mean annual rate of progression is our dataset is comparable with previous studies: the EZ line width reduction was 4.16% per year in accordance with 4.9–10.9% [34, 36, 38] reported in published literature. Similarly to published data, we observed that the progression rate of the EZ band atrophy decreases when the margins of the atrophic retina approach the foveal region [34, 36].
When evaluating the variability of progression rate in our cohort of patients, the highest ICC was attributable to subject grouping: in other terms, being part of the same family cluster, and therefore having the same genetic mutations, does not imply a similar rate of progression (Family ICC = 0.089). Interestingly, we could not find a statistically significant correlation between age and EZ line width among siblings. In other terms, younger siblings in one family cluster could show a more advanced stage of the disease compared to older ones.
However, the lack of any statistically significant interaction between progression rate and mode of inheritance in our study could be explained by the small sample size.
Literature about ERM incidence in RP patients is controversial with range varying from 1 to 64% [41, 45, 46]. In our cohort the percentage of eyes with evidence of ERM ranged between 49% at baseline and 54.7% at the last follow up exam considered.
In conclusion, as our data suggested, RP could be extremely heterogeneous even among siblings: an accurate and sensitive method to follow the progression of the disease is fundamental for future development of clinical trials evaluating new therapeutic strategies for RP and others retinal dystrophies. For this purpose, EZ line width evaluation could be considered one of the markers of disease progression. However, the sources of variability identified in our model had overall low ICCs. A more detailed prospective analysis based on follow up measurements on the same location with retinal tracking, modelling the length of the EZ over time instead of the change in rate, might help better characterize the actual dynamics of the progression and reduce the overall variability of the estimates. Furthermore our data suggest that, due to intra-familiar variability, siblings could not offer evident advantages when used as matched controls in treatment trials.
The limitations of our study are represented by the relatively small sample size, the retrospective nature of the study and, as explained above, the fact that we did not use the same section as reference for all follow up examinations.
Prospective studies evaluating EZ band changes with the same reference and combining other variables as control could be helpful in understanding the precise role of family clustering in determining progression rate of RP.
Autosomal Dominant Retinitis Pigmentosa
autosomal Recessive Retinitis Pigmentosa
Cystoid macular edema
Spectral Domain Optical Coherence Tomography
Hartong DT, Berson EL, Dryja TP. Retinitis pigmentosa. Lancet. 2006;368(9549):1795–809.
Bundey S, Crews SJ. A study of retinitis pigmentosa in the city of Birmingham. I prevalence. J Med Genet. 1984;21(6):417–20.
Haim M, Holm NV, Rosenberg T. Prevalence of retinitis pigmentosa and allied disorders in Denmark. I Main results. Acta Ophthalmol (Copenh). 1992;70(2):178–86.
Bunker C, Berson E, Bromley W, et al. Prevalence of retinitis Pigmentosa in Maine. Am J Ophthalmol. 1984;97(3):357–65.
Farber DB, Heckenlively JR, Sparkes RS, Bateman JB. Molecular genetics of retinitis pigmentosa. West J Med. 1991;155(4):388–99.
Sahel J, Bonnel S, Mrejen S, Paques M. Retinitis pigmentosa and other dystrophies. Dev Ophthalmol. 2010;47:160–7.
Berger W, Kloeckener-Gruissem B, Neidhardt J. The molecular basis of human retinal and vitreoretinal diseases. Prog Retin Eye Res. 2010;29(5):335–75.
Daiger SP, Sullivan LS, Bowne SJ. Genes and mutations causing retinitis pigmentosa. Clin Genet. 2013;84(2):132–41.
Pierrottet CO, Zuntini M, Digiuni M, et al. Syndromic and non-syndromic forms of retinitis pigmentosa: a comprehensive Italian clinical and molecular study reveals new mutations. Genet Mol Res. 2014;13(4):8815–33.
Keats BJ, Savas S. Genetic heterogeneity in usher syndrome. Am J Med Genet A. 2004;130A(1):13–6.
Bessant DA, Ali RR, Bhattacharya SS. Molecular genetics and prospects for therapy of the inherited retinal dystrophies. Curr Opin Genet Dev. 2001;11(3):307–16. Review
JR H, SL Y, Friedman LH, et al. Clinical findings and common symptoms in retinitis pigmentosa. Am J Ophthalmol. 1988;105:504–11.
Strettoi E, Gargini C, Novelli E, et al. Inhibition of ceramide biosynthesis preserves photoreceptor structure and function in a mouse model of retinitis pigmentosa. Proc Natl Acad Sci. 2010;107(43):18706–11.
Falsini B, Iarossi G, Chiaretti A, et al. NGF eye-drops topical administration in patients with retinitis pigmentosa, a pilot study. J Transl Med. 2016;14(1) https://doi.org/10.1186/s12967-015-0750-3.
Léveillard T, Mohand-Saïd S, Lorentz O, et al. Identification and characterization of rod-derived cone viability factor. Nat Genet. 2004;36:755–9.
Léveillard T, Fridlich R, Clérin E, et al. Therapeutic strategy for handling inherited retinal degenerations in a gene-independent manner using rod-derived cone viability factors. C R Biol. 2014;337(3):207–13.
He Y, Zhang Y, Liu X, et al. Recent advances of stem cell therapy for retinitis pigmentosa. Int J Mol Sci. 2014;15(8):14456–74.
Uy HS, Chan PS, Cruz FM. Stem cell therapy: a novel approach for vision restoration in retinitis pigmentosa. Med Hypothesis Discov Innov Ophthalmol. 2013;2(2):52–5.
Acland GM, Aguirre GD, Ray J, et al. Gene therapy restores vision in a canine model of childhood blindness. Nat Genet. 2001;28(1):92–5.
Bainbridge JW, Smith AJ, Barker SS, et al. Effect of gene therapy on visual function in Leber's congenital amaurosis. N Engl J Med. 2008;358(21):2231–9.
Wright AF. Long-term effects of retinal gene therapy in childhood blindness. N Engl J Med. 2015;372(20):1954–5.
Busskamp V, Roska B. Optogenetic approaches to restoring visual function in retinitis pigmentosa. Curr Opin Neurobiol. 2011;21(6):942–6.
Jacobson SG, Sumaroka A, Luo X, Cideciyan AV. Retinal optogenetic therapies: clinical criteria for candidacy. Clin Genet. 2013;84(2):175–82.
Schatz A, Röck T, Naycheva L, et al. Transcorneal electrical stimulation for patients with retinitis Pigmentosa: a prospective, randomized, sham-controlled exploratory study. Invest Opthalmol Vis Sci. 2011;52(7):4485.
Humayun MS, Dorn JD, Dacruz L, et al. Interim results from the international trial of second Sight’s visual prosthesis. Ophthalmol. 2012;119:779–88.
Dacruz L, Bf C, Dorn J, et al. The Argus II epiretinal prosthesis system allows letter and word reading and long-term function in patients with profound vision loss. Br J Ophthalmol. 2013;97:632–6.
Zrenner E, Bartz-Schmidt K, Benav H, et al. Subretinal electronic chips allow blind patients to read letters and combine them to words. Proc Biol Sci. 2011;278:1489–97.
Kitiratschky V, Stingl K, Wilhelm B, et al. Safety evaluation of “retina implant alpha IMS”—a prospective clinical trial. Graefes Arch Clin Exp Ophthalmol. 2014;253(3):381–7.
Gränse L, Ponjavic V, Andréasson S. Full-field ERG, multifocal ERG and multifocal VEP in patients with retinitis pigmentosa and residual central visual fields. Acta Ophthalmol Scand. 2004;82(6):701–6.
Seiple W, Clemens CJ, Greenstein VC, et al. Test-retest reliability of the multifocal electroretinogra and Humphrey visual fields in patients with retinitis pigmentosa. Doc Ophthalmol. 2004;109(3):255–72.
Holopigian K, Greenstein V, Seiple W, Carr RE. Rates of change differ among measures of visual function in patients with retinitis pigmentosa. Ophthalmol. 1996;103(3):398–495.
Grover S, Fishman GA, Birch DG, et al. Variability of full-field electroretinogram responses in subjects without diffuse photoreceptor cell disease. Ophthalmol. 2003;110(6):1159–63.
Fishman GA, Chappelw AV, Anderson RJ, et al. Short term inter-visit variability of ERG amplitudes in normal subjects and patients with retinitis pigmentosa. Retina. 2005;25(8):1014–21.
Birch DG, Locke KG, Wen Y, et al. Spectral-domain optical coherence tomography measures of outer segment layer progression in patients with X-linked retinitis pigmentosa. JAMA Ophthalmol. 2013;131(9):1143–50.
Cai CX, Locke KG, Ramachandran R, et al. A comparison of progressive loss of the ellipsoid zone (EZ) band in autosomal dominant and x-linked retinitis pigmentosa. Invest Ophthalmol Vis Sci. 2014;55(11):7417–22.
Sujirakul T, Lin MK, Duong J, et al. Multimodal Imaging of Central Retinal Disease Progression in a 2-Year Mean Follow-up of Retinitis Pigmentosa. Am J Ophthalmol. 2015;160(4):786–98.
Hariri AH, Zhang HY, Ho A, et al. Trial of oral Valproic acid for retinitis Pigmentosa group. Quantification of Ellipsoid Zone Changes in Retinitis Pigmentosa Using en Face Spectral Domain-Optical Coherence Tomography. JAMA Ophthalmol. 2016;134(6):628–35.
Lima LH, Burke T, Greenstein VC, et al. Progressive constriction of the hyperautofluorescent ring in retinitis pigmentosa. Am J Ophthalmol. 2012;153(4):718–27.
McNabb RP, Grewal DS, Mehta R, et al. Wide field of view swept-source optical coherence tomography for peripheral retinal disease. Br J Ophthalmol. 2016;100(10):1377–82.
Adackapara CA, Sunness JS, Dibernardo CW, et al. Prevalence of cystoid macular edema and stability in oct retinal thickness in eyes with retinitis pigmentosa during a 48-week lutein trial. Retina. 2008;28(1):103–10.
Testa F, Rossi S, Colucci R, et al. Macular abnormalities in Italian patients with retinitis pigmentosa. Br J Ophthalmol. 2014;98(7):946–50.
Hirakawa H, Iijima H, Gohdo T, Tsukahara S. Optical coherence tomography of cystoid macular edema associated with retinitis pigmentosa. Am J Ophthalmol. 1999;128(2):185–91.
Chung H, Hwang JU, Kim JG, Yoon YH. Optical coherence tomography in the diagnosis and monitoring of cystoid macular edema in patients with retinitis pigmentosa. Retina. 2006;26(8):922–7.
Makiyama Y, Oishi A, Otani A, et al. Prevalence and spatial distribution of cystoid spaces in retinitis pigmentosa: investigation with spectral domain optical coherence tomography. Retina. 2014;34(5):981–8.
Grigoropoulos VG, Emfietzoglou J, Nikolaidis P, et al. Optical coherence tomography findings in patients with retinitis pigmentosa and low visual acuity. Ophthalmic Surg Lasers Imaging. 2010;41(1):35–9.
Hagiwara A, Yamamoto S, Ogata K, et al. Macular abnormalities in patients with retinitis pigmentosa: prevalence on OCT examination and outcomes of vitreoretinal surgery. Acta Ophthalmol. 2011;89(2):e122–5.
Availability of data and materials
The dataset used and/or analysed in the study available from the corresponding author upon reasonable request.
Ethics approval and consent to participate
The study was approved by the internal review boards of San Paolo Hospital – University of Milan. Written informed consent was obtained from each participant.
Consent for publication
Prof Luca Rossetti is a section editor for Glaucoma of BMC Ophthalmology.
Other authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Colombo, L., Montesano, G., Sala, B. et al. Comparison of 5-year progression of retinitis pigmentosa involving the posterior pole among siblings by means of SD-OCT: a retrospective study. BMC Ophthalmol 18, 153 (2018). https://doi.org/10.1186/s12886-018-0817-z