Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Detecting retinal neurodegeneration in people with diabetes: Findings from the UK Biobank

  • Roomasa Channa ,

    Roles Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Visualization, Writing – original draft, Writing – review & editing

    rchanna@wisc.edu

    Affiliation Department of Ophthalmology and Visual Sciences, University of Wisconsin, Madison, WI, United States of America

  • Kyungmoo Lee,

    Roles Formal analysis, Software, Writing – review & editing

    Affiliation Department of Ophthalmology and Visual Sciences, University of Iowa, Iowa City, Iowa, United States of America

  • Kristen A. Staggers,

    Roles Formal analysis, Writing – original draft, Writing – review & editing

    Affiliation Institute for Clinical and Translational Research, Baylor College of Medicine, Houston, TX, United States of America

  • Nitish Mehta,

    Roles Formal analysis, Writing – original draft, Writing – review & editing

    Affiliation New York University, New York, NY, United States of America

  • Sidra Zafar,

    Roles Writing – original draft, Writing – review & editing

    Affiliation Wilmer Eye Institute, Johns Hopkins University, Baltimore, MD, United States of America

  • Jie Gao,

    Roles Formal analysis, Investigation, Software, Writing – original draft

    Affiliation Department of Ophthalmology, Baylor College of Medicine, Houston, TX, United States of America

  • Benjamin J. Frankfort,

    Roles Conceptualization, Formal analysis, Resources, Supervision, Writing – original draft, Writing – review & editing

    Affiliation Department of Ophthalmology and Neurosciences, Baylor College of Medicine, Houston, TX, United States of America

  • Sharon Y. L. Chua,

    Roles Formal analysis, Writing – review & editing

    Affiliation NIHR Biomedical Research Centre at Moorfields Eye Hospital NHS Foundation Trust & UCL Institute of Ophthalmology, London, United Kingdom

  • Anthony P. Khawaja,

    Roles Conceptualization, Writing – review & editing

    Affiliation NIHR Biomedical Research Centre at Moorfields Eye Hospital NHS Foundation Trust & UCL Institute of Ophthalmology, London, United Kingdom

  • Paul J. Foster,

    Roles Formal analysis, Funding acquisition, Writing – review & editing

    Affiliation NIHR Biomedical Research Centre at Moorfields Eye Hospital NHS Foundation Trust & UCL Institute of Ophthalmology, London, United Kingdom

  • Praveen J. Patel,

    Roles Funding acquisition, Writing – review & editing

    Affiliation NIHR Biomedical Research Centre at Moorfields Eye Hospital NHS Foundation Trust & UCL Institute of Ophthalmology, London, United Kingdom

  • Charles G. Minard,

    Roles Formal analysis, Writing – review & editing

    Affiliation Institute for Clinical and Translational Research, Baylor College of Medicine, Houston, TX, United States of America

  • Chris Amos,

    Roles Funding acquisition, Resources, Supervision, Writing – review & editing

    Affiliation Institute for Clinical and Translational Research, Baylor College of Medicine, Houston, TX, United States of America

  • Michael D. Abramoff

    Roles Conceptualization, Resources, Supervision, Visualization, Writing – review & editing

    Affiliation Department of Ophthalmology and Visual Sciences, University of Iowa, Iowa City, Iowa, United States of America

Abstract

Importance

Efforts are underway to incorporate retinal neurodegeneration in the diabetic retinopathy severity scale. However, there is no established measure to quantify diabetic retinal neurodegeneration (DRN).

Objective

We compared total retinal, macular retinal nerve fiber layer (mRNFL) and ganglion cell-inner plexiform layer (GC-IPL) thickness among participants with and without diabetes (DM) in a population-based cohort.

Design/setting/participants

Cross-sectional analysis, using the UK Biobank data resource. Separate general linear mixed models (GLMM) were created using DM and glycated hemoglobin as predictor variables for retinal thickness. Sub-analyses included comparing thickness measurements for patients with no/mild diabetic retinopathy (DR) and evaluating factors associated with retinal thickness in participants with and without diabetes. Factors found to be significantly associated with DM or thickness were included in a multiple GLMM.

Exposure

Diagnosis of DM was determined via self-report of diagnosis, medication use, DM-related complications or glycated hemoglobin level of ≥ 6.5%.

Main outcomes and measures

Total retinal, mRNFL and GC-IPL thickness.

Results

74,422 participants (69,985 with no DM; 4,437 with DM) were included. Median age was 59 years, 46% were men and 92% were white. Participants with DM had lower total retinal thickness (-4.57 μm, 95% CI: -5.00, -4.14; p<0.001), GC-IPL thickness (-1.73 μm, 95% CI: -1.86, -1.59; p<0.001) and mRNFL thickness (-0.68 μm, 95% CI: -0.81, -0.54; p<0.001) compared to those without DM. After adjusting for co-variates, in the GLMM, total retinal thickness was 1.99 um lower (95% CI: -2.47, -1.50; p<0.001) and GC-IPL was 1.02 μm lower (95% CI: -1.18, -0.87; p<0.001) among those with DM compared to without. mRNFL was no longer significantly different (p = 0.369). GC-IPL remained significantly lower, after adjusting for co-variates, among those with DM compared to those without DM when including only participants with no/mild DR (-0.80 μm, 95% CI: -0.98, -0.62; p<0.001). Total retinal thickness decreased 0.40 μm (95% CI: -0.61, -0.20; p<0.001), mRNFL thickness increased 0.20 μm (95% CI: 0.14, 0.27; p<0.001) and GC-IPL decreased 0.26 μm (95% CI: -0.33, -0.20; p<0.001) per unit increase in A1c after adjusting for co-variates. Among participants with diabetes, age, DR grade, ethnicity, body mass index, glaucoma, spherical equivalent, and visual acuity were significantly associated with GC-IPL thickness.

Conclusion

GC-IPL was thinner among participants with DM, compared to without DM. This difference persisted after adjusting for confounding variables and when considering only those with no/mild DR. This confirms that GC-IPL thinning occurs early in DM and can serve as a useful marker of DRN.

Introduction

Globally there are an estimated 93 million people with diabetic retinopathy (DR) and 21 million with diabetic macular edema (DME) [1]. Vision loss from diabetic retinal diseases is thought to be primarily due to vascular damage, which manifests as DR or DME. However, there is increasing evidence of an underlying neurodegenerative process [2, 3] and recognition regarding the need to update the DR severity scale to incorporate diabetic retinal neurodegeneration(DRN) [4, 5]. Spectral Domain Optical Coherence Tomography (SD-OCT) is the most widely available imaging technique to quantify structural changes in the neuroretina. While prospective studies have shown that DRN occurs before DR manifests [6], results from cross-sectional studies evaluating retinal thickness parameters to quantify DRN have reported mixed results, have included small samples and have not always controlled for the many factors that can impact thickness measurements including, age, gender, ethnicity, refractive error, cognitive impairment, glaucoma, systemic co-morbidities and medications [710]. Consequently, some studies report increased thickness of inner retinal layers, some show no changes and others suggest that the RNFL and GC-IPL may be thinner in patients with DM compared to controls [11, 12]. Lim et al reported a higher rate of peripapillary RNFL thinning among patients with DM compared to those without DM [13]. Although this information is useful, clinically it is helpful to identify changes in the macula as macula-centered OCT scans are already routinely used to evaluate macular edema in patients with diabetes. A recent study of 112 patients with diabetes and 63 healthy controls reported a higher rate of GC-IPL thinning among patients with diabetes, but did not report changes in the macular RNFL [14]. Also, given the small number of patients in this study it is difficult to fully evaluate factors associated with inner retinal thickness changes in patients with diabetes. If we are to incorporate evaluation of DRN in clinical practice, using OCT scans, it is important to fully understand factors associated with inner retinal thickness measurements in participants with and without diabetes. Data from population-based studies such as the UK Biobank cohort, offer the opportunity to evaluate changes in RNFL and GC-IPL thickness in a large number of participants with and without diabetes. Additionally, these data allow us to evaluate the impact of systemic and ocular conditions on retinal thickness measurements providing important information regarding factors to consider when using OCT-based thickness changes in the evaluation of DRN.

We aimed to identify structural changes to the neuroretina from diabetes, by comparing total retinal thickness and retinal thickness of the neural unit i.e., soma + dendrite, identified as ganglion cell-inner plexiform layer (GC-IPL) and axon as identified by retinal nerve fiber layer (RNFL) on macula centered SD-OCT scans of participants with and without DM in a large population-based cohort. We also evaluated the impact of A1c level on retinal thickness measurements, and factors associated with retinal thickness measurements in participants with diabetes.

Methods

Study population

The UK Biobank is a large population-based cohort of over 500,000 participants [15]. Detailed study protocols are available online (http://www.ukbiobank.ac.uk/resources/ and http://biobank.ctsu.ox.ac.uk/crystal/docs.cgi). Briefly, from 2006–10, all residents of the United Kingdom, between the ages of 40 to 69 years, who were registered with the National Health Service were invited to participate. The study was conducted in accordance with the principles of the declaration of Helsinki and approved by the ethics committee and institutional regulatory boards, and all participants gave informed consent. Participants initially completed a baseline questionnaire related to demographics, medical history, medication use, as well as lifestyle. This was followed by a verbal interview by a nurse, physical exam, blood draw and additional measures [15]. Ophthalmic assessment was added in 2009 for selected assessment centers. This included measurement of visual acuity, intraocular pressure, refractive error, fundus photos and SD-OCT scans. All participants who had SD-OCT measures were considered for inclusion in the study.

Retinal imaging in the UK Biobank

Macular SD-OCT scans plus a one-field 45-degree fundus photograph were acquired using a Topcon camera (Topcon 3D OCT-1000 Mark II; Topcon GB). The SD-OCT image acquisition protocol is described in detail in a previous publication [16]. Participants who had an eye infection or any eye surgery within the last 4 weeks did not undergo any ocular measurements. For this study we used Iowa Reference Algorithm software version 3.8 to determine 1) total retinal thickness 2) the thickness of the mRNFL and 3) GC-IPL. This software has been previously demonstrated to give accurate and reliable thickness measurements [1719]. Original image files were downloaded from the UKBB servers. Two automated measures of image quality were used: 1) undefined region and 2) surface cost. Undefined region quantified the percentage area of a scan that had missing data or insufficient signal. Images with undefined region > 0% were excluded from analysis. The surface cost is calculated using edge-based costs of the dark to bright and bright to dark transition of the retinal sub-layers. This is an inverted Gaussian-smoothed gradient magnitude of the OCT voxel intensities of the retinal sub-layers. A lower value corresponds to more reliable segmentation [20]. We randomly evaluated 100 OCT scans with a range of surface cost values and determined that scans with surface cost values > 62, 000 had unreliable segmentation. Subsequently all scans with surface cost values > 62, 000 were excluded. For the remaining OCT images, any scans with thickness measurements 2 or more standard deviation above or below the mean were manually reviewed. Scans with incorrect algorithm, poor quality, poor signal or presence of pathology such as epiretinal membrane, macular edema, drusen, intra-retinal or sub-retinal fluid, pigment epithelial detachment, decentration that could impact thickness measurements were excluded.

A previously validated procedure was used to determine DM status based on self-reported DM diagnosis, use of DM medications and presence of DM complications [21]. We also used HbA1c ≥ 6.5% as a criterion for identifying DM. Two graders (NM and JG) independently evaluated the fundus photos for quality and graded DR based on the International Clinical Diabetic Retinopathy Scale as no/mild DR (no changes or microaneurysms only) or more than mild DR [22]. Final DR grade was determined by the higher of the two grades. We defined glaucoma as self-reported diagnosis of glaucoma, taking glaucoma medications, or having had glaucoma surgery. We defined cognitive impairment as a failing score on 2 or more tests as has been described in prior publication [23].

Exclusions

Participants were excluded if they withdrew consent, were identified to have pre-diabetes or gestational diabetes; if they had poor quality OCT scans as described above, eyes affected by injury or trauma resulting in loss of vision; eyes affected by retinal pathology as described above which could distort the retinal layers; those with self-reported history of retinal surgery or macular degeneration; extreme IOP values of less than 5mmHg or more than 60mmHg.

Statistical analysis

Patient and eye characteristics were summarized using mean with standard deviation, median with minimum and maximum values, or frequency with percentage. Summary statistics were stratified by diabetes. Patient summary statistics were compared using the Wilcoxon rank sum, Fisher’s exact or Chi-square test. Eye summary statistics were compared using generalized estimating equations (GEE) with a logit link and exchangeable correlation. Independent general linear mixed models (GLMM) with a random intercept for patient and eye were used to test the association between patient/eye characteristics and SD-OCT thickness measurements (RNFL, GC-IPL, and total retinal thickness). Factors found to be significantly associated with DM or thickness, or those that were clinically significant were included in a multiple GLMM. When correlation between predictor variables was > 0.5, only one was included in the GLMM to avoid multicollinearity. Separate GLMMs were created using DM and A1c as predictor variables for thickness. Sub-analyses were done 1) including only patients with no/mild DR in the DM group, and 2) to evaluate factors associated with GC-IPL thickness specifically for DM patients. In addition to the co-variates used in other models, the analysis with only DM patients tested whether thickness differed for patients who had an A1c > = 6.5% or who were on metformin, whether they had DR and the duration of DM. P < 0.05 was considered statistically significant.

Results

74,422 participants (131,555 eyes) met all inclusion and exclusion criteria described above and were included in the analysis. 10,024 participants (35,096 eyes) were excluded from the analysis after quality control and other exclusion criteria described above. Comparison of summary statistics between those included and excluded in the analyses are presented in S1A and S1B Tables and a detailed list of the number of patients/eyes excluded from analysis is shown in S2 Table.

Of the 74,422 participants included in the study: 69,985 did not have DM and 4,437 had DM. Participants with DM were older (median 62.3 years vs 59 years, p<0.001), more likely to be men (63.0% vs 45.2%, p<0.001), more likely to be on cholesterol lowering (45.1% vs 9.5%, p<0.001) or blood pressure lowering medications (39.2% vs 9.9%, p<0.001), and more likely to have impairment on cognitive testing (29.6% vs 19.3%, p<0.001) compared to those without DM (See Tables 1and 2). Participants with DM had on average lower total retinal thickness (4.57 μm lower, 95% CI: -5.00, -4.14; p<0.001), GC-IPL thickness (1.73 μm lower, 95% CI: -1.86, -1.59; p<0.001) and mRNFL thickness (0.68 μm lower, 95% CI: -0.81, -0.54; p<0.001) compared to those without DM. A full list of the associations with total retinal thickness, mRNFL and GC-IPL thickness measures are shown in S3 Table. Systolic blood pressure and diastolic blood pressure were correlated and use of cholesterol lowering and blood-pressure lowering medications were correlated. To avoid collinearity, we chose systolic blood pressure and cholesterol lowering medication for the multivariable GLMM.

thumbnail
Table 1. Patient summary statistics stratified by diabetes (all included patients).

https://doi.org/10.1371/journal.pone.0257836.t001

thumbnail
Table 2. Eye summary statistics stratified by diabetes (all included eyes).

https://doi.org/10.1371/journal.pone.0257836.t002

After adjusting for co-variates in the multivariable GLMM (See Table 3), total retinal thickness was 1.99 um thinner (95% CI: -2.47–1.50, p<0.001) among participants with diabetes compared to those without DM, and GC-IPL was 1.02 um lower (95% CI: -1.18, -0.87; p<0.001) among those with diabetes compared to those without. The difference in mRNFL thickness was no longer significantly different (p = 0.369). Fig 1A–1C show the adjusted differences in average retinal thickness measurements after controlling for co-variates among participants with and without DM for total retinal thickness (Fig 1A), mRNFL (Fig 1B) and GC-IPL thickness (Fig 1C). On average, total retinal thickness decreased 0.40 um (95% CI:-0.61, -0.20; p<0.001), mRNFL thickness increased 0.20 um (95% CI: 0.14, 0.27; p<0.001) and GC-IPL decreased 0.26 um (95% CI: -0.33, -0.20; p<0.001) per unit increase in A1c, after adjusting for co-variates. See Fig 2A–2C.

thumbnail
Fig 1. Difference in retinal thickness measurements between participants with and without diabetes.

Fig 1 shows the difference in average thickness measurements (95% confidence interval) between participants with and without diabetes (DM) after adjusting for co-variates: a) total retinal thickness, b) retinal nerve fiber layer (RNFL) thickness, c) ganglion cell inner plexiform layer (GC-IPL) thickness.

https://doi.org/10.1371/journal.pone.0257836.g001

thumbnail
Fig 2. Change in retinal thickness measurements with change in glycated hemoglobin levels.

Fig 2 shows the change in a) average total retinal thickness, b) average retinal nerve fiber layer (RNFL) thickness, c) average ganglion cell-inner plexiform (GC-IPL) layer thickness after adjusting for co-variates.

https://doi.org/10.1371/journal.pone.0257836.g002

thumbnail
Table 3. Multivariable general linear mixed models for retinal thickness measurements.

https://doi.org/10.1371/journal.pone.0257836.t003

Factors associated with GC-IPL thickness among participants with DM

We further evaluated the factors associated with GC-IPL thickness specifically among participants with DM. Average GC-IPL thickness was not significantly different between patients who did and did not have A1c> = 6.5% (p = 0.120), and was not associated with duration of DM (p = 0.172). Average GC-IPL was 1.24 um thinner among those with more than mild DR compared to those with no/mild DR (95% CI: -1.89, -0.59; p<0.001). See Table 4 for multivariable model of the factors associated with GC-IPL thickness among participants with DM. This sub-analysis included 1,959 participants (3,239 eyes) with no/mild DR and 146 participants (190 eyes) with more than mild DR.

thumbnail
Table 4. Multivariable general linear mixed model for GC-IPL thickness among participants with diabetes.

https://doi.org/10.1371/journal.pone.0257836.t004

Sub-analysis: No diabetes vs diabetes (with no/mild DR). From 7,687 eyes with diabetes, 5,433 eyes were graded as no/mild and 278 as more than mild. 1,976 eyes were excluded from this analysis because either a) image was missing (n = 94) or b) image was graded as poor quality (n = 1,882).

Average total retinal thickness was lower among those with DM (and no/mild DR) compared to those without DM (4.25 um lower, 95% CI: -4.73, -3.76; p<0.001). The average RNFL thickness was lower among participants with DM (and no/mild DR) compared to those without DM (0.60 um lower, 95% CI: -0.76, -0.44; p<0.001). The average GCL-IPL was thinner among participants with DM and no/mild DR compared to those without DM (1.41 um lower, 95% CI: -1.56, -1.25; p<0.001).

After adjusting for the co-variates mentioned above, the total retinal thicknesses of those with DM (and no/mild DR) was 1.89 um lower compared to those without DM (95% CI: -2.43, -1.34; p<0.001); the average RNFL thickness did not differ among those with DM (and no/mild DR) compared to those without DM (p = 0.417); average GC-IPL thickness was 0.80 microns lower among participants with DM (and no/mild DR) compared to those without DM (95% CI: -0.98, -0.62; p<0.001). S1A–S1C Fig show the differences in average retinal thickness measurements after controlling for co-variates among patients with DM (no/mild DR) and without DM.

After controlling for co-variates, the total retinal thickness decreased 0.39 um (95% CI: -0.62, -0.17; p<0.001) per unit increase in A1c; average RNFL thickness increased 0.21 um (95% CI: 0.14, 0.28; p<0.001) per unit increase in A1c and GC-IPL thickness decreased 0.22 um (95% CI: -0.30, -0.15; p<0.001) per unit increase in A1c.

Discussion

In this study we sought to address the question regarding the impact of DM on total retinal thickness and thickness of the neuro-retina (mRNFL and GC-IPL), as measured using macular OCT scans, in a large population-based cohort. Our findings show that GC-IPL is thinner among participants with DM compared to non-DM. Additionally, this study shows that thinning persisted after controlling for the many confounding factors associated with retinal thickness measurements, and after excluding participants with DM who had more than mild DR. These findings are consistent with our prior findings from 45 participants with Type 1 DM where we showed progressive GC-IPL thinning over time and those of Aschauer et al showing GC-IPL thinning but no mRNFL thinning over time among patients with Type 2 DM [6, 24]. A recent meta-analysis on this topic also showed a stronger association with GC-IPL thickness compared to mRNFL thickness between DM patients with no DR and controls [25]. Our findings, together with prior work from others, suggest that GC-IPL thickness may be a good marker for evaluation of DRN and animal studies suggest that thinning of GC-IPL may represent loss of ganglion cells [26]. Direct visualization of ganglion cells in-vivo would help further elucidate this.

Our large population-based study also allows us to study ocular and systemic factors associated with GC-IPL thickness among those with and without diabetes. This information is helpful when using GC-IPL clinically as a measure to identify DRN. Factors associated with GC-IPL thickness among patients with DM included age, ethnicity, glaucoma, refractive error, DR grade and logMAR visual acuity. Independent of DR grade, thinner GC-IPL was associated with worse visual acuity, suggesting that DRN, even in the early stages may have an impact on vision. BMI was significant with a small coefficient of -0.07 suggesting that the impact on thickness is likely small. Interestingly gender was not significantly associated with GC-IPL thickness among participants with DM. Khawaja et al reported a lower GC-IPL thickness among men compared to women, and attributed this to the higher rates of glaucoma reported among men [27]. This finding would need to be explored further in additional studies as we know from prior literature that men tend to have higher total retinal thickness measurements compared to women [28]. Consistent with the prospective study by Sohn et al among Type 1 DM patients, and Lim et al [14], A1c of 6.5% or more was not associated with GC-IPL thickness among those with DM [6]. Unlike this study, duration of DM was significantly associated with GC-IPL thickness in the study by Sohn et al. The difference is likely due to the difficulty in measuring duration of disease in patients with Type 2 DM, who may have sub-clinical disease before clinical diagnosis is made. In the study by Sohn et al, all patients had Type 1 DM, and the duration of diabetes was exactly known [6]. In this study, participants, predominantly had type 2 DM, making it difficult to determine the exact date of diagnosis. Aschauer et al, in their study evaluating changes in retinal thickness over time among patients with Type 2 DM, also did not find an association with duration, likely due to the same reason [24]. Further studies are needed to identify patients who are at the highest risk of progressive DRN as treating them with neuroprotective agents, such as brimonidine and somatostatin may be beneficial, as shown in the EUROCONDOR trial [29].

The changes in mRNFL thickness among those with and without DM seemed to be more complex. mRNFL thickness was on average lower in participants with DM compared to non-DM. When we adjusted for confounders, this difference was no longer significant. When we did a sub-analysis only including DM participants with no/mild DR, mRNFL thickness did not differ between those with and without DM. When A1c was used as a predictor, mRNFL thickness increased with increasing A1c levels. There may be several possible explanations for this: 1) neurodegeneration is a process which involves axonal swelling followed by atrophy, our findings suggest that in early stages of DR, mRNFL thickness may be higher than normal; 2) DM affects GC-IPL differentially than mRNFL, mRNFL also includes Muller cells and activated astrocytes/glia, which when activated can lead to thickening of the mRNFL;3) sub-clinical macular edema may also contribute to thicker mRNFL. Baseline data of participants enrolled in the EUROCONDOR trial showed that about a third of patients with Type 2 DM had early microvascular changes but no multifocal ERG changes, suggesting that while some patients have DRN, others do not [30]. Pooling the findings of those with DRN and those without DRN may be minimizing the overall differences in layer thickness between groups. Irrespective of the exact reasons for this finding, our results as well as that of others [24] suggest that GC-IPL thickness may be a more robust marker of early DRN changes compared to mRNFL.

To the best of our knowledge, our study, is one of the largest study to evaluate retinal thickness measures in patients with and without DM, while controlling for multiple confounders, and to report factors associated with GC-IPL thickness among patients with DM.

The strength of our study includes the large numbers, detailed data available on cofounding factors, as well as the rigorous criteria used to define DM. Accuracy of self-reported DM in population-based studies can be limited [3133]. In order to make our groups as “clean” as possible we also included hemoglobin A1c levels in our determination of DM and excluded patients with pre-diabetes or gestational diabetes.

There are certain limitations to our study, which was retrospective by nature. We were limited by the imaging modalities available in the UKBB dataset. We could not assess peripapillary RNFL thickness and can only comment on the changes in macular RNFL thickness, but as macular centered OCT scans offer the additional advantage of evaluating diabetic macular edema, information regarding thickness measurements from macular scans may be more relevant for clinical practice. Diagnosis of glaucoma was based on self-report and we could have missed some participants with undiagnosed glaucoma or some who incorrectly recalled their diagnosis of glaucoma. However, as the criteria for identifying glaucoma were the same across groups, it is likely that those misclassified, were equally distributed between the two groups. Although, we were limited by one-field fundus photos in our dataset, prior work has shown high agreement between single-field fundus photos and the standard seven field stereoscopic photographs, with a sensitivity and specificity of 78% and 86% respectively, for DR detection [34]. Despite these limitations, our results show, that in a large population-based study GC-IPL thickness is lower among patients with DM compared to those without DM, and this difference persists when considering only those with no/mild DR and after controlling for confounding factors. Our findings also suggest the mRNFL thickness may be more impacted by confounding variables.

Our report adds to a growing body of evidence that DRN is an important component of diabetic eye disease, and that GC-IPL thickness may be a good measure to quantify DRN. Our study further reports factors associated with GC-IPL thickness among those with DM. These include age, ethnicity, refractive error, DR grade, BMI and visual acuity, suggesting that these factors must be considered when evaluating GC-IPL thickness in patients with DM. Additional studies are needed to determine which patients with DM are most likely to develop DRN and to improve our understanding of ganglion cell health by directly visualizing retinal ganglion cells in vivo using newer imaging techniques such as adaptive optics [35]. Our study also showed that that thinner GC-IPL was associated with worse vision independent of DR grade, additional studies are needed to fully evaluate the impact of DRN on functional measures of vision. Finally, as the eye offers unique insights into future risks of disease elsewhere, links between DRN and diabetic peripheral neuropathy and cognitive impairment warrants further study.

Supporting information

S1 Table.

a. Patient summary statistics stratified by included vs excluded. b. Eye summary statistics stratified by included vs excluded.

https://doi.org/10.1371/journal.pone.0257836.s001

(DOCX)

S2 Table. Detailed list of exclusions and reasons for exclusions.

https://doi.org/10.1371/journal.pone.0257836.s002

(DOCX)

S3 Table. Factors association with retinal thickness measurements (Univariate GLMMs).

https://doi.org/10.1371/journal.pone.0257836.s003

(DOCX)

S1 Fig. Difference in retinal thickness measurements between participants with and without diabetes (with no/mild diabetic retinopathy).

S1 Fig shows the difference in average thickness measurements (95% confidence interval) between participants with and without diabetes (DM) after adjusting for co-variates and excluding those with more than mild diabetic retinopathy: a) total retinal thickness, b) retinal nerve fiber layer (RNFL) thickness, c) ganglion cell inner plexiform layer (GC-IPL) thickness.

https://doi.org/10.1371/journal.pone.0257836.s004

(JPG)

Acknowledgments

This analysis was conducted using data from UK Biobank, under project number 41561 under terms of the access MTA.

References

  1. 1. Yau J.W.Y., et al., Global Prevalence and Major Risk Factors of Diabetic Retinopathy. Diabetes Care, 2012: p. DC_111909. pmid:22301125
  2. 2. Lynch S.K. and Abramoff M.D., Diabetic retinopathy is a neurodegenerative disorder. Vision Res, 2017. 139: p. 101–107. pmid:28408138
  3. 3. Bloodworth J.M. Jr., Diabetic retinopathy. Diabetes, 1962. 11: p. 1–22. pmid:13870120
  4. 4. Abramoff M.D., et al., Approach for a clinically useful comprehensive classification of vascular and neural aspects of diabetic retinal disease. Investigative ophthalmology & visual science, 2018. 59(1): p. 519–527. pmid:29372250
  5. 5. Sun J.K., et al., Updating the Staging System for Diabetic Retinal Disease. Ophthalmology, 2020. pmid:33218709
  6. 6. Sohn E.H., et al., Retinal neurodegeneration may precede microvascular changes characteristic of diabetic retinopathy in diabetes mellitus. Proc Natl Acad Sci U S A, 2016. 113(19): p. E2655–64. pmid:27114552
  7. 7. Rauscher F.M., et al., Myopia affects retinal nerve fiber layer measurements as determined by optical coherence tomography. Journal of glaucoma, 2009. 18(7): p. 501. pmid:19745664
  8. 8. Budenz D.L., et al., Determinants of normal retinal nerve fiber layer thickness measured by Stratus OCT. Ophthalmology, 2007. 114(6): p. 1046–1052. pmid:17210181
  9. 9. Oktem E.O., et al., The relationship between the degree of cognitive impairment and retinal nerve fiber layer thickness. Neurological Sciences, 2015. 36(7): p. 1141–1146. pmid:25575807
  10. 10. Bowd C., et al., The retinal nerve fiber layer thickness in ocular hypertensive, normal, and glaucomatous eyes with optical coherence tomography. Archives of ophthalmology, 2000. 118(1): p. 22–26. pmid:10636409
  11. 11. Zafar S., et al., Retinal neurodegeneration as an early manifestation of diabetic eye disease and potential neuroprotective therapies. Current diabetes reports, 2019. 19(4): p. 17. pmid:30806815
  12. 12. Schroder K., et al., German Diabetes Study—Baseline data of retinal layer thickness measured by SD-OCT in early diabetes mellitus. Acta Ophthalmol, 2018. pmid:30238609
  13. 13. Lim H.B., et al., Longitudinal Changes in the Peripapillary Retinal Nerve Fiber Layer Thickness of Patients With Type 2 Diabetes. JAMA Ophthalmology, 2019. 137(10): p. 1125–1132. pmid:31343674
  14. 14. Lim H.B., et al., Ganglion Cell–Inner Plexiform Layer Damage in Diabetic Patients: 3-Year Prospective, Longitudinal, Observational Study. Scientific reports, 2020. 10(1): p. 1–9. pmid:31913322
  15. 15. Bycroft C., et al., The UK Biobank resource with deep phenotyping and genomic data. Nature, 2018. 562(7726): p. 203–209. pmid:30305743
  16. 16. Patel P.J., et al., Spectral-Domain Optical Coherence Tomography Imaging in 67 321 Adults: Associations with Macular Thickness in the UK Biobank Study. Ophthalmology, 2016. 123(4): p. 829–40. pmid:26746598
  17. 17. Sohn E.H., et al., Reproducibility of diabetic macular edema estimates from SD-OCT is affected by the choice of image analysis algorithm. Investigative ophthalmology & visual science, 2013. 54(6): p. 4184–4188. pmid:23696607
  18. 18. Terry L., et al., Automated retinal layer segmentation using spectral domain optical coherence tomography: evaluation of inter-session repeatability and agreement between devices. PloS one, 2016. 11(9): p. e0162001. pmid:27588683
  19. 19. Lee K., et al., Distribution of damage to the entire retinal ganglion cell pathway: quantified using spectral-domain optical coherence tomography analysis in patients with glaucoma. Archives of ophthalmology, 2012. 130(9): p. 1118–1126. pmid:22965586
  20. 20. Lee K., et al., Automated segmentability index for layer segmentation of macular SD-OCT images. Translational vision science & technology, 2016. 5(2): p. 14–14. pmid:27066311
  21. 21. Eastwood S.V., et al., Algorithms for the capture and adjudication of prevalent and incident diabetes in UK Biobank. PloS one, 2016. 11(9): p. e0162388. pmid:27631769
  22. 22. Wilkinson C., et al., Proposed international clinical diabetic retinopathy and diabetic macular edema disease severity scales. Ophthalmology, 2003. 110(9): p. 1677–1682. pmid:13129861
  23. 23. Ko F., et al., Association of Retinal Nerve Fiber Layer Thinning With Current and Future Cognitive Decline: A Study Using Optical Coherence Tomography. JAMA neurology, 2018. pmid:29946685
  24. 24. Aschauer J., et al., Longitudinal analysis of microvascular perfusion and neurodegenerative changes in early type 2 diabetic retinal disease. British Journal of Ophthalmology, 2020. pmid:33293271
  25. 25. Tang Z., et al., Assessment of retinal neurodegeneration with spectral-domain optical coherence tomography: a systematic review and meta-analysis. Eye (Lond), 2020. pmid:32581390
  26. 26. Hedberg-Buenz A., et al., Biological Correlations and Confounding Variables for Quantification of Retinal Ganglion Cells Based on Optical Coherence Tomography using Diversity Outbred Mice. bioRxiv, 2020: p. 2020.12.23.423848.
  27. 27. Khawaja A.P., et al., Comparison of associations with different macular inner retinal thickness parameters in a large cohort: the UK Biobank. Ophthalmology, 2020. 127(1): p. 62–71. pmid:31585827
  28. 28. Chalam K.V., et al., Retinal thickness in people with diabetes and minimal or no diabetic retinopathy: Heidelberg Spectralis optical coherence tomography. Investigative ophthalmology & visual science, 2012. 53(13): p. 8154–8161. pmid:23132803
  29. 29. Simó R., et al., Effects of topically administered neuroprotective drugs in early stages of diabetic retinopathy: results of the EUROCONDOR clinical trial. Diabetes, 2019. 68(2): p. 457–463. pmid:30389750
  30. 30. Santos A.R., et al., Functional and structural findings of neurodegeneration in early stages of diabetic retinopathy: cross-sectional analyses of baseline data of the EUROCONDOR project. Diabetes, 2017. 66(9): p. 2503–2510. pmid:28663190
  31. 31. Goto A., et al., Validity of diabetes self-reports in the Saku diabetes study. J Epidemiol, 2013. 23(4): p. 295–300. pmid:23774288
  32. 32. Jackson J.M., et al., Validity of diabetes self-reports in the Women’s Health Initiative. Menopause, 2014. 21(8): p. 861–8. pmid:24496083
  33. 33. Sluijs I., et al., Ascertainment and verification of diabetes in the EPIC-NL study. Neth J Med, 2010. 68(1): p. 333–9. pmid:20739736
  34. 34. Lin D.Y., et al., The sensitivity and specificity of single-field nonmydriatic monochromatic digital fundus photography with remote image interpretation for diabetic retinopathy screening: a comparison with ophthalmoscopy and standardized mydriatic color photography11InternetAdvance publication at ajo.com. April 12, 2002. American Journal of Ophthalmology, 2002. 134(2): p. 204–213. pmid:12140027
  35. 35. Liu W.W. and Margeta M.A., Imaging retinal ganglion cell death and dysfunction in glaucoma. International ophthalmology clinics, 2019. 59(4): p. 41–54. pmid:31569133