RNFLT2Vec: Artifact-corrected representation learning for retinal nerve fiber layer thickness maps

Optical coherence tomography imaging provides a crucial clinical measurement for diagnosing and monitoring glaucoma through the two-dimensional retinal nerve fiber layer (RNFL) thickness (RNFLT) map. Researchers have been increasingly using neural models to extract meaningful features from the RNFLT map, aiming to identify biomarkers for glaucoma and its progression. However, accurately representing the RNFLT map features relevant to glaucoma is challenging due to significant variations in retinal anatomy among individuals, which confound the pathological thinning of the RNFL. Moreover, the presence of artifacts in the RNFLT map, caused by segmentation errors in the context of degraded image quality and defective imaging procedures, further complicates the task. In this paper, we propose a general framework called RNFLT2Vec for unsupervised learning of vectorized feature representations from RNFLT maps. Our method includes an artifact correction component that learns to rectify RNFLT values at artifact locations, producing a representation reflecting the RNFLT map without artifacts. Additionally, we incorporate two regularization techniques to encourage discriminative representation learning. Firstly, we introduce a contrastive learning-based regularization to capture the similarities and dissimilarities between RNFLT maps. Secondly, we employ a consistency learning-based regularization to align pairwise distances of RNFLT maps with their corresponding thickness distributions. Through extensive experiments on a large-scale real-world dataset, we demonstrate the superiority of RNFLT2Vec in three different clinical tasks: RNFLT pattern discovery, glaucoma detection, and visual field prediction. Our results validate the effectiveness of our framework and its potential to contribute to a better understanding and diagnosis of glaucoma.

Glucose tolerance and insulin resistance/sensitivity associate with retinal layer characteristics: the LIFE-Adult-Study

AIMS/HYPOTHESIS

As the prevalence of insulin resistance and glucose intolerance is increasing throughout the world, diabetes-induced eye diseases are a global health burden. We aim to identify distinct optical bands which are closely related to insulin and glucose metabolism, using non-invasive, high-resolution spectral domain optical coherence tomography (SD-OCT) in a large, population-based dataset.

METHODS

The LIFE-Adult-Study randomly selected 10,000 participants from the population registry of Leipzig, Germany. Cross-sectional, standardised phenotyping included the assessment of various metabolic risk markers and ocular imaging, such as SD-OCT-derived thicknesses of ten optical bands of the retina. Global and Early Treatment Diabetic Retinopathy Study (ETDRS) subfield-specific optical retinal layer thicknesses were investigated in 7384 healthy eyes of 7384 participants from the LIFE-Adult-Study stratified by normal glucose tolerance, prediabetes (impaired fasting glucose and/or impaired glucose tolerance and/or HbA1c 5.7-6.4% [39-47 mmol/mol]) and diabetes. The association of optical retinal band characteristics with different indices of glucose tolerance (e.g. fasting glucose, area under the glucose curve), insulin resistance (e.g. HOMA2-IR, triglyceride glucose index), or insulin sensitivity (e.g. estimated glucose disposal rate [eGDR], Stumvoll metabolic clearance rate) was determined using multivariable linear regression analyses for the individual markers adjusted for age, sex and refraction. Various sensitivity analyses were performed to validate the observed findings.

RESULTS

In the study cohort, nine out of ten optical bands of the retina showed significant sex- and glucose tolerance-dependent differences in band thicknesses. Multivariable linear regression analyses revealed a significant, independent, and inverse association between markers of glucose intolerance and insulin resistance (e.g. HOMA2-IR) with the thickness of the optical bands representing the anatomical retinal outer nuclear layer (ONL, standardised β=-0.096; p<0.001 for HOMA2-IR) and myoid zone (MZ; β=-0.096; p<0.001 for HOMA2-IR) of the photoreceptors. Conversely, markers of insulin sensitivity (e.g. eGDR) positively and independently associated with ONL (β=0.090; p<0.001 for eGDR) and MZ (β=0.133; p<0.001 for eGDR) band thicknesses. These global associations were confirmed in ETDRS subfield-specific analyses. Sensitivity analyses further validated our findings when physical activity, neuroanatomical cell/tissue types and ETDRS subfield categories were investigated after stratifying the cohort by glucose homeostasis.

CONCLUSIONS/INTERPRETATION

An impaired glucose homeostasis associates with a thinning of the optical bands of retinal ONL and photoreceptor MZ. Changes in ONL and MZ thicknesses might predict early metabolic retinal alterations in diabetes.

The Association of Race, Ethnicity, and Insurance Status with the Visual Acuity of Retinoblastoma Survivors in the IRIS® Registry

PURPOSE

To identify sociodemographic factors associated with the visual outcomes of retinoblastoma survivors.

METHODS

Retrospective cohort study using a US-based clinical data registry. All individuals < 18 years of age with a history of retinoblastoma in the Intelligent Research in Sight (IRIS®) Registry (1/1/2013-12/31/2020). The primary outcome was visual acuity below the threshold for legal blindness (20/200 or worse) in at least one eye. Multivariable logistic regression was used to evaluate the association between visual outcomes and age, sex, laterality, race, ethnicity, type of insurance, and geographic location.

RESULTS

This analysis included 1545 children with a history of retinoblastoma. The median length of follow-up was 4.1 years (IQR, 2.2-5.9 years) and the median age at most recent clinical visit was 12 years (IQR, 8-16 years). Retinoblastoma was unilateral in 54% of cases. Poor vision in at least one eye was identified in 78% of all children and poor vision in both eyes in 17% of those with bilateral disease. Poor visual outcomes were associated with unilateral diagnosis (OR, 1.55; 95% CI,1.13-2.12; p = .007), Black race (OR, 2.03; 95% CI, 1.19-3.47; p = .010), Hispanic ethnicity (OR, 1.65; 95% CI, 1.16-2.37; p = .006), and non-private insurance (OR, 1.47; 95% CI, 1.02-2.10; p = .037).

CONCLUSIONS

Poor visual outcomes appear to be more common among Black, Hispanic, and publicly insured children with a history of retinoblastoma, raising concerns regarding healthcare inequities. Primary care physicians should ensure that young children receive red reflex testing during routine visits and consider retinoblastoma in the differential diagnosis of abnormal eye exams.

Polygenic Risk Scores for Glaucoma Onset in the Ocular Hypertension Treatment Study

Importance

Primary open-angle glaucoma (POAG) is a highly heritable disease, with 127 identified risk loci to date. Polygenic risk score (PRS) may provide a clinically useful measure of aggregate genetic burden and improve patient risk stratification.

Objective

To assess whether a PRS improves prediction of POAG onset in patients with ocular hypertension.

Design, Setting, and Participants

This was a post hoc analysis of the Ocular Hypertension Treatment Study. Data were collected from 22 US sites with a mean (SD) follow-up of 14.0 (6.9) years. A total of 1636 participants were followed up from February 1994 to December 2008; 1077 participants were enrolled in an ancillary genetics study, of which 1009 met criteria for this analysis. PRS was calculated using summary statistics from the largest cross-ancestry POAG meta-analysis, with weights trained using 8 813 496 variants from 449 186 cross-ancestry participants in the UK Biobank. Data were analyzed from July 2022 to December 2023.

Exposures

From February 1994 to June 2002, participants were randomized to either topical intraocular pressure-lowering medication or close observation. After June 2002, both groups received medication.

Main Outcomes and Measures

Outcome measures were hazard ratios for POAG onset. Concordance index and time-dependent areas under the receiver operating characteristic curve were used to compare the predictive performance of multivariable Cox proportional hazards models.

Results

Of 1009 included participants, 562 (55.7%) were female, and the mean (SD) age was 55.9 (9.3) years. The mean (SD) PRS was significantly higher for 350 POAG converters (0.24 [0.95]) compared with 659 nonconverters (-0.12 [1.00]) (P < .001). POAG risk increased 1.36% (95% CI, 1.08-1.64) with each higher PRS decile, with conversion ranging from 9.52% (95% CI, 7.09-11.95) in the lowest PRS decile to 21.81% (95% CI, 19.37-24.25) in the highest decile. Comparison of low-risk and high-risk PRS tertiles showed a 2.0-fold increase in 20-year POAG risk for participants of European and African ancestries. In the subgroup randomized to delayed treatment, each increase in PRS decile was associated with a 0.52-year (95% CI, 0.01-1.03) decrease in age at diagnosis (P = .047). No significant linear association between PRS and age at POAG diagnosis was present in the early treatment group. Prediction models significantly improved with the addition of PRS as a covariate (C index = 0.77) compared with the Ocular Hypertension Treatment Study baseline model (C index = 0.75) (P < .001). Each 1-SD higher PRS conferred a mean hazard ratio of 1.25 (95% CI, 1.13-1.44) for POAG onset.

Conclusions and Relevance

Higher PRS was associated with increased risk for POAG in patients with ocular hypertension. The inclusion of a PRS improved the prediction of POAG onset.

Trial Registration

ClinicalTrials.gov Identifier: NCT00000125.

Comparison of Structural and Functional Features in Primary Angle Closure and Open Angle Glaucomas

PRCIS

Using a large data set, we showed structural and functional differences between primary angle closure glaucoma (PACG) and primary open angle glaucoma (POAG). Primary angle closure glaucoma has relative structural preservation and worse functional loss inferiorly.

PURPOSE

To identify structural and functional differences in PACG and POAG.

MATERIALS AND METHODS

In this large cross-sectional study, differences in structural and functional damage were assessed among patients with POAG and PACG with optical coherence tomography and reliable visual field testing.

RESULTS

In all, 283 patients with PACG and 4110 patients with POAG were included. Despite similar mean deviation on visual fields (mean [SD] -7.73 [7.92] vs. -7.53 [6.90] dB, P =0.72), patients with PACG had thicker global retinal nerve fiber layer (RNFL), smaller cup volume, smaller cup-to-disc ratio, and larger rim area than POAG (77 [20] vs. 71 [14] µm, 0.32 [0.28] vs. 0.40 [0.29] mm 3 , 0.6 [0.2] vs. 0.7 [0.1], 1.07 [0.40] vs. 0.89 [0.30] mm 2 , P <0.001 for all), while patients with POAG had more pronounced inferior RNFL thinning (82 [24] vs. 95 [35] µm, P <0.001). In a multivariable analysis, hyperopia [odds ratio (OR): 1.24, confidence interval (CI): 1.13-1.37], smaller cup-to-disc ratio (OR: 0.69, CI: 0.61-0.78), thicker inferior RNFL (OR: 1.15, CI: 1.06-1.26) and worse mean deviation (OR: 0.95, CI: 0.92-0.98) were associated with PACG. Functionally, POAG was associated with superior paracentral loss and PACG with inferior field loss. After adjusting for average RNFL thickness, PACG was associated with more diffuse loss than POAG (total deviation differences 1.26-3.2 dB).

CONCLUSIONS

Patients with PACG had less structural damage than patients with POAG despite similar degrees of functional loss. Regional differences in patterns of functional and structural loss between POAG and PACG may improve disease monitoring for these glaucoma subtypes.

Association of retinal optical coherence tomography metrics and polygenic risk scores with cognitive function and future cognitive decline

PURPOSE

To evaluate the potential of retinal optical coherence tomography (OCT) measurements and polygenic risk scores (PRS) to identify people at risk of cognitive impairment.

METHODS

Using OCT images from 50 342 UK Biobank participants, we examined associations between retinal layer thickness and genetic risk for neurodegenerative disease and combined these metrics with PRS to predict baseline cognitive function and future cognitive deterioration. Multivariate Cox proportional hazard models were used to predict cognitive performance. P values for retinal thickness analyses are false-discovery-rate-adjusted.

RESULTS

Higher Alzheimer's disease PRS was associated with a thicker inner nuclear layer (INL), chorio-scleral interface (CSI) and inner plexiform layer (IPL) (all p<0.05). Higher Parkinson's disease PRS was associated with thinner outer plexiform layer (p<0.001). Worse baseline cognitive performance was associated with thinner retinal nerve fibre layer (RNFL) (aOR=1.038, 95% CI (1.029 to 1.047), p<0.001) and photoreceptor (PR) segment (aOR=1.035, 95% CI (1.019 to 1.051), p<0.001), ganglion cell complex (aOR=1.007, 95% CI (1.002 to 1.013), p=0.004) and thicker ganglion cell layer (aOR=0.981, 95% CI (0.967 to 0.995), p=0.009), IPL (aOR=0.976, 95% CI (0.961 to 0.992), p=0.003), INL (aOR=0.923, 95% CI (0.905 to 0.941), p<0.001) and CSI (aOR=0.998, 95% CI (0.997 to 0.999), p<0.001). Worse future cognitive performance was associated with thicker IPL (aOR=0.945, 95% CI (0.915 to 0.999), p=0.045) and CSI (aOR=0.996, 95% CI (0.993 to 0.999) 95% CI, p=0.014). Prediction of cognitive decline was significantly improved with the addition of PRS and retinal measurements.

CONCLUSIONS AND RELEVANCE

Retinal OCT measurements are significantly associated with genetic risk of neurodegenerative disease and may serve as biomarkers predictive of future cognitive impairment.

Effectiveness of angle based minimally invasive glaucoma surgery after laser trabeculoplasty: an analysis of the IRIS® Registry

OBJECTIVE

Angle-based minimally invasive glaucoma surgery (ab-MIGS) has grown substantially, although long-term efficacy is poorly understood. We analyze ab-MIGS effectiveness with/without preceding laser trabeculoplasty (LTP).

DESIGN

Retrospective cohort study SUBJECTS: Eyes undergoing ab-MIGS (Canaloplasty, Goniotomy, Trabectome, iStent) with/without prior SLT (<2 years preceding MIGS) were identified in the IRIS® Registry (Intelligent Research in Sight) 2013-2018.

METHODS

Propensity score matching (PSM) was undertaken to define 4 cohorts: (i) standalone ab-MIGS, no prior LTP vs (ii) standalone ab-MIGS, with prior LTP; and (iii) ab-MIGS + phacoemulsification, no prior LTP vs (iv) ab-MIGS + phacoemulsification, with prior LTP.

MAIN OUTCOME MEASURES

Failure was defined as subsequent glaucoma reoperation following ab-MIGS (either MIGS or traditional glaucoma surgery). Time to event outcome and incidence rates were calculated using survival analysis, and adjusted hazard ratios (aHR) generated using multivariate Cox proportional hazards models. Medication data was not available for analysis.

RESULTS

164,965 unique MIGS procedures were performed, from 2013-2018. After PSM, we identified 954 eyes undergoing standalone ab-MIGS and 7,522 undergoing ab-MIGS + phacoemulsification. For eyes undergoing standalone ab-MIGS, those with prior LTP (n=477) were more likely to undergo reoperation vs those without LTP (n=477) at 6 and 12 months. In multivariable models, those with prior LTP were more likely to undergo reoperation over the 36-month period vs those without prior LTP (aHR 1.53 (CI 1.15-2.04), p=0.004). For eyes undergoing ab-MIGS + phacoemulsification, those with prior LTP (n=3,761) were more likely to undergo reoperation vs those without LTP (n=3,761) at 12, 24 and 36 months. In multivariable models, those with prior LTP were more likely to undergo reoperation over the 36-month period vs those without prior LTP (aHR 1.53 (CI 1.15-2.04), p=0.004).

CONCLUSIONS

Prior LTP may be associated with a higher chance of subsequent glaucoma surgery following ab-MIGS, either with or without concurrent phacoemulsification. These findings have important implications in understanding who may benefit most from ab-MIGS, and guiding patient and surgeon treatment expectations.

Harvard Glaucoma Fairness: A Retinal Nerve Disease Dataset for Fairness Learning and Fair Identity Normalization

Fairness (also known as equity interchangeably) in machine learning is important for societal wellbeing, but limited public datasets hinder its progress. Currently, no dedicated public medical datasets with imaging data for fairness learning are available, though minority groups suffer from more health issues. To address this gap, we introduce Harvard Glaucoma Fairness (Harvard-GF), a retinal nerve disease dataset including 3,300 subjects with both 2D and 3D imaging data and balanced racial groups for glaucoma detection. Glaucoma is the leading cause of irreversible blindness globally with Blacks having doubled glaucoma prevalence than other races. We also propose a fair identity normalization (FIN) approach to equalize the feature importance between different identity groups. Our FIN approach is compared with various state-of-the-art fairness learning methods with superior performance in the racial, gender, and ethnicity fairness tasks with 2D and 3D imaging data, demonstrating the utilities of our dataset Harvard-GF for fairness learning. To facilitate fairness comparisons between different models, we propose an equity-scaled performance measure, which can be flexibly used to compare all kinds of performance metrics in the context of fairness. The dataset and code are publicly accessible via https://ophai.hms.harvard.edu/harvard-gf3300/.

An Artificial Intelligence Enabled System for Retinal Nerve Fiber Layer Thickness Damage Severity Staging

Purpose

To develop an objective glaucoma damage severity classification system based on OCT-derived retinal nerve fiber layer (RNFL) thickness measurements.

Design

Algorithm development for RNFL damage severity classification based on multicenter OCT data.

Subjects and Participants

A total of 6561 circumpapillary RNFL profiles from 2269 eyes of 1171 subjects to develop models, and 2505 RNFL profiles from 1099 eyes of 900 subjects to validate models.

Methods

We developed an unsupervised k-means model to identify clusters of eyes with similar RNFL thickness profiles. We annotated the clusters based on their respective global RNFL thickness. We computed the optimal global RNFL thickness thresholds that discriminated different severity levels based on Bayes' minimum error principle. We validated the proposed pipeline based on an independent validation dataset with 2505 RNFL profiles from 1099 eyes of 900 subjects.

Main Outcome Measures

Accuracy, area under the receiver operating characteristic curve, and confusion matrix.

Results

The k-means clustering discovered 4 clusters with 1382, 1613, 1727, and 1839 samples with mean (standard deviation) global RNFL thickness of 58.3 (8.9) μm, 78.9 (6.7) μm, 87.7 (8.2) μm, and 101.5 (7.9) μm. The Bayes' minimum error classifier identified optimal global RNFL values of > 95 μ m , 86 to 95 μ m , 70 to 85 μ m , and < 70 μ m for discriminating normal eyes and eyes at the early, moderate, and advanced stages of RNFL thickness loss, respectively. About 4% of normal eyes and 98% of eyes with advanced RNFL loss had either global, or ≥ 1 quadrant, RNFL thickness outside of normal limits provided by the OCT instrument.

Conclusions

Unsupervised machine learning discovered that the optimal RNFL thresholds for separating normal eyes and eyes with early, moderate, and advanced RNFL loss were 95 μ m , 85 μm, and 70 μ m , respectively. This RNFL loss classification system is unbiased as there was no preassumption or human expert intervention in the development process. Additionally, it is objective, easy to use, and consistent, which may augment glaucoma research and day-to-day clinical practice.

Financial Disclosures

Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Risk Factors for Glaucoma Diagnosis and Surgical Intervention following Pediatric Cataract Surgery in the IRIS® Registry

PURPOSE

To compare demographic and clinical factors associated with glaucoma following cataract surgery (GFCS) and glaucoma surgery rates between infants, toddlers, and older children using a large, ophthalmic registry.

DESIGN

Retrospective cohort study.

PARTICIPANTS

Patients in the IRIS® Registry (Intelligent Research in Sight) who underwent cataract surgery at ≤ 17 years old and between January 1, 2013 and December 31, 2020.

METHODS

Glaucoma diagnosis and procedural codes were extracted from the electronic health records of practices participating in the IRIS Registry. Children with glaucoma diagnosis or surgery before cataract removal were excluded. The Kaplan-Meier estimator was used to determine the cumulative probability of GFCS diagnosis and glaucoma surgery after cataract surgery. Multivariable Cox regression was used to identify factors associated with GFCS and glaucoma surgery.

MAIN OUTCOME MEASURES

Cumulative probability of glaucoma diagnosis and surgical intervention within 5 years after cataract surgery.

RESULTS

The study included 6658 children (median age, 10.0 years; 46.2% female). The 5-year cumulative probability of GFCS was 7.1% (95% confidence interval [CI], 6.1%-8.1%) and glaucoma surgery was 2.6% (95% CI, 1.9%-3.2%). The 5-year cumulative probability of GFCS for children aged < 1 year was 22.3% (95% CI, 15.7%-28.4%). Risk factors for GFCS included aphakia (hazard ratio [HR], 2.63; 95% CI, 1.96-3.57), unilateral cataract (HR, 1.48; 95% CI, 1.12-1.96), and Black race (HR, 1.61; 95% CI, 1.12-2.32). The most common surgery was glaucoma drainage device insertion (32.6%), followed by angle surgery (23.3%), cyclophotocoagulation (15.1%), and trabeculectomy (5.8%).

CONCLUSIONS

Glaucoma following cataract surgery diagnosis in children in the IRIS Registry was associated with young age, aphakia, unilateral cataract, and Black race. Glaucoma drainage device surgery was the preferred surgical treatment, consistent with the World Glaucoma Association 2013 consensus recommendations for GFCS management.

FINANCIAL DISCLOSURE(S)

Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

The Incidence of Strabismus After Upper and Lower Blepharoplasty in the United States

PURPOSE

To compare the incidence of strabismus after upper and lower blepharoplasty in the United States.

METHODS

Retrospective cohort study of adults (age ≥18 years) in the IRIS Registry (Intelligent Research in Sight) who underwent blepharoplasty between January 1, 2013 and December 31, 2020. The primary outcome was the Kaplan-Meier estimated cumulative incidence of strabismus diagnosis and surgery within 3 years of blepharoplasty. Multivariable Cox regression was used to determine the association of blepharoplasty type with strabismus diagnosis and surgery, adjusting for patient age, sex, and geographic region.

RESULTS

Blepharoplasty was performed in 368,623 patients (median [interquartile range] age, 69 [63-75] years, and 69% female). Compared with those undergoing upper eyelid blepharoplasty, patients treated with lower eyelid blepharoplasty were slightly younger (median age, 66 vs. 69 years; p < 0.001) and more likely to be female (71% vs. 69%; p < 0.001). There was a greater 3-year incidence of strabismus diagnosis (2.0% vs. 1.5%; p < 0.001) and a greater 3-year incidence of strabismus surgery (0.15% vs. 0.06%; p = 0.003) for individuals undergoing lower vs. upper blepharoplasty. After adjusting for age, sex, and geographic region, lower blepharoplasty was associated with a higher 3-year risk of strabismus diagnosis (HR, 1.49; 95% CI, 1.23-1.81; p < 0.001) and surgery (HR, 2.53; 95% CI, 1.27-5.03; p = 0.008).

CONCLUSIONS

This registry-based analysis found that individuals undergoing lower eyelid blepharoplasty were at higher risk of strabismus compared with those undergoing upper eyelid blepharoplasty. Using large databases to understand the incidence of complications of frequently performed procedures may improve ophthalmologists' ability to provide data-driven counseling on surgical risks prior to intervention.

Factors associated with the use of botulinum toxin injections for adult strabismus in the IRIS Registry

This cross-sectional study used data from a large nationwide registry to describe the factors associated with use of botulinum toxin injections for adults with strabismus in the United States. Botulinum toxin injections were performed on 3.1% of adults undergoing an intervention for strabismus between 2013 and 2020. Adults treated with botulinum toxin injections were more likely to be older and female. Compared to non-Hispanic White patients, non-Hispanic Black patients were three times less likely to receive treatment with botulinum toxin after adjusting for age, sex, geographic region, and type of insurance. Efforts to understand the factors contributing to disparities in the use of botulinum toxin for strabismus may lead to opportunities for more equitable access to this intervention.

Fellow Eyes Conversion Rates in Patients With Unilateral Exudative Age-Related Macular Degeneration: An Academy IRIS® Registry Analysis

BACKGROUND AND OBJECTIVE

This study aimed to examine conversion rates from non-exudative to exudative age-related macular degeneration (AMD) in the fellow eye of patients with unilateral exudative AMD using the Academy IRIS® Registry.

PATIENTS AND METHODS

This study was a retrospective, cohort analysis from 2016 to 2019. Patient and disease characteristics including initial AMD stage were collected. Cox proportional-hazard (PH) and logistic regression modeling were performed.

RESULTS

The risk of conversion was lower for men relative to women and for Asians and Blacks relative to Whites. Compared to never-smokers, active smokers were at increased risk of conversion, and compared to initially early non-exudative AMD eyes, intermediate and advanced non-exudative AMD eyes had higher rates of conversion. Compared to active choroidal neovascularization eyes, eyes with inactive choroidal neovascularization and inactive scars had lower rates of fellow eye conversion.

CONCLUSIONS

In this cohort analysis of unilateral exudative AMD patients, women, Whites, and active smokers had higher rates of non-exudative to exudative AMD conversion in the fellow eye. [Ophthalmic Surg Lasers Imaging Retina 2024;55:xx-xx.].

Machine Learning-Derived Baseline Visual Field Patterns Predict Future Glaucoma Onset in the Ocular Hypertension Treatment Study

Purpose

The Ocular Hypertension Treatment Study (OHTS) identified risk factors for primary open-angle glaucoma (POAG) in patients with ocular hypertension, including pattern standard deviation (PSD). Archetypal analysis, an unsupervised machine learning method, may offer a more interpretable approach to risk stratification by identifying patterns in baseline visual fields (VFs).

Methods

There were 3272 eyes available in the OHTS. Archetypal analysis was applied using 24-2 baseline VFs, and model selection was performed with cross-validation. Decomposition coefficients for archetypes (ATs) were calculated. A penalized Cox proportional hazards model was implemented to select discriminative ATs. The AT model was compared to the OHTS model. Associations were identified between ATs with both POAG onset and VF progression, defined by mean deviation change per year.

Results

We selected 8494 baseline VFs. Optimal AT count was 19. The highest prevalence ATs were AT9, AT11, and AT7. The AT-based prediction model had a C-index of 0.75 for POAG onset. Multivariable models demonstrated that a one-interquartile range increase in the AT5 (hazard ratio [HR] = 1.14; 95% confidence interval [CI], 1.04-1.25), AT8 (HR = 1.22; 95% CI, 1.09-1.37), AT15 (HR = 1.26; 95% CI, 1.12-1.41), and AT17 (HR = 1.17; 95% CI, 1.03-1.31) coefficients conferred increased risk of POAG onset. AT5, AT10, and AT14 were significantly associated with rapid VF progression. In a subgroup analysis by high-risk ATs (>95th percentile or <75th percentile coefficients), PSD lost significance as a predictor of POAG in the low-risk group.

Conclusions

Baseline VFs, prior to detectable glaucomatous damage, contain occult patterns representing early changes that may increase the risk of POAG onset and VF progression in patients with ocular hypertension. The relationship between PSD and POAG is modified by the presence of high-risk patterns at baseline. An AT-based prediction model for POAG may provide more interpretable glaucoma-specific information in a clinical setting.

The Epidemiology and Risk Factors for the Progression of Sympathetic Ophthalmia in the United States: An IRIS Registry Analysis

PURPOSE

To investigate the demographic and clinical characteristics of patients with sympathetic ophthalmia (SO) and define the risk factors for its incidence following trauma and ophthalmic procedures.

DESIGN

Retrospective cohort study.

PARTICIPANTS

Patients in the American Academy of Ophthalmology's (Academy) IRIS Registry (Intelligent Research in Sight) who were (n=1523) or were not diagnosed with SO following a documented procedure or trauma between January 1, 2013, and December 31, 2019.

METHODS

Multiple demographic and clinical factors were collected, descriptive statistics and prevalence were calculated, and multivariate linear regression models were fit to the data.

MAIN OUTCOME MEASURES

Prevalence of SO, demographic and clinical characteristics, and beta coefficient (β) estimates of demographic and clinical characteristics impacting time to SO onset after procedure (Procedure Only cohort) or trauma (Trauma cohort).

RESULTS

Of 65,348,409 distinct IRIS Registry patients, 1523 (0.0023%) were diagnosed with SO between 2013 and 2019, and also had a documented preceding trauma or procedure. Of these, 927 (60.87%) were female, 1336 (87.72%) belonged to the Procedure Only cohort, and 187 (12.28%) belonged to the Trauma cohort. The prevalence of SO after trauma was 0.0207%, whereas after procedures it was 0.0124%. The highest risk of procedure-related SO was seen in patients with history of "other anterior segment" (0.122%) followed by glaucoma (0.066%) procedures, whereas the lowest prevalence was noted with cataract surgeries (0.011%). The average time to onset of SO across both cohorts combined was 527.44 (±715.60) days, with statistically significant differences between the 2 cohorts (P < .001). On average, the time to onset from inciting event to SO was shorter with increasing age, by 9.02 (95% CI: -11.96, -6.08) days for every 1-year increase.

CONCLUSIONS

SO following trauma and ophthalmic procedure is potentially rarer than previously reported, as measured in this large ophthalmic medical record database. Female sex may be a risk factor for SO. Older age may be a risk factor for quicker onset. These findings can guide clinical decision-making and management.

The epidemiology of pediatric dry eye disease in the United States: An IRIS® registry (Intelligent Research in Sight) analysis

PURPOSE

Dry-eye disease (DED) is a chronic progressive ocular surface disorder with limited studies in the pediatric population. The Academy of Ophthalmology's IRIS® Registry was leveraged to investigate the prevalence of DED in the pediatric population (PDED, patients <18 years old) and the demographic differences of DED between pediatric and adult patients (ADED).

METHODS

Retrospective cohort study. Patients with DED between January 1st, 2013 and December 31st, 2019 (N = 4,795,979) were included. Descriptive statistics, Pearson's chi-squared tests and two-sample proportions tests were conducted to compare key demographic distributions between the ADED and PDED cohorts.

RESULTS

The average age at onset for ADED patients was 61.06 (±14.75) years and for PDED patients was 12.51 (±3.86). The overall tests for independence and the individual tests of proportions of each category were statistically significant for all demographic characteristics (p < 0.001). Characteristics with the largest discrepancies between patients of PDED and the IRIS Registry pediatric patient pool (PIRIS) included female sex (58.08 % vs. 50.60 %), male sex (41.58 % vs. 48.78 %) and Asian race (6.02 % vs. 3.11 %) respectively. Within the PDED cohort, females were at higher risk of PDED (58 % vs. 42 %). PDED was more prevalent in children with refractive errors (76 %) and eyelid/conjunctival disorders (41 %). Characteristics with the largest discrepancies between PDED and ADED patients included female sex (58.08 % vs. 68.12 %), male sex (41.58 % vs. 31.55 %) and Caucasian race (50.24 % vs. 67.06 %) respectively.

CONCLUSIONS

Significant differences in the PDED cohort are demonstrated in this study. PDED was more prevalent in the female sex and Caucasian race compared to PIRIS and was more commonly associated with refractive errors and eyelid/conjunctival disorders.

Phenome- and genome-wide analyses of retinal optical coherence tomography images identify links between ocular and systemic health

The human retina is a multilayered tissue that offers a unique window into systemic health. Optical coherence tomography (OCT) is widely used in eye care and allows the noninvasive, rapid capture of retinal anatomy in exquisite detail. We conducted genotypic and phenotypic analyses of retinal layer thicknesses using macular OCT images from 44,823 UK Biobank participants. We performed OCT layer cross-phenotype association analyses (OCT-XWAS), associating retinal thicknesses with 1866 incident conditions (median 10-year follow-up) and 88 quantitative traits and blood biomarkers. We performed genome-wide association studies (GWASs), identifying inherited genetic markers that influence retinal layer thicknesses and replicated our associations among the LIFE-Adult Study (N = 6313). Last, we performed a comparative analysis of phenome- and genome-wide associations to identify putative causal links between retinal layer thicknesses and both ocular and systemic conditions. Independent associations with incident mortality were detected for thinner photoreceptor segments (PSs) and, separately, ganglion cell complex layers. Phenotypic associations were detected between thinner retinal layers and ocular, neuropsychiatric, cardiometabolic, and pulmonary conditions. A GWAS of retinal layer thicknesses yielded 259 unique loci. Consistency between epidemiologic and genetic associations suggested links between a thinner retinal nerve fiber layer with glaucoma, thinner PS with age-related macular degeneration, and poor cardiometabolic and pulmonary function with a thinner PS. In conclusion, we identified multiple inherited genetic loci and acquired systemic cardio-metabolic-pulmonary conditions associated with thinner retinal layers and identify retinal layers wherein thinning is predictive of future ocular and systemic conditions.

Surgical Approach and Reoperation Risk in Intermittent Exotropia in the IRIS Registry

Importance

There is no consensus on the optimal surgical treatment for children with intermittent exotropia (IXT).

Objective

To compare the 5-year reoperation rates for children with IXT treated with horizontal muscle strabismus surgery using bilateral lateral rectus recession (BLR) vs unilateral lateral rectus recession with medial rectus resection (RR).

Design, Setting, and Participants

This cohort study examined data obtained from the Intelligent Research in Sight (IRIS) Registry on 7482 children (age, <18 years) with IXT who underwent horizontal eye muscle strabismus surgery between January 1, 2013, and December 31, 2017. Children undergoing initial surgeries involving 3 or more horizontal muscles, vertical muscles, or reoperations were excluded.

Main Outcomes and Measures

The primary outcome was the adjusted cumulative incidence of repeat horizontal muscle surgery within 5 years after the initial surgery. Reoperation risk was analyzed using adjusted hazard ratios (AHRs) derived from multivariable Cox regression models, adjusting for individual demographic and surgical factors (age, sex, race and ethnicity, US Census region, and surgeon subspecialty). Data were analyzed between January 16 and September 20, 2023.

Results

The study included 7482 children (median [IQR] age at initial surgery, 6 [4-9] years; 3945 females [53%]) with IXT treated with horizontal muscle strabismus surgery. Bilateral lateral rectus recession was performed more frequently than RR (85.3% vs 14.7%, P < .001), especially in younger children (rates of BLR vs RR by age: age 0 to ≤4 years, 88.4% vs 11.6%; age 5 to ≤11 years, 84.7% vs 15.3%; age 12 to ≤17 years, 78.1% vs 21.9%; P < 0.001). After data adjustment, the 5-year cumulative incidence of reoperation was 21.3% (95% CI, 20.1%-22.5%). The adjusted 5-year cumulative incidence of reoperation was higher for BLR than for RR (22.2% vs 17.2%; difference, 4.9%; 95% CI, 1.9%-8.0%). Unilateral lateral rectus recession with medial rectus resection was associated with a lower 5-year reoperation risk compared with BLR (AHR, 0.77; 95% CI, 0.64-0.93). Younger age at time of initial surgery was associated with a higher reoperation risk (AHR per 1-year decrease, 1.09; 95% CI, 1.07-1.11) after adjusting for all other covariates.

Conclusions and Relevance

In this nationwide registry, approximately 1 in 5 children with IXT underwent reoperation within 5 years after the initial surgery. Children treated with RR were less likely to require a reoperation within 5 years compared with those treated with BLR. Further efforts to identify modifiable risk factors for reoperation are needed to reduce the surgical burden and improve outcomes for children with IXT.

Characterizing Macular Neovascularization in Myopic Macular Degeneration and Age-Related Macular Degeneration Using Swept Source OCTA

Purpose

Visual prognosis and treatment burden for macular neovascularization (MNV) can differ between myopic macular degeneration (MMD) and age-related macular degeneration (AMD). We describe and compare MNV associated with MMD and AMD using swept-source (SS)-OCTA.

Patients and Methods

Adult patients with documented MNV associated with MMD or AMD were consecutively recruited. Qualitative and quantitative features were assessed from 6x6mm angiograms, including the MNV area and vessel density (VD). Descriptive statistics and linear regression analyses were carried out.

Results

Out of 75 enrolled eyes with diagnosed MNV (30 MMD-MNV and 45 AMD-MNV; mean age 55±19 and 75±8 years, respectively), 44 eyes had discernible MNV (11 MMD-MNV and 33 AMD-MNV) on SS-OCTA at the time of the study and were included in the analysis. The MMD-MNV group exhibited a three-fold smaller sized MNV (p=0.001), lower greatest linear dimension (p=0.009) and greatest vascular caliber (p<0.001) compared to AMD-MNVs, and had a higher prevalence of tree-in-bud pattern. Eyes with AMD showed a higher prevalence of type 1 MNVs with medusa pattern. There was no difference in the location of the MNV, shape's regularity, margins, presence of core vessel, capillary fringe, peripheral loops, or perilesional dark halo (p>0.05) between both conditions. After adjustment, decreased MNV area and increased VD were associated with the tree-in-bud pattern, whereas the diagnosis did not significantly influence those parameters.

Conclusion

While larger studies are warranted, this study is the first to describe and compare MMD-MNV and AMD-MNV using SS-OCTA, providing relevant clinical insight on MNV secondary to MMD and AMD. These findings also further validate OCTA as a powerful tool to detect and characterize MNV non-invasively.

Structure-function associations between contrast sensitivity and widefield swept-source optical coherence tomography angiography in diabetic macular edema

PURPOSE

To evaluate the relationship between contrast sensitivity (CS) and widefield swept-source optical coherence tomography angiography (WF SS-OCTA) vascular metrics in diabetic macular edema (DME) was the purpose.

METHODS

This prospectively enrolled cross-sectional observational study included 61 eyes of 48 patients that were tested with the quantitative CS function (qCSF) test on the same day as imaging with WF SS-OCTA (PLEX® Elite 9000, Carl Zeiss Meditec) 3 × 3, 6 × 6, and 12 × 12 mm scans. Outcomes included visual acuity (VA) and multiple qCSF metrics. Vascular metrics included vessel density (VD) and vessel skeletonized density (VSD) in the superficial (SCP) and deep capillary plexus (DCP) and whole retina (WR) and foveal avascular zone (FAZ) parameters. Mixed effects multivariable linear regression models controlling for age, lens status, and diabetic retinopathy stage were performed. Standardized beta coefficients were calculated by refitting the standardized data.

RESULTS

SS-OCTA metrics had a significant association with CS and VA. The effect size of OCTA metrics was larger on CS compared to VA. For example, the standardized beta coefficients for VSD and CS at 3 cpd (βSCP = 0.76, βDCP = 0.71, βWR = 0.72, p < 0.001) were larger than those for VA (βSCP =  - 0.55, p < 0.001; βDCP =  - 0.43, p = 0.004; βWR =  - 0.50, p < 0.001). On 6 × 6 mm images, AULCSF, CS at 3 cpd, and CS at 6 cpd were significantly associated with VD and VSD in all three slab types (SCP, DCP, and WR), while VA was not.

CONCLUSION

Structure-function associations in patients with DME leveraging the qCSF device suggest microvascular changes on WF SS-OCTA are associated with larger changes in contrast sensitivity than VA.

Participant Experience Using Novel Perimetry Tests to Monitor Glaucoma Progression

PRCIS

Participant surveys taken after using tablet-based and smart visual function analyzer (SVFA) perimetry tests suggest that patients may prefer novel perimetry tests over traditional visual field machines.

PURPOSE

Compare patient experience using the IMOvifa SVFA and the tablet-based Melbourne Rapid Fields visual field (VF) tests to the Humphrey Field 24-2 Swedish Interactive Threshold Algorithm Standard.

PATIENTS AND METHODS

Prospective observational cohort study on adult participants with diagnoses of glaucoma suspect, ocular hypertension, or glaucoma. Participants attended 2 study visits ~3 months apart. During the first visit, participants were trained to use the 2 novel perimeters, took 1 test on both new devices and the Humphrey Field Analyzer, then were surveyed. Participants received tablets and performed weekly tablet VF tests at home between study visits. At the final study visit, participants re-took the VF tests and completed the same surveys.

RESULTS

Eighty-one participants were surveyed twice. At the baseline survey, participants preferred the SVFA (71.7%) and tablet tests (69.2%) over the Humphrey VF. Most were willing to perform weekly monitoring at home on the SVFA (69.1%) and tablet (75.4%). Participants generally had a "very good" overall experience when testing on the SVFA (71.6%) and tablet (90.1%). At the final visit, fewer participants were willing to test on the tablet daily (23.5% to 9.9%; P = 0.02 for change) and more were willing to test monthly (18.5% to 33.3%; P = 0.03 for change).

CONCLUSION

Users reported a preference for novel VF devices. Overall participant experience using these devices was positive, supporting the feasibility of home monitoring of VFs from an experience perspective.

Deep Learning for Localized Detection of Optic Disc Hemorrhages

PURPOSE

To develop an automated deep learning system for detecting the presence and location of disc hemorrhages in optic disc photographs.

DESIGN

Development and testing of a deep learning algorithm.

METHODS

Optic disc photos (597 images with at least 1 disc hemorrhage and 1075 images without any disc hemorrhage from 1562 eyes) from 5 institutions were classified by expert graders based on the presence or absence of disc hemorrhage. The images were split into training (n = 1340), validation (n = 167), and test (n = 165) datasets. Two state-of-the-art deep learning algorithms based on either object-level detection or image-level classification were trained on the dataset. These models were compared to one another and against 2 independent glaucoma specialists. We evaluated model performance by the area under the receiver operating characteristic curve (AUC). AUCs were compared with the Hanley-McNeil method.

RESULTS

The object detection model achieved an AUC of 0.936 (95% CI = 0.857-0.964) across all held-out images (n = 165 photographs), which was significantly superior to the image classification model (AUC = 0.845, 95% CI = 0.740-0.912; P = .006). At an operating point selected for high specificity, the model achieved a specificity of 94.3% and a sensitivity of 70.0%, which was statistically indistinguishable from an expert clinician (P = .7). At an operating point selected for high sensitivity, the model achieves a sensitivity of 96.7% and a specificity of 73.3%.

CONCLUSIONS

An autonomous object detection model is superior to an image classification model for detecting disc hemorrhages, and performed comparably to 2 clinicians.

Artifact Correction in Retinal Nerve Fiber Layer Thickness Maps Using Deep Learning and Its Clinical Utility in Glaucoma

Purpose

Correcting retinal nerve fiber layer thickness (RNFLT) artifacts in glaucoma with deep learning and evaluate its clinical usefulness.

Methods

We included 24,257 patients with optical coherence tomography and reliable visual field (VF) measurements within 30 days and 3,233 patients with reliable VF series of at least five measurements over ≥4 years. The artifacts are defined as RNFLT less than the known floor value of 50 µm. We selected 27,319 high-quality RNFLT maps with an artifact ratio (AR) of <2% as the ground truth. We created pseudo-artifacts from 21,722 low-quality RNFLT maps with AR of >5% and superimposed them on high-quality RNFLT maps to predict the artifact-free ground truth. We evaluated the impact of artifact correction on the structure-function relationship and progression forecasting.

Results

The mean absolute error and Pearson correlation of the artifact correction were 9.89 µm and 0.90 (P < 0.001), respectively. Artifact correction improved R2 for VF prediction in RNFLT maps with AR of >10% and AR of >20% up to 0.03 and 0.04 (P < 0.001), respectively. Artifact correction improved (P < 0.05) the AUC for progression prediction in RNFLT maps with AR of ≤10%, >10%, and >20%: (1) total deviation pointwise progression: 0.68 to 0.69, 0.62 to 0.63, and 0.62 to 0.64; and (2) mean deviation fast progression: 0.67 to 0.68, 0.54 to 0.60, and 0.45 to 0.56.

Conclusions

Artifact correction for RNFLTs improves VF and progression prediction in glaucoma.

Translational Relevance

Our model improves clinical usability of RNFLT maps with artifacts.

Chatbots Vs. Human Experts: Evaluating Diagnostic Performance of Chatbots in Uveitis and the Perspectives on AI Adoption in Ophthalmology

PURPOSE

To assess the diagnostic performance of two chatbots, ChatGPT and Glass, in uveitis diagnosis compared to renowned uveitis specialists, and evaluate clinicians' perception about utilizing artificial intelligence (AI) in ophthalmology practice.

METHODS

Six cases were presented to uveitis experts, ChatGPT (version 3.5 and 4.0) and Glass 1.0, and diagnostic accuracy was analyzed. Additionally, a survey about the emotions, confidence in utilizing AI-based tools, and the likelihood of incorporating such tools in clinical practice was done.

RESULTS

Uveitis experts accurately diagnosed all cases (100%), while ChatGPT achieved a diagnostic success rate of 66% and Glass 1.0 achieved 33%. Most attendees felt excited or optimistic about utilizing AI in ophthalmology practice. Older age and high level of education were positively correlated with increased inclination to adopt AI-based tools.

CONCLUSIONS

ChatGPT demonstrated promising diagnostic capabilities in uveitis cases and ophthalmologist showed enthusiasm for the integration of AI into clinical practice.

Normative Percentiles of Retinal Nerve Fiber Layer Thickness and Glaucomatous Visual Field Loss

Purpose

Circumpapillary retinal nerve fiber layer thickness (RNFLT) measurement aids in the clinical diagnosis of glaucoma. Spectral domain optical coherence tomography (SD-OCT) machines measure RNFLT and provide normative color-coded plots. In this retrospective study, we investigate whether normative percentiles of RNFLT (pRNFLT) from Spectralis SD-OCT improve prediction of glaucomatous visual field loss over raw RNFLT.

Methods

A longitudinal database containing OCT scans and visual fields from Massachusetts Eye & Ear glaucoma clinic patients was generated. Reliable OCT-visual field pairs were selected. Spectralis OCT normative distributions were extracted from machine printouts. Supervised machine learning models compared predictive performance between pRNFLT and raw RNFLT inputs. Regional structure-function associations were assessed with univariate regression to predict mean deviation (MD). Multivariable classification predicted MD, pattern standard deviation, MD change per year, and glaucoma hemifield test.

Results

There were 3016 OCT-visual field pairs that met the reliability criteria. Spectralis norms were found to be independent of age, sex, and ocular magnification. Regional analysis showed significant decrease in R2 from pRNFLT models compared to raw RNFLT models in inferotemporal sectors, across multiple regressors. In multivariable classification, there were no significant improvements in area under the curve of receiver operating characteristic curve (ROC-AUC) score with pRNFLT models compared to raw RNFLT models.

Conclusions

Our results challenge the assumption that normative percentiles from OCT machines improve prediction of glaucomatous visual field loss. Raw RNFLT alone shows strong prediction, with no models presenting improvement by the manufacturer norms. This may result from insufficient patient stratification in tested norms.

Translational Relevance

Understanding correlation of normative databases to visual function may improve clinical interpretation of OCT data.

Evaluating the Diagnostic Accuracy and Management Recommendations of ChatGPT in Uveitis

INTRODUCTION

Accurate diagnosis and timely management are vital for favorable uveitis outcomes. Artificial Intelligence (AI) holds promise in medical decision-making, particularly in ophthalmology. Yet, the diagnostic precision and management advice from AI-based uveitis chatbots lack assessment.

METHODS

We appraised diagnostic accuracy and management suggestions of an AI-based chatbot, ChatGPT, versus five uveitis-trained ophthalmologists, using 25 standard cases aligned with new Uveitis Nomenclature guidelines. Participants predicted likely diagnoses, two differentials, and next management steps. Comparative success rates were computed.

RESULTS

Ophthalmologists excelled (60-92%) in likely diagnosis, exceeding AI (60%). Considering fully and partially accurate diagnoses, ophthalmologists achieved 76-100% success; AI attained 72%. Despite an 8% AI improvement, its overall performance lagged. Ophthalmologists and AI agreed on diagnosis in 48% cases, with 91.6% exhibiting concurrence in management plans.

CONCLUSIONS

The study underscores AI chatbots' potential in uveitis diagnosis and management, indicating their value in reducing diagnostic errors. Further research is essential to enhance AI chatbot precision in diagnosis and recommendations.

Artifact-Tolerant Clustering-Guided Contrastive Embedding Learning for Ophthalmic Images in Glaucoma

Ophthalmic images, along with their derivatives like retinal nerve fiber layer (RNFL) thickness maps, play a crucial role in detecting and monitoring eye diseases such as glaucoma. For computer-aided diagnosis of eye diseases, the key technique is to automatically extract meaningful features from ophthalmic images that can reveal the biomarkers (e.g., RNFL thinning patterns) associated with functional vision loss. However, representation learning from ophthalmic images that links structural retinal damage with human vision loss is non-trivial mostly due to large anatomical variations between patients. This challenge is further amplified by the presence of image artifacts, commonly resulting from image acquisition and automated segmentation issues. In this paper, we present an artifact-tolerant unsupervised learning framework called EyeLearn for learning ophthalmic image representations in glaucoma cases. EyeLearn includes an artifact correction module to learn representations that optimally predict artifact-free images. In addition, EyeLearn adopts a clustering-guided contrastive learning strategy to explicitly capture the affinities within and between images. During training, images are dynamically organized into clusters to form contrastive samples, which encourage learning similar or dissimilar representations for images in the same or different clusters, respectively. To evaluate EyeLearn, we use the learned representations for visual field prediction and glaucoma detection with a real-world dataset of glaucoma patient ophthalmic images. Extensive experiments and comparisons with state-of-the-art methods confirm the effectiveness of EyeLearn in learning optimal feature representations from ophthalmic images.

Comparison of Perimetric Outcomes from a Tablet Perimeter, Smart Visual Function Analyzer, and Humphrey Field Analyzer

PURPOSE

The tablet-based Melbourne Rapid Fields (MRF) visual field (VF) test and the IMOvifa Smart Visual Function Analyzer (SVFA) are portable perimeters that may allow for at-home monitoring and more frequent testing. We compared tablet and SVFA results with outputs from the Humphrey Field Analyzer (HFA) 24-2 Swedish Interactive Threshold Algorithm Standard program.

DESIGN

Observational cross-sectional study.

SUBJECTS

Adult participants with a diagnosis of glaucoma, suspected glaucoma, or ocular hypertension seen in the Massachusetts Eye and Ear glaucoma clinic were enrolled. All participants were reliable and experienced HFA testers.

METHODS

Participants were tested with the SVFA and HFA. The study staff also trained participants on the MRF tablet with instructions to take weekly tests at home for 3 months. Visual field results from the 3 devices were compared.

MAIN OUTCOME MEASURES

Mean deviation (MD), pattern standard deviation (PSD), reliability parameters, and point sensitivity.

RESULTS

Overall, 79 participants (133 eyes) with a mean age of 61 ± 13 years (range, 26-79 years) were included; 59% of the participants were female, and the mean HFA MD was -2.7 ± 3.9 dB. The global indices of MD and PSD did not significantly vary between HFA and the 2 novel devices, except that the tablet VF reported a 0.6 dB higher PSD compared with HFA. However, tablet and SVFA sensitivities significantly differed from those of the HFA at 36 and 39 locations, respectively, out of 52 locations. Relative to HFA, the tablet overestimated light sensitivity in the nasal field while underestimating the temporal field. The SVFA generally underestimated light sensitivity, but its results were more similar to HFA results compared with the tablet.

CONCLUSIONS

Although average MD values from the 2 novel devices suggest that they provide similar results to the HFA, point-by-point comparisons highlight notable deviations. Differences in specific point sensitivity values were significant, especially between the tablet and the other 2 devices. These differences may in part be explained by differences in the devices' normative databases as well as how MD is calculated. However, the tablet had substantial differences based on location, indicating that the tablet design itself may be responsible for differences in local sensitivities.

FINANCIAL DISCLOSURE(S)

Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Deep Ocular Phenotyping Across Primary Open-Angle Glaucoma Genetic Burden

Importance

Better understanding of primary open-angle glaucoma (POAG) genetics could enable timely screening and promote individualized disease risk prognostication.

Objective

To evaluate phenotypic features across genetic burden for POAG.

Design, Setting, and Participants

This was a cross-sectional, population-based study conducted from 2006 to 2010. Included participants were individuals from the UK Biobank aged 40 to 69 years. Individuals with non-POAG forms of glaucoma were excluded from the analysis. Data were statistically analyzed from October 2022 to January 2023.

Main Outcomes and Measures

POAG prevalence based on structural coding, self-reports, and glaucoma-related traits.

Results

Among 407 667 participants (mean [SD] age, 56.3 [8.1] years; 219 183 majority sex [53.8%]) were 14 171 POAG cases. Area under receiver operating characteristic curve for POAG detection was 0.748 in a model including polygenic risk score (PRS), age, sex, and ancestry. POAG prevalence in the highest decile of PRS was 7.4% (3005 of 40 644) vs 1.3% (544 of 40 795) in lowest decile (P < .001). A 1-SD increase in PRS was associated with 1.74 times higher odds of POAG (95% CI, 1.71-1.77), a 0.61-mm Hg increase in corneal-compensated intraocular pressure (IOP; 95% CI, 0.59-0.64), a -0.09-mm Hg decrease in corneal hysteresis (95% CI, -0.10 to -0.08), a 0.08-mm Hg increase in corneal resistance factor (95% CI, 0.06-0.09), and a -0.08-diopter decrease in spherical equivalent (95% CI, -0.11 to -0.07; P < .001 for all). A 1-SD increase in PRS was associated with a thinning of the macula-region retinal nerve fiber layer (mRNFL) of 0.14 μm and macular ganglion cell complex (GCC) of 0.26 μm (P < .001 for both). In the subset of individuals with fundus photographs, a 1-SD increase in PRS was associated with 1.42 times higher odds of suspicious optic disc features (95% CI, 1.19-1.69) and a 0.013 increase in cup-disc ratio (CDR; 95% CI, 0.012-0.014; P < .001 for both). A total of 22 of 5193 fundus photographs (0.4%) in decile 10 had disc hemorrhages, and 27 of 5257 (0.5%) had suspicious optic disc features compared with 9 of 5158 (0.2%) and 10 of 5219 (0.2%), respectively, in decile 1 (P < .001 for both). CDR in decile 10 was 0.46 compared with 0.41 in decile 1 (P < .001).

Conclusion and Relevance

Results suggest that PRS identified a group of individuals at substantially higher risk for POAG. Higher genetic risk was associated with more advanced disease, namely higher CDR and corneal-compensated IOP, thinner mRNFL, and thinner GCC. Associations with POAG PRS and corneal hysteresis and greater prevalence of disc hemorrhages were identified. These results suggest that genetic risk is an increasingly important parameter for risk stratification to consider in clinical practice.

Automated detection of genetic relatedness from fundus photographs using Siamese Neural Networks

Heritability of common eye diseases and ocular traits are relatively high. Here, we develop an automated algorithm to detect genetic relatedness from color fundus photographs (FPs). We estimated the degree of shared ancestry amongst individuals in the UK Biobank using KING software. A convolutional Siamese neural network-based algorithm was trained to output a measure of genetic relatedness using 7224 pairs (3612 related and 3612 unrelated) of FPs. The model achieved high performance for prediction of genetic relatedness; when computed Euclidean distances were used to determine probability of relatedness, the area under the receiver operating characteristic curve (AUROC) for identifying related FPs reached 0.926. We performed external validation of our model using FPs from the LIFE-Adult study and achieved an AUROC of 0.69. An occlusion map indicates that the optic nerve and its surrounding area may be the most predictive of genetic relatedness. We demonstrate that genetic relatedness can be captured from FP features. This approach may be used to uncover novel biomarkers for common ocular diseases.

Polygenic Risk Score Improves Prediction of Primary Open Angle Glaucoma Onset in the Ocular Hypertension Treatment Study

Objective or Purpose

Primary open-angle glaucoma (POAG) is a highly heritable disease with 127 identified risk loci. Polygenic risks score (PRS) offers a measure of aggregate genetic burden. In this study, we assess whether PRS improves risk stratification in patients with ocular hypertension.

Design

A post-hoc analysis of the Ocular Hypertension Treatment Study (OHTS) data.

Setting Participants and/or Controls

1636 participants were followed from 1994 to 2020 across 22 sites. The PRS was computed for 1009 OHTS participants using summary statistics from largest cross-ancestry POAG metanalysis with weights trained using 8,813,496 variants from 488,395 participants in the UK Biobank.

Methods Interventions or Testing

Survival regression analysis, with endpoint as development of POAG, predicted disease onset from PRS incorporating baseline covariates.

Main Outcomes and Measures

Outcome measures were hazard ratios for POAG onset. Concordance index and time-dependent AUC were used to compare the predictive performance of multivariable Cox-Proportional Hazards models.

Results

Mean PRS was significantly higher for POAG-converters (0.24 ± 0.95) than for non-converters (-0.12 ± 1.00) (p < 0.01). POAG risk increased 1.36% with each higher PRS decile, with conversion ranging from 9.5% in the lowest PRS decile to 21.8% in the highest decile. Comparison of low- and high-risk PRS tertiles showed a 1.8-fold increase in 20-year POAG risk for participants of European and African ancestries (p<0.01). In the subgroup randomized to delayed treatment, each increase in PRS decile was associated with a 0.52-year decrease in age at diagnosis, (p=0.05). No significant linear relationship between PRS and age at POAG diagnosis was present in the early treatment group. Prediction models significantly improved with the addition of PRS as a covariate (C-index = 0.77) compared to OHTS baseline model (C-index=0.75) (p<0.01). One standard deviation higher PRS conferred a mean hazard ratio of 1.25 (CI=[1.13, 1.44]) for POAG onset.

Conclusions

Higher PRS is associated with increased risk for, and earlier development of POAG in patients with ocular hypertension. Early treatment may mitigate the risk from high genetic burden, delaying clinically detectable disease by up to 5.2 years. The inclusion of a PRS improves the prediction of POAG onset.

Metabolite and Lipid Biomarkers Associated With Intraocular Pressure and Inner Retinal Morphology: 1H NMR Spectroscopy Results From the UK Biobank

Purpose

The purpose of this study was to assess metabolites associated with intraocular pressure (IOP) and inner retina structure.

Methods

We cross-sectionally assessed 168 non-fasting plasma metabolites measured by nuclear magnetic resonance (NMR) spectroscopy with IOP (n = 28,195), macular retinal nerve fiber layer thickness (mRNFL; n = 10,584), and macular ganglion cell inner plexiform layer thickness (mGCIPL; n = 10,554) in the UK Biobank. We used multiple linear regression models adjusting for various covariates with probit-transformed metabolite levels as predictors for each outcome. Each estimate represents the difference in outcome variable per standard deviation increase in the probit-transformed metabolite values. We used the number of effective (NEF) tests and false discovery rate (FDR) to adjust for multiple comparisons for metabolites and metabolite classes, respectively.

Results

In individual metabolite analysis, multiple amino acids, especially branched-chain amino acids, were associated with lower IOP (-0.12 mm Hg; 95% confidence interval = -0.16 to -0.07; NEF = 2.7E-05). Albumin, 3 hydroxybutyrate, lactate, and several lipids were associated with higher IOP (range = 0.07 to 0.18 mm Hg, NEF = ≤ 0.039). In IOP-adjusted analyses, five HDL-related metabolites were associated with thinner mRNFL (-0.15 microns for all metabolites, NEF = ≤ 0.027), whereas five LDL-related metabolites were associated with thicker mGCIPL (range = 0.17 to 0.20 microns; NEF = ≤ 0.044). In metabolite class analysis, the lipid components of lipoproteins (cholesterol, triglycerides, etc.) were not associated with our outcomes (FDR > 0.2 for all); yet multiple lipoproteins were significantly (FDR < 0.05) associated with all outcomes.

Conclusions

Branched-chain amino acids were associated with lower IOP, HDL metabolites were associated with thinner mRNFL, and LDL metabolites were associated with thicker mGCIPL.

Impact of the Affordable Care Act on Glaucoma Severity at First Presentation

PURPOSE

To test whether the increase in insurance coverage in Massachusetts due to the Affordable Care Act (ACA) is associated with a decrease in glaucoma severity in patients presenting for the first time at a tertiary health system.

METHODS

Clinical and billing data of first-time glaucoma patients from a tertiary health system in Massachusetts from 2006 to 2021 was used. Pre-ACA is defined as before 2014 and post-ACA is defined as after 2014. Visual field mean deviation was used to define glaucoma severity: greater than -6 dB, less than -6 dB and greater than -12 dB, and less than -12 dB was classified as "mild," "moderate," and "severe" respectively. Ordinal logistic regression models adjusted for age, race, gender, and insurance type were used to determine the odds of presenting with more severe glaucoma.

RESULTS

2,394 pre-ACA and 3,651 post-ACA first-time glaucoma patients were identified. There was no significant difference in the likelihood of more severe glaucoma at first presentation post-ACA compared to pre-ACA (OR=0.96; 95% CI 0.86-1.08; p=0.49) among the entire population. In stratified analyses, patients who utilized Medicaid for insurance had 52% decreased odds for presenting with more severe glaucoma at first presentation post-ACA compared to pre-ACA (OR=0.48; 95% CI 0.33-0.69; p<0.001). This remained significant after adjustment for age, race, and gender (adjusted OR=0.44; 95% CI 0.29-0.65; p<0.001).

CONCLUSION

At a Massachusetts-based tertiary healthcare center, individuals on Medicaid were more likely to have more severe glaucoma at first presentation before the implementation of the ACA, compared to after.

Plasma metabolite profile for primary open-angle glaucoma in three US cohorts and the UK Biobank

Glaucoma is a progressive optic neuropathy and a leading cause of irreversible blindness worldwide. Primary open-angle glaucoma is the most common form, and yet the etiology of this multifactorial disease is poorly understood. We aimed to identify plasma metabolites associated with the risk of developing POAG in a case-control study (599 cases and 599 matched controls) nested within the Nurses' Health Studies, and Health Professionals' Follow-Up Study. Plasma metabolites were measured with LC-MS/MS at the Broad Institute (Cambridge, MA, USA); 369 metabolites from 18 metabolite classes passed quality control analyses. For comparison, in a cross-sectional study in the UK Biobank, 168 metabolites were measured in plasma samples from 2,238 prevalent glaucoma cases and 44,723 controls using NMR spectroscopy (Nightingale, Finland; version 2020). Here we show higher levels of diglycerides and triglycerides are adversely associated with glaucoma in all four cohorts, suggesting that they play an important role in glaucoma pathogenesis.

Insights into human health from phenome- and genome-wide analyses of UK Biobank retinal optical coherence tomography phenotypes

The human retina is a complex multi-layered tissue which offers a unique window into systemic health and disease. Optical coherence tomography (OCT) is widely used in eye care and allows the non-invasive, rapid capture of retinal measurements in exquisite detail. We conducted genome- and phenome-wide analyses of retinal layer thicknesses using macular OCT images from 44,823 UK Biobank participants. We performed phenome-wide association analyses, associating retinal thicknesses with 1,866 incident ICD-based conditions (median 10-year follow-up) and 88 quantitative traits and blood biomarkers. We performed genome-wide association analyses, identifying inherited genetic markers which influence the retina, and replicated our associations among 6,313 individuals from the LIFE-Adult Study. And lastly, we performed comparative association of phenome- and genome- wide associations to identify putative causal links between systemic conditions, retinal layer thicknesses, and ocular disease. Independent associations with incident mortality were detected for photoreceptor thinning and ganglion cell complex thinning. Significant phenotypic associations were detected between retinal layer thinning and ocular, neuropsychiatric, cardiometabolic and pulmonary conditions. Genome-wide association of retinal layer thicknesses yielded 259 loci. Consistency between epidemiologic and genetic associations suggested putative causal links between thinning of the retinal nerve fiber layer with glaucoma, photoreceptor segment with AMD, as well as poor cardiometabolic and pulmonary function with PS thinning, among other findings. In conclusion, retinal layer thinning predicts risk of future ocular and systemic disease. Furthermore, systemic cardio-metabolic-pulmonary conditions promote retinal thinning. Retinal imaging biomarkers, integrated into electronic health records, may inform risk prediction and potential therapeutic strategies.

Archetypal analysis of longitudinal visual fields for idiopathic intracranial hypertension patients presenting in a clinic setting

We previously applied archetypal analysis (AA) using visual fields (VF) from the Idiopathic Intracranial Hypertension Treatment Trial (IIHTT) to derive a model, which quantified patterns (or archetypes [ATs] of VF loss), anticipated recovery, and identified residual VF deficits. We hypothesized that AA could produce similar results using IIH VFs collected in clinical practice. We applied AA to 803 VFs from 235 eyes with IIH from an outpatient neuro-ophthalmology clinic and created a clinic-derived model of ATs, with the relative weight (RW) and average total deviation (TD) for each AT. We also created a combined-derived model from an input dataset containing the clinic VFs and 2862 VFs from the IIHTT. We used both models to decompose clinic VF into ATs of varying percent weight (PW), correlated presentation AT PW with mean deviation (MD), and evaluated final visit VFs considered "normal" by MD ≥ -2.00 dB for residual abnormal ATs. The 14-AT clinic-derived and combined-derived models revealed similar patterns of VF loss previously identified in the IIHTT model. AT1 (a normal pattern) was most prevalent in both models (RW = 51.8% for clinic-derived; 35.4% for combined-derived). Presentation AT1 PW correlated with final visit MD (r = 0.82, p < 0.001 for the clinic-derived model; r = 0.59, p < 0.001 for the combined-derived model). Both models showed ATs with similar patterns of regional VF loss. The most common patterns of VF loss in "normal" final visit VFs using each model were clinic-derived AT2 (mild global depression with enlarged blind spot; 44/125 VFs; 34%) and combined-derived AT2 (near-normal; 93/149 VFs; 62%). AA provides quantitative values for IIH-related patterns of VF loss that can be used to monitor VF changes in a clinic setting. Presentation AT1 PW is associated with the degree of VF recovery. AA identifies residual VF deficits not otherwise indicated by MD.

Smoking Is Associated With a Higher Risk of Surgical Intervention for Thyroid Eye Disease in the IRIS Registry

PURPOSE

To describe the association of smoking status with surgical intervention for thyroid eye disease (TED) at the population-level.

DESIGN

Retrospective cohort study.

METHODS

This study included all adults (aged ≥18 years) with Graves disease in the Intelligent Research in Sight (IRIS) Registry (January 1, 2013, to December 31, 2020). The primary outcome was surgical intervention for TED, stratified into orbital decompression, strabismus surgery, and eyelid recession surgery. The Kaplan-Meier estimated 5-year cumulative probability for each surgical intervention was calculated. Multivariable Cox regression was used to evaluate the association between smoking status and each surgical intervention, adjusting for age, sex, race, ethnicity, and geographic region.

RESULTS

This study included 87,774 patients. Median age was 59 years (IQR, 48-68 years); 81% were female patients. Current smokers had a greater 5-year cumulative probability of orbital decompression (3.7% vs 1.9%; P < .001), strabismus surgery (4.6% vs 2.2%; P < .001), and eyelid recession (4.1% vs 2.6%; P < .001) compared to never smokers. After adjusting for demographic factors, current smokers were at greater risk for orbital decompression (hazard ratio [HR], 2.1; 95% CI, 1.8-2.4; P < .001), strabismus surgery (HR, 2.0; 95% CI, 1.8-2.3; P < .001), and eyelid recession (HR, 1.7; 95% CI, 1.5-1.9; P < .001) than never smokers. Former smokers were at higher risk for each type of surgery for TED, albeit at lower levels than current smokers.

CONCLUSIONS

Smoking was associated with increased risk of surgical intervention for TED in the IRIS Registry. Former smokers were at a lower risk than current smokers, supporting the role of smoking cessation on lowering the burden of surgical disease at the population-level.

Factors Associated With Nasolacrimal Duct Probing Failure Among Children in the Intelligent Research in Sight Registry

Importance

Understanding the factors associated with nasolacrimal duct probing failure in young children may help inform practice patterns.

Objective

To identify factors associated with repeated nasolacrimal duct probing in young children.

Design, Setting, and Participants

This retrospective cohort study analyzed data from the Intelligent Research in Sight (IRIS) Registry for all children who underwent nasolacrimal duct probing before 4 years of age between January 1, 2013, and December 31, 2020.

Main Outcomes and Measures

The Kaplan-Meier estimator was used to assess the cumulative incidence of a repeated procedure within 2 years of the initial procedure. Hazard ratios (HRs) derived from multivariable Cox proportional hazards regression models were used to evaluate the association between repeated probing and patient age, sex, race and ethnicity, geographic region, operative side, laterality of obstruction, type of initial procedure, and surgeon volume.

Results

This study included 19 357 children (9823 [50.7%] male; mean [SD] age, 1.40 [0.74] years) undergoing nasolacrimal duct probing. The cumulative incidence of repeated nasolacrimal duct probing was 7.2% (95% CI, 6.8%-7.5%) within 2 years of the initial procedure. Among 1333 repeated procedures, the second procedure involved silicone intubation in 669 (50.2%) and balloon catheter dilation in 256 (19.2%). Among 12 008 children aged 1 year or younger, office-based simple probing was associated with a slightly higher probability of reoperation compared with facility-based simple probing (9.5% [95% CI, 8.2%-10.8%] vs 7.1% [95% CI, 6.5%-7.7%]; P < .001). In the multivariable model, a greater risk of repeated probing was associated with bilateral obstruction (HR, 1.48; 95% CI, 1.32-1.65; P < .001) and office-based simple probing (HR, 1.33; 95% CI, 1.13-1.55; P < .001), and a lower risk was associated with primary balloon catheter dilation (HR, 0.69; 95% CI, 0.56-0.85; P < .001) and procedures performed by high-volume surgeons (HR, 0.84; 95% CI, 0.73-0.97; P = .02). Age, sex, race and ethnicity, geographic region, and operative side were not associated with reoperation risk in the multivariable model.

Conclusions and Relevance

In this cohort study, most children in the IRIS Registry undergoing nasolacrimal duct probing before 4 years of age did not require any additional intervention. Factors associated with lower risk of reoperation include surgeon experience, probing performed under anesthesia, and primary balloon catheter dilation.

Effectiveness of Microinvasive Glaucoma Surgery in the United States: Intelligent Research in Sight Registry Analysis 2013-2019

PURPOSE

To evaluate the effectiveness of microinvasive glaucoma surgery (MIGS) with and without concurrent phacoemulsification.

DESIGN

Multicenter, retrospective cohort study.

PARTICIPANTS

Patients in the Intelligent Research in Sight (IRIS®) Registry who underwent Xen gel stent (ab interno) implantation, endoscopic cyclophotocoagulation (ECP), or goniotomy or canaloplasty from 2013 through 2019.

METHODS

Kaplan-Meier survival analysis was used to assess reoperation rates. We defined reoperation as any subsequent glaucoma surgery occurring 1 month to 3 years after the initial procedure. Multivariable Cox proportional hazard models were used to determine factors predictive of reoperation.

MAIN OUTCOME MEASURES

Reoperation rate, mean intraocular pressure (IOP) and visual acuity (VA), postoperative complications, predictors of reoperation, and reoperation procedure type.

RESULTS

A total of 79 363 eyes from 57 561 patients were included, with 15 118 eyes (19%) receiving stand-alone MIGS and 64 245 eyes (81%) receiving MIGS concurrent with phacoemulsification. Overall, patients who underwent MIGS concurrently with phacoemulsification showed lower reoperation rates compared with stand-alone MIGS, most pronounced in ECP and goniotomy or canaloplasty. At postoperative year 2, the cumulative reoperation rate for stand-alone procedures was 15% for ECP, 24% for Xen implantation, and 24% for goniotomy or canaloplasty compared with 3% for ECP, 19% for Xen implantation, and 6% for goniotomy or canaloplasty concurrent with phacoemulsification (P < 0.001 for each stand-alone MIGS vs. MIGS with phacoemulsification). Black race, older age, moderate and severe glaucoma, higher baseline IOP, and glaucoma subtype were associated with higher reoperation risk. Although IOP decreased in all groups, stand-alone MIGS showed a more substantial decrease in mean IOP. Complication rates from MIGS were low overall: 1% for ECP, 1% for Xen implantation, and 2% for goniotomy or canaloplasty.

CONCLUSIONS

In current United States clinical practice, MIGS has substantially lower reoperation rates when performed with phacoemulsification, especially for ECP and goniotomy or canaloplasty. Approximately one-sixth of patients undergoing stand-alone ECP and one-quarter of patients undergoing stand-alone Xen implantation or goniotomy or canaloplasty require reoperation by 2 years. Black race, diagnosis coding of moderate to severe glaucoma, and higher baseline IOP were associated with higher risk of reoperation after MIGS procedures.

FINANCIAL DISCLOSURE(S)

Proprietary or commercial disclosure may be found after the references.

PyVisualFields: A Python Package for Visual Field Analysis

Purpose

Artificial intelligence (AI) methods are changing all areas of research and have a variety of capabilities of analysis in ophthalmology, specifically in visual fields (VFs) to detect or predict vision loss progression. Whereas most of the AI algorithms are implemented in Python language, which offers numerous open-source functions and algorithms, the majority of algorithms in VF analysis are offered in the R language. This paper introduces PyVisualFields, a developed package to address this gap and make available VF analysis in the Python language.

Methods

For the first version, the R libraries for VF analysis provided by vfprogression and visualFields packages are analyzed to define the overlaps and distinct functions. Then, we defined and translated this functionality into Python with the help of the wrapper library rpy2. Besides maintaining, the subsequent versions' milestones are established, and the third version will be R-independent.

Results

The developed Python package is available as open-source software via the GitHub repository and is ready to be installed from PyPI. Several Jupyter notebooks are prepared to demonstrate and describe the capabilities of the PyVisualFields package in the categories of data presentation, normalization and deviation analysis, plotting, scoring, and progression analysis.

Conclusions

We developed a Python package and demonstrated its functionality for VF analysis and facilitating ophthalmic research in VF statistical analysis, illustration, and progression prediction.

Translational Relevance

Using this software package, researchers working on VF analysis can more quickly create algorithms for clinical applications using cutting-edge AI techniques.

Effectiveness of Trabeculectomy and Tube Shunt with versus without Concurrent Phacoemulsification: Intelligent Research in Sight Registry Longitudinal Analysis

OBJECTIVE

To determine the effectiveness of trabeculectomy and glaucoma drainage device (GDD) surgery performed with concurrent phacoemulsification compared with stand-alone procedures.

DESIGN

Multicenter retrospective cohort study.

PARTICIPANTS

Patients in the Intelligent Research in Sight Registry who underwent trabeculectomy or GDD from 2013 through 2019.

METHODS

The Kaplan-Meier survival analysis was used to determine reoperation rates. Reoperation was defined as any subsequent glaucoma surgery occurring 1 month to 3 years after the initial procedure. Multivariable Cox proportional hazard models were used to determine reoperation risk factors.

MAIN OUTCOME MEASURES

Reoperation rate, intraocular pressure (IOP), visual acuity, reoperation procedure type, postoperative complications, and predictors of surgical failure.

RESULTS

A total of 117 697 eyes undergoing glaucoma surgery alone and 35 657 eyes undergoing surgery with phacoemulsification were included. The cumulative reoperation rates at postoperative years 1 and 3 were 4.9% and 11.5%, respectively, for trabeculectomy alone and 3.0% and 7.3%, respectively, for trabeculectomy combined with phacoemulsification (P < 0.001). The reoperation rates at postoperative 1 and 3 years were 3.8% and 7.8%, respectively, for GDD alone and 2.1% and 5.4%, respectively, for GDD with phacoemulsification (P < 0.001). Stand-alone procedures achieved greater IOP reduction by percentage change from baseline (trabeculectomy alone, 35.3% vs. trabeculectomy with phacoemulsification, 23.1%, P < 0.001; and GDD alone, 36.0% vs. GDD with phacoemulsification, 29.3%; P < 0.001). Visual acuity improved by 0.12 logarithm of the minimum angle of resolution (logMAR) (95% confidence interval [CI], 0.11-0.12) and 0.10 logMAR (95% CI, 0.08-0.11) after trabeculectomy and GDD with phacoemulsification and declined by 0.15 logMAR (95% CI, 0.14-0.15) and 0.12 logMAR (95% CI, 0.11-0.12) after stand-alone trabeculectomy and GDD. The overall documented complication rate was 2.9% for GDD and 1.4% for trabeculectomy. Age, sex, race, ethnicity, baseline IOP, and glaucoma diagnosis and severity were associated with surgical failure risk. The most common reoperation procedure was GDD.

CONCLUSIONS

Reoperation rates within the first 3 years after trabeculectomy and GDD with and without phacoemulsification were low. Trabeculectomy and GDD with phacoemulsification had lower reoperation rates than those with stand-alone procedures. However, stand-alone procedures resulted in greater IOP reduction compared with combined procedures. Postoperative complications were uncommon overall. Patient age, sex, race, ethnicity, baseline IOP, and glaucoma diagnosis and severity were associated with surgical success.

Wide-field swept-source optical coherence tomography angiography in the assessment of retinal microvasculature and choroidal thickness in patients with myopia

BACKGROUND/AIMS

Pathological myopia (PM) is a leading cause of blindness worldwide. We aimed to evaluate microvascular and chorioretinal changes in different stages of myopia with wide-field (WF) swept-source (SS) optical coherence tomography angiography (OCTA).

METHODS

This prospective cross-sectional observational study included 186 eyes of 122 patients who had undergone imaging between November 2018 and October 2020. Vessel density (VD) and vessel skeletonised density (VSD) of superficial capillary plexus, deep capillary plexus and whole retina, as well as foveal avascular zone parameters, retinal thickness (RT) and choroidal thickness (CT), were calculated.

RESULTS

This study evaluated 75 eyes of 48 patients with high myopia (HM), 43 eyes of 31 patients with mild to moderate myopia and 68 eyes of 53 age-matched controls. Controlling for age and the presence of systemic hypertension, we found that HM was associated with decrease in VD and VSD in all layers on 12×12 mm² scans. Furthermore, HM was associated with a VD and VSD decrease in every Early Treatment Diabetic Retinopathy Study grid, with a larger decrease temporally (β_VD=-0.39, β_VSD=-10.25, p<0.01). HM was associated with decreased RT and CT. Reduction in RT was outside the macular region, while reduction in CT was in the macular region.

CONCLUSION

Using WF SS-OCTA, we identified reduction in microvasculature and structural changes associated with myopia. Decrease in VD and VSD was greater in the temporal quadrant, and reductions in RT and CT were uneven across the retina. Further work may help identify risk factors for the progression of PM and associated vision-threatening complications.

Neurotrophic Keratopathy in the United States: An Intelligent Research in Sight Registry Analysis

PURPOSE

To describe the characteristics of neurotrophic keratopathy (NK) in the United States.

DESIGN

Retrospective database study.

PARTICIPANTS

Thirty-one thousand nine hundred fifteen eyes of 27 483 patients with a diagnosis of NK.

METHODS

Retrospective analysis of visits associated with a diagnosis of NK between 2013 and 2018 using the American Academy of Ophthalmology Intelligent Research in Sight (IRIS®) Registry.

MAIN OUTCOME MEASURES

Demographic information, prevalence, visual acuity (VA), concomitant diagnosis and procedure codes, and risk factors impacting VA most closely after NK onset date.

RESULTS

Mean ± standard deviation (SD) age at initial diagnosis of NK was 68.0 ± 16.0 years, and 58.91% of patients were women (P < 0.0001). Presentation was unilateral in 58.14%, bilateral in 16.13%, and unspecified in 25.73%. Average 6-year prevalence of NK in the IRIS Registry was 21.34 cases per 100 000 patients. Mean ± SD VA was 0.60 ± 0.79 logMAR before diagnosis and 0.88 ± 0.94 logMAR after diagnosis (P < 0.0001). Most common concomitant diagnoses included herpetic keratitis (33.70%), diabetes (31.59%), and corneal dystrophy (14.28%). Common procedures for NK management included the use of amniotic membrane (29.90%), punctal plugs (29.65%), and bandage contact lenses (22.67%). Age, male sex, Black race, Hispanic or Latino ethnicity, unilateral involvement, concomitant diagnoses of diabetes, corneal transplantation, and herpetic keratitis were associated significantly with worse VA.

CONCLUSIONS

Based on the IRIS Registry, the prevalence of NK is 21.34 cases per 100 000 patients. Visual acuity was significantly worse after NK diagnosis compared with other time points. Neurotrophic keratopathy was associated most commonly with herpetic keratitis and diabetes. Worse VA in patients with NK was associated with several demographic characteristics, history of diabetes, corneal transplantation, and herpetic keratitis.

Archetypal analysis of visual fields in optic neuritis reveals functional biomarkers associated with outcome and treatment response

BACKGROUND AND OBJECTIVES

Archetypal analysis (AA), a form of unsupervised machine learning, can identify quantifiable visual field (VF) patterns seen in optic neuritis (ON), known as archetypes (ATs). We hypothesized that AT weight changes over time would reflect the course of recovery and the effects of therapy in ON. We explored whether baseline AT weights would be associated with VF status at the clinical trial outcome and if ATs would indicate residual VF defects in eyes with mean deviation (MD) ≥ -2.00 at six months.

METHODS

We used a published 16-AT model derived from 3892 Optic Neuritis Treatment Trial VFs (456 eyes) for all analyses. We measured AT weight changes over the six-month study period and used asymptotic regression to analyze the rate of change. We compared AT weights at six months between treatment groups. We evaluated associations between baseline AT weight thresholds and VF outcome or treatment effect. We calculated residual AT weights in eyes with MD ≥ -2.00 dB at six months.

RESULTS

Over six months, AT1 (a normal VF pattern) demonstrated the greatest median weight change, increasing from 0.00% (IQR 0.00-0.00%) at baseline to 60.0% (IQR 38.3-70.8%) at six months (p < 0.001). At outcome, the intravenous methylprednisolone (IVMP) group had the highest median AT1 weight (IVMP: 63.3%, IQR 51.3-72.8%; placebo: 56.2%, IQR 35.1-71.6%; prednisone 58.3%, IQR 35.1-71.6%; p = 0.019). Eyes with AT1 weight ≥ 19% at baseline had superior median MD values (-0.91 vs. -2.07 dB, p < 0.001) and AT1 weights (70.8% vs. 57.8% p < 0.001) at six months. Only eyes with AT1 weight < 19% at baseline showed a treatment benefit for IVMP, with a higher six-month median AT1 weight compared to placebo (p = 0.015) and prednisone (p = 0.016), and a higher median MD compared to placebo (p = 0.027). At six months, 182 (80.2%) VFs with MD ≥ -2.00 had at least one abnormal AT.

DISCUSSION

Changes in quantifiable, archetypal patterns of VF loss reflect recovery in ON. Machine learning analysis of the VFs in optic neuritis reveals associations with response to therapy and VF outcome, and uncovers residual deficits, not readily seen with standard evaluations.

Risk Factors Associated With Pterygium Reoperation in the IRIS Registry

This cohort study used data from the IRIS Registry to assess the risk factors associated with pterygium reoperation.

Longitudinal visual field archetypal analysis of optic neuritis treated in a clinical setting

Background/aims

We previously used archetypal analysis (AA) to create a model that quantified patterns (archetypes (ATs)) of visual field (VF) loss that can predict recovery and reveal residual VF deficits from eyes in the Optic Neuritis Treatment Trial (ONTT). We hypothesised that AA could produce similar results for ON VFs collected in clinical practice.

Methods

We applied AA to a retrospective data set of 486 VFs collected in 1 neuro-ophthalmology service from 141 eyes with acute ON and typical VF defects, to create a clinic-derived AT model. We also used the ONTT-derived AT model to analyse this new dataset. We compared the findings of both models by decomposing VFs into component ATs of varying per cent weight (PW), correlating presentation AT PW with mean deviation (MD) at final visits for each eye and identifying residual deficits in VFs considered normal.

Results

Both models, each with 16 ATs, decomposed each presentation VF into 0–6 abnormal ATs representative of known patterns of ON-related VF loss. AT1, the normal pattern in both models, correlated strongly with MD for VFs collected at presentation (r=0.82; p<0.001) and the final visit (r=0.81, p<0.001). The presentation AT1 PW was associated with improvement in MD over time. 67% of VFs considered ‘normal’ at final visit had 1.2±0.4 abnormal ATs, and both models revealed similar patterns of regional VF loss.

Conclusions

AA is a quantitative method to measure change and outcome of ON VFs. Presentation AT features are associated with MD at final visit. AA identifies residual VF deficits not otherwise indicated by MD.