Comparisons between retinal vessel calibers and various optic disc morphologic parameters with different optic disc appearances: The Glaucoma Stereo Analysis Study

Abstract

The Glaucoma Stereo Analysis Study (GSAS) is a multicenter collaborative study of the characteristics of glaucomatous optic disc morphology using a stereo fundus camera. This study evaluated the retinal vessel calibers and correlations using GSAS fundus photographs between retinal vessels and 38 optic nerve head (ONH) morphologic parameters comprehensively. In all 240 eyes, the mean central retinal arteriolar equivalent (CRAE) and central retinal venular equivalent (CRVE) were 138.4 and 216.5 μm, respectively; the CRAE correlated with age, visual field scores and 19 ONH parameters and CRVE correlated with age, intraocular pressure, visual field scores and 11 ONH parameters. Among the different optic disc appearances including focal ischemia (FI) (n = 53, 22%), generalized enlargement (GE) (n = 53, 22%), myopic glaucoma (MY) (n = 112, 47%), and senile sclerosis (SS) (n = 22, 9%), the CRAE did not differ significantly; CRVE was significantly narrower in SS than in FI and MY. In FI, GE, MY, and SS disc types, CRAE correlated with 3, 14, 9, and 2 ONH parameters, respectively, and CRVE corelated with 9, 0, 12, and 6 ONH parameters, respectively. We confirmed previous observations on the effect of retinal vessel narrowing on glaucomatous changes in the ONH and visual field. The associations between retinal vessel caliber and ONH morphologic parameters vary among different optic disc appearances, suggesting different effects of vascular changes in each disc type.

About the Authors:

Kazunobu Sugihara

Roles Data curation, Formal analysis, Writing – original draft

Affiliation: Department of Ophthalmology, Shimane University Faculty of Medicine, Izumo, Japan

Yasuyuki Takai

Roles Data curation, Formal analysis, Writing – original draft

Affiliation: Department of Ophthalmology, Shimane University Faculty of Medicine, Izumo, Japan

Ryo Kawasaki

Roles Data curation, Writing – review & editing

Affiliation: Department of Vision Informatics, Osaka University Graduate School of Medicine, Osaka Japan

Koji Nitta

Roles Data curation, Writing – review & editing

Affiliation: Department of Ophthalmology, Fukui-ken Saiseikai Hospital, Fukui, Japan

ORCID logo https://orcid.org/0000-0001-5139-8539

Maki Katai

Roles Data curation, Writing – review & editing

Affiliation: Department of Ophthalmology, NTT Medical Center Sapporo, Sapporo, Japan

ORCID logo https://orcid.org/0000-0002-7975-2605

Yasushi Kitaoka

Roles Data curation, Writing – review & editing

Affiliation: Department of Ophthalmology, St. Marianna University School of Medicine, Kawasaki, Japan

Yu Yokoyama

Roles Data curation, Writing – review & editing

Affiliation: Department of Ophthalmology, Tohoku University Graduate School of Medicine, Sendai, Japan

Kazuko Omodaka

Roles Data curation, Writing – review & editing

Affiliation: Department of Ophthalmology, Tohoku University Graduate School of Medicine, Sendai, Japan

Tomoko Naito

Roles Data curation, Writing – review & editing

Affiliation: Grace Eye Clinic, Okayama, Japan

Takehiro Yamashita

Roles Data curation, Writing – review & editing

Affiliation: Department of Ophthalmology, Kagoshima University Graduate School of Medical and Dental Sciences, Kagoshima, Japan

Shiro Mizoue

Roles Data curation, Writing – review & editing

Affiliation: Department of Ophthalmology, Ehime University Graduate School of Medicine, Toon, Japan

Aiko Iwase

Roles Data curation, Writing – review & editing

Affiliation: Tajimi Iwase Eye Clinic, Gifu, Japan

Toru Nakazawa

Roles Data curation, Writing – review & editing

Affiliation: Department of Ophthalmology, Tohoku University Graduate School of Medicine, Sendai, Japan

Masaki Tanito

Roles Conceptualization, Data curation, Formal analysis, Writing – original draft

* E-mail: [email protected]

Affiliation: Department of Ophthalmology, Shimane University Faculty of Medicine, Izumo, Japan

ORCID logo https://orcid.org/0000-0001-6512-7203

Introduction

Glaucoma is a leading cause of irreversible blindness in Japan and worldwide [1, 2]. Optic disc damage due to retinal ganglion cell (RGC) and RGC axon loss is characteristic of glaucoma [3], making it essential to detect glaucomatous optic neuropathy (GON) using different modalities for disease diagnosis. GON accompanies alterations in the optic disc morphologic parameters that measure depth, i.e., excavation of optic disc cupping; however, performing ophthalmoscopy or subjectively examining one fundus photograph monoscopically may underestimate optic disc cupping and glaucoma severity [4, 5]. Using a noninvasive noncontact imaging technique to topographically analyze optic discs that does not require dilation of the pupils, such as simultaneous stereo fundus camera, provides excellent reproducibility and interexaminer consistency [6], making it promising to objectively assess the morphologic changes in GON. This technique was used in a multicenter collaborative study, the Glaucoma Stereo Analysis Study (GSAS), to assess various optic disc morphologic parameters in Japanese patients with primary open-angle glaucoma (POAG) [7–11].

There are substantial differences among patients with glaucoma in the patterns and progression of visual field defects and prevalence of risk factors, such as the intraocular pressure (IOP) and refractive error [12, 13]. A previous study suggested classifying various GONs as follows into subgroups based on the appearance of the optic disc: focal ischemic (FI), generalized enlargement (GE), myopic glaucoma (MY), and senile sclerosis (SS) [14]. Eyes of patients with differences in the appearance of the optic disc, selected based on only the assessment of fundus photographs, exhibited different demographic characteristics, prevalence rates of some systemic and ocular risk factors, IOP levels, and patterns of visual field damage [14–21]. Therefore, classifying GONs based on the appearance of the optic disc may facilitate more accurate diagnoses and better disease management because it may reflect different pathologic mechanisms of glaucoma.

It has been suggested that vascular pathologies or blood flow changes may play roles in glaucoma development and progression. In normal-tension glaucoma (NTG), optic disc hemorrhage [22], migraine [23], high frequency of nighttime hypotension [24], lower recovery rate from a cold recovery test [25], lower blood flow velocities, and higher resistive indices in patients with progressive visual field defects [26] may be associated with the pathology. Reduced retrobulbar and retinal blood flow over time can be associated with structural glaucomatous progression [27]; oral calcium channel blockers were reported to maintain the optic disc rim, visual fields, and posterior choroidal circulation in OAG [28]. Jonas et al. reported peripapillary retinal vessel narrowing in glaucomatous eyes [29, 30]. Analyses of optic disc photographs showed that arteriolar narrowing was more common in patients with glaucoma than ocular hypertension [31]. Large-scale population-based studies have reported associations between retinal arteriolar and/or venular narrowing and elevated glaucoma prevalence [32–35] and incidence [36].

Vessel caliber measurements can be obtained using fundus photography [29, 30], retinal vessel analyzer [37], retinal oximetry [38, 39], scanning laser ophthalmoscopy [40], spectral-domain optical coherence tomography (OCT) [41, 42], OCT angiography [43], Doppler OCT [43], and adaptive optics imaging [44]. Analysis of digital color fundus photographs using semi-automated methods such as interactive vessel analysis reproducibly measured arteriolar and venular calibers [45]. The current study measured retinal vessel calibers in fundus photographs obtained during the GSAS and assessed correlations between the retinal vessels and optic nerve head (ONH) morphologic parameters; the assessment was conducted in the entire dataset and for each optic disc appearance.

Methods

Subjects

This study adhered to the tenets of the Declaration of Helsinki. The institutional review boards of Shimane University Hospital, Fukui-ken Saiseikai Hospital, Sapporo Teishin Hospital, St. Marianna University School of Medicine, Tohoku University Graduate School of Medicine, Okayama University Graduate School of Medicine, Kagoshima University Graduate School of Medical and Dental Sciences, Ehime University Graduate School of Medicine, and Tajimi Iwase Eye Clinic reviewed and approved the research. The subjects provided written informed consent; otherwise, based on the regulations of the Japanese Guidelines for Epidemiologic Study issued by the Japanese Government, the study protocols did not require each patient provide written informed consent. The protocol was posted at the outpatient clinic to inform patients about the study. All subjects were recruited from among nine institutions (Shimane University Hospital, Fukui-ken Saiseikai Hospital, Sapporo Teishin Hospital, Hospital of St. Marianna University School of Medicine, Tohoku University Hospital, Okayama University Hospital, Kagoshima University Hospital, Ehime University Hospital, and Tajimi Iwase Eye Clinic) as reported previously [10, 11]. In this study, of the 280 GSAS enrolled subjects, 240 eyes of 240 patients with POAG were included for whom we could estimate the retinal vessel calibers. The inclusion and exclusion criteria, methods of ophthalmic examinations, diagnosis of POAG, and rules for data collection have been reported previously [7, 9]. Dataset underlying the findings described in this manuscript is shown in S1 File.

Demographic and clinical parameters

In the current study, the following demographic and clinical parameters were extracted from the GSAS database: age, sex, spherical equivalent refractive error (SERE), treated IOP at the time of the fundus camera examination, and visual field mean deviation (MD), pattern standard deviation (PSD), and MD slope (dB/year) calculated from at least six visual fields tested during longer than 3 years. The IOP was measured using Goldmann applanation tonometry. The MD and PSD were measured using the Humphrey Visual Field Analyzer Swedish Interactive Thresholding Algorithm Standard central 30–2 or 24–2 program (Carl Zeiss Meditec, Dublin, CA, USA).

Optic disc topography

Stereo images of the ONH were obtained using a stereo fundus camera (nonmyd WX, Kowa Company, Ltd., Aichi, Japan) that produces nonmydriatic fundus stereographs and simultaneous right and left parallactic images using one optical system to handle light paths in two directions [6]. The built-in software (VK-2 WX, prototype version, Kowa Company, Ltd.) automatically calculates the ONH morphologic parameters based on manually set contour lines for the ONH disc and cup, which in this study were determined by one of the authors (M.T.) who viewed the images stereoscopically. According to the recommendations of the Japan Glaucoma Society Guidelines for Glaucoma [46], the disc contour was delineated by the inner margin of Elschnig’s scleral ring, and the cup contour was delineated by the outer cup margin, which was indicated by the bending of the ONH vessels at the rim. The observer determined several points on the contour (typically 8–14), and the contour line then was generated automatically by software spline interpolation. Excellent intra- and inter-observer agreements between the contour delineation and the software-assisted ONH analysis were reported previously [5]. In the GSAS dataset, 38 optic disc parameters were obtained using the VK-2 WX software. The depth and volume values were calculated based on the disparity between the right and left images of the stereo image pair with a stereo matching technique, with correction for magnification by a modification of Littmann’s method using the refractive error and corneal curvature of each eye [7]. Rim decentration can range from -1, defined as rim thinning in the superotemporal area without rim thinning in the inferotemporal area, and 1, defined as rim thinning at the inferotemporal area without rim thinning at the superotemporal area; 0 indicates equally thinned rims at both the supero- and inferotemporal rims. The disc tilt angle was defined as the degree of the angle between the horizontal plane and the line drawn from the temporal to the nasal disc edge, passing through the center of the ONH; plus and minus values indicate disc tilt to the temporal and nasal sides, respectively. The details of the parameters were described previously [7]. The summary of the 38 optic disc parameters is shown in S2 File.

Grader classification of optic disc appearance

As reported previously [8], three independent graders (T.N., K.O., and Y.Y.) classified each optic disc appearance into four different types according to the proposal of Nicolela and Drance [14], i.e., a FI disc with localized tissue loss at the superior or inferior poles and a relatively intact neuroretinal rim elsewhere; a GE disc characterized by a diffusely enlarged round cup and lack of localized defects of the neuroretinal rim; a MY disc that was tilted with temporal crescent peripapillary atrophy (PPA), excluding discs with degenerative myopia; and a SS disc with a saucerized shallow cup and diffuse neuroretinal rim tissue loss accompanied by surrounding PPA and choroidal sclerosis. Discs with features of multiple (mixed) disc types were assigned to the most prominent type. Any disagreements among the three graders were resolved by discussion. Using Cohen’s kappa statistics, the agreement among the graders was considered substantial to excellent [10, 11]. Examples of each optic disc appearance are shown in S3 File.

Retinal vessel caliber measurement

The retinal vessel diameters were estimated on stereo fundus photographs collected during the GSAS using IVAN software developed by the Department of Ophthalmology and Visual Science, University of Wisconsin, Madison, WI, USA. Either one of the photographs from a stereo pair of 20-degree ONH photographs were used for the measurement. The detailed measurement method has been described previously [47, 48]. Briefly, all vessels 25 microns or larger were measured and passed completely through a circumferential zone 0.5 to 1 disc diameter from the optic disc margin [32]. Trained graders (K.S., Y.T.) identified each vessel as an arteriole or a venule. The six widest arteriolar and venular diameters then were defined as the central retinal arteriolar equivalent (CRAE) and the central retinal venular equivalent (CRVE), respectively, using the revised Parr-Hubbard formulas of Knudtson et al. [49]. The magnification of the optic media in an eye was corrected according to the method of Littmann, with consideration of the SERE and corneal curvature [50]. The examples of original stereophotograph (S4 File), identified vessels on IVAN software (S5 File), defined values of CRAE (S6 File) and CRVE (S7 File) by IVAN software, magnification factor calculated by Kowa nonmyd WX Tool (Kowa Company) (S8 File; in this example, the magnification factor is 1.8360 mm/300 pixel), and accompanying optic disc topography results by VK-2 WX software are found in the Supporting information files.

Statistical analysis

All statistical analyses were performed using JMP Pro statistical software version 12.2 (SAS Institute, Inc., Cary, NC, USA). The continuous data are expressed as the means ± standard error. The continuous data were compared among the four optic disc appearances by one-way analysis of variance, which was followed by a comparison between each pair of optic disc appearances using the post-hoc Student t-test. The categorical data were compared among the four optic disc appearances using the G test. Based on Bonferroni’s method for correction of multiple comparisons, p<0.0083 and p<0.0017 were considered to be significance levels of 5% and 1%, respectively, in the post-hoc test. Correlations between the retinal vessel calibers (i.e., CRAE or CRVE) and demographic parameters (i.e., age, IOP, MD, MD slope, or PSD) or 38 optic disc parameters were assessed using the Spearman’s rank correlation coefficient; the analyses were performed in all subjects and each of the optic disc appearances.

Results

The graders classified the 240 study eyes as follows: 53 eyes (22%) with FI, 53 eyes with GE (22%), 112 eyes with MY (47%), and 22 eyes with SS (9%). The demographic data of the subjects are summarized in Table 1. The comparisons among the different disc appearances indicated that the patient age in the MY group was significantly younger and SERE was significantly more myopic compared with the other three groups (Table 1). Other demographic parameters, including gender, visual field MD, PSD, MD slope, and IOP did not differ among the various disc appearances. The comparisons of 38 optic disc morphologic parameters among the disc types are shown in S2 File; as expected, significant differences were detected in multiple optic disc parameters.

[Figure omitted. See PDF.]

Table 1. Comparison of demographic subject data among different optic disc appearances.

https://doi.org/10.1371/journal.pone.0250245.t001

The CRAE and the CRVE measured in the fundus photographs are summarized in Table 2. In all eyes, the mean CRAE and CRVE were 138.4 μm and 216.5 μm, respectively. Among the different disc appearances, the CRAE did not differ significantly; the CRVE was significantly narrower in the SS group than those in the FI and MY groups.

[Figure omitted. See PDF.]

Table 2. Comparison of retinal vessel calibers among different optic disc appearances.

https://doi.org/10.1371/journal.pone.0250245.t002

Possible correlations between retinal vessel parameters and various demographic parameters are shown in Tables 3 (for CRAE) and 4 (for CRVE). In all eyes, the CRAE and CRVE correlated positively with each other (ρ = 0.43, p<0.0001); the correlation was significant for each disc appearance (ρ = 0.41–0.61) except for FI (ρ = 0.24). In all eyes, the CRAE correlated positively with the MD and MD slope and negatively with age and PSD; the CRVE correlated positively with the IOP, MD, and MD slope and negatively with age and SERE. In men and women, the CRAE values were 138.0±15.6 μm and 138.8±14.5 μm, respectively, and the CRVE values were 217.3±22.8 μm and 215.7±20.4 μm, respectively; no gender difference was seen for the CRAE (p = 0.6627) and CRVE (p = 0.5520).

[Figure omitted. See PDF.]

Table 3. Correlation between CRAE and demographic parameters in each optic disc appearance group.

https://doi.org/10.1371/journal.pone.0250245.t003

[Figure omitted. See PDF.]

Table 4. Correlation between CRVE and demographic parameters in each optic disc appearance group.

https://doi.org/10.1371/journal.pone.0250245.t004

Possible correlations between retinal vessel parameters and various optic disc morphologic parameters are shown in S10 (for CRAE) and S11 (for CRVE) Files. The parameters that correlated significantly with CRAE and CRVE were summarized in Tables 5 and 6, respectively. In all eyes, the CRAE correlated significantly with 19 ONH parameters (Table 5), and CRVE correlated significantly with 11 ONH parameters (Table 6). In each optic disc type, multiple parameters correlated with the CRAE and/or CRVE with various patterns. The numbers of disc morphologic parameters that correlated significantly with CRAE or CRVE in each disc appearance group are shown in Table 7. In the GE group, the CRAE correlated with 14 optic disc parameters while CRVE didn’t correlate with any parameters. In contrast, more optic disc parameters correlated with the CRVE than the CRAE in the three other disc types.

[Figure omitted. See PDF.]

Table 5. Parameters significantly correlated with CRAE (sorted by effect size).

https://doi.org/10.1371/journal.pone.0250245.t005

[Figure omitted. See PDF.]

Table 6. Parameters significantly correlated with CRVE (sorted by effect size).

https://doi.org/10.1371/journal.pone.0250245.t006

[Figure omitted. See PDF.]

Table 7. Number of optic disc parameters significantly correlated with the CRAE or CRVE.

https://doi.org/10.1371/journal.pone.0250245.t007

Discussion

In all eyes, a moderate (ρ = 0.43) but significant correlation was seen between the CRAE and CRVE; surprisingly few studies have reported this correlation previously. Narrowing of the CRAE and CRVE was associated with a larger cup-to-disc area ratio and smaller rim-to-disc area ratio or smaller rim volume. An association between narrower retinal vessels and glaucoma was seen in the arterioles and venules of subjects with glaucoma [34], OAG [32], POAG [30], or NTG [51]; and in the arterioles of subjects with glaucoma [33], POAG [38], NTG [38, 42, 52], or primary angle-closure glaucoma [38]. Regarding the morphologic disc parameters, an association between narrower retinal vessels and larger cup-to-disc ratio was found in the arterioles and venules of subjects with glaucoma [34], and in narrower arterioles and larger cup/smaller rim areas in POAG [29]. Other studies have reported an association between narrower retinal vessels and thinner retinal nerve fiber layer (RNFL) thickness in both the arterioles and venules of general populations [53], non-glaucomatous eyes [54], NTG [51], and in the arterioles of subjects with NTG [52]. Collectively, the current study confirmed the previous findings regarding involvement of vascular changes in glaucoma.

The current study showed that narrower CRAE and CRVE were associated with steeper MD slopes, and narrower CRVE were associated with smaller disc areas. Previously, an association between narrower retinal arterioles and worse visual fields was found in subjects with NTG [39, 52]. In subjects with asymmetric glaucoma, the retinal arteriolar diameter was narrower in eyes with advanced disease [55]. We previously reported that smaller disc areas (odd ratio, 0.49/mm2) were associated with faster glaucoma progression in the GSAS dataset [9]. Another group reported a significant association between smaller disc diameters and narrower retinal arteriolar and venular diameters [56]. This association also was confirmed in the current study in that both the CRAE and CRVE correlated with the disc area. Although the exact underlining mechanism is uncertain, our and previous observations have suggested a possible link among ONH size, vessel vascular morphology, and visual field progression in glaucoma.

We have found significant associations between aging and narrowing of the retinal vessel diameters in both CRAE and CRVE. Age-dependent narrowing of the retinal vessels was reported previously in both the arterioles and venules in subjects with OAG [57] and non-glaucomatous subjects older than 81 years [33], and in the arterioles in NTG [52]. Using adaptive optics imaging in POAG, narrowing of the arteriolar caliber was associated with narrowing of the vascular lumen rather than thickening of the vascular wall, suggesting dysregulation of the vascular tone other than atherosclerosis was involved in the mechanism of arteriolar narrowing [44]. Decreases in the nitric oxide (NO) marker and increases in the endothelin level were found in glaucomatous eyes [58]. Disruption of the ONH microcirculation can result in ONH ischemia and ultimately development of GON [59]. Narrowing of the retinal vessels may result from dysregulation of the vascular endothelium by a disproportionate release of endogenous vasodilators such as NO and prostacyclin and/or vasoconstrictors such as endothelin [53]. A 10-year cohort study found that patients in the lowest quartile of retinal arteriolar calibers at baseline had an approximate 4-fold higher risk of developing OAG compared with those in the highest quartile of arteriolar calibers, independent of other systemic and ocular OAG risk factors [36]. This finding was explained by whether retinal vessel caliber changes preceded subtle nerve fiber layer defects (hypothesized in the vascular theory of glaucoma) or were a subsequent change secondary to the RNFL thinning [36]. Although not conclusive, we speculated that aging can associate with vessel narrowing by its effects on tone dysregulation and atherosclerosis mechanisms.

In the GE group, both CRAE and CRVE associated with visual field parameters and CRVE associated with IOP, while only CRAE associated with ONH morphologic parameters. GE tended to be associated with younger age [14] and higher baseline IOP [10]. Previous studies have found an association between narrower retinal vessels and higher IOP in the arterioles and venules of subjects with OAG [57] and NTG [51], whereas no such association was observed in glaucoma in another study [34]. In a population-based study, no significant difference in the arteriolar caliber was found between high- and low-tension glaucoma [32]. IOP-dependent changes in retinal vessels, i.e., dilation of retinal arterioles and venules with reduction of IOP, diminished in subjects 55 years and older [57]. Among the different disc appearances, we found a significant association between vessel parameter and IOP only in the GE and MY groups; thus, a heterogenous effect of IOP on each glaucoma subject and different age distribution may explain the inconsistent results in previous and current studies. A large ONH cup is the hallmark of GE; previously, normal-caliber retinal vessels were seen in eyes with non-glaucomatous large cups, while narrowing of the retinal vessels was found in glaucomatous eyes [60]. In the GE group, the ONH microcirculation was correlated with the RNFL thickness, vertical cup-to-disc ratio, and visual field MD [61]. GE is the predominant disc type among patients with severe glaucoma [62]. Thus, GE can be a good candidate to assess the IOP-dependent glaucomatous changes in the vascular blood flow, ONH/RNFL structures, and visual functions.

Rim decentration, a marker of asymmetrical rim thinning obtained by comparing the superior and inferior rim areas, associated with CRAE in FI, while no association was detected in other disc appearances. Thus, the results might detect the focal changes in the retinal vessels in this optic disc type. The correlations between the circumpapillary RNFL thickness and visual field MD were significant in the MY, GE, and SS groups, while no association was found with FI; thus, the magnitude of the structure-function relationships may differ among different disc appearances [19]. A geographic association between the location of the retinal vessel narrowing and visual field damage was found in arterioles in subjects with POAG (no association was found in the venules) [63] or NTG [52]. Focal narrowing of the retinal arterioles accompanied thinning of the ONH rim width/area in cross-sectional [29] and longitudinal studies [64]. Thus, analyses in each segmented area of the ONH, fundus photographs, and visual field may increase the chance of detecting an association between vessel caliber and other parameters in this type of optic disc.

In the MY group, the CRAE and CRVE associated with decreased rim area and/or volume. The disc damage likelihood scale in the MY group tended to be higher than the scores in the other disc types even though the visual field MD was equivalent [8]. In myopic discs, reduced ONH microcirculation measured by laser speckle flowgraphy associated significantly with losses of RNFL and visual fields in glaucoma [65, 66]. Reduced temporal rim width and RNFL thickness, the hallmarks of MY [10], were the direct risk factors of visual acuity loss in glaucoma [67]. Accordingly, this type of ONH can be a good candidate for assessing the visual acuity and retinal vessel changes in a future study.

Among the different disc appearances, the fewest optic disc morphologic parameters associated with vessel caliber parameters in the SS group, which might associate with older age of this disc group [10, 14] and age-related loss of the vascular response with IOP change [57]. Patients with SS had a higher prevalence of ischemic heart disease and systemic hypertension [14]. A population-based study reported that the major systemic determinant of narrower CRAE was higher blood pressure, while those of wider CRVE were smokers, higher blood pressure, systemic inflammation and obesity [68], suggesting the complex mechanisms of vascular diameter regulation. Using color Doppler imaging, subjects with SS had a lower diastolic velocity and a higher resistive index in their ophthalmic, central retinal, and short posterior ciliary arteries [69]; SS associated with greater circulatory abnormalities in the retrobulbar vessels with higher downstream resistance in these vessels. These circulatory abnormalities may relate to the pathogenesis of glaucoma in SS; thus, SS is potentially a good candidate for assessing a possible association between glaucoma and age-related changes in retinal vessels and systemic diseases. The fewest subjects in this disc group may have resulted in a weaker statistical power than in the other groups; therefore, the results in this study may change if more patients are included in a future study.

Ocular hypotensive medications may have vasodilatory effects [70]; most subjects in the GSAS used medication, so that the association between glaucoma and retinal arteriolar narrowing may have been underestimated in this study. Oral antihypertensive drug may increase the retinal vessel caliber [71]; some subjects in the GSAS presumably were using antihypertensive medication, although we have no data regarding systemic drug use; thus, uncontrol for systemic drug use was a study limitation. Difference in race, environmental factors, and lifestyle habits may contribute to the results [72]; all GSAS participants were Japanese; thus, the generalizability of our results must be confirmed in different datasets.

Using the GSAS dataset, we assessed the correlation between retinal vessel narrowing and topographic changes in ONH. The associations between retinal vessel calibers and optic disc morphologic parameters vary among different optic disc types, suggesting that different etiologies may underlie glaucoma in each subject.

Supporting information

S1 File. Dataset underlying the findings described in this manuscript.

https://doi.org/10.1371/journal.pone.0250245.s001

(PDF)

S2 File. Comparison of optic nerve disc parameters among different optic disc appearances.

https://doi.org/10.1371/journal.pone.0250245.s002

(PDF)

S3 File. Examples of each optic disc appearance.

https://doi.org/10.1371/journal.pone.0250245.s003

(PDF)

S4 File. Examples of original stereophotograph.

https://doi.org/10.1371/journal.pone.0250245.s004

(JPG)

S5 File. Example of identified vessels on IVAN software.

https://doi.org/10.1371/journal.pone.0250245.s005

(PNG)

S6 File. Example of defined value of CRAE by IVAN software.

https://doi.org/10.1371/journal.pone.0250245.s006

(PNG)

S7 File. Example of defined value of CRVE by IVAN software.

https://doi.org/10.1371/journal.pone.0250245.s007

(PNG)

S8 File. Example of magnification factor calculation by Kowa nonmyd WX tool.

https://doi.org/10.1371/journal.pone.0250245.s008

(PNG)

S9 File. Example of optic disc topography results by VK-2 WX software.

https://doi.org/10.1371/journal.pone.0250245.s009

(PNG)

S10 File. Correlation between CRAE and other parameters in each optic disc appearance group.

https://doi.org/10.1371/journal.pone.0250245.s010

(PDF)

S11 File. Correlation between CRVE and other parameters in each optic disc appearance group.

https://doi.org/10.1371/journal.pone.0250245.s011

(PDF)

Acknowledgments

The authors are grateful to the GSAS group members for establishing the dataset and valuable discussion and to Dr. Toshiaki Nakagawa for his technical help with the development of the prototype optic disc analyses software.

Citation: Sugihara K, Takai Y, Kawasaki R, Nitta K, Katai M, Kitaoka Y, et al. (2021) Comparisons between retinal vessel calibers and various optic disc morphologic parameters with different optic disc appearances: The Glaucoma Stereo Analysis Study. PLoS ONE 16(7): e0250245. https://doi.org/10.1371/journal.pone.0250245

References

1. Foster A, Resnikoff S. The impact of Vision 2020 on global blindness. Eye (Lond). 2005;19(10):1133–5. pmid:16304595

2. Iwase A, Araie M, Tomidokoro A, Yamamoto T, Shimizu H, Kitazawa Y. Prevalence and causes of low vision and blindness in a Japanese adult population: the Tajimi Study. Ophthalmology. 2006;113(8):1354–62. pmid:16877074

3. Weinreb RN, Khaw PT. Primary open-angle glaucoma. Lancet. 2004;363(9422):1711–20. pmid:15158634

4. Wolfs RC, Ramrattan RS, Hofman A, de Jong PT. Cup-to-disc ratio: ophthalmoscopy versus automated measurement in a general population: The Rotterdam Study. Ophthalmology. 1999;106(8):1597–601. pmid:10442909

5. Tanito M, Sagara T, Takamatsu M, Kiuchi Y, Nakagawa T, Fujita Y, et al. Intraobserver and interobserver agreement of computer software-assisted optic nerve head photoplanimetry. Jpn J Ophthalmol. 2014;58(1):56–61. pmid:24100930

6. Asakawa K, Kato S, Shoji N, Morita T, Shimizu K. Evaluation of optic nerve head using a newly developed stereo retinal imaging technique by glaucoma specialist and non-expert-certified orthoptist. J Glaucoma. 2013;22(9):698–706. pmid:22936277

7. Yokoyama Y, Tanito M, Nitta K, Katai M, Kitaoka Y, Omodaka K, et al. Stereoscopic analysis of optic nerve head parameters in primary open angle glaucoma: the glaucoma stereo analysis study. PLoS One. 2014;9(6):e99138. pmid:24922327

8. Kitaoka Y, Tanito M, Yokoyama Y, Nitta K, Katai M, Omodaka K, et al. Estimation of the Disc Damage Likelihood Scale in primary open-angle glaucoma: the Glaucoma Stereo Analysis Study. Graefes Arch Clin Exp Ophthalmol. 2016;254(3):523–8. pmid:26666233

9. Kitaoka Y, Tanito M, Yokoyama Y, Nitta K, Katai M, Omodaka K, et al. A Small Disc Area Is a Risk Factor for Visual Field Loss Progression in Primary Open-Angle Glaucoma: The Glaucoma Stereo Analysis Study. J Ophthalmol. 2018;2018:8941489. pmid:29750124

10. Tanito M, Nitta K, Katai M, Kitaoka Y, Yokoyama Y, Omodaka K, et al. Differentiation of glaucomatous optic discs with different appearances using optic disc topography parameters: The Glaucoma Stereo Analysis Study. PLoS One. 2017;12(2):e0169858. pmid:28178303

11. Tanito M, Nitta K, Katai M, Kitaoka Y, Yokoyama Y, Omodaka K, et al. Validation of formula-predicted glaucomatous optic disc appearances: the Glaucoma Stereo Analysis Study. Acta Ophthalmol. 2019;97(1):e42–e9. pmid:30022606

12. Suzuki Y, Iwase A, Araie M, Yamamoto T, Abe H, Shirato S, et al. Risk factors for open-angle glaucoma in a Japanese population: the Tajimi Study. Ophthalmology. 2006;113(9):1613–7. pmid:16828504

13. Yamamoto S, Sawaguchi S, Iwase A, Yamamoto T, Abe H, Tomita G, et al. Primary open-angle glaucoma in a population associated with high prevalence of primary angle-closure glaucoma: the Kumejima Study. Ophthalmology. 2014;121(8):1558–65. pmid:24746386

14. Nicolela MT, Drance SM. Various glaucomatous optic nerve appearances: clinical correlations. Ophthalmology. 1996;103(4):640–9. pmid:8618765

15. Nicolela MT, McCormick TA, Drance SM, Ferrier SN, LeBlanc RP, Chauhan BC. Visual field and optic disc progression in patients with different types of optic disc damage: a longitudinal prospective study. Ophthalmology. 2003;110(11):2178–84. pmid:14597527

16. Broadway DC, Drance SM. Glaucoma and vasospasm. Br J Ophthalmol. 1998;82(8):862–70. pmid:9828767

17. Broadway DC, Drance SM, Parfitt CM, Mikelberg FS. The ability of scanning laser ophthalmoscopy to identify various glaucomatous optic disk appearances. Am J Ophthalmol. 1998;125(5):593–604. pmid:9625542

18. Broadway DC, Nicolela MT, Drance SM. Optic disk appearances in primary open-angle glaucoma. Surv Ophthalmol. 1999;43 Suppl 1:S223–43. pmid:10416767

19. Omodaka K, Takada N, Yamaguchi T, Takahashi H, Araie M, Nakazawa T. Characteristic correlations of the structure-function relationship in different glaucomatous disc types. Jpn J Ophthalmol. 2015;59(4):223–9. pmid:25860862

20. Takada N, Omodaka K, Kikawa T, Takagi A, Matsumoto A, Yokoyama Y, et al. OCT-Based Quantification and Classification of Optic Disc Structure in Glaucoma Patients. PLoS One. 2016;11(8):e0160226. pmid:27557112

21. Omodaka K, An G, Tsuda S, Shiga Y, Takada N, Kikawa T, et al. Classification of optic disc shape in glaucoma using machine learning based on quantified ocular parameters. PLoS One. 2017;12(12):e0190012. pmid:29261773

22. Kitazawa Y, Shirato S, Yamamoto T. Optic disc hemorrhage in low-tension glaucoma. Ophthalmology. 1986;93(6):853–7. pmid:3737131

23. Phelps CD, Corbett JJ. Migraine and low-tension glaucoma. A case-control study. Invest Ophthalmol Vis Sci. 1985;26(8):1105–8. pmid:4019101

24. Graham SL, Drance SM, Wijsman K, Douglas GR, Mikelberg FS. Ambulatory blood pressure monitoring in glaucoma. The nocturnal dip. Ophthalmology. 1995;102(1):61–9. pmid:7831043

25. Ishida K, Yamamoto T, Kitazawa Y. Clinical factors associated with progression of normal-tension glaucoma. J Glaucoma. 1998;7(6):372–7. pmid:9871858

26. Yamazaki Y, Drance SM. The relationship between progression of visual field defects and retrobulbar circulation in patients with glaucoma. Am J Ophthalmol. 1997;124(3):287–95. pmid:9439354

27. Tobe LA, Harris A, Hussain RM, Eckert G, Huck A, Park J, et al. The role of retrobulbar and retinal circulation on optic nerve head and retinal nerve fibre layer structure in patients with open-angle glaucoma over an 18-month period. Br J Ophthalmol. 2015;99(5):609–12. pmid:25467967

28. Koseki N, Araie M, Tomidokoro A, Nagahara M, Hasegawa T, Tamaki Y, et al. A placebo-controlled 3-year study of a calcium blocker on visual field and ocular circulation in glaucoma with low-normal pressure. Ophthalmology. 2008;115(11):2049–57. pmid:18672290

29. Jonas JB, Naumann GO. Parapapillary retinal vessel diameter in normal and glaucoma eyes. II. Correlations. Invest Ophthalmol Vis Sci. 1989;30(7):1604–11. pmid:2745001

30. Jonas JB, Nguyen XN, Naumann GO. Parapapillary retinal vessel diameter in normal and glaucoma eyes. I. Morphometric data. Invest Ophthalmol Vis Sci. 1989;30(7):1599–603. pmid:2745000

31. Rankin SJ, Drance SM. Peripapillary focal retinal arteriolar narrowing in open angle glaucoma. J Glaucoma. 1996;5(1):22–8. pmid:8795730

32. Mitchell P, Leung H, Wang JJ, Rochtchina E, Lee AJ, Wong TY, et al. Retinal vessel diameter and open-angle glaucoma: the Blue Mountains Eye Study. Ophthalmology. 2005;112(2):245–50. pmid:15691558

33. Wang S, Xu L, Wang Y, Wang Y, Jonas JB. Retinal vessel diameter in normal and glaucomatous eyes: the Beijing eye study. Clin Exp Ophthalmol. 2007;35(9):800–7. pmid:18173406

34. Amerasinghe N, Aung T, Cheung N, Fong CW, Wang JJ, Mitchell P, et al. Evidence of retinal vascular narrowing in glaucomatous eyes in an Asian population. Invest Ophthalmol Vis Sci. 2008;49(12):5397–402. pmid:18719076

35. Kiuchi Y, Yanagi M, Itakura K, Takahashi I, Hida A, Ohishi W, et al. Association between radiation, glaucoma subtype, and retinal vessel diameter in atomic bomb survivors. Sci Rep. 2019;9(1):8642. pmid:31201344

36. 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(1):84–90. pmid:23062656

37. Garhofer G, Bek T, Boehm AG, Gherghel D, Grunwald J, Jeppesen P, et al. Use of the retinal vessel analyzer in ocular blood flow research. Acta Ophthalmol. 2010;88(7):717–22. pmid:19681764

38. Cheng CS, Lee YF, Ong C, Yap ZL, Tsai A, Mohla A, et al. Inter-eye comparison of retinal oximetry and vessel caliber between eyes with asymmetrical glaucoma severity in different glaucoma subtypes. Clin Ophthalmol. 2016;10:1315–21. pmid:27486303

39. Yap ZL, Ong C, Lee YF, Tsai A, Cheng C, Nongpiur ME, et al. Retinal Oximetry in Subjects With Glaucomatous Hemifield Asymmetry. J Glaucoma. 2017;26(4):367–72. pmid:28079654

40. Pellegrini E, Robertson G, Trucco E, MacGillivray TJ, Lupascu C, van Hemert J, et al. Blood vessel segmentation and width estimation in ultra-wide field scanning laser ophthalmoscopy. Biomed Opt Express. 2014;5(12):4329–37. pmid:25574441

41. Park SC, De Moraes CG, Teng CC, Tello C, Liebmann JM, Ritch R. Enhanced depth imaging optical coherence tomography of deep optic nerve complex structures in glaucoma. Ophthalmology. 2012;119(1):3–9. pmid:21978593

42. Yabana T, Shiga Y, Kawasaki R, Omodaka K, Takahashi H, Kimura K, et al. Evaluating retinal vessel diameter with optical coherence tomography in normal-tension glaucoma patients. Jpn J Ophthalmol. 2017;61(5):378–87. pmid:28667424

43. Hosseinaee Z, Tan B, Martinez A, Bizheva KK. Comparative Study of Optical Coherence Tomography Angiography and Phase-Resolved Doppler Optical Coherence Tomography for Measurement of Retinal Blood Vessels Caliber. Transl Vis Sci Technol. 2018;7(4):18. pmid:30159211

44. Hugo J, Chavane F, Beylerian M, Callet M, Denis D, Matonti F. Morphologic Analysis of Peripapillary Retinal Arteriole Using Adaptive Optics in Primary Open-angle Glaucoma. J Glaucoma. 2020;29(4):271–5. pmid:31977547

45. Wong TY, Knudtson MD, Klein R, Klein BE, Meuer SM, Hubbard LD. Computer-assisted measurement of retinal vessel diameters in the Beaver Dam Eye Study: methodology, correlation between eyes, and effect of refractive errors. Ophthalmology. 2004;111(6):1183–90. pmid:15177969

46. Gloster J, Parry DG. Use of photographs for measuring cupping in the optic disc. Br J Ophthalmol. 1974;58(10):850–62. pmid:4451651

47. Mitchell P, Smith W, Attebo K, Healey PR. Prevalence of open-angle glaucoma in Australia. The Blue Mountains Eye Study. Ophthalmology. 1996;103(10):1661–9. pmid:8874440

48. Takai Y, Tanito M, Omura T, Kawasaki R, Kawasaki Y, Ohira A. Comparisons of retinal vessel diameter and glaucomatous parameters between both eyes of subjects with clinically unilateral pseudoexfoliation syndrome. PLoS One. 2017;12(6):e0179663. pmid:28644852

49. Knudtson MD, Lee KE, Hubbard LD, Wong TY, Klein R, Klein BE. Revised formulas for summarizing retinal vessel diameters. Curr Eye Res. 2003;27(3):143–9. pmid:14562179

50. Littmann H. [Determination of the real size of an object on the fundus of the living eye]. Klin Monbl Augenheilkd. 1982;180(4):286–9. pmid:7087358

51. Chang M, Yoo C, Kim SW, Kim YY. Retinal vessel diameter, retinal nerve fiber layer thickness, and intraocular pressure in korean patients with normal-tension glaucoma. Am J Ophthalmol. 2011;151(1):100–5.e1. pmid:21094935

52. Rong X, Cai Y, Li M, Fang Y, Tian T, Pan Y. Peripapillary retinal artery in first diagnosed and untreated normal tension glaucoma. BMC Ophthalmol. 2019;19(1):203. pmid:31590635

53. Zheng Y, Cheung N, Aung T, Mitchell P, He M, Wong TY. Relationship of retinal vascular caliber with retinal nerve fiber layer thickness: the singapore malay eye study. Invest Ophthalmol Vis Sci. 2009;50(9):4091–6. pmid:19443726

54. Tham YC, Cheng CY, Zheng Y, Aung T, Wong TY, Cheung CY. Relationship between retinal vascular geometry with retinal nerve fiber layer and ganglion cell-inner plexiform layer in nonglaucomatous eyes. Invest Ophthalmol Vis Sci. 2013;54(12):7309–16. pmid:24135760

55. De Leon JM, Cheung CY, Wong TY, Li X, Hamzah H, Aung T, et al. Retinal vascular caliber between eyes with asymmetric glaucoma. Graefes Arch Clin Exp Ophthalmol. 2015;253(4):583–9. pmid:25795058

56. Lee KE, Klein BE, Klein R, Meuer SM. Association of retinal vessel caliber to optic disc and cup diameters. Invest Ophthalmol Vis Sci. 2007;48(1):63–7. pmid:17197517

57. Shin DH, Tsai CS, Parrow KA, Kim C, Wan JY, Shi DX. Intraocular pressure-dependent retinal vascular change in adult chronic open-angle glaucoma patients. Ophthalmology. 1991;98(7):1087–92. pmid:1679918

58. Nathanson JA, McKee M. Alterations of ocular nitric oxide synthase in human glaucoma. Invest Ophthalmol Vis Sci. 1995;36(9):1774–84. pmid:7543463

59. Flammer J, Orgül S, Costa VP, Orzalesi N, Krieglstein GK, Serra LM, et al. The impact of ocular blood flow in glaucoma. Prog Retin Eye Res. 2002;21(4):359–93. pmid:12150988

60. Jonas JB, Zäch FM, Gusek GC, Naumann GO. Pseudoglaucomatous physiologic large cups. Am J Ophthalmol. 1989;107(2):137–44. pmid:2913807

61. Chiba N, Omodaka K, Yokoyama Y, Aizawa N, Tsuda S, Yasuda M, et al. Association between optic nerve blood flow and objective examinations in glaucoma patients with generalized enlargement disc type. Clin Ophthalmol. 2011;5:1549–56. pmid:22125400

62. Nakazawa T, Fuse N, Omodaka K, Aizawa N, Kuwahara S, Nishida K. Different types of optic disc shape in patients with advanced open-angle glaucoma. Jpn J Ophthalmol. 2010;54(4):291–5. pmid:20700795

63. Hall JK, Andrews AP, Walker R, Piltz-Seymour JR. Association of retinal vessel caliber and visual field defects in glaucoma. Am J Ophthalmol. 2001;132(6):855–9. pmid:11730648

64. Papastathopoulos KI, Jonas JB. Follow up of focal narrowing of retinal arterioles in glaucoma. Br J Ophthalmol. 1999;83(3):285–9. pmid:10365034

65. Yokoyama Y, Aizawa N, Chiba N, Omodaka K, Nakamura M, Otomo T, 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. Clin Ophthalmol. 2011;5:1721–7. pmid:22205831

66. Aizawa N, Kunikata H, Shiga Y, Yokoyama Y, Omodaka K, Nakazawa T. Correlation between structure/function and optic disc microcirculation in myopic glaucoma, measured with laser speckle flowgraphy. BMC Ophthalmol. 2014;14:113. pmid:25252729

67. Omodaka K, Nakazawa T, Yokoyama Y, Doi H, Fuse N, Nishida K. Correlation between peripapillary macular fiber layer thickness and visual acuity in patients with open-angle glaucoma. Clin Ophthalmol. 2010;4:629–35. pmid:20668666

68. Liew G, Sharrett AR, Wang JJ, Klein R, Klein BE, Mitchell P, et al. Relative importance of systemic determinants of retinal arteriolar and venular caliber: the atherosclerosis risk in communities study. Arch Ophthalmol. 2008;126(10):1404–10. pmid:18852419

69. Nicolela MT, Walman BE, Buckley AR, Drance SM. Various glaucomatous optic nerve appearances. A color Doppler imaging study of retrobulbar circulation. Ophthalmology. 1996;103(10):1670–9. pmid:8874441

70. Mayama C, Araie M. Effects of antiglaucoma drugs on blood flow of optic nerve heads and related structures. Jpn J Ophthalmol. 2013;57(2):133–49. pmid:23321913

71. Hughes AD, Stanton AV, Jabbar AS, Chapman N, Martinez-Perez ME, Mc GTSA. Effect of antihypertensive treatment on retinal microvascular changes in hypertension. J Hypertens. 2008;26(8):1703–7. pmid:18622251

72. Li X, Wong WL, Cheung CY, Cheng CY, Ikram MK, Li J, et al. Racial differences in retinal vessel geometric characteristics: a multiethnic study in healthy Asians. Invest Ophthalmol Vis Sci. 2013;54(5):3650–6. pmid:23652487