Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 13 October 2025

Sec. Biomechanics

Volume 13 - 2025 | https://doi.org/10.3389/fbioe.2025.1609639

Retinal features as predictive indicators for high myopia: insights from explainable multi-machine learning models

Haohan Zou,,&#x;Haohan Zou1,2,3Jing Liu&#x;Jing Liu4Shenda ShiShenda Shi5Saiguang LingSaiguang Ling6Qian Fan,Qian Fan1,2Yan Huo,,Yan Huo1,2,3Zhou DongZhou Dong6Guoge Han,Guoge Han1,2Shengjin Wang,Shengjin Wang5,7Yan Wang,,
Yan Wang1,2,3*
  • 1Tianjin Eye Hospital, Tianjin, China
  • 2Tianjin Key Lab of Ophthalmology and Visual Science, Tianjin Eye Institute, Tianjin, China
  • 3Nankai University Eye Institute, Nankai University, Tianjin, China
  • 4Department of Ophthalmology, Beijing Friendship Hospital, Capital Medical University, Beijing, China
  • 5HuaHui Jian AI Tech Ltd., Tianjin, China
  • 6EVision technology (Beijing) Co., Ltd., Beijing, China
  • 7Department of Electronic Engineering, Tsinghua University, Beijing, China

Objectives: To investigate the role of retinal characteristics for high myopia (HM) prediction based on multi-machine learning (ML) and to provide an interpretable framework for the results.

Methods: A total of 2981 patients (2981 eyes) were included, comprising 1191 HM eyes and 1790 non-HM eyes. A deep semantic segmentation network was used to quantify retinal structural parameters. Five ML algorithms were evaluated to develop predictive models, and SHapley Additive exPlanations method was applied to analyze feature contribution to the outcomes.

Results: The eXtreme Gradient Boosting achieved an accuracy of 0.81 (95% confidence interval [CI] 0.78–0.85), and an area under the receiver operating characteristic curve of 0.87 (95% CI 0.84–0.89), outperforming other models. The 12 most important factors affecting prediction included tessellation density, seven vascular parameters, two parapapillary atrophy parameters, and two optic disc parameters. The tessellated density >0.025, width of parapapillary atrophy >400 um, parapapillary atrophy area >0.60 × 106 um2 were associated with an increased risk of HM. The mean curvature of the arteries >0.00063, diameter of vessels >55.2 um, curvature of the veins >0.00128, vertical diameter of the optic disc >1320 um, diameter of veins >58.5 um, and diameter of artery >47.1 um was associated with a decreased risk of HM.

Conclusion: The XGBoost model outperformed other algorithms, and SHAP-derived cut-off values for critical risk factors enhanced clinical interpretability.

1 Introduction

The global prevalence of high myopia (HM) has risen significantly, drawing attention due to its association with increased risks of blindness and substantial social burdens (Baird et al., 2020; Xu et al., 2021; Shinojima et al., 2022). Patients with HM are more than 50 percent more likely to develop cataracts, glaucoma, macular degeneration, retinal detachment, and optic neuropathy (Modjtahedi et al., 2021; Verhoeven et al., 2015). These risks increase further with age, refractive error, and axial length (Haarman et al., 2020). Thus, early detection and diagnosis of HM is crucial.

HM often results in alterations to the fundus structures, which may not be evident in the early stages (Salih et al., 2023). As a non-invasive and straightforward diagnostic tool, retinal imaging provides comprehensive information about anatomical structures and optical characteristics (Li et al., 2021; Chen et al., 2024). However, traditional image analysis methods, such as meta-analyses of pathologic myopia (META-PM) and the atrophy, traction, and neovascularization classification systems (Ohno-Matsui et al., 2016; Ruiz-Medrano et al., 2019), are commonly dependent on the knowledge of ophthalmologists. These methods frequently fail to provide rapid and precise outcomes, and the requisite expertise limits their broader applicability. Hence, automatic identification and quantification of key features, and accurate data interpretation is vital for the diagnosis and prediction of disease progression (Conti et al., 2021). Recent advancements in computer technology and machine learning (ML) algorithms have made radiomics invaluable in clinical decision-making (Sun R. et al., 2018; Sun et al., 2021; Yu et al., 2021). By combining retinal imaging with omics, radiomics can automatically extract data on fundus shape, texture, geometric characteristics, and relative positions (Zhang et al., 2025). This approach enables rapid identification of HM features and provides valuable insights.

This study investigated the role of fundus features in predicting HM using a retinal imaging omics-based analytical method. Additionally, SHapley Additive exPlantations (SHAP) was employed to help ophthalmologists better understand the relationships between these features and HM, thereby supporting accurate screening and prediction of high-risk groups.

2 Methods

2.1 Study approval

This study was registered with the Chinese Clinical Trial Register (ChiCTR2100049885) and approved by the Ethics Committee of Tianjin Eye Hospital (TJYYLL2021018) in accordance with the Declaration of Helsinki. The study followed the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines (Von Elm et al., 2007). The ethics committee waived the requirement for informed consent because the data were de-identified before the study to ensure patient privacy.

2.2 Patients and data

The examination records of 4,048 patients from August 2019 to December 2021 were retrospectively collected using an electronic medical record system. Of these, 53.9% (n = 2,183) were male and 46.1% (n = 1,865) were female. The mean age was 22.72 ± 5.24 (range 16–45) years. The inclusion criteria included age ≥16 years, diagnosis of refractive error, and a best-corrected visual acuity of at least 0.0 (LogMAR). Patients with disease that affect vision acuity, such as cataract, various types of glaucoma, retinal tears, retinal detachment, macular schisis, macular holes, epiretinal membranes, macular neovascularization, were excluded. Patients with systemic diseases such as hypertension, diabetes, and cardiovascular diseases that affected fundus characteristics were excluded. To avoid inter-eye effects, only the left side of the eye was selected for analysis.

In order to obtained the accuracy refraction, all patients underwent both subjective and objective refraction measurements conducted before and after cycloplegia. Spherical power ranged from −0.25 D to −13.00 D, and cylinder power ranged from −0.50 D to −6.00 D. The spherical equivalent (SE) of the refractive error was calculated as spherical error +1/2 cylindrical error. The subjects included in this study were all cases of simple myopia. Data were divided into two groups: the non-HM group (SE > −6.00 D and ≤ −0.50 D) and the HM group (SE ≤ −6.00 D) (Jong et al., 2021). Fundus images were taken using an automated image quality control system, with images centered on the macula and always included the optic disc.

2.3 Feature recognition and quantification

Object detection and semantic segmentation algorithms were used to extract and quantify retinal features, including quantitative information on tessellation, optic discs, optic cups, parapapillary atrophy, and retinal vascular structures. The annotation was performed using a semi-automatic machine-assisted workflow. Two board-certified ophthalmologists (specializing in fundus diseases with ≥5 years of experience) independently validated the annotations. The first physician performed initial corrections, and the second (with ≥8 years of experience) served as adjudicator. Given the second reviewer’s greater clinical experience, their annotations were designated as ground truth. The methodology details are provided in Document S1.

2.4 Model construction and evaluation

Data from the baseline examinations were incorporated into a ML modeling. Five ML algorithms (eXtreme Gradient Boosting (XGBoost), Categorical Boosting (CatBoost), Light Gradient Boosting Machine (LightGBM), Random Forest (RF), and Neural Network (NN)) were utilized to build the model. To mitigate overfitting, the data were randomly split into separate training and testing datasets. Model performance was assessed using 5-fold cross-validation (Document S2). Evaluation metrics included sensitivity, specificity, accuracy, precision, F1 score, and the area under the receiver operating characteristic curve (AUROC). Higher values for these metrics indicate superior model performance. The SHAP method evaluated the importance of each input variable by analyzing the average impact on model output with and without the target input variable. Model training was conducted using Python 3.6 (https://www.python.org).

2.5 Statistical analyses

Data were analyzed using SPSS (version 26.0; IBM Corp., Armonk, NY, United States). Outliers were identified and removed using Tukey’s test with a coefficient of 1.5. Categorical data were expressed as frequencies and percentages and compared using the χ2 test. The normality of continuous parameters was assessed with the Kolmogorov–Smirnov test, with P > 0.05 indicating a normal distribution. Normally distributed data were expressed as mean ± standard deviation and compared using independent samples t-tests. A two-tailed p-value of 0.05 or less was considered statistically significant. To maintain consistency, data from the same eye was selected for analyses to avoid inter-eye variability.

3 Results

3.1 Basic patient information

The study followed the steps shown in Figure 1, which included clinical data collection, image segmentation and quantification, feature engineering, model construction and evaluation, and a SHAP analysis.

Figure 1
Flowchart detailing a study on eye data processing, comprising four sections. Section 1: Data selection from 4048 patients, with exclusions leading to 2981 eyes included. Section 2: Depicts ROI segmentation and feature quantification, comparing original and segmented eye images with analysis charts. Section 3: SHAP values, displaying performance metrics like precision and sensitivity of various models. Section 4: Model construction and evaluation, showing SHAP value distributions across four iterations, included graphs and a confusion matrix.

Figure 1. Flow chart for data selection and analyses in the study.

A total of 2,981 eyes were included. Of these, 1,191 (39.9%) were classified as highly myopic, whereas 1790 (60.1%) were non-highly myopic. The following variables including age, sex, refraction, intraocular pressure, visual acuity, and 29 retinal characteristics were obtained. The baseline analysis of the dataset is shown in Table 1. The average age of 2,981 patients were 22.71 ± 5.24 years, including 1539 males (51.63%) and 1442 females (48.37%).

Table 1
www.frontiersin.org

Table 1. Baseline characteristics of the participants included in the analysis.

3.2 Algorithm evaluation

All participants were divided into five equal subsets, with one subset serving as the validation set and the remaining four subset serving as the training set for 5-fold cross-validation. Table 2 lists the performance metrics for each model. The results showed that XGBoost performed best across the five ML, with an accuracy of 0.81 (95% confidence interval [CI], 0.78–0.85) and an AUROC of 0.87 (95% CI, 0.84–0.89).

Table 2
www.frontiersin.org

Table 2. Performance of different models in the validation set.

Feature importance in the XGBoost model was assessed using feature importance bars based on average absolute SHAP values (Figure 2A). The larger absolute values represented significant influence. We summarized the top 12 features that contributed most significantly to the model predictions. Retinal tessellation density was at the top of the ranking list, followed by the mean arterial curvature, parapapillary atrophy width and area, mean vessel diameter, mean venous curvature, vertical optic disc diameters, mean venous and arterial diameters, retinal vessel density, horizontal optic disc diameters, and vessel fractal dimension. Each point corresponding to the SHAP summary plot (Figure 2B) represents the contribution of each feature of each patient to the model.

Figure 2
Bar and violin plots showcasing the impact of different eye-related features on model output. Left plot shows mean SHAP values, with

Figure 2. Feature Importance Visualization and SHAP Contributions. Left, Feature importance plot of the XGBoost model: this plot highlights the top 16 features influencing HM, identified through cross-validation. Features are ranked by their impact. The y-axis shows the average absolute SHAP values. Importance decreases from top to bottom. Right, SHAP summary plot: each row represents a feature, and each point corresponds to a sample. A wider scatter indicates a stronger effect, whereas a distribution near 0 suggests minimal influence on most samples.

The SHAP dependence plot demonstrates the effect of a single feature on the output of the XGBoost model. We visualize the SHAP values of the nine features which appeared consistently in each iteration and analyze the relationships between them (Figure 3). When the SHAP value of each feature exceeds zero, this indicates an increased risk of HM. The tessellated density value greater than 0.025, width of parapapillary atrophy greater than 400 um, parapapillary atrophy area greater than 0.6 × 106 um2 were associated with an increased risk of HM. The mean curvature of the arteries value over 0.00063, diameter of vessels over 55.2 um, curvature of the veins over 0.00128, vertical diameter of the optic disc over 1320 um, diameter of veins over 58.5 um, and diameter of artery over 47.1 um was associated with a decreased risk of HM. These findings highlight the key retinal features that significantly influenced the prediction of HM.

Figure 3
Nine scatter plots show SHAP values against various eye-related metrics. Top row: density of fundus tessellation, curvature of arteries, and width of parapapillary atrophy. Middle row: area of parapapillary atrophy, diameter of vessels, and curvature of veins. Bottom row: diameter of optic disc, diameter of veins, and diameter of artery. Each plot depicts the relationship between the SHAP values and the respective metric, illustrating patterns and trends.

Figure 3. SHAP dependence plot. The y-axis shows the feature SHAP values. The SHAP values for specific features that exceed 0 represent an increased possibility of HM. The x-axis shows the feature ranges.

4 Discussion

This study presents a novel approach of employing an explainable ML framework based on retinal imaging omics to investigate high myopia-related fundus structural parameters. Unlike previous studies that relied entirely on algorithmic results, our method improved overall understanding by capturing and interpreting the entire retinal structure. We identified and ranked 29 parameters associated with the prediction. Of these, fundus tessellation density emerged as the strongest predictor. Additionally, during validation, we found out the effect of nine consistently present features on the predicted results, including increased tessellation density, width of parapillary atrophy, area of parapillary atrophy, and reduced mean arterial curvature, vessel diameter, vein curvature, vertical optic disc diameter, vein diameter, and artery diameter. The method provides robust evidence for doctors to accurately identify and predict HM.

Based on the META-PM classification and previous research, tessellation is recognized as an early indicator of myopic maculopathy (Ohno-Matsui et al., 2021; Yii et al., 2024). Most studies have used qualitative or semi-quantitative approaches to examine the severity of myopia and the degree of tessellation (Yamashita et al., 2018; Yan et al., 2015). Some studies (Shao et al., 2021; Huang et al., 2025; Li et al., 2023) quantified tessellation density using AI image processing technology and found a correlation relationship with the SE. However, to our knowledge, this was the first quantitative ordering of the structure of the entire retina. We measured tessellation patterns by calculating the average choroidal exposure per unit area to determine the tessellation density value greater than 0.025 increased the risk of HM. This feature consistently comes out on top in five validations underscores its critical role in predicting.

Arterial tortuosity was identified as the second most influential predictive factor after tessellation density. The risk of HM was reduced when the mean arterial curvature value exceeded 0.00063. Reduced arterial curvature is associated with tissue hypoxia (Malek et al., 2014). Hypoxia-induced vascular endothelial growth factor plays a role in the curvature of retinal arteriovenous system (Malek et al., 2015). Extrinsic factors such as shear stress and mechanical pressure contribute to variations in tortuosity (Kalitzeos et al., 2013). In patients with HM, increased axial length exerts lateral mechanical forces on the eye, potentially stretching the retina and making smaller arteries more susceptible to reduced curvature (Witt et al., 2006). These findings are consistent with Lim et al. (2011) observations using computer-assisted methodologies, HM exhibit significantly lower arterial tortuosity compared to non-high myopic individuals (P < 0.001). Importantly, this study also highlights the previously overlooked role of venous tortuosity in model predictions, suggesting that clinical procedures should evaluate both arterial and venous tortuosity to better assess retinal health and disease progression. Identifying the association between retinal vascular geometry and myopic changes is a critical research priority. Although previous studies have suggested that refractive error is independent of retinal vessel caliber (Cheung et al., 2007a; Cheung et al., 2007b), we founded that diameter of vessels less than 55.2 um, veins less than 58.5 um,and artery less than 47.1 um were positively associated with an increased risk of HM. Unlike previous studies that calculated the central retinal arteriolar and venular equivalents by summing individual retinal vessel diameters (Pappelis and Jansonius, 2023), we determined the ratio of image pixel size to the actual distance using a deep learning-recognized region of interest. This method provides a more accurate representation of the retinal vascular system (Long et al., 2022). By quantifying the global retinal vascular geometry using this method, a more accurate representation of the overall vascular condition can be obtained.

In this study, we found that both the width and area of parapapillary atrophy significantly contributed to the prediction of the model. Typically, parapapillary atrophy is located horizontally on the temporal side, followed by the supratemporal and inferotemporal regions, and this distribution affects the measurement of the height of atrophy (Fang et al., 2018). In addition, the gamma region cannot be accurately identified because the border of the Brush’s membrane could not be detected by fundus photography. Notably, the vertical diameter of the optic disc emerged as a significant factor in each iteration, offering new insights for clinical application in HM.

To date, the analysis of myopic fundus characteristics has relied primarily on subjective evaluations (Hu et al., 2020). Despite the ability to manually label a few structures, the lack of objective quantification due to empirical differences and methodological limitations precludes the comprehensive assessment of complex retinal parameters (Sun J. et al., 2018; Guo et al., 2021). Thus, the relationship between these structures and the severity of myopia remains poorly understood. Artificial intelligence (AI) offers the potential to improve the identification and segmentation for these structures. Our study employed the SHAP-XGBoost to quantify the impact of retinal characteristics on HM. XGBoost, a versatile nonparametric model, provides superior accuracy over linear models. Due to its ability to handle heterogeneous and high-dimensional clinical data, and its robustness in dealing with missing values and complex variables of different categories, it is particularly suitable for the structured and multi-factor nature of ophthalmic datasets (Chen and Guestrin, 2016). This method has been applied in previous studies on long-term outcomes of refractive surgery (Kim et al., 2022) and prognosis prediction of ocular trauma (Meng et al., 2024). Although improving predictive accuracy is essential, enhancing model interpretability is equally important for clinical credibility and utility. SHAP summary and dependence plots clarify the contribution of each variable to the model, providing valuable insights for interpretability in AI-based predictions.

4.1 Study limitations

This study has several limitations. First, although the data were collected from a single center, the cross-validation method confirmed that the model demonstrated satisfactory performance. Second, the cohort comprised solely adult patients with myopia. Future studies should include analyses of individuals older than 45 years and adolescents with myopia. Finally, the main purpose of this study is to investigate the role of retinal characteristics in HM, and the SHAP analysis provided valuable insights regarding how various features influenced the prediction model, although the axial length was not fully collected before the operation.

5 Conclusion

This study effectively utilized explainable ML frameworks based on retinal imaging omics to rank the contributions of 29 retinal features. Tessellation density emerged as the most significant marker, and when its value was greater than 0.025, the risk of developing HM increased significantly. In addition, parapillary atrophy width and area were associated with increased risk of HM prediction, mean curvature of the arteries, diameter of vessels, curvature of the veins, vertical diameter of the optic disc, diameter of veins, and diameter of artery was associated with reduced risk predicted by HM. These findings can enhance clinician understanding and confidence in using ML for high myopia prediction, thus expand its potential and applicability in ophthalmology.

Data availability statement

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

Ethics statement

This study was registered with the Chinese Clinical Trial Register and approved by the Ethics Committee of Tianjin Eye Hospital. 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

HZ: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Methodology, Project administration, Writing – original draft, Writing – review and editing. JL: Writing – original draft, Writing – review and editing. SS: Methodology, Visualization, Writing – review and editing. SL: Methodology, Software, Writing – review and editing. QF: Writing – review and editing. YH: Writing – review and editing. ZD: Software, Writing – review and editing. GH: Writing – review and editing. SW: Writing – review and editing. YW: Funding acquisition, Project administration, Resources, Writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by the National Program on Key Research Project of China (2022YFC2404502), Tianjin Health Research Project (No. TJWJ2024QN076), and Tianjin Key Medical Discipline Construction (No. TJYXZDXK-3-004A-3).

Conflict of interest

Authors SS and SW were employed by HuaHui Jian AI Tech Ltd.

Authors SL and ZD were employed by EVision technology (Beijing) Co., Ltd.

The remaining 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.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2025.1609639/full#supplementary-material

References

Baird, P. N., Saw, S. M., Lanca, C., Guggenheim, J. A., Smith III, E. L., Zhou, X., et al. (2020). Myopia. Nat. Rev. Dis. Prim. 6 (1), 99. doi:10.1038/s41572-020-00231-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, T., and Guestrin, C. (2016). “Xgboost: a scalable tree-boosting system,” in Proceedings of the 22nd ACM sigkdd international conference on knowledge discovery and data mining (San Francisco, CA, USA), 785–794.

Google Scholar

Chen, N., Zhu, Z., Yang, W., and Wang, Q. (2024). Progress in clinical research and applications of retinal vessel quantification technology based on fundus imaging. Front. Bioeng. Biotechnol. 12, 1329263. doi:10.3389/fbioe.2024.1329263

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheung, N., Tong, L., Tikellis, G., Saw, S. M., Mitchell, P., Wang, J. J., et al. (2007a). Relationship of retinal vascular caliber with optic disc diameter in children. Invest Ophthalmol. Vis. Sci. 48 (11), 4945–4948. doi:10.1167/iovs.07-0472

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheung, N., Tikellis, G., Saw, S. M., Amirul Islam, F., Mitchell, P., Wang, J. J., et al. (2007b). Relationship of axial length and retinal vascular caliber in children. Am. J. Ophthalmol. 144 (5), 658–662.e1. doi:10.1016/j.ajo.2007.07.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Conti, A., Duggento, A., Indovina, I., Guerrisi, M., and Toschi, N. (2021). Radiomics in breast cancer classification and prediction. Semin. Cancer Biol. 72, 238–250. doi:10.1016/j.semcancer.2020.04.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Fang, Y., Yokoi, T., Nagaoka, N., Shinohara, K., Onishi, Y., Ishida, T., et al. (2018). Progression of myopic maculopathy during 18-Year Follow-up. Ophthalmology 125 (6), 863–877. doi:10.1016/j.ophtha.2017.12.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, X., Chen, X., Li, M., Li, S., You, R., and Wang, Y. (2021). Association between morphological characteristics of the optic disc and other anatomical features of the fundus in highly myopic eyes. Eur. J. Ophthalmol. 31 (5), 2329–2338. doi:10.1177/1120672120945901

PubMed Abstract | CrossRef Full Text | Google Scholar

Haarman, A. E. G., Enthoven, C. A., Tideman, J. W. L., Tedja, M. S., Verhoeven, V. J. M., and Klaver, C. C. W. (2020). The complications of myopia: a review and meta-analysis. Invest Ophthalmol. Vis. Sci. 61 (4), 49. doi:10.1167/iovs.61.4.49

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, G., Chen, Q., Xu, X., Lv, H., Du, Y., Wang, L., et al. (2020). Morphological characteristics of the optic nerve head and choroidal thickness in high myopia. Invest Ophthalmol. Vis. Sci. 61 (4), 46. doi:10.1167/iovs.61.4.46

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, D., Lin, X., Zhu, H., Ling, S., Dong, Z., Ke, X., et al. (2025). The associations between myopia and fundus tessellation in school children: a comparative analysis of macular and peripapillary regions using deep learning. Transl. Vis. Sci. Technol. 14 (1), 4. doi:10.1167/tvst.14.1.4

PubMed Abstract | CrossRef Full Text | Google Scholar

Jong, M., Jonas, J. B., Wolffsohn, J. S., Berntsen, D. A., Cho, P., Clarkson-Townsend, D., et al. (2021). IMI 2021 yearly digest. Invest Ophthalmol. Vis. Sci. 62 (5), 7. doi:10.1167/iovs.62.5.7

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalitzeos, A. A., Lip, G. Y., and Heitmar, R. (2013). Retinal vessel tortuosity measures and their applications. Exp. Eye Res. 106, 40–46. doi:10.1016/j.exer.2012.10.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, J., Ryu, I. H., Kim, J. K., Lee, I. S., Kim, H. K., Han, E., et al. (2022). Machine learning predicting myopic regression after corneal refractive surgery using preoperative data and fundus photography. Graefes Arch. Clin. Exp. Ophthalmol. 260 (11), 3701–3710. doi:10.1007/s00417-022-05738-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, T., Bo, W., Hu, C., Kang, H., Liu, H., Wang, K., et al. (2021). Applications of deep learning in fundus images: a review. Med. Image Anal. 69, 101971. doi:10.1016/j.media.2021.101971

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, R., Guo, X., Zhang, X., Lu, X., Wu, Q., Tian, Q., et al. (2023). Application of artificial intelligence to quantitative assessment of fundus tessellated density in young adults with different refractions. Ophthalmic Res. 66 (1), 706–716. doi:10.1159/000529639

PubMed Abstract | CrossRef Full Text | Google Scholar

Lim, L. S., Cheung, C. Y., Lin, X., Mitchell, P., Wong, T. Y., and Mei-Saw, S. (2011). Influence of refractive error and axial length on retinal vessel geometric characteristics. Invest Ophthalmol. Vis. Sci. 52 (2), 669–678. doi:10.1167/iovs.10-6184

PubMed Abstract | CrossRef Full Text | Google Scholar

Long, T., Xu, Y., Zou, H., Lu, L., Yuan, T., Dong, Z., et al. (2022). A generic pixel pitch calibration method for fundus camera via automated ROI extraction. Sensors (Basel) 22 (21), 8565. doi:10.3390/s22218565

PubMed Abstract | CrossRef Full Text | Google Scholar

Malek, J., Azar, A. T., and Tourki, R. (2014). Impact of retinal vascular tortuosity on retinal circulation. Neural. comput. Appl. 26 (1), 25–40.

CrossRef Full Text | Google Scholar

Malek, J., Azar, A. T., and Tourki, R. (2015). Impact of retinal vascular tortuosity on retinal circulation. Neural Comput. Applic 26, 25–40. doi:10.1007/s00521-014-1657-2

CrossRef Full Text | Google Scholar

Meng, X., Wang, Q., Chen, S., Zhang, S., Yu, J., Li, H., et al. (2024). An interpretable model predicts visual outcomes of no light perception eyes after open globe injury. Br. J. Ophthalmol. 108 (2), 285–293. doi:10.1136/bjo-2022-322753

PubMed Abstract | CrossRef Full Text | Google Scholar

Modjtahedi, B. S., Abbott, R. L., Fong, D. S., Lum, F., and Tan, D.Task Force on Myopia (2021). Reducing the global burden of myopia by delaying the onset of myopia and reducing myopic progression in children: the academy's task force on myopia. Ophthalmology 128 (6), 816–826. doi:10.1016/j.ophtha.2020.10.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Ohno-Matsui, K., Lai, T. Y., Lai, C. C., and Cheung, C. M. (2016). Updates of pathologic myopia. Prog. Retin Eye Res. 52, 156–187. doi:10.1016/j.preteyeres.2015.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Ohno-Matsui, K., Wu, P. C., Yamashiro, K., Vutipongsatorn, K., Fang, Y., Cheung, C. M. G., et al. (2021). IMI pathologic myopia. Invest Ophthalmol. Vis. Sci. 62 (5), 5. doi:10.1167/iovs.62.5.5

PubMed Abstract | CrossRef Full Text | Google Scholar

Pappelis, K., and Jansonius, N. M. (2023). Retinal vessel caliber measurement bias in fundus images in the presence of the central light reflex. Transl. Vis. Sci. Technol. 12 (7), 16. doi:10.1167/tvst.12.7.16

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruiz-Medrano, J., Montero, J. A., Flores-Moreno, I., Arias, L., García-Layana, A., and Ruiz-Moreno, J. M. (2019). Myopic maculopathy: current status and proposal for a new classification and grading system (ATN). Prog. Retin Eye Res. 69, 80–115. doi:10.1016/j.preteyeres.2018.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Salih, A., Boscolo Galazzo, I., Gkontra, P., Lee, A. M., Lekadir, K., Raisi-Estabragh, Z., et al. (2023). Explainable artificial intelligence and cardiac imaging: toward more interpretable models. Circ. Cardiovasc Imaging 16 (4), e014519. doi:10.1161/circimaging.122.014519

PubMed Abstract | CrossRef Full Text | Google Scholar

Shao, L., Zhang, Q. L., Long, T. F., Dong, L., Zhang, C., Da Zhou, W., et al. (2021). Quantitative assessment of fundus tessellated density and associated factors in fundus images using artificial intelligence. Transl. Vis. Sci. Technol. 10 (9), 23. doi:10.1167/tvst.10.9.23

PubMed Abstract | CrossRef Full Text | Google Scholar

Shinojima, A., Negishi, K., Tsubota, K., and Kurihara, T. (2022). Multiple factors causing myopia and the possible treatments: a mini review. Front. Public Health 10, 897600. doi:10.3389/fpubh.2022.897600

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, R., Limkin, E. J., Vakalopoulou, M., Dercle, L., Champiat, S., Han, S. R., et al. (2018). A radiomics approach to assess tumour-infiltrating CD8 cells and response to anti-PD-1 or anti-PD-L1 immunotherapy: an imaging biomarker, retrospective multicohort study. Lancet Oncol. 19 (9), 1180–1191. doi:10.1016/s1470-2045(18)30413-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, J., Wang, J., You, R., and Wang, Y. (2018). Is the retinal vasculature related to β-Peripapillary atrophy in nonpathological high myopia? An optical coherence tomography angiography study in Chinese adults. J. Ophthalmol. 2018, 1–8. doi:10.1155/2018/7895238

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Q., Chen, Y., Liang, C., Zhao, Y., Lv, X., Zou, Y., et al. (2021). Biologic pathways underlying prognostic radiomics phenotypes from paired MRI and RNA sequencing in glioblastoma. Radiology 301 (3), 654–663. doi:10.1148/radiol.2021203281

PubMed Abstract | CrossRef Full Text | Google Scholar

Verhoeven, V. J., Wong, K. T., Buitendijk, G. H., Hofman, A., Vingerling, J. R., and Klaver, C. C. (2015). Visual consequences of refractive errors in the general population. Ophthalmology 122 (1), 101–109. doi:10.1016/j.ophtha.2014.07.030

PubMed Abstract | CrossRef Full Text | Google Scholar

von Elm, E., Altman, D. G., Egger, M., Pocock, S. J., Gøtzsche, P. C., and Vandenbroucke, J. P. (2007). The strengthening the reporting of observational studies in epidemiology (STROBE) statement: guidelines for reporting observational studies. PLoS Med. 4 (10), e296. doi:10.1371/journal.pmed.0040296

PubMed Abstract | CrossRef Full Text | Google Scholar

Witt, N., Wong, T. Y., Hughes, A. D., Chaturvedi, N., Klein, B. E., Evans, R., et al. (2006). Abnormalities of retinal microvascular structure and risk of mortality from ischemic heart disease and stroke. Hypertension 47 (5), 975–981. doi:10.1161/01.hyp.0000216717.72048.6c

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, L., Zhuang, Y., Zhang, G., Ma, Y., Yuan, J., Tu, C., et al. (2021). Design, methodology, and baseline of whole city-million scale children and adolescents myopia survey (CAMS) in wenzhou, China. Eye Vis. (Lond). 8 (1), 31. doi:10.1186/s40662-021-00255-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamashita, T., Iwase, A., Kii, Y., Sakai, H., Terasaki, H., Sakamoto, T., et al. (2018). Location of ocular tessellations in Japanese: population-based kumejima study. Invest Ophthalmol. Vis. Sci. 59 (12), 4963–4967. doi:10.1167/iovs.18-25007

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, Y. N., Wang, Y. X., Xu, L., Xu, J., Wei, W. B., and Jonas, J. B. (2015). Fundus tessellation: prevalence and associated factors. Ophthalmology 122 (9), 1873–1880. doi:10.1016/j.ophtha.2015.05.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Yii, F., Nguyen, L., Strang, N., Bernabeu, M. O., Tatham, A. J., MacGillivray, T., et al. (2024). Factors associated with pathologic myopia onset and progression: a systematic review and meta-analysis. Ophthalmic Physiol. Opt. 44 (5), 963–976. doi:10.1111/opo.13312

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, Y.D., He, Z., Ouyang, J., Tan, Y., Chen, Y., Gu, Y., et al. (2021). Magnetic resonance imaging radiomics predicts preoperative axillary lymph node metastasis to support surgical decisions and is associated with tumor microenvironment in invasive breast cancer: a machine learning, multicenter study. Ebiomedicine 69, 103460. doi:10.1016/j.ebiom.2021.103460

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Zhang, H., Jiang, M., Li, J., Li, J., Zhou, H., et al. (2025). Radiomics in ophthalmology: a systematic review. Eur. Radiol. 35 (1), 542–557. doi:10.1007/s00330-024-10911-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: high myopia, retinal imaging omics, machine learning, deep learning, shapley additive exPlanation

Citation: Zou H, Liu J, Shi S, Ling S, Fan Q, Huo Y, Dong Z, Han G, Wang S and Wang Y (2025) Retinal features as predictive indicators for high myopia: insights from explainable multi-machine learning models. Front. Bioeng. Biotechnol. 13:1609639. doi: 10.3389/fbioe.2025.1609639

Received: 22 April 2025; Accepted: 29 September 2025;
Published: 13 October 2025.

Edited by:

Weihua Yang, Southern Medical University, China

Reviewed by:

Tae Keun Yoo, Hangil Eye Hospital, Republic of Korea
Peng Gao, Tongji University, China

Copyright © 2025 Zou, Liu, Shi, Ling, Fan, Huo, Dong, Han, Wang and Wang. 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: Yan Wang, d2FuZ3lhbjcxNDNAdmlwLnNpbmEuY29t

These authors have contributed equally to this work

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.