Miguel O. Bernabeu
Yang Lu
Omar Abu-Qamar
Lloyd P. Aiello2,3
Jennifer K. Sun- 1Centre for Medical Informatics, Usher Institute, The University of Edinburgh, Edinburgh, United Kingdom
- 2Beetham Eye Institute, Joslin Diabetes Center, Boston, MA, United States
- 3Department of Ophthalmology, Harvard Medical School, Boston, MA, United States
Diabetic retinopathy (DR) is a leading cause of vision loss worldwide. Microaneurysms (MAs), which are abnormal outpouchings of the retinal vessels, are early and hallmark lesions of DR. The presence and severity of MAs are utilized to determine overall DR severity. In addition, MAs can directly contribute to retinal neural pathology by leaking fluid into the surrounding retina, causing abnormal central retinal thickening and thereby frequently leading to vision loss. Vascular perfusion parameters such as shear rate (SR) or wall shear stress (WSS) have been linked to blood clotting and endothelial cell dysfunction, respectively in non-retinal vasculature. However, despite the importance of MAs as a key aspect of diabetic retinal pathology, much remains unknown as to how structural characteristics of individual MAs are associated with these perfusion attributes. MA structural information obtained on high resolution adaptive optics scanning laser ophthalmoscopy (AOSLO) was utilized to estimate perfusion parameters through Computational Fluid Dynamics (CFD) analysis of the AOSLO images. The HemeLB flow solver was used to simulate steady-state and time-dependent fluid flow using both commodity hospital-based and high performance computing resources, depending on the degree of detail required in the simulations. Our results indicate that WSS is lowest in MA regions furthest away from the feeding vessels. Furthermore, areas of low SR are associated with clot location in saccular MAs. These findings suggest that morphology and CFD estimation of perfusion parameters may be useful tools for determining the likelihood of clot presence in individual diabetic MAs.
Introduction
As the worldwide prevalence of diabetes mellitus continues to increase, diabetic retinopathy (DR) remains the most common vascular complication in diabetic patients (Kempen et al., 2004; Klein, 2007; Ko et al., 2012). The chronic hyperglycemic state of diabetes results in pathological alterations of retinal microvascular structures and blood flow (Curtis et al., 2009). Retinal microaneurysms (MAs), which are outpouchings of the retinal capillary walls, are one of the earliest clinical signs in the diabetic eye and are among the key lesions for DR severity classification (ETDRS_No10, 1991; ETDRS_No12, 1991; Wilkinson et al., 2003; Hirai et al., 2007). Whereas some MAs do not appear to affect vision, other MAs can be associated with abnormal vascular leakage caused by the local loss of endothelial barrier function. In some cases, this may lead to subsequent retinal edema and associated vision loss (Nunes et al., 2009; Murakami et al., 2011). MA leakage affecting the local neural retina can often be detected by fluorescein angiography (FA), and treated by intraocular injections of anti vascular endothelial growth factor agents or steroids, as well as macular laser photocoagulation (Duh et al., 2017).
Several studies have evaluated the pathogenesis and natural history of MAs using ex vivo (e.g., transmission electron microscopy and scanning electron microscopy) and in vivo (scanning laser ophthalmoscopy and optical coherence tomography) imaging approaches to characterize pericyte loss, basement membrane thickening, and endothelial proliferation and disruption (Wise, 1957; Cogan et al., 1961; de Oliveira, 1966; Ashton, 1974; Moore et al., 1999). One study (Ezra et al., 2013) proposed using MA-to-vessel radius ratio as a potential marker for assessing risk of leakage, and suggested that shear stress at the MA wall may lead to endothelial dysfunction.
Advances in adaptive optics scanning laser microscopy (AOSLO) have recently enabled non-invasive investigation of the living human retina with single cell level resolution (~2 μm) (Tam et al., 2010; Chui et al., 2012, 2013), allowing detailed characterization of MA features (wall hyper-reflectivity, wall deformability), morphology (saccular, fusiform, focal bulge, irregular) and perfusion status (fully/partially perfused or non-perfused). One recent study (Dubow et al., 2014), which combined high resolution AOSLO with FA to provide a high-resolution and high-contrast view of individual MAs, extended the qualitative morphologic classification into six morphology groups.
Retinal MAs are known to be highly dynamic lesions. Over the course of the disease, some lesions will disappear (possibly due to thrombus formation and revascularization) while others will either stabilize or grow. A series of studies (Goatman et al., 2003; Bernardes et al., 2009; Ribeiro et al., 2013) have characterized MA turnover (defined as the sum of the MA formation and disappearance rates Ribeiro et al., 2013) and found this metric to be a predictor of macular edema progression. However, these studies were limited in their ability to fully characterize MAs and did not include perfusion status or morphological characteristics of individual MAs in their analysis.
In a recent study, we demonstrated the feasibility of computational fluid dynamics (CFD) analysis to characterize the hemodynamic environment of the diabetic eye (Lu et al., 2016). Comparable approaches have been extensively used for the characterization of larger scale vascular lesions, such as intracranial aneurysms (IA) (Dhar et al., 2008; Chien et al., 2011). Morphological parameters, such as aneurysm aspect ratio and non-sphericity index (Chien and Sayre, 2014) have been identified as risk factors for rupture of IA. Perfusion parameters, such as velocity, wall shear stress (Tarbell, 2010), and shear rate have been proposed to study IA progression and resolution. In particular, a relationship between shear rate and IA thrombosis has been established (Ribeiro de Sousa et al., 2016), leading to a better understanding of IA progression.
In this study, morphological and CFD analyses of individual diabetic MAs were performed based on high resolution AOSLO technology. Our aim is to develop a method capable of establishing which MA characteristics are associated with a higher risk of leakage or clotting. We propose two novel morphological indices to quantify MA shape and aspect ratio. In addition, we introduce two CFD-based perfusion parameters to predict areas with higher risks of endothelial dysfunction and blood clotting. Finally, we demonstrate how to account for the pulsatile nature of blood flow in the models development and investigate the previous indices throughout the cardiac cycle.
Methods
Imaging Instrument
The AOSLO used in this study was a modified version of the Indiana system described previously (Burns et al., 2007). A near infrared superluminesent diode (SLD) with a central wavelength of 830 nm (BLM-S-830, Superlum, Ireland) was used for imaging. Another SLD with a central wavelength of 780 nm (BLM-S-780, Superlum, Ireland) was used for wavefront sensing. A micro-electro-mechanical system deformable mirror (DM, Multi-DM, Boston Micromachines Corp., Cambridge, MA, USA) provided wavefront correction. The DM has an active area of 4.95 × 4.95 mm and 12 × 12 actuators with a maximum stroke of 5.5 μm. The system uses doubler mirrors to amplify the usable stroke of the DM (Webb et al., 2004). The maximum beam size at the exit pupil is 6.5 mm. Based on theoretical calculations, this AOSLO system is capable of compensating for over 90% of the optical aberrations from an eye with clear media and a dilated pupil, achieving ~2.5 μm resolution. With such resolution, MA structural and perfusion information can be characterized in much greater detail than previously achievable with standard techniques such as fundus photography or fluorescein angiography (Figure 1).
Figure 1. Three MAs imaged by AOSLO superimposed on a digital fundus photograph from an eye with diabetic retinopathy.
Image Processing and Morphological Analysis
MA Segmentation and Skeletonisation
The body and feeding/draining vessels of the MAs under study were manually segmented from AOSLO images. The MA outline was created by using the Fiji/ImageJ “Polygon Selections” tool to define series of line segments along the MA wall. The outline was adjusted based on both the scattered light images (Figure 2a) and their corresponding “perfusion map” (standard deviation map) images calculated from the AOSLO frames (Figure 2b). The “Create Mask” function was used to turn the segmentation file into a binarized figure file. In the segmentation process, the length of each feeding/draining capillary was taken to be roughly equal to the MA body length along the flow direction axis (see section Hemodynamic Analysis for more details). The region representing the MA body was differentiated from the feeding/draining vessels (Figure 2c). In the subset of MAs in which we could identify blood clots, the perfused versus clotted areas within the MA body were also segmented (Figure 2d). Direction of flow was recorded from the AOSLO videos of each MA. Binary masks defining the two-dimensional projection of the MA body along with feeding/draining capillaries were prepared for further processing. In the cases where clots were present, the clotted area was also included in the binary mask. We employed the methodology described previously (Bernabeu et al., 2014) to calculate the MA centerline and radii along the centerline from the Voronoi diagram of the pixels defining the boundary of each binary mask (Attali and Montanvert, 1997). Briefly, the centerline is the subset of the Voronoi diagram defining the medial axis of the mask. For any point along the centerline, its radius is given by the largest circle centered on that point and inscribed within the mask (Figure 3).
Figure 2. Segmentation of a partially clotted MA, (a) delineation of the MA body and feeding/draining capillaries against AOSLO movie frame, (b) same delineation against perfusion map, (c) MA body (red) and MA feeding/draining capillaries (black) were independently segmented, (d) using AOSLO video files for reference, the perfused (blue) vs. clotted regions (green) within the MA body were also segmented.
Figure 3. BNR is defined as the largest caliber registered along the centerline of the MA body (blue arrows) divided by the narrowest caliber along the feeding/draining vessels (red arrows).
Morphological Analysis
In this work, we propose two novel indices to describe the morphology of a retinal MA: the body-to-neck ratio (BNR) and the asymmetry ratio (AR). BNR provides a measure of how dilated the MA body is in relation to the caliber of the feeding/draining capillaries. BNR is defined as the quotient between the MA body width and the caliber of the feeding/draining vessels (see Figure 3). Chien et al. employed a similar measure to characterize arterial brain aneurysms and found a trend for increases in this index when comparing aneurysms before and after rupture (Chien and Sayre, 2014). BNR is computed based on the skeleton/radii analysis described in section MA Segmentation and Skeletonisation. Briefly, the MA width and feeding/draining vessel caliber are defined to be the largest and smallest radii registered along the skeleton of the MA, respectively.
AR quantifies the degree of asymmetry of the MA body. AR is defined as the ratio between the larger (A1) and smaller (A2) areas in the MA body mask to each side of the centerline (A1 divided by A2 in Figure 4, respectively, where A1>A2). Vorp et al. proposed a comparable measure of asymmetry for idealized aortic abdominal aneurysm (AAA) geometries and used it to characterize mechanical wall stress (Vorp et al., 1998). In subsequent work, Finol et al. studied the impact of AAA asymmetry on their hemodynamics and found that asymmetry tends to increase the maximum wall shear stress at peak flow and to induce the appearance of secondary flows in late diastole in idealized AAA geometries (Finol et al., 2003). AR is computed based on the MA body segmentation and centerline. Briefly, the polygon approximating the MA body is split into two along the MA centerline and the area of each sub-polygon is subsequently calculated. Custom Python scripts were developed to calculate BNR and AR.
Figure 4. AR is defined as the projected area on one side of the centerline (A1) divided by the area on the other side of the centerline (A2), where A1 > A2. This definition applies to both saccular (left) and fusiform (right) MAs as shown in this figure.
Hemodynamic Analysis
Based on the MA skeletonisation previously described and assuming rotational symmetry, we reconstructed the three-dimensional luminal surface of each MA under study (Figure 5). This surface encloses the approximate MA volume including its body and feeding/draining capillaries. The CFD package HemeLB (Bernabeu et al., 2014) was used to simulate both steady-state and time-dependent flow of a shear-thinning fluid modeled with the Carreau-Yasuda rheology model parametrized for human blood (Boyd and James, 2007). HemeLB uses the Lattice Boltzmann Method for the numerical simulation of blood flow. The interested reader can refer to (Aidun and Clausen, 2010; Krüger et al., 2017) for a complete presentation. The velocity field at the inlet was assumed to be parabolic (Poiseuille flow) for a given centerline peak velocity. To define this velocity, we took advantage of recent measurements of blood flow velocities in parafoveal capillaries by de Castro et al. (2016). Figure 4b of de Castro et al. (2016) reports velocity values over 4 cardiac cycles (equivalent to 3.13 s), which we used in the time-dependent flow simulations, with a mean capillary velocity of 1.69 mm/s, which we used in the steady-state flow simulations. Furthermore, no-slip velocity was imposed at the walls and a reference pressure was set at the outlet. To ensure that the flow field in the MAs is not affected by the finite length of the feeding/draining capillaries, we take them to be longer than the entrance length, Le, required for laminar flow to fully develop in a circular straight pipe. This is given by the expression Le = 0.035*D*Re (Bird et al., 2002), where D and Re are the diameter and Reynolds number of the feeding vessel, respectively. In all the MAs studied Le can be shown to be shorter than D. Therefore, the feeding/draining capillaries were segmented to be of length comparable to the MA body length along the flow axis for statistical purposes in the hemodynamic analyses that follow. Steady-state HemeLB simulations were run inexpensively in a four-core commodity hospital-based workstation, while time-dependent simulations made use of ARCHER, the UK National Supercomputing Service (http://www.archer.ac.uk). Typical execution times for the latter ranged between 4 and 10 h using 312 cores. All computational domains were discretized as a regular grid ensuring a minimum of 8 lattices sites across the narrowest point in the domain (Bernabeu et al., 2014) and comprised between 45,000 and 520,000 fluid lattices sites.
Figure 5. The image series shows the process for hemodynamic analysis. (A) The MA is imaged using AOSLO multiply scattered light imaging modality; (B) a perfusion map of the MA is created (calculated based on pixel-by-pixel standard deviation method) highlighting blood flow; (C) a binary mask of the MA is generated using the outline of the MA and its feeding and draining vessels; (D) a 3-D model of the MA is created under the assumption that it is rotationally symmetric; (E) flow streamlines, colored according to velocity magnitude, are plotted to show the paths followed by blood inside the MA, arrow indicates direction of flow.
Our computer simulations generated a description of the velocity, shear rate, and pressure fields in the whole computational domain as well as the wall shear stress on the model surface. In this study, we decided to characterize the changes in shear rate (SR) and wall shear stress (WSS) present in the MAs. Low SR has been associated with blood cell aggregation and clotting (Runyon et al., 2007) and abnormal WSS levels have been linked to endothelial cell dysfunction and changes in permeability (Tarbell, 2010). To reduce the dimensionality of the data and facilitate further statistical analysis we propose two indices for the characterization of the hemodynamic state of an MA: the shear rate mean drop (SRMD) and the wall shear stress mean drop (WSSMD). SRMD reports the ratio between the mean of the SR field in the MA feeding/draining vessels and the same measurement inside the MA body. Similarly, WSSMD indicates the ratio between the mean of the WSS on the MA feeding/draining vessels surface and the same measurement on the surface of the MA body. SRMD and WSSMD are dimensionless quantities. Finally, in the case of MA displaying clots, we also estimated SRMD for the clotted and perfused parts of the MA separately. Custom Python scripts were developed to calculate SRMD and WSSMD.
Study Cohort
In this study, 20 MAs were imaged from 13 eyes of 11 diabetic patients with varying severity of DR. The patient and MA characteristics are given in Supplementary Table 1. In this cohort, 9 of 11 patients had Type 1 diabetes, mean diabetes duration was 25 years and mean HbA1c was 8.1%. Informed written consent was obtained from each subject prior to the performance of any study procedures. This study adhered to the tenets of the Declaration of Helsinki and was approved by the institutional review board of the Joslin Diabetes Center.
Imaging Protocol and Light Safety
Mydriasis and cycloplegia were achieved by instillation of 1 drop each of 1% tropicamide (Akorn, Inc., Lake Forest, IL) and 2.5% phenylephrine hydrochloride (Akorn, Inc., Lake Forest, IL). Prior to AOSLO imaging, eye axial length (IOL Master, Zeiss, Germany) was measured to determine the magnification factor on AOSLO images. Ultrawide field, 200° digital fundus sphotographs (Optos 200Tx and Optos California, UK) were taken to determine MA location. During imaging, the subject's head was placed on a chin rest, and a head rest was used against the forehead for secure positioning. Precise head position adjustment and pupil alignment were achieved using a three-axis motorized stage (MT3-Z8 Thorlabs, NJ). MAs were imaged using AOSLO confocal imaging mode and multiply scattered light imaging mode with 75-frame videos. A 500 μm and 150 μm pinhole was used for forward scattering image and confocal imaging, respectively. Two SLDs were used for for imaging (830 nm) and wavefront sensing (780 nm). Output power at the cornea was 200 μW for the imaging SLD, and 70 μW for the wavefront sensing SLD. The light power was checked periodically to ensure compliance with the ANSI laser safety standard (American National Standards Institute, 2014).
Statistical Analysis
The segmentation of all the MAs and clotted regions are performed by at least 2 trained graders. For agreement between graders, <10% area variation for each MA and sub-region is ensured. All statistical analyses are completed using custom Python scripts and the Statistics package of the SciPy library (https://www.scipy.org). The Wilcoxon rank-sum test is used to test for significance in the comparison between groups. A p < 0.05 is used to reject the null hypothesis that two sets of measurements are drawn from the same distribution. Associations between continuous variables are evaluated using Pearson's correlation coefficient.
Results
Morphological and Hemodynamic Indices
Twenty MAs were imaged from 13 eyes of 11 diabetic subjects as shown in Figure 6, 10 were classified as saccular (5 partially clotted) and 10 as fusiform (none was clotted). For each MA, projected MA body size, asymmetry ratio (AR), body-to-neck ratio (BNR), shear rate mean drop (SRMD), and wall shear stress mean drop (WSSMD) are shown in Supplementary Table 2.
Figure 6. All MAs numbered according to decreasing AR. Images not to scale.
Analysis of Partially Perfused Mas
In the 5 partially perfused MAs (Figure 7), which had evidence of clot within the MA body, we calculated the hemodynamic indices within the perfused and clotted regions of the MA separately (Table 1).
Figure 7. AOSLO images of the 5 MAs where a clot was identified.
Table 1. Perfusion indices of the 5 MA within perfused versus clotted areas.
Among the partially clotted MAs, the SRMD and WSSMD values in the perfused regions were lower (mean ± SD: 63.36 ± 39.66 and 29.02 ± 16.74, respectively) than the values (211.85 ± 118.22 and 82.94 ± 30.91, respectively) in the clotted regions.
Asymmetry Ratio Predicts Manual MA Morphology Classification
All of the MAs in the study were qualitatively classified as saccular or fusiform according to the taxonomy proposed by Dubow et al. (2014). AR was calculated for all MAs and was found to be lower on average in the fusiform group compared to the saccular group (p < 0.001, Figure 8). Our data indicate that an AR threshold of ~1.5 reliably distinguishes fusiform from saccular MAs in this cohort. However, given the degree of overlap between both groups in terms of AR, it may not be advisable to define a unique cutoff value for automatic classification. Instead we propose a semi-automatic approach were MA with an AR below 1.4 and above 1.8 are automatically classified as fusiform and saccular, respectively, while those in the 1.4-1.8 region are labeled for manual classification by graders.
Figure 8. Asymmetry ratio index separated by MA groups.
Association of MA Morphology and Size
The area defined by the MA body segmentation, which is determined from an en face projection (xy plane) of the MA volume, was calculated for all the MAs in the study and used as a surrogate measure of MA size. Saccular MAs were found to be smaller than fusiform MAs (p = 0.004, Figure 9) with some saccular outliers having comparable size to the fusiform group. Moore et al. (1999) measured the extent of saccular and fusiform MAs in the direction perpendicular to the en face projection (z axis) and found no statistically significant difference. Taken together these results could indicate that size variability is more likely to be observed along the en face cross section compared to the transverse direction.
Figure 9. MA body size by MA groups.
Shear Rate Mean Drop Is Higher in MA Regions Likely to Clot
In the current study, all the MAs containing clots were of saccular type. MAs presenting clots appeared to have a higher AR approaching statistical significance in the comparison (p = 0.061, Figure 10).
Figure 10. Asymmetry ratio index in saccular MAs by clot presence.
Clots were always identified in contact with the MA wall. We performed hemodynamics analysis of the MAs to understand the relationship between flow and clot formation. The flow models were defined to include both the perfused and clotted portions of any given MA. In the 5 partially perfused MAs, we found that SRMD and WSSMD were higher in the regions where the clots were present compared to those that had not developed clots (p = 0.028 and p = 0.009, respectively, Figure 11). This speaks in favor of a model where MA thrombosis occurs in regions adjacent to the wall that experience low shear rates (hence high SRMD). In agreement with our results, low SR has been associated with blood clotting in vitro (Runyon et al., 2007) and with thrombus formation in intracranial aneurysms (Ribeiro de Sousa et al., 2016). Indeed, despite the obvious structural and hemodynamic differences between the macro and microcirculation, flow diverters, which rely on the principle of flow reduction from the parent circulation into the aneurysm body (hence SR reduction), are an established treatment for brain aneurysms (Jiang et al., 2016) to promote progressive intra-aneurysmal thrombosis.
Figure 11. Shear rate mean drop (SRMD, left) and wall shear stress mean drop (WSSMD, right) in partially perfused MAs by clotted/unclotted regions.
Body-To-Neck Ratio Correlates With Perfusion Changes in the MA Body
Both saccular and fusiform MAs are characterized by a sudden and non-uniform expansion of the vascular lumen. This change is most asymmetrical in the saccular class of MAs. This abnormal morphological configuration has a profound impact on the hemodynamics of the MA. We propose BNR as a simple metric for the quantification of hemodynamic abnormalities. Our results demonstrate that BNR is a good surrogate marker of SRMD (Pearson's r = 0.9, Figure 12) and WSSMD (Pearson's r = 0.83, Figure 12). Furthermore, mean WSSMD in this cohort was 35.2 with values as high as 78.4 (compared to a theoretical value of ~1 in the absence of MAs) showing the highly abnormal level of WSS experienced by endothelial cells lining the MA body wall compared to those in neighboring vessels.
Figure 12. Scatterplot of body-to-neck ratio vs wall shear stress mean drop index (left) and body-to-neck ratio vs. shear rate mean drop index (right) for each of the MAs studied. Regression lines and associated correlation coefficients are given to demonstrate the good correlation between the pairs of variables. The marginal histograms in each of the plots present the distribution of each of the variables studied.
Hemodynamic Changes Throughout the Cardiac Cycle
Blood flow displays pulsatile characteristics throughout the cardiac cycle. In our flow models, we can account for this property by defining a time-dependent inlet boundary condition based on the velocity traces measured by de Castro et al. (2016). Based on these simulations, we investigate the changes in velocity and shear rate throughout the cardiac cycle and their potential link with MA perfusion status and MA progression.
As expected, we find velocity and shear rate to be largest during systole, with regions that have developed clots experiencing reduced velocity and shear rate. We hypothesize that clots will form in areas of slow flow (i.e., low velocity) due to a sustained reduction in shear rate throughout the cardiac cycle (i.e., a low shear rate threshold). This is in agreement with in vitro studies looking at clot formation and propagation (Runyon et al., 2007). We calculate this threshold for the clotted region of MA1 to be ~1 s−1 on the previously described AOSLO delineation. In Supplementary Movie 1, we show the variation in the velocity field inside MA1 throughout the cardiac cycle and, color-coded in yellow, the regions of the MA experiencing a shear rate smaller or equal to 15 s−1. Interestingly, we observe how MA regions adjacent to the clotted part will fall below the threshold following systole (hence the yellow color disappear/appear in this region) when flow in the MA slows down.
Based on this observation, we postulate that a clot can propagate over time in areas where shear rate remains under threshold. We selected two MAs from the same eye for follow-up, one partially clotted at the time of baseline imaging (MA1) and another fully perfused (MA4). After 15 months of follow-up, the body of MA1 appeared to become non-perfused with persistent blood flow through a central vessel (Figure 13). Interestingly, the shape of MA4 remained unchanged and no clot development was observed.
Figure 13. Saccular (A) and fusiform (B) MAs from the same eye of a patient with severe NPDR. The shape and the perfusion status of the saccular MA changed dramatically, whereas the fusiform MA's shape and perfusion status was maintained during the 15-month non-treatment period.
Discussion and Conclusions
In the current work, we propose 4 novel indices for the classification and study of retinal MAs. Two of them are structural (asymmetry ratio, AR and body-to-neck ratio, BNR), and the other two describe the hemodynamic environment of the MA (shear rate mean drop, SRMD and wall shear stress mean drop, WSSMD). The limitations of the CFD methodology include the assumption of rotational symmetry in the MA surface reconstruction and the use of non-patient-specific boundary conditions. We calculated these indices in a set of 20 retinal MAs imaged with AOSLO. Our aim is to develop a method capable of establishing which MA characteristics are associated with a higher risk of leakage or clotting.
The data demonstrate that the proposed AR index is highly correlated with the qualitative MA classification of being either saccular or fusiform as performed by trained graders. The area calculated from the en face AOSLO projection of the MA body volume was found to be smaller in the saccular MAs studied compared to fusiform MAs.
It remains elusive why only some MAs are associated with retinal edema due to the disruption of endothelial cell barrier function. Previous work has linked abnormal WSS levels to endothelial cell dysfunction and changes in permeability (Tarbell, 2010). In the current work, we have proposed a method for the quantification of the changes in WSS experienced by the cells lining the MAs. Our results show a consistent WSS reduction with up to one order of magnitude difference among all cases (7- vs. 78-fold reduction). In future work, we will investigate the association between WSSMD and clinically observed MA leakage in longitudinal datasets. Furthermore, we shall investigate associations between the changes in WSSMD throughout the cardiac cycle and MA outcomes as changes in hemodynamic frequency have been shown to regulate pathologic phenotypes in endothelial cells (Feaver et al., 2012).
Previous studies have described and quantified the dynamic turnover of MAs in retinal vasculature (Goatman et al., 2003, Bernardes et al., 2009). In the current work, we took advantage of high resolution AOSLO imaging to observe partially clotted MAs. Five out of 20 MAs presented clots. All the partially clotted cases were of saccular type. Therefore asymmetry appeared to play a role in clotting. In one occasion, we could observe thrombosis of the MA body and remodeling of the affected capillary. Based on previous reports of the relationship between hemodynamics and blood clotting (Runyon et al., 2007) and thrombosis of vascular lesions (Ribeiro de Sousa et al., 2016), we studied SRMD and WSSMD in the MAs prior to clot development and identified a statistically significant reduction of both indices in the regions that would subsequently develop clots. Taken together, these results are consistent with the hypothesis that MA asymmetry promotes MA thrombosis through the well-characterized mechanism of blood clotting at low shear stress.
We anticipate that this work will shed light on the assessment of the dynamic processes of retinal MA development, clotting, and regression. We believe the proposed indices can be exploited as biomarker for vascular stability and DR disease progression. In future work, we will quantify this relationship and establish WSSMD/SRMD thresholds that facilitate the prediction of MA progression on a lesion-specific basis, as well as their relationship with MA leakage.
Ethics Statement
This study was carried out in accordance with the recommendations of the institutional review board of the Joslin Diabetes Center with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the institutional review board of the Joslin Diabetes Center.
Author Contributions
MB, YL, LA, and JS designed research. MB, YL, and OA-Q performed research. MB, YL, OA-Q, and JS analyzed data. MB, YL, OA-Q, LA, and JS wrote the manuscript.
Funding
National Eye Institute (1R01EY0-24702-01, 2R44EY-16295-04A1); Juvenile Diabetes Research Foundation (2-SRA-2014-264-M-R, 17-2011-359); Eleanor Chesterman Beatson Childcare Ambassador Program Foundation Grant; Massachusetts Lions Eye Research Fund; Engineering and Physical Sciences Research Council (EP/L00030X/1); Fondation Leducq (17 CVD 03). Funding support for author OA-Q comes from the Dubai Harvard Foundation for Medical Research.
Conflict of Interest Statement
JS received research support from Boston Micromachines Corp.
The remaining 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.
Acknowledgments
BEI research study coordinators (mydriasis), technicians. The authors acknowledge the contributions of the HemeLB development team. This work used the ARCHER UK National Supercomputing Service (http://www.archer.ac.uk).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2018.00989/full#supplementary-material
References
Aidun, C. K., and Clausen, J. R. (2010). Lattice-Boltzmann method for complex flows. Annu. Rev. Fluid Mech. 42, 439–472. doi: 10.1146/annurev-fluid-121108-145519
American National Standards Institute (2014). American National Standard for Safe Use of Lasers. Orlando, FL: American National Standards Institute.
Ashton, N. (1974). Vascular basement membrane changes in diabetic retinopathy. Montgomery lectur. Br. J. Ophthalmol. 58, 344–366.
Attali, D., and Montanvert, A. (1997). Computing and simplifying 2D and 3D continuous skeletons. Comput. Vis. Image Underst. 67, 261–273. doi: 10.1006/cviu.1997.0536
Bernabeu, M. O., Jones, M. L., Nielsen, J. H., Krüger, T., Nash, R. W., Groen, D., et al. (2014). Computer simulations reveal complex distribution of haemodynamic forces in a mouse retina model of angiogenesis. J. R. Soc. Interface 11:20140543. doi: 10.1098/rsif.2014.0543
Bernardes, R., Nunes, S., Pereira, I., Torrent, T., Rosa, A., Coelho, D., et al. (2009). Computer-assisted microaneurysm turnover in the early stages of diabetic retinopathy. Int. J. Ophthalmol. 223, 284–291. doi: 10.1159/000213638
Bird, R. B., Stewart, W. E., and Lightfoot, E. N. (2002). Transport Phenomena,. 2nd Edn New York, NY: John Wiley & Sons Inc.
Boyd, J., and James, M. B. (2007). Analysis of the casson and carreau-yasuda non-newtonian blood models in steady and oscillatory flows using the lattice boltzmann method. Phys. Fluids 19:093103. doi: 10.1063/1.2772250
Burns, S. A., Tumbar, R., Elsner, A. E., Ferguson, D., and Hammer, D. X. (2007). Large-Field-of-View, modular, stabilized, adaptive-optics-based scanning laser ophthalmoscope. J. Opt. Soc. Am. A Opt. Image Sci. Vis. 24, 1313–1326. doi: 10.1364/JOSAA.24.001313
Chien, A., and Sayre, J. (2014). Morphologic and hemodynamic risk factors in ruptured aneurysms imaged before and after rupture. AJNR Am. J. Neuroradiol. 35, 2130–2135. doi: 10.3174/ajnr.A4016
Chien, A., Sayre, J., and Viñuela, F. (2011). Comparative morphological analysis of the geometry of ruptured and unruptured aneurysms. Neurosurgery 69, 349–356. doi: 10.1227/NEU.0b013e31821661c3
Chui, T. Y., Gast, T. J., and Burns, S. A. (2013). Imaging of vascular wall fine structure in human retina using adaptive optics scanning laser ophthalmoscopy. Invest. Ophthalmol. Vis. Sci. 54, 7115–7124. doi: 10.1167/iovs.13-13027
Chui, T. Y., Vannasdale, D. A., and Burns, S. A. (2012). The use of forward scatter to improve retinal vascular imaging with an adaptive optics scanning laser ophthalmoscope. Biomed. Opt. Express 3, 2537–2549. doi: 10.1364/BOE.3.002537
Cogan, D. G., Toussaint, D., and Kuwabara, T. (1961). Retinal vascular patterns. IV. Diabetic retinopathy. Arch. Ophthalmol. 66, 366–378.
Curtis, T. M., Gardiner, T. A., and Stitt, A. W. (2009). Microvascular lesions of diabetic retinopathy: clues towards understanding pathogenesis? Eye (Lond) 23, 1496–1508. doi: 10.1038/eye.2009
de Castro, A., Huang, G., Sawides, L., Luo, T., and Burns, S. (2016). Rapid high resolution imaging with a dual-channel scanning technique. Opt. Lett. 41, 1881–1884. doi: 10.1364/OL.41.001881
Dhar, S., Tremmel, M., Mocco, J., Kim, M., Yamamoto, J., Siddiqui, A. H., et al. (2008). Morphology parameters for intracranial aneurysm rupture risk assessment. Neurosurgery 63, 185–197. doi: 10.1227/01.NEU.0000316847.64140.81
Dubow, M., Pinhas, A., Shah, N., Cooper, R. F., Gan, A., Gentile, R. C., et al. (2014). Classification of human retinal microaneurysms using adaptive optics scanning light ophthalmoscope fluorescein angiographyclassification of retinal microaneurysms. Invest. Ophthalmol. Vis. Sci. 55, 1299–1309. doi: 10.1167/iovs.13-13122
Duh, E. J., Sun, J. K., and Stitt, A. W. (2017). Diabetic retinopathy: current understanding, mechanisms, and treatment strategies. JCI Insight 2:93751. doi: 10.1172/jci.insight.93751
ETDRS_No10 (1991). Grading diabetic retinopathy from stereoscopic color fundus photographs—an extension of the modified Airlie House classification. ETDRS report number 10. Early treatment diabetic retinopathy study research group. Ophthalmology 98, 786–806.
ETDRS_No12 (1991). Fundus photographic risk factors for progression of diabetic retinopathy. ETDRS report number 12. Early treatment diabetic retinopathy study research group. Ophthalmology 98, 823–833.
Ezra, E., Keinan, E., Mandel, Y., Boulton, M. E., and Nahmias, Y. (2013). Non- dimensional analysis of retinal microaneurysms: critical threshold for treatment. Integr. Biol. 5, 474–480. doi: 10.1039/c3ib20259c
Feaver, R. E., Gelfand, B. D., and Blackman, B. R. (2012). Human haemodynamic frequency harmonics regulate the inflammatory phenotype of vascular endothelial cells Nat. Commun. 4:1525. doi: 10.1038/ncomms2530
Finol, E. A., Keyhani, K. K., and Amon, C. H. (2003). The effect of asymmetry in abdominal aortic aneurysms under physiologically realistic pulsatile flow conditions. ASME J. Biomech. Eng. 125, 207–217. doi: 10.1115/1.1543991
Goatman, K. A., Cree, M. J., Olson, J. A., Forrester, J. V., and Sharp, P. F. (2003). Automated measurement of microaneurysm turnover. Invest. Ophthalmol. Vis. Sci. 44, 5335–5341. doi: 10.1167/iovs.02-0951
Hirai, F. E., Moss, S. E., Knudtson, M. D., Klein, B. E., and Klein, R. (2007). Retinopathy and survival in a population without diabetes: the Beaver Dam Eye Study. Am. J. Epidemiol. 166, 724–730. doi: 10.1093/aje/kwm126
Jiang, B., Paff, M., Colby, G. P., Coon, A. L., and Lin, L. M. (2016). Cerebral aneurysm treatment: modern neurovascular techniques. Stroke Vasc. Neurol. 1, 93–100. doi: 10.1136/svn-2016-000027
Kempen, J. H., O'Colmain, B. J., Leske, M. C., Haffner, S. M., Klein, R., Moss, S. E., et al. (2004). The prevalence of diabetic retinopathy among adults in the United States. Arch. Ophthalmol. 122, 552–563. doi: 10.1001/archopht.122.4.552
Klein, B. E. (2007). Overview of epidemiologic studies of diabetic retinopathy. Ophthalmic Epidemiol. 14, 179–183. doi: 10.1080/09286580701396720
Ko, F., Vitale, S., Chou, C. F., Cotch, M. F., Saaddine, J., and Friedman, D. S. (2012). Prevalence of nonrefractive visual impairment in US adults and associated risk factors, 1999-2002 and 2005-2008. JAMA 308, 2361–2368. doi: 10.1001/jama.2012.85685
Krüger, T., Kusumaatmaja, H., Kuzmin, A., Shardt, O., Silva, G., and Viggen, E. M. (2017). The Lattice Boltzmann Method. Principles and Practice. Springer.
Lu, Y., Bernabeu, M. O., Lammer, J., Cai, C. C., Jones, M. L., Franco, C. A. et al. (2016). Computational fluid dynamics assisted characterization of parafoveal hemodynamics in normal and diabetic eyes using adaptive optics scanning laser ophthalmoscopy. Biomed. Opt. Express 7, 4958–4973. doi: 10.1364/BOE.7.004958
Moore, J., Bagley, S., Ireland, G., McLeod, D., and Boulton, M. E. (1999). Three dimensional analysis of microaneurysms in the human diabetic retina. J. Anat. 194, 89–100.
Murakami, T., Nishijima, K., Sakamoto, A., Ota, M., Horii, T., and Yoshimura, N. (2011). Foveal cystoid spaces are associated with enlarged foveal avascular zone and microaneurysms in diabetic macular edema. Ophthalmology 118, 359–367. doi: 10.1016/j.ophtha.2010.03.035
Nunes, S., Pires, I., Rosa, A., Duarte, L., Bernardes, R., and Cunha-Vaz, J. (2009). Microaneurysm turnover is a biomarker for diabetic retinopathy progression to clinically significant macular edema: findings for type 2 diabetics with nonproliferative retinopathy. Int. J. Ophthalmol. 223, 292–297. doi: 10.1159/000213639
Ribeiro de Sousa, D., Vallecilla, C., Chodzynski, K., Corredor Jerez, R., Malaspinas, O., Eker, O. F., et al. (2016). Determination of a shear rate threshold for thrombus formation in intracranial aneurysms J. Neurointerv. Surg. 8, 853–858. doi: 10.1136/neurintsurg-2015-011737
Ribeiro, M. L., Nunes, S. G., and Cunha-Vaz, J. G. (2013). Microaneurysm turnover at the macula predicts risk of development of clinically significant macular edema in persons with mild nonproliferative diabetic retinopathy. Diabetes Care 36, 1254–1259. doi: 10.2337/dc12-1491
Runyon, M. K., Johnson-Kerner, B. L., Kastrup, C. J., Van Ha, T. G., and Ismagilov, R. F. (2007). Propagation of blood clotting in the complex biochemical network of hemostasis is described by a simple mechanism. J. Am. Chem. Soc. 129, 7014–7015. doi: 10.1021/ja072602p
Tam, J., Martin, J. A., and Roorda, A. (2010). Noninvasive visualization and analysis of parafoveal capillaries in humans. Invest. Ophthal- mol. Vis. Sci. 51, 1691–1698. doi: 10.1167/iovs.09-4483
Tarbell, J. M. (2010). Shear stress and the endothelial transport barrier. Cardiovasc. Res. 87, 320–330. doi: 10.1093/cvr/cvq146
Vorp, D. A., Raghavan, M. L., and Webster, M. W. (1998). Mechanical wall stress in abdominal aortic aneurysm: influence of diameter and asymmetry. J. Vasc. Surg. 27, 632–639.
Webb, R. H., Albanese, M. J., Zhou, Y., Bifano, T., and Burns, S. A. (2004). Stroke amplifier for deformable mirrors. Appl. Opt. 43, 5330–5333. doi: 10.1364/AO.43.005330
Wilkinson, C. P., Ferris, F. L., Klein, R. E., Lee, P. P., Agardh, C. D., Davis, M., et al. (2003). Proposed international clinical diabetic retinopathy and diabetic macular edema disease severity scales. Ophthalmology 110, 1677–1682. doi: 10.1016/S0161-6420(03)00475-5
Keywords: diabetic retinopathy, microaneurysm, adaptive optics, blood flow, computational fluid dynamics
Citation: Bernabeu MO, Lu Y, Abu-Qamar O, Aiello LP and Sun JK (2018) Estimation of Diabetic Retinal Microaneurysm Perfusion Parameters Based on Computational Fluid Dynamics Modeling of Adaptive Optics Scanning Laser Ophthalmoscopy. Front. Physiol. 9:989. doi: 10.3389/fphys.2018.00989
Received: 15 December 2017; Accepted: 05 July 2018;
Published: 07 September 2018.
Edited by:
Peter V. Coveney, University College London, United KingdomReviewed by:
Paolo Melillo, Università degli Studi della Campania “Luigi Vanvitelli” Caserta, ItalyMariano Vázquez, Barcelona Supercomputing Center, Spain
Jazmin Aguado-Sierra, Barcelona Supercomputing Center, Spain
Copyright © 2018 Bernabeu, Lu, Abu-Qamar, Aiello and Sun. 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: Jennifer K. Sun, jennifer.sun@joslin.harvard.edu
†These authors have contributed equally to this work