Faculty Bibliography

PURPOSE:To determine whether peripapillary atrophy (PPA) area is an indicator of glaucomatous structural and functional damage and progression. METHODS:In this retrospective longitudinal analysis from ongoing prospective study we qualified 71 eyes (50 subjects) with glaucoma. All subjects had a comprehensive ophthalmic examination, visual field (VF), and spectral-domain optical coherence tomography (OCT) testing in at least three visits. PPA was manually delineated on en face OCT optic nerve head scans, while observing the corresponding cross-sectional images, as the hyper-reflective area contiguous with the optic disc. RESULTS:The mean follow-up duration was 4.4 ± 1.4 years with an average of 6.8 ± 2.2 visits. At baseline, PPA area was significantly associated only with VF's mean deviation (MD; P = 0.041), visual field index (VFI; P = 0.041), superior ganglion cell inner plexiform layer (GCIPL; P = 0.011), and disc area (P = 0.011). Longitudinally, PPA area was negatively and significantly associated with MD (P = 0.015), VFI (P = 0.035), GCIPL (P = 0.009), superior GCIPL (P = 0.034), and disc area (P = 0.007, positive association). CONCLUSIONS:Longitudinal change in PPA area is an indicator of glaucomatous structural and functional progression but PPA area at baseline cannot predict future progression. TRANSLATIONAL RELEVANCE:Longitudinal changes in peripapillary atrophy area measured by OCT can be an indicator of structural and functional glaucoma progression.

PURPOSE/UNASSIGNED:To establish generalizable pointwise spatial relationship between structure and function through occlusion analysis of a deep-learning (DL) model for predicting the visual field (VF) sensitivities from 3-dimensional (3D) OCT scan. DESIGN/UNASSIGNED:Retrospective cross-sectional study. PARTICIPANTS/UNASSIGNED:A total of 2151 eyes from 1129 patients. METHODS/UNASSIGNED:A DL model was trained to predict 52 VF sensitivities of 24-2 standard automated perimetry from 3D spectral-domain OCT images of the optic nerve head (ONH) with 12 915 OCT-VF pairs. Using occlusion analysis, the contribution of each individual cube covering a 240 × 240 × 31.25 μm region of the ONH to the model's prediction was systematically evaluated for each OCT-VF pair in a separate test set that consisted of 996 OCT-VF pairs. After simple translation (shifting in x- and y-axes to match the ONH center), group t-statistic maps were derived to visualize statistically significant ONH regions for each VF test point within a group. This analysis allowed for understanding the importance of each super voxel (240 × 240 × 31.25 μm covering the entire 4.32 × 4.32 × 1.125 mm ONH cube) in predicting VF test points for specific patient groups. MAIN OUTCOME MEASURES/UNASSIGNED:The region at the ONH corresponding to each VF test point and the effect of the former on the latter. RESULTS/UNASSIGNED:The test set was divided to 2 groups, the healthy-to-early-glaucoma group (792 OCT-VF pairs, VF mean deviation [MD]: -1.32 ± 1.90 decibels [dB]) and the moderate-to-advanced-glaucoma group (204 OCT-VF pairs, VF MD: -17.93 ± 7.68 dB). Two-dimensional group t-statistic maps (x, y projection) were generated for both groups, assigning related ONH regions to visual field test points. The identified influential structural locations for VF sensitivity prediction at each test point aligned well with existing knowledge and understanding of structure-function spatial relationships. CONCLUSIONS/UNASSIGNED:This study successfully visualized the global trend of point-by-point spatial relationships between OCT-based structure and VF-based function without the need for prior knowledge or segmentation of OCTs. The revealed spatial correlations were consistent with previously published mappings. This presents possibilities of learning from trained machine learning models without applying any prior knowledge, potentially robust, and free from bias. FINANCIAL DISCLOSURES/UNASSIGNED:Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

PURPOSE:The presence of imbalanced datasets in medical applications can negatively affect deep learning methods. This study aims to investigate how the performance of convolutional neural networks (CNNs) for glaucoma diagnosis can be improved by addressing imbalanced learning issues through utilizing glaucoma suspect samples, which are often excluded from studies because they are a mixture of healthy and preperimetric glaucomatous eyes, in a semi-supervised learning approach. METHODS:A baseline 3D CNN was developed and trained on a real-world glaucoma dataset, which is naturally imbalanced (like many other real-world medical datasets). Then, three methods, including reweighting samples, data resampling to form balanced batches, and semi-supervised learning on glaucoma suspect data were applied to practically assess their impacts on the performances of the trained methods. RESULTS:The proposed method achieved a mean accuracy of 95.24%, an F1 score of 97.42%, and an area under the curve of receiver operating characteristic (AUC ROC) of 95.64%, whereas the corresponding results for the traditional supervised training using weighted cross-entropy loss were 92.88%, 96.12%, and 92.72%, respectively. The obtained results show statistically significant improvements in all metrics. CONCLUSIONS:Exploiting glaucoma suspect eyes in a semi-supervised learning method coupled with resampling can improve glaucoma diagnosis performance by mitigating imbalanced learning issues. TRANSLATIONAL RELEVANCE:Clinical imbalanced datasets may negatively affect medical applications of deep learning. Utilizing data with uncertain diagnosis, such as glaucoma suspects, through a combination of semi-supervised learning and class-imbalanced learning strategies can partially address the problems of having limited data and learning on imbalanced datasets.

PURPOSE/UNASSIGNED:The structural changes measured by optical coherence tomography (OCT) are related to functional changes in visual fields (VFs). This study aims to accurately assess the structure-function relationship and overcome the challenges brought by the minimal measurable level (floor effect) of segmentation-dependent OCT measurements commonly used in prior studies. METHODS/UNASSIGNED:We developed a deep learning model to estimate the functional performance directly from three-dimensional (3D) OCT volumes and compared it to the model trained with segmentation-dependent two-dimensional (2D) OCT thickness maps. Moreover, we proposed a gradient loss to utilize the spatial information of VFs. RESULTS/UNASSIGNED:Our 3D model was significantly better than the 2D model both globally and pointwise regarding both mean absolute error (MAE = 3.11 + 3.54 vs. 3.47 ± 3.75 dB, P < 0.001) and Pearson's correlation coefficient (0.80 vs. 0.75, P < 0.001). On a subset of test data with floor effects, the 3D model showed less influence from floor effects than the 2D model (MAE = 5.24 ± 3.99 vs. 6.34 ± 4.58 dB, P < 0.001, and correlation 0.83 vs. 0.74, P < 0.001). The gradient loss improved the estimation error for low-sensitivity values. Furthermore, our 3D model outperformed all prior studies. CONCLUSIONS/UNASSIGNED:By providing a better quantitative model to encapsulate the structure-function relationship more accurately, our method may help deriving VF test surrogates. TRANSLATIONAL RELEVANCE/UNASSIGNED:DL-based VF surrogates not only benefit patients by reducing the testing time of VFs but also allow clinicians to make clinical judgments without the inherent limitations of VFs.

PURPOSE/UNASSIGNED:Lamina cribrosa (LC) deformation is hypothesized to play a major role in glaucoma pathogenesis. The purpose of this study was to determine in vivo how varying intraocular pressure (IOP) under fixed intracranial pressure (ICP), and vice versa, deforms the pore paths throughout the LC volume. METHODS/UNASSIGNED:Spectral-domain optical coherence tomography scans of the optic nerve head were acquired from healthy adult rhesus monkeys under different pressures. IOP and ICP were controlled with gravity-based perfusion systems into the anterior chamber and lateral ventricle, respectively. IOP and ICP were modulated from baseline to high (19-30 mmHg) and highest (35-50 mmHg) levels while maintaining a fixed ICP of 8 to 12 mmHg and IOP of 15 mmHg, respectively. After three-dimensional registration and segmentation, the paths of pores visible in all settings were tracked based on their geometric centroids. Pore path tortuosity was defined as the measured distance divided by the minimal distance between the most anterior and posterior centroids. RESULTS/UNASSIGNED:The median pore tortuosity at baseline varied among the eyes (range, 1.16-1.68). For the IOP effect under fixed ICP (six eyes, five animals), two eyes showed statistically significant increased tortuosity and one showed a decrease (P < 0.05, mixed-effects model). No significant change was detected in three eyes. When modulating ICP under fixed IOP (five eyes, four animals), a similar response pattern was detected. CONCLUSIONS/UNASSIGNED:Baseline pore tortuosity and the response to acute pressure increase vary substantially across eyes. TRANSLATIONAL RELEVANCE/UNASSIGNED:LC pore path tortuosity could be associated with glaucoma susceptibility.

BACKGROUND:Spectral-domain (SD-) optical coherence tomography (OCT) can reliably measure axonal (peripapillary retinal nerve fiber layer [pRNFL]) and neuronal (macular ganglion cell + inner plexiform layer [GCIPL]) thinning in the retina. Measurements from 2 commonly used SD-OCT devices are often pooled together in multiple sclerosis (MS) studies and clinical trials despite software and segmentation algorithm differences; however, individual pRNFL and GCIPL thickness measurements are not interchangeable between devices. In some circumstances, such as in the absence of a consistent OCT segmentation algorithm across platforms, a conversion equation to transform measurements between devices may be useful to facilitate pooling of data. The availability of normative data for SD-OCT measurements is limited by the lack of a large representative world-wide sample across various ages and ethnicities. Larger international studies that evaluate the effects of age, sex, and race/ethnicity on SD-OCT measurements in healthy control participants are needed to provide normative values that reflect these demographic subgroups to provide comparisons to MS retinal degeneration. METHODS:Participants were part of an 11-site collaboration within the International Multiple Sclerosis Visual System (IMSVISUAL) consortium. SD-OCT was performed by a trained technician for healthy control subjects using Spectralis or Cirrus SD-OCT devices. Peripapillary pRNFL and GCIPL thicknesses were measured on one or both devices. Automated segmentation protocols, in conjunction with manual inspection and correction of lines delineating retinal layers, were used. A conversion equation was developed using structural equation modeling, accounting for clustering, with healthy control data from one site where participants were scanned on both devices on the same day. Normative values were evaluated, with the entire cohort, for pRNFL and GCIPL thicknesses for each decade of age, by sex, and across racial groups using generalized estimating equation (GEE) models, accounting for clustering and adjusting for within-patient, intereye correlations. Change-point analyses were performed to determine at what age pRNFL and GCIPL thicknesses exhibit accelerated rates of decline. RESULTS:The healthy control cohort (n = 546) was 54% male and had a wide distribution of ages, ranging from 18 to 87 years, with a mean (SD) age of 39.3 (14.6) years. Based on 346 control participants at a single site, the conversion equation for pRNFL was Cirrus = -5.0 + (1.0 × Spectralis global value). Based on 228 controls, the equation for GCIPL was Cirrus = -4.5 + (0.9 × Spectralis global value). Standard error was 0.02 for both equations. After the age of 40 years, there was a decline of -2.4 μm per decade in pRNFL thickness ( P < 0.001, GEE models adjusting for sex, race, and country) and -1.4 μm per decade in GCIPL thickness ( P < 0.001). There was a small difference in pRNFL thickness based on sex, with female participants having slightly higher thickness (2.6 μm, P = 0.003). There was no association between GCIPL thickness and sex. Likewise, there was no association between race/ethnicity and pRNFL or GCIPL thicknesses. CONCLUSIONS:A conversion factor may be required when using data that are derived between different SD-OCT platforms in clinical trials and observational studies; this is particularly true for smaller cross-sectional studies or when a consistent segmentation algorithm is not available. The above conversion equations can be used when pooling data from Spectralis and Cirrus SD-OCT devices for pRNFL and GCIPL thicknesses. A faster decline in retinal thickness may occur after the age of 40 years, even in the absence of significant differences across racial groups.