The Oscillatory Potentials of the Dark-Adapted Electroretinogram in Retinopathy of Prematurity

Abstract

purpose. To study the development of the electroretinographic (ERG) oscillatory potentials (OPs) in two rat models of ROP and in human subjects with a history of ROP.

methods. Sprague-Dawley rats (n = 36) were studied longitudinally. Rat models of ROP were induced, either by exposure to alternating 50%/10% oxygen (50/10 model) from postnatal day (P) 0 to P14 or by exposure to 75% oxygen (75 model) from P7 to P14. Control rats were reared in room air. Infant and adult human subjects with and without a history of ROP (n = 91) were also studied. Dark-adapted ERGs were recorded and filtered to demonstrate the OPs. Discreet Fourier transform (DFT) allowed evaluation of the OP power spectrum. OP energy (E), dominant frequency (F _peak), and sensitivity (log i _1/2) were evaluated. results. In 50/10 model rats, E was low compared with that in the 75 model rats and control animals. F _peak (∼95 Hz) did not vary with age or group. Intriguingly, log i _1/2 in 75 model rats was greater than that in controls or 50/10 model rats. Human adults with a history of ROP had lower-energy OPs than did the control adults, but infants with a history of ROP had higher-energy OPs than did the control infants. F _peak was lower (∼120 Hz) in infants than in adults (∼130 Hz). ROP did not affect log i _1/2 in humans.

conclusions. Differences between OPs in healthy rats and healthy humans were substantial, suggesting that OPs in rat models of ROP are unlikely to provide insight into the effects of ROP on human OPs. Indeed, neither ROP model studied showed a pattern of effects similar to that in human ROP.

The electroretinogram (ERG) is the only objective, noninvasive method available for determining the function of photoreceptors and other neural cells in the retina. The ERG recorded at the cornea of the dark-adapted vertebrate eye in response to a high-intensity flash includes two principle components: an initial negative deflection, the a-wave, and a subsequent positive deflection, the b-wave. The a-wave originates mainly in the collapse of the circulating current of the photoreceptors; the b-wave originates mainly in the onset of a radial current circulating around depolarizing bipolar cells. ¹ Superimposed on the leading edge of the b-wave is a series of rhythmic, low-amplitude potentials called the oscillatory potentials (OPs). ^{2 3}  

The OPs represent neural activity distinct from that of the a- and b-waves. Early OPs have been associated with photoreceptors and bipolar cells of the distal retina and late OPs with amacrine and ganglion cells of the proximal retina. ^{4 5 6 7 8} However, possibly excepting the first OP which originates in the photoreceptors, ^{9 10} the OPs can be well modeled as one complex. ^{9 10 11} Thus, a neuron that spans the retinal layers and registers the activity of many classes of cells, such as the interplexiform cell, may be important in the generation of the OPs. ^{2 12 13} The amplitude and timing of the OP wavelets can be measured and studied individually or in sum. Alternatively, the OPs can be evaluated in the frequency domain by discreet Fourier transform (DFT) of the record. The DFT technique is especially advantageous when the signal-to-noise ratio is low. 

The OPs are known to be affected in ocular conditions in which abnormal retinal blood vessels are the primary clinical sign, such as diabetic retinopathy, ¹⁴ glaucoma, ^{15 16} vascular occlusion, ¹⁷ and the developmental disorder retinopathy of prematurity (ROP). ¹⁸ Rat models of all these diseases have been reported. ^{19 20 21 22 23 24 25 26 27 28} These ailments have in common oxidative challenge to the inner retina. 

In ERG studies of patients with a history of ROP, significant abnormalities in the a- and b-waves are found to persist into adulthood. ^{18 29 30} Similarly, ERG studies of rat models of ROP, created by exposing infant rats to periods of relatively high and low oxygen during the first weeks after birth, reveal deficits in a- and b-wave parameters similar to those found in the human condition. ^{13 27 28 31 32 33} Because the OPs can measure retinal sensitivity in ways not specified by the a- and b-waves, they may provide additional insight into disease processes in the inner neural retina. Whereas reports indicate that the OPs are diminished in rat ROP, ^{13 33} a comparison of the OPs in ROP rats and in human patients with ROP has not previously been undertaken. In this study, the OPs in two rat models of ROP and in human subjects with a history of ROP were characterized. Major differences in the effects of ROP on the OPs were found between the two species. 

Methods

Rat Studies

Subjects.

Thirty-six Sprague-Dawley albino rats (Charles River Laboratories, Inc., Worcester, MA) from 10 litters were studied. The light cycle was 12 hours dark and 12 hours light (50–100 lux). Care was taken to minimize the effects of mother and birth order. Pups were weaned at postnatal day (P)25. All tests were performed in all subjects at P20 ± 1, P30 ± 1, and P60 ± 1. Twelve rats were exposed to 75% oxygen from P7 to P14 (the 75 model). Twelve other rats were exposed to alternating 24 hour periods of 50% and 10% oxygen from P0 (the day of birth) to P14 (the 50/10 model). The two models produce a range of effects on the retinal vasculature and on the neural retina, as observed in the broad scope of human ROP. ²⁸ The remaining 12 rats were reared in room air as controls. The a- and b-wave responses and the morphology of the retinal vasculature in these rats have been previously described. ²⁸ At P20, the a- and b-waves of the ROP rats (75 and 50/10 models) were markedly attenuated. ^{27 28} All procedures in this study were approved by the Animal Care and Use Committee at Children’s Hospital Boston and were performed in accordance with the ARVO Statement for the Use of Animals in Ophthalmic and Vision Research. 

Procedures.

Electroretinography.

Dark-adapted ERGs were recorded at the cornea of the anesthetized rat by previously reported techniques. ²⁸ The stimuli were brief (<1 ms), white flashes delivered through a 41-cm integrating sphere. The pupil was dilated with cyclopentolate/phenylephrine (Cyclomydril; Alcon, Fort Worth, TX). Stimulus intensity was controlled by calibrated neutral-density filters (Wratten; Eastman-Kodak, Rochester, NY). The unattenuated flash, measured at the position of the rat’s eye using an integrating radiometer (S350; United Detector Technology, Orlando, FL) produced ∼40,000 μW · cm⁻² at the cornea. This flash was calculated to elicit ∼135,000 photoisomerizations of rhodopsin per rod (R*) in the adult control rat. ^{27 34} Flashes of doubling intensity (0.3-log-unit steps) ranging from ∼6.3 R* to ∼65,000 R* were presented. Half of the subjects in each group were not tested at flash intensities between ∼270 R* and ∼2100 R* to abbreviate recording time. Each stimulus was delivered at an interval that did not attenuate subsequent response amplitudes, 2 seconds for the dimmest and more than 2 minutes for the brightest. Nearly all records were single, unaveraged responses. Responses were amplified (×1000; 1–1000 Hz), digitized (2 kHz), and stored for off-line analysis (UTAS-E 2000, LKC Inc., Gaithersburg, MD). 

Measurement of the Oscillatory Potentials.

Figure 1demonstrates the extraction of the OPs from an intact ERG. To mitigate the effects of a-wave contamination on the early OPs, ^{8 9 10 13 38 39} we fitted a mathematical model of the activation of phototransduction, ^{35 36 37} P3, to the leading edge of the a-wave and digitally subtracted the result from each trace. We digitally filtered the difference, P2, a putatively pure postreceptoral response, by using a second-order Butterworth filter ⁴⁰ (MatLab; The MathWorks) with band-pass 60 to 235 Hz. ¹³  

The amplitude (in microvolts) and implicit time (in milliseconds) of the troughs and peaks in the filtered response were determined. For each OP (OP2, -3, -4, and -5), amplitude was defined as the difference between the peak and the trough immediately preceding it. Implicit time was defined as the time from stimulus onset to the peak. Despite the subtraction of P3, OP1 was occasionally visible but was not measured. The interpeak interval, IPI (in milliseconds), was defined as the time between adjacent OP peaks (OP2–3, OP3–4, OP4–5); mean IPI is a measure of periodicity. ³⁹ For each trace, summed OP amplitude (SOPA) and summed OP implicit time (SOPIT) were used to characterize the OP response. 

The OPs were also evaluated in the frequency domain. The first 256 samples (128 ms) of each record were subjected to a fast Fourier transform (radix-2 FFT algorithm; MATLAB; The MathWorks) to produce a power spectrum. Power in the frequency domain increases as the square of response amplitude in the time domain. However, the output of the FFT varies with the number of samples. To calibrate the output of the FFT, sine waves of fixed amplitude and an integer number of cycles per sample-set were subjected to the algorithm. These produced a single spike at the selected frequency that could be calibrated. ⁴¹ It was found that the equation  

\[P(P{^\prime},s){=}\left(2{\cdot}\frac{\sqrt{P{^\prime}}}{\sqrt{s}}\right)^{2}\]

properly calibrated OP power at each frequency, P (μV²), from the radix-2 FFT power spectrum values, P′, and the number of samples used in the FFT, s. Though the resistance in the recording circuit is unspecified, P is linearly proportional to OP power in Watts. 

For responses with components at frequencies that do not result in an integer number of cycles in the measurement window, overspill occurs. ⁴¹ Thus, to characterize the OPs in the frequency domain, the calibrated data were fit by a Gaussian ¹⁰ of the form:  

\[P(F){=}\mathrm{exp}\left[{-}{\frac{1}{2}}{\cdot}\left(\frac{F{-}F_{\mathrm{peak}}}{S}\right)^{2}\right]{\cdot}P_{\mathrm{peak}},\]

where P is the power in the response at frequency F (Hz), F _peak (Hz) is the frequency at the peak of the Gaussian, S (Hz) is a measure of the width of the distribution, and P _peak is the power at F _peak. F _peak is the dominant frequency of the OPs and is related to periodicity. The area under the Gaussian, E (μV² · s⁻¹), was calculated by integrating equation 2 , such that  

\[E{=}P_{\mathrm{peak}}{\cdot}S{\cdot}\sqrt{2{\pi}}\]

where E indicates the energy (proportional to joules) in the OPs and is related to the square of the summed OP amplitude, SOPA. To minimize any residual line-noise (60 Hz), observations in the frequency range of approximately 50 to 70 Hz were excluded from the fits. 

Preliminary inspection of the data indicated that the relationship of OP energy and stimulus intensity was asymptotic. For each rat at each age, E was normalized to its maximum observed value (E _max), and the equation  

\[\frac{E(i)}{E_{\mathrm{max}}}{=}\frac{i^{n}}{i^{n}{+}i_{1/2}^{n}}\]

fit to OP power; both i _1/2 and n were free to vary. In this equation, i is the stimulus intensity (R*) and i _1/2 is the intensity of the stimulus that produces OPs with half-maximum energy. Although 1/i _1/2 is a measure of OP sensitivity, it is analytically convenient to define the sensitivity of the OPs as log i _1/2. 

Statistical Analyses.

OP time domain parameters (SOPA, SOPIT, mean IPI) and frequency domain parameters (E, F _peak, log i _1/2) were evaluated by respective repeated-measures analyses of variance (ANOVA) with the factors age (P20, P30, P60), group (control, 75 model, 50/10 model), and intensity. Two intensities were selected for analysis because not all subjects were tested at all intensities. The lower intensity produced ∼65 R*, too dim to activate cones but sufficient to saturate the rod-driven b-wave; the higher intensity produced ∼8500 R*, sufficient to elicit a saturating a-wave and a saturating rod+cone b-wave. ⁴² Post hoc analyses were performed with the Tukey honestly significant difference (q) statistical test. The Pearson product moment correlation (r ²) was used to evaluate the strength of the relationship between OP parameters and previously reported ²⁸ parameters of the a- and b-waves. The significance level (α) for all tests was P < 0.05. 

Human Studies

Subjects.

Ninety-one subjects participated: 51 infants and 40 adults. All subjects fell into one of four groups: (1) term-born infants in good general health (n = 38; age, 63–78 days at test); (2) ROP infants, infants with a history of ROP (n = 13; age 60–82 days postterm at test); (3) adult controls with good visual acuity and no high refractive errors (n = 26; age, 8–52 years at test); or (4) ROP adults, adults with a history of ROP (n = 14; age, 7–23 years at test). ROP was diagnosed in the neonatal intensive care unit by a pediatric ophthalmologist expert in ROP. In each ROP subject, the ROP was mild and resolved spontaneously. Rod a- and b-wave responses have been reported previously for nearly all subjects in this study. ^{30 43} In normal subjects, all a- and b-wave response parameters are mature by age 6 months. ⁴³ This study was approved by the Committee on Clinical Investigation at Children’s Hospital Boston and adhered to the tenets of the World Medical Association Declaration of Helsinki. Written informed consent was obtained from all subjects or their parents. 

Procedures.

Electroretinography.

Dark-adapted ERGs were recorded from one eye of each subject, using previously described techniques. ³⁹ The pupil was dilated (Cyclogyl 1%; Alcon, Fort Worth, TX). Brief (<1 ms), full-field blue flashes (Wratten 47B; λ < 510 nm) were presented (model 600 VR, series 2100; Novatron, Dallas, TX) in an integrating sphere. Responses were differentially amplified (×1000; 1–1000 Hz), digitized (2.56 kHz), and stored (Compact 4 system; Nicolet Biomedical, Madison, WI). 

Retinal illuminance varies directly with pupil diameter and the transmissivity of the ocular media and inversely with the square of the posterior-nodal distance. ⁴⁴ Though the density of the ocular media is significantly lower in 10-week-old infants than in adults, ⁴⁵ when pupil area, posterior-nodal distance, and ocular media density are taken into account, a given intensity flash produces approximately equal retinal illuminance in young infants and adults. ^{46 47} Based on the work of Kraft et al., ⁴⁸ Hansen and Fulton ⁴⁹ estimated that the difference in rate of photoisomerization in adults and 10 week-old infants was only ∼0.1 log unit. The unattenuated flash produced ∼3.0 log scot td · s, or ∼8500 R*. Neutral density filters were used to present stimuli of doubling intensity over an ∼2.1-log-unit range (∼65 to ∼8500 R*). Stimuli were delivered at a rate that did not attenuate subsequent response amplitudes. 

Measurement of the Oscillatory Potentials.

Similar to the method for the rat ERGs, for each human response, P3 was fitted to the leading edge of the a-wave and digitally subtracted from the record. To demonstrate the OPs, the resultant P2 was digitally filtered with a second-order Butterworth filter with the ISCEV (International Society for Clinical Electrophysiology of Vision)-recommended bandpass of 75 to 300 Hz. ⁵⁰ In the time domain, the amplitudes and implicit times of OP2, -3, -4, and -5 were measured and the OP2–3, OP3–4, and OP4–5 interpeak intervals calculated. The first 128 samples (50 ms) of each record were then subjected to FFT (MatLab; The MathWorks). This duration is sufficient to capture all OPs. ³⁹ The resultant frequency domain data were fitted by a Gaussian, and the OP frequency, F _peak (equation 2 ); energy, E (equation 3 ); and sensitivity, log i _1/2 (equation 4 ) were calculated. 

Statistical Analyses.

OP frequency domain parameters (E, F _peak, log i _1/2) were evaluated by respective repeated measures ANOVAs, with the factors age (infant, adult), clinical history (ROP, no ROP), and intensity (65–8500 R*). Post hoc analyses were performed with the Tukey honestly significant difference (q) statistical test. The significance level (α) for all tests was P < 0.05. 

Results

Rat Studies

For each group (control, 75 model, 50/10 model) and intensity (65 R*, 8500 R*) mean OP2, -3, -4, and -5 trough and peak amplitudes and implicit times were calculated. Figure 2summarizes the time domain analysis in P20 rats. Time-domain OP parameter values at each age and results of ANOVA on these data are presented in Table 1 . Amplitudes were larger and implicit times were shorter for OP responses to the 8500 R* stimulus than to the 65 R* stimulus. At the higher intensity, the OP amplitudes were larger in the control rats than in either the 50/10 model or the 75 model rats; at the lower (scotopic) intensity, the 50/10 model rats’ OP amplitudes were markedly attenuated, but 75 model rats’ OPs differed little from controls’. ROP did not affect SOPIT, but 50/10 model rats’ IPIs were significantly longer than IPIs in either the control or 75 model groups. 

Figure 3shows the OPs of P20 rats transformed into the frequency domain by DFT. In the left panels, mean OP power for control, 75 model, and 50/10 model rats is plotted against frequency and flash intensity. In the right panels, fits of equation 1to the 65 R* and 8500 R* flashes are shown. Frequency domain OP parameter values at each age and results of ANOVAs on these data are presented in Table 2 . OP power at F _peak, P _peak (equation 2 ), and OP energy, E (equation 3 ), the area under the curve, increased with stimulus intensity. E and F _peak were significantly lower in the 50/10 model rats than in the controls. The 75 model rats did not differ significantly from controls in either parameter. 

Figure 4plots the change in the time domain parameters SOPA and IPI, and in the frequency domain parameters E and F _peak, as a percent of control at each test age for the 8500 R* flash. The group × age interaction was not significant for either SOPA (Table 1)or E (Table 2) , indicating that the OPs in ROP rats remained below normal across ages. Periodicity, indicated by IPI and F _peak, became relatively shorter with age in ROP rats. 

OP sensitivity, indicated by the parameter log i _1/2 (equation 4 ; Fig. 5A ), differed significantly from normal in ROP rats (Table 2) . Surprising, however, was the direction of the abnormality. As indicated in Figures 5B and 5C , 75 model rats at all test ages were significantly more sensitive than both controls and 50/10 model rats. At P20, 50/10 model rats were slightly (but not significantly) less sensitive than controls, but at P30 and P60 were slightly more sensitive than controls. Log i _1/2 was the only parameter for which 75 model rats differed significantly from controls (Tukey HSD test). Across groups, at each test age (P20, P30, and P60) both SOPA and E correlated significantly with the previously reported saturating a- and b-wave amplitudes (Rm _P3 and Vm parameters) in these same rats. ²⁸ The sensitivity of the OPs (log i _1/2), however, did not correlate with the sensitivities of either the a- or b-waves (S and log σ parameters). ²⁸ The tortuosity of the retinal vasculature (IC _A parameter) did not correlate significantly with OP amplitude, energy, or sensitivity. ²⁸  

Human Studies

Figure 6Ademonstrates the derivation of the OPs in an adult control subject. The OPs in human adults were comparable in amplitude to those in control rats. However, as previously reported, the amplitudes of the OPs in the healthy, term-born infants were small (median, <5 μV), similar to root mean squared noise (∼3 μV). In the frequency domain, however, a small peak at the frequency expected from the interpeak interval was always found, and thus the DFT proved a superior OP analysis technique to time domain trough and peak measurements in these data with low signal. 

Data from the adult controls, ROP adults, term-born infants, and ROP infants groups, transformed into the frequency domain, are presented in Figure 6B . Frequency domain OP parameters at each age and results of ANOVAs on these data are presented in Table 3 . In the frequency domain data, the spectra in both control and ROP adults were bimodally distributed with peaks at ∼100 and ∼140 Hz at 8500 R*. The higher frequency peak agreed well with interpeak interval (∼7 ms). Accordingly, equation 3was fit to the higher-frequency peak. In contrast, spectra in both term-born and ROP infants show a single peak at ∼127 and ∼112 Hz, respectively. In both adults and infants, with and without a history of ROP, the dominant OP frequency, F _peak, increased with increasing flash intensity (Table 3) . The OPs in ROP adults had less energy than those in adult controls. Surprisingly, OPs in ROP infants had significantly more energy than the OPs in age-matched, term-born infants. 

Figure 7plots the frequency domain parameters in ROP infants and ROP adults as a proportion of the values in term-born infants and adult controls, respectively. The energy in the OPs, E, was significantly higher (>1 log unit) in ROP than in term-born infants but was significantly lower (∼0.3 log units) in ROP than in control adults. Overall, OP energy was significantly higher in adults than in infants. There was no significant difference in the dominant frequency of the OPs, F _peak, between ROP and healthy subjects in either infancy or adulthood. ROP did not affect the sensitivity of the OPs (log i _1/2) in either infants or adults. Sensitivity was, however, significantly greater in adult than in infant subjects. 

Discussion

The OPs were affected by ROP in both species, rat and human, but in different ways. In the rats, the OPs in the 50/10 model were slower, smaller, and lower energy than in controls and did not improve with age. Sensitivity did not differ from controls in 50/10 model rats but was the only parameter that differed from controls in 75 model rats. Intriguingly, the sensitivity of the OPs was greater in 75 model rats than in controls or in 50/10 model rats. In the human subjects, on the other hand, only OP energy was abnormal. The direction of the abnormality changed with age. The energy in the ROP infants’ OPs was dramatically higher, about a log unit on average, than in the term-born infants’, but only half as high in ROP adults as in adult controls. The significantly greater deficits in both the ROP rat models’ OP amplitude (SOPA, Table 1 ) and energy (E, Table 2 ) at higher stimulus intensities may reflect a greater deficit in the cone pathway. In the human subjects, the history × intensity interaction for E (Table 3)was not significant. 

Why, with the exception of OP sensitivity, were the rat OPs less affected by ROP in the 75 model than in the 50/10 model? Presumably, fully differentiated cell types are less susceptible to oxygen perturbations than those still undergoing extensive developmental alterations. ⁵¹ The 50/10 induction begins at P0, when the retina is extremely immature. Ganglion cells are identifiable, but the remainder of the retina is a single layer of neuroblastic cells. ^{52 53} Thus, the 50/10 exposure may interfere with the differentiation and development of every class of neuron in the retina. The 75 model induction begins later, at P7, after all cell types are differentiated and the developmental elongation of photoreceptor outer segments has begun; however, extensive refinements to neural circuitry are still underway. ⁵⁴ Thus, three possible (not mutually exclusive) explanations for the robust OPs in the 75 model are: (1) By P7, the cells that generate the OPs are relatively insensitive to retinal hyperoxia (2); the cells that generate the OPs are more susceptible to hypoxia than hyperoxia; and (3) the alternation between hyperoxia and hypoxia is more damaging to the OP generator(s) than is sustained hyperoxia. In addition, attenuation of inhibitory retinal circuitry results in enhanced OPs. ⁵⁵ There is some evidence that excitatory pathways to retinal ganglion cells in the rat are present, if not mature, by P5, but inhibitory pathways are not. ⁵⁶ Thus, it is also possible that the observed enhanced sensitivity of the OPs in the 75 model is due to selective disruption of later-developing inhibitory synaptic connections. 

Vascular endothelial growth factor (VEGF) is induced by hypoxia and drives developmental vascularization of the retina as well as the pathologic retinal vascularization characteristic of ROP. ^{57 58 59 60} In addition, amacrine cell progenitors are directed to their cell fate by a number of cues, including activation of the VEGF receptor flk-1. ⁶¹ In murine retina, amacrine cells, including the interplexiform cell, are still differentiating at P0 but are fully differentiated at P7. ⁵⁴ Thus, the number and type of postreceptoral cells may be at least partially dependent on the status of retinal oxygenation. We found that the tortuous vasculature characteristic of ROP remodeled (that is, became more normal) remarkably after P20 in the 50/10 model rats, but only slightly in the 75 model rats. ²⁸ Since the OPs are susceptible to disturbances of retinal circulation, ⁶² their function may be related to the severity of retinal vascular abnormality in ROP. However, the data did little to affirm this. There was no significant correlation between vascular abnormality and E or log i _1/2 at P20, P30, or P60. As in the ROP rats, the adult human subjects with a history of ROP had lower energy in their OPs than the adult controls. However, the ROP infants’ OPs contained more energy than those in term-born infants. Developing retinal circuitry must underpin the maturation of the OPs. Both infant groups were age-matched postterm, but the ROP infants had, on average, been using their eyes for more than twice the duration of the term-born infants (∼3 extra months). If use of vision promotes the development of the inner retinal circuitry that generate the OPs, then the ROP infants’ additional early visual experience may have more than offset the negative effects of mild ROP ³⁰ and resulted in the observed increase in their OP energy. Evidence from the sweep visual evoked potential indicated that prematurely born infants, both with ⁶³ and without ⁶⁴ ROP, have better acuity than term-born infants, though behavioral acuity tests at ∼6 months found that the advantage has vanished if not reversed. ⁶⁵ The visual experience of the rats did not vary with group. 

The shapes of the rat and human OP power spectra also differed. First, the dominant frequency of the OPs (F _peak) was much higher in the human subjects than in the rats (∼130 Hz vs. ∼95 Hz at 8500 R*) across all ages and groups. Second, though recorded and processed nearly analogously, the rats’ OP power spectrum was unimodally distributed at all ages, whereas, in agreement with previous reports, ^{9 10} the adult humans’ OP power spectrum was bimodal. The faster peak matched the interpeak interval (∼7 ms). The origin of the lower-frequency peak remains to be specified. This second, slower peak, was not present in the human infants’ power spectra. 

The DFT permitted reliable evaluation of the OPs in the term-born human infants, even though the signal-to-noise ratio in the time domain was very low. Although use of an analog (Butterworth) filter introduces a slight phase shift, ⁴⁰ the filter has the benefit of removing the b-wave ramp-artifact ⁴¹ (low-frequency, high-amplitude noise) from the record. The use of a Gaussian fit to the power spectrum is advised because periodicity and sampling rate both affect OP power, but the area under the curve (E in equation 3 ) remains constant. The peak of the Gaussian (F _peak in equation 2 ) also provides an unbiased measure of the dominant frequency of the OPs that is less susceptible to artifactual spikes in the data. 

That there were marked interspecies differences between the OP power spectra in healthy rats and humans implies that the neural underpinnings of the OPs may be distinct in the two species. Therefore, it is not surprising that the pattern of OP responses in the 75 and 50/10 models of ROP did not resemble the pattern of OP responses in human ROP. Until the origins of the OPs are better specified, it will be difficult to interpret the physiological significance of unexpected observations such as the heightened sensitivity of the OPs in 75 model ROP rats or the increased OP energy in infants with a history of ROP. 

 Supported by National Eye Institute Grant R01 EY10597 and by grants from the Knights Templar Eye Research Fund, the Fight for Sight, the March of Dimes Birth Defects Foundation, the Pearle Vision Foundation, and the Massachusetts Lions Eye Research Fund.

 Submitted for publication July 13, 2007; revised August 20, 2007; accepted October 15, 2007.

 Disclosure: J.D. Akula, None; J.A. Mocko, None; A. Moskowitz, None; R.M. Hansen, None; A.B. Fulton, None

 The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked “advertisement” in accordance with 18 U.S.C. §1734 solely to indicate this fact.

 

Figure 1.

 

Extraction of the electroretinographic oscillatory potentials (OPs). (A) An intact ERG elicited from a P20 control rat by a flash producing ∼8500 photoisomerizations of rhodopsin per rod (R*). Fitted to the leading edge of the a-wave is a mathematical model of the activation of phototransduction, P3 (dashed line). ^{35 36 37} P3 is digitally subtracted from the trace to produce the putatively pure postreceptoral P2. The OPs are demonstrated by passing P2 through a second-order Butterworth filter with band-pass of 60 to 235 Hz. The amplitude of the OPs (numbered) is measured to the peak of the OP wavelet from the trough immediately preceding it. The implicit time of each OP is measured from stimulus onset to the peak. (B) The OPs are transformed into the frequency domain by DFT of the filtered record. A Gaussian (equation 2 ) is fitted to the power spectrum (dashed line).

Figure 2.

 

Time domain analysis. (A) Mean (±SD) OP2, -3, -4, and -5 trough and peak amplitudes and implicit times plotted for P20 control, 75 model, and 50/10 model rats’ responses to the 8500 R* and 65 R* stimuli. In this display, the points are connected with a spline curve (SigmaPlot 10; Systat Software, Inc., San Jose, CA) to simulate average OPs. The mean implicit time for each OP peak is indicated below each record. (B) Mean SOPA, mean SOPIT, and mean IPI (periodicity) are plotted as a function of flash intensity. Dotted lines: intensities shown in (A). Error bars at the top right of each panel indicate the average ± SEM across groups at 8500 R*. Error at other intensities was a similar proportion of the mean.

Table 1.

 

Parameters of Rat OPs in the Time Domain

Parameter	Intensity	Age	Control	75 Model	50/10 Model	ANOVA*
SOPA (μV)	65 R*	P20	207 (125)	209 (158)	62 (65)	F group = 12.15; P < 0.001
		P30	171 (91)	167 (96)	113 (53)	F age = 3.17; P = 0.048
		P60	132 (83)	110 (78)	73 (42)	F intensity = 110.03; P < 0.001
						F group×age = 2.12; P = 0.088
	8,500 R*	P20	425 (164)	283 (158)	121 (74)	F group×intensity = 11.55; P < 0.001
		P30	322 (99)	255 (112)	172 (91)	F age×intensity = 1.04; P = 0.359
		P60	322 (178)	225 (124)	148 (95)	F group×age×intensity = 0.99; P = 0.418
SOPIT (ms)	65 R*	P20	218 (20)	234 (19)	224 (20)	F group = 0.96; P = 0.392
		P30	214 (19)	207 (26)	213 (16)	F age = 22.53; P < 0.001
		P60	199 (18)	192 (16)	199 (25)	F intensity = 679.13; P < 0.001
						F group×age = 1.68; P = 0.166
	8,500 R*	P20	173 (18)	183 (12)	191 (17)	F group×intensity = 1.79; P = 0.184
		P30	166 (11)	164 (9)	177 (13)	F age×intensity = 2.89; P = 0.063
		P60	165 (20)	158 (13)	160 (19)	F group×age×intensity = 2.28; P = 0.070
Mean IPI (ms)	65 R*	P20	11.1 (2.6)	10.9 (1.9)	11.2 (2.2)	F group = 6.25; P = 0.005
		P30	10.4 (2.8)	9.2 (1.5)	9.9 (1.7)	F age = 52.73; P < 0.001
		P60	8.8 (4.7)	8.8 (3.1)	8.6 (2.1)	F intensity = 11.29; P = 0.002
						F group×age = 4.65; P = 0.002
	8,500 R*	P20	11.3 (2.1)	11.8 (2.8)	14.1 (3.4)	F group×intensity = 1.44; P = 0.252
		P30	9.7 (1.4)	9.8 (1.9)	11.5 (2.4)	F age×intensity = 3.99; P = 0.023
		P60	10.3 (2.7)	8.7 (1.5)	9.0 (2.0)	F group×age×intensity = 2.74; P = 0.036

 Data are the mean (SD).

*  Group factor is between subjects; age and intensity factors are within subjects. Significance was defined as P < 0.05.

Figure 3.

 

Frequency domain analysis. Left: the mean OP power as a function of frequency and intensity for P20 control, 75 model, and 50/10 model rats. Bold lines: 8500 R* and 65 R* stimulus levels. Right: replot the frequency domain responses (circles) to these stimuli with fits of equation 2 . The area under the curve (equation 3 ) indicates OP energy.

Table 2.

 

Parameters of Rat OPs in the Frequency Domain

Parameter	Intensity	Age	Control	75 Model	50/10 Model	ANOVA*
E (μV2 · s−1)	65 R*	P20	3,310 (3,450)	3,730 (5,210)	610 (570)	F group = 7.24; P = 0.002
		P30	2,000 (1,510)	1,730 (1,970)	710 (660)	F age = 7.08; P = 0.002
		P60	1,040 (1,270)	620 (1,070)	210 (240)	F intensity = 37.73; P < 0.001
						F group×age = 1.33; P = 0.269
	8,500 R*	P20	9,100 (6,330)	6,320 (7,500)	1,650 (1,430)	F group×intensity = 7.80; P = 0.002
		P30	4,920 (3,080)	2,930 (2,550)	1,690 (1,820)	F age×intensity = 2.84; P = 0.065
		P60	5,240 (7,500)	1,780 (1,930)	840 (1,260)	F group×age×intensity = 0.86; P = 0.490
F peak (Hz)	65 R*	P20	92 (16)	84 (5)	84 (12)	F group = 3.46; P = 0.043
		P30	101 (14)	105 (13)	97 (12)	F age = 62.13; P < 0.001
		P60	107 (16)	113 (16)	111 (17)	F intensity = 4.87; P = 0.034
						F group×age = 4.30; P = 0.004
	8,500 R*	P20	91 (12)	85 (10)	77 (10)	F group×intensity = 1.51; P = 0.236
		P30	98 (9)	101 (15)	88 (11)	F age×intensity = 0.43; P = 0.651
		P60	99 (10)	117 (9)	104 (11)	F group×age×intensity = 0.84; P = 0.507
Log i 1/2 (log R*)		P20	0.705 (0.676)	0.784 (0.574)	0.756 (0.614)	F group = 9.53; P = 0.001
		P30	0.672 (0.713)	0.792 (0.561)	0.710 (0.670)	F age = 3.60; P = 0.033
		P60	0.738 (0.637)	0.788 (0.568)	0.781 (0.579)	F group×age = 0.94; P = 0.448

 Data are the mean (SD).

*  Group factor is between subjects; age and intensity factors are within subjects. Significance was defined as P < 0.05.

Figure 4.

 

The time domain parameters SOPA and SOPIT and frequency domain parameters OP energy (E) and dominant frequency (F _peak) plotted as a percentage of the control (+SEM) at each test age (P20, P30, and P60) for 75 model and 50/10 model rats. The control mean for each parameter is 100% (dotted line).

Figure 5.

 

The sensitivity of the OPs derived from OP energy (E). (A) For each rat at each age, OP energy at every intensity was divided by the maximum observed energy and fit by equation 4 . Shown are representative energy versus intensity relationships for respective P20 control, 75 model, and 50/10 model rats. (B) The sensitivity of the OPs, indicated by the log of the intensity required to elicit a response with half-maximum OP energy, log i _1/2, plotted for each group at P20, P30, and P60. Horizontal bars: mean for each group at each age. (C) Mean OP sensitivity (from B) for 75 model and 50/10 model rats replotted as a proportion of the control mean (+SEM) at each test age. The control mean for each parameter is 0.0 (dotted line).

Figure 6.

 

The OPs in human subjects. (A) Extraction of the time and frequency domain OPs in a healthy adult subject. The steps are as described in Figure 1 . (B) Mean OP power plotted as a function of frequency and intensity for healthy adult controls, adults with a history of ROP, healthy term-born infants, and age-matched infants with a history of ROP. The adult power spectra are bimodal.

Table 3.

 

Parameters of Human OPs in the Frequency Domain

Parameter	Intensity	Group	Mean (SD)	ANOVA*
E (μV2 · s−1)	65 R*	Term-born infants	80 (200)	F history = 0.01; P = 0.935
		ROP infants	370 (620)	F age = 25.92; P < 0.001
		Adult controls	2,840 (2,540)	F intensity = 20.94; P < 0.001
		ROP adults	2,150 (2,280)	F history×age = 11.99; P = 0.001
				F history×intensity = 0.41; P = 0.900
	8,500 R*	Term-born infants	580 (780)	F age×intensity = 4.82; P < 0.001
		ROP infants	7,180 (10,180)	F history×age×intensity = 5.95; P < 0.001
		Adult controls	17,860 (20,930)
		ROP adults	10,020 (8,110)
F peak (Hz)	65 R*	Term-born infants	137 (35)	F history = 0.35; P = 0.558
		ROP infants	146 (32)	F age = 24.81; P < 0.001
		Adult controls	152 (10)	F intensity = 5.69; P < 0.001
		ROP adults	147 (13)	F history×age = 3.44; P = 0.067
				F history×intensity = 0.52; P = 0.817
	8,500 R*	Term-born infants	127 (28)	F age×intensity = 1.03; P = 0.412
		ROP infants	112 (19)	F history×age×intensity = 2.34; P = 0.023
		Adult controls	134 (16)
		ROP adults	140 (14)
Log i 1/2 (log R*)		Term-born infants	3.049 (0.484)	F history = 0.24; P = 0.624
		ROP infants	3.201 (0.276)	F age = 13.34; P < 0.001
		Adult controls	2.753 (0.465)	F history×age = 0.82; P = 0.368
		ROP adults	2.708 (0.546)

*  History (ROP vs. healthy) and age (adult vs. infant) factors are between subjects; intensity factor is within subjects. Significance was defined as P < 0.05.

Figure 7.

 

Left: the frequency domain parameters OP energy (E) and dominant frequency (F _peak) are plotted for ROP subjects as a percent of normal for age. Right: the mean OP sensitivity (log i _1/2) is plotted as a proportion of normal. Error bars, +1 SEM.

The authors thank Professor Michael Bach for aid in understanding the behavior of the DFT and an anonymous reviewer for contributions to the Discussion.