Associations between macular retinal vasculature and severity of idiopathic epiretinal membrane

Abstract

Background

To demonstrate the associations between the morphology of macular retinal vasculature and disease severity of idiopathic epiretinal membrane (ERM).

Methods

Macular structures were assessed using optical coherence tomography (OCT), and were classified as “with pseudohole” or “without pseudohole”. The 3 × 3 mm macular OCT angiography images were analyzed using the Fiji software to obtain the vessel density, skeleton density, average vessel diameter, vessel tortuosity, fractal dimension, and foveal avascular zone (FAZ)-related parameters. The correlations between these parameters and ERM grading as well as visual acuity were analyzed.

Results

For ERM with or without a pseudohole, increased average vessel diameter, decreased skeleton density, and decreased vessel tortuosity were all associated with inner retinal folding and thickened inner nuclear layer, indicating more severe ERM. In 191 eyes without a pseudohole, the average vessel diameter increased, fractal dimension decreased and vessel tortuosity decreased with increasing ERM severity. The FAZ was not associated with ERM severity. Decreased skeleton density (r = -0.37), vessel tortuosity (r = -0.35), and increased average vessel diameter (r = 0.42) were correlated with worse visual acuity (All P < 0.001). In 58 eyes with pseudoholes, a larger FAZ was associated with a smaller average vessel diameter (r = -0.43, P = 0.015), higher skeleton density (r = 0.49, P < 0.001), and vessel tortuosity (r = 0.32, P = 0.015). However, none of the retinal vasculature parameters correlated with visual acuity and central foveal thickness.

Conclusion

Increased average vessel diameter, decreased skeleton density, decreased fractal dimension and decreased vessel tortuosity were good indicators of ERM severity and associated visual impairment.

Background

Idiopathic epiretinal membrane (ERM) is characterized by fibrocellular proliferation on the inner surface of the retina [1, 2]. It can cause alteration in the retinal structures and may result in loss of visual acuity and metamorphopsia [1, 2]. The advent of optical coherence tomography (OCT) enables the visualization of changes in inner and outer retinal structures. Several specific biomarkers have been identified to depict the features of ERM, including ectopic inner foveal layer [3], disorganization of retinal inner layers [4], cotton ball signs [5, 6], and the gap between the ERM and the retinal surface [7]. However, the OCT images used for analysis were mostly B-scan images, which could only represent one or two cross sections of the whole macula.

In fact, ERM not only causes significant retinal vessel displacement as observed by fundus photography [8, 9], but also leads to changes in macular retinal vasculature, as demonstrated by optical coherence tomography angiography (OCTA). Previous studies have demonstrated that eyes with idiopathic ERM had decreased superficial and deep vessel density at the macula [10, 11] but increased vessel density at the central subfield [12, 13] and a smaller area of the foveal avascular zone (FAZ) [14, 15] compared to healthy control subjects. Greater vascular changes were associated with greater central foveal thickness (CFT) and poor postoperative visual outcomes [10]. Most studies have investigated the OCTA parameters obtained automatically using built-in software including vessel density and FAZ. However, other OCTA parameters such as skeleton density, vessel diameter, vessel tortuosity, and fractal dimension could describe more characteristics of the vascular changes accompanied with the retinal fold and distortion in ERM, and they could only be analyzed after image processing for the en face OCTA images [16, 17]. These parameters have been evaluated in various retinal diseases and found to be associated with disease severity [16,17,18,19,20], but seldom discussed in idiopathic ERM. We thought that the en face OCTA images could reveal the conditions of the central macula which may not be shown in one or two cross-sectional B-scan OCT images.

In the present study, we aimed to describe retinal vasculature changes using OCTA in a large cohort of patients with idiopathic ERM. We analyzed the associations between these vascular parameters and disease severity, as assessed using visual acuity, CFT, and the ERM staging system. Additionally, we compared the features of retinal vasculature changes between patients with and without macular pseudoholes. Furthermore, the OCT biomarkers and clinical features that predicted vascular parameters were analyzed.

Methods

Patients diagnosed with idiopathic ERM by Dr. Yi-Ting Hsieh at the National Taiwan University Hospital between January 2016 and December 2019 were retrospectively enrolled. Patients with retinopathy other than ERM, high myopia, opaque media, visually significant cataract (Lens Opacity Classification System III grades C2–5, P2–5, NC3-6, or NO3-6), and a history of intraocular surgery other than uncomplicated cataract surgery were excluded. Additionally, patients with secondary ERM due to previous surgery, retinal laser therapy, or ocular inflammation were excluded. All patients underwent detailed ophthalmic examinations including best-corrected visual acuity (BCVA), slit-lamp biomicroscopy, dilated fundus examination, OCT, and OCTA. If a patient had bilateral idiopathic ERM, only one eye was randomly selected. Institutional review board approval was obtained from the National Taiwan University Hospital for this study (201803097RINB), and the study adhered to the tenets of the Declaration of Helsinki. The requirement for informed consent was waived owing to the retrospective nature of this study.

OCT measurement

All patients underwent spectral-domain OCT (RTVue RT-100, version 3.5; Optovue, Inc., Fremont, CA, USA) to obtain standard 10 mm vertical and horizontal scans centered at the fovea. The presence of ellipsoid zone disruption, external limiting membrane disruption, inner retinal folding [7], and loss of foveal depression were documented. The CFT was obtained from the retinal thickness map. The thickness of the segmented retinal layers, including the internal limiting membrane-outer plexiform layer (ILM-OPL), ganglion cell layer-inner plexiform layer (GC-IPL), inner nuclear layer (INL), OPL, outer nuclear layer-retinal pigment epithelium (ONL-RPE), and maximal retinal thickness (MRT), was measured using a manual caliber. Eyes were then classified as “with pseudohole” and “without pseudohole,” [21] and macular pseudohole was characterized by a deep foveal pit, verticalized edges, and thickened macula caused by ERM traction [22, 23]. Those without pseudoholes were further classified into three grades of ERM severity according to the staging system of ERM proposed by Govetto et al. [3]. Grade 1 and grade 2 corresponded to the stage 1 and stage 2 ERM, respectively, and grade 3 consisted of the stage 3 and stage 4 ERM. In summary, the grade 1 ERM was characterized by a preserved foveal pit and well-defined retinal layers. The grade 2 ERM was characterized by loss of foveal pit and well-defined retinal layers. The grade 3 ERM was defined as having an ectopic inner foveal layer, with or without disrupted retinal layers.

OCTA Parameters

All patients underwent OCTA using the Optovue RTVue XR Avanti with the AngioVue OCTA system. Only scans with a scan quality indicator of ≥ 5 were eligible for further analysis. A 3 × 3 mm region centered on the fovea was scanned. FAZ parameters, including area, perimeter, and circularity, were measured using a built-in software. Other OCTA parameters were measured using Fiji software (version 2.0.0, National Institutes of Health, Bethesda, MD, USA) to process the en face angiography map of the whole inner retinal slab (ILM-OPL). The image processing procedures were as follows: First, the FAZ was marked and circularity was calculated using the “circularity” function. The original RGB images were converted to 8-bit grayscale images, and the “default” algorithm was used for auto thresholding. The vessel density was calculated as the area occupied by white pixels (blood vessels) divided by the entire area of the binarized image. Next, the “skeletonize” function was used to create a skeletonized image. Skeleton density was calculated as the area occupied by the white pixels (skeletonized blood vessels) divided by the entire area of the skeletonized image. The average vessel diameter was calculated as vessel density divided by skeleton density. The fractal dimension was calculated by the “box-counting” method to quantify the branching complexity [19, 24, 25]. Vessel tortuosity was calculated as an average of the ratio of geodesic distance and the Euclidean distance of each vessel branch [25].

Statistical analysis

All statistical analyses were performed using the R software (version 4.0.3; R Foundation for Statistical Computing, Vienna, Austria). For descriptive statistics, means and standard deviations were calculated for parametric numerical data, and percentages were calculated for categorical variables. BCVA was converted to the logarithm of the minimum angle of resolution (logMAR) for statistical analysis. Analysis of variance was used to compare the clinical characteristics and OCTA parameters between the patients with different ERM grades. Scheffe’s test was used for post-hoc analyses. Pearson’s correlation was used to assess the relationships between vascular parameters and visual acuity, CFT, or FAZ area. Independent t-test and chi-squared test were used to compare the clinical characteristics between patients with and without pseudoholes. Univariate linear regression analysis was performed to analyze the association between OCT characteristics and OCTA biomarkers. P values < 0.05 were considered as statistically significant.

Results

We examined 249 eyes of 249 patients (146 men and 103 women) with idiopathic ERM. Fifty-eight eyes (23.3%) showed macular pseudoholes. The remaining 191 eyes without pseudoholes were further classified into three groups according to ERM severity: 75, grade 1; 50, grade 2; and 66, grade 3.

Pseudohole vs. no pseudohole

Table 1 shows a comparison of clinical characteristics between patients with and without pseudoholes. Patients with pseudoholes were younger and had better BCVA. The retinal thickness and FAZ parameters were not different between the two groups; however, the OCTA parameters were significantly different. Eyes with pseudoholes had a lower skeleton density, larger average vessel diameter, and lower vessel tortuosity. We further analyzed the correlation between the FAZ area and other OCTA parameters in eyes with pseudoholes and found that a larger FAZ area was associated with higher vessel density (P = 0.008, r = 0.34, 95% confidence interval [CI] = 0.09 ‒ 0.55), higher skeleton density (P < 0.001, r = 0.49, 95% CI = 0.26 ‒ 0.67), higher vessel tortuosity (P = 0.015, r = 0.32, 95% CI = 0.07 ‒ 0.53), and smaller average vessel diameter (P = 0.015, r = -0.43, 95% CI = -0.62 ‒ -0.20).

Table 1 Comparisons of clinical characteristics and OCTA parameters between ERM with and without psuedohole

Correlations between OCTA Parameters and ERM Severity, CFT and BCVA

Table 2 shows the characteristics of the eyes without pseudoholes. Patients with more severe ERM grades had increased CFT and MRT. Eyes with grade 3 ERM had worse BCVA than those with grade 1 or 2 ERM. As the ERM severity increased, the skeleton density, fractal dimension, and vessel tortuosity decreased, whereas the average vessel diameter increased. The vessel density increased in grade 2, but then decreased in grade 3 when compared to grade 1 ERM. None of the FAZ parameters was associated with ERM severity. Figure 1 shows some examples of the OCTA and OCT of cases with different ERM severity. Table 3 shows the correlations between OCTA parameters, CFT, and visual acuity. For eyes without pseudoholes, a thicker CFT was associated with lower vessel density, skeleton density, fractal dimension, vessel tortuosity, and a larger average vessel diameter. Regarding BCVA, a larger logMAR, which indicates poorer vision, was associated with lower skeleton density, vessel tortuosity, and higher average vessel diameter. None of the FAZ parameters was associated with CFT or BCVA. For eyes with pseudoholes, none of the OCTA parameters were associated with CFT or BCVA.

Table 2 The clinical characteristics and OCTA parameters in patients with different severity of idiopathic epiretinal membrane

Fig. 1

Optical coherence tomography and en face optical coherence tomography angiography of epiretinal membrane (ERM) with different disease severitiy. A Grade 1 ERM: Preservation of the foveal pit and well-defined retinal layers on the OCT. B Grade 2 ERM: Loss of the foveal pit but preservation of all retinal layers. C Grade 3 ERM: Ectopic inner foveal layer with organized retinal layers. As the grading of ERM increased, the skeleton density appeared to decrease, and the average vessel diameter increased

Table 3 The correlations between OCTA parameters and central foveal thickness and visual acuity

Correlations between OCTA Parameters and OCT Characteristics

Table 4 shows the correlations between the OCT characteristics and OCTA parameters. Among the OCTA parameters, skeleton density, vessel tortuosity, and average vessel diameter were most correlated with OCT characteristics in ERM. Both decreased skeleton density, decreased vessel tortuosity, and increased average vessel diameter were associated with loss of foveal depression, presence of inner retinal folding, and increased inner retinal thickness. Decreased fractal dimension, indicating loss of branching complexity, was associated with increased inner retinal thickness and ELM disruption. Vessel density was negatively associated only with INL thickness. Decreased FAZ circularity was associated with loss of foveal depression and ellipsoid zone disruption. The FAZ area and perimeter were not associated with any of the OCT parameters analyzed.

Table 4 The analysis of correlations between OCT characteristics and OCTA parameters using univariable linear regression

Discussion

Tangential traction of idiopathic ERM can cause significant changes in macular structure and retinal vasculature [10, 11, 13,14,15, 26,27,28]. The superficial capillary plexus (SCP) and deep capillary plexus (DCP) nourish the inner retina and are located in the layer of RNFL and retinal ganglion cell and in the junction between IPL and OPL, respectively [29]. The tractional force of idiopathic ERM from the retinal surface causes significant disorganization of the inner retinal structures and results in changes in these two capillary plexuses [4]. Several previous studies have investigated differences in vascular parameters between eyes with idiopathic ERM and their fellow eyes or healthy controls [10, 11, 13,14,15, 26,27,28]. In the present study, we further analyzed whether these vascular parameters may change with the severity of ERM assessed by the OCT-based staging system proposed by Govetto et al. [3] in which the ectopic inner foveal layer is the critical feature for grading. This ERM staging system exhibited clinical relevance by showing that visual acuity declined and CFT increased in more advanced stages [3]. In the present study, we not only confirmed the correlation between ERM severity and BCVA, but also found that the skeleton density decreased and the average vessel diameter increased proportionally with the severity of ERM. Since the tractional centripetal displacement of the inner retinal layer and Müller cell–driven proliferation simultaneously contribute to the development of the ectopic inner foveal layer in idiopathic ERM [3, 30], the changes in retinal vasculature may be more complicated than the simple result of contraction. Theoretically, the traction force of the ERM may damage retinal capillaries. In addition, as retinal folding and disorganization of retinal layers become more severe, the retinal capillary network becomes more difficult to detect. Furthermore, Kim et al. proposed that diminished microvasculature might be an artifact due to distorted retinal structure or stagnated blood flow caused by the compressive force of the ERM [11]. We hypothesized that the vascular network in eyes with severe ERM would be more difficult to be detected thoroughly, and therefore the vasculature shown in the en face OCTA would also be less complete. As a result, there is a decrease in the skeleton density measured using OCTA. Because large vessels are less affected and easier to detect, the average vessel diameter measured by OCTA will increase as ERM become worse. The findings that skeleton density decreased and average vessel diameter increased in the presence of inner retinal folding, loss of foveal depression, and increased INL thickness and CFT in this study further confirmed our proposal. Furthermore, these changes in the retinal vasculature correlated with worse visual acuity. Therefore, we believe that decreased skeletal density and increased average vessel diameter are good indicators of more severe ERM and functional deterioration.

Previous studies have shown that eyes with idiopathic ERM had lower vessel density in the SCP and DCP than in healthy controls [26], and the changes in vessel density correlated with visual acuity and CFT [10]. The centripetal tractional force led to the displacement of the capillary plexus toward the fovea. Therefore, the vessel density in the foveal region increased, the area of the FAZ decreased, and the vessel density of the parafoveal region decreased [10, 11, 13, 15, 26,27,28, 31]. In addition to the tractional force, pre-existing vascular insufficiency may also attribute to the reduction in blood flow before the development of ERM [26]. However, we found that vessel density was less significantly correlated with ERM severity and CFT. In ERM, large retinal vessels appear thicker in the presence of traction. Compensatory vessel dilatation may also occur in the presence of disturbed perfusion [19]. Kim et al. also proposed that the central displacement of larger retinal vessels and dilatation of stagnated capillaries may both contribute to the increasing vessel diameter and vessel density in more severe ERM [11]. These mixed effects of ERM on changes in vessel density may make it a less important indicator of ERM severity or visual impairment.

Fractal dimension is a measurement of the branching complexity of retinal vessels and usually decreases in myopia, glaucoma, and hypertension; however, it has shown conflicting results in diabetic retinopathy [32]. In the present study, the fractal dimension decreased significantly in grade 3 ERM. Kim et al. reported that the DCP had a decreased fractal dimension and higher lacunarity in more severe ERM stages [11]. They speculated that this finding was an artifact due to distorted retinal structure or stagnated blood flow caused by the compressive force of the ERM [11]. In the present study, decreased fractal dimension was associated with increased INL thickness. INL thickness is a biomarker that reflects the extent of tangential traction from ERM [33]. The abnormal traction force may cause damage to tiny capillaries, resulting in dropout and decreased fractal dimension. Since DCP supplied the INL, the changes in the INL may affect DCP more prominently, which is in line with the findings of Kim et al. [11].

A previous study found that eyes with ERM had more tortuous vessels than fellow eyes [34] and healthy controls [35] and were associated with more significant metamorphopsia [34]. In contrast, the present study demonstrated that the tortuosity of macular retinal vessels decreased in patients with more severe ERM. In addition, decreased vessel tortuosity was correlated with thicker CFT and worse BCVA. Theoretically, tangential traction from the ERM could cause prominent curling of the perifoveal capillary; alternatively, it could straighten and drag the surrounding retinal vessel toward the fovea [8, 36, 37]. The abnormally slow and fluctuated blood flow may cause the vulnerable and pulled capillary wall to become tortuous [16, 38]. However, our study results implied that the straightening effect of the strong tractional force might outweigh the curling effect in more severe ERM. Furthermore, as stated previously, the capillary network may be under detected due to retinal folding and disorganization of retinal layers in severe ERM; therefore, the measured vessel tortuosity could not reflect the presence of torturous capillaries.

Previous studies have shown that eyes with ERM have a smaller FAZ than healthy eyes [10, 11, 13, 15, 26,27,28, 31]. However, the present study found that the FAZ-related parameters had no association with ERM severity, OCT parameters, and visual acuity. Since the FAZ area is highly variable among individuals and could be affected by other factors, including CFT, axial length, sex, and choroidal thickness [28, 31, 39], direct comparison of FAZ among individuals may be biased. Furthermore, FAZ is largely different between fovea center-sparing ERM and fovea-attached ERM. Despite the similar pathogenesis, foveal center-sparing attachment may result in different morphological changes and clinical characteristics from fovea-attached ERM, including pseudoholes. In the present study, we classified patients according to the presence or absence of macular pseudoholes, which appeared as discrete, reddish, round, or oval lesions [23]. We found that eyes with pseudoholes had better visual acuity than those without pseudoholes. Hwang et al. reported that eyes with pseudoholes had similar visual acuity as those with mild ERM, but had better visual acuity than those with significant inner retinal thickening [21]. Since the eyes in our non-pseudohole group had ERM severity ranging from mild to severe, the overall visual acuity in this group may have been worse than that in the pseudohole group. Pierro et al. revealed that the FAZ area was smaller in eyes without pseudoholes than in those with pseudoholes owing to the stronger traction force [40]. In contrast, we did not find a significant difference in FAZ-related parameters between eyes with and without pseudoholes. The conflicting results between the present study and those of Pierro et al. could be explained by the different disease severities and individual variability of the FAZ. Macular pseudoholes with straight edges have the fovea as the epicenter of centripetal traction [22]. A more prominent centripetal force would lead to a shorter distance of the foveal edges and result in a smaller FAZ area. We found that a larger FAZ in eyes with pseudoholes was associated with a higher skeleton density, vessel tortuosity, and smaller vessel diameter. In summary, the FAZ may be a good indicator of ERM severity in comparison to healthy fellow eyes. However, it was less reliable to indicate the severity of ERM among different subjects, as it varied between individuals and with the diverse force direction of ERM. Instead, the skeleton density, vessel tortuosity, and vessel diameter correlated with OCT parameters and visual acuity, and could be regarded as a proxy of fovea-attached ERM severity.

This study has several limitations. First, we did not enroll healthy eyes; therefore, we could not investigate the differences in vascular parameters between healthy eyes and eyes with idiopathic ERM. However, the data provided by previous studies delineated the changes in vascular parameters in eyes with idiopathic ERM. Second, the present study focused on vascular changes in different types and severity of ERM. In the analysis of en face OCTA images, we chose the image from the retinal slab, which included both the SCP and DCP. However, previous studies have suggested that the changes in SCP and DCP may be distinct owing to the different locations and inherent characteristics of capillaries [11]. Third, vessel density is affected by several factors. Although we excluded patients with other types of retinopathy, high myopia, opaque media, systemic disease, and systemic vasoactive medication could have biased our results. Despite of the extensive use of skeletonized metrics from the en face OCTA in the literature, the concerns of low agreement between results obtained using different imaging processing methods and algorithms had been raised recently [41, 42]. Therefore, further studies are necessary to verify whether the findings in this study are consistent with those of other methods. Due to the retrospective nature of the present study, we did not have detailed documentation of metamorphopsia, which may be associated with changes in the macular retinal vasculature. Lastly, the enrollment of severe ERM with disorganized retinal layers was relatively difficult due to the poorer image quality and small case numbers. This may limit the analysis of microvascular changes and their correlations with visual impairment in this subgroup of patients.

In conclusion, the present study investigated the morphology of macular retinal vasculature in idiopathic ERM of different types and severity. Eyes with more severe fovea-attached ERM had decreased skeleton density, a larger average vessel diameter, and decreased vessel tortuosity, which was also associated with worse visual acuity and increased CFT. In eyes with pseudoholes, smaller FAZ were associated with decreased skeleton density, larger vessel diameter, and decreased vessel tortuosity. Both skeleton density and average vessel diameter were significantly correlated with loss of foveal pit depression, inner retinal folding, and a thicker INL, all of which are characteristics of severe ERM. Based on these findings, we propose that skeleton density, vessel tortuosity, and average vessel diameter could be potential indicators of disease severity and strength of the centripetal tractional force in idiopathic ERM. Future studies could focus on the changes in these parameters in post-operative eyes and their role in predicting visual outcomes.