Faculty Bibliography

Excitotoxicity from the impairment of glutamate uptake constitutes an important mechanism in neurodegenerative diseases such as Alzheimer's, multiple sclerosis, and Parkinson's disease. Within the eye, excitotoxicity is thought to play a critical role in retinal ganglion cell death in glaucoma, diabetic retinopathy, retinal ischemia, and optic nerve injury, yet how excitotoxic injury impacts different retinal layers is not well understood. Here, we investigated the longitudinal effects of N-methyl-D-aspartate (NMDA)-induced excitotoxic retinal injury in a rat model using deep learning-assisted retinal layer thickness estimation. Before and after unilateral intravitreal NMDA injection in nine adult Long Evans rats, spectral-domain optical coherence tomography (OCT) was used to acquire volumetric retinal images in both eyes over 4 weeks. Ten retinal layers were automatically segmented from the OCT data using our deep learning-based algorithm. Retinal degeneration was evaluated using layer-specific retinal thickness changes at each time point (before, and at 3, 7, and 28 days after NMDA injection). Within the inner retina, our OCT results showed that retinal thinning occurred first in the inner plexiform layer at 3 days after NMDA injection, followed by the inner nuclear layer at 7 days post-injury. In contrast, the retinal nerve fiber layer exhibited an initial thickening 3 days after NMDA injection, followed by normalization and thinning up to 4 weeks post-injury. Our results demonstrated the pathological cascades of NMDA-induced neurotoxicity across different layers of the retina. The early inner plexiform layer thinning suggests early dendritic shrinkage, whereas the initial retinal nerve fiber layer thickening before subsequent normalization and thinning indicates early inflammation before axonal loss and cell death. These findings implicate the inner plexiform layer as an early imaging biomarker of excitotoxic retinal degeneration, whereas caution is warranted when interpreting the ganglion cell complex combining retinal nerve fiber layer, ganglion cell layer, and inner plexiform layer thicknesses in conventional OCT measures. Deep learning-assisted retinal layer segmentation and longitudinal OCT monitoring can help evaluate the different phases of retinal layer damage upon excitotoxicity.

The cellular-level visualization of retinal microstructures such as blood vessel wall components, not available with other imaging modalities, is provided with unprecedented details by dark-field imaging configurations; however, the interpretation of such images alone is sometimes difficult since multiple structural disturbances may be present in the same time. Particularly in eyes with retinal pathology, microstructures may appear in high-resolution retinal images with a wide range of sizes, sharpnesses, and brightnesses. In this paper we show that motion contrast and phase gradient imaging modalities, as well as the simultaneous acquisition of depth-resolved optical coherence tomography (OCT) images, provide additional insight to help understand the retinal neural and vascular structures seen in dark-field images and may enable improved diagnostic and treatment plans.

PURPOSE/UNASSIGNED:Broken stick analysis is a widely used approach for detecting unknown breakpoints where the association between measurements is nonlinear. We propose LIMBARE, an advanced linear mixed-effects breakpoint analysis with robust estimation, especially designed for longitudinal ophthalmic studies. LIMBARE accommodates repeated measurements from both eyes and over time, and it effectively addresses the presence of outliers. METHODS/UNASSIGNED:The model setup of LIMBARE and the computing algorithm for point and confidence interval estimates of the breakpoint were introduced. The performance of LIMBARE and other competing methods was assessed via comprehensive simulation studies and application to a longitudinal ophthalmic study with 216 eyes (145 subjects) followed for an average of 3.7 ± 1.3 years to examine the longitudinal association between structural and functional measurements. RESULTS/UNASSIGNED:In simulation studies, LIMBARE showed the smallest bias and mean squared error for estimating the breakpoint, with an empirical coverage probability of corresponding confidence interval estimates closest to the nominal level for scenarios with and without outlier data points. In the application to the longitudinal ophthalmic study, LIMBARE detected two breakpoints between visual field mean deviation (MD) and retinal nerve fiber layer thickness and one breakpoint between MD and cup-to-disc ratio, whereas the cross-sectional analysis approach detected only one and none, respectively. CONCLUSIONS/UNASSIGNED:LIMBARE enhances breakpoint estimation accuracy in longitudinal ophthalmic studies, and the cross-sectional analysis approach is not recommended for future studies. TRANSLATIONAL RELEVANCE/UNASSIGNED:Our proposed method and companion R package provide a valuable computational tool for advancing longitudinal ophthalmology research and exploring the association relationships among ophthalmic variables.

The cellular-level visualization of retinal microstructures such as blood vessel wall components, not available with other imaging modalities, is provided with unprecedented details by dark-field imaging configurations; however, the interpretation of such images alone is sometimes difficult since multiple structural disturbances may be present in the same time. Particularly in eyes with retinal pathology, microstructures may appear in high-resolution retinal images with a wide range of sizes, sharpnesses, and brightnesses. In this paper we show that motion contrast and phase gradient imaging modalities, as well as the simultaneous acquisition of depth-resolved optical coherence tomography (OCT) images, provide additional insight to help understand the retinal neural and vascular structures seen in dark-field images and may enable improved diagnostic and treatment plans.

PURPOSE/UNASSIGNED:Race disparities in the healthcare system and the resulting inequality in clinical data among different races hinder the ability to generate equitable prediction results. This study aims to reduce healthcare disparities arising from data imbalance by leveraging advanced transfer learning (TL) methods. METHOD/UNASSIGNED:We examined the ophthalmic healthcare disparities at a population level using electronic medical records data from a study cohort (N = 785) receiving care at an academic institute. Regression-based TL models were usesd, transferring valuable information from the dominant racial group (White) to improve visual field mean deviation (MD) rate of change prediction particularly for data-disadvantaged African American (AA) and Asian racial groups. Prediction results of TL models were compared with two conventional approaches. RESULTS/UNASSIGNED:Disparities in socioeconomic status and baseline disease severity were observed among the AA and Asian racial groups. The TL approach achieved marked to comparable improvement in prediction accuracy compared to the two conventional approaches as evident by smaller mean absolute errors or mean square errors. TL identified distinct key features of visual field MD rate of change for each racial group. CONCLUSIONS/UNASSIGNED:The study introduces a novel application of TL that improved reliability of the analysis in comparison with conventional methods, especially in small sample size groups. This can improve assessment of healthcare disparity and subsequent remedy approach. TRANSLATIONAL RELEVANCE/UNASSIGNED:TL offers an equitable and efficient approach to mitigate healthcare disparities analysis by enhancing prediction performance for data-disadvantaged group.

Glaucoma causes progressive visual field deterioration and is the leading cause of blindness worldwide. Glaucomatous damage is irreversible and greatly impacts quality of life. Therefore, it is critically important to detect glaucoma early and closely monitor progression to preserve functional vision. Glaucoma is routinely monitored in the clinical setting using optical coherence tomography (OCT) for derived measures such as the thickness of important visual structures. There is not a consensus of what measures represent the most relevant biomarkers of glaucoma progression. Further, despite the increasing availability of longitudinal OCT data, a quantitative model of 3D structural change over time associated with glaucoma does not exist. In this paper we present an algorithm that will perform hierarchical geodesic modeling at the imaging level, considering 3D OCT images as observations of structural change over time. Hierarchical modeling includes subject-wise trajectories as geodesics in the space of diffeomorphisms and population level (glaucoma vs control) trajectories are also geodesics which explain subject-wise trajectories as deviations from the mean. Our preliminary experiments demonstrate a greater magnitude of structural change associated with glaucoma compared to normal aging. Our algorithm has the potential application in patient-specific monitoring and analysis of glaucoma progression as well as a statistical model of population trends and population variability.

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.

Glaucoma is an age-related neurodegenerative disease of the visual system, affecting both the eye and the brain. Yet its underlying metabolic mechanisms and neurobehavioral relevance remain largely unclear. Here, using proton magnetic resonance spectroscopy and functional magnetic resonance imaging, we investigated the GABAergic and glutamatergic systems in the visual cortex of glaucoma patients, as well as neural specificity, which is shaped by GABA and glutamate signals and underlies efficient sensory and cognitive functions. Our study shows that among the older adults, both GABA and glutamate levels decrease with increasing glaucoma severity regardless of age. Further, our study shows that the reduction of GABA but not glutamate predicts the neural specificity. This association is independent of the impairments on the retina structure, age, and the gray matter volume of the visual cortex. Our results suggest that glaucoma-specific decline of GABA undermines neural specificity in the visual cortex and that targeting GABA could improve the neural specificity in glaucoma.

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:Rhesus macaques (Macaca mulatta) are the premier nonhuman primate model for studying human health and disease. We investigated if age was associated with clinically relevant ocular features in a large cohort of free-ranging rhesus macaques from Cayo Santiago, Puerto Rico. METHODS:We evaluated 120 rhesus macaques (73 males, 47 females) from 0 to 29 years old (mean ± SD: 12.6 ± 6.4) from September to December 2021. The ophthalmic evaluation included intraocular pressure (IOP) assessment, corneal pachymetry, biomicroscopy, A-scan biometry, automated refraction, and fundus photography after pupil dilation. The associations of age with the outcomes were investigated through multilevel mixed-effects models adjusted for sex and weight. RESULTS:On average, IOP, pachymetry, axial length, and automated refraction spherical equivalent were 18.37 ± 4.68 mmHg, 474.43 ± 32.21 µm, 19.49 ± 1.24 mm, and 0.30 ± 1.70 diopters (D), respectively. Age was significantly associated with pachymetry (β coefficient = -1.20; 95% confidence interval [CI], -2.27 to -0.14; P = 0.026), axial length (β coefficient = 0.03; 95% CI, 0.01 to 0.05; P = 0.002), and spherical equivalent (β coefficient = -0.12; 95% CI, -0.22 to -0.02; P = 0.015). No association was detected between age and IOP. The prevalence of cataracts in either eye was 10.83% (95% CI, 6.34-17.89) and was significantly associated with age (odds ratio [OR] = 1.20; 95% CI, 1.06-1.36; P = 0.004). Retinal drusen in either eye was observed in 15.00% (95% CI, 9.60-22.68) of animals, which was also significantly associated with age (OR = 1.14; 95% CI, 1.02-1.27; P = 0.020). CONCLUSIONS:Rhesus macaques exhibit age-related ocular associations similar to those observed in human aging, including decreased corneal thickness, increased axial length, myopic shift, and higher prevalence of cataract and retinal drusen.