A mathematical model for describing the retinal nerve fiber bundle trajectories in the human eye: average course, variability, and influence of refraction, optic disc size and optic disc position

Exploring the relationship between 24-2 visual field and widefield optical coherence tomography data across healthy, glaucoma suspect and glaucoma eyes

PURPOSE

To utilise ganglion cell-inner plexiform layer (GCIPL) measurements acquired using widefield optical coherence tomography (OCT) scans spanning 55° × 45° to explore the link between co-localised structural parameters and clinical visual field (VF) data.

METHODS

Widefield OCT scans acquired from 311 healthy, 268 glaucoma suspect and 269 glaucoma eyes were segmented to generate GCIPL thickness measurements. Estimated ganglion cell (GC) counts, calculated from GCIPL measurements, were plotted against 24-2 SITA Faster visual field (VF) thresholds, and regression models were computed with data categorised by diagnosis and VF status. Classification of locations as VF defective or non-defective using GCIPL parameters computed across eccentricity- and hemifield-dependent clusters was assessed by analysing areas under receiver operating characteristic curves (AUROCCs). Sensitivities and specificities were calculated per diagnostic category.

RESULTS

Segmented linear regression models between GC counts and VF thresholds demonstrated higher variability in VF defective locations relative to non-defective locations (mean absolute error 6.10-9.93 dB and 1.43-1.91 dB, respectively). AUROCCs from cluster-wide GCIPL parameters were similar across methods centrally (p = 0.06-0.84) but significantly greater peripherally, especially when considering classification of more central locations (p < 0.0001). Across diagnoses, cluster-wide GCIPL parameters demonstrated variable sensitivities and specificities (0.36-0.93 and 0.65-0.98, respectively), with the highest specificities observed across healthy eyes (0.73-0.98).

CONCLUSIONS

Quantitative prediction of VF thresholds from widefield OCT is affected by high variability at VF defective locations. Prediction of VF status based on cluster-wide GCIPL parameters from widefield OCT could become useful to aid clinical decision-making in appropriately targeting VF assessments.

Maculopapillary Bundle Degeneration in Optic Neuropathies

PURPOSE OF REVIEW

Degeneration of the maculopapillary bundle (MPB) is a prominent feature in a spectrum of optic neuropathies. MPB-selective degeneration is seen in specific conditions, such as nutritional and toxic optic neuropathies, Leber hereditary optic neuropathy (LHON), and dominant optic atrophy (DOA). Despite their distinct etiologies and clinical presentations, which encompass variations in age of incidence and monocular or binocular onset, these disorders share a core molecular mechanism: compromised mitochondrial homeostasis. This disruption is characterized by dysfunctions in mitochondrial metabolism, biogenesis, and protein synthesis. This article provides a comprehensive understanding of the MPB's role in optic neuropathies, emphasizing the importance of mitochondrial mechanisms in the pathogenesis of these conditions.

RECENT FINDINGS

Optical coherence tomography studies have characterized the retinal nerve fiber layer changes accompanying mitochondrial-affiliated optic neuropathies. Selective thinning of the temporal optic nerve head is preceded by thickening in early stages of these disorders which correlates with reductions in macular ganglion cell layer thinning and vascular atrophy. A recently proposed mechanism underpinning the selective atrophy of the MPB involves the positive feedback of reactive oxygen species generation as a common consequence of mitochondrial dysfunction. Additionally, new research has revealed that the MPB can undergo degeneration in the early stages of glaucoma, challenging the historically held belief that this area was not involved in this common optic neuropathy. A variety of anatomical risk factors influence the propensity of glaucomatous MPB degeneration, and cases present distinct patterns of ganglion cell degeneration that are distinct from those observed in mitochondria-associated diseases. This review synthesizes clinical and molecular research on primary MPB disorders, highlighting the commonalities and differences in their pathogenesis.

KEY POINTS (BOX)

1. Temporal degeneration of optic nerve fibers accompanied by cecocentral scotoma is a hallmark of maculopapillary bundle (MPB) degeneration. 2. Mechanisms of MPB degeneration commonly implicate mitochondrial dysfunction. 3. Recent research challenges the traditional belief that the MPB is uninvolved in glaucoma by showing degeneration in the early stages of this common optic neuropathy, yet with features distinct from other MPB-selective neuropathies. 4. Reactive oxygen species generation is a mechanism linking mitochondrial mechanisms of MPB-selective optic neuropathies, but in-vivo and in-vitro studies are needed to validate this hypothesis.

Visual field progression patterns in the ocular hypertension treatment study correspond to vulnerability regions of the disc

OBJECTIVES

To determine the locations on the 24-2 visual field (VF) testing grid that are most likely to progress in patients with ocular hypertension (OHTN). Based on a structural model of superior and inferior areas of relative vulnerability at the optic disc, we hypothesized that the nasal and paracentral regions are more prone to show a reduction in sensitivity.

METHODS

Posthoc analysis of data collected in phases 1 and 2 of the Ocular Hypertension Treatment Study (OHTS). A pointwise analysis was applied to determine the progression patterns in the early and delayed treatment groups. Each group's progression rate and frequency were calculated for each of the 52 locations corresponding to the 24-2 VF strategy, using trend- and event-based analyses, respectively.

RESULTS

For the event-based analysis, the events were most commonly found in the nasal and paracentral regions. The same regions, with some modest variation, were found to have the fastest rates of progression (ROP) measured with trend analysis. A similar pattern of progression was observed in both the early and delayed treatment groups. The difference in event rates and ROP between the early and delayed treatment groups was also greatest in the nasal and paracentral regions.

CONCLUSIONS

Development of VF loss in ocular hypertensive eyes appears to be consistent with the vulnerability zones previously described in glaucomatous eyes with established VF loss. Ocular hypotensive treatment likely helps to slow the rate of progression in these regions. This suggests that careful monitoring of these locations may be useful.

Progressive Changes in the Neuroretinal Rim and Retinal Nerve Fiber Layer in Glaucoma: Impact of Baseline Values and Floor Effects

PURPOSE

To determine whether more severe baseline damage impedes measurement of minimum rim width (MRW) and peripapillary retinal nerve fiber layer thickness (RNFLT) change in glaucoma patients because of a floor effect.

DESIGN

Prospective, longitudinal cohort study in a hospital-based setting.

PARTICIPANTS

The study included patients with open-angle glaucoma and healthy control subjects. Participants had at least 5 years of follow-up with OCT every 6 months.

METHODS

Baseline global and sectorial MRW and RNFLT values were classified as within normal limits, borderline, or outside normal limits based on reference normative values. Regression analysis was used to determine the magnitude and significance of MRW and RNFLT change. Additionally, the follow-up period for each participant was divided into 2 equal halves (first and second periods) to determine whether there was attenuation of MRW and RNFLT change with follow-up time.

MAIN OUTCOME MEASURES

Rates of global and sectoral MRW and RNFLT changes (slopes).

RESULTS

A total of 97 patients with glaucoma (median age, 70.3 years) and 42 healthy subjects (median age, 64.8 years) were followed for a median of 6.9 years and 7.0 years, respectively. The median mean deviation of the visual field in glaucoma patients was -4.30 decibels (dB) (interquartile range, -7.81 to -2.06 dB; range, -20.68 to 1.37 dB). Statistically significant changes in global and sectoral MRW and RNFLT were detected across all baseline classifications; however, there was a tendency for less change with increasing baseline damage. In glaucoma patients, RNFLT slopes, but not MRW slopes, were significantly more positive (less change) in the second period compared with the first. There were also no differences in MRW or RNFLT slopes in the first and second periods in healthy subjects.

CONCLUSIONS

Significant MRW and RNFLT changes were detected at all levels of baseline damage. However, an attenuation in the rate of RNFLT change compared with MRW indicates an earlier floor effect in RNFLT measurements globally and in equivalent sectors. Because the axonal component of these measurements should be equivalent, our results suggest important differences in tissue remodeling at the level of the optic nerve head and peripapillary retina.

FINANCIAL DISCLOSURE(S)

Proprietary or commercial disclosure may be found after the references.

Characterisation of the normal human ganglion cell-inner plexiform layer using widefield optical coherence tomography

PURPOSE

To describe variations in ganglion cell-inner plexiform layer (GCIPL) thickness in a healthy cohort from widefield optical coherence tomography (OCT) scans.

METHODS

Widefield OCT scans spanning 55° × 45° were acquired from 470 healthy eyes. The GCIPL was automatically segmented using deep learning methods. Thickness measurements were extracted after correction for warpage and retinal tilt. Multiple linear regression analysis was applied to discern trends between global GCIPL thickness and age, axial length and sex. To further characterise age-related change, hierarchical and two-step cluster algorithms were applied to identify locations sharing similar ageing properties, and rates of change were quantified using regression analyses with data pooled by cluster analysis outcomes.

RESULTS

Declines in widefield GCIPL thickness with age, increasing axial length and female sex were observed (parameter estimates -0.053, -0.436 and -0.464, p-values <0.001, <0.001 and 0.02, respectively). Cluster analyses revealed concentric, slightly nasally displaced, horseshoe patterns of age-related change in the GCIPL, with up to four statistically distinct clusters outside the macula. Linear regression analyses revealed significant ageing decline in GCIPL thickness across all clusters, with faster rates of change observed at central locations when expressed as absolute (slope = -0.19 centrally vs. -0.04 to -0.12 peripherally) and percentage rates of change (slope = -0.001 centrally vs. -0.0005 peripherally).

CONCLUSIONS

Normative variations in GCIPL thickness from widefield OCT with age, axial length and sex were noted, highlighting factors worth considering in further developments. Widefield OCT has promising potential to facilitate quantitative detection of abnormal GCIPL outside standard fields of view.

Training in cortically-blind fields confers patient-specific benefit against retinal thinning after occipital stroke

Purpose

Damage to the adult primary visual cortex (V1) causes vision loss in the contralateral hemifield, initiating a process of trans-synaptic retrograde degeneration (TRD). Here, we examined retinal correlates of TRD using a new metric to account for global changes in inner retinal thickness, and asked if perceptual training in the intact or blind field impacts its progression.

Methods

We performed a meta-analysis of optical coherence tomography (OCT) data in 48 participants with unilateral V1 stroke and homonymous visual defects, who completed clinical trial NCT03350919. After measuring the thickness of the macular ganglion cell and inner plexiform layers (GCL-IPL), and the peripapillary retinal nerve fiber layer (RNFL), we computed individual laterality indices (LI) at baseline and after ~6 months of daily motion discrimination training in the intact- or blind-field. Increasingly positive LI denoted greater layer thinning in retinal regions affected versus unaffected by the cortical damage.

Results

Pre-training, the affected GCL-IPL and RNFL were thinner than their unaffected counterparts, generating LI values positively correlated with time since stroke. Participants trained in their intact-field exhibited increased LIGCL-IPL. Those trained in their blind-field had no significant change in LIGCL-IPL. LIRNFL did not change in either group.

Conclusions

Relative shrinkage of the affected versus unaffected macular GCL-IPL can be reliably measured at an individual level and increases with time post-V1 stroke. Relative thinning progressed during intact-field training, but appeared to be halted by training within the blind field, suggesting a potentially neuroprotective effect of this simple behavioral intervention.

Association between Glaucoma Progression in Macular Ganglion Cell Complex and Disc Hemorrhage: Differences between Superior and Inferior Hemiretinas

Disc hemorrhage (DH) is often associated with glaucoma progression. A vertically asymmetrical pattern is typical of glaucoma progression, but it remains unclear whether the association between DH and glaucoma progression differs between the superior and inferior hemiretinas. We compared the thickness changes of the macular ganglion cell complex (GCC) in the hemiretinas of normal-tension glaucoma patients with or without DH, as well as between hemiretinas positive and negative for DH, during five years. Both the superior and inferior hemiretinas in the DH-positive group had a more negative GCC thickness slope in association with more DH counts compared to those in the DH-negative group. Conversely, only the inferior hemiretina exhibited a significant relationship between GCC thickness slope and DH counts when hemiretinas positive and negative for DH in the DH-positive group were compared. In the superior hemifield, the slope of the total deviation changes in the DH-positive hemifield of the DH-positive group was more negative compared to that of the DH-negative group. The association between DH and glaucoma progression in the macular GCC may be stronger in the inferior hemiretina, suggesting that more attention should be paid to DH in the inferior disc area as a sign of glaucoma progression.

Structure-function assessment in glaucoma based on perimetric sensitivity and en face optical coherence tomography images of retinal nerve fiber bundles

Many studies have assessed structure-function relations in glaucoma, but most without topographical comparison across the central 30°. We present a method for assessing structure-function relations with en face images of retinal nerve fiber layer (RNFL) bundles allowing topographical comparison across much of this retinal area. Forty-four patients with glaucoma (median age 61 years) were recruited and tested with Optical Coherence Tomography (OCT) and perimetry. Six rectangular volume scans were gathered, and then montaged to provide en face views of the RNFL bundles. We calculated the proportion of locations showing a perimetric defect that also showed an en face RNFL defect; and the proportion of locations falling on an RNFL defect that also showed a perimetric defect. A perimetric defect for a location was defined as a total deviation (TD) value equal to or deeper than -4 dB. We found that the median (IQR) number of locations with abnormal RNFL bundle reflectance that also had abnormal TD was 78% (60%) and for locations with abnormal TD that also had abnormal RNFL bundle reflectance was 75% (44%). We demonstrated a potential approach for structure-function assessment in glaucoma by presenting a topographic reflectance map, confirming results of previous studies and including larger retinal regions.

On the protective role of the blood vessels in glaucomatous damage: A transversal study

PURPOSE

To corroborate whether vessels on the surface of the optic nerve head can provide protection against the loss of underlying axons in subjects with manifest glaucoma.

METHODS

In this pilot study, thirty-six glaucomatous eyes with a perimetric defect in the Bjerrum area were included. The retinal nerve fiber layer (RNFL) thickness was measured in each of the sectors of the clock-hour map obtained by Cirrus HD-OCT considering the presence or absence of blood vessels. These sectors were related with their corresponding areas of the retina examined in the visual field using a mathematical model of the retina introduced by Jansonius, in order to determine the values of threshold sensitivity in those areas in the presence or absence of vessels.

RESULTS

We corroborated the protective role of the blood vessel for peripapillary RNFL thickness of clock-hour 12 despite obtaining a p-value (p = 0.023; w = 228.5) close to the acceptance zone (p ≥ 0.05). The mean ± standard deviation with vessel and without vessel were 70.95 ± 24.35 and 88.46 ± 23.96, respectively. No differences were found between the mean values of threshold sensitivity to the presence or absence of blood vessels in each of the sectors considered.

CONCLUSIONS

Our findings do not allow us to affirm that there is an association between the presence of a vessel and protection against glaucomatous damage in subjects with an advanced manifestation of the disease. In the future, more extensive studies are needed to study this relationship in subjects with early glaucoma.

Differing Associations between Optic Nerve Head Strains and Visual Field Loss in Patients with Normal- and High-Tension Glaucoma

PURPOSE

To study the associations between optic nerve head (ONH) strains under intraocular pressure (IOP) elevation with retinal sensitivity in patients with glaucoma.

DESIGN

Clinic-based cross-sectional study.

PARTICIPANTS

Two hundred twenty-nine patients with primary open-angle glaucoma (subdivided into 115 patients with high-tension glaucoma [HTG] and 114 patients with normal-tension glaucoma [NTG]).

METHODS

For 1 eye of each patient, we imaged the ONH using spectral-domain OCT under the following conditions: (1) primary gaze and (2) primary gaze with acute IOP elevation (to approximately 35 mmHg) achieved through ophthalmodynamometry. A 3-dimensional strain-mapping algorithm was applied to quantify IOP-induced ONH tissue strain (i.e., deformation) in each ONH. Strains in the prelaminar tissue (PLT), the retina, the choroid, the sclera, and the lamina cribrosa (LC) were associated (using linear regression) with measures of retinal sensitivity from the 24-2 Humphrey visual field test (Carl Zeiss Meditec). This was performed globally, then locally according to a previously published regionalization scheme.

MAIN OUTCOME MEASURES

Associations between ONH strains and values of retinal sensitivity from visual field testing.

RESULTS

For patients with HTG, we found (1) significant negative linear associations between ONH strains and retinal sensitivity (P < 0.001; on average, a 1% increase in ONH strains corresponded to a decrease in retinal sensitivity of 1.1 decibels [dB]), (2) that high-strain regions colocalized with anatomically mapped regions of high visual field loss, and (3) that the strongest negative associations were observed in the superior region and in the PLT. In contrast, for patients with NTG, no significant associations between strains and retinal sensitivity were observed except in the superotemporal region of the LC.

CONCLUSIONS

We found significant negative associations between IOP-induced ONH strains and retinal sensitivity in a relatively large glaucoma cohort. Specifically, patients with HTG who experienced higher ONH strains were more likely to exhibit lower retinal sensitivities. Interestingly, this trend in general was less pronounced in patients with NTG, which could suggest a distinct pathophysiologic relationship between the two glaucoma subtypes.

Comparison of machine learning approaches for structure-function modeling in glaucoma

To evaluate machine learning (ML) approaches for structure-function modeling to estimate visual field (VF) loss in glaucoma, models from different ML approaches were trained on optical coherence tomography thickness measurements to estimate global VF mean deviation (VF MD) and focal VF loss from 24-2 standard automated perimetry. The models were compared using mean absolute errors (MAEs). Baseline MAEs were obtained from the VF values and their means. Data of 832 eyes from 569 participants were included, with 537 Asian eyes for training, and 148 Asian and 111 Caucasian eyes set aside as the respective test sets. All ML models performed significantly better than baseline. Gradient-boosted trees (XGB) achieved the lowest MAE of 3.01 (95% CI: 2.57, 3.48) dB and 3.04 (95% CI: 2.59, 3.99) dB for VF MD estimation in the Asian and Caucasian test sets, although difference between models was not significant. In focal VF estimation, XGB achieved median MAEs of 4.44 [IQR 3.45-5.17] dB and 3.87 [IQR 3.64-4.22] dB across the 24-2 VF for the Asian and Caucasian test sets and was comparable to VF estimates from support vector regression (SVR) models. VF estimates from both XGB and SVR were significantly better than the other models. These results show that XGB and SVR could potentially be used for both global and focal structure-function modeling in glaucoma.

Test of a Retinal Nerve Fiber Bundle Trajectory Model Using Eyes With Glaucomatous Optic Neuropathy

Purpose

To test a model of retinal nerve fiber bundle trajectories that predicts the arcuate-shaped patterns seen on optical coherence tomography (OCT) retinal nerve fiber layer (RNFL) probability/deviation maps (p-maps) in glaucomatous eyes.

Methods

Thirty-one glaucomatous eyes from a database of 250 eyes had clear arcuate-shaped patterns on RNFL p-maps derived from an OCT cube scan. The borders of the arcuate patterns were extracted from the RNFL p-maps. Next, the trajectories from an arcuate model were compared against these borders via a normalized root-mean-square difference analysis. The model's parameter, β, was varied, and the best-fitting, initial clock-hour position of the trajectory to the border was found for each β. Finally, the regions, as determined by the arcuate border's best-fit, initial clock-hour positions, were compared against the abnormal regions on the circumpapillary retinal nerve fiber layer (cpRNFL) profile.

Results

The arcuate model's mean β_Sup and β_Inf parameters minimized large differences between the trajectories and the arcuate borders on the RNFL p-maps. Furthermore, on average, 68% of the cpRNFL regions defined by the arcuate border's best-fit, initial clock-hour positions were abnormal (i.e., below the ≤5% threshold).

Conclusions

The arcuate model performed well in predicting the borders of arcuate patterns seen on RNFL p-maps. It also predicted the associated abnormal regions of the cpRNFL thickness plots.

Translational Relevance

This model should prove useful in helping clinicians understand topographical comparisons among different OCT representations and should improve structure-structure, as well as structure-function agreement analyses.

Immersive Virtual Reality Simulations of Bionic Vision

Bionic vision uses neuroprostheses to restore useful vision to people living with incurable blindness. However, a major outstanding challenge is predicting what people "see" when they use their devices. The limited field of view of current devices necessitates head movements to scan the scene, which is difficult to simulate on a computer screen. In addition, many computational models of bionic vision lack biological realism. To address these challenges, we present VR-SPV, an open-source virtual reality toolbox for simulated prosthetic vision that uses a psychophysically validated computational model to allow sighted participants to "see through the eyes" of a bionic eye user. To demonstrate its utility, we systematically evaluated how clinically reported visual distortions affect performance in a letter recognition and an immersive obstacle avoidance task. Our results highlight the importance of using an appropriate phosphene model when predicting visual outcomes for bionic vision.

Cluster Formation for Analyses of Glaucomatous Visual Field Defects in Central 10-2 Visual Field in Normal Tension Glaucoma Eyes

Purpose

To develop a cluster system to analyze the retinal sensitivity loss of 68 test points in the central 10 degrees of standard automated perimetry (SAP) in eyes with normal tension glaucoma (NTG).

Patients and Methods

Patients with NTG who met the following criteria were included: visual acuity ≥0.7, SAP-derived mean deviation ≥−15 dB, and pattern deviation probability plots with at least one point with a probability of <0.5% and/or two or more contiguous points with a probability of <1% that did not cross the horizontal meridian in the central 12 points of the 24-2 test points. SAP with the Swedish Interactive Threshold Algorithm Standard (SITA-S) 10-2 program (10-2) was performed within 6 months of the SITA-S 24-2. The averaged total deviation (TD) for each of the 68 test points in the 10-2 was calculated. Hierarchical cluster analyses were performed based on the deviation of the TDs of the test points, and a dendrogram was created. The number of clusters was determined following the Sturges’ rule.

Results

One hundred and twenty-six eyes of 126 patients (61.9±11.4 years) were studied. Hierarchical cluster analysis of the TD values statistically obtained a dendrogram that divided the 68 test points into 7 clusters. Of these 7 clusters, 21 points belonging to the clusters in the papillomacular region included cluster 5. Cluster 5 was distributed above and below the horizontal meridian, which does not agree with the course of the retinal nerve fiber layer (RNFL).

Conclusion

The hierarchical cluster analysis of the TD values stratified the 68 test points of the 10-2 into seven clusters. Considering the course of the RNFL, cluster 5 was divided into clusters of 5a and 5b, and consequently eight clusters were considered to be appropriate for detecting glaucomatous visual field defects in the central 10 degrees in NTG eyes.

Detecting glaucoma with only OCT: Implications for the clinic, research, screening, and AI development

A method for detecting glaucoma based only on optical coherence tomography (OCT) is of potential value for routine clinical decisions, for inclusion criteria for research studies and trials, for large-scale clinical screening, as well as for the development of artificial intelligence (AI) decision models. Recent work suggests that the OCT probability (p-) maps, also known as deviation maps, can play a key role in an OCT-based method. However, artifacts seen on the p-maps of healthy control eyes can resemble patterns of damage due to glaucoma. We document in section 2 that these glaucoma-like artifacts are relatively common and are probably due to normal anatomical variations in healthy eyes. We also introduce a simple anatomical artifact model based upon known anatomical variations to help distinguish these artifacts from actual glaucomatous damage. In section 3, we apply this model to an OCT-based method for detecting glaucoma that starts with an examination of the retinal nerve fiber layer (RNFL) p-map. While this method requires a judgment by the clinician, sections 4 and 5 describe automated methods that do not. In section 4, the simple model helps explain the relatively poor performance of commonly employed summary statistics, including circumpapillary RNFL thickness. In section 5, the model helps account for the success of an AI deep learning model, which in turn validates our focus on the RNFL p-map. Finally, in section 6 we consider the implications of OCT-based methods for the clinic, research, screening, and the development of AI models.

Policy-Driven, Multimodal Deep Learning for Predicting Visual Fields from the Optic Disc and Optical Coherence Tomography Imaging

PURPOSE

To develop and validate a deep learning (DL) system for predicting each point on visual fields (VF) from disc and optical coherence tomography (OCT) imaging and derive a structure-function mapping.

DESIGN

Retrospective, cross-sectional database study PARTICIPANTS: 6437 patients undergoing routine care for glaucoma in three clinical sites in the UK.

METHODS

OCT and infrared reflectance (IR) optic disc imaging was paired with the closest VF within 7 days. Efficient-Net B2 was used to train two single modality DL models to predict each of the 52 sensitivity points on the 24-2 VF pattern. A policy DL model was designed and trained to fuse the two model predictions.

MAIN OUTCOME MEASURES

Pointwise Mean Absolute Error (PMAE) RESULTS: A total of 5078 imaging to VF pairs were used as a held-out test set to measure the final performance. The improvement in PMAE with the policy model was 0.485 [0.438, 0.533] dB compared to the IR image of the disc alone and 0.060 [0.047, 0.073] dB compared to the OCT alone. The improvement with the policy fusion model was statistically significant (p < 0.0001). Occlusion masking shows that the DL models learned the correct structure function mapping in a data-driven, feature agnostic fashion.

CONCLUSIONS

The multimodal, policy DL model performed the best; it provided explainable maps of its confidence in fusing data from single modalities and provides a pathway for probing the structure-function relationship in glaucoma.

Combining vascular and nerve fiber layer thickness measurements to model glaucomatous focal visual field loss

We compare the focal structure-function correlation of structural measurements of peripapillary retinal nerve fiber layer thickness (RNFL-T) using optical coherence tomography (OCT), capillary density (CD) measurements using OCT-angiography (OCT-A), or a combination of both, with visual field deviation (VFD) in early to advanced glaucoma. Primary open angle glaucoma patients (n = 46, mean ± SD age: 67 ± 10 years; VF mean deviation: -10.41 ± 6.76 dB) were included in this cross-sectional study. We performed 30-2 standard automated perimetry OCT (3.5-mm diameter ring scan) and 15°×15° OCT-A (superficial vascular complex slab). Based on a nerve fiber trajectory model, each VF test spot was assigned to an OCT-A wedge and an OCT ring-sector. Two univariate linear models (M_v and M_t ) using either CD-based vascular (M_v ) or RNFL-T-based thickness information (M_t ) and one multivariate model using both (M_v:t ) were compared in their associations with measured focal VFD, which were higher for the multivariate model M_v:t (mean ± SD correlation coefficient: 0.710 ± 0.086) than for either nested model (0.627 ± 0.078 for M_v and 0.578 ± 0.095 for M_t ). Using a focal visual field approach, the combination of RNFL-T and CD showed better structure-function correlations than thickness or vascular information only.

Automatic retinal nerve fiber bundle tracing based on large field of view polarization sensitive OCT data

A technique to accurately estimate trajectories of retinal nerve fiber bundles (RNFB) in a large field of view (FOV) image covering 45° is described. The method utilizes stitched projections of polarization-sensitive optical coherence tomography (PS-OCT) data, as well as a mathematical model of average RNFB trajectories as prior. The fully automatic process was applied to data recorded in healthy subjects and glaucoma patients and automatically detected individual RNFB trajectories are compared to manual traces.

Combining OCT and OCTA for Focal Structure-Function Modeling in Early Primary Open-Angle Glaucoma

Purpose

To investigate modeling of the focal visual field (VF) loss by combining structural measurements and vascular measurements in eyes with early primary open-angle glaucoma (POAG).

Methods

In this cross-sectional study, subjects with early glaucoma (VF mean deviation, ≥−6 dB) underwent optical coherence tomography (OCT) imaging, optical coherence tomography angiography (OCTA) imaging, and Humphrey 24-2 VF tests. Capillary perfusion densities (CPDs) were calculated after the removal of large vessels in the OCTA images. Focal associations between VF losses at the individual VF test locations, circumpapillary retinal nerve fiber layer (RNFL) thickness measurements from OCT, and CPDs were determined using nerve fiber trajectory tracings. Linear mixed models were used to model focal VF losses at each VF test location.

Results

Ninety-seven eyes with early POAG (VF mean deviation, −2.47 ± 1.64 dB) of 71 subjects were included. Focal VF modeling using a combined RNFL–CPD approach resulted in a median adjusted R² value of 0.30 (interquartile range [IQR], 0.13–0.55), whereas the RNFL-only and CPD-only approaches resulted in median values of 0.22 (IQR, 0.10–0.51) and 0.26 (IQR, 0.10–0.52), respectively. Seventeen VF locations with the combined approach had an adjusted R² value greater than 0.50. Likelihood testing at each VF test location showed that the combined approach performed significantly better at the superior nasal VF regions of the eyes compared with the univariate approaches.

Conclusions

Modeling of focal VF losses showed improvements when structural thickness and vascular parameters were included in tandem. Evaluation of VF defects in early glaucoma may benefit from considering both RNFL and OCTA characteristics.

Ganglion cell-inner plexiform layer measurements derived from widefield compared to montaged 9-field optical coherence tomography

CLINICAL RELEVANCE

With equivalent inner retinal thickness measurements compared to a more conventional composite optical coherence tomography (OCT) protocol, Widefield optical coherence tomography (WF-OCT) is a clinically viable, time-saving option facilitating detection of ocular pathologies within the central 55° of the retina.

PURPOSE

To compare ganglion cell-inner plexiform layer (GCIPL) thicknesses obtained using a single WF-OCT scan and standard composite OCT scans acquired in 9 fields of gaze (9F-OCT).

METHODS

Thirteen healthy participants underwent WF-OCT and 9F-OCT using the Spectralis OCT. The GCIPL was automatically segmented with a manual review for 9F-OCT and was manually segmented for WF-OCT. After registration, differences in GCIPL thicknesses were compared using 95% confidence intervals computed from one-sample t-tests and Bland-Altman analyses. Location-specific differences in B-scan tilt were analysed using Spearman correlations and linear regression models. To determine whether B-scan tilt influences GCIPL measurements, regression models of tilt versus differences between perpendicular and axial GCIPL thickness were applied.

RESULTS

While scattered locations demonstrated significant GCIPL thickness differences between WF-OCT and 9F-OCT, most differences did not exceed the axial pixel resolution of the instrument of 3.87 µm. Bland-Altman analyses indicated no notable bias using WF-OCT. Moderate correlations indicating significant location-specific differences in B-scan tilt were observed for temporal, central and inferior B-scans (r = -0.62 to 0.72), with linear regression models predicting a maximum difference in the tilt of 4.65°. The quadratic regression model indicated that at tilts greater than 27.3°, perpendicular GCIPL measurements become increasingly thin relative to axial measurements.

CONCLUSIONS

GCIPL thicknesses and B-scan tilts from WF-OCT are comparable to 9F-OCT, indicating that WF-OCT can be applied clinically to obtain valid inner retinal OCT measurements over 55° of the central retina with relative ease. However, for peripheral locations, B-scan tilt may need to be considered when measuring GCIPL thicknesses.

Developing a normative database for retinal perfusion using optical coherence tomography angiography

Visualizing and characterizing microvascular abnormalities with optical coherence tomography angiography (OCTA) has deepened our understanding of ocular diseases, such as glaucoma, diabetic retinopathy, and age-related macular degeneration. Two types of microvascular defects can be detected by OCTA: focal decrease because of localized absence and collapse of retinal capillaries, which is referred to as the non-perfusion area in OCTA, and diffuse perfusion decrease usually detected by comparing with healthy case-control groups. Wider OCTA allows for insights into peripheral retinal vascularity, but the heterogeneous perfusion distribution from the macula, parapapillary area to periphery hurdles the quantitative assessment. A normative database for OCTA could estimate how much individual's data deviate from the normal range, and where the deviations locate. Here, we acquired OCTA images using a swept-source OCT system and a 12×12 mm protocol in healthy subjects. We automatically segmented the large blood vessels with U-Net, corrected for anatomical factors such as the relative position of fovea and disc, and segmented the capillaries by a moving window scheme. A total of 195 eyes were included and divided into 4 age groups: < 30 (n=24) years old, 30-49 (n=28) years old, 50-69 (n=109) years old and >69 (n=34) years old. This provides an age-dependent normative database for characterizing retinal perfusion abnormalities in 12×12 mm OCTA images. The usefulness of the normative database was tested on two pathological groups: one with diabetic retinopathy; the other with glaucoma.

Improving Personalized Structure to Function Mapping From Optic Nerve Head to Visual Field

Purpose

Maps are required to relate visual field locations to optic nerve head regions. We compare individualized structure-to-function mapping (CUSTOM-MAP) to a population-derived mapping schema (POP-MAP).

Methods

Maps were compared for 118 eyes with glaucomatous field loss, circumpapillary retinal nerve fiber layer (cpRNFL) thickness measured using spectral domain optical coherence tomography (OCT), and two landmarks: the optic nerve head (ONH) position relative to the fovea and the temporal raphe angle. Locations with visual field damage (total deviation < −6 dB) were mapped to 30° ONH sectors centered on the angle given by each mapping schema. The concordance between damaged function and damaged structure was determined per location for various cpRNFL damage probability levels, with the number of concordant locations divided by the total number of damaged field locations providing a concordance ratio per eye.

Results

For the strictest concordance criteria (minimum cpRNFL thickness < 1% of normal), CUSTOM-MAP had higher mean concordance ratio than POP-MAP (60.5% c.f. 57.0% paired Wilcoxon, P = 0.005), with CUSTOM-MAP having a higher ratio in 43 eyes and POP-MAP having a higher ratio in 21 eyes. For all cpRNFL probability levels <20% of normal, more locations concorded for CUSTOM-MAP than POP-MAP. Inspection of the spatial patterns of differences revealed that CUSTOM-MAP often performed better in the arcuate regions, whereas POP-MAP had benefits inferior to the macula.

Conclusions

Anatomic parameters required for individualized structure-function mapping are readily measured with OCT and can provide improved concordance for some eyes.

Translational Relevance

Personalizing structure-function mapping may improve concordance between these measures. We provide a web-based tool for creating customized maps.

Association Between Structure-function Characteristics and Visual Field Outcomes in Glaucoma Subjects With Intraocular Pressure Reduction After Trabeculectomy

PRECIS

Improvements in post-trabeculectomy visual field (VF) outcomes were found to be significantly associated with preoperative nerve fiber layer thickness parameters extracted from the sectorized structure-function relationship, baseline VF, and severity of glaucoma.

OBJECTIVE

To determine whether the preoperative structure-function relationship helps to predict visual outcomes at 1-year post-trabeculectomy.

PATIENTS AND METHODS

In total, 91 eyes from 87 participants who successfully underwent trabeculectomy were included in our study. All eyes received optical coherence tomography imaging and VF assessment using 30-2 standard automated perimetry preoperatively at baseline and postoperatively 1 year after trabeculectomy. The linear mixed-model analysis was used to assess the association of structure and function at baseline, and multivariate analysis to investigate factors associated with postoperative VF outcomes.

RESULTS

Results from multivariate and univariate analysis for VF 1 year after trabeculectomy showed that a positive preoperative retinal nerve fiber layer thickness deviation from the structure-function model was found to be significantly associated with improved postoperative VF outcomes [β=0.06 dB/μm; 95% confidence interval (CI), 0.03-0.09]. Other significant factors included baseline VF MD (β=-0.18; 95% CI, -0.23 to -0.13) and the presence of severe glaucoma (β=-1.69; 95% CI, -2.80 to -0.57). Intraocular pressure was positively associated with improved VF outcomes only in univariate analysis (β=0.06; 95% CI, 0.01-0.11).

CONCLUSIONS AND RELEVANCE

Characteristics derived from the baseline structure-function relationship were found to be strongly associated with postoperative VF outcomes. These findings suggest that the structure-function relationship could potentially have a role in predicting VF progression after trabeculectomy.

Focal Structure-Function Relationships in Primary Open-Angle Glaucoma Using OCT and OCT-A Measurements

Purpose

To evaluate the focal structure-function associations among visual field (VF) loss, optical coherence tomography angiography (OCT-A) vascular measurements, and optical coherence tomography (OCT) structural measurements in glaucoma.

Methods

In this cross-sectional study, subjects underwent standard automated perimetry, OCT-based nerve fiber thickness measurements, and OCT-A imaging. Mappings of focal VF test locations with OCT and OCT-A measurements were defined using anatomically adjusted nerve fiber trajectories and were studied using multivariate mixed-effects analysis. Segmented regression analysis was used to determine the presence of breakpoints in the structure-function associations.

Results

The study included 119 eyes from 86 Chinese subjects with primary open-angle glaucoma (POAG). VF mean deviation was significantly associated with global capillary perfusion density (β = 0.13 ± 0.08) and global retinal nerve fiber layer thickness (β = 0.09 ± 0.02). Focal capillary density (FCD) was significantly associated with VF losses at 34 VF test locations (66.7% of 24-2 VF), with 24 of the 34 locations being within 20° of retinal eccentricity. Focal nerve layer (FNL) thickness was significantly associated with 16 VF test locations (31.4% of 24-2 VF; eight locations within 20° eccentricity). For VF test locations in the central 10° VF, VF losses below the breakpoint were significantly associated with FCD (slope, 0.89 ± 0.12, P < 0.001), but not with FNL thickness (slope, 0.57 ± 0.39, P = 0.15).

Conclusions

Focal capillary densities were significantly associated with a wider range of visual field losses and in a larger proportion of the visual field compared to nerve fiber thickness.

Anatomical Location of the Raphe and Extended Raphe in the Human Retina: Implications for Assessment of the Optic Nerve with OCT

Purpose

To determine the location of (1) the superior–inferior watershed between the fovea and optic disc (extended raphe) at the peripapillary optical coherence tomography (OCT) measurement circle and (2) the raphe, temporal to the fovea.

Methods

We used existing data consisting of 2285 traced retinal nerve fiber bundle trajectories from 83 fundus images. For localization of the extended raphe at the 3.46-mm-diameter OCT measurement circle, trajectories were classified as belonging to the superior or inferior hemiretina, using predefined criteria. For the raphe, we localized the endings of trajectories coming from the superior and inferior arcuate bundles.

Results

At the measurement circle, the extended raphe is located 14° (range, 12°–16°) inferiorly to a horizontal line through the optic disc center. The raphe follows a horizontal line at the latitude of the fovea if the disc is assumed to be located 15° nasal to and 2° above the fovea.

Conclusions

At the measurement circle, OCT brands use either the 9 o'clock location or a straight line connecting the center of the optic disc and the fovea as a reference for separating the hemiretinas. This results, on average, in a 14° and 6° misalignment with respect to the anatomical watershed, respectively. For the macular area, the commonly used line through the center of the optic disc and the fovea fails to describe the raphe adequately.

Translational Relevance

An unbiased asymmetry assessment of the optic nerve requires a detailed knowledge of the shape and location of the (extended) raphe.

Artificial Intelligence Mapping of Structure to Function in Glaucoma

Purpose

To develop an artificial intelligence (AI)-based structure-function (SF) map relating retinal nerve fiber layer (RNFL) damage on spectral domain optical coherence tomography (SDOCT) to functional loss on standard automated perimetry (SAP).

Methods

The study included 26,499 pairs of SAP and SDOCT from 15,173 eyes of 8878 patients with glaucoma or suspected of having the disease extracted from the Duke Glaucoma Registry. The data set was randomly divided at the patient level in training and test sets. A convolutional neural network (CNN) was initially trained and validated to predict the 52 sensitivity threshold points of the 24-2 SAP from the 768 RNFL thickness points of the SDOCT peripapillary scan. Simulated localized RNFL defects of varied locations and depths were created by modifying the normal average peripapillary RNFL profile. The simulated profiles were then fed to the previously trained CNN, and the topographic SF relationships between structural defects and SAP functional losses were investigated.

Results

The CNN predictions had an average correlation coefficient of 0.60 (P < 0.001) with the measured values from SAP and a mean absolute error of 4.25 dB. Simulated RNFL defects led to well-defined arcuate or paracentral visual field losses in the opposite hemifield, which varied according to the location and depth of the simulations.

Conclusions

A CNN was capable of predicting SAP sensitivity thresholds from SDOCT RNFL thickness measurements and generate an SF map from simulated defects.

Translational Relevance

AI-based SF map improves the understanding of how SDOCT losses translate into detectable SAP damage.

An Automated Method for Assessing Topographical Structure-Function Agreement in Abnormal Glaucomatous Regions

Purpose

To develop an automated/objective method for topographically comparing abnormal regions on optical coherence tomography (OCT) and visual field (VF) tests of eyes with early glaucoma.

Methods

A custom R program was developed that allows for both visualization and automatic assessment of the topographical agreement between functional (24-2 and/or 10-2 VF) and structural (widefield OCT retinal nerve fiber layer and/or retinal ganglion cell layer) deviation/probability maps. It was optimized using information from 98 eyes: 53 diagnosed as "definitely glaucoma" (DG) and 45 recruited as healthy (H) controls. Different pairs of abnormal VF (P <1%, <2%, <5%) and abnormal OCT (P <5%, <10%, <15%) criteria were evaluated. The percentages of abnormal structure-abnormal function (aS-aF) agreement found in DG eyes and nonagreement found in H eyes were used to define the optimal criteria and number of aS-aF locations for the detection of aS-aF agreement.

Results

A criterion of two aS-aF locations with "OCT <10% and VF <5%" on VF pattern deviation (PD) probability and OCT deviation/probability maps yielded high overall agreement (92%) with high aS-aF agreement for the DG eyes (89%) and high aS-aF nonagreement for the H eyes (95%). Total deviation probability maps achieved slightly lower performance than PD maps.

Conclusions

The method described here can automatically and objectively evaluate aS-aF agreement with a direct comparison of abnormal regions of function and structure.

Translational Relevance

As glaucoma diagnosis often involves assessing structure-function agreement, this technique can overcome subjectivity in this assessment.

[Structure and function of the visual pathway]

Humans receive information from their environment mainly via the visual system. Signals from the photoreceptors of the retina via bipolar and ganglion cells are projected onto specific neuronal subpopulations in the lateral geniculate body and from there are forwarded to appropriate layers of the primary visual cortex. The most important anatomical and functional features of the visual system are explained. For this purpose, a selective literature search was carried out in the databases PubMed (also in Europe PubMed Central), Psychline, Google Scholar, Cochrane Library and Web of Science as well as additional information in relevant books or websites in the fields of (neuro)anatomy, (neuro)physiology, (neuro)ophthalmology and (neuro)otology, among others with the search terms Sehbahn, visual system, visual pathway, receptors, spatial cognition and visual cognition.

Structure-Function Agreement Is Better Than Commonly Thought in Eyes With Early Glaucoma

Purpose

To assess the agreement between structural (optical coherence tomography [OCT]) and functional (visual field [VF]) glaucomatous damage with an automated method and deviation/probability maps, and to compare this method to a metric method.

Methods

Wide-field spectral-domain OCT scans, including the disc and macula, and 24-2 and 10-2 VFs were obtained from 45 healthy control (H) eyes/individuals, and 53 eyes/patients with 24-2 mean deviation (MD) better than -6 dB diagnosed as "definite glaucoma" (DG) by experts. Abnormal structure-abnormal function (aS-aF) agreement was assessed with an automated topographic (T) method based upon VF pattern deviation and OCT probability maps. Results were compared to a metric (M) method optimized for accuracy, (abnormal 24-2 glaucoma hemifield test [GHT] or pattern standard deviation [PSD], or 10-2 PSD AND abnormal OCT [quadrant]).

Results

For the T-method, 47 (88.7%) of the 53 DG eyes showed aS-aF agreement, compared to 2 (4.5%) of the 45 H eyes. The aS-aF agreement for these two H eyes was easily identified as mistaken, and did not replicate on a subsequent test. Without the 10-2, the aS-aF agreement decreased from 47 to 34 (64.2%) of 53 DG eyes. For the M-method, 37 (69.8%) of the 53 DG eyes showed aS-aF agreement, while omitting the 10-2 VF resulted in agreement in only 33 (62.3%) eyes.

Conclusions

There is good agreement between structural and functional damage, even in eyes with confirmed early glaucomatous damage, if both 24-2 and 10-2 VFs are obtained, and abnormal locations on the VFs are compared to abnormal regions seen on OCT macular and disc scans. This can be done in an objective, automated fashion. (ClinicalTrials.gov number, NCT02547740.).

OCT-Detected Optic Nerve Head Neural Canal Direction, Obliqueness, and Minimum Cross-Sectional Area in Healthy Eyes

PURPOSE

To assess anterior scleral canal opening (ASCO) offset relative to Bruch's membrane opening (BMO) (ASCO/BMO offset) so as to determine neural canal direction, obliqueness, and minimum cross-sectional area (NCMCA) in 362 healthy eyes.

DESIGN

Cross-sectional study.

METHODS

After optical coherence tomography optic nerve head and retinal nerve fiber layer thickness (RNFLT) imaging, BMO and ASCO were manually segmented. Planes, centroids, size, and shape were calculated. Neural canal direction was defined by projecting the neural canal axis vector (connecting BMO and ASCO centroids) onto the BMO plane. Neural canal obliqueness was defined by the angle between the neural canal axis and the BMO plane perpendicular vector. NCMCA was defined by projecting BMO and ASCO points onto a neural canal axis perpendicular plane and measuring the area of overlap. The angular distance between superior and inferior peak RNFLT was measured, and correlations between RFNLT, BMO, ASCO, ASCO/BMO offset, and NCMCA were assessed.

RESULTS

Mean (SD) NCMCA was significantly smaller than either the BMO or ASCO area (1.33 (0.42), 1.82 (0.38), 2.22 (0.43) mm2, respectively), and most closely correlated to RNFLT (P < .001, R2 = 0.158). Neural canal direction was most commonly superior-nasal (55%). Mean neural canal obliqueness was 39.4° (17.3°). The angular distance between superior and inferior peak RNFLT correlated to neural canal direction (P ≤ .008, R2 = 0.093).

CONCLUSIONS

ASCO/BMO offset underlies neural canal direction, obliqueness, and NCMCA. RNFLT is more strongly correlated to NCMCA than to BMO or ASCO, and its peripapillary distribution is influenced by neural canal direction.

Development of a Spatial Model of Age-Related Change in the Macular Ganglion Cell Layer to Predict Function From Structural Changes

PURPOSE

To develop location-specific models of normal, age-related changes in the macular ganglion cell layer (GCL) from optical coherence tomography (OCT). Using these OCT-derived models, we predicted visual field (VF) sensitivities and compared these results to actual VF sensitivities.

DESIGN

Retrospective cohort study.

METHODS

Single eyes of 254 normal participants were retrospectively enrolled from the Centre for Eye Health (Sydney, Australia). Macular GCL measurements were obtained using Spectralis OCT. Cluster algorithms were performed to identify spatial patterns demonstrating similar age-related change. Quadratic and linear regression models were subsequently used to characterize age-related GCL decline. Forty participants underwent additional testing with Humphrey VFs, and 95% prediction intervals were calculated to measure the predictive ability of structure-function models incorporating cluster-based pooling, age correction, and consideration of spatial summation.

RESULTS

Quadratic GCL regression models provided a superior fit (P value <.0001-.0066), establishing that GCL decline commences in the late 30s across the macula. The equivalent linear rates of GCL decline showed eccentricity-dependent variation (0.13 μm/yr centrally vs 0.06 μm/yr peripherally); however, average, normalized GCL loss per year was consistent across the 64 macular measurement locations at 0.26%. The 95% prediction intervals describing predicted VF sensitivities were significantly narrower across all cluster-based structure-function models (3.79-4.99 dB) compared with models without clustering applied (5.66-6.73 dB, P < .0001).

CONCLUSIONS

Combining spatial clustering with age-correction based on regression models allowed the development of robust models describing GCL changes with age. The resultant superior predictive ability of VF sensitivity from ganglion cell measurements may be applied to future models of disease development to improve detection of early macular GCL pathology.

Relating optical coherence tomography to visual fields in glaucoma: structure-function mapping, limitations and future applications

Combining information from optical coherence tomography (OCT) imaging and visual field testing is useful in the clinical assessment and monitoring of patients with glaucoma. Measurements of retinal nerve fibre layer thickness or neuroretinal rim width taken around the optic nerve head may be related to the visual field using a structure-function map. In this review, the structure-function mapping methods in clinical use are discussed. Typical clinical maps provide a population average, 'one size fits all' representation, but in recent years methods for customising structure-function maps to individual eyes have been developed and these are reviewed here. In the macula, visual field stimuli stimulate photoreceptors for which associated retinal ganglion cells are peripherally displaced. Recently developed methods that relate OCT measurements to visual field test locations in the macula are therefore also reviewed. The use of structure-function maps to relate OCT measurements to localised visual field sensitivity in new applications is also explored. These new applications include the selection of visual field test locations and stimulus intensities based on OCT data, and the formal post-test combination of results across modalities. Such applications promise to exploit the structure-function relationship in glaucoma to improve disease diagnosis and monitoring of progression. Limitations in the validation and use of current structure-function mapping techniques are discussed.

Identifying Glaucomatous Damage to the Macula

SIGNIFICANCE

Measurements of the macula have been increasingly used to diagnose and manage patients with glaucoma. Asymmetry analysis was clinically introduced to assess damage to the macular ganglion cells in patients with glaucoma, but its effectiveness is limited by high normal between-subject variability.

PURPOSE

We aimed to reduce the high normal between-subject variability and improve the potential of asymmetry analysis to identify glaucomatous damage to the macula.

METHODS

Twenty patients with glaucoma (aged 57 to 85 years) and 30 age-similar control subjects (aged 53 to 89 years) were recruited from a longitudinal glaucoma study. Participants were imaged with the Spectralis OCT using the posterior pole protocol; measurements of the averaged retinal thickness and ganglion cell layer (GCL) thickness were obtained. We established three zones per hemifield within the central ±9°, based on the lowest between-subject variability that we previously found and the course of retinal nerve fiber layer projections. The criteria for flagging abnormality were at least two contiguous zones when P < 5% or one zone when P < 1% with two-tailed tests.

RESULTS

Between-subject variability of the asymmetry analysis for both retinal and GCL thicknesses remained lower than that of the average thickness across each zone in control subjects (F > 2.52, P < .01). Asymmetry analysis of retinal and GCL thicknesses flagged 16 and 18 of 20 patients, respectively.

CONCLUSIONS

Between-subject variability was reduced in control subjects using the three zones; our criteria identified glaucomatous damage to the macula in most of the patients. We used high-density B-scans to confirm the patterns of the glaucomatous damage we found in this study.

Application of Pattern Recognition Analysis to Optimize Hemifield Asymmetry Patterns for Early Detection of Glaucoma

Purpose

To assess the diagnostic utility of a new hemifield asymmetry analysis derived using pattern recognition contrast sensitivity isocontours (CSIs) within the Humphrey Field Analyzer (HFA) 24-2 visual field (VF) test grid. The performance of an optimal CSI-derived map was compared against a commercially available clustering method (Glaucoma Hemifield Test, GHT).

Methods

Five hundred VF results of 116 healthy subjects were used to determine normative distribution limits for comparisons. Pattern recognition analysis was applied to HFA 24-2 sensitivity data to determine CSI theme maps delineating clusters for hemifield comparisons. Then, 1019 VF results from 228 glaucoma patients were assessed using different clustering methods to determine the true-positive rate. We also assessed additional 354 VF results of 145 healthy subjects to determine the false-positive rate.

Results

The optimum clustering method was the CSI-derived seven-theme class map, which identified more glaucomatous VFs compared with the GHT map. The seven-class theme map also identified more cases compared with the five-, six-, and eight-class maps, suggesting no effect of number of clusters. Integrating information regarding the location of glaucomatous defects to the CSI clusters did not improve detection rate.

Conclusions

A clustering map derived using CSIs improved detection of glaucomatous VFs compared with the currently available GHT. An optimized CSI-derived map may serve as an additional means to aid earlier detection of glaucoma.

Translational Relevance

Pattern recognition–derived theme maps provide a means for guiding test point selection for asymmetry analysis in glaucoma assessment.

Retinal nerve fiber bundle trajectories in Chinese myopic eyes: Comparison with a Caucasian based mathematical model

Previously we developed a mathematical model for describing the retinal nerve fiber bundle (RNFB) trajectories in the human retina. The model was based on Caucasian eyes that were not selected regarding refraction. The aim of this study was to determine the characteristics of the RNFB trajectories in Chinese myopic eyes. We collected high quality red free fundus images from 80 eyes of 80 Chinese myopic subjects (median [interquartile range/range] refraction -3.9 [-6.0 to -2.5/-10 to -1] D). We traced all visible RNFBs (n = 1460) and evaluated their trajectories using the previously published mathematical model. In the superior-temporal region, the RNFB trajectories of the Chinese myopic eyes were similar to that of the Caucasian eyes (86% of trajectories within the 95% central range of the Caucasian model). In the inferior-temporal region, the trajectories of the Chinese low to moderate myopic eyes were also similar to that of the Caucasian eyes (85%); trajectories of the high myopic eyes (spherical equivalent beyond -6.00 D) were clearly less curved (75%). Associations between individual deviations from the model and axial length, retinal vessel course, and optic disc anatomy were studied with multiple linear regression analysis. In the superior-temporal region, the trajectories were associated with retinal vessel course (P = 0.008) and optic disc size (P = 0.016). In the inferior-temporal region, there was a significant association with axial length (P < 0.001), retinal vessel course (P = 0.006), and disc torsion (P = 0.009).

Calculation and plotting of retinal nerve fiber paths based on Jansonius et al. 2009/2012 with an R program

The data presented in this article are related to the research article entitled “Retinal conduction speed analysis reveals different origins of the P50 and N95 components of the (multifocal) pattern electroretinogram” (Bach et al., 2018) [1]. That analysis required the individual length data of the retinal nerve fibers (from ganglion cell body to optic nerve head, depending on the position of the ganglion cell body). Jansonius et al. (2009, 2012) [2,3] mathematically modeled the path morphology of the human retinal nerve fibers. We here present a working implementation with source code (for the free and open-source programming environment “R”) of the Jansonius’ formulas, including all errata. One file defines Jansonius et al.’s “phi” function. This function allows quantitative modelling of paths (and any measures derived from them) of the retinal nerve fibers. As a working demonstration, a second file contains a graph which plots samples of nerve fibers. The included R code runs in base R without the need of any additional packages.

Customizing Perimetric Locations Based on En Face Images of Retinal Nerve Fiber Bundles With Glaucomatous Damage

Purpose

Prior studies suggested the use of customized perimetric locations in glaucoma; these studies were limited by imaging only the superficial depths of the retinal nerve fiber layer (RNFL) and by prolonged perimetric testing. We aimed to develop a rapid perimetric test guided by high-resolution images of RNFL bundles.

Methods

We recruited 10 patients with glaucoma, ages 56 to 80 years, median 68 years, and 10 controls, ages 55 to 77 years, median 68 years. The patients were selected based on discrepancies between locations of glaucomatous damage for perimetric and structural measures. Montaging was used to produce optical coherence tomography en face images of the RNFL covering much of the 24-2 grid locations. In experiment 1, we presented the Goldmann size III stimulus at preselected retinal locations of glaucomatous damage, using just two contrasts. In experiment 2, we developed an elongated sinusoidal stimulus, aligned within the defect, to measure contrast sensitivities; abnormalities were defined based on lower 95% reference limits derived from the controls.

Results

The percentage of predicted locations where size III was not seen at 28 dB ranged from 16% to 80%, with a median of 48%. Contrast sensitivity for the sinusoidal stimulus was below the 95% reference range for 37 of 44 stimuli aligned within the defects.

Conclusions

We developed methods for rapid perimetric testing guided by en face images of the RNFL bundles in patients with glaucoma. Results indicated ganglion cell damage under all of the visible RNFL defects.

Translational Relevance

Customized perimetric locations have potential to improve clinical assessment of glaucoma.

[Structural endpoints for glaucoma studies]

BACKGROUND

Structural endpoints have been discussed as surrogate endpoints for the approval of neuroprotective drugs in glaucoma.

OBJECTIVE

Is the evidence strong enough to establish structural endpoints as surrogate endpoints?

MATERIAL AND METHODS

Review of current understanding between structure and function in glaucoma.

RESULTS

The introduction of optical coherence tomography has revolutionized imaging in glaucoma patients. Clinically either the nerve fiber layer thickness can be measured along a circle centered in the optic nerve head or the ganglion cell layer thickness can be assessed in the macular region, the latter being quantified in combination with other inner retinal layers. On a microscopic level there is a strong correlation between structural and functional loss but this relation can only partially be described with currently available clinical methods. This is particularly true for longitudinal course of the disease in glaucoma patients. Novel imaging techniques that are not yet used clinically may have the potential to increase our understanding between structure and function in glaucoma but further research in this field is required.

CONCLUSION

The current evidence does not allow the establishment of structural endpoints as surrogate endpoints for phase 3 studies in glaucoma. Neuroprotective drugs have to be approved on the basis of visual field data because this is the patient-relevant endpoint. Structural endpoints can, however, play an important role in phase 2 and proof of concept studies.

Nerve Fiber Flux Analysis Using Wide-Field Swept-Source Optical Coherence Tomography

Purpose

To devise a method to quantify nerve fibers over their arcuate courses over an extended peripapillary area using optical coherence tomography (OCT).

Methods

Participants were imaged with 8 × 8-mm volumetric OCT scans centered at the optic disc. A new quantity, nerve fiber flux (NFF), represents the cross-sectional area transected perpendicular to the nerve fibers. The peripapillary area was divided into 64 tracks with equal flux. An iterative algorithm traced the trajectory of the tracks assuming that the relative distribution of the NFF was conserved with compensation for fiber connections to ganglion cells on the macular side. Average trajectory was averaged from normal eyes and use to calculate the NFF maps for glaucomatous eyes. The NFF maps were divided into eight sectors that correspond to visual field regions.

Results

There were 24 healthy and 10 glaucomatous eyes enrolled. The algorithm converged on similar patterns of NFL tracks for all healthy eyes. In glaucomatous eyes, NFF correlated with visual field sensitivity in the arcuate sectors (Spearman ρ = 0.53–0.62). Focal nerve fiber loss in glaucomatous eyes appeared as uniform tracks of NFF defects that followed the expected arcuate fiber trajectory.

Conclusions

Using an algorithm based on the conservation of flux, we derived nerve fiber trajectories in the peripapillary area. The NFF map is useful for the visualization of focal defects and quantification of sector nerve fiber loss from wide-area volumetric OCT scans.

Translational Relevance

NFF provides a cumulative measure of volumetric loss along nerve fiber tracks and could improve the detection of focal glaucoma damage.

Deep Defects Seen on Visual Fields Spatially Correspond Well to Loss of Retinal Nerve Fiber Layer Seen on Circumpapillary OCT Scans

Purpose

To examine the structure-function relationship in glaucoma between deep defects on visual fields (VF) and deep losses in the circumpapillary retinal nerve fiber layer (cpRNFL) on optical coherence tomography (OCT) circle scans.

Methods

Thirty two glaucomatous eyes with deep VF defects, as defined by at least one test location worse than ≤ -15 dB on the 10-2 and/or 24-2 VF pattern deviation (PD) plots, were included from 87 eyes with "early" glaucoma (i.e., 24-2 mean deviation better than -6 dB). Using the location of the deep VF points and a schematic model, the location of local damage on an OCT circle scan was predicted. The thinnest location of cpRNFL (i.e., deepest loss) was also determined.

Results

In 19 of 32 eyes, a region of complete or near complete cpRNFL loss was observed. All 19 of these had deep VF defects on the 24-2 and/or 10-2. All of the 32 eyes with deep VF defects had abnormal cpRNFL regions (red, 1%) and all but 2 had a region of cpRNFL thickness <21 μm. The midpoint of the VF defect and the location of deepest cpRNFL had a 95% limit of agreement within approximately two-thirds of a clock-hour (or 30°) sector (between -22.1° to 25.2°). Individual fovea-to-disc angle (FtoDa) adjustment improved agreement in one eye with an extreme FtoDa.

Conclusions

Although studies relating local structural (OCT) and functional (VF) measures typically show poor to moderate correlations, there is good qualitative agreement between the location of deep cpRNFL loss and deep defects on VFs.

Retinal conduction speed analysis reveals different origins of the P50 and N95 components of the (multifocal) pattern electroretinogram

The pattern electroretinogram (PERG), an indicator of retinal ganglion cell (RGC) function, comprises a P50 and an N95 component. We addressed the question of whether the N95 originates, like the P50, from the RGC bodies or from the change of axon orientation at the optic nerve head (ONH). Thus, we recorded multifocal PERGs for 36 retinal locations in 21 participants. Second-order kernel responses were analyzed for the dependence of peak time topography on retinal fiber lengths to the ONH separately for the positive and negative excursions. We found that peak times were longer for macular [P1 (P50-like): 50 ms; N2 (N95-like): 76)] than for peripheral responses [P1: 43; N2: 66]. For the N2 another factor was necessary to explain the variability: The time difference (deltaT: N2 minus P1) was found to be proportional to fiber length from ganglion cell body to the ONH. We calculated retinal fiber length using an analytical function by Jansonius et al. (2009, 2012) and found that a linear model with factors eccentricity and fiber length explained 82% of the total N2 time variance (p«0.001). The conduction speeds of the retinal axons were estimated from deltaT to range from 0.5 to 3.0 m/s for parafovea and periphery, respectively. The dependence of deltaT on the distance from ganglion cell body to the ONH suggests that the N2 originates at the ONH rather than at the ganglion cell body. While the multifocal N2 peaks earlier (≈76 ms) than the non-multifocal PERG-N95 (≈95 ms), considerations of high-pass filtering and frequency dependence of the mfPERG-N2 suggest that the source separation (P50 = ganglion cell body vs. N95 = ONH) also holds for the non-multifocal PERG.

The ISNT Rule: How Often Does It Apply to Disc Photographs and Retinal Nerve Fiber Layer Measurements in the Normal Population?

PURPOSE

To determine what percentage of normal eyes follow the ISNT rule, and whether ISNT rule variants may be more generalizable to the normal population.

DESIGN

Cross-sectional study.

METHODS

Setting: Institutional setting.

STUDY POPULATION

Total of 110 normal subjects.

OBSERVATION PROCEDURES

Neuroretinal rim assessments from disc photographs and retinal nerve fiber layer (RNFL) thickness measurements from spectral-domain optical coherence tomography.

MAIN OUTCOME MEASURES

The percentages of subjects that obeyed the ISNT rule and its variants.

RESULTS

The ISNT rule is only valid for 37.0% of disc photograph rim assessments and 43.8% of RNFL measurements. Deviation of the nasal sector from the expected ISNT pattern was a major cause for the ISNT rule not being obeyed for both rim and RNFL assessments. Specifically, 10.9% of subjects had wider nasal rims than the inferior rims, 29.4% had wider nasal rims than the superior rims, 14.7% had narrower nasal rims than the temporal rims, and 42.9% had thinner nasal RNFLs compared to the temporal quadrant. Exclusion of the nasal quadrant from the ISNT rule significantly increased the validity of ISNT variant rules, with 70.9% and 76.4% of disc photographs following the IST rule and the IS rule, respectively. Similarly, for RNFL thickness, 70.9% and 71.8% of patients followed the IST and IS rule, respectively.

CONCLUSIONS

The ISNT rule is only valid for about a third of disc photographs and less than half of RNFL measurements in normal patients. ISNT rule variants, such as the IST and IS rule, may be considered, as they are valid in more than 70% of patients.

Glaucoma Diagnostic Capability of Global and Regional Measurements of Isolated Ganglion Cell Layer and Inner Plexiform Layer

PURPOSE

To compare glaucoma diagnostic capability of global/regional macular layer parameters in different-sized grids.

MATERIALS AND METHODS

Serial horizontal spectral-domain optical coherence tomography scans of macula were obtained. Automated macular grids with diameters of 3, 3.45, and 6 mm were used. For each grid, 10 parameters (total volume; average thicknesses in 9 regions) were obtained for 5 layers: macular retinal nerve fiber layer (mRNFL), ganglion cell layer (GCL), inner plexiform layer (IPL), ganglion cell-inner plexiform layer (GCIPL; GCL+IPL), and ganglion cell complex (GCC; mRNFL+GCL+IPL).

RESULTS

Sixty-nine normal eyes (69 subjects) and 87 glaucomatous eyes (87 patients) were included. For the total volume parameter, the area under the receiver operating characteristic curves (AUCs) in 6-mm grid were larger than the AUCs in 3- and 3.45-mm grids for GCL, GCC, GCIPL, and mRNFL (all P<0.020). For the average thickness parameters, the best AUC in 6-mm grid (T2 region for GCL, IPL, and GCIPL; I2 region for mRNFL and GCC) was greater than the best AUC in 3-mm grid for GCL, GCC, and mRNFL (P<0.045). The AUC of GCL volume (0.920) was similar to those of GCC (0.920) and GCIPL (0.909) volume. The AUC of GCL T2 region thickness (0.942) was similar to those of GCC I2 region (0.942) and GCIPL T2 region (0.934) thickness.

CONCLUSIONS

Isolated macular GCL appears to be as good as GCC and GCIPL in glaucoma diagnosis, while IPL does not. Larger macular grids may be better at detecting glaucoma. Each layer has a characteristic region with the best glaucoma diagnostic capability.

Structure and Function Relationship of Activated Retinal Glia in Primary Open-Angle Glaucoma Patients

Purpose

To evaluate clinically activated retinal astrocytes and Müller cells (ARAM) regarding retinal sensitivity and retinal nerve fiber layer (RNFL) thickness in primary open-angle glaucoma (POAG).

Methods

Central visual field (VF; i.e., retinal sensitivity) was measured with a custom-made macular pattern by microperimetry and correlated with the presence (ARAM+) or absence (ARAM−) of ARAM on red-free fundus photography and with the corresponding RNFL by optical coherence tomography (OCT).

Results

In the eyes of POAG patients, ARAM+ had overall a significantly lower retinal sensitivity (ARAM+: 7.34 dB, ARAM−: 11.9 dB; p < 0.001) and lower RNFL thickness in the inferior peripapillary quadrants compared to ARAM− (RNFL superior: ARAM+ 74.2 μm, ARAM− 77.5 μm; RNFL temporal: ARAM+ 46.8 μm, ARAM− 53.0 μm, p < 0.001; and RNFL inferior: ARAM+ 63.2 μm, ARAM− 73.1 μm, p < 0.001). Within the same eye, ARAM+ showed a lower retinal sensitivity compared to ARAM− ([ARAM− (11.13 dB)] − [ARAM+ (9.56 dB) = 1.57 dB; p = 0.25). The proportion of ARAM+ per eye correlated strongly with reduced retinal light sensitivity (p = 0.02), corresponding lower peripapillary RNFL thickness (p = 0.02), and lower RNFL temporal quadrant thickness (p < 0.01), but not with greater age (p = 0.45).

Conclusion

ARAM was more frequently identified in the eyes with a lower retinal sensitivity and peripapillary RNFL thickness and may be a clinical sign in the macula for an advanced stage of POAG.

Improving our understanding, and detection, of glaucomatous damage: An approach based upon optical coherence tomography (OCT)

Although ophthalmologists are becoming increasingly reliant upon optical coherence tomography (OCT), clinicians who care for glaucoma patients are not taking full advantage of the potential of this powerful technology. First, we ask, how would one describe the nature of glaucomatous damage if only OCT scans were available? In particular, a schematic model of glaucomatous damage is developed in section 2, and the nature of glaucomatous damage seen on OCT scans described in the context of this model in section 3. In particular, we illustrate that local thinning of the circumpapillary retinal nerve fiber layer (cpRNFL) around the optic disc can vary in location, depth, and/or width, as well as homogeneity of damage. Second, we seek to better understand the relationship between the thinning of the cpRNFL and the various patterns of sensitivity loss seen on visual fields obtained with standard automated perimetry. In sections 4 and 5, we illustrate why one should expect a wide range of visual field patterns, and iilustrate why they should not be placed into discrete categories. Finally, section 6 describes how the clinician can take better advantage of the information in OCT scans. The approach is summarized in a single-page report, which can be generated from a single wide-field scan. The superiority of this approach, as opposed to the typical reliance on summary metrics, is described.

Between-Subject Variability in Healthy Eyes as a Primary Source of Structural-Functional Discordance in Patients With Glaucoma

Purpose

To test with an independent data set the finding that between-subject variability in healthy eyes is the primary source of structural–functional discordance in patients with glaucoma.

Methods

Neuroretinal rim area, retinal nerve fiber layer thickness, and perimetric data were analyzed for one eye in each of 55 control subjects and for 245 right eyes of patients in the United Kingdom Glaucoma Treatment Study. Data were gathered with the Heidelberg Retina Tomograph (HRT), Stratus Optical Coherence Tomograph (OCT), and Humphrey Field Analyzer (HFA). Discordance was quantified as width of the limits of agreement from a Bland-Altman analysis of depth of defect. The ratio of variances (F test) for the patient and control groups was computed for comparisons of HFA-OCT, HFA-HRT, and OCT-HRT. Bonferroni adjustment required P less than 0.017 for statistical significance. The discordance in the patients was also quantified as the 95% prediction interval computed from the discordance in controls using the Hood-Kardon model for the HFA-OCT comparison.

Results

The F ratio comparing discordance in patients and controls was 0.77, 1.43, and 1.32 for the HFA-OCT, HFA-HRT, and OCT-HRT comparisons with P values 0.88, 0.06, and 0.11, respectively. For the Hood-Kardon model, 4.7% of the patients had discordance outside the 95% prediction interval computed from the discordance in controls. Similar results were obtained when all comparisons were repeated for left eyes of patients.

Conclusions

These results confirm previous findings that between-subject variability in healthy eyes is the primary source of structural–functional discordance in patients with glaucoma, and extends this finding to a structural–structural comparison.

Incorporating Spatial Models in Visual Field Test Procedures

Purpose

To introduce a perimetric algorithm (Spatially Weighted Likelihoods in Zippy Estimation by Sequential Testing [ZEST] [SWeLZ]) that uses spatial information on every presentation to alter visual field (VF) estimates, to reduce test times without affecting output precision and accuracy.

Methods

SWeLZ is a maximum likelihood Bayesian procedure, which updates probability mass functions at VF locations using a spatial model. Spatial models were created from empirical data, computational models, nearest neighbor, random relationships, and interconnecting all locations. SWeLZ was compared to an implementation of the ZEST algorithm for perimetry using computer simulations on 163 glaucomatous and 233 normal VFs (Humphrey Field Analyzer 24-2). Output measures included number of presentations and visual sensitivity estimates.

Results

There was no significant difference in accuracy or precision of SWeLZ for the different spatial models relative to ZEST, either when collated across whole fields or when split by input sensitivity. Inspection of VF maps showed that SWeLZ was able to detect localized VF loss. SWeLZ was faster than ZEST for normal VFs: median number of presentations reduced by 20% to 38%. The number of presentations was equivalent for SWeLZ and ZEST when simulated on glaucomatous VFs.

Conclusions

SWeLZ has the potential to reduce VF test times in people with normal VFs, without detriment to output precision and accuracy in glaucomatous VFs.

Translational Relevance

SWeLZ is a novel perimetric algorithm. Simulations show that SWeLZ can reduce the number of test presentations for people with normal VFs. Since many patients have normal fields, this has the potential for significant time savings in clinical settings.