Comparison of Optical Coherence Tomography Angiography Metrics in Primary Angle-Closure Glaucoma and Normal-Tension Glaucoma

DOI: https://doi.org/10.21203/rs.3.rs-121226/v1

Abstract

Purpose: To investigate the peripapillary vascular metrics in early normal tension glaucoma (NTG) and early primary angle closure glaucoma (PACG) eyes using optical coherence tomography angiography (OCT-A).

Methods: One or both eyes of each subjects were imaged for a 3x3mm peripapillary region by swept-source OCT-A (DRI-OCT Triton, Topcon, Japan) and assessed by an automated MATLAB program. OCT-A metrics including circumpapillary vessel density (cpVD) and fractal dimension (cpFD) were compared. Their association with visual field (VF) parameters and retinal nerve fiber layer (RNFL) thickness were determined.

Results: Sixty-eight eyes of 51 PACG, 68 eyes of 48 NTG, and 68 eyes of 49 control subjects were cross-sectionally analyzed. NTG eyes had significantly lower global cpVD (52.369±0.781%) compared with PACG eyes (55.389±0.721%, P=0.004) that had comparable disease severity and average RNFL thickness. Multivariable analysis reviewed that, for PACG and NTG eyes, decreased cpVD ([PACG] β=-4.242; CI:-8.120, -0.363 vs [NTG] β=-5.531; CI:-9.472, -1.590) and cpFD ([PACG] β=-8.894;CI:-11.925, -5.864 vs [NTG] β=-12.064; CI:-17.095, -6.932) were associated with decreased RNFL thickness (all P≤0.032); with a stronger association between decrease cpFD and decreased RNFL thickness in NTG eyes (Fisher’s Z-test, P=0.045). Decreased

cpVD was associated with decrease mean deviation (MD) in NTG eyes (β=-0.707; CI:-1.090, -0.324; P≤0.001) and not associated with the visual field parameters in PACG eyes.

Conclusions: Early NTG had lower global cpVD compared with early PACG, despite similar disease severity and average RNFL thickness. 

Précis:

Optical coherence tomography angiography observed a lower peripapillary microvascular perfusion in early NTG eyes compared with early PACG eyes, despite similar disease severity and retinal nerve fiber layer thickness. 

Introduction

Glaucoma is a heterogeneous group of progressive optic neuropathy characterized by degeneration of retinal ganglion cells (RGCs), resulting in a characteristic optic neuropathy and visual field (VF) loss.1 Some patients experienced continuous VF deterioration despite adequate intraocular pressure (IOP) control, especially in normal tension glaucoma (NTG).2 Studies suggested that vascular factors, other than the IOP-dependent risk factors, may be an important pathogenic mechanism of NTG.36 For instance, population-based studies demonstrated the association between retinal vascular changes observed on retinal photographs and glaucoma;7, 8 NTG eyes was also observed to have significantly narrower retinal arteriolar caliber compared to eyes with primary open angle glaucoma (POAG) and higher IOP (high-tension glaucoma, HTG).9 In contrast, primary angle closure glaucoma (PACG) is likely to be more IOP-dependent.10 If this were the case, one would expect a difference in the pattern of microvascular changes between the two glaucoma subtypes throughout the spectrum of different disease severity.

Optical coherence tomography angiography (OCT-A) allows direct, three dimensional, non-invasive, and depth-resolved visualization of the retinal and choroidal microvasculature through en face reconstruction of OCT combined with motion contrast processing.11 Hence, peripapillary optic disc perfusion and vessel density (VD) can be quantitatively measured. Lower disc flow index and VD have been observed in POAG eyes compared with normal controls;12, 13 and also correlated with more severe disease.14 Microvascular reduction was also associated with VF defects in a region-specific manner.15 These suggested that OCT-A may be potentially useful in detecting retinal microvascular change and vascular mechanism of glaucoma.13, 16

It is unclear whether the characteristic pattern of vascular changes, whilst associated with glaucoma, are the cause or consequence of RGCs loss. Vascular impairment per se could deplete the RGCs’ blood supply and lead to their loss; RGCs death per se could also lead to reduced metabolic demand and the consequential vascular reduction. Neither the comparison of glaucoma patients with normal subjects,15 nor the comparison of eyes with the same glaucoma subtype of different severities,14 could specifically address this question. Comparing the microvasculature of PACG eyes (a predominantly IOP-dependent glaucoma subtype) and NTG eyes (a predominantly non-IOP-dependent glaucoma subtype) could reveal the role of vascular mechanisms in the disease pathogenesis, especially at the early stage of the disease in which early vascular and retinal nerve fiber layer (RNFL) changes could be observed. The purpose of the present study was to evaluate peripapillary vasculature between early PACG and early NTG, as well as to compare the strength of associations in OCT-A vascular metrics with structural and functional glaucoma parameters in PACG and NTG.

Methods

Subjects

This is a cross-sectional clinical study. The study was conducted in accordance with the ethical standards stated in the 2013 Declaration of Helsinki and approved by Hong Kong Kowloon Central Research Ethics Committee with written informed consent obtained. Patients with PACG, NTG, and normal control were recruited from June 2018 to May 2019 at Hong Kong Eye Hospital and the CUHK Eye Centre. The PACG subjects were recruited from an ongoing population based study, the CUHK PACG Longitudinal (CUPAL) study, which was described in detail previously.17, 18 One or both eyes of a subject could be enrolled.

Clinical Examinations

A full ophthalmic assessment was carried out. This included measurement of best-corrected visual acuity (BCVA), refractive error with an autorefractor (Nidek ARK-510A, Gamagori, Japan), axial length (AL) as well as anterior chamber depth (ACD) with A-scan biometer (AL-100, Tomey, Nagoya, Japan), IOP with Goldmann applanation tonometry, central cornea thickness (CCT) with non-contact tonopachymeter (TONOPACHY™ 530P, Nidek co., Ltd., Gamagori, Japan), dark room gonioscopy, anterior segment slit-lamp biomicroscopy, dilated fundal examination with assessment of the optic disc (20D and 90D lens), and VF examination by static automated white-on-white threshold perimetry (Humphrey Field Analyzer; 24 − 2 Swedish interactive threshold algorithm; Carl Zeiss Meditec, Dublin, CA, USA). IOP was measured twice for each eye. If the two readings differed by ≤ 2 mmHg, the mean was recorded as the IOP measurement. Otherwise, a third reading was performed, and the mean was recorded. Gonioscopy was carried out in a dark room. Spherical equivalent (SE) refraction was calculated as the sum of the spherical value and half of the cylindrical value. Peripapillary RNFL thickness was measured with Spectralis SD-OCT (HRA + OCT, Heidelberg Engineering, Heidelberg, Germany). OCT-A imaging was taken by a swept-source OCT (detailed in the session “OCT-A imaging”).

Definitions

All NTG and PACG eyes had structural and functional evidence of glaucoma, including glaucomatous optic disc cupping, thinning of the RNFL, loss of neuroretinal rim, and minimal criteria for glaucomatous VF defect as per published standard:2 glaucoma hemifield test result outside normal limits, pattern standard deviation (PSD) with P < 0.05 or a cluster of 3 or more points in the pattern deviation plot in a single hemifield with P < 0.05, one of which must have P < 0.01. Any one of the preceding criteria, if repeatable, was considered sufficient evidence of a glaucomatous VF defect. Additional criteria of NTG eyes were adopted from the Collaborative Normal Tension Glaucoma study2 that included: (1) six median untreated IOP readings consistently less than 21 mmHg, with no more than 1 reading equal to 23 or 24 mmHg and no single measurement more than 24 mmHg; at least 2 readings were obtained at a different time of the day from the rest, and (2) open drainage angle (Shaffer grade II or above) on dark room gonioscope. Eyes with PACG had (1) a total of ≥ 180o of angle closure obliterating the trabecular meshwork (synechial or appositional), and (2) untreated IOP of > 21 mmHg.19 Patients with early glaucoma were recruited in this study according to the modified Hodapp-Anderson-Parrish (HAP) staging20 – mean deviation (MD) score of -0.01 to -6.00; any point below 5% on the probability plot (with > 3 contiguous, and > 1 of the point below 1%); PSD at P < 0.05, or “outside normal limits” for the glaucoma hemifield test. Age-matched normal controls had open anterior chamber angle (Shaffer grade II or above), untreated IOP of < 21 mmHg, no abnormalities in the anterior segment and posterior segment on clinical examination, no structural or functional evidence of glaucoma, and no family history of glaucoma or other ocular diseases except visually insignificant cataract or myopia/hyperopia of less than 3 diopters (D). The latter is a measure to avoid inclusion of patients with extreme AL that may affect the image scale, hence the quantitative metrics as well as area over which measurements were made.

Exclusion criteria included age of < 18 years old, BCVA of worse than 20/40, glaucoma stage of moderate or worse according to the HAP staging,20 unreliable VF examination (fixation loss > 20%, false positive > 15%, or false negative > 15%, or other evidence of poor quality including inattention fixation, eyelid or lens rim artifacts, fatigue effects, and abnormal results caused by other factors other than glaucoma), OCT images with a signal strength of less than 20, low quality OCT-A imaging (see the session “OCT-A Quality Control”), evidence of secondary causes of glaucoma (e.g. steroid-induced glaucoma, history of uveitis, angle recession resulting from previous trauma), history of ocular surgery except for uneventful cataract surgery and/or laser iridotomy, and eyes with other retinal diseases that might affect the measurement of the microvasculature (e.g., diabetic retinopathy, epiretinal membrane, or age-related macular degeneration). We also excluded patients who had acute angle closure attack.

OCT-A Imaging

All study subjects underwent OCT-A using the swept-source OCT (Triton DRI-OCT, Topcon, Tokyo, Japan). Volumetric OCT scans centered on the optic disc were obtained with a scan area of 3 × 3 mm containing 320 × 320 A-scans. The latest built-in software (IMAGEnet6) was used to generate OCT angiograms, which provides improved detection sensitivity of low blood flow and reduced motion artifacts without compromising axial resolution.21 An OCT-A quality score ranged from 0 to 100 was automatically given by the software for each volumetric OCT scan. The built-in software can separately detect four horizontal depth-resolved segments (the optic nerve head, vitreous, radial peripapillary capillary, and the choroid). The peripapillary capillaries were analyzed from the segment of the optic nerve head, extending from the internal limiting membrane (ILM) to the boundary of RNFL.

OCT-A Quality Control

Each OCT-A image and OCT cross-sectional B-scan images were evaluated in the CUHK Ocular Reading Center by a single reader (YMW), who was masked to the diagnosis and clinical demographics of the participants. OCT-A images with poor image quality or significant image artifacts were excluded before the quantitative analysis, including: (1) image quality score less than 40, (2) inaccurate segmentation of tissue layers or labs, (3) motion artifacts (e.g., vessel discontinuity), (4) blurry images, (5) poor centration or (6) signal loss (e.g., due to eye blinking).

Quantification of Retinal Microvasculature

All included OCT-A images were imported into a customized automated MATLAB program for quantitative analysis. The detailed process of the retinal microvasculature quantification and reliability assessments have been reported before.22 Quantitative vascular indices – including circumpapillary vessel density (cpVD) and fractal dimension (cpFD) – were obtained after large retinal vessels removal. cpVD was calculated as the percentage of area and was defined as perfusion regions over the total area within the circumpapillary region in the binarized image.23 cpFD is a mathematical index that quantifies complex geometric patterns in structures that are self-similar in their scaling patterns, and is considered as a classifier to discriminate the class of normal structures from abnormal or pathological structures.24 It represents a global measurement of the complexity of the vascular branching pattern;25, 26 a higher cpFD value indicates a more complex vascular branching pattern, whereas a lower cpFD value indicates a sparser vascular network. In our study, the MATLAB program automatically determined cpFD from the skeletonized image using the box-counting method.27 Fig. 1 demonstrates the quantitation of retinal microvasculature from OCT-A images using our customized program in normal, PACG, and NTG eyes.

Statistical Analysis

All statistical analyses were performed using the IBM SPSS statistics version 25.0 (SPSS, Inc., Chicago, IL, USA). Normality of the variables was assessed by histogram. A generalized estimating equation (GEE) was used to compare the parameters among three groups to adjust the inter-eye correlation. A chi-square test was used to compare the categorical variables among the three groups. Multivariable linear regression models adjusting for the age, gender, axial length, IOP and OCT-A quality score were performed to determine the associations of average RNFL thickness and visual field parameters for each standard deviation change with OCT-A metrics. The different linear regression models were compared to determine whether OCT-A vascular metrics reflect stronger associations with functional or structural parameters between the PACG and NTG group using Fisher’s Z-test. A P value of < 0.05 was considered statistically significant.

Results

Two hundred and eighty-nine eyes (90 eyes of PACG, 115 eyes of NTG, and 84 eyes of normal control) were originally selected in this study. Eighty-five eyes (22 eyes of PACG, 47 eyes of NTG and 16 eye of normal control) were excluded due to low quality score (5 eyes of PACG, 6 eyes of NTG), motion artifacts (8 eyes of PACG, 22 eyes of NTG, 10 eye of normal control), blurred images (2 eye of PACG, 8 eyes of NTG, 2 eye of normal control), signal loss (4 eye of PACG, 7 eyes of NTG, 2 eye of normal control), and poor centration (3 eyes of PACG, 4 eyes of NTG, 2 eye of normal control), leaving a total of 204 eyes (68 eyes from 51 PACG subjects, 68 eyes from 48 NTG subjects and 68 eyes from 49 normal subjects) for the final analysis.

Table 1 shows the demographics and clinical characteristics of the participants. There were no significant differences in gender, age, and CCT amongst the three groups. The average BCVA was > 0.8 (Snellen scale) for all the participants, with the PACG eyes having the lowest value amongst the groups (P < 0.001 versus control and P = 0.003 versus NTG eyes). The PACG eyes had shorter ACD, shorter AL, and a higher IOP compared with either NTG eyes (all P ≤ 0.001) or normal controls (all P ≤ 0.001). Eyes with glaucoma had lower MD, higher PSD, lower VFI, and used more glaucoma medications than the normal controls (all P < 0.001), without significant differences between the two groups of glaucoma.

Table 1

Demographics and ocular characteristics among PACG, NTG, and Controls

Variables

PACG

NTG

Controls

P value

PACG vs NTG

PACG vs Controls

NTG vs Controls

Subjects (n)

51

48

49

N/A

N/A

N/A

Eyes (n)

68

68

68

N/A

N/A

N/A

Gender (F/M)

35/16

27/21

33/16

0.218

0.891

0.277

Age at recruitment (years)

59.82 ± 0.91

62.42 ± 1.30

59.42 ± 1.15

0.100

0.783

0.083

SE (D)

0.14 ± 0.22

-1.35 ± 0.32

-1.20 ± 0.27

< 0.001

< 0.001

0.720

AL (mm)

22.43 ± 0.12

24.60 ± 0.17

24.26 ± 0.15

< 0.001

< 0.001

0.138

IOP (mmHg) (range)

16.76 ± 0.52

(7.5–24.0)

14.65 ± 0.38

(10.0–20.0)

14.07 ± 0.33

(10.0–21.0)

0.001

< 0.001

0.259

CCT (𝜇m)

542.97 ± 4.52

531.54 ± 4.18

539.06 ± 5.06

0.063

0.564

0.252

BCVA

0.84 ± 0.03

0.95 ± 0.02

1.02 ± 0.02

0.003

< 0.001

0.003

ACD (mm)

2.48 ± 0.06

3.22 ± 0.08

3.52 ± 0.31

< 0.001

0.001

0.357

SAP MD (dB) (range)

-1.44 ± 0.24

(-5.73–1.17)

-1.55 ± 0.20

(-4.92– 1.65)

0.53 ± 0.08

(-1.82–1.96)

0.745

< 0.001

< 0.001

SAP PSD (dB) (range)

2.69 ± 0.20

(1.04–9.67)

3.10 ± 0.23

(1.39–10.55)

1.54 ± 0.04

(0.89–2.62)

0.172

< 0.001

< 0.001

SAP VFI (%) (range)

96.54 ± 0.56

(81–100)

96.30 ± 0.40

(82–100)

99.49 ± 0.08

(97–100)

0.721

< 0.001

< 0.001

Number of glaucoma medications

1.14 ± 0.13

0.82 ± 0.12

0

0.066

< 0.001

< 0.001

Abbreviations: CCT = central corneal thickness; BCVA = best-corrected visual acuity; ACD = anterior chamber depth; SE = spherical equivalent; AL = axial length; IOP = intraocular pressure; SAP = standard automated perimetry; MD = mean deviation; PSD = pattern standard deviation; VFI = visual field index; PACG = primary angle closure glaucoma; NTG = normal tension glaucoma.
*Statistical significance was tested by generalized estimating equation (for continuous variables) or chi-square test (for categorical variables).

Table 2 shows the comparison of RNFL thickness and OCT-A metrics among the three groups. Whilst eyes in either glaucoma group have thinner RNFL thickness than the control eye (P < 0.001), there was no significant difference between the average RNFL thickness of PACG eyes and NTG eyes. NTG eyes had thinner RNFL thickness in the superotemporal and superonasal regions compared with PACG eyes (all P ≤ 0.038). Compared with normal eyes, eyes in either glaucoma group also showed lower global cpVD (both P ≤ 0.014) and lower cpFD (both P < 0.001). Comparing NTG eyes with PACG eyes, NTG eyes had a significantly lower global cpVD (52.369 ± 0.781% vs 55.389 ± 0.721%), lower VD at the superotemporal (53.410 ± 1.047% vs 57.094 ± 1.106%), inferotemporal (53.971 ± 1.129% vs 58.458 ± 0.909%), inferonasal (48.879 ± 1.130% vs 54.946 ± 1.027%) and superonasal (50.457 ± 1.169% vs 53.365 ± 0.903%) regions (all P ≤ 0.049).

Table 2

Comparison of RNFL thickness and OCT-A metrics among PACG, NTG, and Controls

Variables

PACG

NTG

Controls

P value

PACG vs NTG

PACG vs Controls

NTG vs Controls

RNFL thickness (µm)

           

Average

88.00 ± 1.89

83.34 ± 1.62

102.45 ± 1.44

0.061

< 0.001

< 0.001

Superotemporal

120.93 ± 3.33

111.67 ± 2.97

140.79 ± 3.03

0.038

< 0.001

< 0.001

Temporal

69.48 ± 1.70

73.41 ± 1.82

86.31 ± 3.21

0.116

< 0.001

< 0.001

Inferotemporal

121.09 ± 3.98

112.99 ± 4.24

146.79 ± 4.32

0.163

< 0.001

< 0.001

Inferonasal

97.82 ± 2.99

92.48 ± 2.50

115.27 ± 3.08

0.171

< 0.001

< 0.001

Nasal

62.97 ± 1.90

59.92 ± 2.18

70.92 ± 2.72

0.287

0.016

0.001

Superonasal

98.21 ± 3.29

84.73 ± 3.28

111.24 ± 3.59

0.004

0.007

< 0.001

cpVD (%)

           

Global

55.389 ± 0.721

52.369 ± 0.781

57.700 ± 0.599

0.004

0.014

< 0.001

Superotemporal

57.094 ± 1.106

53.410 ± 1.047

58.686 ± 0.828

0.016

0.249

< 0.001

Temporal

56.939 ± 0.977

56.920 ± 0.924

59.956 ± 0.755

0.989

0.015

0.011

Inferotemporal

58.458 ± 0.909

53.971 ± 1.129

60.817 ± 0.951

0.002

0.073

< 0.001

Inferonasal

54.946 ± 1.027

48.789 ± 1.130

56.658 ± 1.005

< 0.001

0.234

< 0.001

Nasal

51.643 ± 0.896

50.483 ± 1.130

54.547 ± 0.786

0.421

0.015

0.003

Superonasal

53.365 ± 0.903

50.457 ± 1.169

55.309 ± 0.768

0.049

0.101

0.001

cpFD

1.519 ± 0.004

1.516 ± 0.004

1.548 ± 0.003

0.523

< 0.001

< 0.001

Abbreviations: cpRNFL = circumpapillary retinal nerve fiber layer; cpVD = circumpapillary vessel density; cpFD = circumpapillary fractal dimension; PACG = primary angle closure glaucoma; NTG = normal tension glaucoma.
*Statistical significance was tested by generalized estimating equation.

Table 3 shows the relationships of OCT-A metrics with average RNFL thickness in PACG, NTG, and normal controls. For the PACG and NTG group, after adjusting for age, gender, AL, IOP, number of glaucoma medications, and OCT-A image quality score, multivariable analyses showed that decreased cpVD ([PACG] β = -4.242; CI: -8.120, -0.363, and [NTG] β = -5.531; CI: -9.472, -1.590), and cpFD ([PACG] β = -8.894; CI: -11.925, -5.864, and [NTG]: β = -12.064; CI: -17.195, -6.932) were associated with decreased RNFL thickness (all P ≤ 0.032). In contrast, there was no significant association between OCT-A metrics and average RNFL thickness in normal eyes. The relationships of OCT-A metrics with the MD and VFI are shown in Table 4. In the multivariable analyses, decreased cpVD was significantly associated with decreased MD in NTG eyes (β = -0.707; CI: -1.090, -0.324; P < 0.001). However, in PACG eyes, there was no significant association between the OCT-A metrics and VF parameters. We also compared the strength of associations of OCT-A metrics with average RNFL thickness, MD, and VFI between PACG and NTG (Table 5). The Fisher’s Z-test revealed that there was a significantly stronger association between cpFD and average RNFL thickness in the NTG group compared with that in the PACG group (P = 0.045).

Table 3

Relationship of OCT-A metrics with average RNFL thickness in PACG, NTG, and Controls

OCT-A metrics

 

Univariable model

 

Multivariable model

 
 

Average RNFL thickness, µm (95%CI)

P value

Average RNFL thickness, µm (95%CI)

P value

PACG

         

cpVD (%)

per SD decrease

-5.855 (-9.567, -2.143)

0.002

-4.242 (-8.120, -0.363)

0.032

cpFD

per SD decrease

-9.297 (-11.682, -6.913)

< 0.001

-8.894 (-11.925, -5.864)

< 0.001

NTG

         

cpVD (%)

per SD decrease

-4.998 (-7.674, -2.322)

< 0.001

-5.531 (-9.472, -1.590)

0.006

cpFD

per SD decrease

-9.831 (-15.901, -3.761)

0.002

-12.064 (-17.195, -6.932)

< 0.001

Controls

         

cpVD (%)

per SD decrease

1.918 (-0.684, 4.520)

0.148

2.221 (-0.197, 4.638)

0.072

cpFD

per SD decrease

-0.646 (-3.284, 1.992)

0.631

-1.325 (-3.891, 1.241)

0.312

Abbreviations: cpVD = circumpapillary vessel density; cpFD = circumpapillary fractal dimension; PACG = primary angle closure glaucoma; NTG = normal tension glaucoma; RNFL = retinal nerve fiber layer; CI = confidence interval; SD = standard deviation. Multivariable model adjusted for age, gender, spherical equivalent, intraocular pressure, number of glaucoma medications and OCT-A image quality score.

Table 4

Relationship of OCT-A metrics with visual field mean deviation and visual field index in PACG and NTG

OCT-A metrics

 

Univariable model

Multivariable model

 

SAP MD, dB (95%CI)

P value

SAP VFI, % (95%CI)

P value

SAP MD, dB (95%CI)

P value

SAP VFI, % (95%CI)

P value

PACG

                 

cpVD (%)

per SD decrease

-0.096

(-0.576, 0.385)

0.697

-0.714

(-1.715, 0.286)

0.162

-0.073

(-0.601, 0.455)

0.786

-1.304

(-3.070, 0.463)

0.148

cpFD

per SD decrease

-0.345

(-0.914, 0.225)

0.236

-2.126

(-4.002, -0.250)

0.026

-0.49

(-1.713, 0.733)

0.433

-1.874

(-4.678, 0.931)

0.190

NTG

                 

cpVD (%)

per SD decrease

-0.253

(-0.508, 0.002)

0.052

-0.334

(-0.724, 0.056)

0.093

-0.707

(-1.090, -0.324)

< 0.001

-0.338

(-1.100, 0.424)

0.385

cpFD

per SD decrease

-0.221

(-0.600, 0.157)

0.252

-0.867

(-2.139, 0.406)

0.182

-0.032

(-1.106, 1.042)

0.953

-0.424

(-1.621, 0.774)

0.488

Controls

                 

cpVD (%)

per SD decrease

0.01

(-0.162, 0.181)

0.913

0.036

(-0.151, 0.223)

0.706

-0.002

(-0.172, 0.167)

0.977

0.015

(-0.173, 0.202)

0.880

cpFD

per SD decrease

0.126

(-0.076, 0.329)

0.222

-0.098

(-0.331, 0.135)

0.409

0.132

(-0.088, 0.351)

0.239

-0.168

(-0.405, 0.069)

0.165

Abbreviations: cpVD = circumpapillary vessel density; cpFD = circumpapillary fractal dimension; PACG = primary angle closure glaucoma; NTG = normal tension glaucoma; SAP = standard automated perimetry; MD = mean deviation; PSD = pattern standard deviation; VFI = visual field index; CI = confidence interval; SD = standard deviation. Multivariable model adjusted for age, gender, spherical equivalent, intraocular pressure, number of glaucoma medications and OCT-A image quality score.

Table 5

Comparison of linear regression analysis of OCT-A metrics with structural and functional alterations between PACG and NTG

OCT-A metrics

Average RNFL thickness, µm

 

β (95% CI)

R2

cpVD (%)

 

PACG

-4.242 (-8.120, -0.363) *

0.388

NTG

-5.531 (-9.472, -1.590) *

0.341

P value†

0.592

 

cpFD

   

PACG

-8.894 (-11.925, -5.864) **

0.548

NTG

-12.064 (-17.195, -6.932) **

0.342

P value†

0.045

 
*denotes P < 0.05; **denotes P < 0.001; † Fisher's Z test for comparison of linear regression models. OCT-A metrics were analyzed as per SD decrease. Results were adjusted for age, gender, spherical equivalent, intraocular pressure, number of glaucoma medications and OCT-A image quality score.
Abbreviations: cpVD = circumpapillary vessel density; cpFD = circumpapillary fractal dimension; PACG = primary angle closure glaucoma; NTG = normal tension glaucoma; RNFL = retinal nerve fiber layer; β = regression coefficient; CI = confidence interval; MD = mean deviation; VFI = visual field index; R2 = adjusted coefficient of determination.

Discussion

Our findings of reduced cpVD and cpFD of glaucomatous eyes compared with age-matched normal controls were in concordance with previous studies that showed reduction of the microvascular density of the optic disc in glaucoma patients.11, 2831 Rao et al29 reported that VD in optic disc region in PACG eyes was significantly lower than control eyes, whereas VD in primary angle closure (PAC) with high IOP and thinner superotemporal peripapillary RNFL thickness was similar to that of the controls. They suggested that high IOP affects the RNFL measurements earlier than VD in PACG. Scripsema et al11 found that both POAG and NTG patients demonstrated decreased perfused capillary density compared to normal subjects, with POAG patients having a lower perfused capillary density than NTG patients. They attributed the latter finding to be possibly related to medication effect (POAG patients using larger number of drops) and/or the different pathophysiological processes of NTG and HTG.

To our knowledge, this is the first study to compare the peripapillary microvasculature in early PACG and early NTG eyes using OCT-A. We utilized OCT-A to compare the peripapillary microvasculature of two subtypes of early glaucoma which have likely different pathogenic mechanisms – NTG that is less IOP dependent and has possibly a stronger vascular pathogenic component,5, 32 and PACG that has likely a predominantly IOP-dependent mechanism.33 This allowed a “snapshot” comparison of the detailed microvasculature between NTG and PACG at the early stage of the disease. Our results showed a significantly reduced global cpVD in NTG eyes compared with PACG eyes, despite the comparable RNFL thickness and disease severity. We also identified a significant association between the OCT-A metrics and RNFL thickness in both glaucoma groups, with a stronger relationship between the cpFD and RNFL thickness in the NTG group compared with the PACG group. The findings might reflect a more specific and early reduction of microvascular perfusion at the peripapilliary region of NTG eyes at early stage glaucoma. Such a difference was not observed in other studies that compared patients with PACG and POAG (involved both NTG and HTG).34, 35

Given the similar disease severity and the number of medications used, the difference in cpVD between NTG and PACG eyes could be related to the different pathogenic mechanisms of the two glaucoma subtypes. For NTG at the early stage, the impairment of vascular autoregulation – postulated to be an important risk factor for disease progression in NTG5, 6, 3638 – could lead to reduced blood flow but has yet to cause dysfunction, death, or atrophy of the RGC, as well as the consequential RNFL thinning. Therefore, in NTG eyes, there was a reduction of peripapillary microvascular perfusion prior to RGCs loss, and a further delay for the development of a detectable RNFL thinning after the RGCs’ loss. The latter is supported by in vivo study that showed an initially faster decline of RGC soma counts compared with RNFL thickness following optic nerve injury in animal model.39 Our findings of the stronger association between cpFD and RNFL thickness in NTG eyes compared with PACG eyes, as well as the association between decrease in cpVD and decrease MD of VF in NTG (which was not observed in PACG eyes), may provide a modest but concordant support of this theory. Further in vivo studies are needed to validate the role of microvasculature in the pathophysiology in NTG. In PACG eyes, the loss of RGCs was possibly mainly due to elevation of IOP that occurred prior to IOP lowering treatments (including lens extraction and/or laser iridotomy) and the reduction of the OCT-A metrics could be a secondary consequence of RGCs loss; this echoed with the suggestion of Rao et al. that high IOP affects the RNFL measurements earlier than VD in PACG.29 The atrophy of RGCs may lead to a reduced demand of blood supply and blood flow that was reflected as a reduction in cpVD. However, this secondary reduction of cpVD, unlike NTG eyes at a similar stage of disease severity, was not extensive enough to a degree that would cause further RGCs loss. Hence, this possible difference in the pathogenic algorithm of NTG and PACG might lead to the differences in cpVD measurement between the two subgroups, despite similar RNFL thickness. A longitudinal study with larger number of patients may, in the future, verify the causal relationship of vascular-RNFL thickness in these two subtypes of glaucoma. This is not only important in terms of understanding the pathophysiology of the disease, it is also clinically implicative if OCT-A is to be utilized as a diagnostic and monitoring tool in glaucoma management.

In the multivariable analyses, we found a significant association between OCT-A metrics and RNFL thickness in both glaucoma subtypes but only an association between cpVD and MD in the NTG group (Table 4). This differed from previous studies that reported a stronger association between decreased cpVD with the severity of VF damage, compared with the association between RNFL thickness and VF function in PACG34 and POAG eyes.14, 34 However, our results were expected because we only included glaucoma patients with mild severity with the MD score of better than − 6.0 D. Identifying an association between OCT-A metrics within a narrow range of VF parameters (MD of -5.73 to 1.65 dB and VFI of 81 to 100% in the glaucoma groups) is understandably difficult. Furthermore, the association between microvasculature with functional change may not be strong in these patients with early stage disease that have minimal functional loss. Indeed, a study by Shin et al showed that whilst there was a significant relationship between VD and VF function in moderate-to-advance POAG regardless of location, the relationship of VD and VF function was only significant in the superotemporal and inferotemporal regions for early stage POAG40.

The strength of our study was the inclusion of early NTG and PACG patients – that have similar age, disease severity, RNFL thickness and number of medications used – for comparison. This could better reflect the role of microvasculature in the glaucomatous pathogenic mechanisms at the early stage of the disease. We also used an objective, automated MATLAB program to quantitatively measure retinal microvasculature. Limitations of this study included a cross-sectional study design, a relatively small sample size, as well as not taking systemic diseases into consideration (e.g. obstructive sleep apnea, hypertension). We have limited the range of refractive errors (+ 3.0 to -3.0 D) in an attempt to avoid inclusion of eyes with extreme AL. We acknowledge that the AL of PACG eyes were statistically shorter than either the NTG eyes or control eyes in the current study (Table 1). However, the difference in AL between the PACG and NTG group was reasonable in the clinical point of view (22.43 ± 0.12 D vs 24.60 ± 0.17 D; P < 0.001). Nonetheless, our findings provide the basis for future longitudinal study that may review the causal relationship between retinal microvasculature change and RNFL thickness change.

In summary, the cpVD was significantly lower in early NTG eyes when compared to early PACG eyes, despite similar RNFL thickness and VF parameters. Reductions in cpVD and cpFD were associated with average RNFL thickness thinning in both NTG and PACG eyes. Longitudinal study may verify the differences of microvasculature change in different glaucoma subtypes, improve our understanding of the pathogenic mechanisms of the diseases, and also establish a role for OCT-A in the management of glaucoma.

Declarations

Grant support:

Health and Medical Research Fund, Hong Kong (Ref. No.05162836 to C.C.T.). General Research Fund, Hong Kong (Ref. No. 14107516 to C.C.T.). The funding organization had no role in the design or conduct of this research.

Conflict of interest:

nil

References

  1. Weinreb RN, Aung T, Medeiros FA. The pathophysiology and treatment of glaucoma: a review. JAMA 2014;311:1901–11.
  2. The effectiveness of intraocular pressure reduction in the treatment of normal-tension glaucoma. Collaborative Normal-Tension Glaucoma Study Group. Am J Ophthalmol 1998;126:498–505.
  3. Leske MC, Wu SY, Hennis A, Honkanen R, Nemesure B, Group BES. Risk factors for incident open-angle glaucoma: the Barbados Eye Studies. Ophthalmology 2008;115:85–93.
  4. Bonomi L, Marchini G, Marraffa M, Bernardi P, Morbio R, Varotto A. Vascular risk factors for primary open angle glaucoma - The Egna-Neumarkt study. Ophthalmology 2000;107:1287–1293.
  5. Flammer J, Orgul S, Costa VP, et al. The impact of ocular blood flow in glaucoma. Prog Retin Eye Res 2002;21:359–93.
  6. Harris A, Rechtman E, Siesky B, Jonescu-Cuypers C, McCranor L, Garzozi HJ. The role of optic nerve blood flow in the pathogenesis of glaucoma. Ophthalmol Clin North Am 2005;18:345 – 53, v.
  7. Wu R, Cheung CY, Saw SM, Mitchell P, Aung T, Wong TY. Retinal vascular geometry and glaucoma: the Singapore Malay Eye Study. Ophthalmology 2013;120:77–83.
  8. Kawasaki R, Wang JJ, Rochtchina E, Lee AJ, Wong TY, Mitchell P. Retinal vessel caliber is associated with the 10-year incidence of glaucoma: the Blue Mountains Eye Study. Ophthalmology 2013;120:84–90.
  9. Lee JY, Yoo C, Park JH, Kim YY. Retinal vessel diameter in young patients with open-angle glaucoma: comparison between high-tension and normal-tension glaucoma. Acta Ophthalmol 2012;90:e570-1.
  10. Ritch R, Lowe RF. Angle closure glaucoma: mechanisms and epidemiology. In: Ritch R, Shields MB, Krupin T, editors. The Glaucomas. St Louis, Mo: Mosby, 1996:801–819.
  11. Scripsema NK, Garcia PM, Bavier RD, et al. Optical Coherence Tomography Angiography Analysis of Perfused Peripapillary Capillaries in Primary Open-Angle Glaucoma and Normal-Tension Glaucoma. Invest Ophthalmol Vis Sci 2016;57:OCT611-OCT620.
  12. Jia Y, Wei E, Wang X, et al. Optical coherence tomography angiography of optic disc perfusion in glaucoma. Ophthalmology 2014;121:1322–32.
  13. Liu L, Jia Y, Takusagawa HL, et al. Optical Coherence Tomography Angiography of the Peripapillary Retina in Glaucoma. JAMA Ophthalmol 2015;133:1045–52.
  14. Yarmohammadi A, Zangwill LM, Diniz-Filho A, et al. Relationship between Optical Coherence Tomography Angiography Vessel Density and Severity of Visual Field Loss in Glaucoma. Ophthalmology 2016;123:2498–2508.
  15. Akagi T, Iida Y, Nakanishi H, et al. Microvascular Density in Glaucomatous Eyes With Hemifield Visual Field Defects: An Optical Coherence Tomography Angiography Study. Am J Ophthalmol 2016;168:237–249.
  16. Wang X, Jiang C, Ko T, et al. Correlation between optic disc perfusion and glaucomatous severity in patients with open-angle glaucoma: an optical coherence tomography angiography study. Graefes Arch Clin Exp Ophthalmol 2015;253:1557–64.
  17. Cheung CY, Li SL, Chan N, et al. Factors Associated With Long-term Intraocular Pressure Fluctuation in Primary Angle Closure Disease: The CUHK PACG Longitudinal (CUPAL) Study. Journal of glaucoma 2018;27:703–710.
  18. Cheung CY, Li SL, Chan PP, et al. Intraocular pressure control and visual field changes in primary angle closure disease: the CUHK PACG Longitudinal (CUPAL) study. Br J Ophthalmol 2020;104:629–635.
  19. Tham CC, Kwong YY, Leung DY, et al. Phacoemulsification vs phacotrabeculectomy in chronic angle-closure glaucoma with cataract: complications [corrected]. Arch Ophthalmol 2010;128:303–11.
  20. Hodapp E, Parrish RK, Anderson DR. Clinical decision in glaucoma: St Louis: Mosby, 1993.
  21. Stanga PE, Tsamis E, Papayannis A, Stringa F, Cole T, Jalil A. Swept-Source Optical Coherence Tomography Angio (TM) (Topcon Corp, Japan): Technology Review. Dev Ophthalmol 2016;56:13–17.
  22. Tang FY, Ng DS, Lam A, et al. Determinants of Quantitative Optical Coherence Tomography Angiography Metrics in Patients with Diabetes. Sci Rep-Uk 2017;7.
  23. Kim AY, Chu ZD, Shahidzadeh A, Wang RKK, Puliafito CA, Kashani AH. Quantifying Microvascular Density and Morphology in Diabetic Retinopathy Using Spectral-Domain Optical Coherence Tomography Angiography. Invest Ophth Vis Sci 2016;57:Oct362-Oct370.
  24. Masters BR. Fractal analysis of the vascular tree in the human retina. Annu Rev Biomed Eng 2004;6:427–52.
  25. Cheung CY, Tay WT, Mitchell P, et al. Quantitative and qualitative retinal microvascular characteristics and blood pressure. J Hypertens 2011;29:1380–91.
  26. Liew G, Wang JJ, Cheung N, et al. The Retinal Vasculature as a Fractal: Methodology, Reliability, and Relationship to Blood Pressure. Ophthalmology 2008;115:1951–1956.
  27. Bao MH, Guo SY, Tang Q, Zhang F. Optimization of the bwmorph Function in the MATLAB Image Processing Toolbox for Binary Skeleton Computation. Proceedings of the 2009 International Conference on Computational Intelligence and Natural Computing, Vol Ii 2009:273-+.
  28. Zhang S, Wu C, Liu L, et al. Optical Coherence Tomography Angiography of the Peripapillary Retina in Primary Angle-Closure Glaucoma. Am J Ophthalmol 2017;182:194–200.
  29. Rao HL, Pradhan ZS, Weinreb RN, et al. Vessel Density and Structural Measurements of Optical Coherence Tomography in Primary Angle Closure and Primary Angle Closure Glaucoma. Am J Ophthalmol 2017;177:106–115.
  30. Rao HL, Kadambi SV, Weinreb RN, et al. Diagnostic ability of peripapillary vessel density measurements of optical coherence tomography angiography in primary open-angle and angle-closure glaucoma. Br J Ophthalmol 2017;101:1066–1070.
  31. Zhu L, Zong Y, Yu J, et al. Reduced Retinal Vessel Density in Primary Angle Closure Glaucoma: A Quantitative Study Using Optical Coherence Tomography Angiography. J Glaucoma 2018;27:322–327.
  32. Araie M, Sekine M, Suzuki Y, Koseki N. Factors contributing to the progression of visual field damage in eyes with normal-tension glaucoma. Ophthalmology 1994;101:1440–4.
  33. Gazzard G, Foster PJ, Devereux JG, et al. Intraocular pressure and visual field loss in primary angle closure and primary open angle glaucomas. Br J Ophthalmol 2003;87:720–5.
  34. Jo YH, Sung KR, Yun SC. The Relationship Between Peripapillary Vascular Density and Visual Field Sensitivity in Primary Open-Angle and Angle-Closure Glaucoma. Invest Ophthalmol Vis Sci 2018;59:5862–5867.
  35. Hou TY, Kuang TM, Ko YC, Chang YF, Liu CJ, Chen MJ. Optic Disc and Macular Vessel Density Measured by Optical Coherence Tomography Angiography in Open-Angle and Angle-Closure Glaucoma. Sci Rep 2020;10:5608.
  36. Chung HS, Harris A, Evans DW, Kagemann L, Garzozi HJ, Martin B. Vascular aspects in the pathophysiology of glaucomatous optic neuropathy. Surv Ophthalmol 1999;43 Suppl 1:S43-50.
  37. Grieshaber MC, Mozaffarieh M, Flammer J. What is the link between vascular dysregulation and glaucoma? Surv Ophthalmol 2007;52 Suppl 2:S144-54.
  38. Orgul S, Gugleta K, Flammer J. Physiology of perfusion as it relates to the optic nerve head. Surv Ophthalmol 1999;43 Suppl 1:S17-26.
  39. Munguba GC, Galeb S, Liu Y, et al. Nerve fiber layer thinning lags retinal ganglion cell density following crush axonopathy. Invest Ophthalmol Vis Sci 2014;55:6505–13.
  40. Shin JW, Lee J, Kwon J, Choi J, Kook MS. Regional vascular density-visual field sensitivity relationship in glaucoma according to disease severity. Br J Ophthalmol 2017;101:1666–1672.