Molecular Vision 2025; 31:526-537
<http://www.molvis.org/molvis/v31/526>
Received 12 October 2025 |
Accepted 18 December 2025 |
Published 20 December 2025
Sarah Pereira de Freitas Cenachi,1 Maria Frasson,1 Anna Laura Marques Nascentes,2 Anna Luiza Braga Albuquerque,3 Rodrigo Rezende Arantes,2 Virgínia Mares,1 Luiz Armando Cunha De Marco,3 Márcio Bittar Nehemy1
1Department of Ophthalmology, Federal University of Minas Gerais, Belo Horizonte, Brazil; 2Genetics Service, Federal University of Minas Gerais, Belo Horizonte, Brazil; 3Department of Surgery, Federal University of Minas Gerais, Belo Horizonte, Brazil
Correspondence to: Sarah Pereira de Freitas Cenachi, Department of Ophthalmology 190 Alfredo Balena Avenue, Belo Horizonte - MG - Brazil - 30.190-090; Phone: +55 31 99588-4359; FAX: +55 31 3409-964; email: sarahcenachi@gmail.com
Purpose: To identify causative genetic variants and associated clinical phenotypes in patients of a tertiary referral center in Brazil with suspected retinitis pigmentosa (RP).
Methods: RP diagnosis was established based on predefined clinical criteria. The patients underwent detailed ophthalmologic assessments and multimodal retinal imaging. Genomic DNA was analyzed using a next-generation sequencing (NGS) panel targeting 238 genes associated with inherited retinal diseases.
Results: Among 55 patients, the genetic diagnostic yield was 71% (39/55), with 13 novel variants identified. The most frequently implicated genes were RHO, RPGR, and USH2A, accounting for approximately 50% of genetically diagnosed cases. Fundus autofluorescence revealed patchy hypoautofluorescence surrounding the vascular arcades as the most frequent finding. On spectral-domain optical coherence tomography, the pattern of ellipsoid zone presentation in the central macula was significantly correlated with best-corrected visual acuity (p<0.001).
Conclusions: This study delineates the genetic and phenotypic spectrum of RP in a tertiary Brazilian referral center, highlighting the utility of NGS for molecular diagnosis and clinical management.
Retinitis pigmentosa (RP) encompasses a group of progressive inherited retinal diseases (IRDs) characterized by the primary degeneration of rod photoreceptors, followed by the loss of cone photoreceptors [1]. RP is the most prevalent hereditary retinal degenerative disorder, with an estimated global prevalence of 1 in 4,000 individuals and approximately 2.5 million affected worldwide [2,3]. A Brazilian study identified RP as the most common IRD, with the nonsyndromic form accounting for roughly 35% of cases [4]. Nonsyndromic is the main form of RP. In addition, 20% to 30% of patients present with a syndromic form of RP associated with extraocular abnormalities [1].
RP typically manifests during the first or second decade of life, with initial symptoms including nyctalopia, followed by progressive concentric visual field constriction. Central visual function often remains relatively preserved until advanced disease stages. Classic fundus findings include bone spicule pigmentation in the peripheral retina, attenuation of retinal vessels, and optic disc pallor [1].
Retinitis pigmentosa is exceptionally heterogeneous. This includes (1) genetic heterogeneity, with many different genes causing the same disease phenotype; (2) allelic heterogeneity, with many different disease-causing mutations in each gene; (3) phenotypic heterogeneity, with different mutations in the same gene causing different diseases; and (4) clinical heterogeneity, with the same mutation in different individuals producing different clinical consequences, even among members of the same family [5,6]. Nowadays, more than 90 genes have been linked to RP, and this number will likely increase as diagnostic testing techniques improve [7,8]. Furthermore, RP is inherited in a complex manner, involving autosomal dominant inheritance (around 20%), autosomal recessive inheritance (around 30%), and X-linked inheritance (around 10%). Occasionally, digenic inheritance and mitochondrial inheritance also have been reported, highlighting the typical genetic heterogeneity of RP [9].
Despite these complexities, individualized characterization remains essential to inform prognosis, guide genetic counseling, and identify potential therapeutic options. Next-generation sequencing (NGS) has emerged as a powerful diagnostic tool, offering high sensitivity and accuracy in detecting pathogenic variants in IRDs [10-12]. In this study, we applied NGS to a cohort of Brazilian patients with clinically suspected RP, aiming to delineate the genetic landscape and associated clinical features.
To explore the clinical features, we used data from multimodal imaging assessment, particularly optical coherence tomography (OCT) of the macula. Because the disease involves progressive loss of cones over its course, morphological assessment of photoreceptors in the macular region can be useful for estimating residual central visual function in patients with RP, and OCT is a well established method for in vivo examination of retinal architecture [13]. Several OCT studies in patients with RP have examined the correlation between retinal structure and visual function [14-17].
This study adhered to the principles of the Declaration of Helsinki and was approved by the local ethics committee of the Ezequiel Dias Foundation (Protocol No. 53,981,621.1.3001.9507). Informed consent was obtained from all participants or their legal guardians.
Patients with clinically suspected RP were recruited from the Retina Department at São Geraldo Hospital, Federal University of Minas Gerais, a tertiary referral center in the state of Minas Gerais that provides high-complexity ophthalmologic care to patients from Minas Gerais and other states of Brazil. From June 2022 to May 2024, a total of 55 participants were enrolled, encompassing 36 unrelated patients clinically diagnosed with suspected RP based on the following inclusion criteria: (1) history of nyctalopia and progressive visual acuity loss, (2) visual field constriction, and (3) characteristic fundoscopic findings, including peripheral bone spicule pigmentation, vascular attenuation, and optic disc pallor [1,8]. Demographic data, age of symptom onset, clinical symptoms, family history, and consanguinity status were documented. Syndromic RP cases were excluded [1,9].
The patients underwent comprehensive ophthalmologic evaluation, including best-corrected visual acuity (BCVA), slit-lamp examination, and multimodal retinal imaging. Imaging modalities included color fundus photography (FP; Optos Daytona, Dunfermline, UK, or Topcon, Tokyo, Japan), fundus autofluorescence (FAF; Optos Daytona), and spectral-domain OCT (SD-OCT; Heidelberg Engineering, Heidelberg, Germany). BCVA was assessed using the early treatment diabetic retinopathy study chart and converted to logarithm of the minimum angle of resolution units for statistical analysis. For macular OCT analysis, dense scans as well as vertical and horizontal line scans passing through the fovea were used.
Statistical analyses were conducted on two parameters regarding macular OCT: integrity of the ellipsoid zone and presence of hyperreflective foci in the outer retina. These features were classified according to their location on SD-OCT: within the central 3-mm diameter of the macula, corresponding to the inner ring of the early treatment diabetic retinopathy study grid (central macula), and in the outer retinal layers beyond the central 3 mm (perifoveal macula) [17]. To assess the ellipsoid zone, participants were classified based on the state of the ellipsoid zone (EZ), which could be preserved (when the line was intact), disorganized (when a large part of the integrity of the line was not preserved) [18], or not detected when the zone could not be visualized (not necessarily destroyed but not discernible) [13]. Hyperreflective foci were defined as well circumscribed hyperreflective lesions (with reflectivity at least as bright as the retinal pigment epithelium [RPE] band) and a maximum size of 50 μm. For a more detailed analysis of the correlation between these features and visual acuity, eyes with poor fixation due to extremely low visual acuity (<0.01), vitreomacular traction, macular hole, or macroaneurysm were excluded.
Genomic DNA was extracted from saliva samples of all affected individuals. A custom NGS panel targeting 238 genes associated with IRDs was used (Illumina platform). Sequence reads were aligned to the GRCh38 human genome reference, and variants were filtered by proprietary software Abracadabra (Mendelics, Sao Paulo, Brazil) to exclude low-quality calls, common polymorphisms (based on public and internal frequency databases), and variants outside the targeted gene panel, retaining only high-confidence rare variants. DNA from relatives was not available for segregation analysis.
The variants were assessed for pathogenicity according to the American College of Medical Genetics and Genomics guidelines [19]. Annotation used multiple databases, including OMIM, gnomAD, dbSNP, Varsome, and ClinVar.
Variants were classified into five categories: pathogenic, likely pathogenic, variant of uncertain significance, likely benign, and benign. Criteria for classification included minor allele frequency (<0.01% for dominant, <0.5% for recessive variants), previous reports in ClinVar or peer-reviewed literature, in silico deleterious predictions, and consistency with known phenotype and inheritance patterns. Variants not reported in HGMD, ClinVar, or the literature were considered novel [20]. Possible origins of inheritance of each variant were evaluated based on populations in which they had been previously reported, based on gnomAD Exomes.
Data were tabulated and analyzed using Microsoft Excel 365 (version 2505, Build 16.0.18827.20102; Microsoft, Redmond, WA). Data were initially analyzed descriptively. For categorical variables, absolute and relative frequencies were presented, while for numerical variables, summary measures (mean, quartiles, minimum, maximum, and standard deviation) were reported.
To compare mean visual acuity according to specific characteristics (EZ integrity and presence of hyperreflective foci), a linear mixed-effects model was employed to account for possible dependence between observations from both eyes of the same patient. When significant mean differences were identified, groups with distinct means were determined using ad hoc multiple comparisons under the estimated model, with Bonferroni correction applied to maintain the overall significance level. The linear mixed model assumes data normality, which was assessed using the Kolmogorov–Smirnov test. However, as noted by Gelman and Hill, deviations from normality do not bias the estimates [21]. All statistical tests were performed at a 5% significance level. Analyses were conducted using SPSS version 20.0 (SPSS, Chicago, IL) and STATA version 18 (StataCorp LLC, College Station, TX).
Of the 55 patients, no pathogenic variants were detected in nine cases, and two cases carried mutations in CHM, associated with choroideremia rather than RP. The remaining 44 patients were evaluated for RP-related mutations, yielding 67 variants: 44 pathogenic, 16 likely pathogenic, 5 variants of uncertain significance, and 2 likely benign. After excluding the two cases carrying likely benign variants, 42 patients were retained for analysis. Of these, 39 patients had pathogenic and likely pathogenic variants, resulting in a diagnostic yield of 71% (39/55; Figure 1).
Variants were identified in 19 RP-related genes: AHI1, CLN3, CLRN1, CRB1, EYS, FSCN2, OPA1, PCARE, PDE6A, PHYC, PRCD, PRPF31, PRPF8, PRPH2, RHO, RPE65, RPGR, SDCCAG8, and USH2A. The most frequently implicated genes were RHO, RPGR, and USH2A, collectively accounting for nearly 50% of genetically resolved cases. Variants were mapped to known structural domains of the respective proteins (Figure 2).
Autosomal recessive retinitis pigmentosa (ARRP) was the predominant inheritance pattern, present in 11 genes, followed by autosomal dominant RP (ADRP) in 4 genes, mixed ARRP/ADRP in 3 genes, and X-linked RP (XLRP) in 1 gene (Figure 1). Compound heterozygosity was inferred in 14 patients with ARRP, based on clinical presentation suspected of RP, NGS findings of two distinct variants in the same gene, and the known AR inheritance pattern.
A total of 65 genetic variants were identified, comprising 44 distinct types. Of these distinct types, 13 were novel mutations (6 missense, 6 frameshift, and 1 insertion), 27 variants were previously reported (14 missense, 7 frameshift, 5 nonsense, 1 intronic), and 4 copy number variants CNVs could not be conclusively categorized as novel. Novel variants were identified in 13 genes: AHI1, CRB1, EYS, FSCN2, PCARE, PDE6A, PHYC, PRCD, RHO, RPE65, RPGR, SDCCAG8, and USH2A. Comprehensive genetic findings are summarized in Appendix 1.
Two sisters were found to carry the c.940+1G>T mutation in CHM, consistent with X-linked choroideremia. Both had been clinically misdiagnosed with RP based on overlapping features, including nyctalopia, peripheral pigmentation, vascular narrowing, and optic disc pallor. Based on genetic results, their diagnoses were revised to choroideremia.
Among the 42 analyzed patients, 21 were male and 21 were female, with ages ranging from 9 to 90 years (mean: 39±18 years). Age of symptom onset ranged from 5 to 52 years (mean: 25±14 years; Appendix 2). Consanguinity was reported in 7% of patients (P2, P25, and P34), and a family history of visual impairment was present in 52% (22/42; Appendix 1).
All 42 patients reported visual symptoms. Nyctalopia was the most frequent complaint, reported by 86% of patients (36/42), followed by photophobia in 19% (8/42). BCVA ranged from 20/25 to light perception, with a median of 20/50. Myopia was the predominant refractive error, observed in 59% of eyes (47/80) among the 40 patients who underwent refraction. Bilateral cataracts were identified in 33% of patients (14/42), and among them, 36% (n=5 patients) had undergone lensectomy in at least one eye (Appendix 2).
FP was performed in 39 patients, yielding 78 eyes for analysis. Optic disc pallor was observed in 92% of eyes (72/78) and vascular attenuation in 90% (70/78). Retinal pigment epithelium alterations were present in all eyes, with peripheral pigmentation in the form of bone spicules observed in 97% (76/78). Two eyes showed no fundus pigmentation but exhibited RPE rarefaction in the peripheral retina. All findings were bilateral, highlighting the symmetric presentation characteristic of inherited retinal diseases (Appendix 2, Figure 3, and Figure 4).
FAF was performed in 38 patients. Due to poor fixation, the examination could not be completed in 4 eyes, resulting in 72 eyes analyzed. In the macular region, 35% of eyes (25/72) exhibited a hyperautofluorescent ring, consistent with Robson-Holder ring [22], and 14% (10/72) showed patchy hypoautofluorescence surrounding the fovea, indicative of bull’s-eye atrophy [23]. In the peripheral retina, 79% of eyes (57/72) demonstrated patchy hypoautofluorescence surrounding the vascular arcades. Representative FAF findings are shown in Figure 3 and Figure 4.
SD-OCT was performed in 39 of the 42 patients, totaling 78 eyes analyzed. Retinal cysts were identified in 49% of eyes (38/78), with small cysts accounting for 84% (32/38). An epiretinal membrane was present in 47% of eyes (37/78), with bilateral involvement in 17 of the 20 affected patients. Macular holes were observed in three eyes from three different patients: one eye had a full-thickness macular hole with associated sensory shallow macular detachment (P29), while two had lamellar holes with epiretinal proliferation and adjacent epiretinal membrane (P11 and P42).
Central macular thickness was measured by SD-OCT in 74 of the 78 eyes analyzed, with 4 eyes excluded due to poor fixation. Central macular thickness values ranged from 96 to 630 µm, with a mean of 251 µm and a median of 247 µm. SD-OCT findings are shown in Appendix 1 and Figure 3.
Regarding the correlation between visual acuity and EZ integrity features and the presence of hyperreflective foci in the macular region, 65 eyes from 34 patients were included in the statistical analysis. The data indicate that alterations in the EZ in the central macula on OCT are associated with visual acuity (p<0.001). The mean BCVA in the group with a nondetected EZ was significantly lower than that of the preserved EZ and disorganized EZ groups, with no significant difference observed between the latter two. The analysis of other variables in relation to BCVA did not reach statistical significance (Table 1 and Figure 5).
RP is the most prevalent form of IRD, with over 90 genes currently implicated in its pathogenesis [24]. The genetic heterogeneity of RP is profound, encompassing multiple inheritance modes, variable penetrance, and phenotypic diversity—even among individuals carrying the same pathogenic variant [5]. As clinical features often overlap with other IRDs, accurate molecular diagnosis remains critical for definitive classification, prognosis, genetic counseling, and potential therapeutic stratification [9]. In the present cohort, a targeted NGS panel yielded a genetic diagnosis in 71% of patients, corroborating its role as a powerful diagnostic tool, as reported in previous studies [9,25,26]. Moreover, we identified 13 novel variants, expanding the known mutational spectrum associated with RP.
Variants in three genes—RHO, RPGR, and USH2A—accounted for nearly half of the cases resolved by genetics. These genes have also been reported among the most prevalent in epidemiological studies on RP conducted both within and outside Brazil [4,9,27-30] (Table 2). In our cohort, 35% of the patients were related, and this information was used to explore the clinical heterogeneity of RP (Figure 3 and Appendix 1).
RHO, encoding rhodopsin, was the most frequently implicated gene. Rhodopsin is a G protein-coupled receptor localized to rod photoreceptor outer segment discs, where it mediates phototransduction upon activation by light [7]. Mutations in RHO are the predominant cause of ADRP, accounting for approximately 25% to 30% of ADRP cases [31]. In the present study, RHO variants were identified in eight patients: c.601G>T in patient 1 (P1), c.316G>A in P2, c.140T>G in P3, and c.403C>T in P4 to P8, who were related. The onset of visual impairment ranged from 5 to 37 years of age. Among the patients carrying RHO variants, P2 exhibited sectoral peripheral retinal involvement, whereas the others demonstrated global peripheral retinal involvement, showing phenotypic heterogeneity. Moreover, clinical heterogeneity was observed; notably, P8, aged 9 years, carries the same mutation as her older sisters but presents a more severe disease spectrum, with more pronounced macular edema and worse visual acuity (Figure 3). Preservation of the fovea was observed in most of these cases, consistent with the “OCT-based transition zone” model proposed by Hood et al., in which the outer nuclear layer remains partially intact over the central photoreceptors but disappears beyond the foveal region [32]. FAF imaging frequently revealed the presence of the Robson-Holder ring. This hyperautofluorescent circular band surrounding the fovea represents the transition zone between relatively preserved central photoreceptors and degenerating peripheral retina [22].
The second most frequently mutated gene was RPGR, associated with approximately 70% of XLRP cases [33]. RPGR is localized to the connecting cilia of photoreceptors and plays a key role in ciliary trafficking, although its precise function remains incompletely understood. XLRP due to RPGR mutations typically manifests early and progresses rapidly, often with profound visual impairment [34]. In our study, RPGR variants were identified in seven patients—five males and two females. Age of onset ranged from 5 to 29 years. OCT findings in affected individuals frequently showed macular atrophy of photoreceptor layers. Female carriers demonstrated variable phenotypes, consistent with previous reports of XLRP carrier expressivity [35]. One female patient harboring the c.2405_2406delAG variant exhibited significant retinal involvement, including bilateral loss in the EZ, retinal cysts, an epiretinal membrane, and a unilateral macular hole (Figure 5).
The third most common gene was USH2A, previously reported as a major contributor to both syndromic and nonsyndromic RP in various populations [36]. Located on chromosome 1q41, USH2A encodes usherin, a basement membrane protein involved in cochlear and retinal development [37]. Four patients carried USH2A variants—three with compound heterozygous mutations and one with a homozygous variant. Symptom onset was later than in RHO or RPGR carriers (range: 38–49 years). No patient reported hearing loss at the time of evaluation; however, audiological assessment was recommended.
The variants were most commonly reported in European and African populations. Yet, several were also reported in Asian, admixed American, and Ashkenazi Jewish populations, reflecting Brazil’s population genetic diversity. Brazil’s genetic architecture represents the world’s largest recently admixed population, rooted in an intricate interplay between native indigenous populations, colonization effects, assertive mating, and modern immigration [38]. This cohort represents the mosaic of global haplotypes represented in Brazilian genomic diversity.
Widefield imaging technology was employed in this study to document retinal changes, enhancing the characterization of RP subtypes. FAF analysis identified hypoautofluorescent patches surrounding the arcades, corresponding to areas of pigmentation and retinal atrophy that may not be as evident on retinography (Figure 4). Identification of this finding may be useful for early detection of suspected RP cases and for differential diagnosis, such as uveitis. FAF can also facilitate RP subtype characterization in conjunction with genetic testing. In some patients, FAF patterns such as the “infinity sign” (peripapillary and macular sparing) and “broad-banded hyperautofluorescent leading edge” were observed—patterns previously associated with EYS mutations (Figure 4) [39].
SD-OCT provides high-resolution cross-sectional imaging of retinal morphology and is a valuable tool for assessing disease severity and progression in RP [40-42]. This disease is characterized by progressive degeneration of the photoreceptor and RPE layers, initially manifesting in the mid-peripheral retina and subsequently extending toward the macula [1,18].
In the present study, we observed a significant correlation between the degree of EZ disruption in the central 3 mm of the macula and worsening visual acuity, reinforcing the importance of central photoreceptor integrity as a structural marker of visual function [13]. In contrast, the presence of hyperreflective foci in this region did not show statistical significance in our study, although previous reports have identified them as potential biomarkers of RP progression [16,17]. This discrepancy may be related to our sample size, which may still be insufficient for this level of analysis.
Future studies incorporating functional assessments such as visual field testing and microperimetry, along with the application of artificial intelligence—similar to advances achieved in geographic atrophy associated with age-related macular degeneration—may yield important insights into retinal dystrophies [43].
Notably, two patients harboring a pathogenic CHM splice-site variant (c.940+1G>T) were initially misdiagnosed with RP. The CHM gene causes X-linked choroideremia, a chorioretinal dystrophy that clinically mimics RP but is pathophysiologically distinct. FAF imaging in these patients revealed scalloped mid-peripheral chorioretinal atrophy, a hallmark of choroideremia (Figure 4) [9,44]. This underscores the importance of multimodal retinal imaging and molecular testing in achieving an accurate diagnosis.
Additionally, mutations were identified in genes associated with syndromic disorders. These included AHI1, associated with Joubert syndrome [45], and CLN3, linked to Batten disease [46]. Although the patients reported no systemic symptoms beyond ocular manifestations, given the potential for an association with broader systemic diseases, further genetic investigation was recommended.
This study has certain limitations. The NGS panel employed does not detect variants in deep intronic or regulatory regions, nor does it provide robust CNV characterization. Co-segregation analysis was not feasible, limiting the ability to confirm variant phase in recessive cases. Finally, as discussed in the OCT analysis, the relatively small sample size restricts broader genotype–phenotype correlation analyses, although it remains representative of a rare disease.
This study provides a comprehensive molecular and clinical characterization of patients with suspected RP from a tertiary referral center in Brazil. The 71% diagnostic yield highlights the effectiveness of NGS panels in IRD evaluation. Thirteen novel variants were identified, contributing to the expanding mutational landscape of RP. The findings emphasize the importance of integrating genetic testing with clinical and imaging data to refine diagnoses, guide patient management, and enable personalized therapeutic approaches.
Disclosure statement: The authors report no conflicts of interest. The authors alone are responsible for the content and writing of this article. Funding: The genetic tests were partially funded by Novartis AG.