Correction: Stereoscopic Analysis of Optic Nerve Head Parameters in Primary Open Angle Glaucoma: The Glaucoma Stereo Analysis Study

Purpose: The Glaucoma Stereo Analysis Study (GSAS), a cross sectional multicenter collaborative study, used a stereo fundus camera to assess various morphological parameters of the optic nerve head (ONH) in glaucoma patients and investigated the relationships between these parameters and patient characteristics. Subjects and Methods: The study included 187 eyes of 187 subjects with primary open angle glaucoma or normal tension glaucoma (male: female = 100: 87, age = 6169 years). Stereo pairs of ONH photographs were made with a stereo fundus camera (nonmyd WX). ONH morphological parameters were calculated with prototype analysis software. In addition to 35 standard parameters, we defined three novel parameters: disc tilt angle, rim decentering, and the absolute value of rim decentering. The correlation between each parameter and patient characteristics was analyzed with Spearman's rank correlation coefficient. Results: Patient characteristics included refractive error of 23.3863.75 diopters, intraocular pressure (IOP) of 13.662.6 mmHg, and visual field mean deviation (MD) of 24.7163.26 dB. Representative ONH parameters included a horizontal disc width of 1.6660.28 mm, vertical disc width of 1.8660.23 mm, disc area of 2.4260.63 mm 2 , cup area of 1.4560.57 mm 2 , and cup volume of 0.3160.22 mm 3. Correlation analysis revealed significant negative associations between vertical cup-to-disc ratio (0.8260.08) and MD (r = 20.40, P,0.01) and between disc tilt angle (10.5612.5 degrees) and refractive error (r = 20.36, P,0.01). Seventy-five percent of the eyes had a positive value for rim decentering (0.3060.42), indicating that rim thinning manifested more often as an inferior lesion than a superior lesion. Conclusion: We used stereoscopic analysis to establish a database of ONH parameters, which may facilitate future studies of glaucomatous changes in ONH morphology. Copyright: ß 2014 Yokoyama et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. All data are included within the Supporting Information files. Funding: The authors have no support or funding to report.


Introduction
Glaucoma is characterized by visual field defects that correspond to damaged areas of the optic nerve head (ONH). It affects over 70 million people worldwide, and is the second most common cause of blindness [1,2]. Currently, the only standard treatment to prevent progression of the most common form of the disease, open angle glaucoma (OAG), is maintenance of low intraocular pressure (IOP) with the use of medication or surgery. Although elevated IOP is generally recognized as a major risk factor for glaucoma [3]. It is well known that multiple factors are related to the development and progression of glaucoma. The progression of glaucomatous optic neuropathy (GON) is irreversible, making early diagnosis and treatment critical. Furthermore, assessment of a patient's ONH morphology is an essential part of the correct diagnosis and evaluation of glaucoma.
Generally, changes in ONH morphology, including thinning of the neuronal rim and enlargement of the ONH excavation, precede the progress of visual field defects. Morphological changes in the ONH are therefore considered important early biomarkers of GON and GON progression [4][5][6]. When identification of morphological changes is included in mass examinations and screenings, it has a positive effect on the early diagnosis of glaucoma [7,8]. Observation of the ONH remains a key part of follow-up care for glaucoma. Occurrences of rim notch, thinning of the local rim and enlargement of the ONH cup are especially important signs of risk [9]. In a previous study of glaucomatous ONH types using Nicolela's classification method, we found that the progression speed of glaucomatous visual field defects differed with ONH types [10]. Identifying certain patterns of glaucomatous ONH morphology may improve diagnostic accuracy and aid in identifying the effects of structural changes on visual function [11,12]. However, when the ONH is assessed with ophthalmoscopy or a subjective examination, ONH cupping and the severity of glaucoma may be underestimated because of the difficulty of stereoscopic viewing [13].
Topographic analysis with a simultaneous stereo fundus camera (nonmyd WX, Kowa Company, Ltd., Japan) is a noninvasive, noncontact imaging technique that does not require pupillary mydriasis. Its reliability has been demonstrated in a previous paper [14]._ENREF_16 It is therefore a promising tool for the assessment of ONH morphology. The Glaucoma Stereo Analysis Study (GSAS) is a multicenter study using this technique to assess various morphological parameters of the ONH in Japanese glaucoma patients. In this, the first report from GSAS, we established a database of various ONH parameters.

Subjects and Methods
This study (The Glaucoma Stereo Analysis Study: GSAS) was a cross sectional, multicenter collaborative study. It was approved by the Institutional Review Boards of the Tohoku University Graduate School of Medicine, Shimane University Faculty of Medicine, Fukui-ken Saiseikai Hospital, Sapporo Teishin Hospital, and St. Marianna University School of Medicine. All experimental procedures were conducted in accordance with the tenets set forth in the Declaration of Helsinki. All data collected from the institutions was analyzed anonymously.
One hundred and eighty-seven eyes of 187 patients with normal tension glaucoma (NTG) or primary open angle glaucoma (POAG) were recruited into this study from five institutions: Tohoku University Hospital, the Hospital of Shimane University Faculty of Medicine, Fukui-ken Saiseikai Hospital, Sapporo Teishin Hospital and the Hospital of St. Marianna University School of Medicine. The patients, whose ages ranged from 30 to 80 years, underwent full clinical ophthalmologic evaluations, including testing for visual acuity, refractive error, and intraocular pressure (IOP) with Goldmann applanation tonometry, as well as slit lamp and fundus examinations. At least one measurement of pre-treatment IOP (baseline IOP) was obtained retrospectively. Pre-surgical data on refractive error was also collected from eyes that had undergone refractive procedures such as cataract surgery. Visual field examinations with the Humphrey visual field analyzer (HFA; Carl Zeiss Meditec Inc., Dublin, California) were performed on all subjects within 6 months of recruitment (SITA standard, 30-2 or 24-2). Data from at least six HFA examinations performed over at least the previous three years was also collected retrospectively for each patient. Only reliable visual field data were used, i.e., from examinations with less than 20% false positives, less than 20% false negatives and less than 33% fixation losses. The mean deviation (MD) slope was calculated from these data.
Glaucoma diagnosis was based on the finding of glaucomatous visual field defects in reliable data from an HFA examination, with corresponding GON. GON was defined as an enlarged vertical cup-to-disc (C/D) ratio, narrow neuroretinal rim (rim) width, notching, and nerve fiber layer defects.
Additional inclusion criteria included: 1) best corrected visual acuity of 0.155 or better (LogMAR), 2) no congenital ONH anomalies, 3) ONH size within the typical normal range, defined as a disc-macula distance to disc diameter (DM/DD) ratio between approximately 2.4 and 3.0, 4) no clinically apparent secondary cause of glaucoma and no other disease affecting the visual field, 5) no history of intraocular surgery other than cataract or glaucoma surgery, 6) no history of cataract or glaucoma surgery in the previous three years, and 7) glaucomatous visual field loss more than 212 dB MD. If both eyes met the inclusion criteria, the eye with more advanced glaucoma was selected.
A summary of the inclusion criteria is shown in Table 1. Had no other disease affecting the visual field.
Had no conditions that might cause secondary glaucoma.
Had undergone automated perimetry with the Humphrey Field Analyzer (SITA standard 30-2 or 24-2) in the six months before the test day.
Had remaining visual fields of at least MD .212dB in their most recent tests.
Had been reliably tested with the Humphrey Field Analyzer (SITA standard 30-2) at least six times for more than three years (fixation loss ,33%, false positive and false negative ,20%).
Had visual field defects corresponding to glaucomatous optic disc change.
Had no history of intraocular surgery other than for glaucoma or cataracts.
Had no history of surgery for the previous three years.
?If both eyes met the inclusion criteria, the eye with more advanced glaucoma was selected.
BCVA: best corrected visual acuity DM/DD: disc-macula distance to disc diameter ratio MD: mean deviation doi:10.1371/journal.pone.0099138.t001 Analysis of optic nerve head topography Stereo fundus images of the ONH were obtained with a commercially available simultaneous stereo fundus camera (nonmyd WX). The nonmyd WX produces nonmydriatic fundus stereographs, as well as simultaneous right and left parallactic images, by using a single optical system to handle light paths in two directions [14]. _ENREF_16The built-in software (VK-2 WX, prototype version, Kowa Company, Ltd., Japan) automatically calculates ONH morphological 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.) while viewing the images stereoscopically. This determination was made according to the recommendations of the Japan Glaucoma Society Guidelines for Glaucoma, 3rd Edition [15,16]. The contour of the disc was delineated by the inner margin of Elschnig's scleral ring, and the contour of the cup was delineated by the outer margin of the cup, which was indicated by the bending of the ONH vessels at the rim. The observer determined several points on the contour (typically [8][9][10][11][12][13][14], and the contour line was then automatically generated by software spline interpolation. Parameters nonmyd WX included vertical C/D ratio, upper rim width, lower rim width, cup area, disc area, rim area, C/D area ratio, rim-to-disc (R/D) area ratio, sectional R/D ratio, cup volume, disc volume, rim volume, mean cup depth, maximum cup depth, height variation contour and disc damage likelihood scale (DDLS) stage. Measurements for the area and volume of the cup, disc and rim are illustrated in Figure S1. Depth value maps, necessary to determine parameters such as disc volume and mean cup depth, were generated based on the disparity between the right and left images of the stereo image pair with a stereo matching technique. Area and volume were calculated using three-dimensional analysis software (VK-2 WX) with correction for magnification. This correction, based a modification of Littman's method, was performed after entering the refractive error and corneal curvature of each eye into the software. The vertical and horizontal C/D ratios were defined by the ratio of the maximal vertical or horizontal diameter of the cup to maximal vertical or horizontal diameter of the disc. The aspect ratio was calculated by dividing the length of the largest diameter of the disc by the diameter perpendicular ( Figure S2). DDLS stage (9 stages; 0a, 0b, 1, 2, 3, 4, 5, 6, 7) is a diagnostic parameter proposed by Bayer et al. that provides an estimation of glaucomatous disc damage [17,18]. Analysis of the sectorial rim of the ONH was performed as shown in Figure S3.
In this study, we defined three novel parameters, disc tilt angle, rim decentering, and the absolute value of rim decentering. Rim decentering was calculated with the following formula: rim decentering = (superotemporal rim area 2 inferotemporal rim area)/(superotemporal rim area + inferotemporal rim area) ( Figure  S3). For the statistical analysis, absolute values for rim decentering were also determined. 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 ( Figure S4).

Statistical analysis
Statistical analysis was performed with JMP pro 10.02 (SAS Institute Inc.) for Windows. Continuous variables were expressed as mean values 6 standard deviation (SD). The Spearman's rank correlation coefficient was used to determine correlations between the patients' characteristics and the measured ONH parameters. In our analysis, ordinal data were treated as continuous data. The level of significance was set at 0.05 in all statistical tests.

Results
The detailed characteristics of the 187 subjects included in this study are shown in  Figure 1.
Average ONH topographic parameters are shown in Table 3. The average disc horizontal and vertical widths were 1.6660.28 mm and 1.8660.23 mm, respectively. The average disc aspect ratio was 1.1460.18. The average size of the disc cup was large (area: 1.4560.57 mm 2 , volume: 0.3160.22 mm 3 ) and the average rim of the ONH was thin (area: 0.9760.27 mm 2 , volume: 0.1760.10 mm 3 ) in the GON patients. The average DDLS stage (3.7760.95) indicated that the rim of the ONH had thinned and had defects. We found that 75% of subjects had a positive value (0.3060.42) for rim decentering, one of the new parameters in this study, indicating that damage to the ONH occurred at the lower side of the rim. The average disc tilt angle was 10.5612.5 degrees.
To investigate the relationship between ONH morphological parameters, the characteristics of the patients and visual field defects, we used Spearman's rank correlation coefficient (Table 4). Vertical C/D ratio and rim area were significantly correlated to visual field loss (MD) (r = 20.40, P,0.01, r = 0.40, P,0.01). Disc tilt angle was significantly correlated to age and spherical equivalent refractive error (r = 20.30, P,0.01, r = 20.36, P, 0.01). We also performed a multiple regression analysis, assigning tilt angle as the response variable and each of the latter two parameters as explanatory variables, to determine the strongest correlation between disc tilt angle, age and refractive error. The standardized partial regression coefficient (b value) for age was 2 0.17 (P = 0.03), and the variance inflation factor was 1.15. The b value for refractive error was 20.24 (P,0.01). We thus concluded that myopia had a stronger impact on disc tilt angle than age.

Discussion
The results of this study indicated that stereoscopically quantified parameters of ONH morphology (in particular, vertical C/D ratio, rim width and rim area) correlated to both the clinical background of glaucoma patients and to HFA MD quantifications of the degree of visual field damage. The reliability of the ONH measurement technique used in this study has already been amply demonstrated [14,19]. We believe that measurement of ONH morphological parameters with a simultaneous stereo fundus camera therefore holds promise in future hospital-based multicenter studies of disc morphology in patients with glaucoma.
Although previous studies have often used stereoscopic photography to analyze glaucoma, they have shown that photographic analysis of the ONH has a poor level of inter-observer consistency, even when the observers are experts [20,21]. This limitation prompted the development of clinical imaging technologies capable of objectively determining ONH morphology, such as Heidelberg retinal tomography (HRT). Meanwhile, however, the reliability of stereoscopic photography has continued to progress. Fundus cameras have become increasingly automated, and computer assisted digital analysis of the ONH has been introduced. Computer assisted techniques for analysis of ONH have improved intra-observer reproducibility and inter-observer consistency [22,23]. One of the latest devices, the nonmyd WX, can produce stereoscopic images with its built-in software and optional polarized filters. These enable examinations made by different observers to provide highly consistent results._ENRE-F_14_ENREF_25 Previously, histological studies have determined that the area of the ONH generally ranges from 2.48 to 2.75 mm 2 [24][25][26]. Additionally, population-based studies using fundus photography have found that the area of the ONH ranges from 2.09 to 2.94 mm 2 [27][28][29][30][31]. We found that in our patients, the mean disc area was 2.45 mm 2 , a result that was consistent with earlier reports. Past population-based studies using fundus photography have found that vertical C/D ratio ranges from 0.43 to 0.56 [27,28,30,32]. We found that mean cup area was 1.45 mm 2 and mean vertical C/D ratio was 0.82, results which indicated glaucomatous enlargement of the cup and which were correlated with MD, PSD, disc area and IOP. The result of 0.58 for C/D area ratio with average MD of 24.71 dB, however, was somewhat inconsistent with past reports. Nakatsue et al. used HRT to determine that the C/D area ratio in NTG patients with an average MD of 28.73 dB was 0.54, and that in POAG patients with an average MD of 29.70 dB it was 0.55 [33]. Additionally, Nouri-Mahdavi et al. used HRT to determine that in POAG patients with moderate glaucoma (average MD of 25.1 dB), mean C/D area ratio was 0.52. This inconsistency may be because of differences in the clinical backgrounds of the subjects and in methods of ONH assessment. Ocular morphological parameters can also be affected by many other factors, such as subjective bias, choice of measurement instrument, definition of parameters, ethnicity and sex [28,[34][35][36][37]._ENREF_34_ENREF_33 Moreover, the sample distribution of glaucoma types may have influenced the values of the ONH parameters. In a past study, we demonstrated the correlation between the severity of glaucomatous damage and different glaucomatous ONH types as classified by Nicolela's system [11,38]. Our report found that discs with generalized enlargement, i.e., those with a diffusely enlarged round cup without focal defects on the disc rim, had the highest correlation to MD. By contrast, eyes with the focal ischemic (FI) disc type had no correlation. FI discs have highly localized thinning of the neuronal rim, with other areas of the neuronal rim being normal [11]. Therefore, to assess glaucomatous deformation of the ONH, multiple approaches to investigate ONH morphology are needed.
To obtain results reflecting the nature of this morphology, we devised two novel measurement parameters. Rim decentering, one of the new parameters, was calculated with a formula that determined the difference in the rim ratio between the superior and inferior areas. The resulting data suggested that this ratio might help classify eyes with the FI disc type and identify regionally damaged areas in glaucomatous eyes. Seventy-five percent of patients in this study had a positive value for rim decentering, suggesting that the inferior rim of the ONH was more vulnerable to glaucomatous change than the superior. Furthermore, the angle of the inferior minimum rim disc ratio was 284 degrees, which was consistent with rim damage being more common in the inferotemporal disc region [39,40]._ENREF_27 The other new measurement parameter used in this study was the tilt angle of the disc. We found that the group of glaucoma patients included in our study showed characteristically steeply sloping ONHs, with an average tilt of 10.5 degrees and a tilt of more than 30 degrees in seven percent of patients. Existing studies use varying definitions of disc tilt. Tay et al. used the index of tilt, which was calculated as the quotient of the minimum diameter of the disc divided by its maximum. They defined ONHs as significantly tilted when they had an index of tilt less than or equal to 0.80 [41]._ENREF_42 However, other studies have set 0.70 or 0.75 as the criterion for significant tilt [42,43]._EN-REF_44 These previous methods of determining tilted discs relied on alternative measurement indices, but 3D photographs are capable of providing a direct value for tilt angle and can help in judging the severity of deformity. We found that tilt angle was moderately correlated to refractive error and age, and a multivariate analysis revealed that refractive error was a stronger indicator of a tilted disc than age. In Asia, the prevalence of myopia is significantly higher than in Western countries [44]. _ENREF_39Previous studies have reported that myopia and elongated axial length were associated with a higher prevalence of glaucoma [45][46][47][48][49]._ENREF_45_ENREF_45_ENREF_40 _EN-REF_29Morphologically, an elongated axial length laterally stretches the choroid and retina. As a result, the ONH tends to develop a myopic type tilted disc with crescent peripapillary atrophy. As our data did not include axial length, our analysis of myopic glaucomatous factors used refractive error, which is correlated with axial length. Nevertheless, we believe that quantification of disc tilt with the nonmyd WX may be able to identify a heightened risk of glaucoma in patients, and form a useful future part of clinical analyses of the ONH.
The stereo fundus camera technique used in this study had limits that may have affected our results. In a few patients, structures such as peripapillary atrophy can reduce the color contrast of Elschnig's scleral ring, rendering the image fuzzy and creating difficulties in determining the contour line of the disc edge. Nevertheless, even in such cases, a stereoscopic image provides more information than a monoscopic image, and should lead to more accurate prognoses. This study was also limited by being hospital-based and retrospective, although it included a relatively high number of patients.
In conclusion, we determined baseline ONH data with a simultaneous stereo fundus camera and identified a number of factors, including vertical C/D ratio, rim width and rim area, which were associated with glaucoma severity in the early and middle stages of the disease. Additionally, we found that tilted discs were correlated to spherical equivalent refractive error._ENRE-F_40_ENREF_39_ENREF_39_ENREF_40_ENREF_37_EN-REF_37_ENREF_37 We believe that quantitative data on ONH morphology is a powerful tool for clinical research into glaucoma.

Supporting Information
Figure S1 Schemata of measurements for the area and volume of the cup, disc and rim. The reference plane was defined as the average height of the cup contour. The zero-mm plane was defined as the average height of the nasal retinal surface outside the disc area. Contour lines for the disc and cup were based on observer-determined points on the fundus photograph (typically 8-14,) with computer-generated spline interpolation used to generate the final curve. The area and volume of the cup, disc and rim were determined with these two contour lines and the two planes, as illustrated. The superior and the inferior rim widths were measured on the vertical axis of the optic nerve head. In the cup area, the maximum depth and the mean depth were calculated. The height variation contour was calculated by subtracting the minimum height of the disc contour line from its maximum height. Maximum height and minimum height of the depth map were defined as the maximum height and minimum height of the measurement area. (TIF) Figure S2 Schema of the aspect ratio. The aspect ratio was calculated by dividing the length of the largest diameter (the major axis) of the optic disc by the perpendicular diameter (the minor axis). (TIF) Figure S3 Illustration of rim sections and rim decentering. a) Definitions of the superior and inferior sections of the rim area. The superior section for measurements ranged from 60u to 120u and the inferior section ranged from 240u to 300u. These sections were used to calculate the minimum rim-to-disc ratio and its angle. b) Definitions of the six sections of the rim area. The six sections of the rim were defined as follows: section 1 ranged from 0u to 45u and from 315u to 360u, section 2 ranged from 45u to 90u, section 3 ranged from 90u to 135u, section 4 ranged from 135u to 225u, section 5 ranged from 225u to 270u, and section 6 ranged from 270u to 315u. These six sections were used to calculate the average of the sectional minimum rim-to-disc ratios. Rim decentering was determined with the following formula: rim decentering = (superotemporal rim area [2] 2 inferotemporal rim area [6])/(superotemporal rim area [2] + inferotemporal rim area [6]). (TIF) Figure S4 Definition of disc tilt angle. The disc tilt angle was defined as the degree of the angle between the plane horizontal to the observer and the line between 0u and 180u on the disc edge. T: Temporal, N: Nasal, h: disc tilt angle. (TIF)