April 2014
Volume 55, Issue 4
Free
Retina  |   April 2014
Localizing Functional Damage in the Neural Retina of Adolescents and Young Adults With Type 1 Diabetes
Author Affiliations & Notes
  • Wylie Tan
    Department of Ophthalmology and Vision Sciences, The Hospital for Sick Children, Toronto, Canada
    School of Optometry and Vision Science, University of Waterloo, Waterloo, Canada
  • Tom Wright
    Department of Ophthalmology and Vision Sciences, The Hospital for Sick Children, Toronto, Canada
  • Annie Dupuis
    Clinical Research Services, The Hospital for Sick Children, Toronto, Canada
    Dalla Lana School of Public Health, University of Toronto, Toronto, Canada
  • Ekta Lakhani
    Department of Ophthalmology and Vision Sciences, The Hospital for Sick Children, Toronto, Canada
  • Carol Westall
    Department of Ophthalmology and Vision Sciences, The Hospital for Sick Children, Toronto, Canada
    Department of Ophthalmology and Vision Sciences, University of Toronto, Toronto, Canada
  • Correspondence: Carol Westall, Department of Ophthalmology and Vision Sciences, The Hospital for Sick Children, 555 University Avenue, Toronto, ON, Canada, M5G 1X8; [email protected]
Investigative Ophthalmology & Visual Science April 2014, Vol.55, 2432-2441. doi:https://doi.org/10.1167/iovs.13-13232
  • Views
  • PDF
  • Share
  • Tools
    • Alerts
      ×
      This feature is available to authenticated users only.
      Sign In or Create an Account ×
    • Get Citation

      Wylie Tan, Tom Wright, Annie Dupuis, Ekta Lakhani, Carol Westall; Localizing Functional Damage in the Neural Retina of Adolescents and Young Adults With Type 1 Diabetes. Invest. Ophthalmol. Vis. Sci. 2014;55(4):2432-2441. https://doi.org/10.1167/iovs.13-13232.

      Download citation file:


      © ARVO (1962-2015); The Authors (2016-present)

      ×
  • Supplements
Abstract

Purpose.: It is unknown which regions of the retina are most susceptible to damage by diabetes mellitus. We hypothesized that the standard and slow-flash (sf-) multifocal electroretinogram (mfERG) will localize retinal regions of greatest vulnerability.

Methods.: A total of 55 adolescents and young adults with type 1 diabetes and without diabetic retinopathy (DR) or with mild nonproliferative DR and 54 typically-developing, age-similar control participants underwent mfERG and sf-mfERG testing. The amplitude and implicit time of the first order response of the standard mfERG and of three multifocal oscillatory potentials (mfOPs) of the sf-mfERG were compared between groups at the level of hexagons, quadrants, and rings using separate mixed model ANOVAs. Spatial mapping of the P values from post hoc pairwise comparisons illustrated patterns of retinal dysfunction.

Results.: Delays in mfERG implicit times were evident across the tested retinal areas in the diabetes group. Delays in sf-mfERG implicit times were found at different eccentricities for each mfOP in the diabetes group. The greatest delays were noted in the periphery for mfOP1, in the midperiphery for mfOP2, and in the macular region for mfOP3. There were no significant group differences in amplitude for the mfERG and sf-mfERG protocols.

Conclusions.: Delays in mfERG and sf-mfERG responses suggest that the inner retina is particularly vulnerable to diabetes. Localizing regions of early dysfunction will help guide future studies to examine early structural damage associated with DR.

Introduction
Diabetic retinopathy (DR) is a common complication of diabetes mellitus and the leading cause of new blindness in young individuals in developed countries. 1,2 Approximately 93 million people live with some form of DR worldwide. 3 The prevalence of DR has been shown to be higher in individuals with type 1 diabetes (T1D) compared to type 2 diabetes (T2D), with a landmark study reporting 97.5% prevalence of DR among their participants with 15+ years of T1D. 4 The DR prevalence rates have since declined, 5,6 with latest global estimates reporting a prevalence of approximately 77% in this population. 3 The DR will, however, continue to be a global health concern because of the expected increase in the number of people affected by diabetes. 7  
Hallmark diagnostic features of DR consist of clinically visible microvascular lesions identified through fundus photography or ophthalmoscopy. These lesions are characterized based on the modified Airlie House classification system 8 and serve as the standard for diagnosing DR. Retinal dysfunction occurs before these vascular lesions are detected clinically, 911 suggesting that early functional disturbances in the neural retina may serve as effective biomarkers for DR. 
It has long been known that oscillatory potentials of the electroretinogram (ERG) are one of the earliest electroretinographic components to be affected and they have been suggested to be the most sensitive marker of neuroretinal dysfunction in diabetes. 911 The global response properties of the ERG cannot identify localized regions of early functional damage. The standard and slow-flash (sf-) multifocal electroretinogram (mfERG) may be more ideal tools to localize early retinal damage. Several studies have assessed localized retinal dysfunction in participants with diabetes; the majority, however, focused on adult populations with T2D or a mixture of participants with T1D and T2D. 1216 The Wisconsin Epidemiological Study of Diabetic Retinopathy (WESDR) 4 emphasized that duration of diabetes is a strong risk factor for DR and that prevalence of DR is much higher in individuals with T1D compared to T2D. Therefore, it is of interest to identify sensitive markers for DR as early as possible in the disease process, making adolescents and young adults with T1D an ideal study population. 
Few studies have investigated the retina's vulnerability to diabetes topographically. 17 This study will use the standard mfERG and sf-mfERG to assess the functional integrity of the outer and inner neural retina, respectively, in an adolescent and young adult population with T1D. The mfERG and sf-mfERG responses will be compared spatially between participants with and without T1D at the level of hexagons, quadrants, and rings. Findings from this study will identify regions of the retina that are most vulnerable to damage by the effects of diabetes. 
Methods
Subjects
A total of 55 adolescents and young adults with T1D and without DR, or with mild nonproliferative DR (NPDR) was recruited at The Hospital for Sick Children. Inclusion criteria were duration of T1D of ≥5 years and age 10 to 25 years. Participants with moderate or severe NPDR, or proliferative DR were excluded either by fundus examination by an ophthalmologist or based on seven-field, 30° stereoscopic fundus photographs graded by a retinal specialist according to the modified Airlie House classification system. 8 A total of 54 typically-developing, age-similar participants without diabetes acted as control participants. All participants with other eye diseases, hemoglobinopathy, high refractive error (worse than ±5 diopters [D]), poor visual acuity (worse than 0.3 logMAR), neurologic disorders, and those on medications affecting visual or retinal function were excluded. Informed consent was obtained from all participants after the purpose, protocol, and potential harms and benefits of the study were explained. All procedures were approved by the Research Ethics Board at The Hospital for Sick Children and were conducted in compliance with the tenets of the Declaration of Helsinki. 
Data Acquisition
All participants were tested at The Hospital for Sick Children. Since acute changes in ambient blood glucose are known to affect mfERG responses, 18,19 blood glucose readings were measured with a glucose meter (OneTouch Ultra; LifeScan, Inc., Milpitas, CA, USA) at least three times: before psychophysical testing, before mfERG testing, and after mfERG testing. Participants with T1D completed light exercises, or were administered food or insulin, as advised by a registered nurse, to adjust and maintain blood glucose levels within a 4 to 10 mmol/L range. 
One eye was selected randomly for testing from each participant and the untested eye was occluded. All participants were assessed for visual acuity (Early Treatment of Diabetic Retinopathy Study [ETDRS], logMAR) and contrast sensitivity (Pelli-Robson). Color vision was assessed with the Hardy-Rand-Rittler (HRR) pseudoisochromatic plates and the Mollon-Reffin Minimalist Test. The tested eye was anesthetized with a topical corneal anesthetic (0.5% proparacaine) and dilated pharmacologically (2.5% phenylephrine and 1% tropicamide). Refractive error was measured post dilation. Date of T1D diagnosis and hemoglobin A1c (HbA1c) values closest to the testing date were obtained from hospital records. 
Multifocal ERG
Participants were tested on the mfERG and sf-mfERG protocols with the VERIS FMSII Science System (Electro-Diagnostic Imaging, Inc., Redwood City, CA, USA) according to International Society for Clinical Electrophysiology of Vision (ISCEV) guidelines. 20 A red cross was used as the fixation target and it was enlarged as required to ensure that the target was visible during recordings. The FMSII Stimulator was adjusted to ensure that the mfERG hexagons were in focus. 
The standard mfERG protocol consisted of a stimulus array with 103 hexagons that subtended a field diameter of 40° vertical × 45° horizontal. Hexagons were scaled for eccentricity such that all responses had approximately equal amplitudes across a healthy retina. Each hexagon alternated between black (0 cd/m2) and white (200 cd/m2) following an algorithmic pseudorandom m-sequence (215 − 1) with a base rate of 13.3 ms. For every frame, each hexagon had a 50% chance of being illuminated (mean luminance 100 cd/m2). 
The sf-mfERG protocol consisted of a stimulus array of 61 hexagons that subtended a similar field to the mfERG. Hexagons were again scaled for eccentricity. Each step of the m-sequence (212 − 1) consisted of six frames: in the first frame, each hexagon had a 50% chance of being white and in the subsequent five frames, all hexagons were black. The multifocal frames were separated by 79.8 ms, allowing for the development of inner retinal responses. 
A bipolar Bürian-Allen contact lens electrode (Hansen Ophthalmic Development Laboratory, Iowa City, IA, USA) was used for recordings with a gold-plated electrode (Grass Technologies, Warwick, RI) on the forehead to serve as ground. Incoming signals were amplified (×50,000) and filtered with an analog filter (band-pass 10–300 Hz). Additional filtering was accomplished with digital filtering for mfERG (band stop filter at 60 Hz power) and sf-mfERG (75 Hz high pass filter). Two iterations of artifact removal were used and spatial averaging was not used as this option reduced the spatial accuracy of the data. Participant fixation was monitored with a built-in infrared fundus camera. Recording time was divided into 16 segments for participant comfort. Segments with fixation loss or noise artifacts were repeated. Total recording time for the mfERG (8 minutes) and sf-mfERG (∼7 minutes) protocols was approximately 15 minutes. 
Data Analysis
The outcome measures for the mfERG recordings were the amplitude and implicit time of the first order response of each hexagon. The amplitude was measured from trough-to-peak (N1 to P1) and the implicit time was measured from the time of frame presentation to the time of the major peak (Fig. 1A). 
Figure 1
 
Amplitude (dashed arrow) and implicit time (solid arrow) measurements for the standard mfERG (A) and sf-mfERG with three mfOPs (B) and two mfOPs (C) from a control participant. Gray columns represent the 5-ms time span searched for an mfOP peak.
Figure 1
 
Amplitude (dashed arrow) and implicit time (solid arrow) measurements for the standard mfERG (A) and sf-mfERG with three mfOPs (B) and two mfOPs (C) from a control participant. Gray columns represent the 5-ms time span searched for an mfOP peak.
The outcome measures for the sf-mfERG recordings were the amplitude and implicit time of the three multifocal oscillatory potentials (mfOPs). Peaks and troughs were identified with a custom script written in R Statistical Analysis Software Version 2.13.1. 21 The script identified peaks and troughs by comparing each data point with the two preceding data points and the two following data points; a data point was identified as a peak if it had the greatest positivity and a trough if it had the greatest negativity within a span of 5 ms (Fig. 1B). The amplitude of each mfOP was measured from trough-to-peak, and amplitudes smaller than 10 nV were considered noise and omitted manually. Implicit time was measured from the time of stimulus presentation to the time of the peak of each mfOP. Since the timings of the mfOPs varied with hexagon location, all control data were averaged to find the approximate timing of mfOP1, mfOP2, and mfOP3 at each hexagon. These timings acted as guidelines to classify the peaks identified for all participants. Peaks were classified as mfOP1, mfOP2, or mfOP3 based on the proximity of the peak's implicit time to the guidelines. In several cases, participants with and without diabetes had occasional hexagons that presented with two mfOPs rather than three. These two mfOPs were classified according to the guidelines and it was assumed that one of the three mfOPs was missing (Fig. 1C). It usually was mfOP2 or mfOP3 that was missing, and rarely mfOP1. Almost all participants had at least 1 of 61 hexagons with a missing mfOP and the maximum number of hexagons with a missing mfOP in a single participant was approximately 30 hexagons. 
All left eye recordings were converted to the right eye, such that temporal regions of the left eye were compared to temporal regions of the right eye and similarly for the nasal regions. 
Mixed Model ANOVA
To examine spatial differences in retinal function, between-group comparisons were made at the level of hexagons, quadrants, and rings with separate mixed model ANOVAs. The mfERG and sf-mfERG stimulus patterns were divided into four quadrants (along horizontal and vertical meridians from the fovea) and rings (six for mfERG and five for sf-mfERG, Fig. 2). Quadrant and ring amplitude, and implicit time responses were comprised of the average of the amplitude and implicit time of all of the hexagons within each region. For mfERG and sf-mfERG mixed model ANOVAs, the independent variable was group (T1D or control), the dependent variable was amplitude or implicit time, and the repeated measure variable was hexagon, quadrant, or ring. 
Figure 2
 
Stimulus array for mfERG (AC) and sf-mfERG (DF), divided into hexagons (A, D), quadrants (B, E), and rings (C, F).
Figure 2
 
Stimulus array for mfERG (AC) and sf-mfERG (DF), divided into hexagons (A, D), quadrants (B, E), and rings (C, F).
Mixed model analyses were performed using SAS version 9.2 (SAS Institute, Inc., Rockville, MD, USA). Group means were compared using post hoc pairwise comparisons. Bonferroni corrections were applied to account for the multiple comparisons: α = 0.05 was divided by the number of comparisons or levels within the repeated measure variable. The Bonferroni adjustment is most appropriate when hypothesis tests are independent. As observations across hexagons, quadrants, and rings were correlated, the Bonferroni correction in this case is overly conservative. Thus, all results are presented as actual P values for the reader's discretion. 
Results
Demographic data and psychophysical testing results for all participants are shown in the Table. There were no significant between-group differences in age, visual acuity, or contrast sensitivity. Two participants with T1D had mild red-green color vision defects. Most participants with T1D (45/53) were diagnosed before the age of 10 and approximately 70% had been diagnosed with T1D for ≥8 years. Date of diagnosis was not available for two participants. The HbA1c values were not available for six participants because they were over the age of 18 and had left The Hospital for Sick Children. Over 89% of the participants with T1D had HbA1c levels above the recommended HbA1c target of 7% set by the Canadian Diabetes Association. 22 Eight participants with diabetes had early NPDR. After qualitative comparisons between the fundus photographs and mfERG traces of these eight participants, no correlations between morphologic NPDR findings and local mfERG delays were found. 
Table
 
Demographic Data and Psychophysical Testing Results for Participants with T1D and Control Participants
Table
 
Demographic Data and Psychophysical Testing Results for Participants with T1D and Control Participants
T1D Participants, n = 55 Control Participants, n = 54
Sex, male/female 26/29 19/35
NPDR, n 8
Age at testing, y 16.6 ± 2.3 (12.1–22.1) 17.1 ± 3.6 (10.7–25.1)
Age at diagnosis, y 7.2 ± 3.9 (1.4–16.3)
Duration of T1D, y 9.5 ± 3.1 (3.9–15.6)
HbA1c, % 8.6 ± 1.3 (6.4–12.3)
Visual acuity, logMAR −0.01 ± 0.12 (−0.28–0.28) −0.07 ± 0.13 (−0.60–0.24)
Contrast sensitivity 1.68 ± 0.13 (1.05–2.05) 1.69 ± 0.16 (1.05–2.20)
Data from most of the participants, with and without diabetes, have been included in prior reports. 2325 This study uses a different analysis technique to assess the spatial relationships of the mfERG and sf-mfERG data, and uses the peak-picking technique to analyze waveforms rather than waveform fitting techniques. 
The mixed model analyses were conducted for the mfERG and sf-mfERG data. Rerunning the models for the sf-mfERG data after removing the observations with two mfOPs only yielded similar results, possibly because the number of hexagons with only two mfOPs was relatively small (accounting for 15% and 13% of all data for participants with and without diabetes, respectively) compared to the hexagons with three mfOPs. Rerunning the mixed models after removing the participants with mild NPDR also returned the same results. The analyses presented herein included all of the mfERG and sf-mfERG data from participants with and without mild NPDR. 
Mixed Models for mfERG
Mixed model analyses for mfERG implicit time yielded main effects for group at the level of hexagons (P < 0.0001), quadrants (P < 0.0001), and rings (P < 0.0001), but not for amplitude. The P values for the post hoc pairwise comparisons from the hexagon analysis are presented in a spatial color map to illustrate the distribution of dysfunction (Fig. 3). 
Figure 3
 
Spatial color map representing the mfERG implicit time differences between groups at the level of hexagons (illustrated on an unscaled hexagon array).
Figure 3
 
Spatial color map representing the mfERG implicit time differences between groups at the level of hexagons (illustrated on an unscaled hexagon array).
Positive (+) and negative (−) signs denote whether participants with T1D had earlier or delayed implicit times, respectively, compared to control participants. Implicit times were delayed in participants with T1D compared to control participants at most locations (>60%, α = 0.05). Implicit times of participants with T1D were significantly delayed at three hexagons only with the Bonferroni correction (α = 0.0005). It is important to note, however, that the Bonferroni correction is the most conservative correction for multiple comparisons. Spatial mapping helps to visualize the number of hexagons that are approaching significance. Over 50% of the hexagons have a P < 0.02, indicating that implicit times of participants with T1D are delayed compared to controls across the majority of the retina. Figure 4 shows representative mfERG trace arrays from 2 participants with diabetes. 
Figure 4
 
Examples of representative mfERG trace arrays from two participants (A, B) with diabetes. Blue traces represent individual patient data and dashed black traces represent average control data. Shaded hexagons indicate responses that are significantly delayed (>2 SD from controls) (illustrated on unscaled hexagon arrays).
Figure 4
 
Examples of representative mfERG trace arrays from two participants (A, B) with diabetes. Blue traces represent individual patient data and dashed black traces represent average control data. Shaded hexagons indicate responses that are significantly delayed (>2 SD from controls) (illustrated on unscaled hexagon arrays).
Quadrant (Fig. 5) and ring (Fig. 6) analyses showed that implicit time of participants with T1D were significantly delayed compared to control participants in all quadrants and eccentricities except at the fovea. These comparisons were significant after the Bonferroni correction (α = 0.01 for quadrants and α = 0.0008 for rings). 
Figure 5
 
Comparison of mfERG implicit time between groups at the level of quadrants. Data presented as mean ± 95% confidence limits. Bonferroni correction, P < 0.01.
Figure 5
 
Comparison of mfERG implicit time between groups at the level of quadrants. Data presented as mean ± 95% confidence limits. Bonferroni correction, P < 0.01.
Figure 6
 
Comparison of mfERG implicit time between groups at the level of rings. Data presented as mean ± 95% confidence limits. Bonferroni correction, P < 0.0008.
Figure 6
 
Comparison of mfERG implicit time between groups at the level of rings. Data presented as mean ± 95% confidence limits. Bonferroni correction, P < 0.0008.
Mixed Models for sf-mfERG
Mixed model analyses for mfOP1, mfOP2, and mfOP3 of the sf-mfERG yielded significant between-group differences for implicit time, but not for amplitude. Significant main effects were found for mfOP1 at the level of hexagons (P = 0.0003) and quadrants (P < 0.0001). Significant ring differences were found for mfOP2 (P = 0.007) and mfOP3 (P = 0.0008), mfOP1 approached significance (P = 0.06). 
The P values for the post hoc pairwise comparisons at the level of hexagons were mapped spatially in color plots (Fig. 7). The areas of interest are those in which the implicit time of participants with T1D was delayed significantly compared to control participants at α = 0.05 (color toward the red end of the spectrum with negative sign). The distribution of these delayed regions was different for each mfOP: the red hexagons were concentrated in the periphery for mfOP1, in the midperiphery for mfOP2, and in the macular region for mfOP3. Although there are minimal areas with significant differences with the Bonferroni correction (α = 0.0008), it is important to note the pattern in the distribution of P values. Ring analysis showed a similar distribution of delays (Fig. 7): the greatest implicit time delays were observed in Ring 4 for mfOP1 (P < 0.0001), Ring 3 for mfOP2 (P < 0.0001), and Ring 2 for mfOP3 (P < 0.0001). These differences are significant with the Bonferroni correction (α = 0.01). Figure 8 shows representative sf-mfERG trace arrays from two participants with diabetes. 
Figure 7
 
Spatial color maps illustrating the distribution of mfOP implicit time differences between groups: both hexagon and ring analyses show that each mfOP is most delayed at different eccentricities (illustrated on unscaled hexagon arrays).
Figure 7
 
Spatial color maps illustrating the distribution of mfOP implicit time differences between groups: both hexagon and ring analyses show that each mfOP is most delayed at different eccentricities (illustrated on unscaled hexagon arrays).
Figure 8
 
Representative sf-mfERG trace arrays for two participants with diabetes (A, B). Blue traces represent individual patient data and dashed black traces represent average control data. The same trace array is repeated to demonstrate the hexagons with delayed implicit times for mfOP1 (green), mfOP2 (yellow), mfOP3 (red); the leftmost arrays show a composite of all hexagons with implicit time delays (illustrated on unscaled hexagon arrays).
Figure 8
 
Representative sf-mfERG trace arrays for two participants with diabetes (A, B). Blue traces represent individual patient data and dashed black traces represent average control data. The same trace array is repeated to demonstrate the hexagons with delayed implicit times for mfOP1 (green), mfOP2 (yellow), mfOP3 (red); the leftmost arrays show a composite of all hexagons with implicit time delays (illustrated on unscaled hexagon arrays).
Quadrant analysis showed that mfOP1 is consistently more delayed in participants with T1D compared to control participants in all four retinal quadrants: superior temporal (P < 0.0001), superior nasal (P = 0.007), inferior temporal (P = 0.004), inferior nasal (P = 0.001). 
Discussion
This study found no change in mfERG amplitude, but mfERG implicit times were delayed across the retina in the diabetes group. This is consistent with mfERG literature, where there has been no consensus as to whether mfERG amplitude is affected by diabetes, 14,26,27 but multifocal ERG implicit times has been reported consistently as delayed in participants with diabetes, 14,23,27 and delayed to greater extents with increasing DR severity. 26,28  
Quadrant and ring analyses of the mfERG data revealed that the majority of the retina is affected by diabetes, as participants with T1D showed response delays across all quadrants and eccentricities. Bearse et al. 12 similarly studied and mapped retinal dysfunction in participants with diabetes, and found a spatial preference of vulnerability to the effects of diabetes. They used an sf-mfERG protocol, rather than the mfERG protocol, to evaluate first-order N1, P1, and N2 implicit times, and found that P1 implicit time delays are most frequent in the inferior and peripheral retina, and, thus, concluded these regions are most vulnerable to damage. A few key factors differentiate our studies, including differences in the protocol used, the population studied, and the analysis technique. Bearse et al. 12 studied an adult population with a mix of participants with T1D and T2D, whereas our study examined an adolescent and young adult population with T1D only. The advantages of studying biomarkers with a younger rather than an older adult population have been discussed previously. 4 Bearse et al. 12 plotted the frequency of delays at each hexagon and defined the region most frequently associated with abnormalities as the region of greatest vulnerability. In our study, the region of greatest vulnerability was defined as the region in which the mean implicit time of participants with T1D differed most significantly with the mean implicit time of control participants. 
A comprehensive review conducted by Hood 29 examined studies that used the mfERG to assess different retinal diseases. He derived a guide to locate the site of damage within the retina based on characteristic changes to the mfERG waveform. No amplitude changes and a small implicit time delay (<3 ms) are indicative of damage within the inner plexiform layer. 29 In our study, there was no change in amplitude between groups and subtle (<1 ms), but significant implicit time delays in the diabetes group, suggesting dysfunction at the level of the inner plexiform layer in our participants with T1D. 
Inner retinal function was assessed with the sf-mfERG protocol and participants with T1D showed implicit time delays again. Each mfOP was delayed at different eccentricities, with the greatest delays in the periphery for mfOP1, in the midperiphery for mfOP2, and in the macular region for mfOP3. Delays in mfOPs have been reported consistently in studies involving participants with diabetes. 17,30,31 Kurtenbach et al. 17 studied mfOP topography in adolescents with T1D with the sf-mfERG. A few features in their study differentiate it from our study, including differences in quadrant divisions, the identification of two mfOPS in the first order response instead of three mfOPs, a smaller stimulus size, and a faster flash sequence. Also, they did not take multiple comparisons into consideration or control for blood sugar levels in their participants with diabetes, which have been shown to affect the mfERG. 19,32 Despite these differences, their ring analysis yielded results that were similar to our findings: significant implicit time delays were found peripherally for mfOP1 (∼17°–30°) and midperipherally for mfOP2 (∼11°–22°). 
The exact origin of mfOPs is unknown, but past microelectrode studies have isolated the responses to the inner plexiform layer and suggested that each mfOP may have different cellular origins. 33 Our finding that each of the mfOP delays varied with eccentricity supports the idea of different cellular origins. In the primate retina, the mfOPs have been described as high frequency mfOPs and low frequency mfOPs. 34 High frequency mfOPs correspond to early mfOPs (OP1-3) and are generated mainly by ganglion cells. Low frequency mfOPs correspond to early and late mfOPs, and are mainly generated by amacrine cells and hyperpolarizing second-order cells (OFF-bipolar and horizontal cells). Thus, the three mfOPs in the current study may originate from ganglion, amacrine, and interplexiform cells, and the mfOPs may reflect negative feedback activity between these cell types, but further investigation is needed to identify the specific origins of these mfOPs and the mechanism involved in the generation of the observed patterns of dysfunction. 
Our findings from the mfERG and sf-mfERG protocols implicated the inner plexiform layer as a region particularly vulnerable to the effects of diabetes. Many studies have suggested that hypoxia may have a key role in inner retinal dysfunction. 10,35 The inner retina is supplied by the retinal circulation, and experimental disruption of this circulation in animal studies results in hypoxia and reduces or eliminates OPs. 9,36 Altered rod photoreceptor function has been reported in participants with diabetes. 37 It has been postulated that because rods are compromised in diabetes, the oxygen requirements for dark adaptation increase. The outer retina requires more oxygen and imposes additional hypoxic stresses on the inner retina, resulting in inner retinal dysfunction. 37 In addition, early changes in retinal astrocytes have been identified in the inner and peripheral retina of diabetes-induced rats, and these changes are coincident with early retinal dysfunction and vascular changes. 38 Since astrocytes are glial cells that modulate neuronal and vascular function, their alterations may have a key role in inner retinal dysfunction. 38 Altered glial cells, hypoxia, and retinal stress may have contributed to the retinal dysfunction observed in our participants with T1D. 
Certain factors may have influenced our results. The majority of our study population was at or around the age of puberty. Puberty is a risk factor for the development of DR. 22,24,39 Hormonal fluctuations associated with puberty pose challenges for those with diabetes in controlling their blood glucose levels. 22 These fluctuations may affect retinal function. One of the strengths of our study is that ambient blood glucose was controlled for and monitored throughout testing. Other factors that affect the mfERG are age and high refractive errors. Multifocal ERG amplitude decreases with increasing age and refractive errors. 10,40 Our participants with T1D were compared to age-similar control participants and the majority of participants were refracted. For the participants who were not refracted (15/55 patients and 26/54 controls), a refractive error worse than ±5 D was ruled out if the participant had visual acuity better than 0.3 logMAR. 
The main findings of this study offer insight about the specific regions of the retina that are most susceptible to damage from the effects of diabetes. The inner retina is particularly vulnerable and may suffer the effects of hypoxia, leading to early vascular changes. Spatial mapping of local dysfunction is important because identifying the regions of greatest vulnerability greatly increases the sensitivity of finding early structural changes in the retina. 
Acknowledgments
The authors thank Wai Ching Lam, MD, FRCSC, for grading fundus photographs, Marcia Wilson, RN, for titrating and monitoring the blood glucose levels of participants with diabetes during testing, Melissa Cotesta, OC(C), for conducting refractions, Cynthia VandenHoven, BAA, CRA, for fundus photography, Cheng Lin, BASc, for technical assistance with creating spatial color maps, and Carole Panton, DBO(D), CO, OC(C), for editing the manuscript. 
Supported by the Canadian Institute for Health Research (CW), Canadian Foundation for Innovation (CW), Vision Science Research Program Graduate Student Scholarship (WT), Banting and Best Diabetes Novo Nordisk Graduate Scholarship (WT), and SickKids Restracomp Graduate Scholarship (WT). 
Disclosure: W. Tan, None; T. Wright, None; A. Dupuis, None; E. Lakhani, None; C. Westall, None 
References
Cheung N Mitchell P Wong TY. Diabetic retinopathy. Lancet . 2010; 376: 124–136. [CrossRef] [PubMed]
Congdon NG Friedman DS Lietman T. Important causes of visual impairment in the world today. JAMA . 2003; 290: 2057–2060. [CrossRef] [PubMed]
Yau JWY Rogers SL Kawasaki R Global prevalence and major risk factors of diabetic retinopathy. Diabetes Care . 2012; 35: 556–564. [CrossRef] [PubMed]
Klein R Klein BEK Moss SE Davis MD Demets DL. The Wisconsin Epidemiology Study of Diabetic Retinopathy. II. Prevalence and risk of diabetic retinopathy when age at diagnosis is less than 30 years. Arch Ophthalmol . 1984; 102: 520–526. [CrossRef] [PubMed]
LeCaire T Palta M Zhang H Allen C Klein R D'Alessio D. Lower-than-expected prevalence and severity of retinopathy in an incident cohort followed during the first 4–14 years of type 1 diabetes: the Wisconsin Diabetes Registry Study. Am J Epidemiol . 2006; 164: 143–150. [CrossRef] [PubMed]
Downie E Craig ME Hing S Cusumano J Chan AKF Donaghue KC. Continued reduction in the prevalence of retinopathy in adolescents with type 1 diabetes: role of insulin therapy and glycemic control. Diabetes Care . 2011; 34: 2368–2373. [CrossRef] [PubMed]
International Diabetes Federation Annual Report 2011 . Brussels, Belgium: International Diabetes Federation; 2011: 1–36. Available at: http://www.idf.org/files/idf_publications/annual_report_EN/index.html. Accessed June 27, 2012.
Early Treatment Diabetic Retinopathy Study Research Group. Grading diabetic retinopathy from stereoscopic color fundus photographs: an extension of the modified Airlie House classification. ETDRS report number 10. Ophthalmology . 1991; 98: 786–806. [CrossRef] [PubMed]
Bresnick GH Korth K Groo A Palta M. Electroretinographic oscillatory potentials predict progression of diabetic retinopathy preliminary report. Arch Ophthalmol . 1984; 102: 1307–1311. [CrossRef] [PubMed]
Bresnick GH Palta M. Oscillatory potential amplitudes. relation to severity of diabetic retinopathy. Arch Ophthalmol . 1987; 105: 929–933. [CrossRef] [PubMed]
Speros P Price J. Oscillatory potentials. History, techniques and potential use in the evaluation of disturbances of retinal circulation. Surv Ophthalmol . 1981; 25: 237–252. [CrossRef] [PubMed]
Bearse MA Jr Han Y Schneck ME Adams AJ. Retinal function in normal and diabetic eyes mapped with the slow flash multifocal electroretinogram. Invest Ophthalmol Vis Sci . 2004; 45: 296–304. [CrossRef] [PubMed]
Bearse MA Jr Han Y Schneck ME Barez S Jacobsen C Adams AJ. Local multifocal oscillatory potential abnormalities in diabetes and early diabetic retinopathy. Invest Ophthalmol Vis Sci . 2004; 45: 3259–3265. [CrossRef] [PubMed]
Palmowski AM Sutter EE Bearse MA Fung W. Mapping of retinal function in diabetic retinopathy using the multifocal electroretinogram. Invest Ophthalmol Vis Sci . 1997; 38: 2586–2596. [PubMed]
Han Y Schneck ME Bearse MA Jr Formulation and evaluation of a predictive model to identify the sites of future diabetic retinopathy. Invest Ophthalmol Vis Sci . 2004; 45: 4106–4112. [CrossRef] [PubMed]
Harrison WW Bearse MA Ng JS Multifocal electroretinograms predict onset of diabetic retinopathy in adult patients with diabetes. Invest Ophthalmol Vis Sci . 2011; 52: 772–777. [CrossRef] [PubMed]
Kurtenbach A Langrova H Zrenner E. Multifocal oscillatory potentials in type 1 diabetes without retinopathy. Invest Ophthalmol Vis Sci . 2000; 41: 3234–3241. [PubMed]
Khan MI Barlow RB Weinstock RS. Acute hypoglycemia decreases central retinal function in the human eye. Vision Res . 2011; 51: 1623–1626. [CrossRef] [PubMed]
Klemp K Larsen M Sander B Vaag A Brockhoff PB Lund-Andersen H. Effect of short-term hyperglycemia on multifocal electroretinogram in diabetic patients without retinopathy. Invest Ophthalmol Vis Sci . 2004; 45: 3812–3819. [CrossRef] [PubMed]
Hood DC Bach M Brigell M Keating D Kondo M Lyons JS. ISCEV standard for clinical multifocal electroretinography (2011 edition). Doc Ophthalmol . 2012; 124: 1–13. [CrossRef] [PubMed]
R Development Core Team R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2011. Available at: http://www.R-project.org/. Accessed July 25, 2011.
Canadian Diabetes Association. Canadian Diabetes Association 2008 clinical practice guidelines for the prevention and management of diabetes in Canada. Can J Diabetes . 2008; 32 (suppl 1): S1–S201.
Wright T Cortese F Nilsson J Westall C. Analysis of multifocal electroretinograms from a population with type 1 diabetes using partial least squares reveals spatial and temporal distribution of changes to retinal function. Doc Ophthalmol . 2012; 125: 31–42. [CrossRef] [PubMed]
Lakhani E Wright T Abdolell M Westall C. Multifocal ERG defects associated with insufficient long-term glycemic control in adolescents with type 1 diabetes. Invest Ophthalmol Vis Sci . 2010; 51: 5297–5303. [CrossRef] [PubMed]
McFarlane M Wright T Stephens D Nilsson J Westall CA. Blue flash ERG PhNR changes associated with poor long-term glycemic control in adolescents with type 1 diabetes. Invest Ophthalmol Vis Sci . 2012; 53: 741–748. [CrossRef] [PubMed]
Fortune B Schneck ME Adams AJ. Multifocal electroretinogram delays reveal local retinal dysfunction in early diabetic retinopathy. Invest Ophthalmol Vis Sci . 1999; 40: 2638–2651. [PubMed]
Tyrberg M Ponjavic V Lövestam-Adrian M. Multifocal electroretinography (mfERG) in insulin dependent diabetics with and without clinically apparent retinopathy. Doc Ophthalmol . 2005; 110: 137–143. [CrossRef] [PubMed]
Schneck ME Bearse MA Jr Han Y Barez S Jacobsen C Adams AJ. Comparison of mfERG waveform components and implicit time measurement techniques for detecting functional change in early diabetic eye disease. Doc Ophthalmol . 2004; 108: 223–230. [CrossRef] [PubMed]
Hood DC. Assessing retinal function with the multifocal technique. Prog Retinal Eye Res . 2000; 19: 607–646. [CrossRef]
Klemp K Lund-Andersen H Sander B Larsen M. The effect of acute hypoxia and hyperoxia on the slow multifocal electroretinogram in healthy subjects. Invest Ophthalmol Vis Sci . 2007; 48: 3405–3412. [CrossRef] [PubMed]
Onozu H Yamamoto S. Oscillatory potentials of multifocal electroretinogram in diabetic retinopathy. Doc Ophthalmol . 2003; 106: 327–332. [CrossRef] [PubMed]
Klemp K Sander B Brockhoff PB Vaag A Lund-Andersen H Larsen M. The multifocal ERG in diabetic patients without retinopathy during euglycemic clamping. Invest Ophthalmol Vis Sci . 2005; 46: 2620–2626. [CrossRef] [PubMed]
Wachtmeister L Dowling JE. The oscillatory potentials of the mudpuppy retina. Invest Ophthalmol Vis Sci . 1978; 17: 1176–1188. [PubMed]
Zhou W Rangaswamy N Ktonas P Frishman LJ. Oscillatory potentials of the slow-sequence multifocal ERG in primates extracted using the Matching Pursuit method. Vision Res . 2007; 47: 2021–2036. [CrossRef] [PubMed]
Greenstein VC Holopigian K Hood DC Seiple W Carr RE. The nature and extent of retinal dysfunction associated with diabetic macular edema. Invest Ophthalmol Vis Sci . 2000; 41: 3643–3654. [PubMed]
Brown KT. The electroretinogram: its components and their origins. Vision Res . 1968; 8: 633–677. [CrossRef] [PubMed]
Arden GB Wolf JE Tsang Y. Does dark adaptation exacerbate diabetic retinopathy? Evidence and a linking hypothesis. Vision Res . 1998; 38: 1723–1729. [CrossRef] [PubMed]
Ly A Yee P Vessey KA Phipps JA Jobling AI Fletcher EL. Early inner retinal astrocyte dysfunction during diabetes and development of hypoxia, retinal stress, and neuronal functional loss. Invest Ophthalmol Vis Sci . 2011; 52: 9316–9326. [CrossRef] [PubMed]
Donaghue KC Fairchild JM Craig ME Do all prepubertal years of diabetes duration contribute equally to diabetes complications? Diabetes Care . 2003; 26: 1224–1229. [CrossRef] [PubMed]
Palmowski AM Berninger T Allgayer R Andrielis H Heinemann-Vernaleken B Rudolph G. Effects of refractive blur on the multifocal electroretinogram. Doc Ophthalmol . 1999; 99: 41–54. [CrossRef] [PubMed]
Figure 1
 
Amplitude (dashed arrow) and implicit time (solid arrow) measurements for the standard mfERG (A) and sf-mfERG with three mfOPs (B) and two mfOPs (C) from a control participant. Gray columns represent the 5-ms time span searched for an mfOP peak.
Figure 1
 
Amplitude (dashed arrow) and implicit time (solid arrow) measurements for the standard mfERG (A) and sf-mfERG with three mfOPs (B) and two mfOPs (C) from a control participant. Gray columns represent the 5-ms time span searched for an mfOP peak.
Figure 2
 
Stimulus array for mfERG (AC) and sf-mfERG (DF), divided into hexagons (A, D), quadrants (B, E), and rings (C, F).
Figure 2
 
Stimulus array for mfERG (AC) and sf-mfERG (DF), divided into hexagons (A, D), quadrants (B, E), and rings (C, F).
Figure 3
 
Spatial color map representing the mfERG implicit time differences between groups at the level of hexagons (illustrated on an unscaled hexagon array).
Figure 3
 
Spatial color map representing the mfERG implicit time differences between groups at the level of hexagons (illustrated on an unscaled hexagon array).
Figure 4
 
Examples of representative mfERG trace arrays from two participants (A, B) with diabetes. Blue traces represent individual patient data and dashed black traces represent average control data. Shaded hexagons indicate responses that are significantly delayed (>2 SD from controls) (illustrated on unscaled hexagon arrays).
Figure 4
 
Examples of representative mfERG trace arrays from two participants (A, B) with diabetes. Blue traces represent individual patient data and dashed black traces represent average control data. Shaded hexagons indicate responses that are significantly delayed (>2 SD from controls) (illustrated on unscaled hexagon arrays).
Figure 5
 
Comparison of mfERG implicit time between groups at the level of quadrants. Data presented as mean ± 95% confidence limits. Bonferroni correction, P < 0.01.
Figure 5
 
Comparison of mfERG implicit time between groups at the level of quadrants. Data presented as mean ± 95% confidence limits. Bonferroni correction, P < 0.01.
Figure 6
 
Comparison of mfERG implicit time between groups at the level of rings. Data presented as mean ± 95% confidence limits. Bonferroni correction, P < 0.0008.
Figure 6
 
Comparison of mfERG implicit time between groups at the level of rings. Data presented as mean ± 95% confidence limits. Bonferroni correction, P < 0.0008.
Figure 7
 
Spatial color maps illustrating the distribution of mfOP implicit time differences between groups: both hexagon and ring analyses show that each mfOP is most delayed at different eccentricities (illustrated on unscaled hexagon arrays).
Figure 7
 
Spatial color maps illustrating the distribution of mfOP implicit time differences between groups: both hexagon and ring analyses show that each mfOP is most delayed at different eccentricities (illustrated on unscaled hexagon arrays).
Figure 8
 
Representative sf-mfERG trace arrays for two participants with diabetes (A, B). Blue traces represent individual patient data and dashed black traces represent average control data. The same trace array is repeated to demonstrate the hexagons with delayed implicit times for mfOP1 (green), mfOP2 (yellow), mfOP3 (red); the leftmost arrays show a composite of all hexagons with implicit time delays (illustrated on unscaled hexagon arrays).
Figure 8
 
Representative sf-mfERG trace arrays for two participants with diabetes (A, B). Blue traces represent individual patient data and dashed black traces represent average control data. The same trace array is repeated to demonstrate the hexagons with delayed implicit times for mfOP1 (green), mfOP2 (yellow), mfOP3 (red); the leftmost arrays show a composite of all hexagons with implicit time delays (illustrated on unscaled hexagon arrays).
Table
 
Demographic Data and Psychophysical Testing Results for Participants with T1D and Control Participants
Table
 
Demographic Data and Psychophysical Testing Results for Participants with T1D and Control Participants
T1D Participants, n = 55 Control Participants, n = 54
Sex, male/female 26/29 19/35
NPDR, n 8
Age at testing, y 16.6 ± 2.3 (12.1–22.1) 17.1 ± 3.6 (10.7–25.1)
Age at diagnosis, y 7.2 ± 3.9 (1.4–16.3)
Duration of T1D, y 9.5 ± 3.1 (3.9–15.6)
HbA1c, % 8.6 ± 1.3 (6.4–12.3)
Visual acuity, logMAR −0.01 ± 0.12 (−0.28–0.28) −0.07 ± 0.13 (−0.60–0.24)
Contrast sensitivity 1.68 ± 0.13 (1.05–2.05) 1.69 ± 0.16 (1.05–2.20)
×
×

This PDF is available to Subscribers Only

Sign in or purchase a subscription to access this content. ×

You must be signed into an individual account to use this feature.

×