Effectiveness of Angle-Based Minimally Invasive Glaucoma Surgery after Laser Trabeculoplasty: An Analysis of the IRIS® Registry (Intelligent Research in Sight).

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 and without preceding laser trabeculoplasty (LTP). DESIGN: Retrospective cohort study. SUBJECTS: Eyes undergoing ab-MIGS (Canaloplasty, Goniotomy, Trabectome, and iStent) with and without prior laser trabeculoplasty (< 2 years preceding MIGS) were identified in the IRIS® Registry (Intelligent Research in Sight) 2013 to 2018. METHODS: Propensity score matching (PSM) was undertaken to define the following 4 cohorts: (1) standalone ab-MIGS, no prior LTP vs. (2) standalone ab-MIGS, with prior LTP; and (3) ab-MIGS + phacoemulsification, no prior LTP vs. (4) ab-MIGS + phacoemulsification, with prior LTP. MAIN OUTCOME MEASURES: Failure was defined as subsequent glaucoma reoperation after ab-MIGS (either MIGS or traditional glaucoma surgery). Time-to-event outcome and incidence rates were calculated using survival analysis, and adjusted hazard ratios (aHRs) were generated using multivariate Cox proportional hazards models. Medication data were not available for analysis. RESULTS: A total of 164 965 unique MIGS procedures were performed, from 2013 to 2018. After PSM, we identified 954 eyes undergoing standalone ab-MIGS and 7522 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 = 3761) were more likely to undergo reoperation vs. those without LTP (n = 3761) 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 for understanding who may benefit most from ab-MIGS, and for guiding patient and surgeon treatment expectations. FINANCIAL DISCLOSURES: Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Racial Disparities in Glaucoma Vision Outcomes and Eye Care Utilization: An IRIS Registry Analysis.

PURPOSE: To evaluate racial disparities in vision outcomes and eye care utilization among glaucoma patients. DESIGN: Retrospective cohort study. METHODS: In this population-based IRIS Registry (Intelligent Research in Sight) study, we included patients with minimum one diagnosis code for glaucoma at least 6 months prior to January 1, 2015 and at least one eye exam, visual field (VF), optical coherence tomography (OCT), or eye-related inpatient or emergency department (ED) code in 2015. Multivariable logistic and negative binomial regression models were used to assess vision and utilization outcomes, respectively, across race and ethnicity from January 1, 2015 to January 1, 2020. Vision outcomes included cup-to-disc ratio (CDR) progression > 0.80, poor vision (visual acuity 20/200 or worse), low vision codes, and need for glaucoma filtering surgery. Utilization outcomes included outpatient eye exams, OCTs, VFs, inpatient/ED encounters, and lasers/surgeries. RESULTS: Among 996,297 patients, 73% were non-Hispanic White, 15% non-Hispanic Black, 9% Hispanic, 3% Asian/Pacific Islander, and 0.3% Native American/Alaska Native. Compared to White eyes, Black and Hispanic eyes had higher adjusted odds of CDR progression (odds ratio [OR] = 1.12, 95% confidence interval [CI] = 1.08-1.17; OR = 1.28, 95% CI = 1.22-1.34), poor vision (OR = 1.26, 95% CI = 1.22-1.29; OR = 1.26, 95% CI = 1.22-1.31), glaucoma filtering surgery (rate ratio (RR) = 1.47, 95% CI = 1.42-1.51; RR = 1.13, 95% CI = 1.09-1.18). Hispanic eyes also had increased odds of low vision diagnoses (Hispanic OR = 1.18, 95% CI = 1.07-1.30). Black and Hispanic patients were less likely to have eye exams (RR = 0.94, 95% CI = 0.94-0.95; RR = 0.99, 95% CI = 0.99-0.99) and OCTs (RR = 0.86, 95% CI = 0.85-0.86; RR = 0.97, 95% CI = 0.96-0.98), yet Black patients had higher odds of inpatient/ED encounters (RR = 1.64, 95% CI = 1.37-1.96) compared to White patients. Native American patients were more likely to have poor vision (OR = 1.17, 95% CI = 1.01-1.36) and less likely to have outpatient visits (RR = 0.89, 95% CI = 0.86-0.91), OCTs (RR = 0.85, 95% CI = 0.82-0.89), visual fields (RR = 0.91, 95% CI = 0.88-0.94) or lasers/surgeries (RR = 0.87, 95% CI = 0.79-0.96) compared to White patients. CONCLUSIONS: We found that significant disparities in US eye care exist with Black, Hispanic, and Native American patients having worse vision outcomes and less disease monitoring. Glaucoma may be undertreated in these racial and ethnic minority groups, increasing risk for glaucoma-related vision loss.

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.

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/.

Anti-inflammatory medication use after cataract surgery: online survey of practice patterns.

PURPOSE: To determine current prescribing patterns for topical or intraocular/periocular anti-inflammatory medications (AIMs) after routine cataract surgery. SETTING: kera-net online members. DESIGN: Cross-sectional survey. METHODS: An online survey was distributed to subscribers of kera-net, a global online platform sponsored by the Cornea Society. Questions were asked regarding the use of topical or intraocular/periocular AIM after cataract surgery and types of medications prescribed. RESULTS: Of 217 surgeon respondents (23% response rate), 171 (79%) practiced in the United States and 171 (79%) were cornea subspecialists. Most of the respondents (n = 196, 97%) prescribed topical corticosteroids after routine cataract surgery. The most frequently prescribed were prednisolone acetate (n = 162, 83%), followed by dexamethasone (n = 26, 13%), difluprednate (n = 24, 12%), and loteprednol etabonate (n = 13, 7%). Corticosteroids comprised (n = 40, 32%) of total intraocular/periocular injections, with triamcinolone acetonide 10 or 40 mg (n = 19, 47.5%) most commonly used. 23 surgeons (58%) who utilized intraocular/periocular corticosteroids also prescribed topical corticosteroids. Topical nonsteroidal anti-inflammatory drugs were prescribed postoperatively by 148 surgeons (73%). CONCLUSIONS: Most surgeons prescribed topical AIM after routine cataract surgery. Many surgeons injected intraocular or periocular AIM while prescribing topical AIM. The diversity of practice patterns may reflect the lack of clear evidence-based guidelines.

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.

Trends in Anti-Vascular Endothelial Growth Factor Original Medicare Part B Claims in the United States, 2014-2019.

PURPOSE: To characterize trends in use of and expenditure for the intravitreal injection of anti-vascular endothelial growth factor (anti-VEGF) agents aflibercept, ranibizumab, and bevacizumab among the population enrolled in Original Medicare from 2014 to 2019. METHODS: The Centers for Medicare and Medicaid Services Physician and Other Supplier Public Use File was used to extract Medicare Part B fee-for-service outpatient injection claims data submitted by ophthalmologists. Multivariable linear regression models were used to evaluate the association between reimbursement, ophthalmologist availability, and agent administration rate. RESULTS: Between 2014 and 2019, 17,588,995 intravitreal injection claims were filed by 4218 US ophthalmologists. Medicare costs for anti-VEGF injections increased from 2.51 B USD in 2014 to 4.02 B USD in 2019. Increased state-level ophthalmologist availability and incremental increases in average reimbursement amounts were found to be significantly associated with a 6.8-fold variation in 2019 overall anti-VEGF injection rates across states. CONCLUSIONS: Medicare injection rates and costs for anti-VEGF injections have both increased between 2014 and 2019, largely driven by increased aflibercept use. There is a significant association between ophthalmologist availability and anti-VEGF injection rate on the state level, suggesting access to care may contribute to the observed state-level disparities in intravitreal injection rates. Further characterization of factors contributing to the state-level variation in injection rates of individual anti-VEGF agents may help inform interventions promoting equitable access to and use of these drugs.

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.

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.

Wound modulation in glaucoma surgery: The role of anti-scarring agents.

Filtration surgery is one of the most frequently performed surgeries in the management of glaucoma, and trabeculectomy is considered the gold standard surgical technique for the same. Though trabeculectomy has been reported to have an excellent initial success rate, about 30% of them fail in 3 years, and nearly 50% of them fail in 5 years. The most significant risk of failure still seems to be wound scarring, especially episcleral fibrosis, leading to bleb failure. As a result, it is essential to explore the role of anti-scarring agents, including mitomycin C, and 5-fluorouracil in wound modulation and improving the bleb survival rate. Since these agents are widely used in trabeculectomy, it is crucial to understand the various modes of application, advantages, and adverse effects of these agents. On an evidence-based approach, all these points have been highlighted in this review article. In addition, the newer agents available for wound modulation and their scope for practical application are discussed.

TIE1 and TEK signalling, intraocular pressure, and primary open-angle glaucoma: a Mendelian randomization study.

BACKGROUND: In primary open-angle glaucoma (POAG), lowering intraocular pressure (IOP) is the only proven way of slowing vision loss. Schlemm's canal (SC) is a hybrid vascular and lymphatic vessel that mediates aqueous humour drainage from the anterior ocular chamber. Animal studies support the importance of SC endothelial angiopoietin-TEK signalling, and more recently TIE1 signalling, in maintaining normal IOP. However, human genetic support for a causal role of TIE1 and TEK signalling in lowering IOP is currently lacking. METHODS: GWAS summary statistics were obtained for plasma soluble TIE1 (sTIE1) protein levels (N = 35,559), soluble TEK (sTEK) protein levels (N = 35,559), IOP (N = 139,555) and POAG (Ncases = 16,677, Ncontrols = 199,580). Mendelian randomization (MR) was performed to estimate the association of genetically proxied TIE1 and TEK protein levels with IOP and POAG liability. Where significant MR estimates were obtained, genetic colocalization was performed to assess the probability of a shared causal variant (PPshared) versus distinct (PPdistinct) causal variants underlying TIE1/TEK signalling and the outcome. Publicly available single-nucleus RNA-sequencing data were leveraged to investigate differential expression of TIE1 and TEK in the human ocular anterior segment. RESULTS: Increased genetically proxied TIE1 signalling and TEK signalling associated with a reduction in IOP (- 0.21 mmHg per SD increase in sTIE1, 95% CI = - 0.09 to - 0.33 mmHg, P = 6.57 × 10-4, and - 0.14 mmHg per SD decrease in sTEK, 95% CI = - 0.03 to - 0.25 mmHg, P = 0.011), but not with POAG liability. Colocalization analysis found that the probability of a shared causal variant was greater for TIE1 and IOP than for TEK and IOP (PPshared/(PPdistinct + PPshared) = 0.98 for TIE1 and 0.30 for TEK). In the anterior segment, TIE1 and TEK were preferentially expressed in SC, lymphatic, and vascular endothelium. CONCLUSIONS: This study provides novel human genetic support for a causal role of both TIE1 and TEK signalling in regulating IOP. Here, combined evidence from cis-MR and colocalization analyses provide stronger support for TIE1 than TEK as a potential IOP-lowering therapeutic target.

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.

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.

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.

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.

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.

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.

Relationship between Claims-Based Frailty Index and Eye Care Utilization among Medicare Beneficiaries with Glaucoma.

PURPOSE: To determine differences in eye care utilization by frailty levels among Medicare beneficiaries with glaucoma. DESIGN: Retrospective cohort study. PARTICIPANTS: Medicare fee-for-service beneficiaries over 65 years of age with glaucoma, identified using International Classification of Diseases codes before July 1, 2014. METHODS: By using a validated claims-based frailty index (range, 0-1), beneficiaries were classified as nonfrail/prefrail (0-0.19), mildly frail (0.20-0.29), and moderate-to-severely frail (≥ 0.30). Negative binomial regression analyses were used to estimate incident rate ratios (IRRs) of eye care utilization by frailty levels between July 1, 2014, and December 31, 2016. MAIN OUTCOME MEASURES: Current Procedural Terminology codes for eye examinations and eye care-related office visits; eye care-related inpatient and emergency department (ED) encounters; eye care-related nursing facility and home-visit encounters; visual field (VF) and retinal nerve fiber layer (RNFL) OCT tests; and selective laser trabeculoplasties (SLTs) and glaucoma surgeries. RESULTS: Among 76 260 Medicare beneficiaries with glaucoma, the mean age was 78.9 years (standard deviation, 7.8), female beneficiaries constituted 60.5%, and 78.7% of beneficiaries self-identified as non-Hispanic White. According to a claims-based frailty index, 79.5% of beneficiaries were nonfrail/prefrail, 17.1% were mildly frail, and 3.4% were moderate-to-severely frail. Moderate-to-severely frail beneficiaries were less likely than nonfrail/prefrail beneficiaries to have outpatient encounters (IRR, 0.85, 95% confidence interval [CI], 0.83-0.88); VF tests (IRR, 0.64, 95% CI, 0.60-0.67); RNFL OCT tests (IRR, 0.77, 95% CI, 0.73-0.81); SLT (IRR, 0.74, 95% CI, 0.60-0.92); and glaucoma surgery (IRR, 0.74, 95% CI 0.55-0.99), after adjusting for age, gender, glaucoma severity, race, and socioeconomic status. Compared with nonfrail/prefrail beneficiaries, moderate-to-severely frail beneficiaries had higher rates of inpatient/ED encounters (IRR, 5.03, 95% CI, 2.36-10.71) and nursing facility/home-visit encounters (IRR, 34.89, 95% CI, 14.82-82.13). CONCLUSIONS: Compared with nonfrail/prefrail Medicare beneficiaries with glaucoma, beneficiaries with moderate-to-severe frailty had lower rates of eye care utilization in the outpatient setting and higher rates of utilization in acute care settings. This suggests that frail patients may receive less disease monitoring and fewer interventions for their glaucoma management. 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.