A novel approach to event-based progression analysis in glaucomatous visual fields

doi:10.21203/rs.3.rs-446625/v1

PDF

Research Article

A novel approach to event-based progression analysis in glaucomatous visual fields

https://doi.org/10.21203/rs.3.rs-446625/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 11 Aug, 2021

Read the published version in Scientific Reports →

You are reading this latest preprint version

Glaucoma is a progressive optic neuropathy determining characteristic changes to the optic nerve head and the visual field (VF). Detecting progression of VF damage with Standard Automated Perimetry (SAP) is of paramount importance for clinical care. One common approach to detecting progression is to compare each new VF test to a baseline SAP test (event analysis). This comparison is made difficult by the test-retest variability of SAP, which increases with the level of VF damage, and the lower bounds on the range of measurement, meaning that damage cannot be assessed below a certain level.

We performed a prospective international multi-centre data collection of SAP data on 90 eyes from 90 people with glaucoma and different levels of VF damage over a short period of time (6 tests in 60 days). Data were collected using a fundus tracked perimeter (Compass, CenterVue). We used these data (minus the first test) to develop an improved event analysis that accounts for both the change in variability with damage and the lower bound on the measurement imposed by SAP. Using simulations, we show that our approach is more sensitive compared to previously developed methods, especially in the case of advanced glaucoma, while retaining similar specificity.

Ophthalmology

Cognitive Neuroscience

glaucoma

Standard Automated Perimetry (SAP)

visual field (VF)

Glaucoma is a progressive optic neuropathy characterised by typical morphological changes to the optic nerve head (ONH) and the retina. Functional damage from glaucoma manifests as characteristic defects to the visual field (VF). Monitoring the progression of VF damage is of paramount importance for clinical management and the main outcome of many glaucoma trials^1–8. The VF is usually measured with Static Automated Perimetry (SAP). In SAP, circular stimuli of different intensities are projected at various predetermined locations on the retina. A recent development of SAP, called Fundus tracked Automated Perimetry (FAP), uses continuous infrared imaging of the retina to track and compensate for eye movements and has been widely used to test the central VF in patients with age-related macular degeneration⁹. The Compass fundus perimeter (CMP, CenterVue, Padova, Italy) has extended FAP testing outside the macular region, to include conventional perimetric test patterns employed for testing glaucoma patients, such as the 24 − 2 of the Humphrey Field Analyzer (HFA, Carl Zeiss Meditec, Dublin, CA)¹⁰.

Despite advances in perimetric hardware and software, variability remains a major hurdle in assessing VF progression. FAP has proven useful, compared to SAP, in reducing test-retest variability for global indices ¹⁰. However, measurements at individual locations retained the typical increased variability at lower sensitivities observed with SAP¹⁰. Although several aspects can contribute to long-term variability, the systematic relationship between sensitivity and variability is thought to be a consequence of the shallower psychometric function¹¹ in locations with lower sensitivity. As a consequence, methods to assess pointwise progression of VF damage need to account for such differences in variability. Event-based analyses are common methods to assess point-wise progression. An example is the Guided Progression Analysis (GPA) algorithm implemented for the HFA which seeks to identify changes in sensitivity-related metrics that fall outside the expected 95% test-retest limits for a specific location, based on an estimate of baseline damage. Because this analysis is conducted pointwise, there is an increased chance of identifying locations with change exceeding the test-retest limits, even in stable patients, leading to False Positives (FP). Empirical rules have been implemented to curb this phenomenon, requiring a given number of locations to be outside the limits to classify progression. However, these empirical rules do not account for the fact that variability at more damaged locations is so large that the 95% lower limits would fall outside the range testable by the VF device¹². These locations would therefore not contribute to the FPR and accounting for this characteristic could improve the detection of progression in more damaged VFs without compromising specificity.

For this study, we prospectively collected test-retest series with the CMP fundus perimeter in a cohort of 90 stable glaucoma patients, stratified by VF damage. We used these data to evaluate a novel event-based approach that accounts for the number of locations in which progression can effectively be observed (adaptive rule) and compare this to a fixed decision rule and the empirical GPA decision thresholds. Progression was simulated with realistic rates of VF decay from an independent cohort of glaucoma patients.

Sample description

Ninety eyes of 90 subjects (55 males, 35 females) were recruited. The sample was composed of 34 (37.8%) eyes with early glaucomatous loss, 28 (31.1%) eyes with moderate glaucomatous loss and 28 (31.1%) with advanced loss. The characteristics of the study sample are summarized in Table 1. All patients had 6 reliable VF tests but only 5 were used (the first practice test was discarded).

Table 1

**Descriptive statistics of the test-retest sample.** IQR = Interquartile range; MD = Mean Deviation; PSD = Pattern Standard Deviation; BCVA = Best Corrected Visual Acuity; IOP = Intra Ocular Pressure; Early glaucomatous loss: -6 dB ≤ MD ≤ 0 dB; Intermediate glaucomatous loss: -12 dB ≤ MD < -6 dB; Advanced glaucomatous loss: MD < -12 dB.
	Early loss (N = 34)	Intermediate loss (N = 28)	Advanced loss (N = 28)
Sex (Male:Female)	19:16	17:11	19:8
Age (years)	65 [58.5, 72.5]	68 [64, 73.5]	66 [60.5, 75.5]
BCVA (logMAR)	0.0 [0.0, 0.02]	0.0 [0.0, 0.10]	0.0 [0.0, 0.10]
24 − 2 MD (dB)	-3.5 [-5.2, -2.0]	-8.2 [-10.5, -6.1]	-13.6 [-18.9, -11.2]
24 − 2 PSD (dB)	7.2 [4.2, 8.9]	11.1 [9.8, 13.0]	12.4 [10.8, 13.2]

Test-Retest variability limits

Figure 1 shows the test-retest distribution for the sensitivity, the Total deviation (TD) and the Pattern Deviation (PD) values, stratified by baseline sensitivity. The distribution was obtained from all 10800 permutations of the test-retest series, using all possible pairs of VFs to calculate the baseline sensitivity. As expected, the test-retest variability increased at lower sensitivities with all metrics. The floor effect is clearly visible for sensitivity values for the average baseline sensitivity < 15 dB.

Event analysis

The average pointwise progression rate for the matched sample derived from the Longitudinal Glaucomatous VFs data set (LGVFs) dataset (used for the simulations) was − 0.29 +- 1.24 dB/year. Figures 2 and 3 report the results of the event detection with the three algorithms for sensitivities, TD and PD as cumulative percentage of Hit rate (HR) and False Positive Rate (FPR) at each visit. The shaded bands represent the 95% quantiles around the central estimate from different realisations of the simulation (i.e., permutations of the template/field sequences), used to define the variability of the estimate.

The FPR for one event (E1) at first follow-up (6 months) was in general close to the desired 5%. However, it was on average larger when the adaptive rule was applied to the sensitivities (8.2%) and TD maps (8.9%) in the group with intermediate damage and for the PD maps for eyes with early damage (10.7% for the fixed rule and 10.2% for the adaptive rule). As expected, the FPR was much lower than 5% for the advanced cases with the fixed rule (0.4% for the sensitivities, 1.4% for TD and 0.6% for PD). However, the 5% limit was outside the 95% variability range of the estimate only for the PD maps in subjects with early damage and in the advanced group for the fixed rule with any map. Beyond the first follow-up, the FPR for E1 steadily increased, as expected. The HR was generally higher for the adaptive rule. Importantly, the HR for the advanced cases when the adaptive rule was applied to the sensitivities (7.7%) and the TD values (9.8%) was close to that obtained in the overall sample (10.7% for sensitivities, 11.7% for TD). The HR was instead extremely reduced for the fixed rule (1.1% for sensitivities, 2.3% for TD) and the GPA (0% for sensitivities, 0% for TD).

The FPR for two events (E2) at the second follow-up (12 months) was lower or very close to 5% in all cases and was very similar among all detection methods. The only meaningful difference, as expected, was observed for the advanced cases where the FPR was higher with the adaptive rule applied to the sensitivities (1.6% for the adaptive rule, and 0.2% for the fixed rule and GPA) and TD values (1.9% for the adaptive rule, and 0.5% for the fixed rule and GPA). Concordantly, the HR was also better for the adaptive rule in advanced cases (5.6% for the sensitivities and 6.7% for the TD), doubled compared to the fixed rule and the GPA (3% for the sensitivities and 3.3% for the TD for both methods).

For three events (E3) at last follow-up (18 months) the fixed and adaptive rules behaved very similarly across the whole spectrum of damage both in terms of FPR and HR, with marginally higher FPR and HR compared to GPA.

Tables with the full numerical results are provided in the supplementary material. The average (SD) number of locations available for progression was 45.7 (7.9) for patients with early glaucomatous loss, 37 (9) for patients with intermediate glaucomatous loss and 25.5 (11) for patients with advanced glaucomatous loss.

We compared different algorithms to detect event-based progression. We derived our data from test-retest series in 90 glaucoma patients, stratified by VF damage, using a fundus tracking perimeter that compensates for eye movements. Moreover, the use of fundus tracking allowed the same retinal locations to be tested in all follow-up tests. Our simulations relied on a method that accounts for global fluctuations in sensitivity across different locations in the VF, estimated from the data¹³. This is particularly important to accurately estimate the performance of event-based analyses, which instead assume each location is independent.

The observed test-retest distribution for the Compass was very similar to that reported in a previous study with a smaller dataset¹⁰. The variability increased substantially for sensitivities below 25 dB, in accordance with previous literature^12,14,15, reaching the floor for the lower 5% limits at 15 dB. Interestingly, such a value has often been reported as a lower bound for locations to be informative for progression detection^12,16. Indeed, some newly implemented testing strategies have also considered not testing locations with sensitivities below 15 dB ¹⁷ because of their limited usefulness. Our data confirm this view, showing that such locations are not informative for event-based progression. The test-retest variability was very similar between sensitivities, TD and PD. It is important to notice that, differently from the Swedish interactive thresholding algorithm (SITA) employed by the HFA, the CMP does not make use of any spatial smoothing or growth patterns in the testing strategy^10,18,19. Testing each location independently could lead to an increase in test-retest variability at lower sensitivities despite fundus tracking as reported in our previous study¹⁰. However, locations below 15 dB have been shown to be scarcely informative for progression also with SITA strategies^12,20, implying that such differences are unlikely to be relevant for practical applications. Differently from previous analyses and from the traditional GPA¹⁴, the expected lower 5% limits for TD and PD was modelled based on the average baseline sensitivity rather than the average TD or PD values. This reflects the idea that sensitivity, rather than TD or PD, is the main determinant of response variability¹¹. It is important to note, however, that locations with the same sensitivity have been shown to have very different variability, such as in Chauhan et al.²¹. In their analysis, however, sensitivity was still a better parameter than TD to characterise systematic changes in variability. This supports our choice of baseline sensitivity as the best, albeit imperfect, parameter to model the expected test-retest variability.

When comparing the different methods of event-based analysis, we found that they performed similarly in the overall sample. The GPA for E3 offered a lower HR with minor improvement in the FPR compared to the other methods. Compared to previous reports on the performance of the GPA^14,22−24, our analysis yielded similar values for the FPR and HR for “likely progression” detection (0% and 3.7% respectively in the overall sample with PD values). However, for “probable progression” the specificity was higher than previously reported. Artes et al.²³ showed a cumulative FPR close to 5% at 2 and 3 cumulative tests, similar to our results. However, Wu et al.²² reported a FPR of approximately 8% at 2 and 3 cumulative tests. In our simulations, the FPR with the GPA “probable progression” never exceeded the 5% threshold, despite using a similar simulation methodology. Of course, several differences could have contributed to this discrepancy, including the different perimeter employed and the use of a tobit model to estimate and simulate progression instead of a logistic function and the different sample composition (see supplementary material). The use of the same dataset to define the test-retest limits and to test our detection strategies could have influenced the measured performance, but this effect was minimized by the use of a leave-one-out approach in the simulations. One novel aspect of our analysis is that we grouped locations based on the average baseline sensitivity to calculated the lower 5% limits for event detection. This is different from how variability for the PD values was handled for the GPA¹⁴, which was instead stratified according the average of baseline PD values. Our decision was taken in recognition of the fact that the major determinant of perimetric variability is the sensitivity¹¹, independently of the general height or normal ageing.

Similarly to us, Wu et al. showed that the GPA had a smaller HR for more damaged fields, but did not investigate the corresponding change in specificity²². With our analysis, we show that such a reduction in the HR is connected to an unwanted increase in specificity at more advanced stages, derived from the static rules applied to detect progression. Indeed, such a feature was shared by our fixed rule detection. On the contrary, the adaptive rule maintained a more homogeneous FPR across the severity spectrum. This effect was absent for the E3 detection, which was unchanged between the fixed and adaptive rules, but resulted in important improvements for the advanced cases for the E1 and E2 detections. Interestingly, for early and advanced cases, the FPR for E1 detection was below the 5% threshold with the fixed rule and very close to 5% for the adaptive rule in the advanced group. The 5% limit was also always well within the 95% variability range of the FPR estimate form the simulations for all classes of damage. For the advanced class, this was accompanied by an important increase in HR. This has clinical relevance for the advanced cases, in which faster detection of VF deterioration would be desirable to prevent progression to significant visual impairment or blindness. Progression with E1 criteria had low specificity for patients with intermediate damage using sensitivity and TD values. This could be explained by the gradual shift towards more variable sensitivities. Importantly, except for the E1 criterion, both the HR and the FPR were virtually identical between the fixed and adaptive rule for intermediate and early patients. Moreover, both the fixed and adaptive rule generally outperformed the GPA, especially for the 3-Event detection, in terms of HR (higher) with minimal changes to the FPR, which remained < 5%. It must be noted that most of the methodology for the GPA was developed in the context of the Early Manifest Glaucoma Trial (EMGT)²⁵ which was meant to detect progression in patients with early glaucomatous damage^25,26. It is therefore unsurprising that such a methodology would perform poorly in intermediate and advanced cases, proving often too conservative for clinical applications.

The performance of detection with PD was worse for all detection criteria compared to that with sensitivities and TD. The detrimental effect of using PD was particularly evident for patients with early damage, in which specificity was generally lower compared to sensitivities and TD, and in advanced cases, in which the HR was reduced. This is in agreement with previous reports, showing that PD based metrics might hide progression in glaucoma patients^12,24. One aspect to consider is that our dataset did not contain patients with visually significant cataract. However, the reported effect of lens opacity on the VF has been variable, ranging from significant ^27–30 to negligible ^31,32, with a measured effect on global indices, such as the MD, of approximately 1 dB. Therefore, the role of metrics meant to account for generalised VF depression, such as the PD, remains unclear in the context of determining glaucoma progression, also in consideration of the potential masking effect of such indices²⁴ and their limited reliability in advanced glaucoma damage.

One limitation of our analysis is that the length of the test-retest series did not allow for permutations/simulations of more than 5 visual fields per patient. Further differences between the methods might arise in longer series. However, this study was meant to test the performance of event-based analyses for early detection of progression, where they would be more clinically relevant. Moreover, longitudinal VF series collected with the Compass in glaucoma patients are not yet available and simulation of progression had to rely on data collected with an HFA^33,34. However, we based our simulations on noise estimates from the observed test-retest distribution, in an attempt to reliably mimic what would really be observed in a real-life use of the device. A collection of longitudinal series in glaucoma patients using the CMP is underway and will be able to confirm or refine our results.

In conclusion, our proposed adaptive rule has the potential to be more useful in the clinical practice than current implementations of the event analysis, such as the GPA, especially in patients with advanced damage. Sensitivities and TD values offered similar performance, but TD has the potential for wider application by accounting for normal ageing. Fundus-tracked perimetry can help ensure that the same location is tested in successive follow-ups, a core assumption in any event analysis. However, the effective clinical improvement brought by fundus tracking in detecting progression needs to be determined with longitudinal VF series on progressing glaucoma patients.

Data collection for test-retest variability

Data were collected as part of a longitudinal open-label international multicentre study. Five sites took part in the study: IRCCS Fondazione “G. B. Bietti”, Rome, Italy; Polyclinic Città di Udine, Udine, Italy; Ophthalmology Service and Glaucoma Center HIA Begin. Paris, France; Eye Care Centre, Nova Scotia Health, Halifax, Canada; and ASST Santi Paolo e Carlo, Milan, Italy. Recruitment started on August 30, 2018 and ended on December 20, 2019. We enrolled 90 consecutive patients with glaucoma, testing one eye per subject. The study eye was randomly selected if both eyes met the inclusion criteria. Eyes were stratified into three groups according to disease severity, using the 24 − 2 Mean Deviation (MD) of the most recent reliable HFA test: (1) early glaucomatous loss (-6 dB ≤ MD ≤ 0 dB); (2) moderate glaucomatous loss (-12 dB ≤ MD < -6 dB); (3) advanced glaucomatous loss (MD < -12 dB).

Inclusion criteria were as follows: (1) Age between 18 and 80 years; (2) Best corrected visual acuity logMar ≤ + 0.1 (≥ 8/10 decimal); (3) Spherical refraction between − 12D and + 15D (limits of autofocus in the CMP) ; (4) Astigmatism between − 2D and + 2D; (5) Clinically determined glaucomatous optic nerve damage; (6) Abnormal optical coherence tomography (OCT) scan of the retinal nerve fibre layer (RNFL) and optic nerve head (ONH) in at least one sector, with reference to the normative values specific of each device used; (7) At least 5 available reliable VF tests within the last 5 years; and (8) repeatable VF defects in the last 2 reliable VF tests, defined as: (i) MD or pattern standard deviation (PSD) outside 95% confidence limits (p < 0.05); (ii) Glaucoma Hemifield Test (GHT) outside normal limits; and (iii) A cluster of at least 3 points with p < 0.05 in the PSD map, one of each with p < 0.01 affecting the same hemifield, however, the cluster could not be contiguous with the blind spot and could not cross the horizontal midline.³⁵

Exclusion criteria were as follows: (1) Any ocular surgery except uncomplicated cataract surgery and/or glaucoma surgery performed less than 6 months prior to data collection; (2) Any ocular pathology that could affect the VF other than glaucoma; (3) Use of any drug that could hinder the perimetric examination or its results; and (4) Angle-closure glaucoma.

Eligible patients were identified from the clinical registry of the glaucoma clinics in each recruiting centre. All patients gave written informed consent to participate in the study. Ethics Committee approval was obtained for each centre from the local ethical committee (Approval from the coordinating centre: Comitato Etico Centrale IRCCS Lazio, Sezione IRCCS I.F.O., June, 26, 2018, ref: N. 68 /18/FB). This study adhered to the tenets of the Declaration of Helsinki.

At the first visit, each subject underwent the following examinations in both eyes (besides the VF test): autorefraction, best corrected visual acuity, clinical assessment of the optic disc with slit lamp/ophthalmoscopy/fundus picture, OCT of the RNFL and ONH if a recent scan (< 3 months from data collection) was not available, Goldman applanation tonometry and gonioscopy.

Visual field testing protocol

VF examinations were performed with the CMP using a ZEST (Zippy Estimation through Sequential Testing)^18,36 strategy and a 24 − 2 + grid. The 24 − 2 + grid has the same 52 locations of the standard 24 − 2 grid (excluding the two blind spot locations) with 12 additional points in the macular region (Fig. 4). If the patient was eligible, a first test (Practice) and 5 additional retests were performed within a period of about one month (average 31 days, range [16–44] days) with a mean interval of 6 days (range [1–15] days) between tests. All tests were considered reliable if the rate of false positive errors was ≤ 18%, the rate of false negative errors was ≤ 30% and the average pupil diameter during each test was > 2.8 mm (the threshold to ensure reliable imaging and tracking with the CMP). The tests were performed in follow-up mode, so that the same retinal locations were tested in each examination.

The Practice test was not included in the analysis and served to familiarise the patient with the testing procedure.

Simulation of progressing fields

Simulations were performed using the method proposed by Wu et al.¹³ to obtain longitudinal visual field series. The details are reported in the supplementary material. Briefly, the method uses noise templates to simulate variability in the VFs while maintaining correlations across locations within the same test. This is an important factor in VF variability, sometimes referred to as the Global Visit Effect ³⁷, since the fluctuations in sensitivity at different visits are not independent across different locations. The method also accounts for the larger variability at lower sensitivity values and derives the noise distributions from the empirical test-retest data to specifically simulate the noise observed in CMP from this data collection. The pointwise slopes for sensitivity change over time were derived from the LGVFs from the Rotterdam Ophthalmic data repository^33,38, a collection of 278 VF series from glaucoma patients tested with the Full Threshold (FT) algorithm with the HFA. Both FT and the ZEST procedure of CMP derive sensitivity at test locations independently and without reference to neighbouring points. The pointwise slopes were calculated through tobit regression³⁴ in R (R Foundation for statistical computing, Vienna, Austria) using the AER package³⁹ to account for the floor effect at lower sensitivities (see supplementary material).

From the LGVFs, we selected 90 subjects matched with the participants of our test-retest dataset according to their baseline Mean Sensitivity (MS). For the CMP test-retest series, the MS was calculated as the mean of the 52 sensitivity values from the 24 − 2 locations, each estimated as the average value across all the VFs in the test-retest series from the same eye. For the LGVFs, the MS was the mean of the 52 values predicted by the tobit regression at baseline (first visit). The pointwise slopes from the matched LGVFs subjects were used to simulate progression in the eyes from our test-retest dataset. For our simulations, each VF in the original test-retest series was transformed into a noise template (see supplementary material). We simulated all possible combinations of noise templates (i.e. VFs) for each eye, for a total of 10800 simulated VFs series (120 per eye). This allowed us to retain the VF fluctuations actually observed in each real test-retest series. Essentially, the simulated progressing VF series replicated the permutations of the test-retest series, with the additional effect of the simulated decay in sensitivity. We simulated two baseline visits at 0 and 6 months, and then three more VF at 12, 18 and 24 months (total of 5 VFs, corresponding to the 5 noise templates per eye).

Event analysis

An event analysis method was developed to detect significant progression events. Based on the method proposed by Heijl et al.¹⁴, we used the average of two VFs in the test-retest series to define the baseline sensitivity for each location in each eye. The lower 5% limits for the test-retest distribution were calculated for each baseline sensitivity and used to detect significant progression events. The same calculations were repeated using the TD and the PD values, as provided by the CMP using an internal normative database¹⁰. Of note, when calculating the 5% limits, TD and PD values were grouped according to their average baseline sensitivity instead of the average TD and PD baseline values. This choice was preferred since sensitivity is known to be the main determinant of response variability in perimetry¹¹. The limits were calculated for each eye by excluding all data from that same eye, to avoid inflating the performance (leave-one-out). Additional details are provided in the supplementary material.

We compared different criteria to define a VF series as progressed. Progression could be identified by 1 event (E1) or two (E2) or three (E3) consecutive events at the same location. A set of fixed decision rules were used to combine these events at multiple locations:

E1 progression: one event is observed in at least 6 locations
E2 progression: two consecutive events are observed in at least 3 locations
E3 progression: three consecutive events are observed in at least 2 locations

Similar to the GPA, the fixed rule is meant to account for the rate of expected false positives in the detection of the events, considering each location (N = 52) as a separate statistical test with a 5% probability of false positive error. The probability of detecting at least one false event was assumed to follow a binomial distribution (see supplementary material for details). However, the fixed rule does not account for locations where detecting an event is impossible (locations with sensitivity at Baseline below 15 dB), because the lowest 5% quantile of the test-retest distribution coincides with the lower measurement limit (0 dB). These locations do not contribute to the overall false positive rate and are more frequent in eyes with advanced damage, because lower sensitivities have larger test-retest variability (Fig. 5). We therefore proposed a modification of this rule (the second algorithm in our comparison, called adaptive rule) where the minimum number of locations with an observed event required to detect progression is adjusted according to the effective number of locations in which progression can be measured (potentially less than 52). Details are reported in the supplementary material. Finally, we compared our results with the empirical rules reported for GPA^26,40,41: no VF progression is based on E1; E2 and E3 progression is assessed in both cases if two (“possible progression”) or three (“likely progression”) consecutive events respectively are observed in at least 3 locations. This differs from our fixed E3 criterion, where only 2 locations are required.

We applied these algorithms to all permutations of test-retest series, assumed stable, to evaluate the FPR and to all simulated progressing series to measure the rate of detection (Hit Rate, HR). For each series, the first two tests were used as baselines for the event analysis, over a total of 5 tests. Therefore, E1 progression could be detected for test 3 to 5, E2 for tests 4 and 5 and E3 only for test 5. We considered 5% as the ideal benchmark for FPR for clinical use. The variability of our estimates is reported using the 5%-95% quantiles of the results obtained from different realisations of our simulations (i.e., each permutation of the VF series).

Acknowledgements

The Contribution of Fondazione Bietti to this paper was supported by the Italian Ministry of Health and by Fondazione Roma

Author Contributions

C. R, G. M, D. C developed the original research concept. P. B, B. C. C, L. M. R, P. F, J.M.G, J.R. F and F. O performed data collection. C. R and G. M performed the analysis. C. R and G. M wrote the main manuscript text. All authors critically reviewed and edited the manuscript.

Additional Information

Competing Interest: The author(s) have made the following disclosure(s): G.M, D.C, B.C.C, L.M.R, P.F, F.O: Consultants - CenterVue, Padua, Italy. C.R is an employee of Centerve SpA, Padova, Italy. All other authors have declared no competing interests. CenterVue SpA, Padova, Italy, sponsored the study.

1 Kass, M. A. et al. The Ocular Hypertension Treatment Study: a randomized trial determines that topical ocular hypotensive medication delays or prevents the onset of primary open-angle glaucoma. Arch Ophthalmol 120, 701-713; discussion 829-730, doi:10.1001/archopht.120.6.701 (2002).

2 Heijl, A. et al. Reduction of intraocular pressure and glaucoma progression: results from the Early Manifest Glaucoma Trial. Arch Ophthalmol 120, 1268-1279, doi:10.1001/archopht.120.10.1268 (2002).

3 Comparison of glaucomatous progression between untreated patients with normal-tension glaucoma and patients with therapeutically reduced intraocular pressures. Collaborative Normal-Tension Glaucoma Study Group. Am J Ophthalmol 126, 487-497, doi:10.1016/s0002-9394(98)00223-2 (1998).

4 Anderson, D. R. & Normal Tension Glaucoma, S. Collaborative normal tension glaucoma study. Curr Opin Ophthalmol 14, 86-90, doi:10.1097/00055735-200304000-00006 (2003).

5 Miglior, S. et al. Results of the European Glaucoma Prevention Study. Ophthalmology 112, 366-375, doi:10.1016/j.ophtha.2004.11.030 (2005).

6 Musch, D. C. et al. Visual field progression in the Collaborative Initial Glaucoma Treatment Study the impact of treatment and other baseline factors. Ophthalmology 116, 200-207, doi:10.1016/j.ophtha.2008.08.051 (2009).

7 The Advanced Glaucoma Intervention Study (AGIS): 7. The relationship between control of intraocular pressure and visual field deterioration.The AGIS Investigators. Am J Ophthalmol 130, 429-440, doi:10.1016/s0002-9394(00)00538-9 (2000).

8 Garway-Heath, D. F. et al. Latanoprost for open-angle glaucoma (UKGTS): a randomised, multicentre, placebo-controlled trial. Lancet 385, 1295-1304, doi:10.1016/S0140-6736(14)62111-5 (2015).

9 Pfau, M. et al. Fundus-controlled perimetry (microperimetry): Application as outcome measure in clinical trials. Prog Retin Eye Res, 100907, doi:10.1016/j.preteyeres.2020.100907 (2020).

10 Montesano, G. et al. A Comparison between the Compass Fundus Perimeter and the Humphrey Field Analyzer. Ophthalmology 126, 242-251, doi:10.1016/j.ophtha.2018.08.010 (2019).

11 Henson, D. B., Chaudry, S., Artes, P. H., Faragher, E. B. & Ansons, A. Response variability in the visual field: comparison of optic neuritis, glaucoma, ocular hypertension, and normal eyes. Invest Ophthalmol Vis Sci 41, 417-421 (2000).

12 Artes, P. H., Iwase, A., Ohno, Y., Kitazawa, Y. & Chauhan, B. C. Properties of perimetric threshold estimates from Full Threshold, SITA Standard, and SITA Fast strategies. Invest Ophthalmol Vis Sci 43, 2654-2659 (2002).

13 Wu, Z. & Medeiros, F. A. Development of a Visual Field Simulation Model of Longitudinal Point-Wise Sensitivity Changes From a Clinical Glaucoma Cohort. Transl Vis Sci Technol 7, 22, doi:10.1167/tvst.7.3.22 (2018).

14 Heijl, A., Lindgren, A. & Lindgren, G. Test-retest variability in glaucomatous visual fields. Am J Ophthalmol 108, 130-135, doi:10.1016/0002-9394(89)90006-8 (1989).

15 Chauhan, B. C. & Johnson, C. A. Test-retest variability of frequency-doubling perimetry and conventional perimetry in glaucoma patients and normal subjects. Invest Ophthalmol Vis Sci 40, 648-656 (1999).

16 Junoy Montolio, F. G., Wesselink, C. & Jansonius, N. M. Persistence, spatial distribution and implications for progression detection of blind parts of the visual field in glaucoma: a clinical cohort study. PLoS One 7, e41211, doi:10.1371/journal.pone.0041211 (2012).

17 Turpin, A., Morgan, W. H. & McKendrick, A. M. Improving Spatial Resolution and Test Times of Visual Field Testing Using ARREST. Transl Vis Sci Technol 7, 35, doi:10.1167/tvst.7.5.35 (2018).

18 Turpin, A., McKendrick, A. M., Johnson, C. A. & Vingrys, A. J. Properties of perimetric threshold estimates from full threshold, ZEST, and SITA-like strategies, as determined by computer simulation. Invest Ophthalmol Vis Sci 44, 4787-4795, doi:10.1167/iovs.03-0023 (2003).

19 Bengtsson, B., Olsson, J., Heijl, A. & Rootzen, H. A new generation of algorithms for computerized threshold perimetry, SITA. Acta Ophthalmol Scand 75, 368-375, doi:10.1111/j.1600-0420.1997.tb00392.x (1997).

20 Gardiner, S. K., Swanson, W. H., Goren, D., Mansberger, S. L. & Demirel, S. Assessment of the reliability of standard automated perimetry in regions of glaucomatous damage. Ophthalmology 121, 1359-1369, doi:10.1016/j.ophtha.2014.01.020 (2014).

21 Chauhan, B. C., Tompkins, J. D., LeBlanc, R. P. & McCormick, T. A. Characteristics of frequency-of-seeing curves in normal subjects, patients with suspected glaucoma, and patients with glaucoma. Invest Ophthalmol Vis Sci 34, 3534-3540 (1993).

22 Wu, Z. & Medeiros, F. A. Comparison of Visual Field Point-Wise Event-Based and Global Trend-Based Analysis for Detecting Glaucomatous Progression. Transl Vis Sci Technol 7, 20, doi:10.1167/tvst.7.4.20 (2018).

23 Artes, P. H., O'Leary, N., Nicolela, M. T., Chauhan, B. C. & Crabb, D. P. Visual field progression in glaucoma: what is the specificity of the Guided Progression Analysis? Ophthalmology 121, 2023-2027, doi:10.1016/j.ophtha.2014.04.015 (2014).

24 Artes, P. H., Nicolela, M. T., LeBlanc, R. P. & Chauhan, B. C. Visual field progression in glaucoma: total versus pattern deviation analyses. Invest Ophthalmol Vis Sci 46, 4600-4606, doi:10.1167/iovs.05-0827 (2005).

25 Heijl, A. et al. Measuring visual field progression in the Early Manifest Glaucoma Trial. Acta Ophthalmol Scand 81, 286-293, doi:10.1034/j.1600-0420.2003.00070.x (2003).

26 Leske, M. C., Heijl, A., Hyman, L. & Bengtsson, B. Early Manifest Glaucoma Trial: design and baseline data. Ophthalmology 106, 2144-2153, doi:10.1016/s0161-6420(99)90497-9 (1999).

27 Chen, P. P. & Budenz, D. L. The effects of cataract extraction on the visual field of eyes with chronic open-angle glaucoma. Am J Ophthalmol 125, 325-333, doi:10.1016/s0002-9394(99)80142-1 (1998).

28 Hayashi, K., Hayashi, H., Nakao, F. & Hayashi, F. Influence of cataract surgery on automated perimetry in patients with glaucoma. Am J Ophthalmol 132, 41-46, doi:10.1016/s0002-9394(01)00920-5 (2001).

29 Koucheki, B., Nouri-Mahdavi, K., Patel, G., Gaasterland, D. & Caprioli, J. Visual field changes after cataract extraction: the AGIS experience. Am J Ophthalmol 138, 1022-1028, doi:10.1016/j.ajo.2004.08.006 (2004).

30 Rao, H. L., Jonnadula, G. B., Addepalli, U. K., Senthil, S. & Garudadri, C. S. Effect of cataract extraction on Visual Field Index in glaucoma. J Glaucoma 22, 164-168, doi:10.1097/IJG.0b013e31822e8e37 (2013).

31 Stewart, W. C., Rogers, G. M., Crinkley, C. M. & Carlson, A. N. Effect of cataract extraction on automated fields in chronic open-angle glaucoma. Arch Ophthalmol 113, 875-879, doi:10.1001/archopht.1995.01100070049024 (1995).

32 Carrillo, M. M., Artes, P. H., Nicolela, M. T., LeBlanc, R. P. & Chauhan, B. C. Effect of cataract extraction on the visual fields of patients with glaucoma. Arch Ophthalmol 123, 929-932, doi:10.1001/archopht.123.7.929 (2005).

33 Bryan, S. R., Vermeer, K. A., Eilers, P. H., Lemij, H. G. & Lesaffre, E. M. Robust and censored modeling and prediction of progression in glaucomatous visual fields. Invest Ophthalmol Vis Sci 54, 6694-6700, doi:10.1167/iovs.12-11185 (2013).

34 Russell, R. A. & Crabb, D. P. On alternative methods for measuring visual field decay: Tobit linear regression. Invest Ophthalmol Vis Sci 52, 9539-9540, doi:10.1167/iovs.11-8948 (2011).

35 Oddone, F. et al. Influence of disc size on optic nerve head versus retinal nerve fiber layer assessment for diagnosing glaucoma. Ophthalmology 118, 1340-1347, doi:10.1016/j.ophtha.2010.12.017 (2011).

36 King-Smith, P. E., Grigsby, S. S., Vingrys, A. J., Benes, S. C. & Supowit, A. Efficient and unbiased modifications of the QUEST threshold method: theory, simulations, experimental evaluation and practical implementation. Vision Res 34, 885-912, doi:10.1016/0042-6989(94)90039-6 (1994).

37 Bryan, S. R., Eilers, P. H., Lesaffre, E. M., Lemij, H. G. & Vermeer, K. A. Global Visit Effects in Point-Wise Longitudinal Modeling of Glaucomatous Visual Fields. Invest Ophthalmol Vis Sci 56, 4283-4289, doi:10.1167/iovs.15-16691 (2015).

38 Erler, N. S. et al. Optimizing structure-function relationship by maximizing correspondence between glaucomatous visual fields and mathematical retinal nerve fiber models. Invest Ophthalmol Vis Sci 55, 2350-2357, doi:10.1167/iovs.13-12492 (2014).

39 Kleiber, C. & Zeileis, A. Applied Econometrics with R. Springer-Verlag (2008).

40 Arnalich-Montiel, F., Casas-Llera, P., Munoz-Negrete, F. J. & Rebolleda, G. Performance of glaucoma progression analysis software in a glaucoma population. Graefes Arch Clin Exp Ophthalmol 247, 391-397, doi:10.1007/s00417-008-0986-1 (2009).

41 Heijl, A., Patella, V. M. & Bengtsson, B. The Field Analyzer Primer:Effective Perimetry. (2012).

Competing interest reported. The author(s) have made the following disclosure(s): G.M, D.C, B.C.C, L.M.R, P.F, F.O: Consultants - CenterVue, Padua, Italy. C.R is an employee of Centerve SpA, Padova, Italy. All other authors have declared no competing interests. CenterVue SpA, Padova, Italy, sponsored the study.

SupplementaryMaterialFinalProgression.docx

PDF

Journal Publication

published 11 Aug, 2021

Read the published version in Scientific Reports →

Editorial decision: Major revision
25 May, 2021
Reviews received at journal
24 May, 2021
Reviewers agreed at journal
14 May, 2021
Reviewers invited by journal
12 May, 2021
Editor assigned by journal
07 May, 2021
Editor invited by journal
22 Apr, 2021
Submission checks completed at journal
21 Apr, 2021
First submitted to journal
21 Apr, 2021

You are reading this latest preprint version

A novel approach to event-based progression analysis in glaucomatous visual fields

Status:

Journal Publication

Version 1

Abstract

Figures

Introduction

Results

Discussion

Methods

Declarations

References

Competing Interests

Supplementary Files

Status:

Journal Publication

Version 1