Rodent retinal circulation organization and oxygen metabolism revealed by visible-light optical coherence tomography

: Visible light optical coherence tomography (vis-OCT) is an emerging label-free and high-resolution 3-dimensional imaging technique that can provide retinal oximetry, angiography, and flowmetry in one modality. In this paper, we studied the organization of the arterial and venous retinal circulation in rats using vis-OCT. Arterioles were found predominantly in the superficial vascular plexus whereas veins tended to drain capillaries from the deep capillary plexus. After that, we determined the oxygen metabolic rate supported by retinal microcirculation by combining retinal vessel oxygen saturation and blood flow measurements. The ability to visualize and monitor retinal circulation organization and oxygen metabolism by vis-OCT may provide new opportunities for understanding the pathology of ocular diseases.

the retinal angiography and retinal flowmeter by Doppler OCT [13,14], which had been accomplished previously in near infrared OCT. Although the technique has been successfully applied to human retina [15,16] and other tissue [17] imaging, its great potential on retinal imaging by combing all these functional features have not been fully investigated. We have established a fiber-based vis-OCT with high axial resolution [18] and previously performed automated retinal oximetry with good reproducibility [19]. Our hypothesis for this study was that vis-OCT could enable noninvasive visualization of the organization of the retinal circulation and measurement of oxygen metabolism supported by retinal circulation. We proposed to combine a full set of functional features (i.e. retinal oximetry, angiography and flowmetry) of vis-OCT to verify this hypothesis in the Brown Norway rat.

Animal preparation
Brown Norway rats (17 weeks old, n = 8) were initially anesthetized with 5% isoflurane in a sealed box for 10 minutes, followed by 2.5% isoflurane during the imaging session. Pure oxygen (100% O 2 ) with flow rate at 1 L/min was ventilated for inhalation. Under anesthesia, the animal was immobilized in a custom-made imaging stage with multidimensional manipulation for alignment. The pupil was dilated with 1% tropicamide ophthalmic solution before imaging. To keep the cornea moist, sterile irrigating salt solution (Alcon Laboratories Inc.) was applied to the eye every other minute. The animal's body temperature was maintained with a 38.5 °C water-warming blanket. The exhaust gas, including isoflurane, was removed by a vacuum pump to avoid carbon dioxide accumulation and collected by an anesthesia gas filter (OMNICON F/air, Bickford [20]) before releasing to the open air.
All experimental procedures complied with the Association for Research in Vision and Ophthalmology Statement for the Use of Animals in Ophthalmic and Vision Research and were approved by the Institutional Animal Care and Use Committee (IACUC) of the Oregon Health & Science University (OHSU).

Image acquisition
The vis-OCT used in this work is a custom-built prototype [18] housed in the Center for Ophthalmic Optics & Lasers Laboratory (COOL LAB) [21] of OHSU. The spectrum covered the high absorption contrast region of hemoglobin from 510 nm to 610 nm (λ c = 560 nm, full width at half maximum (FWHM) was ~90nm) and was calibrated with a neon calibration light source (NE-1, Ocean Optics [22]). The dispersion mismatch between the two arms was compensated both physically and numerically. An unbalanced 90:10 wideband fiber optical coupler delivered 10% power into the sample arm (power = 0.8 mw). A telescope tube (f 1 = 75 mm, f 2 = 11 mm) ( Fig. 1(A)) guided light into the rat eye, with rat head immobilized on a device allowing oxygen inhalation and carbon dioxide removal ( Fig. 1(B)). The system operated at a flexible axial scan sampling rate and scanning angle, with a 1.2-μm axial resolution and approximately 6-μm lateral resolution with an image depth of 1.8 mm. The maximum sensitivity was measured as 89 dB with a protected silver mirror. Volumetric raster scans were collected at ± 3° scanning angle at 2-D galvo mirror, resulting in a 3.2 mm × 3.2 mm image field on the rat retina (d eyeball = 6.2 mm). In the fast transverse (X) scanning direction, 512 axial scans were sampled to form a B-scan, resulting in the sampling density of 6.25 µm/line. Three consecutive B-scans were captured before moving to the next slow axis (Y) position. A total of 512 Y positions were sampled within 17 seconds at the A-scan rate of 50 kHz and B-scan rate of 88 Hz.
To obtain the wide field reflectance and angiographic images shown here, the optic disc was first placed at the center of the field of view and a volumetric raster scan was taken for reference. Then the rat was tilted and its optic disc was placed at each of the corners of field of view consecutively to take volumetric raster scans on the inferior-nasal, inferior-temporal, superior-temporal, superior-nasal regions.
In addition, two concentric circular scans (r 1 = 0.6 mm, r 2 = 0.8 mm) consisting of 4096 A-lines (~1 µm/line) were also acquired around the optic disc, transecting all major retinal arteries and veins for the purpose of blood flow measurement and arteriovenous identification by flow direction. The scans were repeated 30 times to account for arterial blood flow differences during systole and diastole of the cardiac cycle. Three repeated scans were acquired within minutes to reduce measurement variation.

Oxygen saturation
The spatial-and depth-resolved optical density OD (z, λ) is defined as the logarithm of the ratio of the reflected intensity spectrum I (z, λ) to the source spectrum I 0 (λ). It stands for the reflectivity ability at a certain illumination wavelength for the tissue at a certain depth and can be expressed as Eq. (1) based on a modified Beer's law [10]. Basically, d is the accumulated absorption length for the voxels in depth. The scattering spectrum of the vessel wall r(λ) can be modeled as a power law A• λ -α under the first-order Born approximation [23]. The scattering spectrum at the reference arm R 0 is considered as a wavelength-independent constant. The subscript HbO 2 and Hb indicate the contribution from oxygenated and deoxygenated hemoglobin respectively, with their extinction coefficients ε referring to the literature [12] and concentrations C being calculated by fitting, in order to find the oxygen saturation sO 2 = C HbO2 / (C HbO2 + C Hb ).
The OD (z, λ) can be extracted by short time Fourier transform (STFT) spectroscopic analysis of the interference fringes. Vessel posterior voxels were selected automatically for analysis [19]. A Gaussian window with a full-width at half-maximum of approximately 9 nm and an interval distance around 3 nm were applied, resulting in 21 split spectral bands in total. Only those bands within the contrast region from 527 nm to 582 nm were selected for linear regression fitting.

Angiography
The recorded data were analyzed with the SSADA algorithm [5] for angiography. Briefly, the algorithm splits the full spectrum into multiple bands and calculates the decorrelation value among consequent B-scans at each band. After that, the contrast of the motion-based angiogram is largely enhanced by averaging the results from all bands. Layer segmentation was done with a graph-search technique [24] on structural B-scan images. The laminar vascular/capillary plexuses were then generated by projecting the flow signal within specific slabs. En face angiograms at the four scanning regions were montaged together for a wide field of view (~60 degree). Note that only a circular region with the diameter of 5 mm was highlighted for visualization.
Projected en face angiograms from the superficial vascular plexus was enhanced by a Frangi multiscale vessel filter [25]. Vessel binary masks were obtained by setting a threshold to the enhanced angiogram. After that, morphologic processing was performed to retain only vessel information. After tracing the vascular network for each major vessel, the vessel binary mask could be further isolated into an artery and a vein binary mask according to the measured oxygen saturation in each major vessel. Vessel density was calculated separately for arteries and veins as the ratio of the vascular area to the total area.

Blood flow
The blood flow F (μL/min) (Eq. (2)) for each major vessel was obtained by the product of the averaged vessel blood flow velocity ν (mm/s) and its cross-section area S that is perpendicular to the flow direction [13] F S ν =  (2) The distribution of blood flow velocity can be expressed as Eq. (3) and obtained from the phase-sensitive Doppler OCT analysis to the circular scans. Briefly, the Δφ is the Doppler phase shift between adjacent A-lines in the B-scan with ultra-high sampling density. The sampling rate f s of the charged coupled device (CCD) camera is 50 kHz. The center wavelength of the light spectrum λ 0 is 565 nm. Refractive index n of the retina is taken as 1.4 [26]. The Doppler angle θ [rad] was calculated from the vessel position shifts between the inner circular scan and outer circular scan. Vessel diameter H was measured axially from Bscan images.

Retinal oxygen metabolic rate
Oxygen metabolic rate is defined as the oxygen consumption of the tissue over time, which can be deduced from the difference of oxygen content supplied by arteries and carried away by veins per unit time [27]. In the inner retina, the oxygen demand is mainly supplied by the retinal circulation, where all vessels originate from the central retinal artery and vein, which enter the eye with the optic nerve. Thus the rMRO 2 (g/min) can be quantified as Where W O2 and W HbO2 are the molecular weights of O 2 [32 g/mol] and oxygenated hemoglobin [68000 g/mol] respectively. One oxygenated hemoglobin can bind 4 oxygen molecules. The total hemoglobin concentration in blood C HbT is 150 g/L for rat. The sO 2 and F are the measured oxygen saturation and blood flow in each vessel respectively. The subscripts a and v denote the contribution of oxygen transported in from the artery and transported out by the vein respectively, i indicates vessel index and the difference between the two transport rates is considered to be the amount consumed by retinal metabolism.

Organization of retinal circulation
From the cross-sectional B-scan and the averaged A-line profile (Fig. 2(A)), the distinct retinal structural layers and three vascular plexuses were observed. In Brown Norway rats, the retinal microvasculature originates from the central retinal artery and central retinal vein, which enter the eye just inferior to the optic nerve. They both branch into several major retinal vessels in a "spoke wheel" distribution radiating towards the retinal periphery ( Fig.  2(B)-2(C)). The rat eye lacks a fovea, and therefore has no foveal avascular zone.
The capillaries in the ICP ( Fig. 2(D)) and DCP ( Fig. 2(E)) follow curvilinear and tortuous paths just as they do in humans, and are evenly distributed over the imaging area. Each plexus possesses a unique distribution pattern, indicating a successful illustration of the laminar vascular plexuses with vis-OCT. Compared to the human retina [2], the vascular patterns are less dense and do not demonstrate radial peripapillary capillaries ( Fig. 2(C)).
On the SVP angiogram with isolated vessel trees color-coded with oxygen saturation of major vessels (Fig. 2(F)), we found major arteries (sO 2 > 93%, red) were located on the top of major veins (sO 2 < 93%, green), and they are interleaved. Arterioles and their branches were predominately distributed extensively within the SVP. The minimal vessel calibers in the en face images in this study were about ~12 μm, which is similar to the diameter of smallest capillaries. Branches down to the 5th order could be visualized and, even at this level, relatively few direct connections between arterial and venous branches were observed (e.g. white arrow in Fig. 2(C)). Precapillary arterioles tend to branch in a dichotomous way before disappearing from the SVP slab. By contrast, fewer venules were present in the SVP and they tended to travel only a short distance within this slab. This difference could be demonstrated by measuring vessel density in vessel binary masks of arteries and veins separately. In this study, SVP vessel densities were ~17% and ~6% for arteries and veins, respectively.
The reason veins have less perfusion area than arteries in the SVP is that their branches tend to vertically connect with DCP, which can be better visualized by overlaying the SVP on the DCP (Fig. 2(G)). In this representative case, nineteen remarkable venule interconnections were indicated by arrows, illustrating that veins primarily drain blood from the DCP. Some venules are connected to veins vertically, while others are oblique and travel a certain path before being collected by major veins in the SVP.

Oxygen metabolism
Through the dual circles Doppler scans (Fig. 3(A)-3(C)), we observed that both arteries and veins can be distinguished based on their flow directions (Fig. 3(D)), consistent with the differentiation made by oxygen saturation. It should be noted that the locations for Doppler OCT should be carefully selected, because the angle between the OCT beam and direction of blood flow is needed and the opposite flow direction is required for this purpose. We also found the values of the velocity are higher in the center of vessels than in the regions near to the vessel walls, in agreement with the expected parabolic blood flow velocity distribution across large vessels.
Retinal blood oxygen transport rate (rTRO 2 ) was obtained by determining the amount of oxygen transportation per unit time through sO 2 and blood flow (Fig. 3(E)-3(F)). The retinal oxygen metabolic rate (Fig. 3(G)) is measured as 261.7 ± 49.2 ng/min from n = 8 eyes, which can be converted to 187.9 ± 35.3 nL/min by the ideal gas law [27]. The measured values suggest that rMRO 2 in hyperoxia condition is equivalent to that in normal air and lower than that in hypoxic conditions [28].

Discussion
New functional OCT modalities are allowing researchers to gain better insight into the pathophysiology of ocular diseases. Previous OCTA research suggested that ischemia is part of the causal chain in many leading causes of blindness [29]. However, how the retinal circulation is organized, how oxygenation differs between retinal plexuses and which plexuses are more vulnerable to ischemia are still unclear. Since vis-OCT can provide maps of the retinal circulation in three dimensions together with vessel oxygen saturation with excellent reproducibility [19], it stands as an instrument with a great potential to help address these questions.
One of the current debates in which vis-OCT exhibits potential to elucidate some clarity is that of whether retinal plexuses drain through venous flow independently, or if this drainage is primarily through the DCP. Previously, both vertical inter-plexus and horizontal intraplexus arteriovenous connections had been observed by our group in human retinas with projection-resolved OCTA techniques [30]. Yet a debate on the organization of retinal circulation still exists [31], and a hybrid pattern has gradually become more convincing, with more details recently revealed by OCTA [32]. However, given the very dense vascular network (60%~80% vessel density in SVP) in the human retina, the architecture of its microvasculature remains unclear at limited resolution.
Here in the rat eye, the retinal microvasculature could be visualized more clearly using our vis-OCT technique for four main reasons. First, vis-OCT enables rapid and accurate identification of arterioles and venules by oxygen saturation. This helps differentiate the vascular type so that we can delineate the arterial and venous networks separately. Second, the axial and transverse resolutions of our prototype vis-OCT are improved with respect to conventional infrared OCT [18]. Third, the vascular patterns in the rat retina are relatively sparse (~25% vessel density in SVP) compared to human retina. Fourth, with a greater distance from the optic disc reached by ultra-wide field angiography, more vascular features can be visualized.
Researchers have previously explored the 3D connectivity of rat retina circulation using OCT angiography [33]. They demonstrated the dense connectivity between the ICP and DCP by tracing the capillaries. However, they failed to reveal the full landscape of retina vascular architecture due to lack of 1) arteriovenous differentiation results 2) accurate retina layers segmentation, and 3) wide-field of view. From our results above, three features of the rat retinal circulation organization can be concluded. First, the arterial network predominates in the SVP. Second, direct arteriovenous connections in the SVP were seldom observed (white arrow in Fig. 2(C)). Third, that veins drain blood flow primarily through the DCP, where arborescent patterns connecting directly into major venules could be appreciated. The findings are similar to the previously reported ex vivo observations in pig and mouse retinas by confocal microscopy [31,34]. Although the capillary connections in the ICP were hard to trace, a reasonable conjecture that the capillaries in the ICP connect to the diving end of arterioles in the SVP could be made through a limited field angiogram centered at the optic disc (Fig. 4(A)-4(D)) and a fly through video (Visualization 1). Based on the findings above, we suggest a predominantly blood flow transition pathway through the DCP in the rat retina ( Fig. 4(E)), i.e. the primary flow transition within the microvasculature architecture would be in a way of artery (SVP) → arteriole (SVP) → capillary (ICP and DCP) → venule (DCP) → vein (SVP). Nevertheless, direct connections to venules were occasionally observed for capillaries in the ICP at the same point where underlying DCP capillaries drain. This observation, together with the few direct arteriovenous connections observed in the SVP suggest that although serial organization is predominant, the complex vascular distribution is more accurately described by a hybrid model (Fig. 4(E), flow transition schematic). Anatomic differences between species are to be expected, particularly when considering those with a fovea and parafoveal vasculature, and several models of microvascular connectivity in the human retina have been discussed. Some authors have proposed a serial model where venous drainage occurs primarily through vortex-like vessels in the DCP [35]. Others have identified frequent connections between arterioles and venules at each of the retinal capillary plexuses in human retina. The situation suggests a predominantly parallel organization within a hybrid connectivity model [32], resembling the more complicated function of human retina. Additionally, the existence of focal, plexus specific ischemia such as nerve fiber layer infarcts (cotton wool spots) and paracentral acute middle maculopathy would be difficult to explain with a pure serial model in humans and further supports the hybrid model. The high axial resolution of vis-OCT angiography may help further clarify this anatomy.
In addition to vascular mapping, noninvasive measurement of blood flow with vis-OCT in all major retinal vessels allows for calculation of total retinal blood flow, which was 6.5-8 μL/min per retina in this study with 100% oxygen inhalation. While this value is significantly affected by rat size and anesthesia protocol, our finding fits well with reported value ranges from previous invasive and noninvasive studies [36,37]. There are many potential clinical applications of vis-OCT determined flowmetry. Measuring flow changes over time or between eyes may be useful in the early diagnosis of impending retinal vascular occlusions or ocular ischemic syndrome. A decrease in blood flow in a retinal tumor feeder vessel could be a marker of successful treatment response. Differences in flow between the systolic and diastolic phases of the cardiac cycle may provide information about vessel stiffness and systemic atherosclerosis.
Besides giving access to the arteriovenous flow distribution in the retina, this technology also allowed us to quantify the retinal oxygen metabolic rate, a potential early biomarker for vision-threatening retinal diseases such as diabetic retinopathy, age-related macular degeneration or central retinal vein occlusion. Vis-OCT could also be used to address fundamental pathophysiology questions such as what are the initial hemodynamic dysregulations and early metabolism changes in diabetic retinopathy, where there is current debate due to conflicting data [38,39]. Furthermore, the combination of assessing vascular anatomy, flow rates and oxygen delivery might reveal and explain flow redistributions in diseases of the retinal microvasculature, help explain the dynamics of oxygen supply in the event of unmet metabolic needs, and shed light into the natural history of complications like neovascularization. The outstanding versatility and non-invasiveness of vis-OCT make it a potentially valuable technology for future exploration of these questions.
Safety of exposure to visible light remains a major concern for vis-OCT. Delori et al have shown that the visible light is safe at limited exposure duration under certain radiant power [40]. Recently, human retina imaging using vis-OCT was achieved using a much reduced laser power [41,42]. However, this work was performed on normal subjects only. Further improvements in instrument sensitivity and/or advanced image processing will be required in order to diminish discomfort while preserving a reasonable imaging quality before this technology can be applied to patients. The question of whether the light exposure of vis-OCT will directly alter retinal oxygen metabolism or blood flow must also be addressed when considering clinical applications. Previous work suggests that any effect may be minimized at the high scan frequency of vis-OCT and in eyes that are already light adapted [43,44].
Currently, oxygen saturation is only assessable in major retinal vessels. Though it is sufficient to calculate the total oxygen metabolism, the lack of capillary oxygen saturation limits the usefulness of distribution maps of oxygen metabolism, which might be altered in very early stages of ocular diseases. Successful assessment of capillary oxygen saturation currently faces two main challenges: 1) red blood cells pass through the capillary one by one, causing a scattering spectrum that is highly dependent on their geometry, orientation and size [45]. 2) There is limited hemoglobin absorption along the axial direction in capillaries, which is not sufficient to provide enough contrast for oxygen saturation measurement. With the development of vis-OCT technology, these issues may be solved by a more effective oxygen saturation algorithm and/or better laser sources.

Conclusion
In summary, three functional features, retinal angiography, oximetry and flowmetry have been combined in the vis-OCT system to image retinal circulation of Brown Norway rats. Using oximetry to help differentiate arteries from veins, we observed a predominance of arterioles in the superficial vascular plexus, whereas venules (the branches of major veins) tended to drain blood from the DCP. This vascular organization suggested a predominantly serial pattern of blood flow transition from arteries to veins in rat's retinal circulation. In addition, the retinal oxygen metabolic rate was also measured by combing retinal oximetry and retinal flowmetry. With these capabilities, future developments based on vis-OCT have potential to play a crucial role in better understanding abnormalities of retinal microvasculature anatomy and oxygen metabolism in various ocular diseases and disorders.