Abstract
Inherited retinal degenerative diseases are a prominent cause of blindness. Although mutations causing death of photoreceptors are mostly known, the pathophysiology downstream in the inner retina and along the visual pathway is incompletely characterized in the earliest disease stages. Here, we investigated retinal, midbrain and cortical visual function using electroretinography (ERG), the optomotor response (OMR), visual evoked potentials (VEPs), respectively, and single unit electrophysiology at the primary visual cortex (V1) in light-adapted juvenile (approximately one-month-old) and young adult (three-month-old) RhoP23H/WT mice, representative of early-stage retinitis pigmentosa (RP). Photopic ERG revealed up to ∼30% hypersensitivity to light in RhoP23H/WT mice, as measured by the light intensity required to generate half-maximal b-wave (I50 parameter). RhoP23H/WT mice also showed increased OMRs toward low spatial frequency (SF) drifting gratings, indicative of visual overexcitation at the midbrain level. At the V1 level, VEPs and single-cell recordings revealed prominent hyperexcitability in the juvenile RhoP23H/WT mice. Mean VEP amplitudes for light ON stimuli were nearly doubled in one-month-old RhoP23H/WT mice compared with controls, and more than doubled for light OFF. Single-cell recordings showed a significantly increased spontaneous V1 neuron firing in the RhoP23H/WT mice, and persistent contrast and temporal sensitivities. In contrast, direction selectivity was severely compromised. Our data suggest that during early RP, the visual pathway becomes hyperexcited. This could have both compensatory and deleterious consequences for visual behavior. Further studies on the mechanisms of hyperexcitability are warranted as this could lead to therapeutic interventions for RP.
Significance Statement
Lost retinal function in many blinding retinal degenerative disorders could soon be alleviated by advanced therapies that restore photoreception. However, it is unknown whether a visual system rewired downstream of the photoreceptors can process signals adequately. We studied the functional consequences of early rod death along the visual pathway in young retinitis pigmentosa (RP) mice. Photopic inner retina responses were moderately hypersensitized in the electroretinograms (ERG) of RP mice. Reflex-based visual behavior and visual cortex electrophysiology showed hyperexcitability. Some aspects of complex visual processing were remarkably resistant to degeneration, whereas others were severely impacted. We conclude that the visual system adapts to lost photoreception by increasing sensitivity, but simultaneously becomes detrimentally hyperexcited. Mechanistic understanding could lead to therapeutic preservation and restoration of vision.
Introduction
Consequences of total or near total loss of outer retinal photoreceptors, rods and cones, to the inner retina function have been widely studied using animal models of retinal degenerative diseases (D’Orazi et al., 2014; Strettoi, 2015; Pfeiffer et al., 2020; Strettoi et al., 2022). With respect to clinical correlates, corresponding anatomic data have been obtained from patient retinas postmortem (Jones et al., 2016). Such studies have demonstrated robust rewiring of neural connections, commonly referred to as “remodeling,” leading to dramatically increased spontaneous neural activity, decreased signal-to-noise ratios, and attenuated light responses in the inner retina (Toychiev et al., 2013; Trenholm and Awatramani, 2015; Haselier et al., 2017; Telias et al., 2020). Much concern has been raised that retinal remodeling may preclude restoration of visual function even if augmentation of photoreception by various advanced therapies could be achieved (Marc et al., 2014; Reh, 2016; Foik et al., 2018; Suh et al., 2020). By default, this refractoriness would be particularly problematic in the adult organism with limited capability for neural plasticity and presumed inability to adapt to and properly process the restored sensory signals.
Much less attention has been directed to investigating functional outcomes of early retinal degeneration and milder disease stages, although most patients with retinal degeneration suffer from a partial loss of vision, and some may never become legally blind (Berson et al., 2002; Xu et al., 2020). Studies focusing on early disease stages could help us understand the triggers and mechanisms of retinal remodeling and lead to increased ability to develop vision restoration interventions. Equally important, increased mechanistic understanding could help us design therapies to increase quality of life in patients with partial loss of vision because of retinal degeneration.
Most inherited retinal degenerations are classified under the umbrella term retinitis pigmentosa (RP). In a typical RP case, primarily rod photoreceptors are affected first; consequently the disease first manifests as loss of night vision and shrinking of the visual field (Hartong et al., 2006). Central and color vision impairment manifest only at a much more advanced disease state. RP phenotype has been studied extensively with respect to the primary insult, rod degeneration and dysfunction. Investigation of the cone-associated phenotype has gained somewhat less attention, although the secondary cone degeneration may be a more realistic target for conventional pharmacological interventions. In fact, preservation of cone-pathway function could retain the most important attributes of vision for humans, such as color vision and high visual acuity. Another gap in knowledge exists in how visual signals are modulated in the brain during early RP. This is crucial as the complex circuitry in the inner retina, which ultimately sends the visual signals to the rest of the brain, rewires early in RP (Soto and Kerschensteiner, 2015). In the current study, we addressed these issues by using an established animal model of autosomal dominant RP (heterozygote RhoP23H/WT mice) before cone degeneration, and by recording photopic electroretinography (ERG), behavioral optomotor responses (OMRs), and primary visual cortex (V1) responses to light ON-OFF and multimodal pattern stimuli.
Materials and Methods
Animals
A mouse model of autosomal dominant RP was used in this study. The mouse carries a P23H mutation in the rhodopsin gene, causing a phenotype of extremely fast rod photoreceptor degeneration in homozygous mutants (RhoP23H/P23H), and an intermediately fast progressing rod degeneration in heterozygous mutants (RhoP23H/WT; Sakami et al., 2011, 2014). This report mainly focuses on the heterozygous mutants. As is common for the RP phenotype (Hartong et al., 2006), cone photoreceptor degeneration is much delayed in RhoP23H/WT mice (Sakami et al., 2011, 2014). The RhoP23H/WT mice were generated by crossbreeding RhoP23H/P23H mice with wild-type (WT) C57BL/6J mice (The Jackson Laboratory, stock # 000664). Age-matched WT mice were used as controls. Both male and female mice were used that were group-housed in a standard vivarium using ad libitum feeding. The light-dark cycle was set at 12 h/12 h (lights on 6:30 A.M., lights off 6:30 P.M.). Cone-transducin knock-out mice (Gnat2−/−, a kind gift of Marie Burns, UC Davis), that lack cone-mediated function but do not show retinal degeneration (Ronning et al., 2018), were used for calibration of photopic ERG background (Extended Data Fig. 1-1). In all procedures, animal subjects were treated following the NIH guidelines for the care and use of laboratory animals, and the ARVO Statement for the Use of Animals in Ophthalmic and Vision Research; and under a protocol approved by the Institutional Animal Care and Use Committee of UC Irvine (AUP-18-124).
Extended Data Figure 1-1
Steady red illumination suppresses rod-driven ERG in cone-transducin knock-out (Gnat2–/–) mice. Gnat2–/– mice were dark-adapted overnight before the recording. Mice were first recorded in a dark-adapted state (no background), and then steady red light at 10, 100, 200, and 500 cd/m2 was introduced, and ERGs were repeated in an increasing “red light-adapted” state. Each step had an adaptation period of 60 s. Two stimuli were presented for each step using an interstimulus interval of 90 s. Note that the ERG response in 200 cd/m2 background illumination is less than 10% from the dark-adapted state, and the response in 500 cd/m2 background illumination cannot be distinguished from the noise. Download Figure 1-1, TIF file.
ERG
The ERG was performed under standard laboratory lighting conditions using a Diagnosys Celeris rodent ERG device (Diagnosys), with some modifications from a previous protocol (Orban et al., 2018). The mice were anesthetized with ketamine (100 mg/kg, KetaVed; Bioniche Teoranta, Inverin Co) and xylazine (10 mg/mg, Rompun; Bayer) by intraperitoneal injection, and their pupils were dilated with 1% tropicamide (Tropicamide Ophthalmic Solution USP 1%; Akorn), and thereafter kept moist with 0.3% hypromellose gel (GenTeal; Alcon). Light stimulation was produced by an in-house scripted simulation series in Espion software (version 6; Diagnosys). The eyes were stimulated with a green light-emitting diode (LED; peak 544 nm, bandwidth 160 nm) or with a UV LED (peak emission 370 nm, bandwidth 50 nm). While green stimulation was used, the steady rod-suppressing background light consisted of 200 cd/m2 red (peak 630 nm, bandwidth 100 nm) and 100 cd/m2 UV. The background during UV stimulation consisted of 200 cd/m2 red and 100 cd/m2 green. Steady red light at 200 cd/m2 suppresses the ERG signal practically fully in Gnat2−/− mice (Extended Data Fig. 1-1). However, a steady strong UV or green background at 100 cd/m2 was further added to facilitate rod suppression during monochromatic green or UV flash stimulation, respectively. Seven different green light stimulus intensities between 0.64 and 214 mW/sr/m2, and eight different UV intensities between 1.53 and 509 mW/sr/m2, were used in ascending order. The whole stimulation protocol lasted <10 min. Stimulus intensities in radiance units were obtained from measurements and conversion coefficients provided by Diagnosys LLC. The details of conversion of luminous energy units to radiance units are presented in Table 1.
The ERG signal was acquired at 2 kHz and filtered with a low-frequency cutoff at 0.25 Hz and a high-frequency cutoff at 300 Hz. Espion software automatically detected the ERG a-wave (first negative ERG component) and b-wave (first positive ERG component) amplitudes; a-wave amplitude was measured from the signal baseline, whereas b-wave amplitude was measured as the difference between the negative trough (a-wave) and the highest positive peak. For assessment of retinal light sensitivity (I1/2), b-wave amplitudes were fitted as a function of stimulus intensity by using the Naka–Rushton equation (Naka and Rushton, 1966), wherein the I1/2 parameter describes the light intensity needed to exert half-maximal b-wave.
OMR behavioral vision test
The OMRs were assessed using a commercial OMR platform (Phenosys qOMR, PhenoSys GmbH) that utilizes automated head tracking and behavior analysis, following protocols described elsewhere (Suh et al., 2020). The OMR arena was lit at ∼80 lux, corresponding to the photopic light level. Rotating (12° s−1) vertical sinusoidal grating stimuli at various spatial frequencies (SFs) were presented for 11 min per trial to light-adapted mice. The contrast between the white and black gratings was set at 100%, whereas the SF (0.05, 0.1, 0.15, 0.20, 0.25, 0.30, 0.35, 0.375, 0.40, 0.425, 0.45 cycles per degree of visual angle; CPD) pattern changed every 60 s in a random order, with one exception; each session always started with 0.1 CPD to facilitate acclimatization to the task, as this SF is known to evoke reliable OMR in WT mice. Each mouse was tested in at least four trials. The performances across the trials were averaged for analysis, excluding those 60-s stimulus periods that led to a correct/incorrect ratio smaller than 0.8. All experiments were performed in the morning before noon.
Single-unit and local field potential (LFP) recordings and visual stimulation
Mice were initially anesthetized with 2% isoflurane in a mixture of N2O/O2 (70%/30%), then placed into a stereotaxic apparatus. A small, custom-made plastic chamber was glued (Vetbond) to the exposed skull. After 1 d of recovery, re-anesthetized animals were placed in a custom-made hammock, maintained under isoflurane anesthesia (1–2% in N2O/O2), and multiple single tungsten electrodes were inserted into a small craniotomy above the visual cortex. Once the electrodes were inserted, the chamber was filled with sterile agar and sealed with sterile bone wax. During recording sessions, animals were sedated with chlorprothixene hydrochloride (1 mg/kg, i.m.; Camillo et al., 2018) and kept under light isoflurane anesthesia (0.2–0.4% in 30% O2). EEG and EKG were monitored throughout the experiments, and body temperature was maintained with a heating pad (Harvard Apparatus).
Data were acquired using a 32-channel Scout recording system (Ripple). The LFP from multiple locations was bandpass filtered from 0.1 to 250 Hz and stored together with spiking data on a computer with a 1-kHz sampling rate. The LFP signal was cut according to stimulus time stamps and averaged across trials for each recording location to calculate visually evoked potentials (VEPs; Foik et al., 2015; Kordecka et al., 2020; Suh et al., 2020; Lewandowski et al., 2022). The spike signal was bandpass filtered from 500 Hz to 7 kHz and stored in a computer hard drive at a 30-kHz sampling frequency. Spikes were sorted online in Trellis (Ripple) while performing visual stimulation. Visual stimuli were generated in MATLAB (MathWorks) using Psychophysics Toolbox (Brainard, 1997; Pelli, 1997; Kleiner et al., 2007) and displayed on a γ-corrected LCD monitor (55 inches, 60 Hz; 1920 × 1080 pixels; 52 cd/m2 mean luminance). Stimulus onset times were corrected for LCD monitor delay using a photodiode and microcontroller (in-house design; Foik et al., 2018).
The vision quality was assessed using protocols published in our previous work (Foik et al., 2018, 2020; Kordecka et al., 2020; Suh et al., 2020). For recordings of visually evoked responses, cells were first tested with 100 repetitions of a 500-ms bright flash of light (105 cd/m2). Receptive fields for visually responsive cells were then located using square-wave drifting gratings, after which optimal orientation/direction, SF and temporal frequency (TF) were determined using sine-wave gratings. SFs tested were from 0.001 to 0.5 cycles/°. TFs tested were from 0.1 to 10 cycles/s. With these optimal parameters, size tuning was assessed using apertures of 1–110° at 100% contrast. With the optimal size, temporal and SF, and at high contrast, the orientation tuning of the cell was tested again using eight orientations × two directions each, stepped by 22.5° increments. This was followed by testing contrast. Single units were recorded along the whole depth of the primary visual cortex (V1) as presented in previous studies (Foik et al., 2018; Suh et al., 2020; Frankowski et al., 2021; Choi et al., 2022; Lewandowski et al., 2022).
VEP data analysis
The response amplitude of LFP was calculated as a difference between the peak of the positive and negative components in the VEP wave. The response latency was defined as the time point where maximum response occurred. The maximum of the response was defined as the maximum of either the negative or positive peak.
Tuning curves were calculated based on the average spike rate. Optimal visual parameters were chosen as the maximum response value. Orientation tuning was measured in degrees as the half-width at half-height (HWHH; 1.18 × σ) based on fits to Gaussian distributions (Carandini and Ferster, 2000; Alitto and Usrey, 2004; Foik et al., 2018, 2020) using:
Size tuning curves were fitted by a difference of Gaussian (DoG) function:
The optimal SF and TF were extracted from the data fitted to Gaussian distributions using the following equation (DeAngelis et al., 1993; Van Den Bergh et al., 2010; Foik et al., 2018, 2020):
The contrast tuning was fitted by using the Naka–Rushton equation (Naka and Rushton, 1966; Albrecht and Hamilton, 1982; Van Den Bergh et al., 2010; Przybyszewski et al., 2014; Foik et al., 2018):
Spectral dissection of ERG and VEP responses
The ERG and VEPs were filtered using fourth order Butterworth bandpass filters of 13–29 Hz for β, and 30–120 Hz for γ ranges. Two filtered signals for each animal were then used to calculate amplitude spectra using the fast Fourier transform (FFT) algorithm (MATLAB) and amplitude spectrum area (AMSA), and shown as bar plots. The dominant frequencies in β and γ ranges were picked for each animal separately as the frequency at maximum amplitude in the calculated spectrum. This process was repeated for each animal in each group and for green and UV light intensities to create the frequency versus light intensity and AMSA versus light intensity traces.
Statistical analysis
Data throughout the manuscript are presented as mean ± SEM. The level of statistical significance was set at p < 0.05. The ERG light intensity series were analyzed by two-way repeated measures ANOVA (RM ANOVA), using the light intensity as the within-subjects factor (repeated) and genotype as the between-subjects factor. Sidak’s post hoc tests were applied if ANOVA found a significant between-subjects main effect, or within between-subjects interaction. Retinal light sensitivity (I1/2) was analyzed by Welch’s t test. In the OMR data, the genotypes were compared using regular two-way ANOVA, whereas repeated measurements in the same animals (e.g., WT animals at one vs three months of age) were analyzed using two-way RM ANOVA by age as the repeated factor. GraphPad Prism 9 software was used for ERG and OMR data analysis. Data from cortical recordings and oscillatory activity was analyzed with the Kruskal–Wallis test followed by Dunn–Sidak post hoc tests. Cortical recording and oscillatory activity offline data analysis and statistics were performed in MATLAB (MathWorks).
Results
Increased photopic ERG b-wave sensitivity in RhoP23H/WT mice
We began our investigations by recording photopic ERGs in light-adapted one-month-old mice (Fig. 1A–H). Because mouse S- and M-cones are most sensitive for UV and green light, with peak sensitivities at around 360 and 510 nm (Jacobs et al., 1991; Lyubarsky et al., 1999), respectively, and since S- and M-cones could be differentially affected during retinal degeneration (Greenstein et al., 1989; Schwarz et al., 2018; Hassall et al., 2020), we recorded photopic ERGs using both monochromatic UV light and green light stimuli. The a-wave or b-wave response amplitudes did not show an overall difference between genotypes as determined by RM ANOVA (a-wave, Extended Data Fig. 1-2A,B; b-wave, Fig. 1B,F). However, anecdotally RhoP23H/WT mice seemed to display higher amplitudes at intermediate light intensities (Fig. 1B,F). Therefore, we used the Naka–Rushton function for curve fitting (Naka and Rushton, 1966) and derivatization of light intensity required to generate half-maximal b-wave responses (I1/2) as a measure of retinal light sensitivity. We found decreased I1/2, indicative of increased sensitivity, in RhoP23H/WT mice for both green light (Fig. 1C) and UV light (Fig. 1G) stimulation, which supports previous findings (Leinonen et al., 2020). In addition, b-wave peak latency on average was faster in RhoP23H/WT mice as compared with WT mice (Fig. 1A,H), but this trend toward shorter latency did not reach a statistically significant value (green, p = 0.09; UV, p = 0.20). Instead, the a-wave latency was significantly faster in response to a UV flash in one-month-old RhoP23H/WT mice compared with WT (Extended Data Fig. 1-2F), but not in response to a green flash (Extended Data Fig. 1-2E, p = 0.73).
Extended Data Figure 1-2
Photopic ERG a-wave amplitude and latency graphs. A, B, E, F, Data from one-month-old RhoP23H/WT mice (blue, n = 11) and WT mice (black, n = 8). C, D, G, H, Data from three-month-old RhoP23H/WT mice (red, n = 11) and WT mice (green, n = 11). The RM ANOVA between-subjects effect revealed that the a-wave peak occurred significantly faster in one-month-old RhoP23H/WT mice compared to WT (graph F; RM ANOVA between-subjects effect: p < 0.01). There were no statistically significant differences for any of the other conditions as measured by RM ANOVA. Download Figure 1-2, EPS file.
Extended Data Figure 1-3
Group-averaged ERG waveforms for every stimulus condition for each animal group. A, Green flash responses. B, UV flash responses. Download Figure 1-3, EPS file.
To study whether increased photopic ERG sensitivity is persistent in adulthood, we performed the same recordings at three months of age (Fig. 1I–P). Here, the maximal b-wave amplitudes in RhoP23H/WT mice tended to remain slightly lower compared with WT (Fig. 1J,N), although no significant RM ANOVA effects were found for the amplitude comparisons (green, p = 0.25; UV, p = 0.76). Instead, RM ANOVA revealed a significant light intensity-response amplitude interaction both in response to green light (Fig. 1J, p < 0.05) and UV light (Fig. 1N, p < 0.001) stimulation. The RM ANOVA interaction together with visual inspection of light intensity-response amplitude curves (Fig. 1J,N) indicate that although maximal b-wave amplitudes tended to decrease in RhoP23H/WT mice, responses to intermediate light intensities were above WT level. Indeed, RhoP23H/WT mice still showed significantly decreased I1/2 b-wave sensitivity parameters compared with WT mice (Fig. 1K,O), regardless of intermediate retinal degeneration at this stage (Sakami et al., 2011; Leinonen et al., 2020). We interpret the relative decrease in b-amplitude and increase in b-wave latency in RhoP23H/WT mice between one and three months of age (compare Fig. 1B,D,F,H and J,L,N,P) to reflect a disease stage when cone degeneration is commencing, although there were not yet decreases in a-wave amplitudes at this age range (Extended Data Fig. 1-2). Overall, the modest photopic ERG changes which we observed in RhoP23H/WT mice appeared to be more evident in the S-cone pathway as compared with the M-cone pathway.
To provide a context for how the advanced-disease stage represents a cone-mediated retinal function, we recorded photopic ERGs also in one-month-old homozygous mutant (RhoP23H/P23H) mice with an almost complete loss of photoreceptors by this age (Sakami et al., 2011). These recordings demonstrated only a residual response in the cone-pathway (Fig. 1A,B,E,F; Extended Data Fig. 1-3), which rendered curve fitting, light sensitivity, and response latency assessments impractical.
Prominent oscillatory potentials (OPs) ride on the ascending phase of the rodent ERG b-wave and could have implications on the b-wave analyses. We therefore performed OP amplitude and peak frequency analyses at two distinct ranges: β (13–29 Hz) and γ (30–120 Hz; Extended Data Fig. 1-4). The OP amplitude was significantly attenuated in RhoP23H/P23H mice as expected, because of extreme retinal degeneration and resulting minimal input from photoreceptors. Instead, there were no significant differences between the RhoP23H/WT and WT mice either in amplitude spectrum area (AMSA) or dominant frequency parameter at any of the conditions, except for one: RhoP23H/WT mice displayed a higher dominant frequency than WT mice in response to the strongest UV flash at the β range (Extended Data Fig. 1-4). Overall, OPs cannot explain the increased b-wave sensitivity observed in RhoP23H/WT mice (Fig. 1).
Extended Data Figure 1-4
Spectral dissection of oscillatory potentials from the ERG responses. A, Filtered ERG response at β range in response to 214 mW·sr/m2 green flash. B, Amplitude spectra corresponding to signals from A. C, The AMSA calculated from B. D, The dominant frequency from response in A. E, AMSA at all green flash intensities. F, Dominant response frequencies at all green flash intensities (G) Filtered ERG response at β range in response to 509 mW·sr/m2 UV flash. H, Amplitude spectra corresponding to signals from G. I, The AMSA calculated from H. J, The dominant frequency from the response in G. K, AMSA at all UV flash intensities. L, Dominant response frequencies at all UV flash intensities. M, Filtered ERG response at γ range in response to 214 mW·sr/m2 green flash. N, Amplitude spectra corresponding to signals from M. O, The AMSA calculated from N. P, The dominant frequency from response in M. Q, AMSA at all green flash intensities. R, Dominant response frequencies at all green flash intensities. S, Filtered ERG response at γ range in response to 509 mW·sr/m2 UV flash. T, Amplitude spectra corresponding to signals from S. U, The AMSA calculated from T. V, The dominant frequency from the response in S. W, AMSA at all UV flash intensities. X, Dominant response frequencies at all UV flash intensities. Statistical analysis was performed by the Kruskal–Wallis test followed by Dunn–Sidak post hoc tests: *p < 0.05, **p < 0.01, ***p < 0.001. Download Figure 1-4, EPS file.
RhoP23H/WT mice display increased optomotor movements toward low SF drifting patterns
We next sought to determine whether visual acuity, an elementary visual behavior as measured by OMRs, is changed in light-adapted/photopic conditions in 1.5- and 3-month-old RhoP23H/WT mice. An earlier report using OMR found contrast sensitivity was well-preserved in RhoP23H/WT mice even in scotopic conditions despite prominent rod death (Leinonen et al., 2020). Here, we used a commercial OMR setup coupled with an automated head-tracking system (Kretschmer et al., 2013, 2015), which allowed us to record the mouse optomotor movements with respect to the drifting pattern in its visual field at 30 frames per second. The common OMR studies have typically investigated the SF (visual acuity assessment) or contrast threshold detection limits using subjective assessment. In contrast, our paradigm ranks correct/incorrect optomotor movements at every frame, providing a more sensitive and objective means to document “supernormal” visual behavior. At the highest SF (smallest pattern sizes) the stimuli elicited similar responses in the 1.5-month-old RhoP23H/WT mice and WT mice; however, with the RhoP23H/WT mice, optomotor movements toward stimuli at the lowest SF (largest pattern size) were more frequent than those of WT mice, as demonstrated by their higher OMR index (Fig. 2A). This behavior reflects sensorimotor hypersensitivity to low SF drifting patterns in the visual field of the RhoP23H/WT mice. The same mice performed the same task later at three months of age. Head-tracking behavior in RhoP23H/WT mice slightly decreased between 1.5 and 3 months of age (Fig. 2B), indicating an initial decline in visual acuity in three-month-old RhoP23H/WT mice. Still, RhoP23H/WT mice at three months of age continued with a supernormal OMR at low SF, as compared with same-age WT mice (Fig. 2C; OMR index-SF interaction, p < 0.05). Intact responses in WT mice between 1.5 and 3 months of age document the feasibility of the repeated analysis (Fig. 2D, p = 0.49).
Stimulus-dependent and spontaneous visual cortex hyperexcitability in juvenile RhoP23H/WT mice
Most importantly for central vision, it was important to investigate the consequences of early retinal degeneration for light responses at the V1 level (note, OMR is a subcortically originating reflex and does not require V1 function; Douglas et al., 2005). To this end, we started by recording VEPs in response to light ON/light OFF stimuli. We found an increase in both ON and OFF responses in one-month-old RhoP23H/WT mice compared with same-age WT mice (Fig. 3A–C,E), indicating hyperexcitability. Conversely, the same age homozygous RhoP23H/P23H mutants showed only a residual response as expected, because of severe retinal degeneration and dysfunction (Fig. 3C,E; Sakami et al., 2011). Heterozygous RhoP23H/WT mice at three months old, when more than half of their photoreceptors are gone (Leinonen et al., 2020), still showed VEP amplitude comparable to WT mice (Fig. 3C,E). The response latency to the ON stimulus increased with disease severity (Fig. 3D, p < 0.01), whereas latency to the OFF stimulus remained intact in heterozygous mutants (Fig. 3F).
To investigate whether oscillatory activity is also affected in V1 in response to ON and OFF stimuli, we analyzed AMSA as well as dominant oscillation frequency at β (13–29 Hz), and γ frequency (30–120 Hz) ranges (Fig. 4A,B and E,F, respectively). As with the main VEP component amplitudes, we found on average higher AMSAs in the β (Fig. 4B,C) and γ (Fig. 4F,G) frequency ranges in one-month-old RhoP23H/WT mice, although far from statistical significance (p = 0.33 and p = 0.46, respectively). Instead, by three months of age, the amplitude of oscillatory activity lowered substantially in RhoP23H/WT mice and reached levels similar to WT (Fig. 4B,C; p = 0.79). As expected, AMSAs at both β and γ ranges were significantly lower in the most advanced disease stage studied, in RhoP23H/P23H mice (Fig. 4B,C,F,G). We found no differences in dominant oscillation ranges at either the β (Fig. 4D, p = 0.13) or γ (Fig. 4H, p = 0.3) ranges. Overall, the lack of major differences in dominant oscillatory frequencies suggests no change in large network processing in the brain.
Single-unit recordings corroborated significant visual cortex hyperexcitability in the juvenile RhoP23H/WT mice (Fig. 5). Both light ON and light OFF conversions evoked more spiking activity in one-month-old RhoP23H/WT mice relative to WT (Fig. 5A,B,G,H). One should note, however, that more than half of the V1 neurons in RhoP23H/WT mice did not respond to the light stimulus, whereas in WT mice <10% of units were nonresponsive (Fig. 5E). The cortical hyperexcitability in juvenile RhoP23H/WT mice was not only stimulus-dependent but also continuous as the stimulus-independent background spiking was also significantly elevated in the juvenile RhoP23H/WT mice versus WT (Fig. 5A,B,F). The spiking activity diminished in the mutant mice with disease progression, so that average spiking activity to light ON or OFF did not differ between three-month-old RhoP23H/WT and WT mice (Fig. 5A,C,G,H). Notably, however, light ON responses in RhoP23H/WT mutants were qualitatively abnormal regardless of age, as highlighted by the second peak in spiking activity occurring at ∼200 ms or later from light onset (Fig. 6C), a phenomenon that was absent in the WT mice (Fig. 5A). The spiking pattern at light OFF remained relatively normal in RhoP23H/WT mutants (Fig. 5A–C). As expected, the spiking activity to both light ON or light OFF was significantly dampened in the homozygous RhoP23H/P23H mice (Fig. 5D,G,H). Conversely, the spontaneous spiking activity was increased in late-stage disease (RhoP23H/P23H) compared with the less advanced disease stage (three-month-old RhoP23H/WT; Fig. 5F), indicating that some phenomenon independent of photoreceptors begins to drive spontaneous V1 neuron spiking in late-stage retinal degeneration, similar to what has been proposed previously regarding activity recordings of cone bipolar, amacrine, and retinal ganglion cells (RGCs; Borowska et al., 2011; Menzler and Zeck, 2011).
Compromised tuning properties of V1 neurons to drifting pattern stimuli
We last investigated V1 neuron tuning properties in RhoP23H/WT mice. For this task, we used drifting sinusoidal-pattern stimuli to study several different parameters relevant for visual behavior: orientation/direction tuning, SF tuning, as well as thresholds for SF, contrast, and TF. Several abnormalities were observed in RhoP23H/WT mice. First, the orientation tuning curve was distinctly broader compared with WT mice already at one month of age (Fig. 6A,B). Strikingly, the three-month-old RhoP23H/WT mice practically preferred no directions, demonstrating a dramatic loss of directional selectivity. Next, the pattern size to which mice optimally responded increased with disease progression, indicating enlargement of receptive fields (Fig. 6C,D). The optimal SF was slightly decreased at one month old for RhoP23H/WT mice, and severely decreased by the age of three months (Fig. 6E,F). In contrast, the optimal TF was abnormally high in one-month-old RhoP23H/WT mice (Fig. 6G,H), mirroring hastened kinetics of retinal responses to light stimulation (Fig. 1D,H; Extended Data Fig. 1-2F). Also, the contrast sensitivity in RhoP23H/WT mice remained at similar levels to those for WT mice, even at three months of age (Fig. 6I,J). Altogether, the VEP and V1 single-unit recordings observed in young RhoP23H/WT mice reveal that hastened ERG kinetics (Fig. 1D,H; Extended Data Fig. 1-2F; Pasquale et al., 2021), maintained behavioral contrast and TF thresholds (Leinonen et al., 2020; Pasquale et al., 2021), and increased preference for lower SF stimuli (Fig. 2A,C) are reflected in V1 function. We were unable to record meaningful responses to pattern stimuli in the most advanced disease state, in RhoP23H/P23H mice.
Discussion
Consequences of early rod degeneration to photopic visual function
Previous studies have suggested paradoxical sensitization of cone-mediated retinal function in the earliest stages of rod degeneration in mice (Leinonen et al., 2017, 2020). To investigate this phenomenon further, we recorded photopic ERG responses in RhoP23H/WT mice at one and three months of age, when ∼20% and ∼60% of rods have degenerated, respectively; but cone count remains normal (Leinonen et al., 2020). Maximal photopic ERG responses were intact in one-month-old RhoP23H/WT mice, indicating a lack of cone degeneration at this age (Fig. 1; Extended Data Fig. 1-2). Instead, the sensitivity of the b-wave was enhanced compared with WT controls. At three months of age, the increase in b-wave sensitivity became even more apparent in RhoP23H/WT mice despite their maximal b-wave amplitude showing some trend toward a decline (Fig. 1). In a more advanced stage of RP, the cone function undoubtedly declines (Sakami et al., 2011; Leinonen et al., 2020). Increased b-wave sensitivity in rod-specific function in RhoP23H/WT/Gnat2−/− mice have been recently observed, regardless of prominent rod degeneration (Leinonen et al., 2020). Relatively better-preserved ERG b-waves as compared with a-waves have been reported also in P23H transgenic rats (Machida et al., 2000; Aleman et al., 2001). The reason for increased sensitivity of the ERG b-wave is unknown but may include overcompensation of the gain at the photoreceptor-to-bipolar synapse in response to outer retina degeneration.
The kinetics of photopic ERG b-wave responses showed a trend (not statistically significant) toward hastening in one-month-old RhoP23H/WT mice (Fig. 1). The a-wave response was significantly hastened in RhoP23H/WT mice in response to UV stimulation (Extended Data Fig. 1-2). Pasquale et al., recently reported significantly improved temporal contrast sensitivity (TCS) in RhoP23H/WT mice (Pasquale et al., 2021). The improvement arose from the rod-pathway, as the phenomenon remained when RhoP23H/WT mice were crossed with GNAT2cpfl3/cpfl3 mice that lack cone function. This could be explained by faster photo-response kinetics and decreased density and collecting area of rods in RhoP23H/WT mice (Sakami et al., 2014; Pasquale et al., 2021). The mechanism that explains faster initial rod response is unclear, but could be related to volume changes in the rod outer segment or changes in Ca2+ dynamics (Pasquale et al., 2021). Our recent transcriptomic profiling of RhoP23H/WT retina homogenates supports the latter, as we showed significant pathway enrichments in, e.g., calmodulin binding and ion-channel-activity-related genes, including cation channel activities (Leinonen et al., 2020). However, analysis at single-cell resolution is needed to test whether these changes occur in the rods or elsewhere in the retina.
As no clear deterioration of photopic ERG was observed in RhoP23H/WT mice up to three months of age, and since scotopic contrast sensitivity is remarkably resistant to early rod degeneration (Leinonen et al., 2020), we were curious to see whether any visual acuity change by OMR measures would occur at this disease stage. To this end, we used an OMR system coupled with automated head-tracking (Kretschmer et al., 2015). During the 60 s recording period per SF, the system recorded 1800 frames, and at every frame compared optomotor movements to correct and incorrect directions, rendering an optomotor index. Since the OMR is a reflex-based visual task, we hypothesized that increased tracking toward the drifting patterns (or increased OMR index) by light-adapted RhoP23H/WT mice was possible. This was partially true, as the OMR index at low SF for the RhoP23H/WT mouse was higher compared with that for the WT mouse (Fig. 2). At high SF, the OMR index was equal between genotypes at 1.5 months, indicating intact visual acuity in RhoP23H/WT mice. However, between 1.5 and 3 months of age, the OMR index at high SF slightly declined in RhoP23H/WT mice, indicating an initial decrease in visual acuity by three months in RhoP23H/WT mice. With respect to rod-driven night vision, a strong link between enhanced scotopic ERG b-wave sensitivity and well-maintained OMR in RhoP23H/WT and RhoP23H/WT/Gnat2−/− mice was recently demonstrated (Leinonen et al., 2020).
Consequences of early rod degeneration to higher-order visual processing
The increased ERG b-wave sensitivity and remarkably well-preserved pattern-contrast sensitivity (Leinonen et al., 2020) point to adaptive functional re-arrangements in the RhoP23H/WT mouse retina during rod degeneration. But what happens at the V1 where the “conscious vision” is initially processed? To address this question, we recorded VEPs and single-cell responses in V1 to light ON/light OFF and drifting pattern stimuli in lightly anesthetized mice. First, as recent OMR data have already suggested (Leinonen et al., 2020), pattern contrast sensitivity remained excellent in RhoP23H/WT mice, and strikingly well-maintained even in the homozygous RhoP23H/P23H mice with extreme retinal degeneration (Fig. 6I). A corresponding persistent pattern-contrast sensitivity has been demonstrated by electrophysiology of the dorsal lateral geniculate (dLGN) nucleus in three- to five-week-old Pde6βRd1 mice (Procyk et al., 2019). The Pde6βRd1 mice at that age represent a late-stage retinal degeneration and have only a residual ERG response remaining (Strettoi et al., 2002). Similarly, the temporal-frequency threshold remained excellent in RhoP23H/WT mice, and was in fact elevated at one month of age (Fig. 6G,H). These data correspond well with the ERG (accelerated kinetics, Fig. 1; Extended Data Fig. 1-2; improved TCS by ERGs; Pasquale et al., 2021) and OMR findings in RhoP23H/WT mice (improved TCS by OMR; Pasquale et al., 2021). Nevertheless, we also found detrimental changes for V1 function starting already at one month of age in RhoP23H/WT mice (Figs. 5, 6). The light responses and spontaneous neuron firing were dramatically hyperexcited in the RhoP23H/WT mouse V1 (Figs. 4, 5). The observed hyperexcitability is well in line with previous literature, demonstrating increased neural noise originating from the inner retina (Borowska et al., 2011; Menzler and Zeck, 2011; Trenholm and Awatramani, 2015; Haselier et al., 2017; Telias et al., 2020). Further, as similarly indicated by the OMR behavior (Fig. 2), RhoP23H/WT mouse V1 neurons preferred lower SFs compared with WT mice (Fig. 6E), indicating increasing RGC receptive fields. However, the most striking abnormality we found was the severely debilitated orientation/direction tuning in RhoP23H/WT mice, which advanced to a practical absence of any direction preference by three months of age (Fig. 6A). It is challenging to speculate how these findings translate into clinical settings, because of the lack of corresponding data from RP patients. Hypersensitivity to bright light (called by the debated misnomer “photophobia”) is common in RP (Hamel, 2006; Otsuka et al., 2020), and could in part result from similar hyperexcitation of the primary visual tract as observed in RP model mice. The relatively well-maintained behavioral vision of RhoP23H/WT mice may also be reflected in clinical settings, as many RP patients show minimal disturbances in subjective vision even at an advanced disease stage (Hartong et al., 2006).
Relevance to treatments in retinal degenerations
There are no pharmacotherapies available that would slow retinal degenerative diseases such as RP or dry age-related macular degeneration (Waugh et al., 2018). It is probable, however, that the impaired photoreception in these diseases can be augmented in the near future via novel therapies, such as photoreceptor transplantation (Pearson et al., 2012), chemical photo-switches (Tochitsky et al., 2018), optogenetic tools (Sahel et al., 2021), prosthetics (Lorach et al., 2015), stem cells (da Cruz et al., 2018; Foik et al., 2018), or gene therapy (Suh et al., 2020). But a great deal of concern has been raised whether, after blinding, the more complex inner retina is capable of normal signal processing and whether restored signals are adequate for sight. Second, most patients experience only a partial loss of vision and are not suitable for advanced vision restoration therapies; therefore, more conventional approaches are needed. Using the Pde6βRd10 mouse model of RP, Telias et al. (2019) recently showed that the firing rate of RGCs spontaneously increased simultaneously with decreasing light responses after mouse eye-opening in multi-electrode array recordings. They were able to link the pathology with deteriorated visual behavior in a light avoidance test by decreasing spontaneous firing pharmacologically. Importantly, blocking the retinoic acid receptor decreased spontaneous RGC firing and increased visual behavior in Pde6βRd10 mice (Telias et al., 2019, 2022). Another group used a different pharmacological strategy, a gap junction blockade, but showed similar beneficial effects, wherein a decrease in spontaneous firing of RGCs led to increased light responses (Toychiev et al., 2013). However, the speed of degeneration is faster and the mechanism of rod death is different in the Pde6βRd10 mouse model as compared with the RhoP23H/WT mouse (Chang et al., 2007; Sakami et al., 2011), rendering direct comparisons difficult. It will be important in future experiments to test whether decreasing the noise signal of RGCs would be therapeutic also in other models of RP, including the RhoP23H/WT mouse with hyperexcited V1 (Figs. 3-6).
On the other hand, some recent studies have concluded positive/adaptive aspects of retinal reorganization after injury or during degeneration (Beier et al., 2017, 2018; Johnson et al., 2017; Care et al., 2019; Leinonen et al., 2020; Shen et al., 2020). For example, Care and coworkers ablated 50% of rods in adult mice and uncovered functional compensation at the level of retinal interneurons, which recovered retinal output to near-normal levels (Care et al., 2020). Contextually similar conclusions of homeostatic adaptation to a partial retinal injury were drawn from rod/cone, cone, or rod bipolar cell ablation studies in adult rabbit, ground squirrel, and mouse retinas (Johnson et al., 2017; Beier et al., 2018; Care et al., 2019; Shen et al., 2020). We showed in the RhoP23H/WT mouse, representing a common form of human autosomal dominant RP, increased sensitivity to visual inputs at the photoreceptor-to-bipolar cell synapse and maintained rod-driven contrast sensitivity regardless of rod population loss of >50% (Leinonen et al., 2020). Importantly, functional adaptation to retinal stress has also been observed in glaucoma models, wherein increased excitability of RGCs and dLGN neurons may mitigate the loss of axon function caused by elevated intraocular pressure (Calkins, 2021; Van Hook et al., 2021).
Altogether, accumulating evidence suggests that retinal remodeling on photoreceptor degeneration may not solely possess disruptive elements, but may initially help the visual system by adapting to sensory function loss. As these countering elements seem to temporally overlap, it will be of utmost importance to uncover what are the original molecular triggers for the retinal rewiring and whether these signals divide into more than one downstream pathway. Will it be possible to pharmacologically intervene to promote neural adaptation while simultaneously suppressing the disruptive hyperexcitability? With respect to vision restoration trials in adult humans, neural plasticity will be crucial.
Acknowledgments
Acknowledgments: We thank Dr. Vladimir Kefalov for valuable comments during manuscript preparation. We also thank our colleagues at University of California, Irvine Center for Translational Vision Research and Gavin Herbert Eye Institute for helpful comments regarding this study. Dr. Palczewski is the Donald Bren Professor and the Irving H. Leopold Chair of Ophthalmology.
Footnotes
K.P. is Chief Scientific Officer of Polgenix Inc. All authors declare no competing financial interests.
This work was supported by the National Institute of Health Grant R01EY009339 (to K.P.) and the National Institutes of Health Grant R24EY027283 (to K.P.). H.L. was supported by the Academy of Finland Grant 346295, by a Knights Templar Eye Foundation Career Starter Grant, and by grants from the Finnish Cultural Foundation, the Orion Research Foundation, the Eye and Tissue Bank Foundation (Finland), Retina Registered Association (Finland), and Sokeain Ystävät/De Blindas Vänner Registered Association. This work was also supported by the National Science Center, Poland Grant 2019/34/E/NZ5/00434 (to A.T.F.); The International Centre for Translational Eye Research (MAB/2019/12) project is carried out within the International Research Agendas programme of the Foundation for Polish Science co-financed by the European Union under the European Regional Development Fund (A.T.F.).
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.