# Optic Disc Vascularization in Glaucoma: Value of Spectral-Domain Optical Coherence Tomography Angiography

Pierre-Maxime Lévêque, Pierre Zéboulon, Emmanuelle Brasnu, Christophe Baudouin, Antoine Labbé

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Pierre-Maxime Lévêque, Pierre Zéboulon, Emmanuelle Brasnu, Christophe Baudouin, Antoine Labbé. Optic Disc Vascularization in Glaucoma: Value of Spectral-Domain Optical Coherence Tomography Angiography. Journal of Ophthalmology, 2016, 2016, pp.6956717. 10.1155/2016/6956717. hal-01294395

HAL Id: hal-01294395

https://hal.sorbonne-universite.fr/hal-01294395

Submitted on 29 Mar 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Research Article

# Optic Disc Vascularization in Glaucoma: Value of Spectral-Domain Optical Coherence Tomography Angiography

Pierre-Maxime Lévêque,¹ Pierre Zéboulon,¹ Emmanuelle Brasnu,¹,²
Christophe Baudouin,¹,²,³ and Antoine Labbé¹,²,³,⁴

¹CHNO des Quinze-Vingts, DHU Sight Restore, INSERM-DHOS CIC, 28 rue de Charenton, 75012 Paris, France

²Sorbonne Universités, UPMC, INSERM, CNRS, Institut de la Vision, 17 rue Moreau, 75012 Paris, France

³Department of Ophthalmology, Ambroise Paré Hospital (AP-HP), University of Versailles Saint-Quentin-en-Yvelines,
55 Avenue de Paris, 78000 Versailles, France

⁴Department of Ophthalmology III, Quinze-Vingts National Ophthalmology Hospital, 28 rue de Charenton, 75012 Paris, France

Correspondence should be addressed to Antoine Labbé; dr.antoinelabbe@gmail.com

Received 30 October 2015; Revised 21 January 2016; Accepted 24 January 2016

Academic Editor: Paolo Fogagnolo

Copyright © 2016 Pierre-Maxime Lévêque et al. This is an open access article distributed under the Creative Commons Attribution
License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

*Purpose.* To detect changes in optic nerve head (ONH) vascularization in glaucoma patients using spectral-domain OCT angiography (OCT-A). *Material and Method.* Fifty glaucoma patients and 30 normal subjects were evaluated with OCT-A (AngioVue®, Optovue). The total ONH vessel density and temporal disc vessel density were measured. Clinical data, visual field (VF) parameters, and spectral-domain OCT evaluation (RNFL: retinal nerve fiber layer thickness, GCC: ganglion cell complex thickness, and rim area) were recorded for glaucoma patients. Correlations among total and temporal ONH vessel density and structural and VF parameters were analyzed. *Results.* In the glaucoma group, total and temporal ONH vessel density were reduced by 24.7% (0.412 versus 0.547; *p* < 0.0001) and 22.88% (0.364 versus 0.472; *p* = 0.001), respectively, as compared with the control group. Univariate analysis showed significant correlation between rim area (mm²) and temporal ONH vessel density (*r* = 0.623; *p* < 0.0001) and total ONH vessel density (*r* = 0.609; *p* < 0.0001). Significant correlations were found between temporal and total ONH vessel density and RNFL, GCC, VF mean deviation, and visual field index. *Conclusion.* In glaucoma patients OCT-A might detect reduced ONH blood vessel density that is associated with structural and functional glaucomatous damage. OCT-A might become a useful tool for the evaluation of ONH microcirculation changes in glaucoma.

**1. Introduction**

Although elevated intraocular pressure (IOP) is the main
risk factor for glaucoma, high numbers of patients also
develop glaucoma with normal IOP levels [1]. Therefore,
other risk factors and in particular vascular risk factors have
been implicated in the pathogenesis of glaucoma [2]. These
vascular factors include autoregulation of the blood flow
(BF) in the optic nerve head (ONH), local vasospasm, arte-
rial hypertension, and nocturnal hypotension [3]. Impaired
microcirculation in the ONH may contribute to the initiation
and progression of glaucomatous neuropathy. BF in the ONH
is supplied by two arterial systems: the central retinal artery
in superficial layers and the posterior ciliary artery in deeper
layers (prelaminar and lamina cribrosa and retrolaminar

structures) [4]. It has been proposed that the main pathologic
changes in glaucoma are located in the deep ONH region,
supplied by ciliary arteries [5].

Consequently, several techniques have been proposed
to measure optic disc perfusion [6, 7]. Laser Doppler
velocimetry has been used to identify the maximum BF
velocity present in large retinal vessels [6]. Other methods
such as scanning laser Doppler flowmetry [8] and laser
speckle flowgraphy [9] were also developed to measure disc
perfusion. Doppler optical coherence tomography (OCT) has
been used to measure the total retinal blood flow around
the ONH [10]. Using this imaging technique, BF in large
vessels of the ONH can be quantified, but the study of the
microvascular structure is impossible due to the low velocity
of the BF in small vessels. The microsphere method, a method
based on radioactive solid microspheres injected into the systemic circulation, has also provided evidence that the laser speckle flowgraphy is capable of assaying BF in the deep ONH region but only in experimental glaucoma models [11].

More recently, Jia et al. developed a new method to study in vivo the ONH microcirculation in glaucoma patients: the split spectrum amplitude-decorrelation angiography (SSADA) algorithm [12]. It provides high-quality three-dimensional (3D) angiography using ultra-high-speed swept-source OCT [13]. These authors found reduced disc perfusion in a group of patients with early glaucoma and established a link between ONH vessel density and visual field pattern standard deviation (PSD) using a custom swept-source OCT device [14]. The OCT angiography with SSADA may be more reliable than LSFG or LDF, offering better intravisit repeatability and intervisit reproducibility [14–16].

The objective of the present study was to detect changes in ONH vascularization by comparing normal subjects and glaucoma patients using a commercial spectral-domain OCT angiography (OCT-A) device.

## 2. Methods

2.1. *Study Population.* This cross-sectional observational study was conducted at the Quinze-Vingts National Ophthalmology Hospital, Paris, France, from January to September 2015 after approval of our Institutional Review Board. The study adhered to the tenets of the Declaration of Helsinki and informed consent was obtained from all subjects. Patients followed up at the Quinze-Vingts National Ophthalmology Hospital for glaucoma were consecutively included. Inclusion criteria were chronic glaucoma defined as glaucomatous optic disc neuropathy and characteristic visual field (VF) loss based on the criteria of the Ocular Hypertension Treatment Study [17]. The glaucomatous VF loss was defined as a glaucoma hemifield test graded “outside normal limits” and a cluster of three contiguous points at the 5% level on the pattern deviation plot, using the threshold test strategy with the 24-2 test pattern of the Humphrey Field Analyzer II [18]. Patients with diabetic retinopathy, other diseases that may cause visual field loss or optic disc abnormalities, and inability to clinically view or photograph the optic discs due to media opacity or poorly dilating pupil were excluded.

For all glaucoma patients, the following data were recorded: age, gender, best corrected visual acuity (BCVA) on a logarithmic scale, IOP (mmHg), central corneal thickness (CCT), type of glaucoma, and cup/disc ratio. IOP was measured using the Goldmann applanation tonometer. For each eye, IOP value represented the mean value of two IOP readings. CCT was measured with an A-scan ultrasound (A-Scan Pachymeter, Ultrasonic, Exton, PA, USA). The mean retinal nerve fiber layer (RNFL) thickness, mean ganglion cell complex (GCC) thickness, and disc rim area were also evaluated (Cirrus® spectral-domain OCT, Carl Zeiss Meditec Inc., Dublin, IE). VF tests were performed with a Humphrey Field Analyzer II (Carl Zeiss Meditec Inc.) set for 24-2 to collect mean deviation (MD), pattern standard deviation (PSD), and the visual field index (VFI). The Bascom Palmer (Hodapp-Anderson-Parrish) glaucoma staging system (GSS)

based on VF results was selected to classify patients as having mild, moderate, and severe glaucoma [19].

Healthy age-matched controls were also included. Normal subjects had IOP < 21 mmHg with clinically normal ONH and no history of ocular or systemic disease. For control subjects, demographic data and IOP were recorded.

2.2. *Optical Coherence Tomography.* We used the commercially available spectral-domain OCT RT XR Avanti with the AngioVue software (Optovue, Inc., Fremont, CA, USA). It detects BF in an acquired volume using the SSADA [12]. The instrument used for OCT-A images is based on the AngioVue Imaging System (Optovue, Inc., Fremont, CA, USA) to obtain amplitude decorrelation angiography images. This instrument has an A-scan rate of 70,000 scans per second, using a light source centered on 840 nm and a bandwidth of 50 nm. Each OCT-A volume contains 304 x 304 A-scans with two consecutive B-scans captured at each fixed position before proceeding to the next sampling location. Split spectrum amplitude-decorrelation angiography was used to extract the OCT angiography information. Each OCT-A volume is acquired in three seconds and two orthogonal OCT-A volumes were acquired in order to perform motion correction to minimize motion artifacts arising from microsaccades and fixation changes.

The software analyzes the amplitude of variation of the OCT signal over time for every location acquired and calculates decorrelation. Static tissue yields low decorrelation values as the signal varies a little over time, and BF yields high decorrelation values because moving red blood cells cross the OCT beam and cause the signal amplitude to vary rapidly over time. A threshold decorrelation value is therefore used to discriminate BF from static tissue [20, 21]. The angiography information displayed was the average of the decorrelation values when viewed perpendicularly through the thickness evaluated. This maps all the detected vessels in the acquired volume. The AngioVue software (Optovue) can then calculate a quantitative variable: the total surface of the flowing vessel on the en face projection. The right eye from each normal subject was scanned. For glaucoma patients, we decided to analyze the most severe eye as determined by the MD 24-2.

2.3. *Image Acquisition and Processing.* First, each subject underwent one set of two scans on the macular region. The SSADA algorithm requires the user to define a threshold of detection of circulating vessels [13]. Jia et al. described a method to calculate the appropriate threshold. The central avascular zone was defined as a noise region and the mean decorrelation value of the whole avascular zone was calculated. A threshold one standard deviation above the mean avascular decorrelation value was defined as circulating vessels.

Then we performed one set of two 3 x 3 mm scans on the ONH. We used the neural canal opening, which is the termination of the retinal pigment epithelium/Bruch membrane complex, to define the ONH margin [22]. The volumes analyzed were defined by both a circular surface on the en face projection and the anteroposterior boundaries, thus creating a cylinder. A circle on the en face projection
<figure><img src="image_3.png" /><figcaption>FIGURE 1: Total (a) and temporal (b) ONH acquisition in a normal patient. Total (c) and temporal (d) ONH acquisition in a glaucoma patient.</figcaption></figure>

and anteroposterior limits defines each cylinder. We studied the total ONH flow area (mm²), representing the number of vessels, detected from all the depths of the disc, from the top of the ONH (inner surface) to the lamina cribrosa. Next, the total ONH vessel density (%) was calculated by dividing flow area by the optic disc surface (πR²). The ONH temporal region (determined by the temporal disc between the temporal arteries and veins) was also studied to avoid BF due to major superior and inferior branches of retinal vessels in order to focus on the ONH laminar microvasculature (Figure 1). Then we measured the temporal disc flow area (mm²) and calculated the temporal disc vessel density by dividing the temporal disc flow area by the temporal disc surface (πR² temporal). The total and temporal optic disc vessel densities, obtained as the ratio of the optic disc flow area to the optic disc surface, were calculated to account for the intereye differences in total and temporal optic disc surface areas. Image quality was determined for all OCT-A scans. Poor quality images with a signal strength index (SSI) < 40 or image sets with remaining motion artifacts were excluded from the study analysis [23].

2.4. *Reproducibility and Repeatability.* Measurement repeatability of the disc flow area was evaluated from three sets of acquisitions from a subset of 6 normal subjects within a single visit. The intravisit coefficient of variation (CV) was calculated by comparing these three measurements. Intervisit reproducibility was obtained from three sets of scans performed on three separate visits from the same subset of normal subjects.

2.5. *Statistical Analysis.* Mann-Whitney tests were used to compare the average values of the measurements between normal and glaucomatous eyes. Linear regression analysis was used in the control group to investigate whether the measurement of ONH vessel density was affected by age, gender, and IOP. Univariate analysis with the Pearson correlation test and multivariate analysis with the ANCOVA test were performed to study the correlation between ONH vessel density and other variables: age, RNFL thickness, GCC, rim area, MD 24-2, and PSD 24-2.

**3. Results**

**3.1. Demographic and Clinical Data.** Fifty glaucoma patients and 30 normal subjects were included in this study. There was
<table><caption>TABLE 1: Qualitative and quantitative demographic data. IOP: intraocular pressure; F: female, M: male.</caption><thead><tr><th></th><th colspan="2">Normal</th><th colspan="2">Glaucoma</th><th rowspan="2">p value</th></tr><tr><th></th><th>(n = 30)</th><th>%</th><th>(n = 50)</th><th>%</th></tr></thead><tbody><tr><td><strong>Qualitative demographic data</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Gender</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>F</td><td>17</td><td>56.7</td><td>23</td><td>46</td><td>NS</td></tr><tr><td>M</td><td>13</td><td>43.3</td><td>27</td><td>54</td><td></td></tr><tr><td>Family history of glaucoma</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>0</td><td>27</td><td>90</td><td>24</td><td>48</td><td>p &lt; 0.0001</td></tr><tr><td>1</td><td>3</td><td>10</td><td>26</td><td>52</td><td></td></tr><tr><td>Ethnic origin</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>African</td><td>8</td><td>26.7</td><td>11</td><td>22</td><td>—</td></tr><tr><td>Caucasian</td><td>20</td><td>66.7</td><td>33</td><td>66</td><td></td></tr><tr><td>Asian</td><td>0</td><td>0</td><td>2</td><td>4</td><td></td></tr><tr><td>Hispanic</td><td>0</td><td>0</td><td>2</td><td>4</td><td></td></tr><tr><td>Maghrebi</td><td>2</td><td>6.7</td><td>2</td><td>4</td><td></td></tr><tr><td><strong>Quantitative demographic data</strong></td><td>Mean</td><td>SD</td><td>Mean</td><td>SD</td><td></td></tr><tr><td>Age (years)</td><td>51.1</td><td>18.1</td><td>58.4</td><td>15.9</td><td>NS</td></tr><tr><td>IOP (mmHg)</td><td>14.4</td><td>2.59</td><td>17.6</td><td>6.41</td><td>0.022</td></tr></tbody></table>

no difference between the two groups concerning age, gender, and ethnic origin. A family history of glaucoma was found in 52% of glaucoma patients and 10% of control subjects (p < 0.0001). IOP was 17.6 mmHg and 14.4 mmHg in the glaucoma and control group, respectively (p = 0.022). The demographic data are summarized in Table 1.

Regarding the type of glaucoma, 43 patients (86%) had open-angle glaucoma (OAG) and seven patients (14%) had angle-closure glaucoma (ACG). Twenty-three patients (46%) had mild glaucoma, eight (16%) had moderate glaucoma, and 19 (38%) had severe glaucoma. Mean RNFL thickness, GCC thickness, and rim area were 66.48 ± 12.6 µm, 63.6 ± 9.07 µm, and 0.786 ± 0.291 mm² in the glaucoma group, respectively. The results of the glaucoma patients' clinical data are summarized in Table 2.

### 3.2. ONH Blood Flow Evaluation.
The en face projection angiograms showed that normal discs had a denser microvascular network, especially in the temporal region, compared with glaucomatous discs (Figure 1). Total ONH vessel density was reduced by 24.7% (0.412 ± 0.117 versus 0.547 ± 0.077, p < 0.0001) in the glaucoma group as compared with the control group. Temporal disc vessel density was significantly reduced by 22.88% (0.364 ± 0.150 versus 0.472 ± 0.105, p = 0.001) in the glaucoma group as compared with the control group. Moreover, SSI and the decorrelation threshold were higher in the control group (72.3 versus 61.8, p < 0.0001, and 0.095 versus 0.085, p = 0.001, resp.). The intravisit CV was 6.17%. The intervisit CV was 6.48%. The results of the ONH BF evaluation are summarized in Table 3.

In the glaucoma group, univariate analysis showed that total ONH vessel density was significantly associated with structural variables including rim area (r = 0.609; p < 0.0001), RNFL thickness (r = 0.406; p < 0.003), and GCC thickness (r = 0.337; p = 0.017) (Figure 2). Total ONH vessel density was also associated with the functional variables VFI

<table><caption>TABLE 2: Quantitative and qualitative clinical data. OAG: open-angle glaucoma; ACG: angle-closure glaucoma. Visual field and OCT spectral-domain testing. MD: mean deviation; PSD: pattern standard deviation; VFI: visual field index; RNFL: retinal nerve fiber layers; GCC: ganglion cell complex.</caption><thead><tr><th></th><th colspan="2">Glaucoma</th></tr><tr><th></th><th>(n = 50)</th><th>%</th></tr></thead><tbody><tr><td><strong>Qualitative clinical data</strong></td><td></td><td></td></tr><tr><td>Type</td><td></td><td></td></tr><tr><td>OAG</td><td>43</td><td>86</td></tr><tr><td>ACG</td><td>7</td><td>14</td></tr><tr><td>Stage</td><td></td><td></td></tr><tr><td>Mild</td><td>23</td><td>46</td></tr><tr><td>Moderate</td><td>8</td><td>16</td></tr><tr><td>Severe</td><td>19</td><td>38</td></tr><tr><td><strong>Quantitative clinical data</strong></td><td>Mean</td><td>SD</td></tr><tr><td>log AV</td><td>0.142</td><td>0.239</td></tr><tr><td>Corneal central thickness (µm)</td><td>525.8</td><td>46</td></tr><tr><td>Cup/disc ratio</td><td>0.821</td><td>0.153</td></tr><tr><td>MD 24-2 (dB)</td><td>-10.5</td><td>8.91</td></tr><tr><td>PSD 24-2 (dB)</td><td>6.81</td><td>4.7</td></tr><tr><td>VFI</td><td>70.87</td><td>27.6</td></tr><tr><td>RNFL thickness (µm)</td><td>66.48</td><td>12.6</td></tr><tr><td>GCC thickness (µm)</td><td>63.6</td><td>9.07</td></tr><tr><td>Rim area (mm²)</td><td>0.786</td><td>0.29</td></tr></tbody></table>

(r = 0.339; p = 0.017) and MD 24-2 (r = 0.301; p = 0.038). The results of these correlations are summarized in Table 4. In the multivariate analysis where total ONH vessel density was considered as a dependent variable, rim area was the dominant explanatory variable, accounting for 38.9% of the variance (R²). In this multivariate model, RNFL thickness, GCC thickness, and VF parameters (MD 24-2 and VFI)
<table><caption>TABLE 3: OCT angiography measurements. SSI: strength signal index.</caption><thead><tr><th rowspan="2">OCT angiography measurements</th><th colspan="2">Normal</th><th colspan="2">Glaucoma</th><th rowspan="2">p value</th></tr><tr><th>Mean</th><th>SD</th><th>Mean</th><th>SD</th></tr></thead><tbody><tr><td>Threshold</td><td>0.095</td><td>0.012</td><td>0.085</td><td>0.012</td><td>0.001</td></tr><tr><td>SSI</td><td>72.3</td><td>8.41</td><td>61.8</td><td>7</td><td>&lt;0.0001</td></tr><tr><td>Total ONH vessel density</td><td>0.547</td><td>0.154</td><td>0.412</td><td>0.117</td><td>&lt;0.0001</td></tr><tr><td>Temporal disc vessel density</td><td>0.472</td><td>0.105</td><td>0.364</td><td>0.15</td><td>0.0001</td></tr><tr><th rowspan="2">Glaucoma stage</th><th colspan="2">Total ONH vessel density</th><th colspan="2">Temporal disc vessel density</th><th rowspan="2"></th></tr><tr><th>Mean</th><th>SD</th><th>Mean</th><th>SD</th></tr><tr><td>Mild</td><td>0.441</td><td>0.145</td><td>0.404</td><td>0.178</td><td></td></tr><tr><td>Moderate</td><td>0.439</td><td>0.067</td><td>0.414</td><td>0.110</td><td></td></tr><tr><td>Severe</td><td>0.365</td><td>0.076</td><td>0.293</td><td>0.096</td><td></td></tr></tbody></table>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /><figcaption>FIGURE 2: Linear regression analysis for temporal disc vessel density, total ONH vessel density, and rim area.</figcaption></figure>

were not statistically associated with ONH vascular density (Table 5).

Considering temporal disc vessel density, in univariate analysis, a correlation was observed with rim area ($r = 0.624$; $p < 0.0001$) (Figure 2), RNFL thickness ($r = 0.448$; $p = 0.001$), and GCC thickness ($r = 0.395$; $p = 0.004$). Temporal disc vessel density was also associated with VFI ($r = 0.423$; $p = 0.002$), MD 24-2 ($r = 0.385$; $p = 0.007$), but not with PSD 24-2 (Table 6). In a multivariate model studying temporal disc vessel density as a dependent variable, rim area

was also the dominant explanatory variable, accounting for 39% of the variance ($R^2$). RNFL thickness, GCC thickness, and VF parameters (MD 24-2, VFI) were not associated with temporal disc vessel density (Table 7).

## 4. Discussion

In this study, we report the use of spectral-domain OCT angiography to measure ONH perfusion in glaucoma. We targeted two regions, the whole disc and a temporal disc area, to compare the two groups. We showed a significant difference in ONH vascular density between glaucoma and normal patients. The difference was similar for the total and temporal vascular area, with a decrease of ONH perfusion as compared with controls by 24.7% and 22.9%, respectively. We observed a significant correlation between ONH perfusion and structural and functional glaucoma damage. Moreover, the multivariate regression models confirmed the relation between disc rim area and ONH perfusion (total and temporal vessel density). These results emphasized the relation between ONH BF and structural loss in glaucoma.

In a previous similar study using a prototype of swept-source angio-OCT, Jia et al. studied ONH total BF with the SSADA. These authors found a correlation between the disc flow index and PSD, but they did not observe a correlation with structural variables such as the C/D ratio, rim area, or RNFL thickness [14]. Although altered microcirculation is reported to be associated with optic nerve head changes in glaucoma [24], these authors suggested that the flow index was independent information on glaucoma severity that is not explained by structural variables alone. The differences between our two studies can be explained by random variations due to the small sample size of the study reported by Jia et al. (24 normal subjects and 11 glaucoma patients). Moreover, most of their glaucoma patients had early-stage disease. In our study, 38% of the patients had severe glaucoma with advanced structural changes. The correlations we observed between ONH perfusion and structural parameters such as rim area, RNFL thickness, and GCC thickness could be explained by the more advanced glaucoma cases included in the present study. Interestingly, in these more advanced glaucoma patients, we also found a low but significant
<table><caption>TABLE 4: Pearson correlation coefficient matrix for visual field, structural variables, and ONH total vessel density in subjects with glaucoma. MD: mean deviation; PSD: pattern standard deviation; VFI: visual field index; RNFL: retinal nerve fiber layers; GCC: ganglion cell complex; total ONH density: total ONH vessel density; SSI: strength signal index. Statistically significant correlations with total ONH vessel density are in bold.</caption><thead><tr><th>Variable</th><th></th><th>Total ONH density</th><th>SSI</th><th>VFI</th><th>PSD 24-2</th><th>MD 24-2</th><th>GCC</th><th>Rim area</th><th>RNFL</th><th>Age</th></tr></thead><tbody><tr><td rowspan="2">Total ONH density</td><td>r</td><td><strong>1</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>p</td><td><strong>0</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">SSI</td><td>r</td><td><strong>0.318</strong></td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>p</td><td><strong>0.026</strong></td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">VFI</td><td>r</td><td><strong>0.339</strong></td><td>0.248</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>p</td><td><strong>0.017</strong></td><td>0.089</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">PSD 24-2</td><td>r</td><td>-0.192</td><td>-0.090</td><td>-0.595</td><td>1</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>p</td><td>0.185</td><td>0.544</td><td>&lt;0.0001</td><td>0</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">MD 24-2</td><td>r</td><td><strong>0.301</strong></td><td>0.202</td><td>0.986</td><td>-0.685</td><td>1</td><td></td><td></td><td></td><td></td></tr><tr><td>p</td><td><strong>0.038</strong></td><td>0.173</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>0</td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">GCC thickness</td><td>r</td><td><strong>0.337</strong></td><td>0.230</td><td>0.580</td><td>-0.472</td><td>0.532</td><td>1</td><td></td><td></td><td></td></tr><tr><td>p</td><td><strong>0.017</strong></td><td>0.111</td><td>&lt;0.0001</td><td>0.001</td><td>0.000</td><td>0</td><td></td><td></td><td></td></tr><tr><td rowspan="2">Rim area</td><td>r</td><td><strong>0.609</strong></td><td>0.187</td><td>0.622</td><td>-0.581</td><td>0.628</td><td>0.598</td><td>1</td><td></td><td></td></tr><tr><td>p</td><td>&lt;0.0001</td><td>0.198</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>0</td><td></td><td></td></tr><tr><td rowspan="2">RNFL thickness</td><td>r</td><td><strong>0.406</strong></td><td>0.201</td><td>0.649</td><td>-0.533</td><td>0.664</td><td>0.789</td><td>0.661</td><td>1</td><td></td></tr><tr><td>p</td><td><strong>0.003</strong></td><td>0.166</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>0</td><td></td></tr><tr><td rowspan="2">Age</td><td>r</td><td>0.166</td><td>-0.308</td><td>0.268</td><td>-0.285</td><td>0.350</td><td>0.047</td><td>0.267</td><td>0.244</td><td>1</td></tr><tr><td>p</td><td>0.248</td><td>0.031</td><td>0.063</td><td>0.047</td><td>0.015</td><td>0.746</td><td>0.061</td><td>0.088</td><td>0</td></tr></tbody></table>

<table><caption>TABLE 5: Multivariate model for total ONH vessel density R² = 0.387.</caption><thead><tr><th>Variable</th><th>Slope</th><th>p value</th></tr></thead><tbody><tr><td>MD 24-2</td><td>-0.002</td><td>0.870</td></tr><tr><td>VFI</td><td>-0.001</td><td>0.955</td></tr><tr><td>GCC thickness</td><td>-0.002</td><td>0.484</td></tr><tr><td>Rim area</td><td><strong>0.269</strong></td><td><strong>0.000</strong></td></tr><tr><td>RNFL thickness</td><td>0.002</td><td>0.445</td></tr></tbody></table>

MD: mean deviation; VFI: visual field index; RNFL: retinal nerve fiber layers; GCC: ganglion cell complex. Statistically significant correlations with total ONH vessel density are in bold.

correlation with VF variables in univariate analysis, as did Jia et al. These results are also in accordance with the study by Liu et al. using OCT-A (RTVue-XR; Optovue, Inc.) that showed reduced peripapillary vessel density in glaucoma patients correlated with VF indexes [23]. In a recent study, Wang et al. comparing optic disc perfusion with OCT-A between 62 glaucoma patients and 20 control subjects found similar results with reduced disc vessel density in the glaucoma group correlated with structural parameters such as the GCC [25].

The SSADA algorithm requires the user to define a threshold of circulating vessel detection. Jia et al. described a method to calculate the appropriate threshold using the central avascular zone as a noise region and calculating the mean decorrelation value of the whole avascular zone. With their OCT-A prototype, they set the threshold two standard deviations above the mean avascular zone decorrelation. In the present study, using the commercially available spectral-domain OCT RT XR Avanti with the AngioVue software (Optovue, Inc. Fremont, CA, USA), we could not detect

obvious circulating vessels using the same threshold as Jia et al. Therefore we set our threshold one standard deviation above the mean avascular decorrelation value. This difference might be explained by the difference in OCT-A hardware and software used in the two studies.

The results obtained with OCT-A confirmed the findings of several other techniques that reported ONH BF changes in glaucoma. Harris et al. measured ONH BF with different imaging technologies including color Doppler imaging, confocal scanning laser ophthalmoscopic angiography, laser blood flowmetry, or scanning laser Doppler flowmetry and demonstrated reduced ONH perfusion in glaucoma patients [6, 7]. Disc BF was also investigated by fluorescein angiography (FA) in several studies showing diffuse disc hypoperfusion and ONH fluorescein filling defects in glaucomatous patients [26]. Nevertheless, FA is not commonly used for monitoring glaucoma because of quantification problems and its potential side effects as an invasive technique [27]. Decreased ONH BF in glaucoma patients versus controls and increased viscosity in the papillary microvascular network were reported with the use of the laser Doppler velocimeter [28]. Using scanning laser Doppler flowmetry, Michelson et al. showed decreased peripapillary retinal BF and neuroretinal BF comparing glaucoma patients to control subjects [29]. Hafez et al. also found lower ONH BF in glaucoma patients and suggested a correlation between visual field defects and ONH perfusion [30]. Using laser speckle flowgraphy, other authors reported reduced BF in glaucoma patients and found a correlation between BF and MD and RNFL thickness [31]. Wang et al. confirmed a high correlation between the BF reduction measured by laser speckle flowgraphy and the microsphere method, providing evidence that this technique is capable of assaying BF for a critical deep ONH region [11].
<table><caption>TABLE 6: Pearson correlation coefficient matrix for visual field, structural variables, and temporal disc vessel density in subjects with glaucoma. MD: mean deviation; PSD: pattern standard deviation; VFI: visual field index; RNFL: retinal nerve fiber layers; GCC: ganglion cell complex; SSI: strength signal index. Statistically significant correlations with temporal disc vessel density are in bold.</caption><thead><tr><th>Variable</th><th>Temporal disc density</th><th>SSI</th><th>VFI</th><th>PSD 24-2</th><th>MD 24-2</th><th>GCC</th><th>Rim area</th><th>RNFL</th><th>Age</th></tr></thead><tbody><tr><td rowspan="2">Temporal disc density</td><td>r</td><td><b>1</b></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>p</td><td><b>0</b></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">SSI</td><td>r</td><td>0.235</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>p</td><td>0.105</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">VFI</td><td>r</td><td><b>0.423</b></td><td>0.248</td><td>1</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>p</td><td><b>0.002</b></td><td>0.089</td><td>0</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">PSD 24-2</td><td>r</td><td>-0.171</td><td>-0.090</td><td>-0.595</td><td>1</td><td></td><td></td><td></td><td></td></tr><tr><td>p</td><td>0.241</td><td>0.544</td><td>&lt;0.0001</td><td>0</td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">MD 24-2</td><td>r</td><td><b>0.385</b></td><td>0.202</td><td>0.986</td><td>-0.685</td><td>1</td><td></td><td></td><td></td></tr><tr><td>p</td><td><b>0.007</b></td><td>0.173</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>0</td><td></td><td></td><td></td></tr><tr><td rowspan="2">GCC thickness</td><td>r</td><td><b>0.396</b></td><td>0.230</td><td>0.580</td><td>-0.472</td><td>0.532</td><td>1</td><td></td><td></td></tr><tr><td>p</td><td><b>0.004</b></td><td>0.111</td><td>&lt;0.0001</td><td>0.001</td><td>0.000</td><td>0</td><td></td><td></td></tr><tr><td rowspan="2">Rim area</td><td>r</td><td><b>0.624</b></td><td>0.187</td><td>0.622</td><td>-0.581</td><td>0.628</td><td>0.598</td><td>1</td><td></td></tr><tr><td>p</td><td>&lt;0.0001</td><td>0.198</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>0</td><td></td></tr><tr><td rowspan="2">RNFL thickness</td><td>r</td><td><b>0.448</b></td><td>0.201</td><td>0.649</td><td>-0.533</td><td>0.664</td><td>0.789</td><td>0.661</td><td>1</td></tr><tr><td>p</td><td><b>0.001</b></td><td>0.198</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>0</td></tr><tr><td rowspan="2">Age</td><td>r</td><td>0.185</td><td>-0.308</td><td>0.268</td><td>-0.285</td><td>0.350</td><td>0.047</td><td>0.267</td><td>0.244</td><td>1</td></tr><tr><td>p</td><td>0.199</td><td>0.031</td><td>0.063</td><td>0.047</td><td>0.015</td><td>0.746</td><td>0.061</td><td>0.088</td><td>0</td></tr></tbody></table>

<table><caption>TABLE 7: Multivariate model for temporal disc vessel density $R^2$ = 0.390.</caption><thead><tr><th>Variable</th><th>Slope</th><th>p value</th></tr></thead><tbody><tr><td>MD 24-2</td><td>0.005</td><td>0.717</td></tr><tr><td>VFI</td><td>-0.002</td><td>0.680</td></tr><tr><td>GCC thickness</td><td>-0.001</td><td>0.690</td></tr><tr><td>Rim area</td><td><b>0.311</b></td><td><b>0.001</b></td></tr><tr><td>RNFL thickness</td><td>0.001</td><td>0.648</td></tr></tbody></table>

MD: mean deviation; VFI: visual field index; RNFL: retinal nerve fiber layers; GCC: ganglion cell complex. Statistically significant correlations with temporal disc vessel density are in bold.

The use of OCT-A to evaluate ONH BF may provide a new tool to investigate the role of the microvasculature in the pathophysiology of glaucoma and its changes with glaucoma progression. One of the longstanding questions is whether reduced BF is a preexisting pathologic change or a consequence of glaucomatous optic atrophy [32]. This is difficult to determine in part because a substantial portion of retinal ganglion cells/axons may have already become degenerated, even at an early stage of glaucoma [33]. Compromised blood supply and therefore reduced ocular BF due to increased IOP or vascular dysregulation has been implicated in the pathogenesis of the disease [7, 32]. In terms of laminar capillary volume flow, Burgoyne et al. hypothesized that, within the lamina, connective tissue IOP-related strain has both direct and indirect effects on axonal nutrition and that axonal ischemia can result from either IOP-induced occlusion of the laminar capillaries (direct effect) or decreased diffusion of nutrients (indirect effect), or both. In terms of only the retrolaminar effects on BF, the level of IOP-related strain within the peripapillary sclera may significantly affect volume

flow through the scleral branches of the short posterior ciliary arteries. This hypothesis suggests that vascular mechanisms damaging the ONH are not necessarily IOP-independent and proposes a logic for understanding the complicated interaction between these three important factors (IOP, volume BF, and nutrient delivery) within the tissues of the ONH [34].

Our study has limitations that need to be considered when interpreting the results. Normal and glaucomatous eyes have been classified according to the visualization the ONH and no visual field was performed to control subjects. As OCT-A was used to evaluate ONH vascularization, we may have increased the chance of finding significant differences between groups. Similarly, by selecting the worst eye of glaucomatous patients we may have influenced the results. Nevertheless, this limitation has not influenced the analysis of correlation of OCT-A results and optic nerve parameters for glaucoma patients as this group was classified according to the results of VF.

The results of the present study suggest that spectral-domain OCT angiography might detect reduced ONH blood flow in glaucoma patients. These ONH blood flow changes were associated with structural and functional glaucomatous alterations. New studies are needed to evaluate the exact potential of OCT-A in glaucoma.

## Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

## References

[1] D. C. Broadway and S. M. Drance, "Glaucoma and vasospasm," *British Journal of Ophthalmology*, vol. 82, no. 8, pp. 862-870, 1998.
[2] B. L. Petrig, C. E. Riva, and S. S. Hayreh, “Laser Doppler flowmetry and optic nerve head blood flow,” *American Journal of Ophthalmology*, vol. 127, no. 4, pp. 413–425, 1999.

[3] C. Prünte, S. Orgül, and J. Flammer, “Abnormalities of micro-circulation in glaucoma: facts and hints,” *Current Opinion in Ophthalmology*, vol. 9, no. 2, pp. 50–55, 1998.

[4] O.-Y. Tektas, E. Lütjen-Drecoll, and M. Scholz, “Qualitative and quantitative morphologic changes in the vasculature and extracellular matrix of the prelaminar optic nerve head in eyes with POAG,” *Investigative Ophthalmology and Visual Science*, vol. 51, no. 10, pp. 5083–5091, 2010.

[5] S. S. Hayreh, “Blood supply of the optic nerve head and its role in optic atrophy, glaucoma, and oedema of the optic disc,” *British Journal of Ophthalmology*, vol. 53, no. 11, pp. 721–748, 1969.

[6] E. Rechtman, A. Harris, R. Kumar et al., “An update on retinal circulation assessment technologies,” *Current Eye Research*, vol. 27, no. 6, pp. 329–343, 2003.

[7] A. Harris, L. Kagemann, R. Ehrlich, C. Rospigliosi, D. Moore, and B. Siesky, “Measuring and interpreting ocular blood flow and metabolism in glaucoma,” *Canadian Journal of Ophthalmology*, vol. 43, no. 3, pp. 328–336, 2008.

[8] C. P. Jonescu-Cuypers, H. S. Chung, L. Kagemann, Y. Ishii, D. Zarfati, and A. Harris, “New neuroretinal rim blood flow evaluation method combining Heidelberg retina flowmetry and tomography,” *British Journal of Ophthalmology*, vol. 85, no. 3, pp. 304–309, 2001.

[9] T. Sugiyama, M. Araie, C. E. Riva, L. Schmetterer, and S. Orgul, “Use of laser speckle flowgraphy in ocular blood flow research,” *Acta Ophthalmologica*, vol. 88, no. 7, pp. 723–729, 2010.

[10] Y. Wang, B. A. Bower, J. A. Izatt, O. Tan, and D. Huang, “Retinal blood flow measurement by circumpapillary Fourier domain Doppler optical coherence tomography,” *Journal of Biomedical Optics*, vol. 13, no. 6, Article ID 064003, 2008.

[11] L. Wang, G. A. Cull, C. Piper, C. F. Burgoyne, and B. Fortune, “Anterior and posterior optic nerve head blood flow in nonhuman primate experimental glaucoma model measured by laser speckle imaging technique and microsphere method,” *Investigative Ophthalmology and Visual Science*, vol. 53, no. 13, pp. 8303–8309, 2012.

[12] Y. Jia, O. Tan, J. Tokayer et al., “Split-spectrum amplitude-decorrelation angiography with optical coherence tomography,” *Optics Express*, vol. 20, no. 4, pp. 4710–4725, 2012.

[13] Y. Jia, J. C. Morrison, J. Tokayer et al., “Quantitative OCT angiography of optic nerve head blood flow,” *Biomedical Optics Express*, vol. 3, no. 12, pp. 3127–3137, 2012.

[14] Y. Jia, E. Wei, X. Wang et al., “Optical coherence tomography angiography of optic disc perfusion in glaucoma,” *Ophthalmology*, vol. 121, no. 7, pp. 1322–1332, 2014.

[15] Y. Tamaki, M. Araie, K. Tomita, M. Nagahara, A. Tomidokoro, and H. Fujii, “Real-time measurement of human optic nerve head and choroid circulation, using the laser speckle phenomenon,” *Japanese Journal of Ophthalmology*, vol. 41, no. 1, pp. 49–54, 1997.

[16] A. S. Hafez, R. L. G. Bizzarro, M. Rivard et al., “Reproducibility of retinal and optic nerve head perfusion measurements using scanning laser Doppler flowmetry,” *Ophthalmic Surgery, Lasers & Imaging*, vol. 34, no. 5, pp. 422–432, 2003.

[17] M. O. Gordon, J. A. Beiser, J. D. Brandt et al., “The Ocular Hypertension Treatment Study: baseline factors that predict the onset of primary open-angle glaucoma,” *Archives of Ophthalmology*, vol. 120, no. 6, pp. 714–730, 2002.

[18] J. L. Keltner, C. A. Johnson, K. E. Cello et al., “Classification of visual field abnormalities in the ocular hypertension treatment study,” *Archives of Ophthalmology*, vol. 121, no. 5, pp. 643–650, 2003.

[19] E. Hodapp, R. K. Parrish II, and D. R. Anderson, *Clinical Decisions in Glaucoma*, The CV Mosby, St Louis, Miss, USA, 1993.

[20] R. F. Spaide, J. M. Klancnik, and M. J. Cooney, “Retinal vascular layers imaged by fluorescein angiography and optical coherence tomography angiography,” *JAMA Ophthalmology*, vol. 133, no. 1, pp. 45–50, 2015.

[21] R. F. Spaide, J. M. Klancnik, and M. J. Cooney, “Retinal vascular layers in macular telangiectasia type 2 imaged by optical coherence tomographic angiography,” *JAMA Ophthalmology*, vol. 133, no. 1, pp. 66–73, 2015.

[22] N. G. Strouthidis, H. Yang, J. F. Reynaud et al., “Comparison of clinical and spectral domain optical coherence tomography optic disc margin anatomy,” *Investigative Ophthalmology and Visual Science*, vol. 50, no. 10, pp. 4709–4718, 2009.

[23] L. Liu, Y. Jia, H. L. Takusagawa et al., “Optical coherence tomography angiography of the peripapillary retina in glaucoma,” *JAMA Ophthalmology*, vol. 133, no. 9, pp. 1045–1052, 2015.

[24] N. Chiba, K. Omodaka, Y. Yokoyama et al., “Association between optic nerve blood flow and objective examinations in glaucoma patients with generalized enlargement disc type,” *Clinical Ophthalmology*, vol. 5, no. 1, pp. 1549–1556, 2011.

[25] X. Wang, C. Jiang, T. Ko et al., “Correlation between optic disc perfusion and glaucomatous severity in patients with open-angle glaucoma: an optical coherence tomography angiography study,” *Graefe’s Archive for Clinical and Experimental Ophthalmology*, vol. 253, no. 9, pp. 1557–1564, 2015.

[26] R. A. Hitchings and G. L. Spaeth, “Fluorescein angiography in chronic simple and low tension glaucoma,” *British Journal of Ophthalmology*, vol. 61, no. 2, pp. 126–132, 1977.

[27] M. R. Stein and C. W. Parker, “Reactions following intravenous fluorescein,” *American Journal of Ophthalmology*, vol. 72, no. 5, pp. 861–868, 1971.

[28] P. Hamard, H. Hamard, J. Dufaux, and S. Quesnot, “Optic nerve head blood flow using a laser Doppler velocimeter and haemorheology in primary open angle glaucoma and normal pressure glaucoma,” *British Journal of Ophthalmology*, vol. 78, no. 6, pp. 449–453, 1994.

[29] G. Michelson, M. J. Langhans, and M. J. Groh, “Perfusion of the juxtapapillary retina and the neuroretinal rim area in primary open angle glaucoma,” *Journal of Glaucoma*, vol. 5, no. 2, pp. 91–98, 1996.

[30] A. S. Hafez, R. L. G. Bizzarro, and M. R. Lesk, “Evaluation of optic nerve head and peripapillary retinal blood flow in glaucoma patients, ocular hypertensives, and normal subjects,” *American Journal of Ophthalmology*, vol. 136, no. 6, pp. 1022–1031, 2003.

[31] Y. Yokoyama, N. Aizawa, N. Chiba et al., “Significant correlations between optic nerve head microcirculation and visual field defects and nerve fiber layer loss in glaucoma patients with myopic glaucomatous disk,” *Clinical Ophthalmology*, vol. 5, no. 1, pp. 1721–1727, 2011.

[32] J. Flammer, S. Orgül, V. P. Costa et al., “The impact of ocular blood flow in glaucoma,” *Progress in Retinal and Eye Research*, vol. 21, no. 4, pp. 359–393, 2002.

[33] L. A. Kerrigan-Baumrind, H. A. Quigley, M. E. Pease, D. F. Kerrigan, and R. S. Mitchell, “Number of ganglion cells in
glaucoma eyes compared with threshold visual field tests in the
same persons,” *Investigative Ophthalmology and Visual Science*,
vol. 41, no. 3, pp. 741–748, 2000.

[34] C. F. Burgoyne, J. Crawford Downs, A. J. Bellezza, J.-K. Francis Suh, and R. T. Hart, “The optic nerve head as a biomechanical structure: a new paradigm for understanding the role of IOP-related stress and strain in the pathophysiology of glaucomatous optic nerve head damage,” *Progress in Retinal and Eye Research*, vol. 24, no. 1, pp. 39–73, 2005.
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

Hindawi
Submit your manuscripts at
http://www.hindawi.com

<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /></figure>

<figure><img src="image_20.png" /></figure>

<figure><img src="image_21.png" /></figure>

<figure><img src="image_22.png" /></figure>