BioMed Research International
Volume 2016 (2016), Article ID 2790194, 12 pages
http://dx.doi.org/10.1155/2016/2790194

Research Article

Assessing the Contribution of the Oscillatory Potentials to the Genesis of the Photopic ERG with the Discrete Wavelet Transform

Mathieu Gauvin,1,2 Allison L. Dorfman,¹ Nataly Trang,¹ Mercedes Gauthier,¹ John M. Little,¹ Jean-Marc Lina,3,4 and Pierre Lachapelle1,2

¹Department of Ophthalmology, Research Institute of the McGill University Health Centre/Montreal Children’s Hospital, Montréal, QC, Canada
²Department of Neurology-Neurosurgery, Research Institute of the McGill University Health Centre/Montreal Children’s Hospital, Montréal, QC, Canada
³Département de Génie Électrique, École de Technologie Supérieure, Montréal, QC, Canada
⁴Centre de Recherches Mathématiques, Montréal, QC, Canada

Received 31 August 2016; Revised 24 October 2016; Accepted 2 November 2016

Academic Editor: Atsushi Mizota

Copyright © 2016 Mathieu Gauvin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The electroretinogram (ERG) is composed of slow (i.e., a-, b-waves) and fast (i.e., oscillatory potentials: OPs) components. OPs have been shown to be preferably affected in some diseases (such as diabetic retinopathy), while the a- and b-waves remain relatively intact. The purpose of this study was to determine the contribution of OPs to the building of the ERG and to examine whether a signal mostly composed of OPs could also exist. DWT analyses were performed on photopic ERGs (flash intensities: −2.23 to 2.64 log cd·s·m⁻² in 21 steps) obtained from normal subjects () and patients () affected with a retinopathy. In controls, the %OP value (i.e., OPs energy/ERG energy) is stimulus- and amplitude-independent (range: 56.6–61.6%; CV = 6.3%). In contrast, the %OPs measured from the ERGs of our patients varied significantly more (range: 35.4%–89.2%; ) depending on the pathology, some presenting with ERGs that are almost solely composed of OPs. In conclusion, patients may present with a wide range of %OP values. Findings herein also support the hypothesis that, in certain conditions, the photopic ERG can be mostly composed of high-frequency components.

1. Introduction

The electroretinogram (ERG) waveform is characterized with a negative a-wave that is generated by the photoreceptors followed by a larger positive b-wave, which originates from the bipolar and Müller cells [1]. Low-voltage high-frequency oscillations, known as the oscillatory potentials (OPs), are also often seen riding on the ascending limb of the b-wave [2, 3]. Although their origin remains debated, it has been suggested that they would represent inner retinal potentials generated by neuronal interactions that might involve bipolar, amacrine, and/or ganglion cells [3–6]. Of interest, while previous studies have shown that, in some retinopathies (such as in diabetic retinopathy or central retinal vein occlusion), the OPs appeared to be selectively more affected compared to the relatively better preserved a- and b-waves [7–9], to date no study has reported ERG responses which seemed to be solely composed of OPs or where the OPs appeared to be selectively better preserved than the slower components of the ERG (i.e., a- and b-waves).

Notwithstanding the above, it must be remembered that traditionally the OPs are extracted using a bandpass filtering technique in order to remove the low-frequency components (i.e., a- and b-waves) from the broadband ERG signal. Unfortunately, bandpass filtering can generate signal distortion such as phase lag, ringing artefacts, and/or attenuation of OP amplitude which can lead to erroneous measures and even, in some instances, “create” artifactual OPs [10, 11]. As a remedy to the latter, it was suggested to quantify the OPs in the frequency domain with the use of the fast Fourier transform (FFT) [6, 12]. Unfortunately, given that the FFT quantifies the power level of all the frequency components contained within a signal (such as the ERG), whether they are time-locked to the stimulus or not (i.e., no temporal resolution), its use can lead to erroneous interpretations [10, 13–15]. Fortunately, the latter limitation can be easily overcome with the use of the discrete wavelet transform (DWT), which is somewhat of an improved FFT since it includes both temporal and frequency resolutions [10, 13–18].

With the above in mind, we sought to determine the contribution of OPs to the building of the photopic ERG waveform obtained from normal subjects and patients and to examine whether a signal mostly composed of OPs could also exist.

2. Methods

2.1. Selection of ERG Responses Analysed

Analysis was performed on photopic ERGs (bandwidth: 1–300 Hz; flash intensities: −2.23 to 2.64 log cd·sec·m⁻² in 21 steps of ~0.2 log-unit; background: 30 cd·m⁻²; averages of 10 to 300 flashes per response) obtained from 40 normal subjects. Results were compared to photopic ERG responses obtained from patients (bandwidth: 1–300 Hz; flash intensities: 0.64 log cd·sec·m⁻²; background: 30 cd·m⁻²; averages of 10 to 300 flashes per response) affected with retinopathies known to selectively abolish the OPs [i.e., diabetic retinopathy (DR) and central retinal vein occlusion (CRVO)] [7–9]. These patients were selected on the basis of the clinical findings that were characteristic of the disease condition (i.e., mostly fundus appearance and in the case of CRVO, unilateral presentation). Results were also compared to photopic ERG response obtained from patients (bandwidth: 1–300 Hz; flash intensities: 0.64 log cd·sec·m⁻²; background: 30 cd·m⁻²; averages of 10 to 300 flashes per response) affected with an advanced retinal degeneration (mostly retinitis pigmentosa) and where, on visual inspection (by a naïve observer: CS), fast oscillations, time-locked to the stimulus and in the frequency range of the OPs (as estimated using a template of a normal OP response) appeared to be the most prominent features of the response. An informed consent form was signed by each subject and the protocol was previously approved by the Institutional Review Board and was conducted in accordance with the declaration of Helsinki. Additional details regarding the ERG setup and recording procedures were previously published by us [19–23].

2.2. Analysis of ERG Responses

The DWT of each ERG signal was computed using the fast wavelet transform algorithm of Mallat [35] implemented with Matlab R2015a. The DWT generates scalograms (Figure 1(a)) which display the energy (-axis) of the signal (maximal values are shown in red; lowest values in blue) as a function of time (-axis) and frequency (-axis). As previously demonstrated by us and others [10, 13, 15–18], this time-frequency approach allows for the identification of energy descriptors, each defined with their respective time and frequency coordinates. The DWT yields measurements of the photopic b-wave (found in the 20 and 40 Hz bands) and OPs (found in the 80 and 160 Hz bands) through the quantification of their respective associated wavelet coefficients [10, 16, 17]. As exemplified in the scalogram of Figure 1(a), two DWT descriptors were used to quantify the 20 Hz and 40 Hz b-wave energy (identified as 20b and 40b) and another two were used to quantify the 80 Hz and 160 Hz OPs energy (identified as 80ops and 160ops), each computed by summating values outlined by the white boxes, respectively. Similar to bandpass filtering, the inverse DWT (IDWT; see [10, 17]) can be used to specifically reconstruct the low-frequency (20 and 40 Hz) and high-frequency (80 and 160 Hz) bands of the signals, which are specific to the slow and fast waves of the ERG. As shown in Figures 1(c) and 1(d), the IDWT confirms that the 20 and 40 Hz descriptors quantify the slow waves (a- and b-waves, as identified in Figure 1(c)) of the ERG, while the 80 Hz and 160 Hz descriptors quantify the fast waves (OPs, identified as 2, 3, and 4 in Figure 1(d)). Consequently, these DWT descriptors were used to quantify the percent contribution of the OPs (%OPs) to the ERG according to the following equation:

Figure 1: Computation of the percentage of oscillatory potentials (%OPs) contributing to the ERG response. The DWT scalogram (a) is shown above the broadband ERG (b). Temporal sections of the scalogram that were taken into consideration when calculating the %OPs are outlined by the dotted gray lines in order to ease appreciation of the time-frequency region in question. The DWT descriptors (20b, 40b, 80ops, and 160ops) that were used to compute the %OPs are identified on the scalogram. The inverse DWT (IDWT) reconstructions of the slow (a- and b-waves, indicated by the corresponding letters) and fast (OPs, indicated by corresponding numbers) waves are shown as the black and red traces (panels (c) and (d)), respectively.

2.3. Statistical Analysis

Descriptive statistics (mean and standard deviation [SD]) of the %OPs were obtained from our selected patients and from our normal subject cohort at each of the 21 stimulus intensities. The coefficient of variation (CV) of selected groups was computed as the SD divided by the mean and multiplied by 100. To further confirm that each of our selected pathological ERGs had a significantly () lower (DR and CRVO patients) or higher (end-stage retinal degeneration patients) than normal %OPs, individual one-way -tests were used with a critical value set at 1.645 [36]. Thus -scores of 1.645 will indicate a significantly higher/lower () than normal contribution of the OPs to the broadband ERG, respectively. Finally, an unpaired -test (one-tail) was also used to validate whether our group of advanced retinal degeneration patients had a significantly () higher percentage of OPs (%OPs) contributing to their ERGs.

3. Results

3.1. %OPs Contribution to the Normal ERG Response

Representative ERG waveforms of the photopic luminance-response (LR) function obtained at four different stimulus intensities are shown in Figures 2(a) to 2(d) and the positions of these responses on the photopic LR function are indicated with a red diamond in Figure 2(e). As shown in Figure 2(a), no OPs could be detected on the rising phase of the b-wave evoked to the dimmest stimulus (i.e., −1.41 log cd·m⁻²) as well as in the response (Figure 2(b)) evoked to a stimulus one log-unit brighter (−0.41 log cd·m⁻²), despite a 10-fold increase in b-wave amplitude (i.e., from ~3 to ~30 μV). The latter contrasts with the ERG evoked to the 0.64 log cd·s·m⁻² stimulus (Figure 2(c)) where two prominent OPs are seen on the ascending phase of the b-wave. A further increase in stimulus intensity to 2.39 log cd·s·m⁻² (Figure 2(d)) will yield a b-wave of reduced amplitude (approximately 30 μV) compared to the 0.64 log cd·s·m⁻² response but with more OPs. The above suggests, at least from visual inspection, that the prominence of OPs increases with stimulus intensity, but not necessarily with b-wave amplitude. However, as indicated at the top of the four representative scalograms (Figures 2(a)–2(d)), when we measure the %OPs content of these ERGs, nearly identical values are found (%OPs = 58.7; 57.8; 56.5; and 59.4 for intensities of −1.41; −0.41; 0.64; and 2.39 log cd·s·m⁻², resp.). This is best exemplified with the graphs of Figures 2(e) and 2(f) reporting the luminance-response functions of DWT descriptors of the b-wave (i.e., 20b + 40b; gray trace; diamond markers) and OPs (i.e., 80ops + 160ops; black traced; round markers) are nearly identical (once normalized). Therefore, as reported in Table 1 and illustrated in Figure 2(f), in normal subjects, across all of the 21 stimulus intensities (close to 5 log-unit range), the %OPs varied between 56.6 (smallest value) and 61.6 (highest value) and was thus significantly less variable than the b-wave amplitude which varied between 1.2 μV (lowest value) and 115.52 μV (highest value). The above demonstrate not only that the %OPs is nearly stimulus-independent but also that it varies very little (see Table 1) as the mean %OPs coefficient of variation (CV) of all the intensities considered was of 6.34% compared to 87.37% for the b-wave amplitude.

Table 1: Normal variation of the %OPs as a function of the 21 stimulus intensities. Column 1: stimulus intensities. Column 2: mean %OPs values obtained at each stimulus intensity. Column 3: coefficient of variation (CV) of the %OPs at each stimulus intensity. Column 4: mean b-wave amplitude obtained at each stimulus intensity. Column 5: CV of the b-wave amplitude at each stimulus intensity. Range, mean (±1 SD).

Intensity (log cd·s·m⁻²)
Mean %OPs (%)
%OPs CV (%)
b-wave amplitude (μV)
b-wave amplitude CV (%)

60.27
12.64
1.23
33.74

61.56
9.96
1.54
30.73

61.27
8.41
1.93
26.79

60.75
7.56
2.67
20.10

58.19
7.15
3.70
20.28

59.23
8.35
5.02
18.27

60.08
9.13
10.86
17.60

61.41
9.40
18.24
27.41

60.47
6.66
28.61
21.25

58.67
7.56
39.02
21.39

56.60
6.12
51.26
23.68

57.76
5.13
76.81
25.69

0.17
59.02
4.56
102.48
22.73

0.39
58.86
5.51
115.42
22.04

0.64
56.64
4.42
112.78
17.25

0.90
59.96
3.26
100.49
18.36

1.40
59.89
3.68
89.71
20.87

1.63
61.38
5.98
57.92
18.55

2.15
60.19
6.61
47.64
20.92

2.39
59.56
5.23
46.35
17.96

2.64
56.99
6.35
49.35
24.36

Range
56.6–61.6
3.3–12.6
1.2–115.4
17.3–33.7

Mean ± SD

Total CV
6.34

87.37

Figure 2: Normal ERGs (obtained at four different stimulus intensities) which show a relatively constant %OPs contribution (from (a) to (d): 58.7%; 57.8%; 56.5%; 59.4%). The four examples are shown in increasing order of stimulus intensities from panels (a) to (d). The DWT scalograms are shown at the top of each broadband ERG. Sections of the scalogram that were not considered in the %OPs computation were shaded to facilitate comparison. The descriptors (20b, 40b, 80ops, and 160ops) that were used to compute the %OPs are identified on each scalogram. (e) Normalized luminance-response functions of the DWT descriptors of the b-wave (20b + 40b) and OPs (80ops + 160ops). The red markers refer to the intensities shown from panels (a) to (d). (f) Mean ± SD %OPs values obtained at each stimulus intensity.

3.2. Pathological ERG Responses Presenting with a Reduced %OPs Contribution

As predicted from previous studies [7–9], ERGs recorded from patients affected with diabetic retinopathy (DR, Figure 3(b)) or central retinal vein occlusion (CRVO, Figure 3(c)) presented no evidence of OPs on the rising phase of the b-wave. These conditions allowed us to investigate if our method of determining the OP content with the DWT (i.e., the %OPs) would be reduced in situations where the OPs are selectively attenuated. The latter is clearly reflected with the significantly reduced %OPs of 35.4% (-score: −8.48; ) and 36.7% (-score: −7.96; ) that we found for DR and CRVO, respectively.

(a)

(b)

(c)

Figure 3: (a) %OP contribution to the photopic ERG response obtained from a normal subject and patients with diabetic retinopathy (DR) (b) and CRVO (c). All photopic ERG responses were evoked to the 0.64 log cd·s·m⁻² stimulus. Sections of the scalogram that were not considered in the %OPs computation were shaded to facilitate comparison. The descriptors (20b, 40b, 80ops, and 160ops) that were used to compute the %OPs are localized on each scalogram. The %OPs is shown at the top of each DWT scalogram.

3.3. Pathological ERG Responses Presenting with a Higher %OPs Contribution

Analysis of ERGs obtained from our databank revealed that the patients in this group were mostly affected with retinitis pigmentosa (RP), accounting for a total of 15 cases. The clinical findings, such as visual fields, visual acuities, and rod ERG amplitudes, which are reported in Table 2, are indicative of end-stage retinopathies. Using the above-defined DWT technique, we determined that the overall percentage of OPs contribution was significantly greater than that obtained from control (%OPs of versus , resp., ) using the same suprathreshold photopic ERG waveform (0.64 log cd·s·m⁻²). Representative examples of enhanced %OPs ERGs are shown in Figure 4 and data from patients are reported in Table 2. For example, as shown in Figure 4(a), patient 1 (also referred to as patient 1 in Table 2) presented with a %OPs of 61.1%, a value that is slightly but nonetheless significantly higher than control (-score: 1.78; ). In contrast, the ERGs of patients 2 and 3 (Figures 4(b) and 4(c)) presented with a more pronounced enhancement of the %OPs parameter to 75.4% and 78.4%, respectively (-scores: 7.46 and 8.67; ). Finally, the ERG tracing and scalogram obtained from patient 4 (Figure 4(d)) disclosed a waveform that is almost solely composed of OPs with a %OPs of 89.2% (-score: 13.14; ), the highest %OPs value obtained from our patient cohort.

Table 2: Ophthalmological and electrophysiological findings of the end-stage patient cohort included in this study. Patient number, age, visual acuity, visual field findings, rod- and cone-mediated ERG responses, diagnosis, %OPs values, and associated -scores are each indicated in individual columns. The %OPs contributing to the cone ERG response and associated -scores are shown in the last two columns, respectively. The critical -score was set to a value of 1.645 and therefore any value > 56.6 + (i.e., 60.7) is considered to be significant. LP: light perception; flat: extinguished ERG; RP: retinitis pigmentosa.

Patient
Age
Visual acuity
Visual field
Rod ERG
Cone ERG
Diagnosis
%OPs
-scores

1
62
20/50
<20°
Flat
<10%
RP
61.1
1.78

2
61
20/400
<10°
Flat
<5%
RP
75.4
7.46

3
51
20/200

<2%
<10%
RP
78.4
8.67

4
9
20/20
<15°
Flat
<10%
RP
89.2
13.14

5
44
20/100
<10°
Flat
<10%
RP
82.0
10.30

6
35
LP
<10°
Flat
<2%
RP
86.5
11.92

7
75
20/400
<10°
Flat
<2%
RP
74.9
7.46

8
6
20/30
N/A
Flat
<5%
RP
73.4
6.65

9
35
20/20
N/A
Flat
<5%
RP
67.1
4.21

10
73
20/400
<10°
Flat
<5%
RP
78.7
9.08

11
39
20/20
<15°
<1%
<10%
RP
64.3
3.00

12
51
20/50
<15°
Flat
<10%
RP
68.6
5.02

13
44
20/20

Flat
<10%
RP
65.4
3.40

14
58
20/60
<10°
Flat
<5%
RP
78.5
8.67

15
53
20/100
<10°
Flat
<10%
RP
71.8
6.24

16
24
20/25
<
<5%
<15%
Choroideremia
65.2
3.40

17
19
20/30
<
<5%
<10%
Choroideremia
79.9
9.49

of an annular scotoma

temporal scotoma with isopter IVe and <10% with isopter IIe

IIe.

(a)

(b)

(c)

(d)

Figure 4: Pathological ERGs (obtained from four different end-stage RP patients) which show a significantly higher %OPs contribution. The four cases are shown in increasing order of %OPs from panels (a) to (d). In each case, the DWT scalogram is shown on the top of the broadband ERGs. Temporal sections of the scalogram that were not considered in the %OPs computation were shaded to ease appreciation for the region in question. The descriptors (20b, 40b, 80ops, and 160ops) that were used to compute the %OPs are shown on each scalogram. The %OPs is shown at the top of each DWT scalogram.

3.4. Example of %OPs Progression with Disease Process

Figure 5(a) shows the photopic ERGs obtained from two brothers (aged 12 and 17 years at initial exam and identified as patients 16 and 17, resp., in Table 2) affected with choroideremia. On initial examination, both brothers presented with severely attenuated ERGs, albeit with normal looking morphologies, that contrast with the more oscillatory ERG waveforms that we obtained seven years later. The corresponding DWT scalograms indicate that while, at the initial visit, the %OPs was within the normal limits (i.e., 56.2% and 55.1%), seven years later it had increased to 65.2% and 79.9%, suggesting that disease progression was most detrimental to the slow frequency generators of the ERG (e.g., b-wave).

(a)

(b)

Figure 5: (a) A control ERG (black tracing) is shown with the ERG follow-up of patients #16 (blue tracings) and 17 (red tracings), at first visit (12 and 17 years of age, resp.) and at second visit (24 and 19 years of age, resp.). The a-, b-, and i-waves are indicated by their corresponding letters and OP2, OP3, and OP4 by their corresponding numbers. Calibration: horizontal: 50 ms; vertical: as shown to the right of each tracing. (b) Corresponding DWT scalograms and associated %OPs. Descriptors (20b, 40b, 80ops, and 160ops) that were used to compute the %OPs are localized on each scalogram. Note that, in the scalogram corresponding to the 2nd visit of Brother 2 (bottom tracing), the 40b descriptor was identified using the value of the left rectangle (yellow) as the higher energy red rectangle was time-locked to the i-wave.

4. Discussion

The purpose of this study was to determine the energy contribution of the oscillatory potentials to the building of the photopic ERG response and to investigate if an ERG could be mostly composed of OPs. Our study demonstrated that, in normal ERGs, the %OP value is relatively constant (overall CV = 6.3%), not stimulus-dependent (range 56.6% to 61.6%), and consequently not influenced by the absolute amplitude of the ERG b-wave (which varied between 1.2 and 115.4 μV; overall CV = 87.4%). This nearly stimulus-independent %OPs value can be explained by the fact that the b-wave and OPs covary across the luminance-response function of the normal photopic hill, increasing in the initial portion, reaching a plateau, and then decreasing with brighter intensities. The latter is supported by previous studies which showed that OPs and b-wave parameters as well as their luminance-response (LR) functions were highly covariant [17, 37–41]. For example, a correlation coefficient of 0.78 () was previously reported between the b-wave amplitude and the sum of OPs (SOP) amplitude [39] and an even higher correlation coefficient of 0.98 () was previously found between LR functions of the b-wave and OPs energy [17]. However, the nearly stimulus-independent %OPs value can be surprising to some given that, on visual inspection, waveforms evoked to progressively brighter stimuli appear to be gradually more oscillatory as exemplified in Figure 2. Nevertheless, DWT analysis suggests that more prominent and more numerous OPs on the ascending limb of the b-wave do not necessarily indicate a larger contribution of the OPs to the building of the ERG since the %OPs value was shown to be nearly stimulus-independent.

The above range of normal %OPs (56.6% to 61.6%) also demonstrates that, in normal subjects, the summed OPs energy contributes to more than half (%OPs > 50%) of the total energy contained within the ERG waveform. Based on normal data included in previously published studies where OPs were extracted using bandpass filtering, we derived (see computation details in Table 3 caption) time domain %OPs (TD%OPs) values (reported in Table 3) as low as 26% [24] or as high as 112% [31]. As aforementioned, bandpass filtering can generate signal distortions such as phase lag, ringing artefacts, attenuation of OPs amplitude, and/or artefactual enhancement of OPs. The latter most probably accounts for the very large variation of the time domain (TD) derivation of the %OPs reported in Table 3. Intuitively, the “real” %OPs value should lie somewhere between the two extremes (i.e., close to the mean of all values). Of interest, the average normal TD%OPs derived from all previous studies (Table 3) was of , a value that is not significantly different () from the %OPs of the normal photopic ERG waveform (irrespective of stimulus intensity) that we quantified with the DWT ( as per Table 1). The latter would suggest that the %OPs computed with the DWT represents an accurate and highly reproducible (CV as low as 3.26% for the 0.9 log cd·s·m⁻² in Table 1) method to estimate the %OPs contribution to the building of the ERG.

Table 3: Time domain (TD) derivation of the TD%OPs contribution to the photopic ERG from previous studies (column 1). Column 2: TD%OPs values are presented in increasing order and were computed as the sum of OPs (SOP) amplitude divided by the associated b-wave amplitude and multiplied by 100. Column 3: the type of filters that were used. N/A indicates unspecified filter type. Columns 4 and 5: ERG and OPs cutoff frequencies that were used for bandpass filtering. Column 6: the attenuation that was used at the specified cutoff frequencies. N/A indicates unspecified attenuation.

Papers
TD%OPs (%)
Filter type
ERG cutoff (Hz)
OPs cutoff (Hz)
Attenuation (dB)

Miyake, 1990 [24]
26
N/A
5.3–100
100–300
N/A

Lachapelle et al., 1993 [25]
28
N/A
1–1000
100–1000
6

Peachey et al., 1991 [26]
36
N/A
1–1000
100–1000
3

Kizawa et al., 2006 [7]
38
N/A
0.5–1000
100–1000
N/A

Lachapelle et al., 1998 [27]
44
N/A
1–1000
100–1000
6

Holopigian et al., 1992 [28]
48
N/A
1–300
100–1000
N/A

Rufiange et al., 2005 [22]
69
N/A
0.3–500
75–500
6

Fortune et al., 2004 [29]
75
Blackman filter
0.3–3000
70–280
3

Kergoat et al., 2001 [30]
99
N/A
1–1500
50–1500
N/A

Peachey et al., 1987 [31]
112
N/A
1–1000
100–1000
3

Contrasting with the above findings obtained from normal subjects, data obtained from our patients did not present a constant %OPs value and showed a wide range of %OPs, some presenting with a relatively suppressed %OPs (such as 35.4%, CRVO ERG) and others with a relatively enhanced %OPs (such as 89.2%, end-stage RP ERG). Our analysis revealed that pathological ERGs presenting with OPs more prominent than normal (on visual inspection) disclosed a higher than normal %OPs, with some almost solely composed of OPs. However, ERGs from RP patients (see Figure 4) were sometimes of very low amplitude (and some with poor signal-to-noise ratio (SNR), e.g., Figure 4(b)) and, as such, one wonders if the %OPs value was repeatable within subjects. To investigate the latter, a representative patient was recorded twice to test for reproducibility (Figure 6). Despite the very low SNRs of 1.3 and 1.6 obtained on test and retest, respectively, the response morphology was found to be reproducible and so was the resulting %OPs (test: %OPs = 70.1%; retest: %OPs = 69.1%). The latter suggest that a repeatable measurement of the %OPs can be obtained, even in low-SNR ERGs. To our knowledge, the findings reported herein demonstrate, for the first time, that in certain conditions the photopic ERG can be mostly (or solely) composed of high-frequency components, known as the OPs.

Figure 6: Representative test-retest repeatability of the %OPs. Tracings (top and bottom panels) were recorded twice in the same RP patient. The %OPs is shown at the top of each tracing.

Interestingly, studies reporting on (slow sequence) multifocal ERGs of normal macaque monkeys revealed that responses obtained from the central retina were more oscillatory in nature and characterized with highly prominent high-frequency waves [6]. Similarly, focal macular ERG responses evoked from patients affected with RP also presented OPs that were relatively better preserved compared to the a- and b-waves, thus suggesting that macular responses are more oscillatory in nature [42], a finding also confirmed with focal macular ERGs obtained in a rabbit model of RP [43]. Given that most of our patients (as summarized in Table 2) that presented with highly prominent OPs had a very constricted visual field (which is a common finding in end stage of rod/cone dystrophies such as RP), the more oscillatory ERG responses that we found could find an explanation in the relatively better preserved central (macular) function. The above would therefore suggest that the macular region would contain a greater proportion of inner retinal cells which are suggested to be involved in OPs generation (i.e., amacrine, bipolar, horizontal, or ganglion cells) compared to retinal cells involved in the genesis of the b-wave (i.e., bipolar and Müller cells). This claim would find support from previously published findings on retinal cell distribution in nonhuman primates [32–34]. Based on these studies, we computed (Table 4) a significantly greater bipolar-to-Müller cell ratio (BMR) at the fovea compared to more peripheral regions (i.e., 24, 40, and 65 degrees from fovea), thus supporting our claim.

Table 4: Derivation of the bipolar-to-Müller cell ratio at different retinal eccentricities from previous nonprimate studies.

Eccentricity
Cones density
Cone bipolar cells density
Müller cells density
Bipolar-to- Müller cell ratio

Fovea
220,000
244444
30,000
8.15

24 deg.
35,000
4730
25,000
0.19

40 deg.
22,000
4889
20,000
0.24

65 deg.
17,000
1809
15,000
0.12

from Wilder et al. [32].

from Chan et al. [33].

from Distler and Dreher [34].

In conclusion, normal subjects disclosed a relatively constant and highly reproducible %OPs values while patients presented with a significantly larger range of %OPs values. We postulated that, in some circumstances, a pathological ERG signal could be mostly (or solely) composed of OPs. Our analyses did indeed reveal that some ERGs were almost solely composed of OPs. Findings herein therefore support the hypothesis that, in certain conditions, the photopic ERG can be mostly composed of high-frequency components, known as the OPs. Furthermore, the significantly wider range of %OPs values measured in patients (with variable b-wave amplitudes) would suggest that one of the two ERG components (b-wave or OPs) might be relatively more affected than the other and the DWT therefore represents a worthwhile approach to help in the segregation of pathological ERGs.

Competing Interests

The authors declare that they have no competing interests.

Authors’ Contributions

Mathieu Gauvin and Allison L. Dorfman equally contributed to this study and should therefore be considered as equal-first authors.

Acknowledgments

This work was supported by grant-in-aid from the Canadian Institutes for Health Research (CIHR MOP-126082), the CIHR (ERA-132932), and the Fonds de Recherche du Québec-Santé (FRQ-S JTC 2013), under the frame of E-Rare-2, the ERA-Net for Research on Rare Diseases, and a doctoral award from the FRQ-S and its thematic research network (Vision Network). The authors are grateful to Cléa Simard (CS) for her help with the chart review process.

References

R. F. Miller and J. E. Dowling, “Intracellular responses of the Müller (glial) cells of mudpuppy retina: their relation to b-wave of the electroretinogram,” Journal of Neurophysiology, vol. 33, no. 3, pp. 323–341, 1970.

L. Wachtmeister, “Basic research and clinical aspects of the oscillatory potentials of the electroretinogram,” Documenta Ophthalmologica, vol. 66, no. 3, pp. 187–194, 1987.

L. Wachtmeister, “Oscillatory potentials in the retina: what do they reveal,” Progress in Retinal and Eye Research, vol. 17, no. 4, pp. 485–521, 1998.

H. Heynen, L. Wachtmeister, and D. van Norren, “Origin of the oscillatory potentials in the primate retina,” Vision Research, vol. 25, no. 10, pp. 1365–1373, 1985.

L. Wachtmeister and J. E. Dowling, “The oscillatory potentials of the mudpuppy retina,” Investigative Ophthalmology and Visual Science, vol. 17, no. 12, pp. 1176–1188, 1978.

W. Zhou, N. Rangaswamy, P. Ktonas, and L. J. Frishman, “Oscillatory potentials of the slow-sequence multifocal ERG in primates extracted using the Matching Pursuit method,” Vision Research, vol. 47, no. 15, pp. 2021–2036, 2007.

J. Kizawa, S. Machida, T. Kobayashi, Y. Gotoh, and D. Kurosaka, “Changes of oscillatory potentials and photopic negative response in patients with early diabetic retinopathy,” Japanese Journal of Ophthalmology, vol. 50, no. 4, pp. 367–373, 2006.

K. Takei, T. Sato, T. Nonoyama, T. Miyauchi, K. Goto, and S. Hommura, “A new model of transient complete obstruction of retinal vessels induced by endothelin-1 injection into the posterior vitreous body in rabbits,” Graefe's Archive for Clinical and Experimental Ophthalmology, vol. 231, no. 8, pp. 476–481, 1993.

D. Yonemura, T. Aoki, and K. Tsuzuki, “Electroretinogram in diabetic retinopathy,” Archives of ophthalmology, vol. 68, pp. 19–24, 1962.

M. Gauvin, J.-M. Lina, and P. Lachapelle, “Advance in ERG analysis: From peak time and amplitude to frequency, power, and energy,” BioMed Research International, vol. 2014, Article ID 246096, 11 pages, 2014.

M. F. Marmor, A. B. Fulton, G. E. Holder, Y. Miyake, M. Brigell, and M. Bach, “ISCEV Standard for full-field clinical electroretinography (2008 update),” Documenta Ophthalmologica, vol. 118, no. 1, pp. 69–77, 2009.

M. Gur and Y. Zeevi, “Frequency-domain analysis of the human electroretinogram,” Journal of the Optical Society of America, vol. 70, no. 1, pp. 53–59, 1980.

H. Kundra, Measurement of the photopic negative response of the electroretinogram by discrete wavelet analysis [Ph.D. thesis], University of Illinois at Chicago, Chicago, Ill, USA, 2015.

S. S. Nair and K. P. Joseph, “Chaotic analysis of the electroretinographic signal for diagnosis,” BioMed Research International, vol. 2014, Article ID 503920, 8 pages, 2014.

S. S. Nair and K. Paul Joseph, “Wavelet based electroretinographic signal analysis for diagnosis,” Biomedical Signal Processing and Control, vol. 9, no. 1, pp. 37–44, 2014.

M. Gauvin, H. Chakor, R. K. Koenekoop, J. M. Little, J.-M. Lina, and P. Lachapelle, “Witnessing the first sign of retinitis pigmentosa onset in the allegedly normal eye of a case of unilateral RP: a 30-year follow-up,” Documenta Ophthalmologica, vol. 132, no. 3, pp. 213–229, 2016.

M. Gauvin, J. M. Little, J.-M. Lina, and P. Lachapelle, “Functional decomposition of the human ERG based on the discrete wavelet transform,” Journal of Vision, vol. 15, no. 16, article 14, 2015.

H. Kundra, J. C. Park, and J. J. McAnany, “Comparison of photopic negative response measurements in the time and time-frequency domains,” Documenta Ophthalmologica. Advances in Ophthalmology, vol. 133, no. 2, pp. 91–98, 2016.

M.-L. Garon, A. L. Dorfman, J. Racine, R. K. Koenekoop, J. M. Little, and P. Lachapelle, “Estimating on and off contributions to the photopic hill: normative data and clinical applications,” Documenta Ophthalmologica, vol. 129, no. 1, pp. 9–16, 2014.

M.-L. Garon, M. Rufiange, R. Hamilton, D. L. McCulloch, and P. Lachapelle, “Asymmetrical growth of the photopic hill during the light adaptation effect,” Documenta Ophthalmologica, vol. 121, no. 3, pp. 177–187, 2010.

M. Rufiange, J. Dassa, O. Dembinska et al., “The photopic ERG luminance-response function (photopic hill): method of analysis and clinical application,” Vision Research, vol. 43, no. 12, pp. 1405–1412, 2003.

M. Rufiange, M. Dumont, and P. Lachapelle, “Modulation of the human photopic ERG luminance-response function with the use of chromatic stimuli,” Vision Research, vol. 45, no. 17, pp. 2321–2330, 2005.

M. Rufiange, S. Rousseau, O. Dembinska, and P. Lachapelle, “Cone-dominated ERG luminance-response function: the Photopic Hill revisited,” Documenta Ophthalmologica, vol. 104, no. 3, pp. 231–248, 2002.

Y. Miyake, “Macular oscillatory potentials in humans. Macular OPs,” Documenta Ophthalmologica, vol. 75, no. 2, pp. 111–124, 1990.

P. Lachapelle, J. Benoit, J. M. Little, and B. Lachapelle, “Recording the oscillatory potentials of the electroretinogram with the DTL electrode,” Documenta Ophthalmologica, vol. 83, no. 2, pp. 119–130, 1993.

N. S. Peachey, K. R. Alexander, D. J. Derlacki, P. Bobak, and G. A. Fishman, “Effects of light adaptation on the response characteristics of human oscillatory potentials,” Electroencephalography and Clinical Neurophysiology, vol. 78, no. 1, pp. 27–34, 1991.

P. Lachapelle, S. Rousseau, M. Mckerral et al., “Evidence supportive of a functional discrimination between photopic oscillatory potentials as revealed with cone and rod mediated retinopathies,” Documenta Ophthalmologica, vol. 95, no. 1, pp. 35–54, 1998.

K. Holopigian, W. Seiple, M. Lorenzo, and R. Carr, “A comparison of photopic and scotopic electroretinographic changes in early diabetic retinopathy,” Investigative Ophthalmology and Visual Science, vol. 33, no. 10, pp. 2773–2780, 1992.

B. Fortune, B. V. Bui, G. Cull, L. Wang, and G. A. Cioffi, “Inter-ocular and inter-session reliability of the electroretinogram photopic negative response (PhNR) in non-human primates,” Experimental Eye Research, vol. 78, no. 1, pp. 83–93, 2004.

H. Kergoat, M.-J. Kergoat, and L. Justino, “Age-related changes in the flash electroretinogram and oscillatory potentials in individuals age 75 and older,” Journal of the American Geriatrics Society, vol. 49, no. 9, pp. 1212–1217, 2001.

N. S. Peachey, H. C. Charles, C. M. Lee, G. A. Fishman, J. G. Cunha-Vaz, and R. T. Smith, “Electroretinographic findings in sickle cell retinopathy,” Archives of Ophthalmology, vol. 105, no. 7, pp. 934–938, 1987.

H. D. Wilder, U. Grunert, B. B. Lee, and P. R. Martin, “Topography of ganglion cells and photoreceptors in the retina of a new world monkey: the marmoset Callithrix jacchus,” Visual Neuroscience, vol. 13, no. 2, pp. 335–352, 1996.

T. L. Chan, P. R. Martin, N. Clunas, and U. Grnert, “Bipolar cell diversity in the primate retina: morphologic and immunocytochemical analysis of a new world monkey, the marmoset Callithrix jacchus,” Journal of Comparative Neurology, vol. 437, no. 2, pp. 219–239, 2001.

C. Distler and Z. Dreher, “Glia cells of the monkey retina—II. Muller cells,” Vision Research, vol. 36, no. 16, pp. 2381–2394, 1996.

S. G. Mallat, A Wavelet Tour of Signal Processing: The Sparse Way, Elsevier/Academic Press, Boston, Mass, USA, 3rd edition, 2009.

R. C. Sprinthall, Basic Statistical Analysis, Pearson Allyn & Bacon, Boston, Mass, USA, 9th edition, 2012.

P. Guité and P. Lachapelle, “The effect of 2-amino-4-phosphonobutyric acid on the oscillatory potentials of the electroretinogram,” Documenta Ophthalmologica, vol. 75, no. 2, pp. 125–133, 1990.

P. Lachapelle, “Analysis of the photopic electroretinogram recorded before and after dark adaptation,” Canadian Journal of Ophthalmology, vol. 22, no. 7, pp. 354–361, 1988.

P. Lachapelle, “Oscillatory potentials as predictors to amplitude and peak time of the photopic b-wave of the human electroretinogram,” Documenta Ophthalmologica, vol. 75, no. 1, pp. 73–82, 1990.

P. Lachapelle and J. Benoit, “Interpretation of the filtered 100- to 1000-Hz electroretinogram,” Documenta Ophthalmologica, vol. 86, no. 1, pp. 33–46, 1994.

P. Lachapelle and S. Molotchnikoff, “Components of the electroretinogram: a reappraisal,” Documenta Ophthalmologica, vol. 63, no. 4, pp. 337–348, 1986.

K. Ikenoya, M. Kondo, C.-H. Piao, S. Kachi, Y. Miyake, and H. Terasaki, “Preservation of macular oscillatory potentials in eyes of patients with retinitis pigmentosa and normal visual acuity,” Investigative Ophthalmology and Visual Science, vol. 48, no. 7, pp. 3312–3317, 2007.

T. Sakai, M. Kondo, S. Ueno, T. Koyasu, K. Komeima, and H. Terasaki, “Supernormal ERG oscillatory potentials in transgenic rabbit with rhodopsin P347L mutation and retinal degeneration,” Investigative Ophthalmology and Visual Science, vol. 50, no. 9, pp. 4402–4409, 2009.