Skip to main content

ORIGINAL RESEARCH article

Front. Cell Dev. Biol., 27 August 2024
Sec. Molecular and Cellular Pathology
This article is part of the Research Topic Artificial Intelligence Applications in Chronic Ocular Diseases, Volume II View all 3 articles

Unveiling the hidden: a deep learning approach to unraveling subzone-specific changes in peripapillary atrophy in type 2 diabetes

Yingying Li,,,&#x;Yingying Li1,2,3,4Xinxin Hu,,,&#x;Xinxin Hu1,2,3,4Xinyu GuoXinyu Guo5Xueqiong Ye,,,Xueqiong Ye1,2,3,4Dandan Wang,,,Dandan Wang1,2,3,4Juntao Zhang,,,Juntao Zhang1,2,3,4Weina Ren,,,Weina Ren1,2,3,4Na Zhao,,,
Na Zhao1,2,3,4*Yitian Zhao
Yitian Zhao5*Qinkang Lu,,,
Qinkang Lu1,2,3,4*
  • 1Department of Ophthalmology, The Affiliated People’s Hospital of Ningbo University, Ningbo University, Ningbo, China
  • 2Ningbo Clinical Research Center for Ophthalmology, Ningbo, China
  • 3Ningbo Key Laboratory for Neuroretinopathy Medical Research, Ningbo, China
  • 4Eye Hospital of Wenzhou Medical University (Ningbo Branch), Ningbo, China
  • 5Laboratory of Advanced Theranostic Materials and Technology, Ningbo Institute of Materials Technology and Engineering, Chinese Academy of Sciences, Ningbo, China

Purpose: This study aimed to evaluate the optical coherence tomography angiography (OCTA) changes in subzones of peripapillary atrophy (PPA) among type 2 diabetic patients (T2DM) with or without diabetic retinopathy (DR) using well-designed deep learning models.

Methods: A multi-task joint deep-learning model was trained and validated on 2,820 images to automate the determination and quantification of the microstructure and corresponding microcirculation of beta zone and gamma zone PPA. This model was then applied in the cross-sectional study encompassing 44 eyes affected by non-proliferative diabetic retinopathy (NPDR) and 46 eyes without DR (NDR). OCTA was utilized to image the peripapillary area in four layers: superficial capillary plexus (SCP), deep capillary plexus (DCP), choroidal capillary (CC) and middle-to-large choroidal vessel (MLCV).

Results: The patients in both groups were matched for age, sex, BMI, and axial length. The width and area of the gamma zone were significantly smaller in NPDR group compared to the NDR group. Multiple linear regression analysis revealed a negative association between the diagnosis of DR and the width and area of the gamma zone. The gamma zone exhibited higher SCP, DCP and MLCV density than the beta zone, while the beta zone showed higher CC density than the gamma zone. In comparison to the NDR group, the MLCV density of gamma zone was significantly lower in NPDR group, and this density was positively correlated with the width and area of the gamma zone.

Discussion: DR-induced peripapillary vascular changes primarily occur in gamma zone PPA. After eliminating the influence of axial length, our study demonstrated a negative correlation between DR and the gamma zone PPA. Longitudinal studies are required to further elucidate the role of the gamma zone in the development and progression of DR.

1 Introduction

The global incidence of diabetes has been increasing in recent years. According to the Diabetes Atlas 10th edition, the prevalence of diabetes is expected to reach 783 million by 2045 (Sun et al., 2023). Diabetic retinopathy (DR) is one of the major microvascular complications in diabetic patients and has become a leading cause of irreversible visual impairment and blindness among the adult working population and elderly individuals worldwide (Yau et al., 2012). Recent studies have suggested that retinal nerve fiber layer and retinal ganglion cells are damaged at an early stage of DR, indicating that DR is a vascular neuropathy (Chhablani et al., 2015; Vujosevic et al., 2018).

Although several studies have demonstrated diabetic macular vascular changes (Sun et al., 2019; You et al., 2020), the evidence regarding the associations between peripapillary vascular changes and the pathogenesis of DR has been limited. Previous studies found a negative correlation between the presence of peripapillary atrophy (PPA)-β and the prevalence of DR, independent of axial length (AL) and refractive status (Tan et al., 2018). However, that study did not quantitatively determine the microstructure and microcirculation of PPA-β, which limited further exploration of the association between PPA-β and DR. Additionally, PPA-β has been found to be associated with the incidence and progression of glaucoma (Jonas and Naumann, 1989; Lee et al., 2011; Uchida et al., 1998). Considering that diabetes is a recognized high-risk factor for glaucoma, and altered ocular microcirculation is one of the common causes in both diseases, it is reasonable to speculate that PPA-β and corresponding vascular metrics are associated with DR.

PPA-β, characterized by visible large choroidal vessels and sclera, has recently been subdivided into a new beta zone and a gamma zone based on the Bruch’s membrane opening (BMO) (Jonas et al., 2012; Dai et al., 2013). The gamma zone represents the peripapillary area free of Bruch’s membrane (BM), while the new beta zone is defined as the continued presence of BM and absence of retinal pigment epithelium (RPE). Previous studies have reported the microvascular features of PPA subzones in myopic eyes using spectral domain optical coherence tomography angiography (SD-OCTA) (Hu et al., 2021). To our knowledge, there have been no studies exploring the changes in PPA subzones secondary to DR. In this study, we introduce a multi-task joint deep-learning framework based on state-of-the-art models to achieve automated determination and quantification of the microstructure and microcirculation of the beta zone and gamma zone PPA in diabetes with or without DR.

2 Methods

2.1 Study design and participants

In this cross-sectional study, we adhered to the principles outlined in the Declaration of Helsinki and obtained approval from Ethics Committee of the Ophthalmology Center of The Affiliated People’s Hospital of Ningbo University, China. Written informed consent was obtained from all participants.

The inclusion criteria for all subjects were as follows: 1) individuals with type 2 diabetes mellitus (T2DM) without diabetic retinopathy (NDR) or with untreated non-proliferative diabetic retinopathy (NPDR). The diagnosis of T2DM followed WHO standards (Alberti and Zimmet, 1998), while the grading of NPDR was based on the International Clinical DR Disease Severity Scale (Wilkinson et al., 2003). 2) Presence of PPA (beta zone, gamma zone, or both).

The exclusion criteria were as follows: 1) intraocular pressure (IOP) > 21 mmHg; 2) presence of other chorioretinal diseases, glaucoma, optic neuropathy, or systemic diseases that may affect ocular perfusion; 3) history of intraocular surgery including cataract surgery, retina laser and intravitreal injections. 4) presence of mental illness or inability to cooperate with the examination; 5) poor quality of OCTA images and structural OCT images. Once the binocular parameters met the above criteria, one eye per participant was randomly selected.

2.2 Data acquisition

A detailed questionnaire collected personal information including age, gender, medical history, and history of eye disease/surgery. General examinations measured height, weight, systolic blood pressure (SBP) and diastolic blood pressure (DBP) for each participant. DR was diagnosed and graded using digital fundus photography (Canon CR-2; Tokyo, Japan). The fundus camera captured two images, one centered on the macula and the other on the optic disc. Axial length and IOP were measured using the IOLMaster (Carl Zeiss Meditec AG, Jena, Germany), and the noncontact tonometer (Auto Tonometer TX-F; Topcon, Tokyo, Japan), respectively. To calculate the patients’ body mass index (BMI), their weight was divided by their height in meters squared (kg/m2). The mean arterial pressure (MAP) and ocular perfusion pressure (OPP) were calculated using the following formulas: MAP = 1/3SBP+2/3DBP and OPP = 2/3MAP-IOP (Longo et al., 2004).

2.3 OCT angiography image acquisition

For the acquisition of peripapillary vascular imaging, we utilized SS-OCT angiography (SS-OCTA) with the BM-400K BMizar from TowardPi Medical Technology Co., Ltd. in Beijing, China. The imaging camera employed a vertical cavity laser with a wavelength of 1,060 nm and a scanning speed of 400 KHz A-scan/s. This setup allowed for a lateral resolution of 10 μm and an axial resolution of 3.8 μm (optical), enabling the visualization of en face images of the retinal and choroidal vascular perfusion system. To evaluate the peripapillary area of each subject, B-scans were performed, covering a 6 mm2 × 6 mm2 area horizontally and vertically.

To stratify the retinal and choroidal vasculature, we utilized the built-in software in the OCTA system. This software allowed us to segment the vasculature into several sub-layers, including the superficial capillary plexus (SCP), deep capillary plexus (DCP), choroidal capillary (CC), and middle-to-large choroidal vessel (MLCV). In our study, we excluded images with significant motion artifacts or those graded below 8. Manual correction was only performed when automatic stratification was inaccurate.

All OCTA images were acquired by an experienced operator (YL) and reviewed by two other authors (XH and NZ).

2.4 Multi-task joint deep-learning framework

2.4.1 PPA segmentation

The beta zone PPA is defined as the presence of BM without the RPE. On the other hand, the gamma zone PPA is characterized by the absence of BM and is located between the optic papilla border and beta zone. To accurately delineate the boundaries of the optic disc, beta zone, and gamma zone PPA, both en face and serial B-scan OCT images of the optic disc were acquired at the same position, as shown in Figure 1. For the segmentation of the optic papilla border and RPE atrophy, we employed an advanced image segmentation algorithm called nnU-Net (Isensee et al., 2021). We chose the Stochastic Gradient Descent (SGD) with Nesterov momentum of 0.99 as optimizer and the initial learning rate was set to 0.001 in our experiment. The total epochs reached 1,000, with each epoch consisting of 250 iterations. To train and validate this model, we utilized an additional dataset of 200 en face images that were manually annotated. The training and validation sets were partitioned in a 7:3 ration, and the segmentation performance demonstrated its robustness and effectiveness, with a Dice coefficient greater than 0.95.

Figure 1
www.frontiersin.org

Figure 1. The workflow of the beta zone and gamma zone PPA automated detection. Firstly, the nnU-Net was used for segmenting the border of optic papilla and RPE atrophy. Then the YoLo-V8 was employed to detect the Bruch’s membrane opening (BMO) on B-scan series. Lastly, the BM boundary curve was obtained by reconstructing the BMO points from B-Scans onto enface images and fitting them with a polynomial function.

2.4.2 BMO detection

To accurately detect the BMO on B-scan series, we utilized YoLo-V8 (Wang et al., 2024), an end-to-end object detection model that excels in real-time processing of high-resolution images. This model has been widely applied in the field of medical image detection, enabling precise target identification. Prior to implementation, YoLo-V8 was trained and validated on 20 subjects, consisting of 2,500 manually annotated B-scan images. We chose the Stochastic Gradient Descent (SGD) with Nesterov momentum of 0.937 as optimizer and the initial learning rate was set to 0.001. The maximum number of epochs reached 300. The BM boundary curve was obtained by reconstructing the BMO points from the B-Scans onto enface images and fitting them with a polynomial function. Subsequently, the width and area of PPA subzone were calculated using MATLAB software.

2.4.3 Vessel segmentation and quantification

We utilized an enhanced U-Net model for the segmentation of vessel in the SCP, DCP, CC and MLCV of OCTA images (Figure 2). This model incorporated a U-Net (Ronneberger et al., 2015) as the underlying architecture, with the original encoder was replaced by a pre-trained ResNet18 (He et al., 2016). We chose the AdamW as optimizer and the initial learning rate was set to 0.0002. The maximum number of epochs reached 200. To train the network, we employed a proprietary OCTA vessel segmentation dataset called OCTA-Z. OCTA-Z comprising 126 OCTA en face angiograms from 42 subjects, with a resolution of 512 × 512 pixels. The vessel annotation at the pixel-level was performed by four proficient image experts, followed by a review and refinement of the initial labeling by two senior ophthalmologists. For this study, we allocated 120 OCTA images from 30 subjects for training the improved U-Net, while the remaining images were used for testing. The network achieved a Dice coefficient of 81.86% and an accuracy is 90.13% on the testing set. Vascular density (%) was determined as the proportion of the measured area occupied by flowing blood vessels within a specific region. We measured the vascular density of four retinal and choroidal sub-layers in beta zone and gamma zone PPA, respectively.

Figure 2
www.frontiersin.org

Figure 2. Illustration of the SS-OCTA imaging and vascular segmentation of peripapillary vasculature parameters in each retina and choroidal layer. (A) Superficial capillary plexus (SCP). (B) Deep capillary plexus (DCP). (C) Choroidal capillary (CC). (D) Middle-to-large choroidal vessel (MLCV). (A1–D1) SS-OCTA images of participant without DR. (A3–D3) SS-OCTA images of participant with NPDR. (A2–D2, A4–D4) The corresponding automated segmented vessel map of A1–D1 and A3–D3, respectively.

2.5 Statistical analysis

All statistical analyses were performed using SPSS version 29.0 (SPSS, Chicago, IL, United States). Continuous data were presented as mean ± standard deviation, while categorical data were reported as frequency and percentage (%). Parametric data were analyzed using Student’s t-test, while non-parametric data were analyzed using the Mann–Whitney U test was for group comparisons. The Chi-squared test was employed to compare categorical variables. The gamma zone indicators that exhibited significant differences between the two groups were selected as dependent variables, and stepwise linear regression analysis was conducted to identify independent risk factors. Additionally, partial correlation coefficient analyses were carried out to assess the effects of PPA subzone parameters on ocular parameters after adjusting for axial length. All P-values were two-sided, and statistical significance was determined when the P-values were <0.05.

3 Results

The study consisted of 46 eyes from 46 subjects without DR (NDR group) and 44 eyes from 44 subjects with NPDR (NPDR group). Among the participants, 54 (60.0%) were male and 36 (40.0%) were female. The mean age was 60.3 ± 9.8 years (median, 60 years; range, 33–85 years), and the mean axial length was 24.1 ± 1.1 mm (median, 23.9 mm; range, 22.3–27.4 mm). All subjects had a diagnosis of T2DM, and their eyes were phakic.

3.1 Comparison of demographic and clinical characteristics between the groups

Table 1 provides a summary of the comparison of demographic and clinical characteristics between the NDR group and NPDR group. The two groups were matched for age, sex, BMI and axial length. Significant differences were observed in fasting plasma glucose (FPG), glycated hemoglobin (HbA1c), and urine albumin-to-creatinine ratio (UACR) between the two groups.

Table 1
www.frontiersin.org

Table 1. Demographic and clinical characteristics of the participants enrolled in this study.

The FPG of the NPDR group (8.1 ± 1.9 mmol/L; range, 4.9–14.6 mmol/L) was significantly higher (P < 0.001) than that of the NDR group (7.0 ± 2.3 mmol/L; range, 4.2–14.7 mmol/L). Similarly, mean HbA1c (7.5% ± 1.5%) and UACR (139.3 ± 275.3 mg/g) in the NPDR group were significantly higher (all P < 0.001) than those in the NDR group (6.5% ± 1.4% and 14.0 ± 18.4 mg/g). No significant differences were observed between the two groups in terms of IOP, OPP, DBP, SBP, MAP, total cholesterol (TC), triglycerides (TG), high-density lipoprotein (HDL), and low-density lipoprotein (LDL) levels (all P > 0.05).

3.2 Comparison of beta and gamma zone characteristics between the groups

In this study, a total of 60 eyes had a beta zone, with 30 (50.0%) eyes in the NDR group and 30 (50.0%) eyes in the NPDR group (P = 0.77). Additionally, 60 eyes had a gamma zone, with 32 (53.3%) eyes in the NDR group and 28 (46.7%) eyes in the NPDR group (P = 0.55) (Table 2). The mean width of the gamma zone was significantly smaller in the NPDR group (0.39 ± 0.14 mm) compared to the NDR group (0.52 ± 0.18 mm) (P < 0.001). Similar statistically significant differences were observed in the gamma zone area, where individuals in the NPDR group exhibited a smaller gamma zone area (0.57 ± 0.17 mm2) compared to those in the NDR group (0.96 ± 0.44 mm2) (P < 0.001). Furthermore, the vascular density of MLCV in the gamma zone was significantly lower in the NPDR group (22.90% ± 10.54%) compared to the NDR group (29.36% ± 8.11%) (P = 0.005).

Table 2
www.frontiersin.org

Table 2. Comparison of beta and gamma zone PPA between the NDR and NPDR groups.

After conducting a multiple forward stepwise linear regression analysis, we examined the relationship between the width and area of the gamma zone PPA and several independent variables, including age, sex, BMI, FPG, HbA1c, UACR, TC, diagnosis of DR, and axial length (Table 3). The findings revealed significant associations between these variables and the width and area of gamma zone. Specifically, the diagnosis of DR was found to have a negative association with both the width and area of the gamma zone (P = 0.005 and P < 0.001, respectively). This suggests that individuals with a diagnosis of DR tend to have narrower and smaller gamma zones. On the other hand, longer axial length was positively associated with the width and area of gamma zone (P = 0.005 and P = 0.007, respectively) (Figure 3). This indicates that individuals with a greater axial length tend to have wider and larger gamma zones. These findings provide valuable insights into the relationships between the diagnosis of DR, axial length, and the characteristics of the gamma zone.

Table 3
www.frontiersin.org

Table 3. Stepwise linear regression analysis of predictors for width and area of gamma zone PPA in all subjects.

Figure 3
www.frontiersin.org

Figure 3. Scatterplot showing the correlation between width and area of gamma zone PPA and axial length.

3.3 Analyses of parapapillary vascular parameters in beta and gamma zone

Table 4 presents the parapapillary vascular differences between beta zone and gamma zone. The SCP density was 31.36% ± 5.55% in the beta zone and 34.29% ± 6.37% in the gamma zone, showing a statistically significant difference (P = 0.003). Similarly, the DCP density in the gamma zone (39.78% ± 6.16%) was significantly higher than that in the beta zone (33.72% ± 5.10%) (P < 0.001). The DCP density of the gamma zone PPA was negatively correlated with the width (P = 0.003) and area (P < 0.001) of the gamma zone (Table 5).

Table 4
www.frontiersin.org

Table 4. Comparison of parapapillary vascular density in beta zone and gamma zone PPA.

Table 5
www.frontiersin.org

Table 5. Correlation with the width and area of beta and gamma zone PPA.

Regarding the choroidal metrics, the beta zone exhibited a higher CC density (50.96% ± 6.93%) compared to the gamma zone (38.72% ± 8.68%), with a statistically significant difference (P < 0.001). The CC density in the beta zone PPA was negatively correlated with the width (P = 0.029) and area (P < 0.001) of the beta zone. However, the MLCV density in the gamma zone (26.35% ± 9.80%) was significantly higher than that in the beta zone (20.34% ± 9.10%) (P < 0.001). Notably, the MLCV density of gamma zone PPA was positively correlated with the width (P = 0.004) and the area (P < 0.001) of the gamma zone.

4 Discussion

This study represents the first investigation into the impact of DR on beta and gamma zone PPA. To achieve automated determination and quantification of the microstructure and microcirculation of the corresponding beta zone and gamma zone PPA, we employed a state-of-the-art deep learning approach. Our findings indicate that patients with DR had a significantly smaller gamma zone PPA compared to those without DR, but no significant difference was observed in the beta zone PPA. Additionally, we observed topographic differences in the SCP, DCP, CC and MLCV between the beta zone and gamma zone PPA. The gamma zone PPA exhibited higher densities of SCP, DCP and MLCV, while the beta zone displayed higher CC density compared to the gamma zone.

The gamma zone PPA is characterized by atrophy of the RPE and BM. Previous studies have suggested a close relationship between the formation and development of the gamma zone PPA and AL (Jonas et al., 2016; Zhang et al., 2018). Our study confirmed a similar result, showing that longer axial length was associated with the width and area of the gamma zone. Furthermore, our multiple forward stepwise linear regression analysis revealed that the presence of DR independently increased the risk factor of a narrower and smaller gamma zone PPA. This aligns with a study by Lin et al. (2022), who reported a significantly smaller PPA-β area in children with diabetes compared to healthy controls. However, their study did not differentiate between the beta zone and gamma zone PPA, limiting the further exploration. Considering that the choroidal circulation serves as the primary source of oxygen and nutrients for the choroid and outer layers of the retina, one may speculate that the effect of DR on the gamma zone PPA could be associated with hyperglycemia-induced choroidal alterations. Previous studies have indicated larger peripapillary choroidal thickness (pCT) in individuals with diabetes compared to healthy individuals (Li et al., 2020). It is possible that the increased supply of peripapillary blood flow during the early stage of DR inhibits atrophic changes to BM, thus slowing down the expansion of the gamma zone PPA.

The beta zone PPA is defined as absence of RPE with continued presence of BM. Previous studies have shown that beta zone PPA is associated with age and the presence of glaucoma (Dai et al., 2013; Kim et al., 2014). Sung et al. (2018) found that myopic eyes with PPA + BM (indicating beta zone PPA in this study) had the lowest superficial radial peripapillary capillary (RPC) and deep vessel density, which were closely negatively correlated with the width of the beta zone. Similarly, we previously reported that the beta zone PPA in myopic eyes had the lowest superficial RPC and choroidal microvasculature compared to the other subzones of PPA, and these parameters were negatively associated with the width and area of the beta zone PPA (Hu et al., 2021). This study is partially consistent with these findings, as we observed lower SCP, DCP and MLCV densities in the beta zone compared to the gamma zone. Both SCP and DCP densities were negatively associated with the width and area of the beta zone PPA, indicating a microcirculatory deficiency in the retinal layer of the beta zone PPA. However, the density of CC measured in our study was significantly lower than that reported in previous studies. This difference could be attributed to the fact that the participants in those studies were normal myopic eyes with a mean age of no more than 30 years. Furthermore, diabetes and its associated ocular complication, DR, are recognized risk factors for increased peripapillary CC flow deficit percentage (CC FD%) (Guo et al., 2023; Wang et al., 2023). It is important to note that the SS-OCTA we utilized in our study showed higher sensitivity and image resolution in distinguishing between CC and MLCV compared to conventional SD-OCTA, which may explain the differences in outcome.

The gamma zone PPA is primarily formed by oblique border tissue in the absence of the BM. The underlying choroidal structure in the gamma zone PPA is distinct from its original status without gamma zone formation (Kim et al., 2018). Our data revealed that the gamma zone had a lower CC density compared to the beta zone PPA, suggesting that the CC in gamma zone would decrease first following the absence of RPE and BM. Additionally, we found that MLCV density in the gamma zone PPA was significantly lower in NPDR group compared to the NDR group. Previous studies have indicated that the peripheral choroid is more susceptible to diabetes-induced injury than the central area (Xu et al., 2023), which could explain the decrease of MLCV density adjacent to the optic disc preceding changes in the macular area. As the gamma zone PPA enlarges, the self-regulatory mechanism of MLCV compensates for the decrease in CC perfusion, maintaining the basic oxygen and nutrient requirements of the retina. This may elucidate the positive correlation observed between MLCV density and the width and area of the gamma zone PPA.

There are several limitations in the present study. Firstly, the participants enrolled in our study primarily consisted of older individuals with T2DM, therefore, our findings may not applicable to children with Type 1 Diabetes Mellitus (T1DM). Secondly, our study excluded healthy subjects without diabetes as well as subjects with proliferative diabetic retinopathy (PDR). Including these two groups in future studies could provide further information. Thirdly, all subjects were from the same center and the sample size was relatively small. Lastly, due to the inherent limitations of the cross-sectional design, we were unable to observe longitudinal changes in the microstructure and microvasculature of the beta and gamma zone PPA in diabetes with or without DR. Therefore, further exploration of the relationship between gamma zone PPA and the development of DR over time should be conducted through larger multi-center longitudinal studies.

In conclusion, our findings demonstrated that DR-induced peripapillary vascular changes primarily occur in the gamma zone PPA. After accounting for AL, our study reveals a negative correlation between DR and gamma zone PPA. Longitudinal studies are necessary to elucidate the role of the gamma zone in the development and progression of DR.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by the Affiliated People’s Hospital of Ningbo University. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

YL: Data curation, Formal Analysis, Methodology, Writing–original draft, Writing–review and editing. XH: Data curation, Formal Analysis, Methodology, Writing–original draft, Writing–review and editing. XG: Data curation, Formal Analysis, Investigation, Software, Writing–original draft, Writing–review and editing. XY: Data curation, Methodology, Writing–original draft, Writing–review and editing. DW: Data curation, Methodology, Writing–original draft, Writing–review and editing. JZ: Methodology, Resources, Validation, Writing–review and editing. WR: Methodology, Resources, Validation, Writing–review and editing. NZ: Conceptualization, Funding acquisition, Supervision, Validation, Writing–original draft, Writing–review and editing. YZ: Conceptualization, Funding acquisition, Supervision, Validation, Writing–original draft, Writing–review and editing. QL: Conceptualization, Funding acquisition, Supervision, Writing–original draft, Writing–review and editing.

Funding

The authors declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by grants from the National Natural Science Foundation of China (No. 82000887), Zhejiang Provincial Natural Science Foundation of China (No. LY20C090001), Ningbo Major Scientific and Technological Research and “Unveiling and Commanding” Project (No. 2021Z054), Zhejiang Medical and Health Science and Technology Plan Project (Nos. 2023RC090, 2024KY1598 and 2023KY301), Ningbo Natural Science Foundation (Nos. 2019A610352, 2022J029 and 2023J385), Ningbo Science and Technology Planning Project in the Public Interest (Nos. 2022S052 and 2022S054), and Ningbo “Innovation Yongjiang 2035” Key Technology Breakthrough program (2024Z233). The authors were supported by Ningbo Science and Technology Project (Nos. 2019C50085 and 2019C50053), Ningbo Medicine and Health Science and Technology Project (No. 2020Y46), Ningbo Clinical Research Center for Ophthalmology (2022L003), Ningbo Key Laboratory for neuroretinopathy medical research, and the Project of NINGBO Leading Medical & Health Discipline (No. 2016-S05).

Acknowledgments

We thank Liwen Home for Researchers (http://www.home-for-researchers.com) for basic language editing of a draft of this manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Alberti, K. G., and Zimmet, P. Z. (1998). Definition, diagnosis and classification of diabetes mellitus and its complications. Part 1: diagnosis and classification of diabetes mellitus provisional report of a WHO consultation. Diabet. Med. 15 (7), 539–553. doi:10.1002/(SICI)1096-9136(199807)15:7<539::AID-DIA668>3.0.CO;2-S

PubMed Abstract | CrossRef Full Text | Google Scholar

Chhablani, J., Sharma, A., Goud, A., Peguda, H. K., Rao, H. L., Begum, V. U., et al. (2015). Neurodegeneration in type 2 diabetes: evidence from spectral-domain optical coherence tomography. Investig. Ophthalmol. Vis. Sci. 56 (11), 6333–6338. doi:10.1167/iovs.15-17334

CrossRef Full Text | Google Scholar

Dai, Y., Jonas, J. B., Huang, H., Wang, M., and Sun, X. (2013). Microstructure of parapapillary atrophy: beta zone and gamma zone. Investig. Ophthalmol. Vis. Sci. 54 (3), 2013–2018. doi:10.1167/iovs.12-11255

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, X., Chen, Y., Bulloch, G., Xiong, K., Chen, Y., Li, Y., et al. (2023). Parapapillary choroidal microvasculature predicts diabetic retinopathy progression and diabetic macular edema development: a three-year prospective study. Am. J. Ophthalmol. 245, 164–173. doi:10.1016/j.ajo.2022.07.008

PubMed Abstract | CrossRef Full Text | Google Scholar

He, K., Zhang, X., Ren, S., and Sun, J. (2016). “Deep residual learning for image recognition,” in 2016 IEEE Conference on Computer Vision and Pattern Recognition, Las Vegas, NV, June 27–30, 2016, 770–778. doi:10.1109/CVPR.2016.90

CrossRef Full Text | Google Scholar

Hu, X., Shang, K., Chen, X., Sun, X., and Dai, Y. (2021). Clinical features of microvasculature in subzones of parapapillary atrophy in myopic eyes: an OCT-angiography study. Eye (Lond). 35 (2), 455–463. doi:10.1038/s41433-020-0872-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Isensee, F., Jaeger, P. F., Kohl, S. A. A., Petersen, J., and Maier-Hein, K. H. (2021). nnU-Net: a self-configuring method for deep learning-based biomedical image segmentation. Nat. Methods 18 (2), 203–211. doi:10.1038/s41592-020-01008-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Jonas, J. B., Jonas, S. B., Jonas, R. A., Holbach, L., Dai, Y., Sun, X., et al. (2012). Parapapillary atrophy: histological gamma zone and delta zone. PLoS One 7 (10), e47237. doi:10.1371/journal.pone.0047237

PubMed Abstract | CrossRef Full Text | Google Scholar

Jonas, J. B., and Naumann, G. O. (1989). Parapapillary chorioretinal atrophy in normal and glaucoma eyes. II. Correlations. Investig. Ophthalmol. Vis. Sci. 30 (5), 919–926.

PubMed Abstract | Google Scholar

Jonas, J. B., Wang, Y. X., Zhang, Q., Fan, Y. Y., Xu, L., Wei, W. B., et al. (2016). Parapapillary gamma zone and axial elongation-associated optic disc rotation: the Beijing eye study. Investig. Ophthalmol. Vis. Sci. 57 (2), 396–402. doi:10.1167/iovs.15-18263

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, M., Choung, H. K., Lee, K. M., Oh, S., and Kim, S. H. (2018). Longitudinal changes of optic nerve head and peripapillary structure during childhood myopia progression on OCT: boramae myopia cohort study report 1. Ophthalmology 125 (8), 1215–1223. doi:10.1016/j.ophtha.2018.01.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, Y. W., Lee, E. J., Kim, T. W., Kim, M., and Kim, H. (2014). Microstructure of β-zone parapapillary atrophy and rate of retinal nerve fiber layer thinning in primary open-angle glaucoma. Ophthalmology 121 (7), 1341–1349. doi:10.1016/j.ophtha.2014.01.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, E. J., Kim, T. W., Weinreb, R. N., Park, K. H., Kim, S. H., and Kim, D. M. (2011). β-Zone parapapillary atrophy and the rate of retinal nerve fiber layer thinning in glaucoma. Investig. Ophthalmol. Vis. Sci. 52 (7), 4422–4427. doi:10.1167/iovs.10-6818

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, T., Jia, Y., Wang, S., Xu, Y., Yin, Y., Wang, A., et al. (2020). Change in peripapillary and macular choroidal thickness change in children with type 1 diabetes mellitus without visual impairment or diabetic retinopathy. Acta Ophthalmol. 98 (2), e203–e211. doi:10.1111/aos.14225

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, Q., Jia, Y., Li, T., Wang, S., Xu, X., Xu, Y., et al. (2022). Optic disc morphology and peripapillary atrophic changes in diabetic children and adults without diabetic retinopathy or visual impairment. Acta Ophthalmol. 100 (1), e157–e166. doi:10.1111/aos.14885

PubMed Abstract | CrossRef Full Text | Google Scholar

Longo, A., Geiser, M. H., and Riva, C. E. (2004). Posture changes and subfoveal choroidal blood flow. Investig. Ophthalmol. Vis. Sci. 45 (2), 546–551. doi:10.1167/iovs.03-0757

PubMed Abstract | CrossRef Full Text | Google Scholar

Ronneberger, O., Fischer, P., and Brox, T. (2015). “U-net: convolutional networks for biomedical image segmentation,” in Medical Image Computing and Computer-Assisted Intervention—MICCAI 2015: 18th International Conference, Munich, Germany, October 5–9, 2015, proceedings, part III 18 (Berlin: Springer), 234–241. doi:10.1007/978-3-319-24574-4_28

CrossRef Full Text | Google Scholar

Sun, H., Saeedi, P., Karuranga, S., Pinkepank, M., Ogurtsova, K., Duncan, B. B., et al. (2023). IDF Diabetes Atlas: global, regional and country-level diabetes prevalence estimates for 2021 and projections for 2045. Diabetes Res. Clin. Pract. 204, 110945. doi:10.1016/j.diabres.2021.109119

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Z., Tang, F., Wong, R., Lok, J., Szeto, S. K. H., Chan, J. C. K., et al. (2019). OCT angiography metrics predict progression of diabetic retinopathy and development of diabetic macular edema: a prospective study. Ophthalmology 126 (12), 1675–1684. doi:10.1016/j.ophtha.2019.06.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Sung, M. S., Heo, H., and Park, S. W. (2018). Microstructure of parapapillary atrophy is associated with parapapillary microvasculature in myopic eyes. Am. J. Ophthalmol. 192, 157–168. doi:10.1016/j.ajo.2018.05.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, N. Y. Q., Tham, Y. C., Ding, Y., Yasuda, M., Sabanayagam, C., Saw, S. M., et al. (2018). Associations of peripapillary atrophy and fundus tessellation with diabetic retinopathy. Ophthalmol. Retina 2 (6), 574–581. doi:10.1016/j.oret.2017.09.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Uchida, H., Ugurlu, S., and Caprioli, J. (1998). Increasing peripapillary atrophy is associated with progressive glaucoma. Ophthalmology 105 (8), 1541–1545. doi:10.1016/S0161-6420(98)98044-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Vujosevic, S., Muraca, A., Gatti, V., Masoero, L., Brambilla, M., Cannillo, B., et al. (2018). Peripapillary microvascular and neural changes in diabetes mellitus: an OCT-angiography study. Investig. Ophthalmol. Vis. Sci. 59 (12), 5074–5081. doi:10.1167/iovs.18-24891

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W., Cheng, W., Yang, S., Chen, Y., Zhu, Z., and Huang, W. (2023). Choriocapillaris flow deficit and the risk of referable diabetic retinopathy: a longitudinal SS-OCTA study. Br. J. Ophthalmol. 107 (9), 1319–1323. doi:10.1136/bjophthalmol-2021-320704

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z., Hua, Z., Wen, Y., Zhang, S., Xu, X., and Song, H. (2024). E-YOLO: recognition of estrus cow based on improved YOLOv8n model. Expert Syst. Appl. 238, 122212. doi:10.1016/j.eswa.2023.122212

CrossRef Full Text | Google Scholar

Wilkinson, C. P., Ferris, F. L., Klein, R. E., Lee, P. P., Agardh, C. D., Davis, M., et al. (2003). Proposed international clinical diabetic retinopathy and diabetic macular edema disease severity scales. Ophthalmology 110 (9), 1677–1682. doi:10.1016/S0161-6420(03)00475-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, F., Li, Z., Yang, X., Gao, Y., Li, Z., Li, G., et al. (2023). Assessment of choroidal structural changes in patients with pre- and early-stage clinical diabetic retinopathy using wide-field SS-OCTA. Front. Endocrinol. (Lausanne) 13, 1036625. doi:10.3389/fendo.2022.1036625

PubMed Abstract | CrossRef Full Text | Google Scholar

Yau, J. W., Rogers, S. L., Kawasaki, R., Lamoureux, E. L., Kowalski, J. W., Bek, T., et al. (2012). Global prevalence and major risk factors of diabetic retinopathy. Diabetes Care 35 (3), 556–564. doi:10.2337/dc11-1909

PubMed Abstract | CrossRef Full Text | Google Scholar

You, Q. S., Wang, J., Guo, Y., Pi, S., Flaxel, C. J., Bailey, S. T., et al. (2020). Optical coherence tomography angiography avascular area association with 1-year treatment requirement and disease progression in diabetic retinopathy. Am. J. Ophthalmol. 217, 268–277. doi:10.1016/j.ajo.2020.04.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Q., Wang, Y. X., Wei, W. B., Xu, L., and Jonas, J. B. (2018). Parapapillary beta zone and gamma zone in a healthy population: the Beijing eye study 2011. Investig. Ophthalmol. Vis. Sci. 59 (8), 3320–3329. doi:10.1167/iovs.18-24141

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: deep learning, peripapillary atrophy, diabetic retinopathy, gamma zone, beta zone

Citation: Li Y, Hu X, Guo X, Ye X, Wang D, Zhang J, Ren W, Zhao N, Zhao Y and Lu Q (2024) Unveiling the hidden: a deep learning approach to unraveling subzone-specific changes in peripapillary atrophy in type 2 diabetes. Front. Cell Dev. Biol. 12:1459040. doi: 10.3389/fcell.2024.1459040

Received: 03 July 2024; Accepted: 16 August 2024;
Published: 27 August 2024.

Edited by:

Huihui Fang, Nanyang Technological University, Singapore

Reviewed by:

Kai Jin, Zhejiang University, China
Yi Zhou, Singapore Eye Research Institute (SERI), Singapore
Xujiao Zhou, Fudan University, China
Meng Wang, Soochow University, China

Copyright © 2024 Li, Hu, Guo, Ye, Wang, Zhang, Ren, Zhao, Zhao and Lu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Na Zhao, rmzhaona@nbu.edu.cn; Yitian Zhao, yitian.zhao@nimte.ac.cn; Qinkang Lu, lqktyyx@163.com

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.