Circumpapillary retinal nerve fiber layer thickness (cpRNFLT) as assessed by spectral domain optical coherence tomography (SD-OCT) is a new technique used for the detection and evaluation of glaucoma and other optic neuropathies. Before translating cpRNFLT into clinics, it is crucially important to investigate anthropometric, biochemical, and clinical parameters potentially affecting cpRNFLT in a large population-based dataset.
The population-based LIFE-Adult Study randomly selected 10,000 participants from the population registry of Leipzig, Germany. All participants underwent standardized systemic assessment of various cardiometabolic risk markers and ocular imaging, including cpRNFLT measurement using SD-OCT (Spectralis, Heidelberg Engineering). After employing strict SD-OCT quality criteria, 8952 individuals were analyzed. Multivariable linear regression analyses were used to evaluate the independent associations of various cardiometabolic risk markers with sector-specific cpRNFLT. For significant markers, the relative strength of the observed associations was compared to each other to identify the most relevant factors influencing cpRNFLT. In all analyses, the false discovery rate method for multiple comparisons was applied.
In the entire cohort, female subjects had significantly thicker global and also sectoral cpRNFLT compared to male subjects (p < 0.05). Multivariable linear regression analyses revealed a significant and independent association between global and sectoral cpRNFLT with biomarkers of renal function and lipid profile. Thus, thinner cpRNFLT was associated with worse renal function as assessed by cystatin C and estimated glomerular filtration rate. Furthermore, an adverse lipid profile (i.e., low high-density lipoprotein (HDL) cholesterol, as well as high total, high non-HDL, high low-density lipoprotein cholesterol, and high apolipoprotein B) was independently and statistically significantly related to thicker cpRNFLT. In contrast, we do not observe a significant association between cpRNFLT and markers of inflammation, glucose homeostasis, liver function, blood pressure, or obesity in our sector-specific analysis and globally.
Markers of renal function and lipid metabolism are predictors of sectoral cpRNFLT in a large and deeply phenotyped population-based study independently of previously established covariates. Future studies on cpRNFLT should include these biomarkers and need to investigate whether incorporation will improve the diagnosis of early eye diseases based on cpRNFLT.
Retinal nerve fiber layer (RNFL) defects are early signs of glaucoma and optic disc deformation . RNFL thickness is, therefore, a major tool in the evaluation of glaucoma and other optic neuropathies . Spectral domain optical coherence tomography (SD-OCT) is an appropriate , non-invasive, in vivo technique for the analysis of the optic nerve, and recent advances have allowed improved image quality for circumpapillary RNFL thickness (cpRNFLT) . Very recently, cpRNFLT has been associated with distinct, basic anthropometric, and biochemical measures in different studies. For instance, Ho et al.  demonstrate a positive association of global cpRNFLT with low-density lipoprotein (LDL) cholesterol and a negative correlation with diabetes prevalence in three different Asian ethnic cohorts. Furthermore, age and a history of stroke or hypertension were negatively, whereas smoking status was positively, related to global cpRNFLT in a cross-sectional meta-analysis of eight European, population-based studies . In contrast, Lamparter et al.  did not find an independent association between global cpRNFLT and cardiovascular disease in multivariable analyses in the Gutenberg Health Study. Taken together, the association of global cpRNFLT shows conflicting results with the presence of cardiometabolic disease states. Before translating the cpRNFLT method into clinics, it is important to investigate anthropometric, biochemical, and clinical parameters potentially affecting cpRNFLT independent of other well-established predictors, i.e., age, sex, and scan radius. Furthermore, other factors influencing cpRNFLT need to be carefully defined to aid the early diagnosis of eye diseases and to prevent misclassification of impaired cpRNFLT due to other clinical and biochemical biomarkers. However, previous studies on cpRNFLT show the following limitations: they (a) included cohorts of smaller sample size; (b) have analyzed global cpRNFLT but not sector-specific data; (c) excluded subjects with different cardiometabolic disease states, e.g., type 2 diabetes or hypertension; (d) did not include a wide range of anthropometric, biochemical, and cardiometabolic markers and other patient-level data; and (e) did not use thoroughly adjusted multivariable models to investigate the independent predictors of cpRNFLT.
We, therefore, investigated a large panel of different anthropometric and cardiometabolic biomarkers and a wide range of clinical phenotypes and their associations with the sector-specific cpRNFLT profile measured by SD-OCT in a large (N = 8952 subjects), unselected, and deeply phenotyped population-based study in Germany. We have applied a highly standardized ophthalmologic and non-ophthalmologic investigation procedure and statistical adjustment with correction for multiple testing.
This analysis is part of the population-based LIFE-Adult Study conducted by the Leipzig Research Centre for Civilization Diseases at Leipzig University between August 2011 and November 2014 . The LIFE-Adult Study includes 10,000 randomly selected participants from the population registry of just over half a million inhabitants of Leipzig, a city located in the east of Germany.
The LIFE-Adult Study recruitment was performed in an age- and sex-stratified manner mainly focusing on subjects with an age between 40 and 79 years . For this purpose, the overall population consisted of 9600 subjects between 40 and 79 years of age, as well as 400 subjects between 19 and 39 years of age. Each age interval (by decade) was balanced with respect to the number of subjects and sex. The study was approved by the Ethical Committee at the Medical Faculty of Leipzig University (approval number: 263-2009-14122009) and adheres to the Declaration of Helsinki and all federal and state laws. Prior to inclusion, informed written consent was obtained from all participants.
Data collection/inclusion and exclusion criteria
During the baseline examination, study participants were deeply phenotyped, including ophthalmological image data, structured interviews, questionnaires, physical examinations, and blood and urine tests . As part of the ophthalmic assessment, SD-OCT imaging (Spectralis, Heidelberg Engineering, Heidelberg, Germany) was performed, yielding cpRNFLT scans around the optic nerve head. The location of the cpRNFLT circle and the coordinate system have been described previously . We excluded subjects with missing SD-OCT scans (excluded N = 931) or SD-OCT scans using the following quality criteria: (1) B-scan number per location < 50, (2) signal to noise ratio < 20 dB, and (3) missing or unreliable RNFLT A-scans > 5% (excluded N = 117). For the remaining 8952 subjects, one eye was randomly selected if both eyes of an included subject were reliable . For validation analyses, we classified optic nerve head (ONH) abnormalities if any of the following were present: excavation (suspected glaucoma [i.e., violation of the inferior-superior-nasal-temporal rule, vertically oval with cup-to-disc ratio > 0.7], optic disc pit, or coloboma of the optic disc), optic disc hemorrhage, neovascularization, optic atrophy, sectoral paleness, ONH swelling, papilledema, or optic disc drusen . Furthermore, patient information on a previous diagnosis of glaucoma, as well as glaucoma medication, were collected.
Anthropometric and biochemical markers
Classical anthropometric (e.g., body mass index [BMI], waist and hip circumferences, blood pressure) measurements were assessed according to standardized procedures by trained study nurses. In all subjects, fasting blood samples were drawn routinely and a panel of laboratory tests was performed on the day of sample collection . The biomarkers of the panel have been described previously  and included measurements of total cholesterol, high-density lipoprotein (HDL), low-density lipoprotein (LDL) cholesterol, triglycerides (TG), apolipoprotein (apo) B, apoA1, lipoprotein (Lp) (a), glucose, insulin, glycated hemoglobin (HbA1c), liver enzymes, interleukin-6, high-sensitivity C-reactive protein (hsCRP), cystatin C, and urinary albumin and creatinine, all being quantified in a central lab by standard methods . In all subjects, the estimated glomerular filtration rate (eGFR) was calculated using the cystatin C-based chronic kidney disease (CKD) epidemiology collaboration equation . As the aim of the current study was to investigate the associations between several cardiometabolic risk markers and cpRNFLT, we have used only the cystatin C-based equation which has been recently reported as the best equation for the assessment of cardiovascular risk . CKD status was defined as a urinary albumin/creatinine ratio ≥ 30 mg/g and/or a decreased eGFR < 60 ml/min/1.73 m2, and the cohort was divided into five eGFR categories (i.e., G1-G5 combining G3a and G3b into one G3 category), as well as four CKD risk categories (i.e., low, moderately increased, high, and very high risk), according to KDIGO .
All statistical analyses were performed in R environment using version 3.5 (R Foundation for Statistical Computing, Vienna, Austria). For comparisons between female and male subjects, the unpaired Student t-test (for continuous variables) or chi-squared test (for categorical variables) were used, respectively.
As a next step, we examined the associations of various anthropometric, as well as cardiometabolic, biomarkers on the sectoral RNFLT. For this purpose, multivariable linear regression analyses were carried out for individual markers adjusted for age, sex, and scanning circle radius in all models. Using these covariates as independent variables in the respective models, the association of each marker with sectoral cpRNFLT (dependent variable) were separately calculated for the temporal (T), supero-temporal (TS), supero-nasal (NS), nasal (N), infero-nasal (NI), and infero-temporal (TI) sectors, and globally (G).
Scanning circle radius was included as an independent variable in all models, since eye size and optical characteristics of the human lens confound the cpRNFLT measurement [12,13,14]. The true scanning circle radius (mm) is estimated from the focus settings used by the Spectralis machine, according to a widely used model . As sex  and age  were shown to affect cpRNFLT, these markers were also included as independent covariates in each model.
We next sought to compare the relative strength of the associations of all biomarkers with sectoral cpRNFLT. Therefore, a sectoral and global heatmap of standardized β values from multivariable analyses for each sectoral cpRNFLT was produced for all biomarkers, and standardized β values were employed in the figure’s color code representing the strength of each association.
As a sensitivity analysis, we further validated the results of the linear regression analyses for lipid markers with sectoral cpRNFLT by stratifying the cohort into subjects on statin treatment compared to non-statin users. For this purpose, we used the Anatomical Therapeutic Chemical (ATC) classification codes to extract participants treated with 3-hydroxy-3-methylglutaryl coenzyme A reductase inhibitors (i.e., statins), thereby reducing cholesterol synthesis. To investigate the potential mediating effects of smoking status on the association between cpRNFLT and the lipid profile, the Bayesian information criterion difference (ΔBIC) was computed for model comparisons. For this purpose, two different linear regression models were calculated with age, sex, measurement radius, and the respective lipid marker, as regressors (model A), as well as an additional model comprising of model A + smoking status (model B). The BIC difference (ΔBIC) was calculated by ΔBIC = BICmodel A − BICmodel B. A ΔBIC > 2 was regarded as statistically relevant according to Madrigal-González et al. , as well as Kass and Raftery .
In all other analyses, a p value < 0.05 was considered as statistically significant. The false discovery rate (FDR) method was applied to correct all p values for multiple comparisons.
Baseline characteristics of the entire study population (N = 8952)
Baseline characteristics of cohort stratified by sex are shown in Table 1. The mean ± standard deviation age of the total population was 57.8 ± 12.4 years. In the entire cohort, about 1243 (13.9%) patients had a diabetes. Hypertension was present in 4444 (49.6%) subjects, and 1130 (12.6 %) participants received statin treatment. Moreover, 1875 (20.9%) subjects were current smokers. Female subjects had lower markers of obesity (i.e., BMI, waist-to-hip ratio), glucose homeostasis (i.e., fasting glucose, fasting insulin, HbA1c), blood pressure, albuminuria, interleukin-6, and liver enzymes compared to male participants (all p < 0.05; Table 1). In contrast, women had higher total cholesterol, HDL cholesterol, apoA1, Lp(a), and hsCRP, compared to men (all p < 0.05; Table 1). Renal function (eGFR), LDL cholesterol, and alkaline phosphatase did not depend on sex (all p > 0.05; Table 1). The mean thicknesses for the respective sectoral cpRNFLT and global cpRNFLT are depicted in Table 2. Female subjects had significantly thicker global cpRNFLT and also regional cpRNFLT in all sectors compared to male participants (all p < 0.05; Table 2), except NS sector (p = 0.067; Table 2).
Associations between sector-specific cpRNFLT and cardiometabolic biomarkers
Visual inspection of selected clinical parameters of renal function (i.e., eGFR) and lipid profile (i.e., HDL cholesterol, non-HDL cholesterol) revealed that there is a linear relationship between these biomarkers and global cpRNFLT (Additional file 1: Fig. S1). Therefore, linear regression models were used in all subsequent analyses. Multivariable associations with cpRNFLT averaged across each of the six respective sectors, as well as the global mean, were separately calculated for each biomarker and sector with adjustment for potential covariates found in previous analyses, i.e., age [6, 14], sex , and scanning circle radius [12,13,14]. A significant and independent association between global cpRNFLT and BMI, cystatin C, eGFR, eGFR category, and presence of CKD, as well as total cholesterol, HDL cholesterol, non-HDL cholesterol, LDL cholesterol, and ApoB, was found, respectively (Table 3, Fig. 1). These results indicate that renal function and lipid profile contribute to cpRNFLT. As correction for multiple comparisons potentially increases type II errors, analysis was repeated without FDR correction and revealed comparable results (data not shown), which demonstrates that the presence and absence of effects in our results are distinct enough from each other that shifting significance criteria within the limits of multiple comparison adjustments does not change any conclusions. When the entire analysis was further adjusted for SD-OCT-derived clinical and subclinical ONH abnormalities (including glaucoma and other ONH diseases), as well as patient-reported glaucoma diagnosis and medication, the associations of all investigated biomarkers and cpRNFLT were virtually unchanged in terms of effect size and direction (Additional file 2: Table S1).
Renal function and sector-specific cpRNFLT
The associations of cpRNFLT with markers of renal function remained virtually the same in all three temporal sectors (i.e., T, TS, TI), as well as NI (Table 3, Fig. 1). In more detail, global, temporal, and NI cpRNFLT were significantly and independently associated with cystatin C, eGFR, eGFR category, and CKD status (Fig. 1). Thus, eGFR was positively associated with cpRNFLT in the respective sectors indicating that better renal function is linked to thicker cpRNFLT (Fig. 1). Conversely, cystatin C, eGFR category, and CKD status were inversely related to cpRNFLT in the respective sectors with the strongest associations observed in TI and globally (Fig. 1). Per 1 mg/l increase in cystatin C, there was a decrease of − 2.2 μm in global cpRNFLT (Table 3).
Lipid profile and sector-specific cpRNFLT
Lipid markers were statistically significantly and independently associated with sector-specific cpRNFLT for total cholesterol, HDL cholesterol, non-HDL cholesterol, LDL cholesterol, and ApoB, in sectors TS, TI, NS, NI, and globally (Table 3, Fig. 1).
For ApoB-containing lipid particles (i.e., non-HDL cholesterol, LDL cholesterol, total cholesterol), a positive correlation between each lipid measure and sector-specific cpRNFLT was found except for T and N, and the strongest associations were observed for non-HDL cholesterol and LDL cholesterol (Fig. 1). In contrast, HDL cholesterol was negatively and statistically significantly related to cpRNFLT in TS, NS, TI, and globally (Fig. 1) indicating that HDL cholesterol is inversely related to cpRNFLT compared to ApoB-containing lipid particles. Thus, an adverse lipid profile (i.e., high total cholesterol, high non-HDL cholesterol, high LDL cholesterol, high ApoB, low HDL cholesterol) is independently and statistically significantly related to thicker cpRNFL (Fig. 1). Clinically, per 1 mmol/l increase in non-HDL cholesterol, there was an increase of 0.5 μm in global cpRNFLT (Table 3).
Other biomarkers and sector-specific cpRNFLT
Current smoking was significantly and independently related to cpRNFLT in sector T and all nasal sectors (i.e., NS, NI, N) (Table 3, Fig. 1). The positive association of global cpRNFLT and BMI was also confirmed in sectors NS and NI (Table 3, Fig. 1). Furthermore, diabetes status and gamma-glutamyltransferase were significantly associated with TI and NI of the cpRNFLT, respectively (Table 3, Fig. 1).
Although we have not found an association between cpRNFLT and inflammatory markers (Table 3, Fig. 1), inflammation might mediate the observed results. Therefore, all analyses were additionally adjusted for hsCRP, and results remained virtually unchanged with additional adjustment for inflammation (Additional file 3: Table S2). Furthermore, when using interleukin-6 instead of hsCRP, comparable associations were observed (data not shown).
Sensitivity analyses—CKD risk groups, statin treatment, smoking status, and glaucoma status
To verify our results obtained from multivariable analyses for markers of renal function, we have stratified the entire cohort based on the CKD risk groups according to KDIGO . Subjects with advanced CKD (i.e., moderate/high/very high risk) had a significantly thinner global and sectoral cpRNFL compared to participants without CKD (low risk; Additional file 4: Table S3) supporting our findings from regression models (Table 3).
Using a similar approach, patients on statin treatment showed a significantly thinner global, as well as temporal, cpRNFL compared to non-statin users (Additional file 5: Table S4). When multivariable analyses of cpRNFLT and lipid profile were stratified by statin usage, associations between lipid profile and cpRNFLT remained similar for non-statin users compared to the entire cohort (Additional file 6: Table S5). In contrast, patients on statin treatment only showed a positive, significant, and independent association of non-HDL cholesterol with cpRNFLT in sector NS, as well as globally (Additional file 6: Table S5).
We further investigated whether smoking status mediates our observed results of cpRNFLT with lipid profile as both smoking and lipid profile were counterintuitively positively associated with cpRNFLT. Using Bayesian information criterion difference (ΔBIC), smoking status does not have an additive effect on linear regression models for cpRNFLT with adverse lipid markers (all ΔBIC < 2; Additional file 7: Table S6).
If including several possible confounders of cpRNFLT (i.e., age, sex, scan radius, smoking status, eGFR, non-HDL cholesterol, and hsCRP) in one regression model for each sector and globally, age, sex, scan radius, eGFR, and non-HDL cholesterol were found to be the most important factors independently contributing to cpRNFLT changes (data now shown). Collectively, these data support our findings of independent associations of renal function, as well as lipid profile, with cpRNFLT (Table 3).
As cpRNFL thinning is associated with glaucoma , we next sought to investigate whether an adverse lipid profile is similarly associated with cpRNFLT in patients with ONH abnormalities, self-reported glaucoma diagnosis, and/or glaucoma medication (N = 1180). In this subset, the associations of lipid markers and cpRNFLT were virtually unchanged (Additional file 8: Table S7) compared to the entire cohort (Table 3).
In the current study, we identify potential anthropometric and cardiometabolic markers which are independently and significantly associated with sector-specific cpRNFLT using a large dataset comprising of 10,000 deeply phenotyped subjects . As we find that different renal and metabolic biomarkers are associated with cpRNFLT, translational implementation of cpRNFLT profiles to clinics and research requires adjustment for these confounders in future studies.
We report that a cystatin C-based calculation of eGFR, i.e., currently the best equation for assessment of cardiovascular risk , was independently associated with cpRNFLT. To the best of our knowledge, no previous large population-based study has investigated the effect of renal function on sectoral cpRNFLT. A small study of 60 patients with diabetic retinopathy vs. 20 healthy controls also showed that thinner RNFLT is correlated with the increase in serum urea and creatinine, i.e., two other surrogate markers of renal dysfunction . Furthermore, data from a Chinese cohort of 1408 patients with type 2 diabetes also revealed a positive relation between eGFR and cpRNFLT . Despite some similar correlations , we refrain from comparing our SD-OCT-derived data to some previous studies evaluating cpRNFLT by fundus photography. In summary, thinner cpRNFLT is a feature of CKD, an established risk factor for the progression to end-stage kidney disease (ESKD) and mortality [10, 22]. As our data show an independent association between renal function and cpRNFLT, future studies need to investigate whether cpRNFLT could predict the risk of progressive CKD and/or mortality, as well as mechanistically identify whether impaired renal function impairs RNFL characteristics similar to brain lesions, e.g., white matter lesions, silent cerebral infarction, microbleeds, and brain atrophy .
A distinct set of adverse lipid markers was significantly and positively associated with cpRNFLT in our cohort. Thus, total cholesterol, non-HDL cholesterol, LDL cholesterol, and circulating apoB are positive predictors of sectoral cpRNFLT in most sectors except for T. It is interesting to note that apoB is the primary apolipoprotein of chylomicrons, very low-density lipoproteins, intermediate-density lipoproteins, and LDL particles, collectively summarized as non-HDL cholesterol. Thus, similar results in our models for non-HDL cholesterol and apoB validate our findings. Furthermore, LDL cholesterol is the main component of total cholesterol, and therefore, the association of total cholesterol with cpRNFLT is weaker but still comparable to apoB, non-HDL cholesterol, and LDL cholesterol. The aforementioned positive associations of adverse lipid particles with cpRNFLT in our large cohort are of interest, as these lipid particles, for instance, total cholesterol  and LDL cholesterol , also predict increased vascular and overall mortality. Furthermore, HDL cholesterol is associated with a beneficial cardiovascular risk profile and reduced mortality . In accordance with our data on the positive correlation of cpRNFLT and adverse apoB-containing lipid particles, HDL cholesterol is negatively associated with sectoral cpRNFLT. Collectively, our data indicate that an adverse lipid profile as assessed by high apoB, high non-HDL cholesterol, high LDL cholesterol, and low HDL cholesterol is independently and statistically significantly related to a thicker cpRNFLT. Importantly, the presence of SD-OCT-derived clinical and subclinical ONH abnormalities (including glaucoma and other ONH diseases), as well as patient-reported glaucoma diagnosis and medication, did not influence the observed results (Additional file 8: Table S7). Pathophysiologically, the retina is capable of the rapid uptake of distinct cholesterol particles from the circulation . Furthermore, many relevant proteins and receptors necessary for uptake, transport, metabolism, synthesis, and efflux of cholesterol and other lipid molecules  are expressed in cells of the human RNFL (i.e., ganglion cells and glia cell types such as astrocytes and Müller cells) [29, 30]. Moreover, lipid-lowering drugs, for instance, statins, are permeable to the blood-retinal barrier , and chronic (i.e., for 6 weeks) simvastatin treatment in mice decreases total retinal cholesterol content by 24% . In more detail, simvastatin reduced retinal cholesterol biosynthesis and to a lesser extent increased retinal uptake of serum cholesterol leading to a reduced retinal cholesterol content . Collectively, retinal lipid synthesis is likely to be the main driver of retinal cholesterol content, but circulating lipids and lipid-lowering compounds also significantly contribute to retinal cholesterol levels [31, 32]. Thus, increased circulating lipid levels hypothetically might induce lipid trafficking from the blood stream into the peripapillary retinal tissue and, concomitantly, increase the volume of cpRNFL by lipid accumulation. Additionally, dysregulated lipid particles could contribute to volume changes and accumulation of cholesterol deposits in the endothelial cells and pericytes of the vasculature around the optic nerve head. Furthermore, retention of lipoproteins in retinal layers potentially induces lipid modifications (i.e., oxidized lipoproteins, different forms of cholesterol) and adverse pro-inflammatory, pro-angiogenic downstream effects—a phenotype similar to atherosclerotic coronary artery disease . It is important to note in this context that an adverse lipid profile cannot prevent glaucoma as both processes differ pathophysiologically. As a consequence, the decay of nerve fibers in glaucoma potentially could be clinically obscured in patients with an adverse lipid profile due to a net effect of “normal” cpRNFLT.
Our observed reduced lipid levels in the circulation, as well as decreased sectoral cpRNFLT, in statin-treated subjects compared to non-statin users (Additional file 5: Table S4), further support the hypothesis of a concentration-difference-driven accumulation of lipids in peripapillary retinal tissue and blood vessels.
Importantly, after stratifying the cohort into statin-treated subjects compared to non-statin users, the direction of the associations remained virtually the same (although statistical significance was partly lost in the statin-treated group) (Additional file 7: Table S6). Potential reasons for the observed effects include an altered lipid metabolism due to statin treatment; disassociation of the cpRNFLT connection to lipid profile by statin treatment, which cannot reverse the cpRNFLT phenotype (in contrast to the statin effects on lipid markers); and/or a reduced number (N = 1130) of included subjects in the subanalysis. As cpRNFL thinning is associated with neuropathies like glaucoma [34, 35] and systemic diseases like diabetes mellitus, even prior to the development of the typical retinal defects related to diabetic retinopathy [36,37,38], the lipid data in the entire cohort may seem counterintuitive. On the other hand, our lipid results are in analogy to the counterintuitive effect of smoking status on cpRNFLT. Here, smoking is also positively associated with the thickening of cpRNFLT in nasal sectors N, NS, and NI in our cohort, similar to an adverse lipid profile. Interestingly, these data are in accordance with the results from a recent meta-analysis  where a positive correlation between current or former smoking and global cpRNFLT was also found in the Rotterdam II and in the Montrachet Study . A further small case-control study by Teberik  also shows numerically thicker RNFL in the three nasal sectors (but the statistical significance has been achieved only in the NI sector). As smoking potentially interferes with lipid profile , we investigated the mediating effects of smoking status and lipid profile. Here, smoking status did not interfere with the observed associations of cpRNFLT with the lipid profile (Additional file 7: Table S6) suggesting separate effects of both variables, i.e., smoking status and lipid profile, on cpRNFLT.
In contrast to the significant associations observed for renal function and dyslipidemia in relation to cpRNFLT, we did not observe significant correlations between cpRNFLT and markers of inflammation, glucose homeostasis, liver function, blood pressure, and obesity in our sector-specific analysis investigating six sectors and globally. While we are aware of the association of diabetes status with thinning of cpRNFLT in some [5, 41] but not all  studies, our large population study suggests that clinicians do not need to consider these cardiometabolic risk markers as potential confounders in sector-specific cpRNFLT examination.
Based on our results, it is tempting to speculate whether incorporating our novel identified RNFLT-linked biomarkers (i.e., renal function, lipid profile, smoking status) additional to conventional (i.e., age, sex, scanning circle radius) parameters will improve the diagnosis of early eye diseases based on cpRNFLT.
This study has several limitations: First, the study population predominantly consisted of European subjects, and therefore, the findings may not be generalizable to populations of different ethnicities. Secondly, the cross-sectional design of this study does not permit causal conclusions. However, strengths of the current study include a large number of deeply phenotyped subjects at a very high level of standardization, as well as a thorough statistical approach accounting for several important covariates.
In conclusion, markers of renal function and lipid metabolism are independent predictors of sectoral cpRNFLT in a large and deeply phenotyped population-based study and should be included as important covariates in future studies on cpRNFLT.
Availability of data and materials
Raw data cannot be shared publicly because of consent restrictions of LIFE-Adult participants. Data are available after an approved project agreement from the LIFE Leipzig Research Center for Civilization Diseases. Please contact Dr. Matthias Nüchter (Head of Managing Office, contact via matthias.nuechter[at]life.uni-leipzig.de) for data access requests.
Lu AT-H, Wang M, Varma R, Schuman JS, Greenfield DS, Smith SD, et al. Combining nerve fiber layer parameters to optimize glaucoma diagnosis with optical coherence tomography. Ophthalmology. 2008;115:1352–1357.e2.
Li D, Rauscher FG, Choi EY, Wang M, Baniasadi N, Wirkner K, et al. Sex-specific differences in circumpapillary retinal nerve fiber layer thickness. Ophthalmology. 2020;127(3):357–68. https://doi.org/10.1016/j.ophtha.2019.09.019.
Mauschitz MM, Bonnemaijer PWM, Diers K, Rauscher FG, Elze T, Engel C, et al. Systemic and ocular determinants of peripapillary retinal nerve fiber layer thickness measurements in the European Eye Epidemiology (E3) population. Ophthalmology. 2018;125(10):1526–36. https://doi.org/10.1016/j.ophtha.2018.03.026.
Lamparter J, Schmidtmann I, Schuster AK, Siouli A, Wasielica-Poslednik J, Mirshahi A, et al. Association of ocular, cardiovascular, morphometric and lifestyle parameters with retinal nerve fibre layer thickness. PLoS ONE. 2018;13(5):e0197682. https://doi.org/10.1371/journal.pone.0197682.
Loeffler M, Engel C, Ahnert P, Alfermann D, Arelin K, Baber R, et al. The LIFE-Adult-study: objectives and design of a population-based cohort study with 10,000 deeply phenotyped adults in Germany. BMC Public Health. 2015;15(1):691. https://doi.org/10.1186/s12889-015-1983-z.
Kang SH, Hong SW, Im SK, Lee SH, Ahn MD. Effect of myopia on the thickness of the retinal nerve fiber layer measured by Cirrus HD optical coherence tomography. Invest Ophthalmol Vis Sci. 2010;51(8):4075–83. https://doi.org/10.1167/iovs.09-4737.
Hirasawa K, Shoji N, Yoshii Y, Haraguchi S. Determination of axial length requiring adjustment of measured circumpapillary retinal nerve fiber layer thickness for ocular magnification. PLoS ONE. 2014. https://doi.org/10.1371/journal.pone.0107553
Wang M, Elze T, Li D, Baniasadi N, Wirkner K, Kirsten T, et al. Age, ocular magnification, and circumpapillary retinal nerve fiber layer thickness. J Biomed Opt. 2017;22:121718.
Madrigal-González J, Calatayud J, Ballesteros-Cánovas JA, Escudero A, Cayuela L, Rueda M, et al. Climate reverses directionality in the richness–abundance relationship across the world’s main forest biomes. Nat Commun. 2020;11(1):5635. https://doi.org/10.1038/s41467-020-19460-y.
Schuman JS, Hee MR, Puliafito CA, Wong C, Pedut-Kloizman T, Lin CP, et al. Quantification of nerve fiber layer thickness in normal and glaucomatous eyes using optical coherence tomography: a pilot study. Arch Ophthalmol. 1995;113(5):586–96. https://doi.org/10.1001/archopht.1995.01100050054031.
Liu S, Wang W, Tan Y, He M, Wang L, Li Y, et al. Correlation between renal function and peripapillary choroidal thickness in treatment-naïve diabetic eyes using swept-source optical coherence tomography. Curr Eye Res. 2020;45(12):1526–33. https://doi.org/10.1080/02713683.2020.1753213.
Choi JA, Ko S-H, Park YR, Jee D-H, Ko S-H, Park CK. Retinal nerve fiber layer loss is associated with urinary albumin excretion in patients with type 2 diabetes. Ophthalmology. 2015;122(5):976–81. https://doi.org/10.1016/j.ophtha.2015.01.001.
Fox CS, Matsushita K, Woodward M, Bilo HJ, Chalmers J, Heerspink HJL, et al. Associations of kidney disease measures with mortality and end-stage renal disease in individuals with and without diabetes: a meta-analysis. Lancet. 2012;380(9854):1662–73. https://doi.org/10.1016/S0140-6736(12)61350-6.
Vogels SCM, Emmelot-Vonk MH, Verhaar HJJ, Koek H (Dineke) L. The association of chronic kidney disease with brain lesions on MRI or CT: a systematic review. Maturitas. 2012;71:331–336, 4, DOI: https://doi.org/10.1016/j.maturitas.2012.01.008.
Prospective Studies Collaboration, Lewington S, Whitlock G, Clarke R, Sherliker P, Emberson J, et al. Blood cholesterol and vascular mortality by age, sex, and blood pressure: a meta-analysis of individual data from 61 prospective studies with 55,000 vascular deaths. Lancet Lond Engl. 2007;370:1829–39.
Ference BA, Ginsberg HN, Graham I, Ray KK, Packard CJ, Bruckert E, et al. Low-density lipoproteins cause atherosclerotic cardiovascular disease. 1. Evidence from genetic, epidemiologic, and clinical studies. A consensus statement from the European Atherosclerosis Society Consensus Panel. Eur Heart J. 2017;38(32):2459–72. https://doi.org/10.1093/eurheartj/ehx144.
Emerging Risk Factors Collaboration, Di Angelantonio E, Sarwar N, Perry P, Kaptoge S, Ray KK, et al. Major lipids, apolipoproteins, and risk of vascular disease. JAMA. 2009;302(18):1993–2000. https://doi.org/10.1001/jama.2009.1619.
Tserentsoodol N, Sztein J, Campos M, Gordiyenko NV, Fariss RN, Lee JW, et al. Uptake of cholesterol by the retina occurs primarily via a low density lipoprotein receptor-mediated process. Mol Vis. 2006;12:1306–18.
Zheng W, Reem RE, Omarova S, Huang S, DiPatre PL, Charvet CD, et al. Spatial distribution of the pathways of cholesterol homeostasis in human retina. PLoS ONE. 2012;7(5):e37926. https://doi.org/10.1371/journal.pone.0037926.
Mast N, Bederman IR, Pikuleva IA. Retinal cholesterol content is reduced in simvastatin-treated mice due to inhibited local biosynthesis albeit increased uptake of serum cholesterol. Drug Metab Dispos. 2018;46(11):1528–37. https://doi.org/10.1124/dmd.118.083345.
Zheng W, Mast N, Saadane A, Pikuleva IA. Pathways of cholesterol homeostasis in mouse retina responsive to dietary and pharmacologic treatments [S]. J Lipid Res. 2015;56(1):81–97. https://doi.org/10.1194/jlr.M053439.
Curcio CA, Johnson M, Huang J-D, Rudolf M. Apolipoprotein B-containing lipoproteins in retinal aging and age-related macular degeneration. J Lipid Res. 2010;51(3):451–67. https://doi.org/10.1194/jlr.R002238.
Sohn EH, van Dijk HW, Jiao C, Kok PHB, Jeong W, Demirkaya N, et al. Retinal neurodegeneration may precede microvascular changes characteristic of diabetic retinopathy in diabetes mellitus. Proc Natl Acad Sci. 2016;113(19):E2655–64. https://doi.org/10.1073/pnas.1522014113.
Clerck EEBD, Schouten JSAG, Berendschot TTJM, Kessels AGH, Nuijts RMMA, Beckers HJM, et al. New ophthalmologic imaging techniques for detection and monitoring of neurodegenerative changes in diabetes: a systematic review. Lancet Diabetes Endocrinol. 2015;3(8):653–63. https://doi.org/10.1016/S2213-8587(15)00136-9.
Chatziralli I, Karamaounas A, Dimitriou E, Kazantzis D, Theodossiadis G, Kozobolis V, et al. Peripapillary retinal nerve fiber layer changes in patients with diabetes mellitus: a case-control study. Semin Ophthalmol. 2020;35(4):257–60. https://doi.org/10.1080/08820538.2020.1810289.
The authors wish to thank the LIFE-Adult Study participants for their time, and furthermore, we gratefully acknowledge the LIFE-Adult Study team for their commitment to the eye investigation and corresponding exams to make this analysis possible.
This research was supported by the LIFE Leipzig Research Center for Civilization Diseases, Leipzig University (LIFE is funded by the EU, the European Social Fund, the European Regional Development Fund, and Free State Saxony’s excellence initiative; project numbers: 713-241202, 14505/2470, 14575/2470); Lions Foundation; Grimshaw-Gudewicz Foundation; Research to Prevent Blindness; BrightFocus Foundation; Alice Adler Fellowship; NIH K99EY028631 to MW; NEI Core Grant P30EYE003790; NIH R21 EY030142; R21EY030631; R01EY030575; German Federal Ministry of Education and Research i:DSem - Integrative data semantics in systems medicine (031L0026) to FGR. ThE was supported by a Novo Nordisk postdoctoral fellowship run in partnership with Karolinska Institutet, Stockholm, Sweden, a Karolinska Institutet Research Foundation grant, as well as by the Swedish Kidney Foundation (Njurfonden) and by the Stiftelsen Stig och Gunborg Westman. ThE was further funded through the EFSD Mentorship Programme supported by AstraZeneca. PS and ThE benefited from the support from the Strategic Research Program in Diabetes at Karolinska Institutet (Swedish Research Council grant No 2009-1068).
The sponsor or funding organization had no role in the design or conduct of this research. Open Access funding provided by Karolinska Institute.
Franziska G. Rauscher, Mengyu Wang, Tobias Elze and Thomas Ebert contributed equally to this work.
Authors and Affiliations
Leipzig Research Centre for Civilization Diseases (LIFE), Leipzig University, Leipzig, Germany
Franziska G. Rauscher, Mengyu Wang, Kerstin Wirkner, Christoph Engel, Joachim Thiery, Markus Loeffler & Tobias Elze
Institute for Medical Informatics, Statistics, and Epidemiology (IMISE), Leipzig University, Leipzig, Germany
Franziska G. Rauscher, Mike Francke, Kerstin Wirkner, Christoph Engel & Markus Loeffler
Schepens Eye Research Institute, Harvard Medical School, Boston, MA, USA
Mengyu Wang & Tobias Elze
Medical Department III – Endocrinology, Nephrology, Rheumatology, University of Leipzig Medical Center, Leipzig, Germany
Anke Tönjes, Michael Stumvoll & Thomas Ebert
Institute of Laboratory Medicine, Clinical Chemistry and Molecular Diagnostics, Leipzig University, Leipzig, Germany
Department of Clinical Science, Intervention and Technology, Division of Renal Medicine, Karolinska Institutet, Stockholm, Sweden
FGR, MW, TE, and ThE designed the study, analyzed data, and wrote the manuscript. MF, KW, AT, CE, JT, PS, MS, and ML analyzed the data and contributed to the interpretation of the data. All authors reviewed, edited, and approved the final manuscript.
The study was approved by the Ethical Committee at the Medical Faculty of Leipzig University (approval number: 263-2009-14122009) and adheres to the Declaration of Helsinki and all federal and state laws. Prior to inclusion, informed written consent was obtained from all participants.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Sectoral multivariable linear regression analyses for cardiometabolic biomarkers and cpRNFLT in all subjects with further adjustment for SD-OCT-derived clinical and subclinical ONH abnormalities (including glaucoma and other ONH diseases), as well as patient-reported glaucoma diagnosis and medication (N = 8952).
Sectoral analyses derived from cpRNFLT in patients with optic nerve head abnormalities, self-reported glaucoma diagnosis, and/or glaucoma medication (N = 1,180).
Rights and permissions
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/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
Rauscher, F.G., Wang, M., Francke, M. et al. Renal function and lipid metabolism are major predictors of circumpapillary retinal nerve fiber layer thickness—the LIFE-Adult Study.
BMC Med19, 202 (2021). https://doi.org/10.1186/s12916-021-02064-8