Faculty Bibliography

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.

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.

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

This study was registered with ClinicalTrials.gov (ID: NCT03715231). A total of 20 participants (37 eyes) who were 18 or older and had glaucoma or were glaucoma suspects were enrolled from the NYU Langone Eye Center and Bellevue Hospital. During their usual ophthalmology visit, they were consented for the study and underwent 360-degree goniophotography using the NIDEK Gonioscope GS-1. Afterwards, the three ophthalmologists separately examined the images obtained and determined the status of the iridocorneal angle in four quadrants using the Shaffer grading system. Physicians were masked to patient names and diagnoses. Inter-observer reproducibility was determined using Fleiss' kappa statistics. The interobserver reliability using Fleiss' statistics was shown to be significant between three glaucoma specialists with fair overall agreement (Fleiss' kappa: 0.266, p < .0001) in the interpretation of 360-degree goniophotos. Automated 360-degree goniophotography using the NIDEK Gonioscope GS-1 have quality such that they are interpreted similarly by independent expert observers. This indicates that angle investigation may be performed using this automated device and that interpretation by expert observers is likely to be similar. Images produced from automated 360-degree goniophotography using the NIDEK Gonioscope GS-1 are similarly interpreted amongst glaucoma specialists, thus supporting use of this technique to document and assess the anterior chamber angle in patients with, or suspected of, glaucoma and iridocorneal angle abnormalities.

Purpose : Multiple sublayers of retina can be visualized with visible light (vis-) OCT.However, image quality can be compromised due to patient movement, cataracts, small pupil size, and light scattering causing haziness and variability in signal to noise ratio in individual A-scans and in entire B-scans.The purpose of this study was to examine the effect of conventional and deep neural network dehazing techniques on the visibility and quantitative assessment of retinal sub-layers on vis-OCT images. Methods : 9 healthy and 5 glaucoma subjects were scanned 3 times during one session.Scanning was done on the superior nasal side of para-foveal region,1.5 mm from the fovea with a 3D speckle reduction raster scanning protocol(3x3x1.6 mm with 8192x16x1024 samplings) using a prototype vis-OCT system.16 A-scan lines were averaged to reduce speckle noise.Gray-scale image dehazing guided by depth information and pretrained Dehazenet deep model following deep convolutional neural network with residual learning(DnCNN) were applied on original B-scans.Quality improvement were evaluated using quality index(QI) and contrast to noise ratio(CNR) on dehazed B-scans.For each subject, the dehazed B-scan of Dehazenet and DnCNN from a fixed location adjacent to the fovea were selected.The distances between each of 3 bright inner plexiform layers(IPL) and retinal pigment epithelium(RPE) sublayers were segmented manually for thickness measurements using a 8 A-scan averaged profile(Fig.).Coefficient of variations (CVs) were calculated to assess the measurement repeatability of the sublayers on original and dehazed B-scans. Results : Healthy and glaucoma subjects were age 45.67+/-11.7and 59.60+/-13.4(p=0.07,t-test),visual field mean deviation(MD)-1.55 to1.20 dB,and from -26.42 to -7.70dB(p= 0.003,Wilcoxon),global mean circumpapillary retinal nerve fiber layer(RNFL)thickness 96.33+/-12.20 and 59.80+/-9.09mm(p<0.001,Wilcoxon),respectively.Dehazed B-scans obtained by deep models have statistically significant better QI and CNR(Table1).Overall intra-subject CVs showed significantly improved reproducibility on all measured sub-layers of dehazed B-scans compared to original scans for all subjects(Tables 2,3). Conclusions : Vis-OCT image quality can be improved using deep neural network dehazing model resulting in higher reproducible thickness measurements of retinal sublayers within subjects in dehazed B-scans

Purpose : Automated segmentation of in-vivo lamina cribrosa (LC) has been challenging, owing to the complex 3D structure and decreased visibility in the lamina depth. Frangi's vesselness filter, which was originally developed for angiogram segmentation, have been successfully demonstrated in segmenting the ex-vivo LC from micro-CT and second harmonic generation microscopy images. In this project we are proposing a new approach of segmenting the in vivo LC from OCT scans, incorporating the Frangi's vesselness principle to facilitate in vivo LC image analysis in much greater detail compared to our previously described 3D analysis method. Methods : In-vivo spectral-domain OCT scans (Leica, Chicago, IL) were acquired from healthy non-human primates. Scans of varying degree of image quality were selected for the analysis and underwent automated brightness and local contrast enhancement. 3D Frangi's vesselness filter was applied using a fixed setting for scans of all qualities. Our previously described segmentation algorithm was then used to quantify the LC microstructure. The measurements generated from the Frangi analysis and from our own conventional method were compared with a standard reference (manually segmented LC by an expert). Paired t tests were performed to compare if the differences between standard reference and conventional method are greater than the differences between standard reference and Frangi analysis. The visibility of analyzable lamina and dice coefficient were also compared to the conventional method using the same test. Results : In vivo scans acquired from 5 rhesus macaques (3 males, 1 female, aged 4.3-10.7 yrs) were used for the analysis. No significant difference was detected for LC microstructure parameters between Frangi's approach and conventional method with respect to the standard reference, except for significantly higher pore count in Frangi's method (p=0.003; Table). Furthermore, visibility (Figure) was significantly higher for the Frangi method compared to the conventional approach (p<0.001) with no difference detected for the semantic segmentation, as reflected by the dice coefficient. Conclusions : The use of Frangi analysis substantially increase the analyzable lamina while providing similar quantification of the LC microstructure compared to our previous 3D analysis method. This improves the potential for automated and thorough volumetric analysis of in vivo OCT LC image

Purpose : The ability to detect progression in eyes with advanced glaucoma is challenging because of known limitations of commonly used structural and functional parameters reaching their minimal measurable limit (floor effect) or increased measurement variability. We examined the ability of inner nuclear layer (INL) thickness measurements to demonstrate change longitudinally in eyes with early and advanced severity glaucoma. Methods : Subjects with glaucoma and >=4 visits were included in the study. Subjects in the ?Early/Moderate? group (EG) had average circumpapillary retinal nerve fiber layer (cRNFL) thicknesses >=60mum and subjects in the ?Advanced? group (AG) had average cRNFL thicknesses <=60mum. All subjects had comprehensive ophthalmic examination, Humphrey visual field (Zeiss, Dublin, CA) testing, and spectral-domain OCT (Cirrus HD-OCT; Zeiss) optic nerve head (ONH) and macula scans. Segmentation of the INL was performed using the Iowa Reference Algorithms (Retinal Image Analysis Lab, Iowa Institute for Biomedical Imaging, Iowa City, IA) and segmentation errors were manually corrected by a trained grader. Overall INL thickness along with the superior and inferior hemifields were used for analysis. Rates of progression were estimated from longitudinal OCT and visual field (VF) data using mixed effects models adjusting for baseline age, follow-up duration, and signal strength at each visit. Results : 23 eyes (23 subjects), 12 with EG and 11 with AG, were included in the study. At baseline, a statistically significant difference between groups was detected in MD, cRNFL, and GCIPL thicknesses (Table 1). In EG eyes, the rate of change was significantly different than a zero slope for cRNFL thickness, C:D ratio, and GCIPL thickness (Table 2). Inferior INL thickness was the only INL parameter showing significant rate of change. However, in the advanced group, all parameters (including both global and sectoral INL thicknesses) showed significant rate of change except for the cRNFL. Conclusions : Longitudinal measurements of INL thickness may be useful for following disease progression in subjects with advanced-stage glaucoma where cRNFL thickness is no longer useful

Purpose : The lamina cribrosa (LC) is hypothesized to be the site of initial axonal damage in glaucoma with the peripapillary retinal nerve fiber layer (RNFL) thickness is widely used as a standard metric for quantifying this damage. The purpose of this study was to determine in vivo changes in the microstructure of the LC in eyes of non-human primates (NHP) with naturally occurring RNFL thinning. Methods : Spectral-domain OCT scans (Leica, Chicago, IL) of the optic nerve head (ONH) were acquired in vivo from a colony of 50 adult rhesus monkeys, suspected of having high prevalence of naturally occurring glaucoma. The circumpapillary global and quadrant RNFL thickness was analyzed using a custom automated segmentation software. From the set of 100 eyes, the 10 eyes with the thinnest global RNFL values were selected as the study group, while 10 eyes with RNFL values around the 50 percentile were used as the control group. A previously described automated segmentation algorithm was used for LC microstructure analysis. The LC microstructure was analyzed globally and in the th following volumetric sectors: quadrants, central and peripheral lamina, and 3 depth slabs (anterior, middle, posterior; Figure). Beam thickness/pore diameter ratio (BPR) and connective tissue volume fraction (CTVF: beam volume/total volume) were calculated globally and in sectors. Results : 20 eyes (15 animals) were analyzed (Table 1). While no significant difference was detected between groups for age, weight or disc size, the study group had significantly thinner RNFL than the control group (p<0.01). The study group had significantly larger BPR and CTVF compared with the control group (Table 2). Significant sectoral differences between study and control group RNFL thickness were noted for BPR and CTVF in the nasal and temporal quadrants, central LC, and in LC depth. Across eyes, the global RNFL thickness was moderately negatively correlated only with the global CTVF (lower RNFL thickness associated with higher CTVF; r² =0.63, p=0.045). Conclusions : Eyes with thinner circumpapillary RNFL had thicker LC BPR and CTVF globally and in various sectors when compared to eyes with normal RNFL thickness. Whether these LC changes are the cause of RNFL damage or the result of remodeling of the LC requires further investigation. (Figure Presented)