Faculty Bibliography

Glaucoma is a neurodegenerative disease of the visual system and is the leading cause of irreversible blindness worldwide. To date, its pathophysiological mechanisms remain unclear. This study evaluated the feasibility of advanced diffusion magnetic resonance imaging techniques for examining the microstructural environment of the visual pathway in glaucoma. While conventional diffusion tensor imaging (DTI) showed lower fractional anisotropy and higher directional diffusivities in the optic tracts of glaucoma patients than healthy controls, diffusion kurtosis imaging (DKI) and the extended white matter tract integrity (WMTI) model indicated lower radial kurtosis, higher axial and radial diffusivities in the extra-axonal space, lower axonal water fraction, and lower tortuosity in the same regions in glaucoma patients. These findings suggest glial involvements apart from compromised axonal integrity in glaucoma. In addition, DKI and WMTI but not DTI parameters significantly correlated with clinical ophthalmic measures via optical coherence tomography and visual field perimetry testing. Taken together, DKI and WMTI provided sensitive and comprehensive imaging biomarkers for quantifying glaucomatous damage in the white matter tract across clinical severity complementary to DTI.

Purpose: Rigorous clinical testing has established that Schlemm"™s canal cross-sec-tional area (SC-CSA) is reduced in glaucomatous eyes. However, to date, it is unclear whether trabecular bypass procedures impact the morphology of the proximal aqueous outflow tract, or if the introduction of a local region of low outflow resistance adversely affects SC-CSA elsewhere, specifically presenting as SC diminution. This study quantifies changes in the morphology of the distal outflow pathway after iStent Trabecular Micro-Bypass stent (Glaukos Corp, Laguna Hills, CA, USA) implantation in living eyes by anterior segment optical coherence tomography (OCT). Design: This was a prospective observational study. Subjects: This study included six patients (eight eyes) with primary-open angle glaucoma. Methods: Patients underwent iStent placement in the nasal anterior chamber angle quadrant. OCT imaging was obtained of both nasal and temporal eye quadrants before and after surgery. For each SC parameter, an average of ten consecutive, evenly spaced measurements were manually obtained over a 1 mm segment of SC on FIJI ImageJ. Linear mixed effects modeling quantified the effect of the iStent on these parameters. Main outcome measures: Main outcome measures were changes in SC-CSA, inner-to-outer wall distance (IOD), and trabecular meshwork (TM) thickness following iStent placement. Results: Following iStent placement, total SC-CSA increased an average of 1,039.12 µm² (P = 0.05). Individually, there were no significant changes in SC-CSA in the nasal or temporal quadrants. Total SC-IOD and nasal SC-IOD increased an average of 2.35 µm (P = 0.01) and 2.96 µm (P = 0.04), respectively. There were no significant changes in temporal quadrant SC-IOD. There were no significant changes in TM thickness in either quadrant. Conclusions: Implantation of the iStent Trabecular Micro-Bypass stent significantly increases SC-IOD in the nasal quadrant at the location of implant, with no evidence of SC diminution in the temporal quadrant. It remains unclear how these observa-tions relate to the surgical efficacy of trabecular bypass procedures.

The direct analysis of 3D Optical Coherence Tomography (OCT) volumes enables deep learning models (DL) to learn spatial structural information and discover new bio-markers that are relevant to glaucoma. Downsampling 3D input volumes is the state-of-art solution to accommodate for the limited number of training volumes as well as the available computing resources. However, this limits the network's ability to learn from small retinal structures in OCT volumes. In this paper, our goal is to improve the performance by providing guidance to DL model during training in order to learn from finer ocular structures in 3D OCT volumes. Therefore, we propose an end-to-end attention guided 3D DL model for glaucoma detection and estimating visual function from retinal structures. The model consists of three pathways with the same network architecture but different inputs. One input is the original 3D-OCT cube and the other two are computed during training guided by the 3D gradient class activation heatmaps. Each pathway outputs the class-label and the whole model is trained concurrently to minimize the sum of losses from three pathways. The final output is obtained by fusing the predictions of the three pathways. Also, to explore the robustness and generalizability of the proposed model, we apply the model on a classification task for glaucoma detection as well as a regression task to estimate visual field index (VFI) (a value between 0 and 100). A 5-fold cross-validation with a total of 3782 and 10,370 OCT scans is used to train and evaluate the classification and regression models, respectively. The glaucoma detection model achieved an area under the curve (AUC) of 93.8% compared with 86.8% for a baseline model without the attention-guided component. The model also outperformed six different feature based machine learning approaches that use scanner computed measurements for training. Further, we also assessed the contribution of different retinal layers that are relevant to glaucoma. The VFI estimation model achieved a Pearson correlation and median absolute error of 0.75 and 3.6%, respectively, for a test set of size 3100 cubes.

Purpose/UNASSIGNED:The prevailing theory about the function of lamina cribrosa (LC) connective tissues is that they provide structural support to adjacent neural tissues. Missing connective tissues would compromise this support and therefore are regarded as "LC defects", despite scarce actual evidence of their role. We examined how so-called LC defects alter IOP-related mechanical insult to the LC neural tissues. Methods/UNASSIGNED:We built numerical models incorporating LC microstructure from polarized light microscopy images. To simulate LC defects of varying sizes, individual beams were progressively removed. We then compared intraocular pressure (IOP)-induced neural tissue deformations between models with and without defects. To better understand the consequences of defect development, we also compared neural tissue deformations between models with partial and complete loss of a beam. Results/UNASSIGNED:The maximum stretch of neural tissues decreased non-monotonically with defect size. Maximum stretch in the model with the largest defect decreased by 40% in comparison to the model with no defects. Partial loss of a beam increased the maximum stretch of neural tissues in its adjacent pores by 162%, compared with 63% in the model with complete loss of a beam. Conclusions/UNASSIGNED:Missing LC connective tissues can mitigate IOP-induced neural tissue insult, suggesting that the role of the LC connective tissues is more complex than simply fortifying against IOP. The numerical models further predict that partial loss of a beam is biomechanically considerably worse than complete loss of a beam, perhaps explaining why defects have been reported clinically but partial beams have not.

Purpose : Predicting progression of primary open angle glaucoma (POAG) continues to be a challenge. Recent studies have shown that macular parameters may predict glaucoma progression and disease onset equally as well or even better than the traditionally used optic nerve head (ONH) retinal nerve fiber layer (RNFL) thickness. We performed a retrospective longitudinal study to identify structural parameters which best predicted visual field (VF) loss. Methods : Subjects with POAG with at least 5 qualified Cirrus OCT (Zeiss, Dublin, CA) macular and optic nerve head (ONH) scans and 5 qualified 24-2 Humphrey VF tests (Zeiss) were enrolled to this study. All parameters from the OCT's report of the macula and ONH scans were used in the analysis. We identified subjects who were OCT ONH progressors or VF progressors using both event and trend based definitions of progression based on the Guided Progression Analyses. Students t-test was used to assess differences in baseline characteristics between groups, and mixed-effects model was used in longitudinal analysis to compare the rate of parameter change between groups. Results : 263 eyes (180 subjects) with a mean follow-up time of 2.4+/-1.8 years were included in the study. Twenty-three eyes were identified as ONH progressors, 25 eyes were identified as VF progressors, and of these, 6 eyes were both VF and ONH progressors. At baseline, only macular average RNFL and macular inferior RNFL were significantly thinner in ONH progressors compared to non-progressors. Between the VF progressors and non-progressor groups, all OCT parameters were significantly thinner at baseline in the progressors except for four focal macular measurements and disc area. Mixed-effects modeling showed that both focal macular and ONH parameters thinned at a significantly greater rate in VF progressors compared to non-progressors (Table 1). Conclusions : Focal macular retinal segmentations and focal ONH parameters thin significantly more per year in VF progressors, and VF progressors tend to have thinner baseline structural parameters

Purpose : Gradient class activation maps (grad-CAM) generated by convolutional neural networks (CNN) have qualitatively indicated that these networks are able to identify important regions in OCT scans. Here, we quantitatively analyse these regions to improve our understanding of the CNN decision making process when detecting glaucoma in OCT volumes. Methods : A total of 1110 OCT (Cirrus HD-OCT, Zeiss, Dublin, Ca) scans from both eyes of 624 subjects (139 healthy and 485 glaucomatous patients (POAG)). An end-to-end 3D-CNN network was trained directly on 3D-volumes for glaucoma detection. Grad-CAM was implemented to highlight structures in the volumes that the network relied on. Grad-CAM heatmaps were generated for 3 different convolutional layers and quantitatively validated by occluding the regions with the highest grad-CAM weights (12.5% of original input volumes) and then evaluating the performance drop. Further, 8-retinal layers segmentation method was used to compute the average heatmap weights for each segmented layer separately, and used to identify the layers that were deemed as important for the task. Results : The model achieved an AUC of 0.97 for the test set (110 scans). Occlusion resulted in a 40% drop in performance (Fig.1). The RNFL and photoreceptors showed the highest median weights for grad-CAM heatmaps (0.1 and 0.2, respectively). The retinal pigment epithelium (RPE) and photoreceptors showed higher weights in the glaucomatous scans (Fig.2-a). RNFL had wider range of weights in healthy cases versus POAG ones. Analysis of the B-scans showed that central part around the optic disc (# 85-135) had the highest contribution to the network decision and the heatmap weights were much higher in glaucoma cases than healthy ones across all B-scans (Fig.2-b). Conclusions : The occlusion experiment indicates that the regions identified by the grad-CAMs are in fact pertinent to the glaucoma detection task. The increased emphasis on the photoreceptors in the glaucoma cases may be attributed to the atrophy in the superficial layers which in turn increased the brightness of this structure. This technique can be used to identify new biomarkers learned for other ocular diseases

Purpose : In mild glaucoma, RNFL thinning and visual field (VF) loss are often localized, but structure-function modeling is impeded by variability due to individual eye anatomy. We perform high-resolution spatial correlations of RNFLT maps for each VF location to identify relevant areas and study further improvements by geometrically aligning RNFLT maps based on artery trajectories. Methods : In 419 SITA Standard 24-2 Humphrey VFs with at most mild glaucoma (mean deviation >=-3dB) with accompanying circumpapillary Cirrus HD-OCT RNFLT maps, we computed pixel-wise correlations (52 VF locations x 40401 pixels). We then performed an alignment operation, ensuring that the two major retinal arteries follow the same lines in all scans. We piecewise linearly approximated the trajectories of the arteries on 4 concentric circles around ONH (Fig. 1a), determined the necessary rotation for each pixel, and morphed the images accordingly (Fig. 1b). Results : For the pre-alignment RNFLT (correlation maps Fig. 2 top) we observed: (1) relatively high correlations (max 0.29); (2) most of the high-correlation regions are highly localized around the median trajectories of the major arteries at most VF locations, possibly due to the stacked character of the fiber bundles close to ONH, which impedes precise spatial mapping to the VF. This observation suggests general retinal vulnerability zones rather than highly VF location-specific areas as assumed by many previous structure-function models. Accordingly, morphing the RNFLT maps by aligning the eyespecific artery locations increased the maximal correlations on 25 of the 52 VF locations (Fig. 2 bottom, marked in green), particularly in nasal and inferior VF, with improvements of up to 0.1 (inferior arcuate region of VF). At many locations, aligned vulnerability areas become substantially more conspicuous (e.g. the location enlarged on the top left) and might have been missed without aligning. Conclusions : High-resolution structure-function correlations reveal retinal vulnerability zones in mild glaucoma. At many VF locations, these zones become better correlated with VF regions when RNFLT maps are aligned along the arteries. Specific attention to RNFL thinning in these zones in glaucoma suspects may improve the detection of initial VF loss glaucoma

Purpose : Development and pre-clinical testing of glaucoma neurotherapeutics have been obfuscated by limited experimental models that provide chronic elevation of intraocular pressure (IOP) while preserving optical media clarity for structural and functional assessments over time. In this study, we developed an in vivo model system involving the use of non-invasive tonometry, optical coherence tomography (OCT) and optokinetics to characterize retinal integrity and visual behavior in a novel hydrogel-induced chronic IOP elevation model. Methods : Six adult C57BL/6J mice underwent unilateral intracameral injection of an optically clear, chemically cross-linked hydrogel composed of hyaluronic acid functionalized with vinyl sulfone and thiol groups. IOP was measured with a rebound tonometer at baseline and 1, 3, 7, 10 and 14 days after hydrogel injection. The optic nerve head (ONH) region was scanned for each eye using OCT at baseline and 2 weeks after injection, and total retinal thickness (TRT) was measured within a 0.26-0.36 mm radius ring centered on the ONH using custom-written software (Fig 1). Visual acuity (VA) was measured for each eye using an optokinetic virtual-reality system at baseline and 2 weeks after injection. Data are presented as mean+/-SEM. Results : Intracameral hydrogel injection resulted in mild-to-moderate IOP elevation throughout the 2-week experimental period (Fig 2a). TRT in the hydrogel-injected eye was 10.06+/-3.61% thinner at 2 weeks post-injection compared to baseline (p<0.01) (Fig 2b). IOP elevation also led to a decline in VA by 58.12+/-7.22% at 2 weeks post-injection compared to baseline (p<0.001) (Fig 2c). Interestingly, among the hydrogel-injected eyes, cumulative IOP measured from 0 to 14 days post-injection did not correlate with TRT or VA (p>0.05) (Fig 2d-e), whereas TRT was positively associated with VA at 2 weeks post-injection (r=0.824, p<0.05) (Fig 2f). No significant change in IOP, TRT or VA was found in the non-injected eye. Conclusions : An in vivo glaucoma model system was developed that showed a positive correlation between retinal thinning and visual behavioral deficits after chronic IOP elevation. The weak association between cumulative IOP and TRT or VA suggests additional factors apart from IOP level in contributing to glaucomatous damage after chronic IOP elevation

Purpose : Oxygen concentration in retinal blood vessels (sO ) can be critical biomarkers for diabetic retinopathy and glaucoma, leading causes of blindness worldwide. We previously demonstrated sO2measurements in rodent and human retinas with spectroscopic visible-light optical coherence tomography (vis-OCT). However, reliable measurements of sO2in a clinical setting remains an open challenge due to constraints on light exposure, imaging time, patient motion, and variation in eye geometry. Spectral calibration to optimize sO2measurements under these non-ideal imaging conditions is needed. Here, we investigate, develop, and implement such calibration. Methods : We developed vis-OCT processing software to optimize sO2measurements in humans. First, we identified an optimal spectral range for spectral measurement in which sO2was most stable. Next, we developed methods to account for alterations induced by the imaging system and eye optics. Specifically, we accounted for depth-dependent variations in the measured spectrum, such as absorption contrast, spectrally-dependent roll-off, chromatic aberrations, and eye morphology. We then imaged the retinas of 12 healthy subjects aged 22 to 60 at Northwestern Medical Hospital in Chicago, IL, and Langone Medical Center in New York, NY. All imaging was approved by the respective IRBs and strictly adhered to the Declaration of Helsinki. Light exposure in the eye was no higher than 250 muW and imaging time was no longer than 5 s. We extracted sO2from vessels larger than 50 mum in diameter using an automated version of our vis-OCT processing software. Results : We measured the sO2in 89 vessels (53 arteries and 36 veins). We found the mean sO2in arteries was 97.70 +/-4.75 % in arteries and mean sO2in veins was 53.11 +/-6.85 %. Conclusions : We developed analytical methods for depth-dependent alterations to the measured spectrum in vis-OCT retinal oximetry. Our measurements yielded spectra that are highly consistent with those reported in literature, despite variations in imaging conditions. Our results indicate a clear path forward for clinical adoption of vis-OCT