Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
Scientific Reports volume 12, Article number: 17750 (2022)
Metrics details
Retinopathy and neuropathy in type 2 diabetes are preceded by retinal nerve fibre layer (RNFL) thinning, an index of neurodegeneration. We investigated whether glucose metabolism status (GMS), measures of glycaemia, and daily glucose variability (GV) are associated with RNFL thickness over the entire range of glucose tolerance. We used cross-sectional data from The Maastricht Study (up to 5455 participants, 48.9% men, mean age 59.5 years and 22.7% with type 2 diabetes) to investigate the associations of GMS, measures of glycaemia (fasting plasma glucose [FPG], 2-h post-load glucose [2-h PG], HbA1c, advanced glycation endproducts [AGEs] assessed as skin autofluorescence [SAF]) and indices of daily GV (incremental glucose peak [IGP] and continuous glucose monitoring [CGM]-assessed standard deviation [SD]) with mean RNFL thickness. We used linear regression analyses and, for GMS, P for trend analyses. We adjusted associations for demographic, cardiovascular risk and lifestyle factors, and, only for measures of GV, for indices of mean glycaemia. After full adjustment, type 2 diabetes and prediabetes (versus normal glucose metabolism) were associated with lower RNFL thickness (standardized beta [95% CI], respectively − 0.16 [− 0.25; − 0.08]; − 0.05 [− 0.13; 0.03]; Ptrend = 0.001). Greater FPG, 2-h PG, HbA1c, SAF, IGP, but not CGM-assessed SD, were also associated with lower RNFL thickness (per SD, respectively − 0.05 [− 0.08; − 0.01]; − 0.06 [− 0.09; − 0.02]; − 0.05 [− 0.08; − 0.02]; − 0.04 [− 0.07; − 0.01]; − 0.06 [− 0.12; − 0.01]; and − 0.07 [− 0.21; 0.07]). In this population-based study, a more adverse GMS and, over the entire range of glucose tolerance, greater glycaemia and daily GV were associated with lower RNFL thickness. Hence, early identification of individuals with hyperglycaemia, early glucose-lowering treatment, and early monitoring of daily GV may contribute to the prevention of RNFL thinning, an index of neurodegeneration and precursor of retinopathy and neuropathy.
Retinopathy and neuropathy, both hallmark microvascular complications of type 2 diabetes, are preceded by subtle neurodegenerative changes. Such changes include retinal nerve fibre layer (RNFL) thinning, which can be non-invasively assessed1,2,3. RNFL thinning reflects a gradual loss of retinal ganglion cell axons, which transmit visual information from the retina to the brain1,2. Mechanistically, elevated glucose concentrations are thought to be toxic for retinal ganglion cells as well as for retinal endothelial and glia cells, which contribute to local metabolic regulation4,5. Thus, hyperglycaemia can lead to retinal ganglion cell apoptosis and loss of retinal ganglion cell axons both directly and through impairment of endothelial and glia cell function6.
As postulated in the “ticking clock hypothesis”, hyperglycaemia-mediated damage is thought to be a continuous (i.e., linear) process that already starts before the onset of type 2 diabetes7,8. Indeed, our group observed linear associations of more adverse glucose metabolism status (GMS) and higher glycaemia (estimated by various measures) with lower heart rate variability9, more structural brain abnormalities10,11, and worse peripheral nerve function12, all of which are measures of neurodegeneration.
The current literature on the associations of hyperglycaemia with RNFL thickness has some important limitations. First, in previous population-based studies, several important confounders were not included, such as age13,14,15, sex14,15, socioeconomic status13,14,15,16,17,18,19,20,21,22,23, cardiovascular risk factors13,14,15,19,21,22,23,24, and lifestyle factors (e.g., alcohol use13,16,17,18,19,21, diet13,14,15,16,17,18,19,20,21,22,23,24, and physical activity)13,14,15,16,17,18,19,21,22,23,24. Second, no population-based studies have yet investigated whether advanced glycation endproducts (AGEs) or daily glucose variability are associated with RNFL thickness. It is important to establish to what extent hyperglycaemia (including AGEs) and daily glucose variability contribute to RNFL thinning over the entire range of glucose tolerance because early diagnosis and treatment of hyperglycaemia as well as novel strategies to monitor glycaemic exposure may contribute to the early prevention of hyperglycaemia-mediated RNFL thinning5,25,26.
In view of above, we investigated, using a large, well-characterized population-based cohort study, whether more adverse GMS, greater glycaemia, and greater daily glucose variability are associated with lower RNFL thickness.
We used data from The Maastricht Study, a prospectively designed, population-based observational cohort study. The rationale and methodology have been described previously27. In brief, the study focuses on the aetiology, pathophysiology, complications and comorbidities of type 2 diabetes mellitus and is characterized by an extensive phenotyping approach. Eligible for participation were all individuals aged between 40 and 75 years and living in the southern part of the Netherlands. Participants were recruited through mass media campaigns and from the municipal registries and the regional Diabetes Patient Registry via mailings. Recruitment was stratified according to known type 2 diabetes status, with an oversampling of individuals with type 2 diabetes, for reasons of efficiency27. The present report includes cross-sectional data of 8005 participants who completed the baseline survey between November 2010 and September 2018. The examinations of each participant were performed within a time window of three months. From 19 September 2016 until 13 September 2018, participants were invited to also undergo CGM26. During this period, a selected group of recently included participants was invited to return for CGM (‘catch-up visit’). For these participants only there was a median time interval (‘visit interval’) of 2.1 years between CGM and all other measurements (more details are provided in the Supplemental Material).
After an overnight fast, participants underwent a standardized seven-point oral glucose tolerance test (OGTT) as part of which venous samples were collected at 15, 30, 45, 60, 90, and 120 min post ingestion of a 75 g glucose drink. All participants underwent an OGTT except for the participants who used insulin or had a fasting plasma glucose (FPG) concentration above 11.0 mmol/L. Based on FPG and 2-h post-load glucose, GMS was determined as normal glucose metabolism (NGM), prediabetes (impaired fasting glucose, impaired glucose tolerance, or both), type 2 diabetes, or other types of diabetes (including type 1 diabetes) in accordance with the World Health Organization 2006 criteria28.
FPG (mmol/L) and haemoglobin A1c (HbA1c; mmol/mol; %) were determined in venous plasma samples collected after an overnight fast. Two-hour post-load glucose (mmol/L) was determined in venous plasma collected at 120 min post glucose drink ingestion. AGEs were assessed with the AGE Reader (DiagnOptics Technologies BV, Groningen, the Netherlands). In brief, the AGE Reader is a desktop device that uses the characteristic fluorescent properties of certain AGEs to quantify their accumulation in the skin as skin autofluorescence (SAF; arbitrary units [AU])29. The AGE Reader illuminates a skin surface of 4 cm2, shielded from other light, and uses the ratio of the reflection of fluorescent light (wavelength 420 to 600 nm) to non-fluorescent light (300–420 nm) to calculate SAF (more details are provided in the Supplemental Materials).
The incremental glucose peak (IGP; mmol/L), a recently validated OGTT-based index of daily glucose variability, was calculated by subtracting FPG from the maximum glucose peak value measured during a complete seven-point OGTT30. At the time of analysis, data on IGP were available in a subset of participants (n = 2407). We used CGM (iPro2 and Enlite Glucose Sensor; Medtronic, Tolochenaz, Switzerland) to assess daily glucose variability during a 1-week period, which was expressed as standard deviation (mmol/L) of mean sensor glucose (mmol/L)26. CGM-assessed data were available in a subset of participants (n = 622). More details on the assessment of glycemic indices with OGTT and CGM are provided in the Supplemental Methods.
We assessed peripapillary RNFL thickness (μm) in both eyes using optical coherence tomography (OCT) (Spectralis unit and Eye Explorer version 5.7.5.0 software; Heidelberg Engineering, Heidelberg, Germany; 3.45-mm-diameter circle scan, manually centred on optic nerve head, 12°, 768 voxels, 100 automatic real-time tracking). Intra- and interindividual reliability, expressed as intraclass correlation coefficients, are 0.97 and 0.96 respectively31. At least 15 min before the examination, pupils were dilated with topical 0.5% tropicamide and 2.5% phenylephrine. Experienced graders masked to clinical information on the participants reviewed the OCT scans and graded their quality. OCT images were excluded if one of the following criteria was present: scan error (i.e., incomplete scan, poor centring of the circular scan on the optic nerve head, RNFL layer incorrectly defined, or technical problem with the OCT device) or poor imaging quality (signal-to-noise ratio < 15 dB)23. If data from both eyes were available (n = 2796 participants) we averaged RNFL thickness in order to reduce measurement error. If data from only one eye were available (n = 2755 participants), we used the RNFL thickness of that eye in the analyses. More details, including on quality criteria, are shown in the Supplemental Methods.
As described previously27, we assessed educational level (low, intermediate, high), socio-economic status (income level and occupational status) 32, smoking status (never, former, current), alcohol use (none, low, high), history of cardiovascular disease, and duration of diabetes by questionnaire; assessed dietary habits (“dietary intake”) with the Dutch Healthy Diet index sum score, a measure of adherence to the Dutch dietary guidelines 201533, based on a validated food frequency questionnaire34; assessed lipid-modifying, antihypertensive, intraocular pressure-lowering, and glucose-lowering medication use as part of a medication interview; assessed weight, height, and waist circumference during a physical examination; calculated body mass index (BMI) based on body weight and height; measured office and 24-h ambulatory blood pressure (BP); measured total daily physical activity (hours/day) with an accelerometer35; measured lipid profile and plasma biomarkers of low-grade inflammation36 (i.e., high-sensitive C-reactive protein, serum amyloid A, interleukin-6, interleukin-8 and tumour necrosis factor alpha) in fasting venous blood samples; measured urinary albumin excretion in two 24-h urine collections; calculated the estimated glomerular filtration rate (eGFR) based on serum creatinine only, since cystatin C was not presently available in all study participants; assessed presence of retinopathy in both eyes via fundus photography (more details in the Supplemental Material); and used an automated refractor and noncontact tonometer (Tonoref II; Nidek, Gamagordi, Japan) to assess spherical equivalent and intraocular pressure in both eyes. Glaucoma was defined as use of intraocular pressure-lowering medication, intraocular pressure higher than 21 mm Hg in any eye (91.3% of all participants had data on intraocular pressure available for at least 1 eye), or both. Spherical equivalent was defined as the mean spherical equivalent of both eyes (available for 91.1% of all participants) or as the spherical equivalent of the eye for which data were available.
We used multivariable linear regression analyses to investigate the associations of GMS (entered as dummy variables of prediabetes, type 2 diabetes, or other types of diabetes versus NGM) and standardized FPG, 2-h post-load glucose, HbA1c, SAF, IGP, and CGM-assessed standard deviation (determinants) with standardized mean RNFL thickness (outcome).
We performed P for trend analyses to test for linear trend with more adverse GMS. To test for trend, we entered GMS into the model as an ordinal variable (i.e., GMS was coded as NGM = 0, prediabetes = 1, type 2 diabetes = 2).In P for trend analyses we excluded participants with other types of diabetes because other types of diabetes (such as type 1 diabetes) do not constitute part of the spectrum of deterioration of GMS from NGM to prediabetes and type 2 diabetes. Then, we checked whether we could assume a linear trend by comparing the statistical variance explained by the statistical model in which GMS was entered as dummy variables to the statistical model in which GMS was entered as an ordinal variable. We used a likelihood ratio test to assess whether the amount of variance explained by both models differed statistically significantly. A P-value > 0.05 indicates that both models are not different and, thus, that a linear trend can be assumed. For all analyses under study, the P-value for the likelihood ratio test was > 0.05 (data not shown) and therefore a linear trend could be assumed.
Model 1 shows crude results. In model 2, we adjusted for age, sex, educational status (low, medium, high). We chose these variables because they are key potential confounders17. In model 3, we additionally adjusted for variables of which their status as confounder has been less firmly established (office systolic BP, use of antihypertensive medication [yes/no], waist circumference, total cholesterol / HDL cholesterol ratio, lipid-modifying medication [yes/no], smoking status [current, former, never], and alcohol consumption status [none, low, high])17. Then, and only for IGP and CGM-assessed standard deviation, we additionally adjusted for HbA1c or mean sensor glucose, respectively (model 3 + HbA1c; model 3 + mean sensor glucose), so that we could differentiate between daily glucose variability and average glycaemia, both of which are strongly related37. Additionally, and only for CGM-assessed indices, we entered ‘visit interval’ in model 1. The associations were expressed as standardized regression coefficient (stβ) and corresponding 95% confidence interval (95%CI). Collinearity diagnostics (i.e., tolerance < 0.10 and/or variance inflation factor > 10) were used to detect excessive multicollinearity between covariates.
We tested for interaction to assess whether associations differed by GMS (i.e., between individuals with type 2 diabetes, individuals with prediabetes, and individuals with NGM) or by sex. For interaction analyses with GMS, we excluded participants with other types of diabetes from the interaction analyses because the number of these participants was small.
To assess the robustness of our findings we performed several sensitivity analyses. In brief, we adjusted for potential confounders which we did not include in the main model (e.g. for lifestyle factors [dietary intake, physical activity] or ocular variables [spherical equivalent, intraocular pressure])38,39,40,41; we performed analyses with other indices of CGM-assessed daily glucose variability that were not included in the main manuscript (i.e. coefficient of variation [CV], time-in-range [TIR], time below range [TBR]; and time above range [TAR])42,43; and we repeated analyses after replacement of covariates which other covariates which reflect a similar underlying construct (e.g. we replaced waist circumference with BMI). A detailed overview of all sensitivity analyses is presented in the Supplemental Methods section.
All analyses were performed with Statistical Package for Social Sciences version 25.0 (IBM SPSS, IBM Corp, Armonk, NY, USA). For all analyses, a P-value < 0.05 was considered statistically significant.
The study has been approved by the institutional medical ethical committee (NL31329.068.10) and the Minister of Health, Welfare and Sports of the Netherlands (Permit 131088-105234-PG). All participants gave written informed consent27.
Figure 1 gives an overview of the study population selection. Participants in whom OCT data were missing or of insufficient quality were excluded first (n = 2454). Next, individuals with missing data on confounders were excluded (n = 96). The sample size of the final study populations depended on the availability of data on the main determinant (n = 5455 for GMS, n = 982 to 5454 for measures of glycaemia, and n = 622 to 2407 for indices of daily glucose variability).
RNFL study population selection. * Not mutually exclusive. OCT, optical coherence tomography; RNFL, retinal nerve fiber layer; HDL, high density lipoprotein; NGM, normal glucose metabolism; PreD, prediabetes; T2D, type 2 diabetes; FPG, fasting plasma glucose; HbA1c, hemoglobin A1C; SAF, skin autofluorescence; IGP, incremental glucose peak; CGM, continuous glucose monitoring.
Table 1 and Supplemental Table 1 show general participant characteristics according to tertiles of RNFL thickness. Overall, participants with a thinner RNFL were more often men, were older, and had a more adverse cardiovascular risk profile. General characteristics of participants included in the study were comparable to those of participants with missing data (Supplemental Table S2).
After full adjustment (model 3), a more adverse GMS was associated with lower RNFL thickness (standardized beta [95%CI], type 2 diabetes versus NGM − 0.16 [− 0.25; − 0.08]; prediabetes versus NGM − 0.05 [− 0.13; 0.03], P for trend = 0.001; Table 2 and Fig. 2).
Associations between glucose metabolism status, measures of glycaemia and indices of daily glucose variability with RNFL thickness (per SD). Standardized regression coefficient (stβ) represents the difference in RNFL thickness in SD for type 2 diabetes or prediabetes status versus NGM status, or per SD greater measure of glycaemia or index of daily glucose variability. In the GMS, fasting plasma glucose, 2-h post-load glucose, HbA1c, and SAF study populations, 1 SD corresponds with 10.8 μm for RNFL thickness, 1.5 mmol/L for fasting plasma glucose, 4.0 mmol/L for 2-h post-load glucose, 0.8% or 9.1 mmol/mol for HbA1c, and 0.5 AU for SAF. For incremental glucose peak, 1 SD corresponds with 11.1 μm for RNFL thickness and 2.9 mmol/L for incremental glucose peak. For CGM-assessed standard deviation, 1 SD corresponds with 10.8 μm for RNFL thickness and 0.58 mmol/L for CGM-assessed standard deviation. Variables entered in the models in addition to glucose metabolism status, measures of glycaemia, or indices of daily glucose variability: age, sex, and educational status (low, medium, high), office systolic blood pressure, total cholesterol to HDL cholesterol ratio, use of antihypertensive or lipid-modifying medication (yes/no), waist circumference, smoking status (current, ever, never), and alcohol consumption status (none, low, high). * Indicates statistically significant (P < 0.05). 1The associations of indices of daily glucose variability with RNFL thickness were additionally adjusted for HbA1c (IGP) or mean sensor glucose (CGM-assessed SD). 2Additionally, and only for CGM-assessed SD, we entered ‘visit interval’ in model 1. SD, standard deviation; CI, confidence interval; NGM, normal glucose metabolism; T2D, type 2 diabetes; FPG, fasting plasma glucose; HbA1c, hemoglobin A1c; SAF, skin autofluorescence; IGP, incremental glucose peak; CGM, continuous glucose monitoring.
After full adjustment (model 3), greater FPG, 2-h post-load glucose, HbA1c, and SAF were statistically significantly associated with lower RNFL thickness (per SD, respectively − 0.05 [− 0.08; − 0.01]; − 0.06 [− 0.09; − 0.02]; − 0.05 [− 0.08; − 0.02]; and − 0.04 [− 0.07; − 0.01]; Table 2 and Fig. 2).
After full adjustment (model 3), greater IGP was statistically significantly associated with lower RNFL thickness (per SD, − 0.06 [− 0.11; − 0.01]; Table 2 and Fig. 2). The association remained statistically significant after additional adjustment for HbA1c (per SD, − 0.06 [− 0.12; − 0.01]). After full adjustment (model 3), CGM-assessed standard deviation was also negatively associated with RNFL thickness, but not statistically significantly (per SD, − 0.08 [− 0.17; 0.01]). The association was similar after further adjustment for mean sensor glucose.
GMS and sex did not modify any of the associations under study. All P-values for interaction are shown in Supplemental Table S4.
Quantitatively similar results were observed in a range of sensitivity analyses and are presented in the Supplemental Results section.
The present population-based study had two main findings. First, a more adverse GMS, greater glycaemia (estimated from FPG, 2-h post-load glucose, HbA1c, and SAF) and greater daily glucose variability (estimated from IGP and CGM-assessed standard deviation) were all linearly and –except for CGM-assessed standard deviation– statistically significantly associated with lower RNFL thickness. Second, the associations between indices of daily glucose variability and RNFL thickness did not materially change after additional adjustment for measures of glycaemia.
Our findings are in line with and extend observations from most previous studies13,19,20,21,22,23,24. The present study is the first large population-based study to comprehensively report associations of GMS, measures of glycaemia, and indices of daily glucose variability with RNFL thickness, and also adjust for an extensive set of potential confounders. Additionally, the present study is the first to present associations of SAF, duration of diabetes, and indices of daily glucose variability with RNFL thickness.
Mechanistically, the linearity of the associations of GMS, measures of glycaemia, and indices of daily glucose variability with RNFL thickness may reflect the increasing loss of retinal ganglion cells due to both hyperglycaemia-induced neurotoxicity and impairment of functioning of retinal cells that contribute to metabolic regulation4,5. Such impairment of metabolic regulation can predispose retinal ganglion cells to ischemia5,44. Importantly, retinal ganglion cells are thought to be highly susceptible to ischemia, since they are highly active and have an energy demand that exceeds that of brain cells44.
These findings extend our previous work on the “ticking clock hypothesis”5,9,10,11,12, which postulates that hyperglycaemia-induced microvascular and neuronal deterioration is a continuous, gradual process that starts in prediabetes, progresses with the onset of type 2 diabetes, and continues during type 2 diabetes8. Indeed, we observed that the regression estimate for prediabetes was in between the estimate for type 2 diabetes and the reference category (i.e., NGM), and was directionally and numerically comparable to our previous findings9,10,11,12. However, the association between prediabetes and RNFL thickness was not statistically significant, which is most likely due to insufficient statistical power. We, therefore, additionally tested for a linear trend with GMS deterioration by using the statistically more powerful P for trend analysis45, which was consistent with a linear decrease in RNFL thickness with more adverse GMS. In support, all measures of glycaemia, regardless of whether they reflect shorter (i.e., FPG, 2-h PG, and HbA1c) or longer (i.e., SAF and duration of diabetes) exposure, were consistently linearly associated with RNFL thickness.
Similarly, a likely explanation why the association between CGM-assessed standard deviation and RNFL thickness was not statistically significant is that the statistical power to detect any such association was too low46. Indeed, we observed that the association between IGP and RNFL thickness, which included almost fourfold the number of participants (n = 2407 versus n = 622), was statistically significant. Moreover, the strength of the associations of IGP and CGM-assessed standard deviation with RNFL thickness were numerically analogous.
A probable explanation why the association between daily glucose variability and RNFL thickness was not materially altered after additional adjustment for measures of average glycaemia is that daily glucose variability, measures of glycaemia, and GMS represent different underlying constructs26. While daily glucose variability reflects oscillating glucose levels, other measures under study reflect exposure to average chronic levels of glycaemia. Mechanistically, substantial glucose fluctuations entail hyperglycaemic peaks, hypoglycaemic nadirs (in individuals with type 2 diabetes treated with agents that can induce hypoglycaemia), or both, which are thought to be potent inducers of retinal ganglion cell apoptosis26,44. Whereas hyperglycaemic peaks may be highly neurotoxic, hypoglycaemic nadirs likely hamper retinal ganglion cell metabolism as their key nutrient is glucose44.
Our findings can have several implications for clinical practice. First, the strength of the association between type 2 diabetes and RNFL thickness corresponds with 15 years of aging and, thus, indicates that with respect to neurodegeneration substantial “additional aging” occurs in individuals with type 2 diabetes (Supplemental Table S14 shows how this comparison was calculated). Second, RNFL thickness may be a biomarker for the identification of individuals at risk of retinopathy and neuropathy. Use of RNFL thickness measurement is feasible because RNFL thickness assessment is non-invasive2, relatively inexpensive2 and easier to perform than other tests of early neuronal dysfunction such as 24-h electrocardiogram9, magnetic resonance imaging10,11, or electromyography12. Indeed, RNFL thickness has been found to be a promising early biomarker for other neurodegenerative diseases (e.g., multiple sclerosis)47. Third, early glycaemic control, possibly already in prediabetes, is likely crucial in the early prevention of microvascular and neuronal complications5. Last, our findings add to growing evidence that control of daily glucose variability besides mean glucose concentrations may be important to prevent microvascular complications48.
Strengths of this study are (1) the large size of this population-based cohort with oversampling of individuals with type 2 diabetes, which enabled accurate comparison of individuals with and without diabetes; (2) the extensive number of potential confounders that were considered; (3) the use of state-of-the-art and novel methods (e.g., CGM)26 to assess all variables included in this study; and (4) the considerable number of additional analyses, which generally yielded consistent findings.
The study has certain limitations. First, due to the cross-sectional nature of the study, causal inferences should be made with caution49. Mechanistically, hyperglycaemia may not only lead to neurodegeneration but the reverse may also be true, thus causing a vicious cycle. Intact neurovascular interaction is required for normal microvascular function and impaired microvascular function may aggravate hyperglycaemia5,50. Second, we may have underestimated the strength of the associations of GMS, measures of glycaemia, and daily glucose variability with RNFL thickness if such an association was similar or stronger in participants that were excluded from the study population (who generally tend to be less healthy)51. The 2-h post-load glucose and IGP results are most susceptible to this form of selection bias, as no data was available in individuals with the most therapy-intensive diabetes because they were excluded from undergoing an OGTT. Such range restriction may lead to underestimated associations51. Third, a single OGTT may misclassify GMS, especially in individuals with prediabetes. Because individuals classified with prediabetes based on their first OGTT are relatively more prone to receive a NGM classification based on their second OGTT52, this would likely lead to an underestimation of the association with RNFL thickness in the prediabetes group. Fourth, although we took an extensive set of confounders into account, we cannot fully exclude bias due to unmeasured confounding (e.g., environmental factors such as air pollution)53. Fifth, due to the relatively low numbers of participants with data on CGM-based glycaemic indices (n = 622)), and—to a lesser extent—IGP (n = 2407), statistical power of analyses with these determinants was reduced compared to statistical power of analyses with GMS and measures of glycaemia (n = 5132 to n = 5455)46. Last, we studied Caucasian individuals aged 40–75 years with access to high-quality diabetes care. Therefore, the generalizability of our results to other populations requires further study.
In summary, the present population-based study demonstrated that adverse GMS, greater glycaemia, and greater daily glucose variability are associated with lower RNFL thickness, i.e., neurodegeneration, independent of demographics, cardiovascular risk factors, and lifestyle risk factors. Hence, these results suggest that RNFL thickness may be an early biomarker for the identification of individuals at risk of retinopathy and neuropathy. Additionally, the combination of early glycaemic monitoring and early glucose-lowering treatment, possibly already in prediabetes, may contribute to the prevention of RNFL thinning and, ultimately, retinopathy and neuropathy.
Data are available from The Maastricht Study for any researcher who meets the criteria for access to confidential data; the corresponding author may be contacted to request data.
Retinal nerve fibre layer thickness
Glucose metabolism status
Oral glucose tolerance test
Normal glucose metabolism
Fasting plasma glucose
Skin autofluorescence
Continuous glucose monitoring
Glucose variability
Incremental glucose peak
Optical coherence tomography
High-density lipid
Body-mass index
Blood pressure
Estimated glomerular filtration rate
Zakin, E., Abrams, R. & Simpson, D. M. Diabetic neuropathy. Semin. Neurol. 39, 560–569 (2019).
PubMed Google Scholar
Srinivasan, S. et al. Optical coherence tomography predicts 4-year incident diabetic neuropathy. Ophthalmic. Physiol. Opt. 37, 451–459 (2017).
PubMed Google Scholar
De Clerck, E. E. et al. New ophthalmologic imaging techniques for detection and monitoring of neurodegenerative changes in diabetes: A systematic review. Lancet Diabetes Endocrinol. 3, 653–663 (2015).
PubMed Google Scholar
Tomlinson, D. R. & Gardiner, N. J. Glucose neurotoxicity. Nat. Rev. Neurosci. 9, 36–45 (2008).
CAS PubMed Google Scholar
Stehouwer, C. D. A. Microvascular dysfunction and hyperglycemia: A vicious cycle with widespread consequences. Diabetes 67, 1729–1741 (2018).
CAS PubMed Google Scholar
Brownlee, M. The pathobiology of diabetic complications: A unifying mechanism. Diabetes 54, 1615–1625 (2005).
CAS PubMed Google Scholar
Haffner, S. M., Stern, M. P., Hazuda, H. P., Mitchell, B. D. & Patterson, J. K. Cardiovascular risk factors in confirmed prediabetic individuals. Does the clock for coronary heart disease start ticking before the onset of clinical diabetes?. JAMA 263, 2893–2898 (1990).
CAS PubMed Google Scholar
Wong, M. S. et al. The Singapore impaired glucose tolerance follow-up study: Does the ticking clock go backward as well as forward?. Diabetes Care 26, 3024–3030 (2003).
PubMed Google Scholar
Coopmans, C. et al. Both prediabetes and type 2 diabetes are associated with lower heart rate variability: The Maastricht Study. Diabetes Care 43, 1126–1133 (2020).
CAS PubMed Google Scholar
van Agtmaal, M. J. M. et al. Prediabetes Is associated with structural brain abnormalities: The Maastricht Study. Diabetes Care 41, 2535–2543 (2018).
PubMed Google Scholar
Vergoossen, L. W. et al. White matter connectivity abnormalities in prediabetes and type 2 diabetes: The Maastricht Study. Diabetes Care 43, 201–208 (2020).
CAS PubMed Google Scholar
van der Velde, J. et al. Cardiometabolic risk factors as determinants of peripheral nerve function: The Maastricht Study. Diabetologia 63, 1648–1658 (2020).
PubMed PubMed Central Google Scholar
Huru, J., Leiviska, I., Saarela, V. & Liinamaa, M. J. Prediabetes influences the structure of the macula: Thinning of the macula in the Northern Finland Birth Cohort. Br. J. Ophthalmol. (2020).
Hashimoto, S. et al. Foveal and macular thickness in a Japanese population: The Hisayama Study. Ophthal. Epidemiol. 23, 202–208 (2016).
Google Scholar
Cheung, C. Y. et al. Determinants of quantitative optic nerve measurements using spectral domain optical coherence tomography in a population-based sample of non-glaucomatous subjects. Invest. Ophthalmol. Vis. Sci. 52, 9629–9635 (2011).
PubMed Google Scholar
Mauschitz, M. M., Holz, F. G., Finger, R. P. & Breteler, M. M. B. Determinants of macular layers and optic disc characteristics on SD-OCT: The Rhineland Study. Transl. Vis. Sci. Technol. 8, 34 (2019).
PubMed PubMed Central Google Scholar
Mauschitz, M. M. et al. Systemic and ocular determinants of peripapillary retinal nerve fiber layer thickness measurements in the European eye epidemiology (E3) population. Ophthalmology 125, 1526–1536 (2018).
PubMed Google Scholar
Lamparter, J. et al. Association of ocular, cardiovascular, morphometric and lifestyle parameters with retinal nerve fibre layer thickness. PLoS ONE 13, e0197682 (2018).
PubMed PubMed Central Google Scholar
Ho, H. et al. Retinal nerve fiber layer thickness in a multiethnic normal asian population: The Singapore Epidemiology of Eye Diseases Study. Ophthalmology 126, 702–711 (2019).
PubMed Google Scholar
Fujiwara, K. et al. Glucose tolerance levels and circumpapillary retinal nerve fiber layer thickness in a general Japanese population: The Hisayama Study. Am. J. Ophthalmol. 205, 140–146 (2019).
CAS PubMed Google Scholar
Duan, X. R. et al. Normal macular thickness measurements using optical coherence tomography in healthy eyes of adult Chinese persons: The Handan Eye Study. Ophthalmology 117, 1585–1594 (2010).
PubMed Google Scholar
De Clerck, E. E. B. et al. Macular thinning in prediabetes or type 2 diabetes without diabetic retinopathy: The Maastricht Study. Acta Ophthalmol. 96, 174–182 (2018).
PubMed Google Scholar
De Clerck, E. E. B. et al. Loss of temporal peripapillary retinal nerve fibers in prediabetes or type 2 diabetes without diabetic retinopathy: The Maastricht Study. Invest. Ophth. Vis. Sci. 58, 1017–1027 (2017).
Google Scholar
Khawaja, A. P. et al. Comparison of associations with different macular inner retinal thickness parameters in a large cohort: The UK Biobank. Ophthalmology 127, 62–71 (2020).
PubMed Google Scholar
Jud, P. & Sourij, H. Therapeutic options to reduce advanced glycation end products in patients with diabetes mellitus: A review. Diabetes Res. Clin. Pract. 148, 54–63 (2019).
CAS PubMed Google Scholar
Danne, T. et al. International consensus on use of continuous glucose monitoring. Diabetes Care 40, 1631–1640 (2017).
PubMed PubMed Central Google Scholar
Schram, M. T. et al. The Maastricht Study: An extensive phenotyping study on determinants of type 2 diabetes, its complications and its comorbidities. Eur. J. Epidemiol. 29, 439–451 (2014).
PubMed Google Scholar
World Health Organization. Definition and Diagnosis of Diabetes Mellitus and Intermediate Hyperglycemia: Report of a WHO/IDF Consultation. Geneva, Switzerland (2006).
van Eupen, M. G. et al. Skin autofluorescence and pentosidine are associated with aortic stiffening: The Maastricht Study. Hypertension 68, 956–963 (2016).
PubMed Google Scholar
Foreman, Y. D. et al. Glucose variability assessed with continuous glucose monitoring: Reliability, reference values, and correlations with established glycemic indices-The Maastricht Study. Diabetes Technol. Ther. 22, 395–403 (2020).
CAS PubMed Google Scholar
Hong, J. T. et al. Retinal nerve fiber layer measurement variability with spectral domain optical coherence tomography. Korean J. Ophthalmol. 26, 32–38 (2012).
PubMed PubMed Central Google Scholar
Qi, Y. et al. Adulthood socioeconomic position and type 2 diabetes mellitus-a comparison of education, occupation, income, and material deprivation: The Maastricht Study. Int. J. Environ. Res. Public Health 16, 1435 (2019).
PubMed Central Google Scholar
Looman, M. et al. Development and evaluation of the Dutch Healthy Diet index 2015. Public Health Nutr. 20, 2289–2299 (2017).
PubMed Google Scholar
van Dongen, M. C. et al. The Maastricht FFQ: Development and validation of a comprehensive food frequency questionnaire for the Maastricht Study. Nutrition 62, 39–46 (2019).
PubMed Google Scholar
van der Berg, J. D. et al. Identifying waking time in 24-h accelerometry data in adults using an automated algorithm. J. Sports Sci. 34, 1867–1873 (2016).
PubMed Google Scholar
van Dooren, F. E. et al. Associations of low grade inflammation and endothelial dysfunction with depression—The Maastricht Study. Brain Behav. Immun. 56, 390–396 (2016).
PubMed Google Scholar
Rodbard, D. The challenges of measuring glycemic variability. J. Diabetes Sci. Technol. 6, 712–715 (2012).
PubMed PubMed Central Google Scholar
Bikbov, M. M. et al. Axial length and its associations in a Russian population: The Ural Eye and Medical Study. PLoS ONE 14, e0211186 (2019).
CAS PubMed PubMed Central Google Scholar
Haeger, A., Costa, A. S., Schulz, J. B. & Reetz, K. Cerebral changes improved by physical activity during cognitive decline: A systematic review on MRI studies. Neuroimage Clin. 23, 101933 (2019).
PubMed PubMed Central Google Scholar
Poulose, S. M., Miller, M. G., Scott, T. & Shukitt-Hale, B. Nutritional factors affecting adult neurogenesis and cognitive function. Adv. Nutr. 8, 804–811 (2017).
CAS PubMed PubMed Central Google Scholar
Schisterman, E. F., Cole, S. R. & Platt, R. W. Overadjustment bias and unnecessary adjustment in epidemiologic studies. Epidemiology 20, 488–495 (2009).
PubMed PubMed Central Google Scholar
Battelino, T. et al. Clinical targets for continuous glucose monitoring data interpretation: Recommendations from the international consensus on time in range. Diabetes Care 42, 1593–1603 (2019).
PubMed PubMed Central Google Scholar
Vatcheva, K. P., Lee, M., McCormick, J. B. & Rahbar, M. H. Multicollinearity in regression analyses conducted in epidemiologic studies. Epidemiology (Sunnyvale). 6, 2 (2016).
Google Scholar
Yu, D. Y. et al. Retinal ganglion cells: Energetics, compartmentation, axonal transport, cytoskeletons and vulnerability. Prog. Retin. Eye Res. 36, 217–246 (2013).
CAS PubMed Google Scholar
Greenland, S. Avoiding power loss associated with categorization and ordinal scores in dose–response and trend analysis. Epidemiology 6, 450–454 (1995).
CAS PubMed Google Scholar
Jones, S. R., Carley, S. & Harrison, M. An introduction to power and sample size estimation. Emerg. Med. J. 20, 453–458 (2003).
CAS PubMed PubMed Central Google Scholar
Britze, J. & Frederiksen, J. L. Optical coherence tomography in multiple sclerosis. Eye (Lond) 32, 884–888 (2018).
Google Scholar
Ceriello, A. Glucose variability and diabetic complications: Is it time to treat?. Diabetes Care 43, 1169–1171 (2020).
PubMed Google Scholar
Rothman, K. J. & Greenland, S. Causation and causal inference in epidemiology. Am. J. Public Health 95(Suppl 1), S144–S150 (2005).
PubMed Google Scholar
Toth, P., Tarantini, S., Csiszar, A. & Ungvari, Z. Functional vascular contributions to cognitive impairment and dementia: Mechanisms and consequences of cerebral autoregulatory dysfunction, endothelial impairment, and neurovascular uncoupling in aging. Am. J. Physiol. Heart Circ. Physiol. 312, H1–H20 (2017).
PubMed Google Scholar
Bland, J. M. & Altman, D. G. Correlation in restricted ranges of data. BMJ 342, d556 (2011).
PubMed Google Scholar
Mooy, J. M. et al. Intra-individual variation of glucose, specific insulin and proinsulin concentrations measured by two oral glucose tolerance tests in a general Caucasian population: The Hoorn Study. Diabetologia 39, 298–305 (1996).
CAS PubMed Google Scholar
Sharma, S. et al. Effect of environmental toxicants on neuronal functions. Environ. Sci. Pollut. Res. Int. 27, 44906–44921 (2020).
CAS PubMed Google Scholar
Download references
The authors would like to acknowledge ZIO foundation (Vereniging Regionale HuisartsenZorg Heuvelland) for their contribution to The Maastricht Study. The researchers are indebted to the participants for their willingness to participate in the study.
This study was supported by the European Regional Development Fund via OP-Zuid, the Province of Limburg, the Dutch Ministry of Economic Affairs (grant 31O.041), Stichting De Weijerhorst (Maastricht, the Netherlands), the Pearl String Initiative Diabetes (Amsterdam, the Netherlands), the Cardiovascular Center (CVC, Maastricht, the Netherlands), CARIM School for Cardiovascular Diseases (Maastricht, the Netherlands), CAPHRI School for Public Health and Primary Care (Maastricht, the Netherlands), NUTRIM School for Nutrition and Translational Research in Metabolism (Maastricht, the Netherlands), Stichting Annadal (Maastricht, the Netherlands), Health Foundation Limburg (Maastricht, the Netherlands), Perimed (Järfälla, Sweden), Diabetes fonds grant 2016.22.1878 (Amersfoort, The Netherlands), Oogfonds (Utrecht, The Netherlands) and by unrestricted grants from Janssen-Cilag B.V. (Tilburg, the Netherlands), Novo Nordisk Farma B.V. (Alphen aan den Rijn, the Netherlands), Sanofi-Aventis Netherlands B.V. (Gouda, the Netherlands), and Medtronic (Tolochenaz, Switzerland).
These authors contributed equally: Frank C. T. van der Heide and Yuri D. Foreman.
CARIM School for Cardiovascular Diseases, Maastricht University (UM), Maastricht, The Netherlands
Frank C. T. van der Heide, Yuri D. Foreman, Iris W. M. Franken, Ronald M. A. Henry, Abraham A. Kroon, Pieter C. Dagnelie, Simone J. P. M. Eussen, Miranda T. Schram, Carla J. H. van der Kallen, Marleen M. J. van Greevenbroek, Casper G. Schalkwijk, Nicolaas C. Schaper, Martijn C. G. J. Brouwers & Coen D. A. Stehouwer
Department of Internal Medicine, Maastricht University Medical Center+ (MUMC+), P. Debyelaan 25, P.O. Box 5800, 6202AZ, Maastricht, The Netherlands
Frank C. T. van der Heide, Yuri D. Foreman, Iris W. M. Franken, Ronald M. A. Henry, Abraham A. Kroon, Pieter C. Dagnelie, Miranda T. Schram, Carla J. H. van der Kallen, Marleen M. J. van Greevenbroek, Casper G. Schalkwijk & Coen D. A. Stehouwer
Heart and Vascular Center, MUMC+ Maastricht, Maastricht, The Netherlands
Ronald M. A. Henry, Abraham A. Kroon & Miranda T. Schram
Department of Epidemiology, UM, Maastricht, The Netherlands
Simone J. P. M. Eussen
University Eye Clinic Maastricht, MUMC+, Maastricht, The Netherlands
Tos T. J. M. Berendschot, Jan S. A. G. Schouten & Carroll A. B. Webers
Department of Ophthalmology, Canisius Wilhelmina Hospital, Nijmegen, The Netherlands
Jan S. A. G. Schouten
Department of Complex Genetics and Epidemiology, NUTRIM School for Nutrition and Translational Research in Metabolism, UM, Maastricht, The Netherlands
Anke Wesselius
CAPHRI Care and Public Health Research Institute, UM, Maastricht, The Netherlands
Nicolaas C. Schaper
Division of Endocrinology and Metabolic Disease, Department of Internal Medicine, Maastricht University Medical Centre+, Maastricht, The Netherlands
Nicolaas C. Schaper & Martijn C. G. J. Brouwers
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
F.C.T.vd.H. and Y.D.F. contributed to conception and design, participated in acquisition of data, analyzed and interpreted data, drafted the manuscript (with C.D.A.S., N.C.S., M.C.G.J.B., R.M.A.H., and J.S.A.G.S.), revised the manuscript critically for important intellectual content, and provided final approval of the version to be published. F.C.T. vd.H and Y.D.F. also are the guarantors of this work and, as such, had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. I.W.M.F., A.J.H.M.H, A.A.K., P.C.D., S.J.P.M.E., T.J.M.B., C.A.B.W., M.T.S., C.J.H.vd.K., M.M.J.v.G., A.W., and C.G.S. contributed to conception and design, revised the manuscript critically for important intellectual content, and provided final approval of the version to be published.
Correspondence to Frank C. T. van der Heide.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
van der Heide, F.C.T., Foreman, Y.D., Franken, I.W.M. et al. (Pre)diabetes, glycemia, and daily glucose variability are associated with retinal nerve fiber layer thickness in The Maastricht Study. Sci Rep 12, 17750 (2022). https://doi.org/10.1038/s41598-022-22748-2
Download citation
Received: 31 March 2022
Accepted: 19 October 2022
Published: 22 October 2022
DOI: https://doi.org/10.1038/s41598-022-22748-2
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Advertisement
© 2022 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.