Abstract
Background:
The retinal microvascular density changes have been identified in thyroid-associated ophthalmopathy (TAO) patients. Whereas a lack of research has been done on the diagnostic ability of optical coherence tomography (OCT) combined with optical coherence tomography angiography (OCTA) parameters.
Objectives:
This study aims to evaluate the retina perfusion variations in eyes with active and stable TAO and its diagnostic abilities using OCT and OCTA.
Design:
This is cohort longitudinal retrospective study.
Methods:
A total of 51 patients with TAO and 39 healthy controls (HCs) were recruited. The TAO eyes were divided into active and stable stage groups. The foveal avascular zone (FAZ), macular perfusion density (mPD), and peripapillary PD were measured by OCTA. The peripapillary retinal nerve fiber layer (RNFL), central retinal thickness (CRT), and whole macular volume (wMV) were measured by OCT. Visual evoked potential (VEP) and visual field (VF) were also assessed.
Results:
The mPD of the superficial retinal capillary plexus (SRCP) was significantly different in all subfields among active, stable, and HC groups (p < 0.05) except for the temporal inner (p = 0.137), and the active group achieved the lowest PD. The FAZ size increased significantly in the active and stable groups compared with the HC group (p < 0.001). Significant difference was observed in mPD of deep retinal capillary plexus (DRCP) in all quadrants among three groups (p < 0.05). Moreover, PD parameters of optic nerve head (ONH) and radial peripapillary capillary plexus (RPCP) showed a different trend among three groups (p < 0.05). The r-value of visual field-mean deviation (VF-MD) of TAO with DRCP-whole PD (wPD) and RPCP-wPD was 0.421 and 0.299, respectively (p < 0.05). The DRCP-wPD in OCTA and RNFL in OCT were significantly higher in area under the receiver operating characteristic curve (AUC) than that of HC eyes.
Conclusion:
OCT and OCTA can noninvasively detect the peripapillary and macular changes in various stages of TAO patients, and it might be a high diagnostic value tool to monitor the TAO progression.
Keywords
Introduction
Thyroid-associated ophthalmopathy (TAO) is an autoimmune mediated inflammatory disorder and affects orbital tissue. The condition shows in patients with hyperthyroidism, euthyroidism, or hypothyroidism. 1 It shows a range of clinical manifestations, including upper eyelid retraction, eyeball proptosis, cornea exposure, limited eye movement, and apical compression of the optic nerve. 2 The proliferation of orbital fibroblasts, inflammatory cellular infiltration with plasma cells, lymphocytes, mast cells and macrophages, and the overproduction of glycosaminoglycans and collagen were the pathogenesis features of this disease. 3 All these contribute to the expansion of extraocular muscles and intraorbital fat which lead to the increased volume of orbital contents and retrobulbar pressure.3,4 It is the leading theory of the severe consequence of compressive optic neuropathy (ON) which requires prompt evaluation and treatment. As early optic nerve variations in clinic are inapparent, it is difficult to detect optic nerve involvement in TAO. Retinal and optic disk structure and function, however, may be acquired alterations before the appearance of the ON;5–7 for earlier detection of ON, previous studies has been reported that latency of P100 in the visual evoked potential (VEP) was an early indicator of ON to prevent visual loss.8,9 Moreover, visual field (VF) test, contrast sensitivity test, and color sensation test have been used to evaluate the visual function and prevent irreversible visual-threatening.9–12
Various techniques have been used to evaluate the orbital blood flow of TAO patients. Alimgil et al. 13 revealed that in Graves ophthalmopathy (GO) patients, elevated intraorbital pressure led to the increase in venous pressure and choroidal vessel resistance, and then, ocular blood flow was significantly decreased. As an advanced technique, optical coherence tomography (OCT) and optical coherence tomography angiography (OCTA) enable visualize detailed structure and microvasculature with high-quality retina image. Several studies have evaluated variations of vessel density in TAO patients,14–16 but obtained conflicting results. In addition, microvasculature changes in the retina and the foveal avascular zone (FAZ) size are considered an indicator of the visual acuity, which facilitates further study of early abnormalities in retina of TAO patients. 17 Meanwhile, perfusion density (PD) focuses on differentiating the movable and stable structures, which can be used to detect changes in retinal blood vessels at different clinical stages and obtain more timely and accurate retinal circulation flow maps. Besides, the determination of the treatment of diverse clinic stages of TAO may take into account the combining of OCTA and OCT parameters.17,18 In this study, we aim to evaluate the retina perfusion variations in eyes with active and stable TAO and its diagnostic abilities using OCT and OCTA and its correlation with the visual function.
Materials and methods
Enrollment of participants
All subjects were recruited from the Department of Ophthalmology at Xiangya Hospital, Central South University from December 2019 to June 2021. A total of 51 random eyes from TAO patients were divided into two groups based on the clinical activity scores (CASs): 19 (1) spontaneous orbital pain, (2) staring induced orbital pain, (3) eyelid swelling, (4) eyelid erythema, (5) conjunctival redness, (6) suppurative necrosis, (7) sarcoma inflammation, or plica. Patients with CAS ⩾3/7 were defined as active TAO, and those with CAS ⩽2/7 were stable TAO. All patients with TAO were newly diagnosed. At the same time, 39 healthy volunteers of matched age were recruited for physical examination as the control group. According to the design of randomized controlled trial, taking deep retinal capillary plexus–whole perfusion density (DRCP-wPD) as the evaluation index, referring to previous research and presurvey data, the DRCP-wPD of the control group was 0.299, the DRCP-wPD of the disease group was 0.421, and the common standard deviation of the two groups was 0.18. Under the condition of α = 0.05 and β = 0.2, assuming the homogeneity of variance between the two groups, PASS software was used to calculate 36 cases in each group, which could achieve 80% statistical efficiency. The cases in the control group and the disease group were 39 and 51, respectively, which met the requirements of sample size. Exclusion criteria for patients and healthy participants included (1) thyroid treatment within 3 months, like radioactive iodine therapy, immunosuppressor agents, and thyroid surgery; (2) hormonotherapy within 6 months; (3) ophthalmic surgery history; (4) concomitant eye diseases, like glaucoma, retinal vein occlusion, maculopathy, dysthyroid optic neuropathy (DON), and ON; refractive error [> ±3 diopters (D) sphere and ±2 D cylinder]; (5) complicated with infection or severe systemic diseases; systolic blood pressure greater than 150 mmHg or less than 100 mmHg, diastolic blood pressure greater than 90 mmHg or less than 60 mmHg; and (6) decreased visual acuity [the best-corrected visual acuity (BCVA) in logMAR visual chart, >0.2], apparent VF defect [mean deviation (MD) in Humphrey perimetry, <−10 dB], and apical congestion is evident on orbital CT and magnetic resonance imaging. OCT and OCTA images with the following characteristics were excluded: (1) signal strength index <7; (2) motor artifacts visible on the en face angiogram; (3) local weak signal; and (4) images off-center on fovea or disk.
Ophthalmic and systemic examination
All participants received a complete ophthalmic examination including BCVA, intraocular pressure (IOP) measurement with the Goldmann applanation tonometry, slit-lamp examination, and axial length (AL). The proptosis was determined by the same examiner using the Hertel exophthalmometry. B-scan ultrasonography was used to evaluate the ocular and orbital structure. Perimetry was conducted with the standard 24-2 Swedish Interactive Thresholding Algorithm on the Humphrey Visual Field Analyzer (Carl Zeiss Meditec, Inc., Dublin, CA, USA), 20 and visual field-mean deviation (VF-MD) was recorded. Pattern VEP (GT-2008 V-VI; GOTEC, Chongqing, China) was only performed for all TAO subjects.
OCTA and OCT imaging
Retina microvasculature imaging based on the Nidek RS-3000 Advance device (Nidek, Gamagori, Japan) 21 of all participants was obtained. The PD (total area of perfusion vessels per unit area in the measurement area) of macular and peripapillary data at different retinal layers was recorded. Spectral domain optical coherence tomography (SD-OCT) and OCTA images were collected and analyzed by the updated AngioScan software. The optic nerve head (ONH) and fovea were artificially entered, and each eye was scanned in 6 mm × 6 mm mode. The PD of macular and peripapillary parameters was measured using the inbuilt software. Retina was automatically segmented into superficial retinal capillary plexus (SRCP) and deep retinal capillary plexus (DRCP), including center area, inner and outer rings of diameters 6 mm, and the inner and outer rings segmented into four sectors, that is, superior, nasal, inferior, and temporal. Automatic segmentation defines the panel (determined by automatic segmentation) where the SRCP extends 12 μm below the core layer from the inner limiting membranes. The panel of the DRCP extends from 8 μm below the inner nuclear layer to 12 μm below the outer nuclear layer. The retinal thickness and macular perfusion density (mPD) were obtained. Central retinal thickness (CRT) was represented as the average thickness of the 1 mm ring in the center of the retina, whole macular volume (wMV), and peripapillary retinal nerve fiber layer (RNFL) were measured by an automatic analysis algorithm of OCT. For optic nerve OCTA, the PD of radial peripapillary capillary plexus (RPCP) network (located in the RNFL) and ONH layer was imaged, including the superior area and inferior area of diameters 6 mm. The FAZ area was calculated using the manufacturer’s angiometric software.
Statistics analysis
SPSS 25.0 was used for statistical analysis. Continuous variables with a normal distribution were represented by mean and standard deviation. The independent samples t-test was used to make comparisons between two groups and one-way analysis of variance (ANOVA) (if the variances were equal) and the Brown–Forsythe test (if the variances were not equal) for multiple groups. The median (interquartile range) was used to analyze continuous variables that did not follow a normal distribution. The Wilcoxon rank sum test was recruited to make comparisons between two groups and the Kruskal–Wallis H test for multiple groups. Categorical variables were represented as frequencies (percentage), and the chi-square test was employed to make comparisons between multiple groups. Correlation analysis was performed using the Pearson and Spearman correlations. A difference of p < 0.05 was considered statistically significant. The area under the receiver operating characteristic curve (AUC) was conducted to evaluate the diagnostic capability of different parameters of OCT and OCTA in TAO.
Results
A total of 90 participants were recruited here, including 39 healthy controls (HCs), 28 active patients, and 23 stable patients. No statistical difference was identified in sex distribution (p = 0.279), age (p = 0.126), IOP (p = 0.538), BCVA (p = 0.242), CRT (p = 0.160), and mWV (p = 0.078) among the active, stable, and HC groups (Table 1). VF-MD and P100 latency were obviously different between the active and stable groups (Table 1, p < 0.001). Moreover, proptosis and RNFL showed significant differences among three groups (Table 1, p < 0.001), the values of RNFL were lower in stable group than in HC, while it reached the highest in active group.
Clinical characteristics of all included eyes.
BCVA, best-corrected visual acuity; CAS, clinical activity scores; CRT, central retinal thickness; HC, health control; IOP, intraocular pressure; PROP, proptosis; RNFL, retinal nerve fiber layer; VF-MD, visual field-mean deviation; wMV, whole macular volume.
Chi-square test.
One-way ANOVA.
Kruskal–Wallis H test.
Brown–Forsythe test.
Independent samples t-test.
MPD parameters
The results demonstrated that mPD of the SRCP was significantly different in all subfields except temporal inner (p = 0.137) among three groups (p < 0.05) (Table 2), the active group obtained the lowest PD. The FAZ size increased significantly in the active and stable groups in comparison with the HC group (p < 0.001) (Table 2). Similarly, significant difference was found mPD of DRCP in all quadrants among three groups (p < 0.05) (Table 3). Post hoc test revealed that the difference of mPD between the active and stable groups was statistically significant only in DRCP except fovea and temporal inner (p < 0.05), active and stable groups significant differs from HC in almost all quadrants of DRCP (Table 3), but only several areas in SRCP (Table 2).
Superficial retinal capillary plexus (SRCP) perfusion densities in the active, stable, and HC eyes.
FAZ, foveal avascular zone; HC, healthy control; II, inferior inner; IO, inferior outer; NI, nasal inner; NO, nasal outer; SI, superior inner; SO, superior outer; SRCP-wPD, superficial retinal capillary plexus-whole perfusion density; TI, temporal inner; TO, temporal outer.
One-way ANOVA.
Kruskal–Wallis H test.
Brown–Forsythe test.
Deep retinal capillary plexus (DRCP) perfusion densities in the active, stable, and HC eyes.
DRCP-wPD, deep retinal capillary plexus-whole perfusion density; HC, healthy control; II, inferior inner; IO, inferior outer; NI, nasal inner; NO, nasal outer; SI, superior inner; SO, superior outer; TI, temporal inner; TO, temporal outer.
Brown–Forsythe test.
Peripapillary PD parameters
The PD parameters of ONH and RPCP in all grids were significantly different among three groups (p < 0.05) (Table 4). The average PD values of the ONH in inferior, superior, and whole of the active group were lower than those of the stable group without statistically significant difference (p > 0.05) (Table 4). PD of the RPCP in inferior and superior has the same tendency with ONH, but the difference between active and stable groups was significant (p < 0.001) (Table 4).
Optic nerve head (ONH) and radial peripapillary capillary layer (RPCP) densities in the active, stable, and HC eyes.
I, inferior; ONH, optic nerve head; RPCP, radial peripapillary capillary plexus; S, superior.
One-way ANOVA.
Brown–Forsythe test.
Correlation between vessel density and other parameters
The Pearson correlation coefficients and their corresponding p-value are presented in Table 5. The r-values of VF-MD of TAO patients with DRCP-wPD and RPCP-wPD were 0.421 and 0.299, respectively (p < 0.05) (Table 5). CRT was negatively correlated with VF-MD and positively correlated with P100 latency (p < 0.05) (Table 5).
Pearson correlation coefficient matrix on parameters of optical coherence tomography (OCT) and optical coherence tomography angiography (OCTA) and visual functions.
BCVA, best-corrected visual acuity; CRT, central retinal thickness; DRCP-wPD, deep retinal capillary plexus-whole perfusion density; FAZ, foveal avascular zone; ONH-wPD, optic nerve head-whole perfusion density; RNFL, peripapillary retinal fiber layer; RPCP-wPD, radial peripapillary capillary plexus-whole perfusion density; SRCP-wPD, superficial retinal capillary plexus-whole perfusion density; VF-MD, visual field-mean deviation.
Diagnostic abilities of OCTA and OCT
The AUC for distinguishing active and HC eyes in OCTA was the highest for DRCP-wPD (0.998), followed by RPCP-wPD (0.939), FAZ (0.823), SRCP-wPD (0.755), and ONH-wPD (0.708). In OCT, RNFL was the highest for active/HC (0.731), followed by ONH-wPD (0.708), wMV (0.640), and CRT (0.577). In the pairwise comparison, the DRCP-wPD in OCTA and RNFL in OCT owned a higher AUC than that of differentiating active from HC eyes (Figure 1).

The area under the receiver operator characteristic curves (AUROCs) for differentiating active from healthy controls. (a) Optical coherence tomography (OCT) measurements: CRT thickness (0.577), RNFL thickness (0.731), and wMV (0.640). (b) OCT angiography (OCTA) measurements: SRCP-wPD (0.755), DRCP-wPD (0.998), ONH-wPD (0.708), RPCP-wPD (0.939), and FAZ (0.823).
Discussion
A lack of research has been done on the diagnostic ability of OCT combined with OCTA parameters and whether they can serve as early signs of visual loss at different stages of TAO. The retinal microvascular density changes has been identified in TAO patients in previous studies, 15 but the results were inconsistent. In this study, we excluded the patients with elevated systemic pressure, which may cause increased cardiac output and orbital blood flow. 22 It might partially explain the discrepancy. In addition, the differences in techniques and disease staging, such as inactive or stable stage, were also factors contributing to this difference. 17 Finally, the TAO participants involved in this study were uncomplicated with DON and without any appearance change. We explored the retinal hemodynamic and structure change in TAO patients at active and stable stages according to CAS. Our research studied PD which symbolized the total area of the perfused vasculature per unit area within an OCTA measurement area. In addition, multidepth assessment of retinal perfusion was performed on cross-sectional images. OCTA offers qualitative and quantitative data of the microcirculation and the PD.
We found in OCT and OCTA parameters in Table 2, except for the PD of the temporal inner quadrant in SRCP, the other quadrants were significantly different among three groups. We observed that PD decreased in active TAO patients. Moreover, as revealed from the post hoc analysis p-value, active and stable groups significantly differ from HC in almost all quadrants of DRCP, but only several areas in SRCP. The trend of ONH and RPCP was similar, and the difference of PD in all fields was significant among groups. But no significant difference was observed in ONH (superior) and ONH (inferior) between stable and HC groups. Retinal perfusion depends mainly upon the orbital blood supply. The Doppler ultrasonography have been used to study the hemodynamic variation and concluded that superior ophthalmic vein had lower flow velocity and stasis in TAO patients. 23 It is known that inflammatory changes recurrent in orbital spaces could lead to chronic atrophic. Then, the orbit compression and macular thinning were found in stable patients. 24 Meanwhile, decreased consumption of nutrients and oxygen occurred in atrophic tissue causes lowered blood flow. 23 Compressive ON might be the severe consequence of orbit compression. Besides, orbital fibroblast mediator, stimulated by orbital fibroblasts and turned into adipocytes and fibroblasts, is thought to be an important cause of orbit compression.25,26 The compression of orbit puts pressure on the venous system, causing retro blood flow. In addition, the inflammatory changes can increase orbital arterial blood flow velocities. 27 Meanwhile, inflammatory variations could cause intraorbital congestion and an increased episcleral venous pressure, and then the central retinal vein is obstructed. 28 It is a possible explanation for reduced macular and ONH perfusion in patients compared with HC. The temporal retinal ring is usually the thinnest, making the difference between HC and affected eyes smaller and therefore hard to detect. Exophthalmos is a result of the expansion of orbital fat and muscle tissue into the bone structure. In this study, the Hertel exophthalmometry values were higher in the active group, but not statistically significant. Marked inflammation was observed during the active phase, followed by a period of stability. Recurrent inflammatory in inactive patients leads to chronic atrophic changes and macular thinning.
DRCP is a more sensitive and early indicator in our investigation, which could be caused by various factors, like the different circulation systems. Retinal circulation produced in endothelial cells is regulated by local factors and only supplies the SRCP. 29 Meanwhile, both the DRCP and macular area are feeded by choroidal circulations with less autoregulation, it is known that perfusion pressure and sympathetic innervation are the main factors influencing choroidal blood flow. 17 Chiara Del Noce et al. 30 found that vascular flow of choriocapillaris was reduced in TAO patients compared with HC, which implied that the decreased PD in DRCP of macular caused by the high congestion of ophthalmic veins reduced choroid drainage. Hence, it is speculated that orbital blood determined the fundus perfusion, and the pathological changes in orbital tissues lead to the perfusion variations in the orbit.
Here, we further found FAZ area enlarged in active TAO patients, and the difference was significant when compared with the HC group. FAZ is the most sensitive in retina and the capillary free zone in central macula. Previous research demonstrated that FAZ enlargement promoted the appearance of capillary perfusion more objectively.21,31
The thickening of the peripapillary RNFL in active stage is predicted by an elevated retrobulbar pressure-based etiology and an inflammatory process. McKeag et al. 32 revealed that disk edema appeared in 56% eyes with diagnosed ON, and it was an effective marker for ON. Intro-orbital inflammation and edema may give rise to optic disk swelling in TAO patients.
A normal-looking appearance of optic disk is, however, thought to be the coexistence of optic disk atrophy with axon loss and optic disk swelling and inflammation combined with inflammatory mechanisms and compressive mechanisms. The compression mechanism is found in glaucomatous eyes which finally results in macular thinning and RNFL loss. 33 This is consistent with our OCT results. Therefore, early changes in TAO can be better monitored and understand with the objective parameters of OCT and OCTA.
Furthermore, we also found that DRCP-wPD and RPCP-wPD positively correlated with VF-MD. According to the previous studies, VF-MD defects and VEP were considered as the early sign of DON even in the presence of normal visual acuity,12,34 suggesting that they can be applied to the early assessment of VF changes.
The AUC showed that the DRCP-wPD was of a high value to distinguish the active eyes from HC in OCTA. Combined with the fact that there was no significant difference among three groups in structural parameters (CRT and wMV) in this research, the diagnostic value of OCTA and OCT was confirmed. We inferred the peripapillary and mPD changes might precede structural changes.
Some limitations, however, remained in our research. First, the sample size was small; we can only see the tendency of the parameters reflecting the stage of the disease. In addition, the use of medicine was not considered in TAO patients, which is a biased evaluation of the severity of the disease. Besides, this study was cross-sectional and failed to capture PD changes over time. Finally, further longitudinal and larger sample study should be performed in the following experiment.
Conclusion
In conclusion, the results suggested the peripapillary and mPD changes might precede structural changes. Furthermore, DRCP-wPD and RNFL have a high value in distinguishing the active eyes from HC which confirmed the diagnose value of OCTA and OCT. The PD in DRCP is promising early sign to detect the variations in various stages of TAO more sensitively.
