Latest developments of generative artificial intelligence and applications in ophthalmology.

The emergence of generative artificial intelligence (AI) has revolutionized various fields. In ophthalmology, generative AI has the potential to enhance efficiency, accuracy, personalization and innovation in clinical practice and medical research, through processing data, streamlining medical documentation, facilitating patient-doctor communication, aiding in clinical decision-making, and simulating clinical trials. This review focuses on the development and integration of generative AI models into clinical workflows and scientific research of ophthalmology. It outlines the need for development of a standard framework for comprehensive assessments, robust evidence, and exploration of the potential of multimodal capabilities and intelligent agents. Additionally, the review addresses the risks in AI model development and application in clinical service and research of ophthalmology, including data privacy, data bias, adaptation friction, over interdependence, and job replacement, based on which we summarized a risk management framework to mitigate these concerns. This review highlights the transformative potential of generative AI in enhancing patient care, improving operational efficiency in the clinical service and research in ophthalmology. It also advocates for a balanced approach to its adoption.

Retinal Microvasculature Causally Affects the Brain Cortical Structure: A Mendelian Randomization Study.

Purpose: To reveal the causality between retinal vascular density (VD), fractal dimension (FD), and brain cortex structure using Mendelian randomization (MR). Design: Cross-sectional study. Participants: Genome-wide association studies of VD and FD involving 54 813 participants from the United Kingdom Biobank were used. The brain cortical features, including the cortical thickness (TH) and surface area (SA), were extracted from 51 665 patients across 60 cohorts. Surface area and TH were measured globally and in 34 functional regions using magnetic resonance imaging. Methods: Bidirectional univariable MR (UVMR) was used to detect the causality between FD, VD, and brain cortex structure. Multivariable MR (MVMR) was used to adjust for confounding factors, including body mass index and blood pressure. Main Outcome Measures: The global and regional measurements of brain cortical SA and TH. Results: At the global level, higher VD is related to decreased TH (ß = -0.0140 mm, 95% confidence interval: -0.0269 mm to -0.0011 mm, P = 0.0339). At the functional level, retinal FD is related to the TH of banks of the superior temporal sulcus and transverse temporal region without global weighted, as well as the SA of the posterior cingulate after adjustment. Vascular density is correlated with the SA of subregions of the frontal lobe and temporal lobe, in addition to the TH of the inferior temporal, entorhinal, and pars opercularis regions in both UVMR and MVMR. Bidirectional MR studies showed a causation between the SA of the parahippocampal and cauda middle frontal gyrus and retinal VD. No pleiotropy was detected. Conclusions: Fractal dimension and VD causally influence the cortical structure and vice versa, indicating that the retinal microvasculature may serve as a biomarker for cortex structural changes. Our study provides insights into utilizing noninvasive fundus images to predict cortical structural deteriorations and neuropsychiatric disorders. Financial Disclosures: The author(s) have no proprietary or commercial interest in any materials discussed in this article.

Development and Evaluation of a Retrieval-Augmented Large Language Model Framework for Ophthalmology.

Importance: Although augmenting large language models (LLMs) with knowledge bases may improve medical domain-specific performance, practical methods are needed for local implementation of LLMs that address privacy concerns and enhance accessibility for health care professionals. Objective: To develop an accurate, cost-effective local implementation of an LLM to mitigate privacy concerns and support their practical deployment in health care settings. Design, Setting, and Participants: ChatZOC (Sun Yat-Sen University Zhongshan Ophthalmology Center), a retrieval-augmented LLM framework, was developed by enhancing a baseline LLM with a comprehensive ophthalmic dataset and evaluation framework (CODE), which includes over 30 000 pieces of ophthalmic knowledge. This LLM was benchmarked against 10 representative LLMs, including GPT-4 and GPT-3.5 Turbo (OpenAI), across 300 clinical questions in ophthalmology. The evaluation, involving a panel of medical experts and biomedical researchers, focused on accuracy, utility, and safety. A double-masked approach was used to try to minimize bias assessment across all models. The study used a comprehensive knowledge base derived from ophthalmic clinical practice, without directly involving clinical patients. Exposures: LLM response to clinical questions. Main Outcomes and Measures: Accuracy, utility, and safety of LLMs in responding to clinical questions. Results: The baseline model achieved a human ranking score of 0.48. The retrieval-augmented LLM had a score of 0.60, a difference of 0.12 (95% CI, 0.02-0.22; P = .02) from baseline and not different from GPT-4 with a score of 0.61 (difference = 0.01; 95% CI, -0.11 to 0.13; P = .89). For scientific consensus, the retrieval-augmented LLM was 84.0% compared with the baseline model of 46.5% (difference = 37.5%; 95% CI, 29.0%-46.0%; P < .001) and not different from GPT-4 with a value of 79.2% (difference = 4.8%; 95% CI, -0.3% to 10.0%; P = .06). Conclusions and Relevance: Results of this quality improvement study suggest that the integration of high-quality knowledge bases improved the LLM's performance in medical domains. This study highlights the transformative potential of augmented LLMs in clinical practice by providing reliable, safe, and practical clinical information. Further research is needed to explore the broader application of such frameworks in the real world.

Shaping school for childhood myopia: the association between floor area ratio of school environment and myopia in China.

AIM: To investigate the association of floor area ratio (FAR), an indicator of built environments, and myopia onset. METHODS: This prospective cohort study recruited 136 753 children aged 6-10 years from 108 schools in Shenzhen, China at baseline (2016-2017). Refractive power was measured with non-cycloplegic autorefraction over a 2-year follow-up period. FAR was objectively evaluated using geographical information system technology. Mixed-effects logistic regression models were constructed to examine the association of FAR with a 2-year cumulative incidence of myopia among individuals without baseline myopia; multiple linear regression model, with a 2-year cumulative incidence rate of myopia at each school. RESULTS: Of 101 624 non-myopic children (56.3% boys; mean (SE) age, 7.657±1.182 years) included in the study, 26 391 (26.0%) of them developed myopia after 2 years. In the individual-level analysis adjusting for demographic, socioeconomic and greenness factors, an IQR in FAR was associated with a decreased risk of 2-year myopia incidence (OR 0.898, 95% CI 0.866 to 0.932, p<0.001). Similar findings were observed in the analysis additionally adjusted for genetic and behavioural factors (OR 0.821, 95% CI 0.766 to 0.880, p<0.001). In the school-level, an IQR increase in FAR was found to be associated with a 2.0% reduction in the 2-year incidence rate of myopia (95% CI 1.3% to 2.6%, p<0.001). CONCLUSIONS: Exposure to higher FAR was associated with a decreased myopia incidence, providing insights into myopia prevention through school built environments in China.

Digital ray: enhancing cataractous fundus images using style transfer generative adversarial networks to improve retinopathy detection.

BACKGROUND/AIMS: The aim of this study was to develop and evaluate digital ray, based on preoperative and postoperative image pairs using style transfer generative adversarial networks (GANs), to enhance cataractous fundus images for improved retinopathy detection. METHODS: For eligible cataract patients, preoperative and postoperative colour fundus photographs (CFP) and ultra-wide field (UWF) images were captured. Then, both the original CycleGAN and a modified CycleGAN (C2ycleGAN) framework were adopted for image generation and quantitatively compared using Frechet Inception Distance (FID) and Kernel Inception Distance (KID). Additionally, CFP and UWF images from another cataract cohort were used to test model performances. Different panels of ophthalmologists evaluated the quality, authenticity and diagnostic efficacy of the generated images. RESULTS: A total of 959 CFP and 1009 UWF image pairs were included in model development. FID and KID indicated that images generated by C2ycleGAN presented significantly improved quality. Based on ophthalmologists' average ratings, the percentages of inadequate-quality images decreased from 32% to 18.8% for CFP, and from 18.7% to 14.7% for UWF. Only 24.8% and 13.8% of generated CFP and UWF images could be recognised as synthetic. The accuracy of retinopathy detection significantly increased from 78% to 91% for CFP and from 91% to 93% for UWF. For retinopathy subtype diagnosis, the accuracies also increased from 87%-94% to 91%-100% for CFP and from 87%-95% to 93%-97% for UWF. CONCLUSION: Digital ray could generate realistic postoperative CFP and UWF images with enhanced quality and accuracy for overall detection and subtype diagnosis of retinopathies, especially for CFP.\ TRIAL REGISTRATION NUMBER: This study was registered with ClinicalTrials.gov (NCT05491798).

Protocol for performing deep learning-based fundus fluorescein angiography image analysis with classification and segmentation tasks.

Fundus fluorescein angiography (FFA) examinations are widely used in the evaluation of fundus disease conditions to facilitate further treatment suggestions. Here, we present a protocol for performing deep learning-based FFA image analytics with classification and segmentation tasks. We describe steps for data preparation, model implementation, statistical analysis, and heatmap visualization. The protocol is applicable in Python using customized data and can achieve the whole process from diagnosis to treatment suggestion of ischemic retinal diseases. For complete details on the use and execution of this protocol, please refer to Zhao et al.1.

Classification of congenital cataracts based on multidimensional phenotypes and its association with visual outcomes.

AIM: To establish a classification for congenital cataracts that can facilitate individualized treatment and help identify individuals with a high likelihood of different visual outcomes. METHODS: Consecutive patients diagnosed with congenital cataracts and undergoing surgery between January 2005 and November 2021 were recruited. Data on visual outcomes and the phenotypic characteristics of ocular biometry and the anterior and posterior segments were extracted from the patients' medical records. A hierarchical cluster analysis was performed. The main outcome measure was the identification of distinct clusters of eyes with congenital cataracts. RESULTS: A total of 164 children (299 eyes) were divided into two clusters based on their ocular features. Cluster 1 (96 eyes) had a shorter axial length (mean±SD, 19.44±1.68 mm), a low prevalence of macular abnormalities (1.04%), and no retinal abnormalities or posterior cataracts. Cluster 2 (203 eyes) had a greater axial length (mean±SD, 20.42±2.10 mm) and a higher prevalence of macular abnormalities (8.37%), retinal abnormalities (98.52%), and posterior cataracts (4.93%). Compared with the eyes in Cluster 2 (57.14%), those in Cluster 1 (71.88%) had a 2.2 times higher chance of good best-corrected visual acuity [<0.7 logMAR; OR (95%CI), 2.20 (1.25-3.81); P=0.006]. CONCLUSION: This retrospective study categorizes congenital cataracts into two distinct clusters, each associated with a different likelihood of visual outcomes. This innovative classification may enable the personalization and prioritization of early interventions for patients who may gain the greatest benefit, thereby making strides toward precision medicine in the field of congenital cataracts.

Association Between Preoperative Ocular Parameters and Myopic Shift in Children Undergoing Primary Intraocular Lens Implantation.

Purpose: To evaluate the association between preoperative ocular parameters and myopic shift following primary intraocular lens (IOL) implantation in pediatric cataracts. Methods: Eyes from pediatric patients undergoing bilateral cataract surgery with primary IOL implantation were included. Eyes were grouped by age at surgery and subdivided into three axial length (AL) subgroups and three keratometry subgroups. Mixed-effects linear regression was utilized to assess the trend in myopic shift among subgroups. Multivariable analysis was performed to determine factors associated with myopic shift. Results: A total of 222 eyes were included. The median age at surgery was 4.36 years (interquartile range [IQR], 3.16-6.00 years) and the median follow-up was 4.18 years (IQR, 3.48-4.64 years). As preoperative AL increased, a decreased trend was observed in myopic shift and rate of myopic shift (P = 0.008 and P = 0.003, respectively, in the 4 to <6 years old group; P = 0.002 and P < 0.001, respectively, in the ≥6 years old group). Greater myopic shift and rate of myopic shift were associated with younger age at surgery (P = 0.008 and P = 0.008, respectively). Both myopic shift and rate of myopic shift were negatively associated with AL. Conclusions: Age at surgery and preoperative AL were associated with myopic shift in pediatric cataracts following primary IOL implantation. Adjusting the target refraction based on preoperative AL could potentially improve patients' long-term refractive outcome. Translational Relevance: This study may help to guide the selection of postoperative target refraction according to age at surgery and preoperative ocular parameters for pediatric cataracts.

Cost-effectiveness and cost-utility of a digital technology-driven hierarchical healthcare screening pattern in China.

Utilization of digital technologies for cataract screening in primary care is a potential solution for addressing the dilemma between the growing aging population and unequally distributed resources. Here, we propose a digital technology-driven hierarchical screening (DH screening) pattern implemented in China to promote the equity and accessibility of healthcare. It consists of home-based mobile artificial intelligence (AI) screening, community-based AI diagnosis, and referral to hospitals. We utilize decision-analytic Markov models to evaluate the cost-effectiveness and cost-utility of different cataract screening strategies (no screening, telescreening, AI screening and DH screening). A simulated cohort of 100,000 individuals from age 50 is built through a total of 30 1-year Markov cycles. The primary outcomes are incremental cost-effectiveness ratio and incremental cost-utility ratio. The results show that DH screening dominates no screening, telescreening and AI screening in urban and rural China. Annual DH screening emerges as the most economically effective strategy with 341 (338 to 344) and 1326 (1312 to 1340) years of blindness avoided compared with telescreening, and 37 (35 to 39) and 140 (131 to 148) years compared with AI screening in urban and rural settings, respectively. The findings remain robust across all sensitivity analyses conducted. Here, we report that DH screening is cost-effective in urban and rural China, and the annual screening proves to be the most cost-effective option, providing an economic rationale for policymakers promoting public eye health in low- and middle-income countries.

Metabolomic profiling of the aqueous humor in patients with pediatric cataract.

Pediatric cataract, including congenital and developmental cataract, is a kind of pediatric vision-threatening disease with extensive phenotypic heterogeneity and multiple mechanisms. We aimed to investigate the metabolite profile of aqueous humor (AH) in patients with pediatric cataracts, and identify underlying mutual correlations between differential metabolites. Metabolomic profiles of AH were analyzed and compared between pediatric cataract patients (n = 33) and age-related cataract patients without metabolic diseases (n = 29), using global untargeted metabolomics with ultra-high-performance liquid chromatography tandem mass spectrometry. Principal component analysis, partial least squares discriminant analysis and heat map were applied. Enriched pathway analysis was conducted using Kyoto Encyclopedia of Genes and Genomes. Receiver-operating characteristic (ROC) analyses were employed to select potential biomarkers. A total of 318 metabolites were identified, of which 54 differential metabolites (25 upregulated and 29 downregulated) were detected in pediatric cataract group compared with controls (variable importance of projection >1.0, fold change ≥1.5 or ≤ 0.667 and P < 0.05). A significant accumulation of N-Acetyl-Dl-glutamic acid was observed in pediatric cataract group. The differential metabolites were mainly enriched in histidine metabolism (increased L-Histidine and decreased 1-Methylhistamine) and the tryptophan metabolism (increased N-Formylkynurenine and L-Kynurenine). 5-Aminosalicylic acid showed strong positive mutual inter-correlation with L-Tyrosinemethylester and N,N-Diethylethanolamine, both of which were down-regulated in pediatric cataract group. The ROC analysis implied 11 metabolites served as potential biomarkers for pediatric cataract patients (all area under the ROC curve ≥0.900). These results illustrated novel potential metabolites and metabolic pathways in pediatric cataract, which provides new insights into the pathophysiology of pediatric cataract.

Combination of serum markers with optical coherence tomography angiography for evaluating neuromyelitis optica spectrum disorders and multiple sclerosis.

BACKGROUND: Neuromyelitis optica spectrum disorder (NMOSD) and multiple sclerosis (MS), autoimmune inflammatory diseases of the central nervous system, affect the optic nerve and brain. A lumbar puncture to obtain biomarkers is highly invasive. Serum biomarkers and optical coherence tomography angiography (OCTA) are more accessible and less expensive than magnetic resonance imaging and provide reliable, reproducible measures of neuroaxonal damage. This study investigated the association between serum neurofilament light chain (sNfL), serum glial fibrillary acidic protein (sGFAP), and OCTA metrics. Serum sNfL and sGFAP levels, OCTA values, and clinical characteristics were compared among 91 patients with NMOSD, 81 patients with MS, and 34 healthy controls (HCs) at baseline and 1-year follow-up. RESULTS: sNfL and sGFAP levels were higher while the sGFAP/sNfL quotients were significantly lower in NMOSD and MS patients than those in HCs. At baseline, the average thicknesses of the peripapillary retinal nerve fibre layer (pRNFL) and macular ganglion cell-inner plexiform layer (mGC-IPL) were significantly smaller in NMOSD and MS patients than those in HCs (pRNFL: MS 92.0 [80.2; 101] µm, NMOSD 80.0 [59.0; 95.8] µm, vs HC 99.0 [92.0; 104] µm, p < 0.001; mGC-IPL: MS 74.5 [64.2; 81.0] µm, NMOSD 68.0 [56.0; 81.0] µm, vs HC 83.5 [78.0; 88.0] µm, p < 0.001). The vessel density (VD) and perfusion density (PD) were increased in MS patients without optic neuritis compared to HCs (VD: MS 16.7 [15.6; 17.9] HC 15.3 [13.4; 16.9], p = 0.008; PD: MS 0.41 [0.38; 0.43], HC 0.37 [0.32; 0.41], p = 0.017). In NMOSD patients without optic neuritis, sNfL was significantly associated with PD at baseline (r = 0.329, q = 0.041). The baseline and follow-up values of the sNfL level and average pRNFL and mGC-IPL thicknesses in MS patients showed significant differences. NMOSD patients showed significant differences between baseline and follow-up sNfL and sGFAP levels but not OCTA metrics. CONCLUSION: Changes in retinal microvasculature might occur earlier than those in retinal structure and may therefore serve as a promising diagnostic marker for early NMOSD. The combination of serum markers and OCTA metrics could be used to evaluate and differentiate between MS and NMOSD.

CRB1-associated retinal degeneration is dependent on bacterial translocation from the gut.

The Crumbs homolog 1 (CRB1) gene is associated with retinal degeneration, most commonly Leber congenital amaurosis (LCA) and retinitis pigmentosa (RP). Here, we demonstrate that murine retinas bearing the Rd8 mutation of Crb1 are characterized by the presence of intralesional bacteria. While normal CRB1 expression was enriched in the apical junctional complexes of retinal pigment epithelium and colonic enterocytes, Crb1 mutations dampened its expression at both sites. Consequent impairment of the outer blood retinal barrier and colonic intestinal epithelial barrier in Rd8 mice led to the translocation of intestinal bacteria from the lower gastrointestinal (GI) tract to the retina, resulting in secondary retinal degeneration. Either the depletion of bacteria systemically or the reintroduction of normal Crb1 expression colonically rescued Rd8-mutation-associated retinal degeneration without reversing the retinal barrier breach. Our data elucidate the pathogenesis of Crb1-mutation-associated retinal degenerations and suggest that antimicrobial agents have the potential to treat this devastating blinding disease.

Early Diagnosis of Syndromic Congenital Cataracts in a Large Cohort of Congenital Cataracts.

PURPOSE: To explore the factors related to the diagnosis yield of syndromic congenital cataracts and describe the phenotype-genotype correlation in congenital cataract patients. DESIGN: Prospective cohort study. METHODS: Setting: the participants from underwent clinical examinations between 2021 and 2022. Facial and anterior eye segment photographs, pre- and postoperative ocular parameters, and medical and family histories were recorded. Bioinformatics analysis was performed using whole-exome sequencing data. Statistical and correlation analyses were performed using the basic characteristics, deep phenotype, and genotype data. PARTICIPANTS: 115 patients with unrelated congenital cataract. INTERVENTIONS: performing clinical examinations, whole-exome sequencing, and bioinformatics analysis for all participants. MAIN OUTCOMES AND MEASURES: factors related to the genetic diagnosis yield of syndromic congenital cataracts. RESULTS: Bilaterally asymmetrical cataracts were identified to be associated with syndromic congenital cataracts. The overall genetic diagnostic yield in the cohort was 72.2%. In total, 34.8% of the probands were early diagnosed with various syndromes with the help of genetic information. A phenotype-genotype correlation was detected for some genes and deep phenotypes. CONCLUSIONS: We highlight the importance of screening syndromic diseases in the patients with asymmetrical congenital cataracts. Application of whole-exome sequencing helps provide early diagnosis and treatment for the patients with syndromic congenital cataracts. This study also achieved a high genetic diagnostic yield, expanded the genotypic spectrum, and found phenotype-genotype correlations. A comprehensive analysis of cataract symmetricity, family history, and deep phenotypes makes the genotype prediction of some congenital cataract patients possible.

Predicting the risk of glaucoma-related adverse events following secondary intraocular lens implantation in paediatric eyes: a 3-year study.

AIMS: To establish and evaluate predictive models for glaucoma-related adverse events (GRAEs) following secondary intraocular lens (IOL) implantation in paediatric eyes. METHODS: 205 children (356 aphakic eyes) receiving secondary IOL implantation at Zhongshan Ophthalmic Center with a 3-year follow-up were enrolled. Cox proportional hazard model was used to identify predictors of GRAEs and developed nomograms. Model performance was evaluated with time-dependent receiver operating characteristic (ROC) curves, decision curve analysis, Kaplan-Meier curves and validated internally through C-statistics and calibration plot of the bootstrap samples. RESULTS: Older age at secondary IOL implantation (HR=1.5, 95% CI: 1.03 to 2.19), transient intraocular hypertension (HR=9.06, 95% CI: 2.97 to 27.67) and ciliary sulcus implantation (HR=14.55, 95% CI: 2.11 to 100.57) were identified as risk factors for GRAEs (all p<0.05). Two nomograms were established. At postoperatively 1, 2 and 3 years, model 1 achieved area under the ROC curves (AUCs) of 0.747 (95% CI: 0.776 to 0.935), 0.765 (95% CI: 0.804 to 0.936) and 0.748 (95% CI: 0.736 to 0.918), and the AUCs of model 2 were 0.881 (95% CI: 0.836 to 0.926), 0.895 (95% CI: 0.852 to 0.938) and 0.848 (95% CI: 0.752 to 0.945). Both models demonstrated fine clinical net benefit and performance in the interval validation. The Kaplan-Meier curves showing two distinct risk groups were well discriminated and robust in both models. An online risk calculator was constructed. CONCLUSION: Two nomograms could sensitively and accurately identify children at high risk of GRAEs after secondary IOL implantation to help early identification and timely intervention.

Green Space Morphology and School Myopia in China.

Importance: China has experienced both rapid urbanization and major increases in myopia prevalence. Previous studies suggest that green space exposure reduces the risk of myopia, but the association between myopia risk and specific geometry and distribution characteristics of green space has yet to be explored. These must be understood to craft effective interventions to reduce myopia. Objective: To evaluate the associations between myopia and specific green space morphology using novel quantitative data from high-resolution satellite imaging. Design, Setting, and Participants: This prospective cohort study included students grades 1 to 4 (aged 6 to 9 years) in Shenzhen, China. Baseline data were collected in 2016-2017, and students were followed up in 2018-2019. Data were analyzed from September 2020 to January 2022. Exposures: Eight landscape metrics were calculated using land cover data from high-resolution Gaofen-2 satellite images to measure area, aggregation, and shape of green space. Main Outcome and Measures: The 2-year cumulative change in myopia prevalence at each school and incidence of myopia at the student level after 2 years were calculated as main outcomes. The associations between landscape metrics and school myopia were assessed, controlling for geographical, demographic, and socioeconomic factors. Principal component analyses were performed to further assess the joint effect of landscape metrics at the school and individual level. Results: A total of 138 735 students were assessed at baseline. Higher proportion, aggregation, and better connectivity of green space were correlated with slower increases in myopia prevalence. In the principal component regression, a 1-unit increase in the myopia-related green space morphology index (the first principal component) was negatively associated with a 1.7% (95% CI, -2.7 to -0.6) decrease in myopia prevalence change at the school level (P = .002). At the individual level, a 1-unit increase in myopia-related green space morphology index was associated with a 9.8% (95% CI, 4.1 to 15.1) reduction in the risk of incident myopia (P < .001), and the association remained after further adjustment for outdoor time, screen time, reading time, and parental myopia (adjusted odds ratio, 0.88; 95% CI, 0.80 to 0.97; P = .009). Conclusions and Relevance: Structure of green space was associated with a decreased relative risk of myopia, which may provide guidance for construction and renovation of schools. Since risk estimates only indicate correlations rather than causation, further interventional studies are needed to assess the effect on school myopia of urban planning and environmental designs, especially size and aggregation metrics of green space, on school myopia.

Urban greenspace and visual acuity in schoolchildren: a large prospective cohort study in China.

BACKGROUND: Greenspace is known to have a positive impact on human health and well-being, but its potential effects on visual acuity have not been extensively studied. OBJECTIVES: Our aim was to examine the relationship between long-term greenspace exposure and visual acuity in children, while also exploring the potential mechanisms in this association. METHODS: We conducted this prospective cohort study based on the Children's growth environment, lifestyle, physical, and mental health development project (COHERENCE), which screened 286,801 schoolchildren in Guangzhou, China, starting in the 2016/17 academic year and followed them up for three academic years (2017/18-2019/20). Visual acuity was measured using a standardized logarithmic chart, and visual impairment was defined as visual acuity worse than 0.0 logarithm of the minimum angle of resolution (LogMAR) units in the better eye. We used the Normalized Difference Vegetation Index (NDVI), the Soil-Adjusted Vegetation Index (SAVI), and the Enhanced Vegetation Index (EVI) to assess the greenspace surrounding child's geocoded home and school at each visit. RESULTS: Our analysis indicated that higher greenspace exposure was associated with greater visual acuity z-score at baseline and with slower decline in visual acuity z-score during the 3-year follow-up. An interquartile range increase in home-school-based NDVI 300m was associated with a 7% decrease [hazard ratios (HRs): 0.93, 95% confidence interval (CI): 0.92, 0.94] in the risk of visual impairment. We also found that air pollution, physical activity, outdoor time, and recreational screen time partially mediated the greenspace-visual acuity association. CONCLUSION: Our findings suggest that increasing greenspace exposure could benefit children's visual acuity development and reduce the risk of visual impairment by reducing air pollution and recreational screen time while increasing physical activity and outdoor time. All results could have potential policy implications, given the individual and societal burdens associated with visual impairment.