Open Access
Glaucoma  |   September 2017
Vertical Macular Asymmetry Measures Derived From SD-OCT for Detection of Early Glaucoma
Author Affiliations & Notes
  • Farideh Sharifipour
    Glaucoma Division, Stein Eye Institute, David Geffen School of Medicine, University of California Los Angeles, California, United States
    Department of Ophthalmology, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, Iran
  • Esteban Morales
    Glaucoma Division, Stein Eye Institute, David Geffen School of Medicine, University of California Los Angeles, California, United States
  • Ji Woong Lee
    Glaucoma Division, Stein Eye Institute, David Geffen School of Medicine, University of California Los Angeles, California, United States
    Department of Ophthalmology, Pusan National University College of Medicine, Busan, South Korea
  • JoAnn Giaconi
    Glaucoma Division, Stein Eye Institute, David Geffen School of Medicine, University of California Los Angeles, California, United States
  • Abdelmonem A. Afifi
    Department of Biostatistics, Jonathan and Karin Fielding School of Public Health, University of California Los Angeles, Los Angeles, California, United States
  • Fei Yu
    Glaucoma Division, Stein Eye Institute, David Geffen School of Medicine, University of California Los Angeles, California, United States
    Department of Biostatistics, Jonathan and Karin Fielding School of Public Health, University of California Los Angeles, Los Angeles, California, United States
  • Joseph Caprioli
    Glaucoma Division, Stein Eye Institute, David Geffen School of Medicine, University of California Los Angeles, California, United States
  • Kouros Nouri-Mahdavi
    Glaucoma Division, Stein Eye Institute, David Geffen School of Medicine, University of California Los Angeles, California, United States
  • Correspondence: Kouros Nouri-Mahdavi, Stein Eye Institute, 100 Stein Plaza, Los Angeles, CA 90095, USA; nouri-mahdavi@jsei.ucla.edu
Investigative Ophthalmology & Visual Science September 2017, Vol.58, 4310-4317. doi:https://doi.org/10.1167/iovs.17-21961
  • Views
  • PDF
  • Share
  • Tools
    • Alerts
      ×
      This feature is available to authenticated users only.
      Sign In or Create an Account ×
    • Get Citation

      Farideh Sharifipour, Esteban Morales, Ji Woong Lee, JoAnn Giaconi, Abdelmonem A. Afifi, Fei Yu, Joseph Caprioli, Kouros Nouri-Mahdavi; Vertical Macular Asymmetry Measures Derived From SD-OCT for Detection of Early Glaucoma. Invest. Ophthalmol. Vis. Sci. 2017;58(10):4310-4317. https://doi.org/10.1167/iovs.17-21961.

      Download citation file:


      © ARVO (1962-2015); The Authors (2016-present)

      ×
  • Supplements
Abstract

Purpose: To test the hypothesis that vertical asymmetry in macular ganglion cell/inner plexiform layer (GCIPL) thickness can improve detection of early glaucoma.

Methods: Sixty-nine normal eyes and 101 glaucoma eyes had macular imaging with spectral-domain optical coherence tomography (OCT; 200 × 200 cube). The resulting GCIPL thickness matrix was grouped into a 20 × 20 superpixel array and superior superpixels were compared to their inferior counterparts. A global asymmetry index (AI) was defined as the grand mean of the asymmetry ratios. To measure local asymmetry, the corresponding thickness measurements of three rows above and below the horizontal raphe were compared individually and in combinations. Global and local AIs were compared to the best-performing GCIPL thickness parameters with area under the receiver operating curves (AUC) and sensitivity/specificities.

Results: Age or axial length did not influence AIs in normal subjects (P ≥ 0.08). Global and local AIs were significantly higher in the glaucoma group compared to normal eyes. Minimum (AUC = 0.962, 95% confidence interval [CI]: 0.936–0.989) and inferotemporal thickness (AUC = 0.944, 95% CI: 0.910–0.977; P = 0.122) performed best for detection of early glaucoma. The AUC for global AI was 0.851 (95% CI: 0.792–0.909) compared to 0.916 (95% CI: 0.874–0.958) for the best local AI. Combining minimum or inferotemporal GCIPL thickness and the best local AI led to higher partial AUCs (0.088 and 0.085, 90% specificity, P = 0.120 and 0.130, respectively) than GCIPL thickness measures.

Conclusions: Macular vertical thickness asymmetry measures did not perform better than sectoral or minimum GCIPL thickness for detection of early glaucoma. Combining local asymmetry parameters with the best sectoral GCIPL thickness measures enhanced this task.

Glaucoma is characterized by progressive loss of retinal ganglion cells and their axons. There is a large body of evidence suggesting that structural changes can be detected several years before measurable visual loss.18 Improvements in the resolution of OCT has enabled us to measure thickness of retinal layers with a precision of 3 to 6 μm.9 Performance of regional macular ganglion cell/inner plexiform layer (GCIPL) thickness measurements is comparable to that of peripapillary retinal nerve fiber layer (ppRNFL) thickness for the detection of early glaucoma.1012 
Early attenuation of retinal nerve fiber layer or GCIPL thinning in glaucoma is usually greater on one side of the horizontal meridian, most frequently inferiorly.13 The asymmetry between threshold sensitivities of superior and inferior hemifields on standard automated perimetry is the basis for the glaucoma hemifield test (GHT)14 used for detection of early glaucomatous field loss. Several studies have evaluated if the vertical asymmetry in full macular or GCIPL thicknesses measured by OCT can serve as a biomarker for early glaucoma.1523 Most have used OCT posterior pole asymmetry analysis (PPAA; Spectralis; Heidelberg Engineering, Dossenheim, Germany) program or global or local measures of vertical asymmetry and showed significant vertical asymmetry in glaucoma patients compared to normal individuals.1722 
The goal of the current study was to test the hypothesis that biomarkers of vertical macular thickness asymmetry based on GCIPL measurements can improve discrimination of early perimetric glaucoma from normal subjects. 
Materials and Methods
We recruited normal subjects and glaucoma patients from ongoing studies at Stein Eye Institute that were all approved by the institutional review board at the University of California Los Angeles (UCLA). This cross-sectional case-control study adhered to the tenets of Declaration of Helsinki. Patients with open-angle glaucoma were enrolled in the study if they met the following criteria: age ≥30 years; open angles; visual acuity equal or better than 20/80; refractive error less than 8.0 diopters (D); and astigmatism ≤3 D. Clinical diagnosis of glaucoma was initially made by a glaucoma specialist at Stein Eye Institute's Glaucoma Division based on the presence of characteristic optic nerve head changes detected on slit lamp funduscopy and/or 3-D disc photographs, presence of RNFL defects on dilated funduscopy, and/or characteristic visual field defects on SAP. The enrolled eyes were later confirmed to have changes consistent with glaucomatous optic neuropathy on review of stereoscopic disc photographs; all eyes had evidence of early glaucomatous field loss (see below). Eyes with evidence of retinal or neurologic diseases were excluded. Normal subjects were recruited by advertising on UCLA's campus-wide website, placing fliers in the clinics, and soliciting spouses and friends of patients examined at the Stein Eye Institute's Glaucoma Clinic. The enrolled normal subjects were required to have open angles, corrected visual acuity of 20/25 or better, normal eye examination, normal visual fields, and lacked definitive evidence of glaucomatous damage at the level of the optic nerve head on review of stereoscopic optic disc photographs by two clinicians (KNM, JAG). 
All subjects underwent a thorough eye examination on the day of imaging, including visual acuity, automated refraction, measurement of intraocular pressure (IOP); gonioscopy; slit-lamp exam; dilated fundus exam; and standard achromatic perimetry (SAP). Axial length and keratometry were measured with an optical biometer (IOLMaster; Carl Zeiss Meditec, Inc., Jena, Germany). Stereoscopic disc photographs and 200 × 200 macular cubes (Cirrus high-definition [HD]-OCT; Carl Zeiss Meditec, Inc.) were obtained after pupillary dilation. All the images were reviewed afterward by one of the investigators, and images with signal strength <7, obvious motion artifact, or incorrect segmentation were excluded. 
A total of 69 eyes from 69 normal subjects and 101 eyes from 101 early glaucoma patients were included in the study. Eyes with visual field mean deviation (MD) equal to or better than −6 dB were considered to have early glaucoma.24 
Imaging Methods
The 200 × 200 macular cube (Carl Zeiss Meditec, Inc.) consists of 40,000 A-scans in a 6 × 6 mm grid (in an emmetropic eye with axial length of 24.46 mm) centered on the fovea. The raw images were exported and the macular GCIPL layers were segmented by custom software at Carl Zeiss Meditec headquarters and the GCIPL thickness measurements were provided to the investigators in 200 × 200 matrices. The data matrices were centered in relation to the foveal center. Thickness data of the left eyes were flipped to convert them to the right eye format. The 40,000 A-scan measurements were then clustered into a 20 × 20 grid of superpixels (Fig. 1), and various approaches were explored to compare the superior superpixels to their inferior counterparts. The global and sectoral GCIPL thickness measurements were also exported as extensible markup language files. 
Figure 1
 
The 200 × 200 macular cube (Carl Zeiss Meditec, Inc.) consists of 40,000 A-scan measurements that were clustered into a 20 × 20 grid of superpixels centered on the fovea. Local asymmetry measures were defined based on comparison of the GCIPL thickness on the three corresponding pairs of rows above and below the temporal horizontal raphe including four or eight temporal columns, both individually and in combination. The best performing local asymmetry parameter was the one incorporating the three rows above and below the temporal raphe consisting of eight superpixels.
Figure 1
 
The 200 × 200 macular cube (Carl Zeiss Meditec, Inc.) consists of 40,000 A-scan measurements that were clustered into a 20 × 20 grid of superpixels centered on the fovea. Local asymmetry measures were defined based on comparison of the GCIPL thickness on the three corresponding pairs of rows above and below the temporal horizontal raphe including four or eight temporal columns, both individually and in combination. The best performing local asymmetry parameter was the one incorporating the three rows above and below the temporal raphe consisting of eight superpixels.
We used four different methods to compare each pair of superior and inferior superpixels and to obtain an asymmetry index (AI) as below: 
  1.  
    superior/inferior thickness
  2.  
    |log (superior/inferior thickness)|
  3.  
    |superior-inferior thickness/max (superior, inferior)|
  4.  
    |superior- inferior thickness|
Global AI was defined as the grand mean of the asymmetry values in the central circular area (6-mm diameter).22 To evaluate local asymmetry, we compared AIs based on three corresponding pairs of rows above and below temporal horizontal line singly and in combinations using the above-mentioned methods including eight or four superpixels on each row temporal to the fovea (Fig. 1). This area is where early glaucomatous damage commonly occurs and roughly corresponds to nasal visual field loss.25 
Statistical Analyses
Sidak or Mann-Whitney U, and χ2 tests were used for descriptive analyses and for parameter comparisons between groups. Univariate linear regression analysis was used to explore the relationship between AI and age, sex, and axial length. In addition to global asymmetry, performance of local asymmetry measures (with the aforementioned four methods) was explored comparing the three corresponding pairs of rows above and below the temporal horizontal line. Areas under the receiver operating curve (AUC) after logistic regression were estimated to compare performance of global and local asymmetry biomarkers to the best-performing GCIPL parameters for glaucoma detection. Partial AUC (pAUC) at 90% specificity was also compared between groups. All global and local asymmetry analyses used absolute values of the asymmetry indices with methods 2 through 4. For analyses using method 1, if the inferior superpixel was thinner than the corresponding superior superpixel, 1 / AI was used. 
Statistical analyses were performed with the commercial software (Stata, version 12.1; StataCorp, College Station TX, USA) and R software (version 3.3.2). All analyses included took into account the difference in age between the two groups. Values of P less than 0.05 were considered as significant. 
Results
We enrolled 69 normal eyes (69 subjects) and 101 early glaucoma eyes (101 patients). Table 1 summarizes the demographic and clinical characteristics of the study sample. The normal group was significantly younger than the glaucoma group (P < 0.001). Regression analysis showed no significant relationship between global or local asymmetry measured by any method and age (P = 0.445 and 0.192, respectively for method 2); sex (P = 0.652 and 0.840 for method 2); or axial length (P = 0.08 and 0.08 for method 2) in normal individuals. 
Table 1
 
Demographic and Clinical Findings of Subjects Enrolled in the Study
Table 1
 
Demographic and Clinical Findings of Subjects Enrolled in the Study
The glaucoma asymmetry index with method 2 (|log (superior/inferior thickness)|) performed best among all asymmetry indices and therefore, all analyses presented here used this method. The global AI was significantly higher in glaucoma group compared to normal eyes (Table 2). Among the local asymmetry methods, the index combining the three rows above versus the 3 rows below the temporal horizontal line demonstrated the best performance (Table 3). All local asymmetry methods using eight temporal superpixels in each row performed significantly better than methods including only four temporal superpixels in the analysis (P < 0.001). 
Table 2
 
Comparison of Sectoral GCIPL Thicknesses and Local and Global AIs Between Study Groups
Table 2
 
Comparison of Sectoral GCIPL Thicknesses and Local and Global AIs Between Study Groups
Table 3
 
Values of AUC for GCIPL Thickness and Local and Global AIs
Table 3
 
Values of AUC for GCIPL Thickness and Local and Global AIs
The best-performing local GCIPL thickness parameters for detection of early glaucoma were the minimum GCIPL (AUC = 0.962, 95% confidence interval [CI]: 0.936–0.989) and inferotemporal (AUC = 0.944, 95% CI: 0.910–0.977; P = 0.122) for the difference between the two; Table 3). The calculated AUCs for global AIs by all four methods were significantly lower than that of local AIs and sectoral GCIPL thicknesses, (P = 0.007, Table 3; Fig. 2). The area under the receiver operating curve for the best local AI (i.e., the one combining 3 rows above and below the horizontal raphe) was significantly smaller than that of the minimum GCIPL (P = 0.008), but was not significantly different from the AUC for inferotemporal GCIPL thickness (P = 0.158, Table 3). 
Figure 2
 
Receiver operating characteristic curves comparing glaucoma diagnostic capabilities of the best local asymmetry index, global asymmetry index, inferotemporal, and minimal GCIPL thicknesses, and combination of the best local asymmetry and regional GCIPL thickness measures.
Figure 2
 
Receiver operating characteristic curves comparing glaucoma diagnostic capabilities of the best local asymmetry index, global asymmetry index, inferotemporal, and minimal GCIPL thicknesses, and combination of the best local asymmetry and regional GCIPL thickness measures.
Combining inferotemporal sector and the best local AI significantly increased the diagnostic performance of inferotemporal GCIPL thickness in early glaucoma (AUC: 0.960 compared to AUC = 0.944 for inferotemporal sector, P = 0.04, Table 3) and improved the sensitivity of the combined parameter at 90% specificity from 71% to 93%. Adding the best local asymmetry to minimum GCIPL thickness led only to a small nonsignificant increase in the AUC from 0.962 to 0.965 (P = 0.392, Table 3). 
Global AI also had the lowest pAUC as compared to minimum or inferotemporal GCIPL, and best local AI (pAUC = 0.024 vs. 0.085, 0.081, and 0.065, respectively, P < 0.001 for all pairwise differences between global AI and other parameters). A combination of minimum or inferotemporal GCIPL thickness and the best local AI had higher pAUCs (0.088 and 0.085, respectively) than all individual parameters including minimum GCIPL; however, the difference did not reach statistical significance (P = 0.12 and 0.13, respectively). 
Discussion
Our study showed that the performance of vertical GCIPL thickness asymmetry along the horizontal temporal line (based on various global or local criteria) was not as good as the best regional GCIPL thickness measures (minimum GCIPL or inferotemporal GCIPL sector) for detection of early glaucoma; however, combining the best-performing local asymmetry parameter with the inferotemporal sector GCIPL thickness improved performance of macular OCT measures based on AUCs. However, based on the pAUCs at ≥90% specificity cutoff (where it matters most clinically), the combination of local asymmetry and either the inferotemporal or minimal GCIPL thickness had a better ability to discriminate glaucomatous from normal eyes than either thickness measures alone although this difference did not reach statistical significance. 
The advent of SD-OCT technology has resulted in significant research endeavors to help diagnose glaucoma at an earlier stage before functional loss occurs. In line with this strategy, ppRNFL thickness and various inner macular thickness measurements have been used to improve glaucoma detection. Several studies have shown that measurements of inner retinal layer (IRL) thickness in localized sectors have a better diagnostic performance compared to average macular IRL or full thickness and are comparable to average ppRNFL thickness for detection of early glaucoma.10,11,2633 
Our study showed that regional GCIPL thickness measures outperformed global and local vertical asymmetry outcome measures. In agreement with some previous reports, our study revealed that a global asymmetry index, which is the average of all asymmetry values across the horizontal meridian, had the worst performance for discrimination of early glaucoma from normal controls.20 Glaucoma frequently starts as a localized thinning in retinal tissue, which might go undetected by averaging asymmetry measures across the macula. Local asymmetry measures could actually detect regional changes across the temporal horizontal line and performed significantly better than global asymmetry measures in this study. Neither the global index nor local asymmetry indices were affected by age, sex, or axial length and were significantly higher in glaucoma group compared to normal individuals. 
There are some reports indicating that intraeye asymmetry between superior and inferior total macular or inner layer thicknesses has equal or greater ability than ppRNFL and ganglion cell complex (GCC) measurements for detection of early glaucoma,1522 although the method used for estimation of asymmetry has varied among such studies. Asrani et al.17 were the first to customize the macular thickness map of the Spectralis SD-OCT to an 8 × 8 grid centered on the fovea and positioned symmetrically to the fovea-disc axis. Intra- and intereye cell-to-cell comparisons were made using a cutoff value of ≥30 μm and the difference presented as a grayscale grid. This was later implemented as the PPAA program on Spectralis SD-OCT.17 Seo et al.18 showed that the diagnostic ability of PPAA for localized RNFL defects was comparable to ppRNFL thickness measurements. In a study by Um et al.19 exploring the Spectralis posterior pole retinal thickness map, hemifield asymmetry in macular full thickness measurements was evaluated in five corresponding superior and inferior zones similar to GHT in perimetry. The macular hemifield test had significantly higher sensitivity and similar specificity compared to average ppRNFL thickness in early glaucoma, but showed comparable sensitivities and specificities in glaucoma suspects and advanced glaucoma. Diagnostic performance of the macular hemifield test was comparable to sectoral ppRNFL thickness measurements. Sullivan-Mee and colleagues20 studied inter- and intraeye asymmetry of ppRNFL and macular thickness measurements on Spectralis PPAA for detection of early glaucoma. Intereye asymmetry parameters had the highest diagnostic sensitivities, although RNFL thickness had better overall performance than total macular thickness. Intra-eye macular thickness asymmetry had the lowest sensitivity and lowest AUC of all studied parameters. The authors reported optimal cutoff values of 5 and 9 μm for inter-eye and intra-eye macular thickness asymmetry, respectively.20 
Yamashita and colleagues21 demonstrated that the symmetry of retinal thickness between upper and lower cells was high in the central and temporal-macular areas but not in the peripheral and nasal macular areas with a range of difference of 3.1 to 23.2 μm. This might be due to lower reproducibility of thickness measurements temporal to the optic disc.34 In a study of preperimetric glaucoma eyes with only inferior RNFL defects, significant asymmetry in macular thickness was reported compared to normal eyes.35 Macular symmetry was measured as (inferior thickness/ superior thickness) ×100 on images from Topcon SD-OCT. The authors also reported significantly thinner inferior macular thickness in normal eyes.35 Yamada et al.22 showed that asymmetry between upper and lower macular hemifield ganglion cell layer thickness had excellent diagnostic performance for early glaucoma compared to full thickness or inner retinal layer thicknesses in Japanese patients. They devised an asymmetry index based on 10 vertically oriented B-scans performed on the Spectralis OCT according to this formula:  
\(\def\upalpha{\unicode[Times]{x3B1}}\)\(\def\upbeta{\unicode[Times]{x3B2}}\)\(\def\upgamma{\unicode[Times]{x3B3}}\)\(\def\updelta{\unicode[Times]{x3B4}}\)\(\def\upvarepsilon{\unicode[Times]{x3B5}}\)\(\def\upzeta{\unicode[Times]{x3B6}}\)\(\def\upeta{\unicode[Times]{x3B7}}\)\(\def\uptheta{\unicode[Times]{x3B8}}\)\(\def\upiota{\unicode[Times]{x3B9}}\)\(\def\upkappa{\unicode[Times]{x3BA}}\)\(\def\uplambda{\unicode[Times]{x3BB}}\)\(\def\upmu{\unicode[Times]{x3BC}}\)\(\def\upnu{\unicode[Times]{x3BD}}\)\(\def\upxi{\unicode[Times]{x3BE}}\)\(\def\upomicron{\unicode[Times]{x3BF}}\)\(\def\uppi{\unicode[Times]{x3C0}}\)\(\def\uprho{\unicode[Times]{x3C1}}\)\(\def\upsigma{\unicode[Times]{x3C3}}\)\(\def\uptau{\unicode[Times]{x3C4}}\)\(\def\upupsilon{\unicode[Times]{x3C5}}\)\(\def\upphi{\unicode[Times]{x3C6}}\)\(\def\upchi{\unicode[Times]{x3C7}}\)\(\def\uppsy{\unicode[Times]{x3C8}}\)\(\def\upomega{\unicode[Times]{x3C9}}\)\(\def\bialpha{\boldsymbol{\alpha}}\)\(\def\bibeta{\boldsymbol{\beta}}\)\(\def\bigamma{\boldsymbol{\gamma}}\)\(\def\bidelta{\boldsymbol{\delta}}\)\(\def\bivarepsilon{\boldsymbol{\varepsilon}}\)\(\def\bizeta{\boldsymbol{\zeta}}\)\(\def\bieta{\boldsymbol{\eta}}\)\(\def\bitheta{\boldsymbol{\theta}}\)\(\def\biiota{\boldsymbol{\iota}}\)\(\def\bikappa{\boldsymbol{\kappa}}\)\(\def\bilambda{\boldsymbol{\lambda}}\)\(\def\bimu{\boldsymbol{\mu}}\)\(\def\binu{\boldsymbol{\nu}}\)\(\def\bixi{\boldsymbol{\xi}}\)\(\def\biomicron{\boldsymbol{\micron}}\)\(\def\bipi{\boldsymbol{\pi}}\)\(\def\birho{\boldsymbol{\rho}}\)\(\def\bisigma{\boldsymbol{\sigma}}\)\(\def\bitau{\boldsymbol{\tau}}\)\(\def\biupsilon{\boldsymbol{\upsilon}}\)\(\def\biphi{\boldsymbol{\phi}}\)\(\def\bichi{\boldsymbol{\chi}}\)\(\def\bipsy{\boldsymbol{\psy}}\)\(\def\biomega{\boldsymbol{\omega}}\)\(\def\bupalpha{\bf{\alpha}}\)\(\def\bupbeta{\bf{\beta}}\)\(\def\bupgamma{\bf{\gamma}}\)\(\def\bupdelta{\bf{\delta}}\)\(\def\bupvarepsilon{\bf{\varepsilon}}\)\(\def\bupzeta{\bf{\zeta}}\)\(\def\bupeta{\bf{\eta}}\)\(\def\buptheta{\bf{\theta}}\)\(\def\bupiota{\bf{\iota}}\)\(\def\bupkappa{\bf{\kappa}}\)\(\def\buplambda{\bf{\lambda}}\)\(\def\bupmu{\bf{\mu}}\)\(\def\bupnu{\bf{\nu}}\)\(\def\bupxi{\bf{\xi}}\)\(\def\bupomicron{\bf{\micron}}\)\(\def\buppi{\bf{\pi}}\)\(\def\buprho{\bf{\rho}}\)\(\def\bupsigma{\bf{\sigma}}\)\(\def\buptau{\bf{\tau}}\)\(\def\bupupsilon{\bf{\upsilon}}\)\(\def\bupphi{\bf{\phi}}\)\(\def\bupchi{\bf{\chi}}\)\(\def\buppsy{\bf{\psy}}\)\(\def\bupomega{\bf{\omega}}\)\(\def\bGamma{\bf{\Gamma}}\)\(\def\bDelta{\bf{\Delta}}\)\(\def\bTheta{\bf{\Theta}}\)\(\def\bLambda{\bf{\Lambda}}\)\(\def\bXi{\bf{\Xi}}\)\(\def\bPi{\bf{\Pi}}\)\(\def\bSigma{\bf{\Sigma}}\)\(\def\bPhi{\bf{\Phi}}\)\(\def\bPsi{\bf{\Psi}}\)\(\def\bOmega{\bf{\Omega}}\)\begin{equation}{\rm{asymmetry\ index}} = \left| {{\rm{log}}_{10}\left( {{\rm{lower\ hemiretinal\ thickness}}\over {\rm{upper\ hemiretinal\ thickness}}} \right)} \right|{\rm {.}}\end{equation}
Interestingly, the performance of this asymmetry index did not change with severity of glaucoma. Similarly, Kim et al.23 reported that a local asymmetry approach utilizing differences in GCIPL thickness within 10 pixels above and below temporal horizontal raphe performed better than the best sectoral GCIPL thickness for detection of early glaucoma in a group of Korean patients. In a similar approach to our study, Hwang et al.36 reported good diagnostic ability of asymmetry analysis in early glaucoma.  
We explored various definitions of vertical asymmetry both globally and regionally and compared their diagnostic performance to each other and also to the best sectoral GCIPL thickness. Interestingly, we observed that the simple difference between superior and inferior macular thickness (which is currently in use by Spectralis PPAA) had the worst diagnostic performance among the four methods of asymmetry explored. Other methods had similar AUCs for detection of early glaucoma, although method 2 (similar to the method presented by Yamada et al.) was the best. The contradictory results reported by different studies may be explained by the limitations inherent to imaging devices. Spectralis SD-OCT provides higher resolution and possibly more precise segmentation compared to 200 × 200 macular cube images of Cirrus HD-OCT, but uses the fovea-disc axis as the reference line for vertical asymmetry in PPAA, as opposed to the horizontal meridian in Cirrus OCT; some investigators such as Yamada and colleagues22 used a specific imaging protocol based on vertical lines perpendicular to the horizontal meridian. The differences in study populations in terms of ethnicity and type of glaucoma may also contribute to disagreement among studies since Asian studies include more normal-tension glaucoma (NTG). 
Our study is the first to explore the performance of various global and local GCIPL-based temporal asymmetry measures in a diverse ethnic group of patients in the United States. One reason for the less optimal performance of local or global asymmetry measures in our study sample could be the lower prevalence of NTG in the US population. Only 15% of our study patients were diagnosed as NTG. Normal-tension glaucoma, which results in deep localized defects in central visual field and macular region early in the disease course,37 could result in significant vertical asymmetry before marked retinal thickness changes are detectable. Therefore, it is possible that the performance of vertical asymmetry measures may depend on the underlying diagnosis and/or population from which the study sample originates. This may explain why Yamada et al.22 and Kim and colleagues23 found asymmetry measures to perform much better compared to total or inner retinal thickness measurements for glaucoma detection. 
We observed that combining inferotemporal GCIPL thickness and local asymmetry measures significantly improved the diagnostic performance compared to the inferotemporal sector alone (P = 0.04, Table 3; Fig. 2). Despite the fact that both measures already had large AUCs, their combination led to a significant increase in glaucoma diagnostic capability similar to that of minimum GCIPL. Not every device provides the minimum GCIPL parameter and hence, combining inferotemporal inner macular thickness and local asymmetry measures would actually be an equivalent surrogate for this parameter. Additionally, when pAUCs were used to compare the performance of various thickness parameters above 90% specificity, not only did the best local vertical asymmetry measure have an equal performance to minimum GCIPL thickness, but also combining the asymmetry measures with either minimum or GCIPL thickness further enhanced performance at high specificity where it is most relevant for early detection of glaucoma—although this difference did not reach statistical significance. 
One of the limitations of our study is that the macular images were not specifically designed to identify vertical asymmetry across the temporal horizontal line and hence the resolution of the images was likely less than optimal. It is possible that the performance of macular asymmetry indices will improve as customized algorithms are more widely adopted and the quality of OCT images improves. We assumed that the temporal raphe is approximately aligned with the horizontal line extending temporally from the center of the fovea. It has been shown that temporal raphe greatly varies among individuals and that the raphe-fovea-disc line looks more like a broken line with an average angle of 170°. However, currently, there is no easy way to exactly define the location of the temporal raphe on commercially available OCT images. We were able to investigate only GCIPL thickness asymmetry; however, other than the report by Yamada et al.,22 there is no strong evidence in the literature for the superiority of the ganglion cell layer or GCC measurements to GCIPL thickness measures. In a recent study, Miraftabi and colleagues38 found that GCIPL actually demonstrated the strongest structure-function relationship among all macular parameter with central visual fields. Our control group was significantly younger than the glaucoma group. To address this age was entered as a covariate in all models. Our study did not include angle closure glaucoma eyes. This potentially limits the generalizability of our study. 
Conclusions
Our study demonstrated that temporal vertical asymmetry in macular GCIPL thickness was comparable to sectoral inferotemporal GCIPL thickness for detection of early glaucoma although minimal GCIPL had the best nominal AUC. The performance of the global asymmetry index was worse than all local vertical asymmetry indices and the best sectoral GCIPL thickness measures. A combination of sectoral GCIPL thickness and the best local vertical asymmetry parameter significantly improved the diagnostic performance of macular images and may be implemented in OCT devices to enhance detection of early glaucoma. 
Acknowledgments
Presented as a poster at the annual meeting of the Association for Research in Vision and Ophthalmology, Denver, Colorado, United States, May 2015. 
Supported by an National Institutes of Health Mentored Patient-oriented Research Career Development Award (5K23EY022659; KNM) and an unrestricted Departmental Grant from Research to Prevent Blindness. 
Disclosure: F. Sharifipour, None; E. Morales, None; J.W. Lee, None; J. Giaconi, None; A.A. Afifi, None; F. Yu, None; J. Caprioli, None; K. Nouri-Mahdavi, Carl Zeiss Meditec (F), Heidelberg Engineering (F) 
References
Sommer A, Miller NR, Pollack I, Maumenee AE, George T. The nerve fiber layer in the diagnosis of glaucoma. Arch Ophthalmol. 1977; 95: 2149–2156.
Sommer A, Katz J, Quigley HA, et al. Clinically detectable nerve fiber atrophy precedes the onset of glaucomatous field loss. Arch Ophthalmol. 1991; 109: 77–83.
Motolko M, Drance SM. Features of the optic disc in preglaucomatous eyes. Arch Ophthalmol. 1981; 99: 1992–1994.
Tuulonen A, Lehtola J, Airaksinen PJ. Nerve fiber layer defects with normal visual fields. Do normal optic disc and normal visual field indicate absence of glaucomatous abnormality? Ophthalmology. 1993; 100: 587–597; discussion 597–598.
Kass MA, Heuer DK, Higginbotham EJ, et al. The Ocular Hypertension Treatment study: a randomized trial determines that topical ocular hypotensive medication delays or prevents the onset of primary open-angle glaucoma. Arch Ophthalmol 2002; 120: 701–713; discussion 829–830.
Subbiah S, Sankarnarayanan S, Thomas PA, Nelson Jesudasan CA . Comparative evaluation of optical coherence tomography in glaucomatous, ocular hypertensive and normal eyes. Indian J Ophthalmol. 2007; 55: 283–287.
Inuzuka H, Kawase K, Sawada A, Aoyama Y, Yamamoto T. Macular retinal thickness in glaucoma with superior or inferior visual hemifield defects. J Glaucoma. 2013; 22: 60–64.
Kuang TM, Zhang C, Zangwill LM, Weinreb RN, Medeiros FA. Estimating lead time gained by optical coherence tomography in detecting glaucoma before development of visual field defects. Ophthalmology. 2015; 122: 2002–2009.
Sung KR, Kim JS, Wollstein G, Folio L, Kook MS, Schuman JS. Imaging of the retinal nerve fibre layer with spectral domain optical coherence tomography for glaucoma diagnosis. Br J Ophthalmol. 2011; 95: 909–914.
Mwanza JC, Durbin MK, Budenz DL, et al. Glaucoma diagnostic accuracy of ganglion cell-inner plexiform layer thickness: comparison with nerve fiber layer and optic nerve head. Ophthalmology. 2012; 119: 1151–1158.
Rolle T, Briamonte C, Curto D, Grignolo FM. Ganglion cell complex and retinal nerve fiber layer measured by Fourier-domain optical coherence tomography for early detection of structural damage in patients with preperimetric glaucoma. Clin Ophthalmol. 2011; 5: 961–969.
Nouri-Mahdavi K, Nowroozizadeh S, Nassiri N, et al. Macular ganglion cell/inner plexiform layer measurements by spectral domain optical coherence tomography for detection of early glaucoma and comparison to retinal nerve fiber layer measurements. Am J Ophthalmol. 2013; 156: 1297–1307.e2.
Hood DC, Raza AS, de Moraes CG, Liebmann JM, Ritch R. Glaucomatous damage of the macula. Prog Retin Eye Res. 2013; 32: 1–21.
Sommer A, Duggan C, Auer C, Abbey H. Analytic approaches to the interpretation of automated threshold perimetric data for the diagnosis of early glaucoma. Trans Am Ophthalmol Soc. 1985; 83: 250–267.
Bagga H, Greenfield DS, Knighton RW. Macular symmetry testing for glaucoma detection. J Glaucoma. 2005; 14: 358–363.
Salgarello T, Colotto A, Valente P, et al. Posterior pole retinal thickness in ocular hypertension and glaucoma: early changes detected by hemispheric asymmetries. J Glaucoma. 2005; 14: 375–383.
Asrani S, Rosdahl JA, Allingham RR. Novel software strategy for glaucoma diagnosis: asymmetry analysis of retinal thickness. Arch Ophthalmol. 2011; 129: 1205–1211.
Seo JH, Kim TW, Weinreb RN, Park KH, Kim SH, Kim DM. Detection of localized retinal nerve fiber layer defects with posterior pole asymmetry analysis of spectral domain optical coherence tomography. Invest Ophthalmol Vis Sci. 2012; 53: 4347–4353.
Um TW, Sung KR, Wollstein G, Yun SC, Na JH, Schuman JS. Asymmetry in hemifield macular thickness as an early indicator of glaucomatous change. Invest Ophthalmol Vis Sci. 2012; 53: 1139–1144.
Sullivan-Mee M, Ruegg CC, Pensyl D, Halverson K, Qualls C. Diagnostic precision of retinal nerve fiber layer and macular thickness asymmetry parameters for identifying early primary open-angle glaucoma. Am J Ophthalmol. 2013; 156: 567–577.e1.
Yamashita T, Sakamoto T, Kakiuchi N, Tanaka M, Kii Y, Nakao K. Posterior pole asymmetry analyses of retinal thickness of upper and lower sectors and their association with peak retinal nerve fiber layer thickness in healthy young eyes. Invest Ophthalmol Vis Sci. 2014; 55: 5673–5678.
Yamada H, Hangai M, Nakano N, et al. Asymmetry analysis of macular inner retinal layers for glaucoma diagnosis. Am J Ophthalmol. 2014; 158: 1318–1329.e3.
Kim YK, Yoo BW, Kim HC, Park KH. Automated detection of hemifield difference across horizontal raphe on ganglion cell--inner plexiform layer thickness map. Ophthalmology. 2015; 122: 2252–2260.
Hodapp E PRI, Anderson DR. Hodapp E, Parrish RKII, Anderson DR. Clinical Decisions in Glaucoma. St Louis, MO: The CV Mosby Co. 1993: 52–61.
Nascimento VC, Kasahara N, Cohen R, Almeida GV, Mandia C, Paolera MD. Location and frequency of visual field defects as measured by SITA (Swedish Interactive Threshold Algorithm) strategy in primary open angle glaucoma [in Portuguese]. Arq Bras Oftalmol. 2005; 68: 661–665.
Kotera Y, Hangai M, Hirose F, Mori S, Yoshimura N. Three-dimensional imaging of macular inner structures in glaucoma by using spectral-domain optical coherence tomography. Invest Ophthalmol Vis Sci. 2011; 52: 1412–1421.
Takayama K, Hangai M, Durbin M, et al. A novel method to detect local ganglion cell loss in early glaucoma using spectral-domain optical coherence tomography. Invest Ophthalmol Vis Sci. 2012; 53: 6904–6913.
Tan O, Chopra V, Lu AT, et al. Detection of macular ganglion cell loss in glaucoma by Fourier-domain optical coherence tomography. Ophthalmology. 2009; 116: 2305–2314.e2.
Kim NR, Lee ES, Seong GJ, Kim JH, An HG, Kim CY. Structure-function relationship and diagnostic value of macular ganglion cell complex measurement using Fourier-domain OCT in glaucoma. Invest Ophthalmol Vis Sci. 2010; 51: 4646–4651.
Seong M, Sung KR, Choi EH, et al. Macular and peripapillary retinal nerve fiber layer measurements by spectral domain optical coherence tomography in normal-tension glaucoma. Invest Ophthalmol Vis Sci. 2010; 51: 1446–1452.
Mori S, Hangai M, Sakamoto A, Yoshimura N. Spectral-domain optical coherence tomography measurement of macular volume for diagnosing glaucoma. J Glaucoma. 2010; 19: 528–534.
Hirashima T, Hangai M, Nukada M, et al. Frequency-doubling technology and retinal measurements with spectral-domain optical coherence tomography in preperimetric glaucoma. Graefes Arch Clin Exp Ophthalmol. 2013; 251: 129–137.
Lisboa R, Paranhos A, Weinreb RN, Zangwill LM, Leite MT, Medeiros FA. Comparison of different spectral domain OCT scanning protocols for diagnosing preperimetric glaucoma. Invest Ophthalmol Vis Sci. 2013; 54: 3417–3425.
Miraftabi A, Amini N, Gornbein J, et al. Local variability of macular thickness measurements with sd-oct and influencing factors. Trans Vis Sci Tech. 2016; 5 (4): 5.
Kawaguchi C, Nakatani Y, Ohkubo S, Higashide T, Kawaguchi I, Sugiyama K. Structural and functional assessment by hemispheric asymmetry testing of the macular region in preperimetric glaucoma. Jpn J Ophthalmol. 2014; 58: 197–204.
Hwang YH, Ahn SI, Ko SJ. Diagnostic ability of macular ganglion cell asymmetry for glaucoma. Clin Exp Ophthalmol. 2015; 43: 720–726.
Thonginnetra O, Greenstein VC, Chu D, Liebmann JM, Ritch R, Hood DC. Normal versus high tension glaucoma: a comparison of functional and structural defects. J Glaucoma 2010; 19: 151–157.
Miraftabi A, Amini N, Morales E, et al. Macular SD-OCT outcome measures: comparison of local structure-function relationships and dynamic range. Invest Ophthalmol Vis Sci. 2016; 57: 4815–4823.
Figure 1
 
The 200 × 200 macular cube (Carl Zeiss Meditec, Inc.) consists of 40,000 A-scan measurements that were clustered into a 20 × 20 grid of superpixels centered on the fovea. Local asymmetry measures were defined based on comparison of the GCIPL thickness on the three corresponding pairs of rows above and below the temporal horizontal raphe including four or eight temporal columns, both individually and in combination. The best performing local asymmetry parameter was the one incorporating the three rows above and below the temporal raphe consisting of eight superpixels.
Figure 1
 
The 200 × 200 macular cube (Carl Zeiss Meditec, Inc.) consists of 40,000 A-scan measurements that were clustered into a 20 × 20 grid of superpixels centered on the fovea. Local asymmetry measures were defined based on comparison of the GCIPL thickness on the three corresponding pairs of rows above and below the temporal horizontal raphe including four or eight temporal columns, both individually and in combination. The best performing local asymmetry parameter was the one incorporating the three rows above and below the temporal raphe consisting of eight superpixels.
Figure 2
 
Receiver operating characteristic curves comparing glaucoma diagnostic capabilities of the best local asymmetry index, global asymmetry index, inferotemporal, and minimal GCIPL thicknesses, and combination of the best local asymmetry and regional GCIPL thickness measures.
Figure 2
 
Receiver operating characteristic curves comparing glaucoma diagnostic capabilities of the best local asymmetry index, global asymmetry index, inferotemporal, and minimal GCIPL thicknesses, and combination of the best local asymmetry and regional GCIPL thickness measures.
Table 1
 
Demographic and Clinical Findings of Subjects Enrolled in the Study
Table 1
 
Demographic and Clinical Findings of Subjects Enrolled in the Study
Table 2
 
Comparison of Sectoral GCIPL Thicknesses and Local and Global AIs Between Study Groups
Table 2
 
Comparison of Sectoral GCIPL Thicknesses and Local and Global AIs Between Study Groups
Table 3
 
Values of AUC for GCIPL Thickness and Local and Global AIs
Table 3
 
Values of AUC for GCIPL Thickness and Local and Global AIs
×
×

This PDF is available to Subscribers Only

Sign in or purchase a subscription to access this content. ×

You must be signed into an individual account to use this feature.

×