Abstract

The role of attentional mechanisms in peripheral vision loss remains an outstanding question. Our study was aimed at determining the effect of genetically determined peripheral retinal dystrophy caused by Retinitis Pigmentosa (RP) on visual cortical function and tested the recruitment of attentional mechanisms using functional magnetic resonance imaging (fMRI). We included thirteen patients and twenty-two age- and gender-matched controls. We analyzed cortical responses under attentional demands and passive viewing conditions while presenting a visual stimulus covering the central and paracentral visual field. Brain activity was studied in visual areas V1, V2, and V3 as well as in cortical regions of interest corresponding to the preserved and the damaged visual field. The influence of visual field extent and age of disease onset were also investigated. Cortical thickness of visual areas was also measured. We found that cortical visual responses under attentional demands were increased in patients with larger degeneration of visual field, as demonstrated by significant interaction effects between group and task conditions. Moreover, activation during the task condition was increased for patients in two cortical regions of interest corresponding to the preserved and damaged visual field, specifically in patients with severe visual field loss. These findings were observed in the presence of preserved visual cortical structure. We conclude that RP patients have enhanced visual attention recruitment despite their retinal degeneration, while cortical structure and overall response levels remain intact. The unmasking of feedback signals from higher level visual regions involved in attentional processes may explain the increased cortical responses. These findings are relevant for the design of strategies for treating retinal diseases, based on attentional cuing.

1. Introduction

Human studies concerning the effects of peripheral retinal loss on adult visual cortical structure and function are scarce. Previous studies have mainly addressed central retinal disorders such as macular degeneration or other hereditary retinal dystrophies and diseases such as glaucoma [13]. Peripheral and central visual information is differentially routed in the brain [4, 5]. Thus, neural adaptation mechanisms might differ when central or peripheral visual degeneration occurs. Previously, we found evidence for visual retinotopic reorganization in RP (peripheral regions responding to more central representations) [6]. It is nevertheless also important to assess the impact of peripheral retinal loss on visual attentional mechanisms. This may provide useful information in the context of low and high level strategies for treating different retinal diseases. Additionally, there is still an ongoing debate on the nature of adult brain functional reorganization induced by retinal diseases [710].

Our study is aimed at determining the interaction between attentional mechanisms and peripheral retinal dystrophy caused by Retinitis Pigmentosa (RP) on brain function using magnetic resonance imaging (MRI). RP is an inherited degeneration of photoreceptors that initially affects the peripheral retina and later advancing towards the central retina. The onset age varies from infancy to adulthood. The disease manifestations comprise night blindness, tunnel vision, and possibly blindness in severe stages [11, 12]. In a functional magnetic resonance imaging (fMRI) case series study of three RP patients [13], authors reported task-dependent changes in cortical responses in the lesion projection zone (LPZ—the cortical region that no longer receives input due to a bilateral retinal lesion or scotomata [1]). They suggested the unmasking of feedback signals from higher-order visual areas under attentional demands when retinal input signals are lost. Another fMRI report with one RP patient found no evidence of functional alterations in the LPZ [14]. Our recent study in a relatively large cohort showed clear topological evidence for reorganization, which was dependent on the long-term extent of visual loss [6].

Contrasting with RP, central vision is primarily affected in macular degeneration. Some authors claimed that the deafferented cortical neurons in the primary visual cortex become responsive to inputs from the peripheral retina in this pathology [1520]. However, other studies have questioned such visual cortical alterations by reporting the existence of a silent LPZ [14, 2125]. Some of the previous studies showed that visual cortical alterations in macular degeneration are associated with the severity of retinal function loss, arguing that large-scale reorganization only occurs when there is a complete foveal visual loss [15, 18, 20]. However, other researchers did not find any signs for cortical reorganization in a large cohort of patients without foveal sparing [22]. Moreover, the influence of age of disease onset on the degree of cortical alterations is not clear [15, 22], although there is some evidence for larger reorganization in macular degeneration patients with earlier forms of the disease [19]. The reduced numbers of participants [21, 25] and the difference in stimuli and tasks used [26] may have also contributed to the controversy in the reported macular degeneration studies ([1], [7].

Our study was aimed at investigating the effect of peripheral retinal dystrophy caused by RP on brain attentional mechanisms using fMRI taking into account the effect of visual field extent and age of disease onset. Our hypothesis stated that, in addition to the previously demonstrated reorganization, visual cortical responses were also altered as a function of attentional demands in RP patients due to the lack of peripheral retinal bottom-up input. Moreover, we hypothesized that these alterations were more prominent for RP patients with more constricted visual fields and earlier disease onset.

2. Methods

2.1. Participants

The participants selected for this study were also included in a previous work from our group on visual retinotopy [6]. We included 13 RP individuals (8 males and 5 females; mean age years; age range 20 - 66 years; 12 right-handed and 1 left-handed; self-reported symptomatic age of onset range 2 - 39 years, resulting in symptomatic duration range 6 - 42 years) and 22 control subjects (11 males and 11 females; mean age years; age range 23 - 66 years; 21 right-handed and 1 left-handed). Both groups were matched for age (, nonsignificant (NS)), gender (, NS), and handedness (, NS) ratio. Patients were recruited at the Centro Hospitalar e Universitário de Coimbra. The control group participants were local volunteers. The study was conducted in accordance with the Declaration of Helsinki and was approved by the Ethics Committee of Faculty of Medicine of the University of Coimbra. Written informed consent was obtained from all participants. Exclusion criteria were intracranial abnormalities, movement during MRI acquisitions, fixation instability, visual alterations in control subjects, or visual alterations other than RP for patients (e.g., diabetic retinopathy or glaucoma).

2.2. Ophthalmological Assessment

For each participant, we measured visual acuity with a decimal chart (converted to logarithm of Minimum Angle of Resolution (logMAR) scale), average cortical thickness, and retinal nerve fiber layer (RNFL) thickness with Frequency Domain Cirrus Ocular Coherence Tomography (OCT, software version 5.1.1.6, Carl Zeiss Meditec AG, USA), and static visual fields with a MonCv3 multifunction perimeter (Metrovision, France) (Figure 1). A detailed description of the methodology used was described in our previous study [6] (Table 1).

2.3. Brain Imaging Procedures

Scanning was performed on a 3 T scanner (Magneton TrioTim, Siemens AG, Germany) at the Portuguese Brain Imaging Network, using a 12-channel birdcage head coil. Two anatomical T1-weighted Magnetization-Prepared Rapid Acquisition with Gradient Echo (MPRAGE) sequences with voxel size, Repetition Time (TR) 2.53 s, Echo Time (TE) 3.42 ms, Flip Angle (FA) 7°, Field Of View (FOV) , and 176 slices were acquired from each participant. Functional sequences consisted of a single shot Echo-Planar Imaging (EPI) acquired in the axial plane parallel to the Anterior Commissure (AC)-Posterior Commissure (PC) plane with voxel size, TR 2 s, TE 39 ms, interslice time (TI) 76 ms, FA 90°, FOV , 26 slices, and imaging matrix.

Stimuli were presented using MRI compatible goggles with refractive correction (VisualSystem, NordicNeurolab, Norway). One eye was covered with a cotton patch while the other received the visual input (the dominant eye, except if it was the eye with the worst visual acuity). The RP and control group were matched for the selected eye ratio (, NS) and for the dominance of the selected eye ratio (, NS; 1 patient with missing data). The maximum field of view was (resolution of ).

2.4. MRI Stimuli

Stimuli were designed using Matlab 2011b (the MathWorks Inc., USA) with Psychophysics Toolbox 3 extensions (http://psychtoolbox.org/). A central red-colored cross with 0.78 deg of diameter was used for fixation [6].

2.4.1. Retinotopic Mapping Stimuli

Polar angle and eccentricity stimuli were employed to delineate the cortical visual areas (V1, V2, and V3) using the traveling-wave approach from the standard phase-encoded retinotopic mapping [27]. Polar angle maps were obtained using a black and white flickering checkerboard wedge with 45 deg rotating in an anticlockwise direction (initial angle of 22.50 deg with horizontal axis). Eccentricity maps were obtained using a black and white flickering checkerboard expanding ring (for additional details see [6]). Checkerboard size varied with cortical magnification factor from the center to the periphery [28, 29]. Stimuli flickering frequency was 8 Hz and contrast was ~100%. Each run comprised 2 baseline blocks (~0% contrast; 12 s) at the beginning and end of the run with 4 cycles of polar angle or eccentricity stimuli (48 s each; total duration of a run 216 s). Two runs of polar angle and 2 runs of eccentricity were acquired for each subject.

2.4.2. Attentional Task Stimuli

The set of stimuli of the main experimental task of this study consisted of a random sequence of 2 different sized checkerboard rings pseudorandomly presented during either passive viewing (fixation only) or a one-back visual memory task condition. During the task condition, participants were instructed to press a button every time a ring was the same size as the immediately preceding one. Ring1 was presented at a foveal location (diameter between 0.78 and 1.90 deg), and Ring2 at a parafoveal location of the visual field (diameter between 6.74 and 9.52 deg). Ring thickness varied with the cortical magnification factor from the center to the periphery [28, 29]. Rings appeared during 0.50 s in random intervals of 1.50, 3.50, or 5.50 s within each block. Four passive viewing blocks and 4 task blocks (~100% contrast; flickering frequency of 8 Hz; 10 rings; 36 s) were alternately presented, intercalated with 9 baseline blocks (~0% contrast; 12 s; each run began and ended with one baseline block). Auditory instructions were provided to the subject before each block, depending on the condition: “Rest” for passive viewing or “Answer” for the one-back task. The average response time and percentage of errors during the task were recorded. Participants who answered during passive viewing blocks were excluded. Figure 2 shows a representation of the one-back task stimulus. We aimed at performance matching between patients and controls. We preferred this choice as compared to a possible 2-back task, because we believe that this would lead to quite large executive load.

2.5. MRI Data Processing

Brain imaging analysis was performed with BrainVoyager QX 2.6.1 (Brain Innovation B. V., Netherlands). The two anatomical images of each participant were averaged, reoriented to AC–PC plane, and transformed to Talairach (TAL) space. The image was segmented into cerebral spinal fluid, gray matter, and white matter to create inflated mesh representations of each hemisphere. A cut was manually drawn along the calcarine sulcus, and meshes were flattened for retinotopic maps projection [6, 30].

The preprocessing of functional sequences consisted of scan time correction, temporal high-pass filtering (2 cycles per run), spatial smoothing (FWHM 2 mm), and a correction for small interscan head movements. Participants were excluded if within-run movements exceeded 4 mm (-2 to 2 mm). Polar angle and eccentricity maps were obtained from the average of the two runs, created based on linear regression analysis, and projected onto the flattened surfaces of each subject (statistical maps with ; Figure 1) [6]. Field Sign Maps were automatically created using polar angle and eccentricity Look-Up Table maps. Retinotopic areas V1 dorsal (V1d), V1 ventral (V1v), V2 dorsal (V2d), V2 ventral (V2v), V3 dorsal (V3d), and V3 ventral (V3v) were manually defined for each subject in each hemisphere on flattened meshes.

Statistical analyses were performed on individual data in TAL space using a general linear model (GLM) (z-transformation, False Discovery Rate (FDR) , correction for temporal serial correlations AR(2)) within the retinotopically defined visual areas. Response predictors for the visual memory task were obtained, and beta values evoked by each stimulus conditions were retrieved: Ring1 Passive Viewing, Ring2 Passive Viewing, Ring1 Task, and Ring2 Task. Figure 3 represents the visual cortical responses for all predictors during the visual memory task for RP patients and control participants.

A multistudy GLM (random fixed effects, -transformation, FDR , AR(2)) was also run in two different regions of interest along the calcarine sulcus (V1). The functional projection zone (FPZ) represented the preserved visual field region, and the LPZ represented either the visual field scotomata in patients or the unstimulated visual field in controls. These cortical regions were manually defined considering the retinotopic eccentricity maps of each participant. Figure 4 illustrates the location of the two regions of interest, the FPZ and the LPZ, in one control participant. Response predictors for the visual working memory task were obtained, and beta values evoked by each stimulus conditions were retrieved inside these regions of interest: Ring1 Passive Viewing, Ring2 Passive Viewing, Ring1 Task, and Ring2 Task.

Cortical thickness was calculated on the retinotopic areas using the standard procedure of BrainVoyager (see [30] for a complete description). To allow an accurate segmentation of white matter–gray matter and gray matter–cerebral spinal fluid boundary, TAL anatomical data were converted to high-resolution . The subcortical structures and the ventricles were filled as white matter. After computation, cortical thickness maps were superimposed on cortical meshes, and mean cortical thickness values of all visual areas were extracted using Matlab BVQXtools toolbox extensions (http://support.brainvoyager.com/available-tools/52-matlab-tools-bvxqtools.html).

2.6. Subgroup Analysis

To understand the influence of the level of peripheral degeneration and disease onset age on visual cortex response, the RP group was divided accordingly to these factors. The RP group was divided accordingly to the extent of visual field measured by the static perimetry test in two subgroups: RPsvf (small visual field) (, patients RP 1 to RP 6) with bilateral visual field diameter under 9.52 deg and RPlvf (large visual field) (, patients RP 7 to RP 13) with bilateral visual field diameter over 9.52 deg (see Table 1). In this way, the RPlvf patients were expected to see the complete visual working memory task stimuli (Ring1 and Ring2, because the maximum diameter of Ring2 was 9.52 deg), whereas most RPsvf patients would only partially see Ring2. Both subgroups and the control group were matched for age (, NS), gender (, NS), handedness (, NS), selected eye (, NS), and selected eye dominance ratio (, NS; 1 patient with missing data). Additionally, the disease onset age (, NS) and the disease symptomatic duration (, NS) were not different between the subgroups RPsvf and RPlvf.

The RP group was also divided accordingly to the disease onset age into two subgroups: RPeo (early onset) (, patients RP 5, 7, 8, 9, 11, and 13) with an onset age lower than 14 years and RPlo (late onset) (, patients RP 1, 2, 3, 4, 6, 10, and 12) with an onset age greater than or equal to 14 years (see Table 1). Patients from subgroup RPeo with earlier disease onset age were expected to have more prominent alterations in visual cortex responses than RPlo patients since some visual plasticity is thought to remain until 14 to 16 years of age [8, 31]. Both subgroups and the control group were matched for age (, NS), gender (, NS), handedness (, NS), selected eye (, NS), and selected eye dominance ratio (, NS; 1 patient with missing data). As expected, the disease onset age (, ) was different between the two subgroups of patients RPeo and RPlo, but the disease symptom duration (, NS) was not different. Table 2 summarizes the information of visual field extent and onset age for all subgroups.

2.7. Statistical Analysis

Statistical analysis was performed with IBM SPSS Statistics, Version 22 (IBM Corporation, USA). Normality and homogeneity of variance were tested using Shapiro-Wilk’s test and Levene’s test, respectively. For data in accordance with these assumptions, statistical parametric tests were performed. Otherwise, nonparametric methods were applied. Bonferroni correction was applied for multiple comparisons ( values presented as corrected ). The epsilon value was used for correction of nonspherical data (Huynh-Feldt for epsilon higher than 0.75 and Greenhouse-Geisser for epsilon lower than 0.75). The significance level was 0.05, and the statistical power was higher than 0.80 for all presented results.

3. Results

3.1. Visual Assessment

The first two parts of the results section (Visual Assessment and Behavioral Data) are focused on ophthalmological features—visual acuity, visual field extent, average retinal thickness, RNFL thickness-, and performance in visual memory task in patient and clinical groups and subgroups. The subsequent parts are focused on the main hypotheses of the article.

3.1.1. All Groups (No Stratification according to Visual Field Extent or Age of Onset)

Visual acuity (left eye (LE) , and right eye (RE) , ) and average retinal thickness (LE , and RE , ) were reduced in both eyes for patients as compared to control participants. Visual field deficit volume (LE , and RE , ) was higher in both eyes for patients as compared to the control group. No differences were found for RNFL thickness in both eyes between the two groups (LE , NS and RE , NS). Moreover, no statistically significant differences were found between the left and right eyes within groups for visual acuity (RP , NS and control , NS), average retinal thickness (RP , NS and control , NS), RNFL thickness (RP , NS and control , NS), visual field deficit volume (RP , NS and control , NS), and visual field extent (RP , NS; controls have a constant visual field extent equal to 48 deg corresponding to the maximum diameter covered by our static perimetry) [6].

3.1.2. Visual Field (Large vs. Small Visual Field) Subgroup Analysis

Visual acuity (LE , and RE , ) and average retinal thickness (LE , and RE , ) were reduced in both eyes for the two RP subgroups as compared to control participants. No differences were found for RNFL thickness in both eyes between the three groups (LE , NS and RE , NS). Visual field deficit volume (LE , and RE , ) was higher in both eyes for both patients’ subgroups as compared to the control group. According to the subgroup division, visual field extent was different between the two subgroups (LE , and RE , ), with RPsvf (small visual field) patients presenting smaller visual fields than RPlvf (large visual field/more preserved). Moreover, no statistically significant differences were found between the left and right eyes within groups for visual acuity (RPsvf , NS, RPlvf , NS, and control , NS), average retinal thickness (RPsvf , NS, RPlvf , NS, and control , NS), RNFL thickness (RPsvf , NS, RPlvf , NS, and control , NS), visual field deficit volume (RPsvf , NS, RPlvf , NS, and control , NS), and visual field extent (RPsvf , NS and RPlvf , NS; controls have a constant visual field extent equal to 48 deg corresponding to the maximum diameter covered by static perimetry).

3.1.3. Age of Onset Subgroup Analysis

Visual acuity (LE , and RE , ) and average retinal thickness (LE , and RE , ) were reduced in both eyes for the two RP subgroups (RPeo, early onset; RPlo, late onset) as compared to control participants. No differences were found for RNFL thickness in both eyes between the three groups (LE , NS and RE , NS). Visual field deficit volume (LE , and RE , ) was higher in both eyes for both patients’ subgroups as compared to the control group. The visual field extent was not different between the two subgroups of patients (LE , NS and RE , NS). Moreover, no statistically significant differences were found between the left and right eyes within groups for visual acuity (RPeo , NS, RPlo , NS, and control , NS), average retinal thickness (RPeo , NS, RPlo , NS, and control , NS), RNFL thickness (RPeo , NS, RPlo , NS, and control , NS), visual field deficit volume (RPeo , NS, RPlo , NS, and control , NS), and visual field extent (RPeo , NS and RPlo , NS; controls have a constant visual field extent equal to 48 deg corresponding to the maximum diameter covered by static perimetry).

In sum, ophthalmological tests showed decreased patients’ visual acuity, visual field extent, and average retinal thickness, while RNFL thickness was preserved. Importantly, subgroups defined by visual field extent and age of onset only differ in the visual field extent or in the age of onset of the disease, respectively. Table 2 shows a summary of participants’ ophthalmologic characterization, presenting the visual parameter values for each group and subgroups.

3.2. Visual Memory Task: Behavioral Data

The mean response time (, NS) and response error (, NS) were not different between the RP and control groups during the performance of the one-back visual memory task, and both groups were actually near ceiling levels. Concerning the visual field subgroups, we did not find differences among the RPsvf and RPlvf subgroups and the control group in response time (, NS) and error (, NS). For the onset age subgroups, the analysis did not show differences among the RPeo and RPlo subgroups and the control group in response time (, NS) and error (, NS). Groups were therefore behaviorally matched.

3.3. Visual Memory Task: Responses in Functional (FPZ) and Lesion Projection Zone (LPZ)

Here, we tested the hypothesis that compensatory allocation of visual attention mechanisms do occur in patients.

3.3.1. All Group

The beta values of the task predictors (Ring1 Passive Viewing, Ring2 Passive Viewing, Ring1 Task, and Ring2 Task) were analyzed between groups for each condition (Condition Passive Viewing and Task), each ring (Ring Ring1 and Ring2), each cortical zone (Zone FPZ and LPZ in V1, see Figure 4), and each hemisphere (Hemisphere Left and Right) with repeated measures ANOVA. We found a significant effect of the predictor beta values between groups (F(1,33) = 4.47, ). Importantly, within-subject effects for Condition (F(1,33) = 97.18, ) and Condition×Group (F(1,33) = 9.66, ) were present. To analyze the effect of the interaction Condition×Group, we used ANOVA between groups for each condition (Passive Viewing and Task). Results showed that the cortical responses during the task were higher for RP patients when compared to the control group (), whereas no differences were found for passive viewing condition between groups, thereby corroborating the main hypothesis.

3.3.2. Analyses by Subgroups Defined by Extent of Visual Field Loss

The statistical analysis described above was conducted to study differences among the two subgroups of patients according to the extent of visual field loss and the control group. Importantly, we found the effects for Condition (F(1,32) = 135.11, ), Condition×Group (F(2,32) = 13.86, p = 4.600 × 10-5), Zone×Ring (F(1,32) = 11.52, ), and Zone×Ring×Group (F(2,32) = 6.86, ). ANOVA for each condition (Passive Viewing and Task) was applied to analyze the interaction effect of Condition×Group. Results showed that the cortical responses for task were higher for RPsvf (small visual field) patients when compared to the RPlvf (large visual field) patients (corrected ) and the control group (corrected ), whereas no differences were found for passive viewing condition among the three groups. In this way, patients with smaller intact visual fields presented higher cortical response during the task than patients with larger visual field extent and controls. To study the effect of the interaction Zone×Ring×Group, we used ANOVA among groups for values of the ring (Ring1 and Ring2) in each cortical zone (FPZ and LPZ). We found increased cortical responses for RPsvf patients for Ring2 in the FPZ when compared to controls (corrected ). In this way, patients with smaller visual fields presented higher cortical response () than controls () in the FPZ during the stimulation of paracentral visual field, while patients with larger visual field extent () had similar activation to controls. These observations suggest that larger damage (present in small visual field patients) lead to larger allocation of visual attention mechanisms.

3.3.3. Analyses Based on Age of Onset Defined Subgroups

The statistical analysis described above was conducted to study differences among the two subgroups of patients according to the age of disease onset and the control group. Effects for Condition (F(1,32) = 90.76, ) and Condition×Group (F(2,32) = 5.35, ) were present. ANOVA for each condition (Passive Viewing and Task) was applied to analyze the interaction effect of Condition×Group. Results showed that the cortical responses for the task were surprisingly higher for RPlo (later onset) patients when compared to the control group (corrected ), whereas no differences were found for passive viewing condition among the three groups. In sum, patients with later disease onset presented surprisingly higher cortical response during task than controls, suggesting that lower activation may represent more “efficient” brain activity patterns in early onset patients (with more prolonged disease evolution and more time to recruit mechanisms where “efficiency” dominates).

3.4. Visual Memory Task: Responses in Retinotopic Regions (V1, V2, and V3)

In this section, we investigate how different retinotopic areas contribute to the observed patterns of activity (overall and according to the subgroups defined above).

3.4.1. All Group

The beta values of the task predictors (Ring1 Passive Viewing, Ring2 Passive Viewing, Ring1 Task, and Ring2 Task) were analyzed between groups for each condition (Condition Passive Viewing and Task), each ring (Ring Ring1 and Ring2), each visual area (Area V1, V2, and V3), each cortical visual region (Region Ventral and Dorsal), and each hemisphere (Hemisphere Left and Right). We found the main effects of Area (F(1.74,57.42) = 50.45, ), Condition (F(1,33) = 44.27, ), and Ring (F(1,33) = 53.81, ). We also found the effects of Area×Condition (F(1.76,58.20) = 25.22, ), Area×Region×Ring (F(1.72,56.92) = 9.12, ), and Area×Region×Ring×Group (F(1.72,56.92) = 5.98, ). We can summarize these interactions by the change in visual cortical responses from V1 to V3 (Area) between passive viewing and the one-back task condition (Condition) and between Ring1 (central visual field) and Ring2 (paracentral visual field) (Ring) for both groups.

We then hypothesized that age of onset and extent of visual field loss could influence the effects in patients.

3.4.2. Visual Field Subgroup

The statistical analysis described above was conducted to study differences between the two subgroups of patients accordingly to the extent of visual field loss and the control group. We found the effects for Area (F(1.82,58.10) = 33.96, ), Condition (F(1,32) = 58.91, ), Ring (F(1,32) = 48.37, ), Condition×Group (F(2,32) = 8.53, ), Area×Condition (F(1.80,57.65) = 20.06, ), Area×Hemisphere×Group (F(3.95,63.20) = 3.71, ), and Area×Region×Ring (F(1.78,57.00) = 11.04, ). An ANOVA for each condition (Passive Viewing and Task) was applied to analyze the critical interaction effect of Condition×Group. Results showed that the cortical responses for task were higher for RPsvf (small visual field) patients when compared to the control group (corrected ), whereas no differences were found for the passive viewing condition among groups. In this way, patients with smaller visual fields presented higher cortical response during task than controls, while patients with larger visual field extent showed responses that were similar to the control participants. To study the effect of the interaction Area×Hemisphere×Group, we used ANOVA between groups for each visual area (V1, V2, and V3) in each hemisphere (Left and Right). We found increased cortical responses for RPsvf patients in the right V1 when compared to controls (corrected ), but the remaining visual areas presented similar cortical activation. In this way, patients with smaller visual fields presented higher cortical response in the right primary visual cortex () when compared to controls (), while patients with larger visual field extent () showed similar activation to controls.

3.4.3. Subgroup Defined by Distinct Age of Onset

The statistical analysis described above was conducted to study differences among the two subgroups of patients according to the age of onset and the control group. We found the effects for Area (F(1.89,60.64) = 42.83, ), Condition (F(1,32) = 40.44, ), Ring (F(1,32) = 51.74, ), Area×Group (F(3.79,60.64) = 4.46, ), Area×Condition (F(1.82,58.365) = 20.28, ), Area×Hemisphere×Group (F(3.96,63.40) = 3.99, ), Region×Ring×Group (F(2,32) = 5.66, ), and Area×Region×Ring (F(1.79,57.38) = 11.51, ).

In sum, this analysis focused on retinotopic areas that essentially mimics the findings observed for the LPZ and FPZ zones, corroborating the main hypothesis of preferential attentional allocation in patients, and distinct effects of the visual field lesion extent and age of onset.

3.5. Cortical Thickness of Visual Areas (V1, V2, and V3)

Finally, we investigated whether functional changes were associated with structural alterations.

Visual cortical thickness differences were evaluated using repeated measures ANOVA with three within-subject factors (Area V1-V3, Region Ventral vs Dorsal, and Hemisphere Left Vs Right), one between-subject factor Group (RP vs. Control), with the average brain cortical thickness as a covariate to account for variability across participants. Cortical thickness of the individually defined visual areas was not different between the two groups (F(1,32) = 0.047, NS), and no within-subject effects or interactions were found. This was also true for subgroup analyses.

4. Discussion

We investigated whether visual cortical responses in a disorder of peripheral vision are related to recruitment of attentional mechanisms. To test this hypothesis, we used a visual one-back task and passive viewing conditions with a visual stimulus covering the central (Ring1) and paracentral (Ring2) visual field in a group of RP patients () with peripheral retinal loss and matched healthy controls (). Cortical responses were studied in visual retinotopic areas (V1, V2, and V3) and in two different regions of interest in V1: the FPZ representing the preserved visual field and the LPZ representing the visual field scotomata. To understand the influence of the level of peripheral degeneration and the disease age of onset, the analysis was further conducted for two distinct RP subgroups: subcategories defined by the extent of visual field loss (RPsvf, remaining small field, and RPlvf, remaining large field patients with bilateral visual field diameters under or over 9.50 deg, respectively), and subgroups defined by distinct ages of disease onset (RPeo, early onset, and RPlo, late onset—patients with age of onset of the disease lower or equal/greater than 14 years, respectively).

Our results demonstrated that RP patients have overall preserved visual cortical responses under central and paracentral visual field stimulation. Visual cortical responses (V1, V2, and V3) to the visual memory task stimuli were also overall preserved. A critical interaction with task condition was however found: RP patients presented higher overall cortical responses during the task condition than control participants in FPZ and LPZ regions. This was further highlighted when extent of visual loss was taken into account.

Concerning the role of the extent of visual loss, RPsvf patients, with smaller visual fields, presented higher overall visual cortical responses during the task condition than control participants, while responses for RPlvf patients, with larger visual field extent, were similar to control participants. Additionally, RPsvf patients had significantly higher cortical activation in the right V1 when compared to controls, while responses for RPlvf patients were similar to healthy participants. In line with these results, RPsvf patients presented higher overall cortical responses during the task condition than control participants and RPlvf patients in LPZ and FPZ regions, while responses for RPlvf patients were similar to control participants. Additionally, the cortical activation in the FPZ was higher for RPsvf patients during the stimulation of the paracentral visual field (Ring2) when compared to the healthy subjects. Because visual stimuli were the same in both task and passive viewing conditions, these responses for RPsvf patients seem to be related to attentional demands during the one-back task. Masuda et al. [13] found a similar increase in striate cortical responses of three RP participants related to changes in task demands and suggested that unmasking of feedback signals from the extrastriate cortex occurs when retinal signals are absent. Such unmasking might come from activation of previously silent synapses [32]. These feedback signals might be associated with attention, visual imagery, and task-related visual processing [26]. In our work, enhanced attentional top-down modulation may compensate for the lack of retinal input from the peripheral visual field in RPsvf patients with greater visual field loss [9, 24, 33]. A recent study presented evidence for increased functional connectivity between afferent early visual areas and cortical regions involved in visual processing (middle occipital gyrus and superior temporal gyrus/sulcus) in RP patients, suggesting a possible compensatory mechanism for peripheral visual loss. These authors also found enhanced functional connectivity between the deafferented visual cortex and higher-order regions (inferior parietal lobe/sulcus and middle frontal gyrus) involved in top-down control, attentional processes, and multisensory integration [34].

However, in the three patients reported in Masuda et al.’s work, the increased responses during task were found in the V1 LPZ [13]. Here, we report an overall increase in cortical responses under task demands while analyzing V1, V2, V3, and V1 FPZ and LPZ for patients with more severe visual field degeneration. In this way, feedback signals might influence both striate and extrastriate visual cortex when there is a severe lack of peripheral retinal input, not being restricted to the LPZ. A recent work from our group provided evidence of functional remapping in V1 in the same group of patients studied here. This functional reorganization was also more prominent in RP patient with larger visual field damage [6].

Previous works with macular degeneration also showed increased cortical responses for V1 LPZ while patients performed a one-back task with peripheral stimulation [1519, 24], contrary to passive viewing stimulation [19, 2125]. Some of the authors showed that these V1 responses were higher for more severe central retinal loss without foveal sparing [15, 18], in accordance with our work. Recently, Plank et al. [26] reported that patients with central scotomata presented enhanced cortical activation in areas beyond the retinotopic cortex for complex images with naturalist scenes, supporting an increased top-down modulation of the deprived visual cortex. This result was further supported by the work of Sabbah et al. [34] showing increased functional connectivity between the LPZ and high-level regions in central retinal disease patients.

Several studies with glaucoma patients found reduced amplitude of cortical responses in V1 during passive viewing stimulation associated with structural damage of the optic disk, the RNFL thickness, and/or the visual field scotomata [3539]. A recent study demonstrated reduced cortical activity within the LPZ in V1 and V2 in glaucoma patients under passive viewing [40]. Here, we did not find decreased activity during passive viewing in the cortical regions studied, which might indicate that RP patients have preserved visual cortical responses even for severe visual field damage, possibly due to a similar compensation mechanism by increased attentional modulation.

The cortical responses in visual areas were not globally significantly different among the RPeo and RPlo patients’ subgroups with different onset ages and the control participants. Thus, overall visual cortical responses do not seem to be influenced by disease onset age, in contrast to the extent of visual loss. However, RPlo patients with later disease onset presented higher overall cortical responses during the task condition than control participants in FPZ and LPZ regions, while responses for RPeo patients with earlier forms of the disease were similar to control participants. This result was unexpected considering our initial hypothesis that earlier onset ages would lead to larger brain alterations. Two sorts of mechanisms might be operating: the first requiring long-term circuit modifications and, more present in early-onset patients, the second entailing stronger frontoparietal recruitment which tends to manifest more in patients with more recent changes in visual experience. There is indeed evidence showing that higher top-down modulation may indeed be stronger in participants with more recent changes in visual experience [41, 42]. Moreover, it may reflect the fact that longer disease durations may lead to efficient compensatory mechanisms and decrease of frontoparietal activation. A second study from Rosa et al. [41, 42] shows that frontoparietal activation decreases over time as patients’ vision becomes more adapted. Less fMRI activation might actually indicate “more efficient” compensation [41, 42]. Studies in macular degeneration patients did not report the effects of the age of onset on visual cortical responses [15, 22], while others showed that juvenile-macular degeneration patients with earlier disease forms have stronger cortical activation than age-related macular degeneration patients with later onset age [19]. Future studies should address the discrimination between age-dependency of neuroplasticity and disease-duration effects, which are separable. Disease onset age is often difficult to determine, which can make this an imprecise measure of RP severity [11, 43].

Given the evidence for cortical reorganization in our prior study [6], remodeling at the retinal level is unlikely. This issue can be further clarified in the future by explicitly computing population receptive fields or alternatively running experiments with artificial scotomata.

Our results did not find evidence for visual cortical structural alterations in this cohort of RP patients, showing that the visual loss level was not sufficient to produce significant cortical atrophy in the visual areas studied (V1 to V3). To our knowledge, few structural MRI studies have been conducted with low vision RP patients [3, 44]. In our study, the RP patients did not present RNFL thickness atrophy which is in line with the preservation of visual cortical thickness. Nonetheless, in more advanced stages of RP disease with larger photoreceptor loss and retinal ganglion cell degeneration, disuse-driven mechanisms may lead to the visual cortical atrophy pattern that is often seen in macular degeneration, glaucoma, and also late-blindness.

5. Conclusion

We found that cortical visual areas (V1, V2, and V3) responses under attentional demands were increased in patients with larger degeneration of visual field. Moreover, activation during the task condition was increased for patients in both cortical regions corresponding to the preserved (FPZ) and the damaged visual field (LPZ), specifically for patients with severe visual field loss. These findings were identified in the presence of preserved visual cortical structure. The age of onset of the disease did not seem to be associated with visual cortical alterations. We conclude that RP patients may have relatively preserved visual cortical responses due to feedback attentional modulation in the absence of cortical atrophy, despite their retinal degeneration. The unmasking of corticocortical feedback signals from higher level visual regions involved in attentional processes might explain the increased cortical responses [1]. Such unmasking might lead to activation of previously silent synapses.

These results might be considered in the context of strategies for treating retinal diseases [21, 45, 46]. This is quite relevant given previous evidence that attentional cueing improves vision restoration therapy in patients with visual field loss [47, 48]. The role of higher-level neuronal networks [49] and their functional connectivity [50] cannot be underestimated in this context. This suggests that visual responses can be dynamically adapted as a function of flexible mechanisms requiring the interaction between high-level regions that implement attentional control.

Abbreviations

RP:Retinitis Pigmentosa
MRI:Magnetic resonance imaging
fMRI:Functional magnetic resonance imaging
RNFL:Retinal nerve Fiber layer
LPZ:Lesion projection zone
FPZ:Function projection zone
deg:Degrees.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Disclosure

This work was presented in an abstract form in the 2017 European Association for Vision and Eye Research Conference.

Conflicts of Interest

The authors declare that no conflicts of interest exist.

Acknowledgments

This work was supported by the Portuguese Funding Agency for Science and Technology (FCT) grants E-Rare2-SAU/0001/2008, E-Rare4/0001/2012, COMPETE, POCI-01-0145-FEDER-007440, FCT. UID/NEU/04539/2013–2020, MEDPERSYST, POCI-01-0145-FEDER-016428, and BIGDATIMAGE, CENTRO-01-0145-FEDER-000016 financed by Centro 2020 FEDER, COMPETE. We would like to thank Carlos Ferreira, João Marques, and Sónia Afonso for technical assistance in the magnetic resonance acquisitions. We express our thanks to the patients and their families and the other participants for their collaboration.