Original Research ARTICLE
Investigating the Fractal Dimension of the Foveal Microvasculature in Relation to the Morphology of the Foveal Avascular Zone and to the Macular Circulation in Patients With Type 2 Diabetes Mellitus
- 1Miller School of Medicine, Bascom Palmer Eye Institute, University of Miami, Miami, FL, United States
- 2Department of Biophysics, University of Belgrade, Belgrade, Serbia
- 3Retinology Unit, Pallas Kliniken, Olten, Switzerland
- 4Department of Ophthalmology, Semmelweis University, Budapest, Hungary
In this study, we examined the relationship between the fractal dimension (FD), the morphology of the foveal avascular zone (FAZ) and the macular circulation in healthy controls and patients with type 2 diabetes mellitus (T2DM) with and with no diabetic retinopathy (DR). Cross-sectional data of 47 subjects were analyzed from a 5-year longitudinal study using a multimodal optical imaging approach. Healthy eyes from nondiabetic volunteers (n = 12) were selected as controls. Eyes from patients with T2DM were selected and divided into two groups: diabetic subjects with mild DR (MDR group, n = 15) and subjects with DM but without DR (DM group, n = 20). Our results demonstrated a higher FD in the healthy group (mean, 1.42 ± 0.03) than in the DM and MDR groups (1.39 ± 0.02 and 1.35 ± 0.03, respectively). Also, a bigger perimeter, area, and roundness of the FAZ were found in MDR eyes. A significant difference in area and perimeter (p ≤ 0.005) was observed for the MDR group supporting the enlargement of the FAZ due to diabetic complications in the eye. A moderate positive correlation (p = 0.014, R2 = 43.8%) between the FD and blood flow rate (BFR) was only found in the healthy control group. The BFR calculations revealed the lowest values in the MDR group (0.98 ± 0.27 μl/s vs. 1.36 ± 0.86 μl/s and 1.36 ± 0.57 μl/sec in the MDR, DM, and healthy groups, respectively, p = 0.2). Our study suggests that the FD of the foveal vessel arborization could provide useful information to identify early morphological changes in the retina of patients with T2DM. Our results also indicate that the enlargement and asymmetry of the FAZ might be related to a lower BFR because of the DR onset and progression. Interestingly, due to the lack of FAZ symmetry observed in the DM and MDR eyes, it appears that the distribution of flow within the retinal vessels loses complexity as the vascular structures distributing the flow are not well described by fractal branching. Further research could determine how our approach may be used to aid the diagnosis of retinal neurodegeneration and vascular impairment at the early stage of DR.
Diabetes mellitus is one of the leading causes of vision loss (Diabetes Fact Sheet, 2017). Central vision loss, which is the most common impairment related to diabetic macular edema (DME), has a big impact on quality of life (Bourne et al., 2013; Korobelnik et al., 2014). The worldwide prevalence of DM is predicted to grow to 430 million patients by 2030, according to the World Health Organization (Korobelnik et al., 2014). Diabetes prevalence has also increased at a faster rate in lower income countries when compared to wealthier nations countries (Diabetes Fact Sheet, 2017). Also, rates of retinopathy are higher among people with type 1 diabetes mellitus (T1DM), individuals with longer duration of diabetes, and Caucasian populations (Yau et al., 2012; Bourne et al., 2013; Korobelnik et al., 2014). There is also a correlation between lower socioeconomic status and higher rates of retinopathy (Yau et al., 2012; Bourne et al., 2013).
Microaneurysms, capillary nonperfusion, and ischemia within the retina are the characteristic pathological features of DR (Elman et al., 2010; Hwang et al., 2015, 2016; Ishibazawa et al., 2015). These pathophysiologic changes can be associated with several complications, such as DME and diabetic macular ischemia (Beltramo and Porta, 2013; American Academy of Ophthalmology., 2014; Varma et al., 2014; Couturier et al., 2015; Das et al., 2015; Bradley et al., 2016). These complications compound impaired blood and oxygen supply of the neuroglial tissues of the retina. The expression of vascular endothelial growth factor (VEGF), which acts as an angiogenic agent and increases vascular permeability, is enhanced in the hypoxic environment (Choi et al., 2013). Diabetic maculopathy is caused by a combination of both VEGF-mediated factors and mediators (Agemy et al., 2015; Mastropasqua et al., 2015; Gorczynska et al., 2016). High glucose levels lead to microvasculopathy with alterations in the blood-retinal barrier, causing pericyte loss and endothelial cell-cell junction breakdown (Kaur et al., 2008). This capillary disruption increases vascular permeability and the pooling of fluid within the plexiform layers of the retina and the subretinal space (Ishibazawa et al., 2015). These pathophysiological changes can result in a gradual loss of visual acuity.
The FAZ is the macular region that is most susceptible to retinal changes in individuals with diabetes (Choi et al., 2013; Varma et al., 2014). The FAZ is the central part of the macula and is surrounded by interconnected capillary beds. This vascular network terminates in the central macula forming a vascular ring-shaped border with an average diameter of 500–600 μm (de Carlo et al., 2015a,c). In some patients with DME, central visual loss may be due to edema as well as ischemia occurring due to capillary dropout sufficient to increase the FAZ area (Choi et al., 2013). Most DR patients will not experience vision changes until late-stage disease, therefore, early detection and immediate intervention may better preserve vision (Aiello, 2003). Consequently, early detection and accurate staging are critical for determining optimal management. DME treatments have been developed, including focal or grid photocoagulation and anti-VEGF therapy, both of which have been demonstrated to be effective (American Academy of Ophthalmology., 2014; Couturier et al., 2015; Wiley et al., 2016). Anti-VEGF therapy has also been shown recently to slow the progression and, in some cases, reverse the degree of ophthalmoscopically observed nonproliferative retinopathy (Ip et al., 2012, 2015). However, the pathology involved in the retinal hemodynamics throughout the course of DR is not completely understood (Pemp and Schmetterer, 2008; Kur et al., 2012; Stitt et al., 2016).
Various diagnostic techniques exist to assess the structure and function of the retina. Ophthalmoscopy and fundus photographs are standard techniques used in DR, with the fundus photographs being comparable to ophthalmoscopy under dilated pupils while remaining cost effective for screening in diabetes clinics (Lee et al., 1993). Fluoroscein angiography is also recognized as a useful tool for healthcare professionals diagnosing and treating DR (Agemy et al., 2015). However, it requires venipuncture which may lead to allergic reactions and, in rare cases, death due to anaphylaxis (Yannuzzi et al., 1986). In addition, the technique is costly and time-consuming, requiring up to 30 min only for image acquisition itself (Matsunaga et al., 2014; Sim et al., 2014; Di et al., 2016). Nevertheless, it has historically been the standard in the assessment of DR and DME. Retinal hemodynamic abnormalities and retinal oxygen metabolism have also been investigated in patients with DR (Kohner et al., 1995; Bursell et al., 1996; Calles-Escandon and Cipolla, 2001; Schram et al., 2005). For example, the retinal function imager (RFI) (Optical Imaging Ltd, Rehovot, Israel) is a noninvasive imaging technique that has been used to investigate the microcirculation in the retina of patients with DM (Grinvald et al., 2004; Nelson et al., 2005).
Research in different areas of complications associated to DM is constantly evolving (Campagnoli et al., 2017; Somfai et al., 2018). The search for risk biomarkers characterizing preclinical abnormalities is fundamental to fast-track the discovery of novel treatments. For example, FD is one of the vascular architectural parameters commonly used to quantify changes in the retinal branching pattern and vascular density due to disease progression (Lim et al., 2009; Cosatto et al., 2010). Fractal geometry studies of the retinal vasculature can be performed by fractal analysis, a mathematical method used to measure complexity in natural phenomenon (Mandelbrot, 1982). The concept of fractal geometry was first described by Mandelbrot in 1989 (Smith et al., 1996; Fernández and Jelinek, 2001; Di Ieva et al., 2015). Later, Family et al. introduced this method in ophthalmology and since then, interest in studying the association between the FD of the retinal vasculature and disease severity and progression has dramatically increased (Family et al., 1989; Fractals medicine., 1991; Cheung et al., 2009; Grauslund et al., 2010; Yau et al., 2010; Aliahmad et al., 2014; Broe et al., 2014). Although, the retinal vasculature tree could be quantified with various methods of fractal analysis (Stosic and Stosic, 2006; Macgillivray et al., 2007), the digital retinal images could be investigated through complexity, space-filling, shape, and tortuosity of retinal blood vessels. These characteristics could be quantified by the box-counting method of fractal analysis (Milošević, 2015).
In this study, we examined the relationship between the fractal dimension (FD), the morphology of the foveal avascular zone (FAZ) and the macular circulation in healthy controls and patients with type 2 diabetes mellitus (T2DM) with and with no diabetic retinopathy (DR). Our results suggest that the FD of the foveal vessel arborization could provide useful information to identify early morphological changes in the retina of patients with T2DM.
The study was approved by the Institutional Review Board (University of Miami, Miami, FL, USA). The research adhered to the tenets outlined in the Declaration of Helsinki and written informed consent was obtained from each study subject. In this prospective study, enrollment was offered to patients with DM referred to the comprehensive ophthalmology clinic that had DR up to early treatment diabetic retinopathy study (ETDRS) level 35 and without macular edema, as well as diabetic patients with no retinopathy and healthy individuals (Group ETDRSR, 1981).
Patients with proliferative disease, clinically significant macular edema (CSME), and anatomic abnormalities that might confound the evaluation of macular architecture (such as glaucoma, vitreoretinal traction, and epiretinal membranes) were excluded. Patients with medical conditions that might affect visual function, those taking medications that might affect retinal thickness (e.g., chloroquine or anti-cholesterol agents containing niacin), recent cataract surgery, previous vitrectomy, or unstable blood sugars were also excluded.
The routine ophthalmic examination was carried out with dilated fundoscopy and patients were divided into two groups based on the absence of DR (DM group) and presence of mild DR (MDR group). Any eyes with more severe DR (i.e., greater than ETDRS level 35) were excluded from the study. Study subjects (age-matched) were selected from a 5-year longitudinal study based on the quality of the overall imaging data required to perform all analyses. Table 1 shows the demographics of the study population. A total of 47 study participants (58 eyes) were identified with good quality images from the RFI system.
This analysis was performed by the RFI system, which is based on a standard fundus camera extended by a customized stroboscopic flash lamp system. A green (“red-free”) light with a spectrum of 548 ± 75 nm is used for illumination and the interval between consecutive flashes is typically 17.5 ms. One session of RFI data consists of 8 images with a resolution of 1,024 × 1,024 pixels in an area of 4.3 × 4.3 mm or 7.2 × 7.2 depending on the choice of field of view (20° or 35°) during the imaging acquisition. In this study, all the images were captured with the setting of 20° field of view (FOV) at a resolution of 4.3 microns/pixel. The heartbeats of the patient were monitored with a finger probe sensor and the image acquisition was synchronized with the cardiac cycle to neutralize the effects of pulsation on arterial blood flow velocity (BFV) (Grinvald et al., 2004; Nelson et al., 2005; Tian et al., 2016).
Once the image was acquired, the RFI built-in software generated (a) the flow movie (a.k.a. “ratio video”) through differential processing so that the motion of individual clusters of red blood cells can be followed by the human eye; and (b) a non-invasive capillary perfusion map (nCPM) was generated through analyzing the difference of pixel intensities in adjacent frames (Nelson et al., 2005; Izhaky et al., 2009). A good quality scanning session is characterized by sharp vessel borders on the raw fundus images, clear red blood cell movement along the vessels on ratio videos and a visible capillary network on the nCPM. Therefore, images were evaluated for optical quality, exposure and focus. We obtained 3 or more good-quality sessions for each eye by the same experienced photographer with at least 5 good images per session that were selected for further analysis using a custom-built software (Tian et al., 2016). In our method, the BFV is calculated by maximizing cross-correlation of intensity profiles between adjacent frames and the BFR is computed by multiplying the BFV with the cross-sectional area (Tian et al., 2016). All BFR measurements were obtained for the overall arteries, overall veins and overall vessels (i.e., arteries and veins) for each study group.
The fractal analysis of the retinal vascular network was performed using the box-counting method (Smith et al., 1996; Fernández and Jelinek, 2001; Milošević, 2015, 2016; Rajkovic et al., 2017). The RFI images were imported in Image J (National Institutes of Health, Bethesda, MD) and used to calculate the FD after grayscale format conversion (Figure 1). The box counting method generates data by “covering” the object with a rectangular coordinate grid and breaking the data into boxes and then analyzing the subsets by counting the number of boxes (Smith et al., 1996; Fernández and Jelinek, 2001). The lower and upper box-dimensions of a subset F⊂Rn are respectively defined by
and if lower and upper values are equal, then the common value is referred to as the box-counting dimension of F and is denoted by
Where Nδ(F) can be the smallest number of cubes of side δ (naturally, in 3D) that covers F or the largest number of disjoint cubes of side δ with centers in F (Falconer, 1989).
When plane-projection of the object is analyzed, this method measures FD by covering the image with a rectangular coordinate grid with cell size r and counts the number of boxes, where the cell size is expressed as the number of foreground pixels. The number of squares N(r) needed to cover the image is given by a power law
Where D is the box dimension, calculated as an absolute value of the slope of the log-log relationship between N(r) and r as previously described in detail (Milošević, 2015).
Morphology Analysis of the FAZ
The foveal morphology analysis was performed on the nCPM composite images, which were generated from the individual scans taken with the RFI (Tian et al., 2016). These capillary maps were then used to gather information regarding the FAZ, which can be identified in the image as the central area where no vasculature is present. An active contouring model was used to identify and outline the FAZ (Nikolay, 2016; Bates et al., 2018). The active contouring program is run by constant user supervision; therefore, each run was reviewed ensuring that it accurately represented the FAZ. If it appeared to be inaccurate, the simulation was rerun with different parameters, specifically with a different region of interest drawn to capture the area better. If it was still not representative of the actual FAZ, the image was excluded from the analysis due to poor quality. This process was done blindly, in that the group of the patient was unknown. Following identification of the FAZ using active contouring, the images were analyzed using ImageJ (National Institutes of Health, Bethesda, MD) (Smith et al., 1996; Schneider et al., 2012). This software allows for the easy acquisition of data regarding the region of interest, including area, circumference, perimeter, and the maximum/minimum Feret diameter. Area, perimeter, and roundness are parameters that we used for this study.
The Human Research Ethics Committee of the University of Miami, Miami, FL, USA approved all protocols and methods described in this study. The research adhered to the tenets outlined in the Declaration of Helsinki. Informed consent was obtained from all participants following a thorough explanation of all test procedures.
Linear regression was used to determine the relationship between the BFR, FD, and FAZ parameters among the three groups. A one-way ANOVA test was used to determine if there was a difference present in any group, and then a post hoc Kruskal-Wallis test was used to identify these individual group differences. In all cases, a p-value of 0.05 was used to define significance (Katz and McSweeney, 1980; Armstrong et al., 2000).
FD was calculated for the three study groups. ANOVA showed a statistically significant difference (p < 0.001) between all study groups. The highest FD values were obtained for the healthy group (1.42 ± 0.03) compared to those calculated for the DM and MDR groups (1.39 ± 0.02 and 1.35 ± 0.03, respectively) (see Figure 2). An example of the nCPM and corresponding skeleton images used in the fractal analysis is shown in Figure 3.
Figure 2. Box plots displaying the FD results. The middle 50% (mean and 95% CI) of the data groups are as follows: DM (1.3765–1.4105), MDR (1.3255–1.365), and Healthy (1.4015–1.447).
Figure 3. Sample images used in the fractal analysis. Images in the left column are the composite nCPM images generated from the RFI system, while those in the right are their respective skeleton images that were used in the fractal analysis. Row (A) is from the healthy group, Row (B) from the DM group, and Row (C) from the MDR group. We note that the blood vessel types used in the analysis were the major arterial and venous trees.
The intergroup differences measured with the Kruskal-Wallis analysis revealed a statistically significantly larger FD in the healthy group compared to the one obtained from the DM and MDR groups (p = 0.002). Also, a statistically significant difference (p < 0.001) was found between the healthy and the MDR group, as well as between the DM and MDR groups (p < 0.001).
The area, perimeter and roundness of the FAZ were calculated to characterize the morphology of the FAZ region (Bates et al., 2018) A summary of the results for the three study groups is presented in Table 2. ANOVA showed significance in all parameters with multiple significant differences in the pairwise comparisons as seen in Table 2.
Table 2. Results of the FAZ characterization for all three groups (first three columns show the morphological parameters).
Blood Flow Rate
Figure 4 shows the BFR results obtained for the overall blood vessels (i.e., both arteries and veins) per study group. Our results showed 0.98 ± 0.27 μl/s, 1.36 ± 0.86 μl/s, and 1.36 ± 0.57 μl/s in the MDR, DM, and healthy groups, respectively. No significant difference (p = 0.2) was found in the ANOVA analysis.
Figure 4. Box Plot displaying the BFR results obtained for the overall blood vessels in the three study groups. The middle 50% of the data groups are as follows: DM (0.89–1.42), MDR (0.809–1.17), and Healthy (0.85–1.868).
Relationship Between the Vascular FD, the Morphology of the FAZ, and Macular Circulation
Linear regression results showed significant differences for the healthy group when analyzing the correlations between the vascular FD and BFR parameters (see Table 3). For all three groups of patients, FD was calculated and compared with BFR data for overall arteries, veins, and vessels. A summary of the results for all three study groups is presented in Table 3. Only the healthy control group showed a significant linear correlation between the FD and BFR parameters.
Table 3. Linear regression results obtained for all three groups after analyzing the correlations between the FD and BFR. Significant results (p < 0.05) are highlighted in bold.
In comparing the BFR with the FAZ characteristics, only the DM group showed a statistically significant linear relationship. These relationships are shown in Figure 5, where the R2 value and the equation of the line of best fit is displayed on their respective graphs.
Figure 5. (A–E) Linear regression results obtained for the DM group of patients, which is the only group that displayed significance when comparing the BFR to the FAZ characteristics. Each figure is labeled with measured parameters. The line represents thel ine of best fit, of which the equation is displayed in the graph region. In all cases, the “y-axis” represents BFR in μl/s and the “x-axis” represents FAZ area in mm2. (A) BFR vs. FAZ area in veins (B) BFR vs. FAZ area in all vessels (C–E) BFR vs. FAZ Perimeter for arteries, veins, and all vessels, respectively. The R2 is the coefficient of determination, which compares the ratio of the average difference of the data point and the line of the best fit with that of the data point and the average “y” value.
In this study, we investigated the FD of the foveal vessel arborization in relation to the morphology of the FAZ and the macular circulation in patients with T2DM. Our results demonstrated the highest FD in the healthy group compared to the DM and MDR groups, with a significant difference between the three groups. Our results are similar to the outcomes presented by Grauslund et al., where the FD of 94 patients with T1DM without proliferative retinopathy were compared with 79 T1DM patients with proliferative retinopathy and found that the group with the most severe disease had lower FD (Grauslund et al., 2010). Similarly, Aliahmad et al. found that the healthy group had higher FD compared with the diabetic group and suggested that a low FD could be a result of the retina being less effective, which may lead to an increased risk of complications like proliferative retinopathy (Mainster, 1990; Aliahmad et al., 2014).
Numerous studies have assessed the vascular complexity (FD) in patients with different types of DM and stages of DR with contradicting findings. A 16-years study that monitored 180 patients with T1DM found that FD of the retinal vasculature generally decreased in these population. It was also found that lower FD could predict neuropathy (Broe et al., 2014). On the other hand, Cheung et al. and Yau et al. found that an increase in vascular FD was associated with an increased incidence of retinopathy (Cheung et al., 2009; Yau et al., 2010). However, all these studies had different study designs and used fundus machines with different resolutions and illumination settings. The quality of the images acquired from the study subjects was also dissimilar. Also, the region of interest for the FD calculation was not the same for all studies. For example, our study analyzed the FD in a 20° FOV while the other studies looked at the vessels at a larger scale. This difference in methodologies may cause contradicting findings among studies. Huang et al. found that the FD must be calculated under very rigorous settings after investigating the reliability of the vascular FD calculated from retinal images acquired with 5 different fundus cameras (Huang et al., 2016). Therefore, when comparing different studies, it is of great importance to consider the study design, the settings used for the analysis of the acquired images and image resolution among other key factors.
The results of our study showed that there were significant differences in FAZ parameters that describe different morphology characteristics between all 3 groups. Results of the FAZ characterization showed no significant differences in the DM group in comparison to the healthy group, while the MDR group had significantly higher FAZ area and perimeter in comparison to the healthy and DM groups. Also, a significant difference in roundness was only observed between DM and MDR groups. These findings might indicate that the FAZ roundness might be a good indicator of DR onset and progression. Intriguingly, as reported in our previous study, the fact that the FAZ area was larger in the MDR group but more asymmetric in the DM group suggest a possible anisotropy in the mechanical properties of the diabetic retina with no retinopathy. This anisotropy may trigger the FAZ elongation in a preferred direction as probably a result of autoregulation (Bates et al., 2018).
The findings in the current study support the loss of symmetry in the FAZ expansion with worsening of the retinopathy condition. Our results are in accordance with the fact that the enlargement of the FAZ area is an indicator of DR onset and progression besides being an indicator of visual prognosis in patients with DME (Bates et al., 2018). It is well known that larger FAZ area is associated with progression of DR, compared to healthy individuals, as reported after using FA (Sakata et al., 2006). Additionally, Optical Coherence Tomography-Angiography (OCTA) studies have found a negative correlation between FAZ area and visual acuity. However, FAZ size is not fixed among individuals; therefore, its normal variation makes assessment of retinal pathology in terms of FAZ size very difficult (de Carlo et al., 2015b; Freiberg et al., 2015; Mammo et al., 2015; Spaide et al., 2015; Takase et al., 2015; Al-Sheikh et al., 2016; Bhanushali et al., 2016; Samara et al., 2016; Tan et al., 2016).
There are many discrepancies in studies targeting retinal blood flow, mainly due to differences in the device used and in the study design. Particularly, the variability in these studies is affected by factors including the recruitment of patients with T1DM diabetes vs. T2DM, controlled vs. poorly controlled diabetes, age, presence or absence of other co-morbidities such as systemic hypertension (Kohner, 1975; Grunwald et al., 1993, 1996; Schmetterer and Wolzt, 1999; Pournaras et al., 2008). Our study population consisted of age-matched groups composed by healthy controls and T2DM subjects with and with no mild DR. The diabetes disease condition in T2DM patients was in control but most of the T2DM patients were hypertensive (under control). The fastest BFR for overall vessels (that is, arteries and veins taken together) was found in healthy subjects, compared with the DM and MDR groups, whereas the slowest BFR was observed in the MDR group. This is in correlation with previous investigations where reduced retinal blood flow and vessel stiffness seemed to be associated with nonperfusion in the vasculature network (Kohner, 1975; Grunwald et al., 1996; Schmetterer and Wolzt, 1999; Sakata et al., 2006; Pournaras et al., 2008).
In comparing the correlation of the BFR and FD for the three groups, a significant correlation between FD and BFR was obtained for the healthy group in all three groupings of the vessels (overall arteries, overall veins and overall vessels). Nevertheless, we did not find any significant correlation for the DM and MDR groups. Also, the BFR and vascular FD among healthy subjects and diabetic patients with and with no DR have not been compared in previous studies. According to our findings, it seems that the flow distribution of the retinal vessels loses complexity while fractal branching does not adequately describe the vascular structures involved in this process.
Although the BFR was lower in the MDR group than in the DM and healthy groups, intriguingly, the correlation between BFR and FAZ was statistically significant with low R-squared values only for the DM group. It is possible that our data contain an intrinsically higher amount of unexplainable variability. Also, there was a negative correlation between FAZ perimeter and BFR for the DM group when compared in overall arteries, overall veins and overall vessels. As well as there was a negative correlation between FAZ area and BFR for overall veins and overall vessels. Interestingly, we saw no statistically significant correlation between BFR and FAZ roundness in the DM group. Also, there was no statistically significant correlation between BFR and FAZ in the MDR and healthy groups.
The reduction in BFR in patients with DM may be due to morphological changes of the vascular bed in combination with lack of capability of vascular autoregulation and the decrease of blood fluidity (Ashton, 1974; Sinclair et al., 1982; McMillan, 1989). Our results could be explained by the fact that an increase in FAZ area and perimeter can result in decreased BFV, and consequently a decreased BFR, as demonstrated by the negative correlation observed between BFR and FAZ size. However, our study does not have the full power to prove this outcome as there are many factors that have an influence on microcirculation in diabetic patients with DR, such as duration of disease, changes in ocular biomechanics, and the presence of other concomitant systemic diseases (e.g., controlled hypertension in our diabetic groups) that can contribute to a slower BFR in patients with DR. A statistically significant negative correlation between capillary BFV and FAZ size in T2DM was found in a previous study that used FA as the measurement method (Sakata et al., 2006). However, a correlation between FAZ parameters and BFR in patients with DR or a healthy control group seems to be unavailable in the current literature. Our fractal dimension results also revealed the potential use of this method to quantify the progressive change in DR between the increased and decreased vessel complexity stages. The complexity index of the retinal vascular pattern characterized by fractal analysis may uncover potential regulation of specific markers of disease status.
There are important limitations to the present study that need to be addressed. There was a different sample size for each one of the 3 groups due to the exclusion of eyes because of the poor image quality obtained with the RFI unit in some eyes with media opacities. Also, only subjects with T2DM were included. Consequently, the degree to which our outcomes can be generalized to individuals with Type 1 diabetes is uncertain. Especially, retinal structure and function may be affected by factors such as hyperlipidemia, older age, and hypertension in Type 2 diabetes. Furthermore, we conducted our study in a relatively small sample of patients. A better understanding of the correlation between BFR, FD, and FAZ measurements in a bigger number of subjects and longitudinal studies is devised as a future study.
There are many structural and hemodynamic parameters that can play a role in the development of DR. Particularly, the development of advanced imaging will facilitate that these quantitative measurements can help with the identification of changes in the retinal structure affected by various diseases with greater precision and detail (Tian et al., 2016). Our results suggest that the FD of the foveal vessel arborization in conjunction with other functional and structural parameters could provide useful information to identify early morphological changes in the retinal tissue of patients with T2DM. The data also lead us to believe that the enlargement and asymmetry of the FAZ area might be related to a lower BFR associated with the onset and progression of DR (Krawitz et al., 2017). Notably, due to the lack of FAZ symmetry observed in DM and MDR eyes, it appears that the distribution of flow within the retinal vessels loses complexity as the vascular structures distributing the flow are not well described by fractal branching. In addition, despite the availability of many different studies about FAZ and DM or FD and FAZ, to our knowledge, there are no known studies that correlate FD with hemodynamic and structural parameters. However, further longitudinal research is warranted to determine how our approach may be used to aid diagnosis of retinal neurodegeneration and vascular impairment at the early stage of DR. There is no doubt that there is a need for a better understanding of structural and hemodynamic parameters which are integrally interdependent. Specifically, the multimodal measurements in our future work would not only provide details of retinal pathophysiology but could possibly contribute as a biomarker in disease staging.
DC conceived and designed the study. DC, JT, GS, WS, W-HL, AK, NG, SP, and MK performed the study; NB, NM, JT, GS, WF, JS, DC, and MK analyzed the data. JT, DC, WF, JS, SP, and NM: contributed reagents, materials, and analysis tools. MK, NB, NM, JT, DC, WS, and GS contributed to the writing of the manuscript.
This study was supported in part by a NIH Grant No. NIH R01EY020607, a NIH Center Grant No. P30-EY014801, by an unrestricted grant to the University of Miami from Research to Prevent Blindness, Inc., and by an Eotvos Scholarship of the Hungarian Scholarship Fund.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Agemy, S. A., Scripsema, N. K., Chirag, M. S., Chui, T., Garcia, P. M., Lee, J. G., et al. (2015). Retinal vascular perfusion density mapping using optical coherence tomography angiography in normals and diabetic retinopathy patients. Retina 35, 2253–2263. doi: 10.1097/IAE.0000000000000862
Aliahmad, B., Kumar, D. K., Sarossy, M. G., and Jain, R. (2014). Relationship between diabetes and grayscale fractal dimensions of retinal vasculature in the Indian population. BMC Ophthalmol. 14:152 doi: 10.1186/1471-2415-14-152
Al-Sheikh, M., Akil, H., Pfau, M., and Sadda, S. R. (2016). Swept-source OCT angiography imaging of the foveal avascular zone and macular capillary network density in diabetic retinopathy. Invest. Ophthalmol. Vis. Sci. 57:3907–3913. doi: 10.1167/iovs.16-19570
American Academy of Ophthalmology (2014). American Academy of Ophthalmology Retina and Vitreous Panel. Preferred Practice Pattern® guidelines. Diabetic retinopathy. San Francisco, CA: American Academy of Ophthalmology. Available online at: http://www.aao.org/ppp
Armstrong, R. A., Slade, S. V., and Esperjesi, F. (2000). An introduction to analysis of variance (ANOVA) with special reference to data from clinical experiments in optometry. Ophthalmic Physiol. Opt. 20, 235–241. doi: 10.1016/S0275-5408(99)00064-2
Bates, N. M., Tian, J., Smiddy, W., Lee, W. H., Somfai, G., Feuer, W., et al. (2018). Relationship between the morphology of the foveal avascular zone, retinal structure, and macular circulation in patients with diabetes mellitus. Sci. Rep. 8:5355. doi: 10.1038/s41598-018-23604-y
Bhanushali, D., Anegondi, N., Gadde, S. G., Srinivasan, P., Chidambara, L., Yadav, N. K., et al. (2016). Linking retinal microvasculature features with severity of diabetic retinopathy using optical coherence tomography angiography. Invest. Ophthalmol. Vis. Sci. 57, 519–525. doi: 10.1167/iovs.15-18901
Bourne, R. R., Stevens, G. A., White, R. A., Smith, J. L., Flaxman, S. R., Price, H., et al. (2013). Causes of vision loss worldwide, 1990–2010: a systematic analysis. Lancet Global Health 1:e339–49. doi: 10.1016/S2214-109X(13)70113-X
Bradley, P. D., Sim, D. A., Keane, P. A., Cardoso, J., Agrawal, R., Tufail, A., et al. (2016). The evaluation of diabetic macular ischemia using optical coherence tomography angiography. Invest. Ophthalmol. Vis. Sci. 57, 626–631. doi: 10.1167/iovs.15-18034
Broe, R., Rasmussen, M. L., Frydkjaer-Olsen, U., Olsen, B. S., Mortensen, H. B., Peto, T., et al. (2014). Retinal vascular fractals predict long-term microvascular complications in type 1 diabetes mellitus: the Danish Cohort of Pediatric Diabetes (DCPD1987). Diabetologia 57, 2215–2221. doi: 10.1007/s00125-014-3317-6
Bursell, S. E., Clermont, A. C., Kinsley, B. T., Simonson, D. C., Aiello, L. M., and Wolpert, H. A. (1996). Retinal blood flow changes in patients with insulin-dependent diabetes mellitus and no diabetic retinopathy. Invest. Ophthalmol. Vis. Sci. 37, 886–897.
Campagnoli, T. R., Somfai, G. M., Tian, J., DeBuc, D. C., and Smiddy, W. E. (2017). Noninvasive, high-resolution functional macular imaging in subjects with retinal vein occlusion. Ophthalmic Surg. Lasers Imag. Retina 48, 799–809. doi: 10.3928/23258160-20170928-04
Cheung, N., Donaghue, K. C., Liew, G., Rogers, S. L., Wang, J. J., Lim, S. W., et al. (2009). Quantitative assessment of early diabetic retinopathy using fractal analysis. Diabetes Care 32, 106–110. doi: 10.2337/dc08-1233
Choi, W., Mohler, K. J., Potsaid, B., Lu, C. D., Liu, J. J., Jayaraman, V., et al. (2013). Choriocapillaris and choroidal microvasculature imaging with ultrahigh speed OCT angiography. PLoS ONE 8:e0081499. doi: 10.1371/journal.pone.0081499
Cosatto, V. F., Liew, G., Rochtchina, E., Wainwright, A., Zhang, Y., Hsu, W., et al. (2010). Retinal vascular fractal dimension measurement and its influence from imaging variation: results of two segmentation methods. Curr. Eye Res. 35, 850–856. doi: 10.3109/02713683.2010.490628
Couturier, A., Mané, V., Bonnin, S., Erginay, A., Massin, P., Gaudric, A., et al. (2015). Capillary plexus anomalies in diabetic retinopathy on optical coherence tomography angiography. Retina 35, 2384–2391. doi: 10.1097/IAE.0000000000000859
de Carlo, T. E., Bonini Filho, M. A., Chin, A. T., Adhi, M., Ferrara, D., Baumal, C. R., et al. (2015a). Spectral-domain optical coherence tomography angiography of choroidal neovascularization. Ophthalmol 122, 1228–1238. doi: 10.1016/j.ophtha.2015.01.029
de Carlo, T. E., Chin, A. T., Bonini Filho, M. A., Adhi, M., Branchini, L., Salz, D. A., et al. (2015b). Detection of microvascular changes in eyes of patients with diabetes but not clinical diabetic retinopathy using optical coherence tomography angiography. Retina 35, 2364–2370. doi: 10.1097/IAE.0000000000000882
Di Ieva, A., Esteban, F. J., Grizzi, F., Klonowski, W., and Martin-Landrove, M. (2015). Fractals in the neurosciences, part II: clinical applications and future perspectives. Neuroscientist 21, 30–43. doi: 10.1177/1073858413513928
Di, G., Weihong, Y., Xiao, Z., Zhikun, Y., Xuan, Z., Yi, Q., et al. (2016). A morphological study of the foveal avascular zone in patients with diabetes mellitus using optical coherence tomography angiography. Graefes Arch. Clin. Exp. Ophthalmol. 254, 873–879. doi: 10.1007/s00417-015-3143-7
Diabetes Fact Sheet (2017). World Health Organization. Available online at: http://www.who.int/mediacentre/factsheets/fs312/en/
Elman, M. J., Aiello, L. P., Beck, R. W., Bressler, N. M., Bressler, S. B., Edwards, A. R., et al. (2010). Randomized trial evaluating ranibizumab plus prompt or deferred laser or triamcinolone plus prompt laser for diabetic macular edema. Ophthalmology 117, 1064–1077. doi: 10.1016/j.ophtha.2010.02.031
Freiberg, F. J., Pfau, M., Wons, J., Wirth, M. A., Becker, M. D., and Michels, S. (2015). Optical coherence tomography angiography of the foveal avascular zone in diabetic retinopathy. Graefes Arch. Clin. Exp. Ophthalmol. 254, 1051–1058. doi: 10.1007/s00417-015-3148-2
Gorczynska, I., Migacz, J. V., Zawadzki, R. J., Capps, A. G., and Werner, J. S. (2016). Comparison of amplitude-decorrelation, speckle-variance and phase-variance OCT angiography methods for imaging the human retina and choroid. Biomed. Opt. Express 7, 911–942. doi: 10.1364/BOE.7.000911
Grauslund, J., Green, A., Kawasaki, R., Hodgson, L., Sjølie, A. K., and Wong, T. Y. (2010). Retinal vascular fractals and microvascular and macrovascular complications in type 1 diabetes. Ophthalmol 117, 1400–1405. doi: 10.1016/j.ophtha.2009.10.047
Grinvald, A., Bonhoeffer, T., Vanzetta, I., Pollack, A., Aloni, E., Ofri, R., et al. (2004). High-resolution functional optical imaging: from the neocortex to the eye. Opthalmol. Clin. North Am. 17, 53–67. doi: 10.1016/j.ohc.2003.12.003
Group ETDRSR (1981). Diabetic retinopathy study. report number 6. design, methods, and baseline results. report number 7. a modification of the airlie house classification of diabetic retinopathy. Invest. Ophthalmol. Vis. Sci. 21, 1–226.
Grunwald, J. E., Brucker, A. J., Grunwald, S. E., and Riva, C. E. (1993). Retinal hemodynamics in proliferative diabetic retinopathy. a laser Doppler velocimetry study. Invest. Ophthalmol. Vis. Sci. 34, 66–71.
Huang, F., Dashtbozorg, B., Zhang, J., Bekkers, E., Abbasi-Sureshjani, S., Berendenschot, T. T., et al. (2016). Reliability of using retinal vascular fractal dimension as a biomarker in the diabetic retinopathy detection. J. Ophthalmol. 2016:6259047. doi: 10.1155/2016/6259047
Hwang, T. S., Gao, S. S., Liu, L., Lauer, A. K., Bailey, S. T., Flaxel, C. J., et al. (2016). Automated quantification of capillary nonperfusion using optical coherence tomography angiography in diabetic retinopathy. JAMA Ophthalmol. 134, 367–373. doi: 10.1001/jamaophthalmol.2015.5658
Hwang, T. S., Jia, Y., Gao, S. S., Bailey, S. T., Lauer, A. K., Flaxel, C. J., et al. (2015). Optical coherence tomography angiography features of diabetic retinopathy. Retina 35, 2371–2376. doi: 10.1097/IAE.0000000000000716
Ip, M. S., Domalpally, A., Hopkins, J. J., Wong, P., and Ehrlich, J. (2012). Long-term effects of ranibizumad on diabetic retinopathy severity and progression. Arch. Ophthalmol. 130, 1145–1152. doi: 10.1001/archophthalmol.2012.1043
Ip, M. S., Domalpally, A., Sun, J. K., and Ehrlich, J. S. (2015). Long-term effects of therapy with ranibizumab on diabetic retinopathy severity and baseline risk factors for worsening retinopathy. Ophthalmol 122, 367–374. doi: 10.1016/j.ophtha.2014.08.048
Ishibazawa, A., Nagaoka, T., Takahashi, A., Omae, T., Tani, T., Sogawa, K., et al. (2015). Optical coherence tomography angiography in diabetic retinopathy: a prospective pilot study. Am. J. Ophthalmol. 160, 35–44. doi: 10.1016/j.ajo.2015.04.021
Izhaky, D., Nelson, D. A., Burgansky-Eliash, Z., and Grinvald, A. (2009). Functional imaging using the retinal function imager: direct imaging of blood velocity, achieving fluorescein angiography-like images without any contrast agent, qualitative oximetry, and functional metabolic signs. Jpn. J. Ophthalmol. 53, 345–351. doi: 10.1007/s10384-009-0689-0
Kaur, C., Foulds, W. S., and Ling, E. A. (2008). Blood-retinal barrier in hypoxic ischaemic conditions: Basic concepts, clinical features and management. Prog. Retin. Eye Res. 27, 622–647. doi: 10.1016/j.preteyeres.2008.09.003
Korobelnik, J. F., Do, D. V., Schmidt-Erfurth, U., Boyer, D. S., Holz, F. G., Heier, J. S., et al. (2014). Intravitreal aflibercept for diabetic macular edema. Ophthalmol 121, 2247–2254. doi: 10.1016/j.ophtha.2014.05.006
Krawitz, B. D., Mo, S., Geyman, L. S., Agemy, S. A., Scripsema, N. K., Garcia, P. M., et al. (2017). Acircularity index and axis ratio of the foveal avascular zone in diabetic eyes and healthy controls measured by optical coherence tomography angiography. Vision Res. 139, 177–186. doi: 10.1016/j.visres.2016.09.019
Kur, J., Newman, E. A., and Chan-Ling, T. (2012). Cellular and physiological mechanisms underlying blood flow regulation in the retina and choroid in health and disease. Prog. Retin. Eye Res. 31, 377–406. doi: 10.1016/j.preteyeres.2012.04.004
Lee, V. S., Kingsley, R. M., Lee, E. T., Lu, M., Russell, D., Asal, N., et al. (1993). The diagnosis of diabetic retinopathy. ophthalmoscopy versus fundus photography. Ophthalmol 100, 1504–1512. doi: 10.1016/S0161-6420(93)31449-1
Lim, S. W., Cheung, N., Wang, J. J., Donaghue, K. C., Liew, G., Amirul Islam, D. M., et al. (2009). Retinal vascular fractal dimension and risk of early diabetic retinopathy: a prospective study of children and adolescents with type 1 diabetes. Diabetes Care 32, 2081–2083. doi: 10.2337/dc09-0719
Macgillivray, T. J., Patton, N., Doubal, F. N., Graham, C., and Wardlaw, J. M. (2007). Fractal analysis of the retinal vascular network in fundus images. Conf. Proc. IEEE Eng. Med. Biol. Soc. 1, 6455–6458. doi: 10.1109/IEMBS.2007.4353837
Mammo, Z., Balaratnasingam, C., Yu, P., Xu, J., Heisler, M., Mackenzie, P., et al. (2015). Quantitative noninvasive angiography of the fovea centralis using speckle variance optical coherence tomography. Invest. Ophthalmol. Vis. Sci. 56, 5074–5086. doi: 10.1167/iovs.15-16773
Mastropasqua, R., Di Antonio, L., Di Staso, S., Agnifili, L., Di Gregorio, A., Ciancaglini, M., et al. (2015). Optical coherence tomography angiography in retinal vascular diseases and choroidal neovascularization. J. Ophthalmol. 2015:343515. doi: 10.1155/2015/343515
Milošević N. T. (2015). “Fractal analysis of two dimensional images: parameters of the space-filling and shape,” in Proceedings of 20th International Conference on Control Systems and Computer Science Vol. 2 (Bucharest), 539–544.
Milošević N. T. (2016). “The morphology of the brain neurons: box-counting method in quantitative analysis of 2D image,” in The Fractal Geometry of the Brain, ed Di A. Ieva (New York, NY: Springer), 109–126.
Nelson, D. A., Krupsky, S., Pollack, A., Aloni, E., Belkin, M., Vanzetta, I., et al. (2005). Special report: noninvasive multi-parametric functional optical imaging of the eye. Ophthalmic Surg. Lasers Imaging 36, 57–66.
Nikolay, S. (2016). Active Contours Implementation and Test Platform GUI. Mathworks File Exchange. Available online at: https://www.mathworks.com/matlabcentral/fileexchange/30284-active-contours-implementation—test-platform-gui?
Pournaras, C. J., Rungger-Brändle, E., Riva, C. E., Hardarson, S. H., and Stefansson, E. (2008). Regulation of retinal blood flow in health and disease. Prog. Retin. Eye Res. 27, 284–330. doi: 10.1016/j.preteyeres.2008.02.002
Rajkovic, N., Krstonošic, B., and Milošević, N. T (2017). Box-counting method of 2D neuronal image: method modification and quantitiative analysis demonstrated on images from the monkey and human brain. Comput. Math. Methods Med. 2017:8967902. doi: 10.1155/2017/8967902
Sakata, K., Funatsu, H., Harino, S., Noma, H., and Hori, S. (2006). Relationship between macular microcirculation and progression of diabetic macular edema. Ophthalmol 113, 1385–1391. doi: 10.1016/j.ophtha.2006.04.023
Samara, W. A., Shahlaee, A., Adam, M. K., Khan, M. A., Chiang, A., Maguire, J. I., et al. (2016). Quantification of diabetic macular ischemia using optical coherence tomography angiography and its relationship with visual acuity. Ophthalmology 124, 235–244. doi: 10.1016/j.ophtha.2016.10.008
Schram, M. T., Chaturverdi, N., Schalwijk, C. G., Fuller, J. H., and Stehouwer, C. D. (2005). Markers of inflammation are cross-sectionally associated with microvaacular complications and cardiovascular disease in type 1 diabetes – the EURODIAB prospective complications study. Diabetologia 48, 370–378. doi: 10.1007/s00125-004-1628-8
Sim, D. A., Keane, P. A., Fung, S., Karampelas, M., Sadda, S. R., Fruttiger, M., et al. (2014). Quantitative analysis of diabetic macular ischemia using optical coherence tomography. Invest. Ophthalmol. Vis. Sci. 55, 417–423. doi: 10.1167/iovs.13-12677
Sinclair, S. H., Grunwald, J. E., Riva, C. E., Braunstein, S. N., Nichols, C. W., and Schwarz, S. S. (1982). Retinal vascular autoregulation in diabetes mellitus. Ophthalmology 89, 748–750. doi: 10.1016/S0161-6420(82)34720-X
Smith, T. G. Jr., Lange, G. D., and Marks, W. B. (1996). Fractal methods and results in cellular morphology – dimensions, lacunarity, and multifractals. J. Neurosci. Methods 69, 123–136. doi: 10.1016/S0165-0270(96)00080-5
Somfai, G. M., Gerding, H., and DeBuc, D. C. (2018). The use of optical coherence tomography for detection of early diabetic retinopathy. Klin Monatsbl Augenheilkd 234, 1–8. doi: 10.1055/s-0044-101827
Stitt, A. W., Curtis, T. M., Chen, M., Medina, R. J., McKay, G. J., Jenkins, A., et al. (2016). The progress in understanding and treatment of diabetic retinopathy. Prog. Retin Eye Res. 51, 156–186. doi: 10.1016/j.preteyeres.2015.08.001
Takase, N., Nozaki, M., Kato, A., Ozeki, H., Yoshida, M., and Ogura, Y. (2015). Enlargement of foveal avascular zone in diabetic eyes evaluated by en face optical coherence tomography angiography. Retina 35, 2377–2383. doi: 10.1097/IAE.0000000000000849
Tan, C. S., Lim, L. W., Chow, V. S., Chay, I. W., Tan, S., and Cheong, K. X. (2016). Optical coherence tomography angiography evaluation of the parafoveal vasculature and its relationship with ocular factors. Invest. Ophthalmol. Vis. Sci. 57, 224–234. doi: 10.1167/iovs.15-18869
Tian, J., Somfai, G. M., Campagnoli, T. R., Smiddy, W. E., and DeBuc, D. C. (2016). Interactive retinal blood flow analysis of the macular region. Microvasc. Res. 104, 1–10. doi: 10.1016/j.mvr.2015.11.003
Varma, R., Bressler, N. M., Doan, Q. V., Gleeson, M., Danese, M., Bower, J. K., et al. (2014). Prevalence of and risk factors for diabetic macular edema in the United States. JAMA Ophthalmol. 132, 1334–1340. doi: 10.1001/jamaophthalmol.2014.2854
Wiley, H. E., Thompson, D. J., Bailey, C., Chew, E. Y., Cukras, C. A., Jaffe, G. J., et al. (2016). A crossover design for comparative efficacy: a 36-week randomized trial of bevacizumab and ranibizumab for diabetic macular edema. Ophthalmology 123, 841–849. doi: 10.1016/j.ophtha.2015.11.021
Yannuzzi, L. A., Rohrer, K. T., Tindel, L. J., Sobel, R. S., Costanza, M. A., Shields, W., et al. (1986). Fluorescein angiography complication survey. Ophthalmology 93, 611–617. doi: 10.1016/S0161-6420(86)33697-2
Yau, J. W., Kawasaki, R., Islam, F. M., Shaw, J., Zimmet, P., Wang, J. J., et al. (2010). Retinal fractal dimension is increased in persons with diabetes but not impaired glucose metabolism: the Australian Diabetes, Obesity and Lifestyle (AusDiab) study. Diabetologia 53, 2042–2045. doi: 10.1007/s00125-010-1811-z
Keywords: diabetic retinopathy, fractal analysis, foveal avascular zone, blood flow rate, fractal dimension
Citation: Kostic M, Bates NM, Milosevic NT, Tian J, Smiddy WE, Lee W-H, Somfai GM, Feuer WJ, Shiffman JC, Kuriyan AE, Gregori NZ, Pineda S and Cabrera DeBuc D (2018) Investigating the Fractal Dimension of the Foveal Microvasculature in Relation to the Morphology of the Foveal Avascular Zone and to the Macular Circulation in Patients With Type 2 Diabetes Mellitus. Front. Physiol. 9:1233. doi: 10.3389/fphys.2018.01233
Received: 30 June 2018; Accepted: 15 August 2018;
Published: 05 September 2018.
Edited by:Sladjana Z. Spasić, Institute for Multidisciplinary Research, University of Belgrade, Serbia
Reviewed by:Daniel Goldman, University of Western Ontario, Canada
Josephine Louise Vaughan, University of Newcastle, Australia
Copyright © 2018 Kostic, Bates, Milosevic, Tian, Smiddy, Lee, Somfai, Feuer, Shiffman, Kuriyan, Gregori, Pineda and Cabrera DeBuc. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Delia Cabrera DeBuc, email@example.com
†These authors have contributed equally to this work.