Measuring oxygen saturation in retinal and choroidal circulations in rats using visible light optical coherence tomography angiography

Siyu Chen; Ji Yi; Hao F. Zhang

doi:10.1364/BOE.6.002840

1. Introduction

Optical coherence tomography (OCT) demonstrated tremendous success in imaging three dimensional (3D) tissue anatomy with micrometer-level spatial resolution on a macroscopic length scale [1]. In the past two decades, the technique has been translated into various clinical practices [2]. In particular, OCT has joined the standard of care in ophthalmology, providing rich 3D anatomical images from, for example, various retinal layers, choroid, and optic nerve head. Besides the technical advances on high-speed and high-resolution retinal imaging, functional imaging of various hemodynamic parameters has emerged as an invaluable supplement to anatomical OCT imaging. Perhaps the most important developments of functional OCT in retinal imaging are the measurement of retinal blood flow [3–9], and label-free retinal OCT microangiography (OMAG) [10–12]. Because the neuroretina has a high metabolic demand, the blood flow in the retinal circulation can be an excellent indicator of the viability of the retina and has been associated to various retinal diseases [13–15].

Another important hemodynamic parameter of the retinal circulation is the hemoglobin oxygenation saturation or sO₂. The sO₂ difference between the arterial and venous blood vessels indicates oxygen extraction from the circulation. When combined with blood flow measurements, oxygen metabolism in the retina can be quantified by gas volume of oxygen consumed per unit time [16, 17]. The attempts of using OCT to measure blood oxygenation dated back more than a decade ago, where the spectroscopic contrast around the isosbestic point around 800 nm from oxy- and deoxy-hemoglobin was characterized [18–21]. Since then, several groups also explored the oxygenation dependent scattering coefficient to create contrast for oxygenated or deoxygenated whole blood. However, none of the efforts has progressed into successful measurements of absolute sO₂ in vivo. The reason is that the absorption contrast around 800 nm was overwhelmed by the wavelength-dependent scattering coefficient of red blood cells/whole blood, which have later been theoretically and numerically demonstrated [22, 23]. Knowing that hemoglobin has much higher absorbing contrast within visible-light spectral range than within NIR range, Yi et. al. and Robles et. al. explored the feasibility of using visible-light OCT (vis-OCT) to quantify sO₂ from ex vivo and in vitro samples, respectively [24, 25]. They further showed in vivo measurement of blood oxygenation [26], and the first demonstration of in vivo retinal oximetry in rodents using vis-OCT combing with the spatially-resolved spectroscopic analysis [27]. Later, to isolate the blood signal from the static tissue in OCT signal, we proposed to combine the speckle contrast and the spectroscopic contrast of blood simultaneously [28] to calculate sO₂ in microvasculature. We demonstrated that sO₂ could be directly calculated from wavelength-dependent OCT angiography without the cumbersome image processing previously needed to locate blood vessels. A similar strategy was recently adapted for both exposed cortical imaging and optical-power-cancelled retinal imaging in rodents [29]. However, the theoretical treatment did not consider the complication of blood flow to the angiography signal, which could confound the blood oxygenation calculation in vivo.

Here, we expanded our theoretical formulation by taking blood flow contrast into account in OCT angiography-based oximetry. Experimental verification was conducted by pumping anticoagulated whole bovine blood through a capillary tube. We tested the influence of different flow rates to the spectral derivative (i.e. the first-order partial derivative of spatially resolved OCT spectrum against wavelength). Further, we proposed a new data processing algorithm to normalize different flow rate and yield a robust sO₂ calculation. We tested our method in rat retina in vivo and detected sO₂ variation upon inhaling different gas mixtures. Moreover we demonstrated, for the first time, a non-invasive measurement of the relative sO₂ changes in choriocapillaris using vis-OCT.

2. Theory

In this section, we described the extended theoretical formulation for the spectroscopic analysis on vis-OCT angiography-based oximetry. The following motion-enhanced dynamic scattering model is to demonstrate that a linear relationship can be established between sO₂ and OCT angiographic derived spectrum measurements. Without loss of generality, we can start the formulation from a wavelength-dependent A-line signal from OCT angiogram $A A (z, λ)$ , where z defines the axial coordinate and λ denotes the wavelength. Note that we will omit the finite axial resolution and sensitivity roll-off for simplicity purpose. At any given point along the A-line,

A A (z, λ) = I_{0} (λ) r (z, λ) η (z, λ) G (Δ φ),

where

I_{0} (λ)

is incidence intensity.

r (λ)

is the wavelength-dependent back scattering coefficient of local tissue, which can be further expressed by a power law

r (λ) = B λ^{- α}

under the first Born approximation.

B

is a scaling constant while the value of

α

depends on the tissue type and has been previously characterized to be ~1 with variations around ± 0.2 [30, 31]. The term

G (Δ φ)

is a non-negative scalar representing the dynamic contrast enhancement caused by random phase variation due to motion. The value of G is assumed to be 0 only if the tissue is static.

η (z, λ)

is the cumulative attenuation when light double passes the tissue,

η (z, λ) = e^{- \int_{0}^{z} 2 μ_{t} (z^{'}, λ) d z'},

where

μ_{t}

is the total attenuation coefficient of local tissue attributed to both optical absorption and scattering. For whole blood, it is a linear combination of absorption coefficient

μ_{a}

, scattering coefficient

μ_{s}

, and a packing factor W,

μ_{t} = μ_{a} + W μ_{s}

[22, 32]. For other retinal tissue, reduced total attenuation coefficient,

μ_{t} = μ_{a} + (1 - g) μ_{s}

may be used, where g denotes the anisotropy factor of the scattering process.

After normalizing $A A (z, λ)$ by the source spectrum $I_{0} (λ)$ , the resulted image is analogous to the wavelength-dependent optical density (OD) mapping used in traditional fundus retinal oximetry, except for one additional motion enhancement term $G (Δ φ)$ . Hence we use the same notation and denote the normalized $A A (z, λ)$ as $O D (z, λ)$ ,

O D (z, λ) = \ln \frac{A A (z, λ)}{I_{0} (λ)} = \ln r (λ) + \ln η (μ_{a}, μ_{s}, λ) + \ln G (Δ φ),

or

O D (z, λ) = α λ + \ln B - 2 \int_{0}^{z} μ_{t} (z^{'}, λ) d z' + \ln G (Δ φ) .

As we take the partial derivative of OD against wavelength (spectral derivative), we have

\frac{\partial O D (z, λ)}{\partial λ} = - \frac{2 \partial \int_{0}^{z} μ_{t} (z^{'}, λ) d z'}{\partial λ} + α + \frac{\partial \ln G (Δ φ)}{\partial λ} .

If we assume that the wavelength is much smaller than the motion displacement of the blood cells, the phase variation becomes random and the enhancement factor

G (Δ φ)

loses wavelength dependence. The spectral derivative of

G (Δ φ)

could be considered to be zero and the equation becomes

\frac{\partial O D (z, λ)}{\partial λ} = - \frac{2 \partial \int_{0}^{z} μ_{t} (z^{'}, λ) d z'}{\partial λ} + α .

Now let us consider the attenuation contribution from the whole blood and the static tissues. As light propagates in tissue, the cumulative attenuation is

\int_{0}^{z} μ_{t} (z^{'}, λ) d z' = \int_{0}^{l_{b}} μ_{t} (z^{'}, λ) d z'_{b l o o d} + \int_{0}^{l_{e l s e}} μ_{t} (z^{'}, λ) d z'_{s t a t i c t i s s u e}

or

\int_{0}^{z} μ_{t} (z^{'}, λ) d z' = μ_{t_b} l_{b} + \bar{μ_{t_e l s e}} l_{e l s e},

where

μ_{t_b}

and

\bar{μ_{t_e l s e}}

denote the total attenuation coefficients of whole blood and averaged attenuation coefficient of static retinal tissue, respectively;

l_{b}

and

l_{e l s e}

denote the cumulative optical path length within the blood vessel and in other retinal tissues, respectively.

Based on the previous characterization of $μ_{t_b}$ and $\bar{μ_{t_e l s e}}$ , we found that the spectral derivative of $\bar{μ_{t_e l s e}}$ was close to zero [32, 33]. As a result, the spectral derivative of Eq. (7) is

\frac{\partial \int_{0}^{z} μ_{t} (z^{'}, λ) d z'}{\partial λ} = \frac{\partial μ_{t_b}}{\partial λ} l_{b} + \frac{\partial \bar{μ_{t_e l s e}}}{\partial λ} l_{e l s e} ≅ \frac{\partial μ_{t_{b}}}{\partial λ} l_{b} .

We can further express

μ_{t_b}

with respect to sO₂ as

μ_{t_b} = s O_{2} \cdot μ_{t_{O x y}} + (1 - s O_{2}) \cdot μ_{t_{D e o x y}},

where

μ_{t_O x y}

and

μ_{t_D e o x y}

are the total attenuation coefficients for fully oxygenated and deoxygenated whole blood, respectively. For simplicity, we define the following two constants,

k_{O x y} = \frac{\partial μ_{t_O x y}}{\partial λ}

and

k_{D e o x y} = \frac{\partial μ_{t_D e o x y}}{\partial λ}

and plug Eqs. (8)-(10) into Eq. (6) so that

\frac{\partial O D (z, λ)}{\partial λ} = - s O_{2} \cdot (k_{O x y} - k_{D e o x y}) l_{b} + k_{D e o x y} l_{b} + α .

It is clear that Eq. (11) establishes a linear relationship between the spectral derivative of the OD function and sO₂, while

(k_{O x y} - k_{D e o x y}) l_{b}

serves as a constant scaling factor. During our experiment, l_b is a constant determined by the choice of ROI. Its exact value can be estimated by measuring the mean optical path on the OCT angiograph.

3. Methods

3.1 Animal preparation

We anesthetized wild-type Long Evans rats using Ketamine/Xylazine cocktail solution (11.45 mg Ketamine and 1.71 mg Xylazine per milliliter of solution, respectively). The solution was administrated via intraperitoneal injection (IP) at a dosage of 87 mg Ketamine and 13 mg Xylazine per kilogram of body weight. A supplementary dose (30% of the original volume) was administrated after 30 minutes of the initial injection to maintain the animal at deep anesthesia. We applied 0.5% Tetracaine hydrochloride ophthalmic solution for local eye anesthesia and 1% Tropicamide ophthalmic solution for pupil dilation. During the entire period of anesthesia, we maintained the core body temperature at 37 °C using an electronic heating pad with feedback control (FHC Inc). The animal was held on a custom-made imaging stage and properly immobilized. A pulse oximeter attached to the rear limb of the animal provided real-time monitoring of its heart rate and systemic arterial oxygen saturation.

The animal was first ventilated with normal air (21% O₂, 79% N₂) at 4.2 Standard Liter per Minute (SLPM). Then we supplied the following gas mixtures in the order of: (1) 5% carbon dioxide (21% O₂, 74% N₂, and 5% CO₂); (2) pure oxygen (100% O₂); and (3) 10% oxygen (10% O₂ and 90% N₂). For each gas inhalation, the animal was allowed to stabilize for at least 3 minutes before the measurement was taken. Additional time was allowed if the pulse oximeter reading showed large fluctuation (>2% deviation over 10 seconds). We supplied normal air for 3 minutes between changes. We took one additional measurement after the sequence with normal air ventilation. For each OCT measurement, the corresponding oximeter readings were also recorded simultaneously.

All the experimental procedures were approved by the Northwestern University IACUC and conformed to the Association for Research in Vision and Ophthalmology Statement on Animal Research.

3.2 Image acquisition

We developed a spectral-domain vis-OCT system to image rodent retina in vivo. A supercontinuum laser source (SuperK EXTREME EXW-6, NKT Photonics) provided broadband illumination. The output beam was delivered to the input arm of a free-space Michaelson interferometer via fiber optics. The beam splitter (CM1-BS013, Thorlabs) of the interferometer divided and redirected the beam to the reference arm and the sample arm. In the reference arm, a set of carefully chosen glass plates compensated for dispersions created by the unbalanced optical components. A continuously adjustable ND filter (NDC-50C-2M, Thorlabs) attenuated the beam intensity to enhance the dynamic range of the interference signal. The sample arm consists of a pair of scanning galvanometer mirror (Nutfield Technology) and a Keplerian telescope (0.2 × magnification) to project rectangular raster scanning pattern on the rat fundus. The spectrometer covered a spectrum range from 520 nm to 630 nm. We estimated the transverse and axial resolutions of our vis-OCT system to be 15 μm and 1.7 μm, respectively. A fast CMOS line scan camera (spL2048-140km, Basler) digitized the interference spectrum for image reconstruction. All acquisition sequence and recording were synchronized by a customized LabView (National Instrument) virtual instrument program.

We adopted uni-directional scans in our OCT angiography. Here we denote the fast-scan direction as the x axis and slow scan direction as the y axis. During one B-scan, we translated the laser spot along the positive x axis at 50-kHz A-line rate and a total of 400 A-lines were recorded. Immediately after each B-scan, we returned (flied back) the laser spot to the x axis origin without acquiring any signal. The returning speed was limited to four times the forward scanning speed, which was constrained by the mechanical properties of the galvanometer. We repeated the process five times at each y location before moving to the next one. For one vis-OCT angiography image stack, we had 5 × 512 = 2560 B-scans at 512 different y locations, each B-scan consisted of 400 A-lines. The entire OCT image stack took 25.6 seconds to acquire.

3.3 Vis-OCT imaging of retinal microvasculature

We generated vis-OCT angiography using phase-sensitive decorrelation algorithm [11, 28, 34]. The contrast originated from the flowing blood cells inside vessels as they randomly change the phase and amplitude of back-scattered light. The random phase shift caused decorrelation in OCT signal and can be enhanced via high-pass filtering. However, such enhancement may be confounded by bulk motion, which also contributed to the signal via the same mechanism. These motion originated from normal physiological movements (e.g. heart beat and breathing) and the rotating galvanometer mirrors, both of which were inevitable. We must take measures to minimize their effects so as not to interfere with vis-OCT angiography. We adopted a three-step approach to correct the aforementioned bulk motion decorrelation within each y location: (1) Bulk image shift; (2) Axial global phase fluctuations (AGPF) correction; and (3) Lateral global phase fluctuations (LGPF) correction [35]. Without loss of generality, we used the notation $A (x, z, t)$ to represent the recovered complex vis-OCT signal of the five consecutive B-scans at a given y location; x and z represents the lateral and axial location, respectively; t indicates the specific B-scan in time order, whose value ranges from 1 to 5 in our case.

The first step was to co-register the five real-valued OCT structure image at each y location, so that motion shift greater than the imaging resolution can be corrected. We calculated the cross-correlation function between the two adjacent B-scan images (t + 1 and t). We identified the maximum value of the cross-correlation function along x and z directions, which indicates when the two images were properly aligned. The second image was then shifted back so the two images were properly registered.

The second step was to correct the global phase fluctuations with two phase modifiers AGPF and LGPF on the co-registered OCT images. The AGPF can be estimated by integrating the phase variation along the A-line axial direction (z axis) as

A G F P (x) = - a r g [\int^{​} A^{*} (x, z, t_{2}) A (x, z, t_{1}) d z],

while the LGPF can be estimated by integrating the phase variation along the x axis as

L G F P (z) = - a r g [\int^{​} A^{*} (x, z, t_{2}) A (x, z, t_{1}) d x] .

We multiplied the two phase modifiers AGFP and LGFP to the corresponding complex OCT signal pixel-wise. The corrected complex valued vis-OCT signal is denoted $A' (x, z, t)$ . The B-scan angiography cross-section $A G (x, z, t)$ at location y can then be calculated by

A G (x, z) = E | A' (x, z, t_{τ + 1}) - A' (x, z, t_{τ}) |,

where

E

represents the expected value; and

| \cdot |

represents the modulus operation. We iterated the process at each y location and generated a three dimensional image stack

A G (x, y, z)

.

At last, en face angiogram was generated by integrating the first 15 pixels of the highest motion contrast along z direction.

3.4 sO₂ estimation

We performed spectroscopic angiography to obtain the spectral profile of motion enhanced contrast of blood vessels [10, 11, 28], which was used to calculate the absolute blood sO₂. Instead of using the full spectrum, we employed short-time Fourier transform (STFT) to perform spatial-spectral dual-domain analysis. We empirically chose four Gaussian-shaped STFT windows centered at 565 nm, 568 nm, 571 nm, and 574 nm, each with 10 nm full-width-half-maximum (FWHM) bandwidth. The full spectrum was windowed using these narrow-band STFT windows, respectively. For each center wavelength chosen, we generated the corresponding narrow-band OCT angiography using the method described in the above section. The calculated wavelength-dependent OCT angiography is denoted $A G (x, y, z, λ)$ . For any given vessel at depth z, we first performed image segmentation to isolate a thin layer containing the vessel. To project the image stack into a two-dimensional en face angiograph $A G (x, y, z, λ)$ , we integrated the first 15 superficial pixels within the blood vessel. The integration interval of 15 was empirically chosen to optimize signal-to-noise ratio. We then calculated the mean value of $A G (x, y, z, λ)$ within the x-y plane inside the vessel region. The value was then assigned to $A A (z, λ)$ to characterize one designated vessel, as in Eq. (1). The process is repeated for all four center wavelength windows.

Next, we normalized $A A (z, λ)$ using the corresponding source spectrum I₀(λ) to calculate the localized $O D (z, λ)$ . We performed a linear regression to calculate the slope of the straight line that went through the center of the four OD points. We used this slope to estimate the partial derivative of the OCT angiograph spectrum against wavelength.

To convert the spectral derivative into absolute sO₂ values, we performed a calibration based on pulse oximetry spO₂ readings where we assumed that the major retinal arteries had the same sO₂ as systemic arterial sO₂. Based on the calibration on each animal, we can measure the absolute sO₂ values in the major retinal arteries and venules around the optic disk.

We summarized the data processing steps in Table 1.

Table 1. Summary of data processing procedures to extract sO_2.

View Table

4. Results

4.1 Phantom experiment verification

To demonstrate that the spectral derivative method is not affected by flow velocity, we designed and performed the following in vitro experiment. We fabricated a capillary phantom using a polystyrene tube (inner diameter: 125 µm) filled with anticoagulated bovine whole blood (Quad Five). We fully oxygenated the blood by exposing it to normal air at room temperature for 1 hour prior to the experiment. The phantom was embedded beneath 2.5% agarose gel to mimic a blood vessel in the retina. We controlled the blood flow rate in the tube using a syringe pump (PHD 2000 Infusion, Harvard Apparatus) at various pre-set values within physiological range.

In order to image the in vitro sample, we slightly modified our OCT system. An objective lens (LSM03-VIS, Thorlabs) replaced the telescope system to create a focusing spot about 5.5 µm in diameter. The same data acquisition and processing protocols were applied.

Figure 1(a) compares the cross-sectional view of OCT structural images and motion-enhanced images. It can be clearly seen that surrounding agarose gel has strong scattering background in conventional OCT. In comparison, the motion-enhanced images largely suppressed the static background, while the moving blood within the tube was enhanced due to dynamic scattering. Even when the flow rate was set to zero, intrinsic random Brownian motion was sufficient to enhance the contrast (the left most images in Fig. 1(a)). On the sectioned en face view of the OCT angiograph showing the capillary phantom, we took four rectangular region of interest (ROI’s) of equal size along its longitudinal axis (Fig. 1(b)). We extracted the wavelength-dependent motion contrast spectrum within the selected ROI’s to evaluate the influence of flow velocity. As expected, the shapes of these spectra were comparable, showing characteristic absorption spectrum of oxygenated whole blood. The spectrum calculated from zero-flow phantom is shown in Fig. 1(c).

Fig. 1 Spectroscopic analysis of in vitro phantom experiment. (a) Cross-sectional images of the tube based on direct (top row) and motion-enhanced dynamic scattering (bottom row) contrasts; (b) En face OCT angiography image of the tube. The yellow rectangles indicate the ROI used; (c) Representative motion contrast spectrum calculated from selected ROI’s, flow rate = 0.00 mL/min; (d) Averaged intensity I_m taken from four ROI’s within the capillary tube; (d) The spectral derivative of motion contrast spectrum obtained within the capillary tube with varying flow rates from 560 to 580 nm. Error bar = S.D.

Download Full Size | PDF

We extracted two parameters to characterize the obtained motion contrast spectrum, the averaged intensity (I_m) and spectral derivative (∂I/∂λ) from 565 nm to 574 nm. Their values were plotted against increasing flow rate in Figs. 1(d) and 1(e), respectively. We performed liner regressions to estimate the slopes of both curves. I_m showed a clear decreasing trend with a negative slope, which is statistically significant (p < 0.001). On the contrary, the slope of the ∂I/∂λ curve did not deviate from zero (p = 0.74). The result showed that ∂I/∂λ was flow-independent across the tested range, suggesting that sO₂ estimation based on ∂I/∂λ would not be affected by different flow rate.

4.2 Measuring rat retinal sO₂ in vivo

Figure 2 compares the cross-section view of the OCT images in rat retina in vivo. Figure 2(a) is the conventional structural OCT image. The anatomical features of the retina can be clearly distinguished. The cross sectional OCT angiogram is shown in Fig. 2(b). We further overlaid the structure image with the angiography in Fig. 2(c), showing locations of each vessel within the retina. Besides the vasculatures in retinal circulation, we also found strong blood motion contrast directly beneath the RPE, which suggested the existence of blood flows in choriocapillaris. Beyond that, visible light attenuated rapidly.

Fig. 2 (a) Example of cross-sectional images of conventional OCT in rat retina in vivo using vis-OCT. Anatomical layers are labeled including, NFL/GCL-nerve fiber layer/ganglion cell layer, IPL-inner plexiform layer, INL-inner nuclear layer, OPL-outer plexiform layer, ONL-outer nuclear layer, ELM-external limiting membrane, IS/OS-inner segments/outer segments, and RPE-retinal pigment epithelium. (b) The corresponding OCT microangiograph generated using motion contrast. Vessels in retinal circulation and choroicapillaries are identified. Artifacts created by the hyper-reflective layers in IS/OS and RPE are indicated with a blue arrow. (c) Overlaid image from (a) and (b). The choroicapillaris lay beneath RPE layer.

Download Full Size | PDF

We then compared the en face view of rat fundus obtained by two OCT modality. Figure 3(a) illustrates the conventional structural OCT image. Although the superficial vessels can be clearly observed, we cannot distinguish the microvasculature. In comparison, Fig. 3(b) showed motion-enhanced dynamic scattering contrast of the same data set, where the hyper-reflective signals from the static tissues are largely suppressed. Due to limited lateral resolution, choriocapillaris appeared as a single bright layer and was excluded through segmentation. Further, we calculated the sO₂ of the major retinal vessels in the retinal circulation during normal oxygen inhalation, as the color-encoded angiographic image in Fig. 3(b). We can clearly distinguish the arteries and veins based on their individual sO₂. Finally, we rendered the 3D vascular structure in Fig. 3(c). Besides the dominating brunches of central retinal arteries and veins, laminated structures of micro-vasculatures within nerve fiber layer (NFL) and outer plexiform layer (OPL) are also distinguishable. The fly-through en face microangiogram is shown in the movie.

Fig. 3 In vivo imaging of a rat retina using vis-OCT. (a) and (b) are en face projection of amplitude and motion enhanced microangiography, respectively. In the microangiography map, each major retinal arterioles and venules are color-encoded based on the measured sO₂ value. (c) A three dimensional rendering of rat retinal microvasculature. See Visualization 1 for the fly-through en face microangiogram.

Download Full Size | PDF

4.3 Absolute sO₂ measurement in retinal circulation

We next performed in vivo verification of our spectral derivative method. We calibrated the absolute sO₂ values from the major retinal arteries against a pulse oximeter readings under different ventilation conditions (normal air, 5% CO₂, pure O₂, and 10% O₂), as shown in Fig. 4(a). The calibration process corresponded to the Step 5 in Table 1. The result shows a strong positive correlation (R² = 0.87) between the two independent measurements, indicating that the calibration offered reasonably good estimation of the true sO₂.

Fig. 4 Measurement of absolute sO₂ values in the retinal circulation in vivo. (a) Measured arterial sO₂ calibrated against pulse oximetry readings (R² = 0.87) from n = 8 rats. Black line is the linear regression of the data (b) Measured sO₂ of major retinal arteries and venues upon changes of inhalation gas mixtures (n = 8 for each bar). Error bar = S.E.M. *: p < 0.05. **: p < 0.01.

Download Full Size | PDF

Based on the established calibration curve, we measured the major retinal arterial and venous sO₂ in respond to the change of inhalation gases. In the course of the experiment, we changed the inhalation gases in the order of normal air, 5% CO₂, pure oxygen, and 10% oxygen. As shown in Fig. 4(b), the arterial and venous sO₂ were significantly different (p < 0.05), except when pure oxygen was inhaled and the shunting between arterioles and venules elevated the venous sO₂. Besides the arteriovenous sO₂ difference, we observed a moderate increase (p = 0.06) in both arterial and venous sO₂ when we supplied 5% CO₂ gas. This increased sO₂ is expected because 5% CO₂ has the effect of vasodilation, which improves tissue oxygenation [36, 37]. Both the arterial and venous sO₂ increased significantly (p < 0.05) under pure O₂ and decreased significantly (p<0.01) under 10% O₂.

4.4 Relative sO₂ change in choriocapillaris

Given the fact that we had strong flow-contrast enhanced signals from choriocapillaris directly underneath the RPE, we attempted to measure the relative changes of sO₂ in choriocapillaris. Due to limited lateral resolution, we could not distinguish individual choriocapillaris; instead, we averaged the signal over the entire area of choriocapillaris. We performed image segmentation to remove any area shadowed by superficial blood vessels, and yielded one averaged spectrum $A A (z_{c}, λ)$ , where z_c was the axial position in the choroid.

In order to eliminate the contamination from the retinal circulation in the inner retina, we normalized the choriocapillaris OCT angiograph spectrum by the intensity spectrum from the layer directly above RPE and calculated the spectral derivative. We used the same calibration curve obtained from retinal circulation to obtain the relative change of sO₂ with respect to normal air breathing, as in Fig. 5. When we supplied 5% CO₂ gas mixture, the average sO₂ in choriocapillaris was observed to increase moderately by 2% (p = 0.29). When we supplied pure O₂ and 10% O₂ gas mixture, the average sO₂ was observed to increase by 3.5% (p < 0.001) and decreased by 5% (p < 0.001), respectively.

Fig. 5 Measured relative choriocapillaris sO₂ change comparing to normal air breathing (n = 6). ***: p < 0.0001.

Download Full Size | PDF

4. Discussion

In this work, we presented an extended theoretical formulation for OCT angiography-based oximetry. This method used the flow contrast in OCT angiography, based on which we normalized the wavelength-dependent angiography and provided a robust measurement of sO₂. The approach differs from our previous vis-OCT oximetry publication [27], where only static OCT amplitude information was used. The introduction of dynamic contrast eliminates the need for cumbersome vessel segmentation, and has arguably higher sensitivity and noise resistance than our previous method. More importantly, the new approach also enables us to measure the relative change of sO₂ in chororiocapillaries, which are generally inaccessible in structural vis-OCT images. We observed the change of sO₂ in both retinal circulation and choriocapillaris under different inhaling gas mixtures, which include normal air, 5% CO₂, 100% O₂, and 10% O₂. Our in vivo measurements were verified by independent and concurrent pulse oximetry.

In our dynamic scattering model, we assumed that the motion enhancement term (G) is wavelength-independent. The assumption requires that the phase fluctuation is random between 0 to 2π without any preference. The condition is automatically satisfied when the displacement of scattering particles during the scanning dwell time are much larger than the wavelength. In our case, the red blood cell travels about 0.13 mm during the 0.01 s dwell time, which suffices the requirement.

We calculated the first derivative of the spectral range from 560 nm to 580 nm to establish its correlation with sO₂. It differs from the least square fitting method over much larger wavelength range (530 nm-600 nm) as previously reported by our and other groups [28, 29]. The major benefit is that the permitted incidence power in the current study could be reduced by about 2/3, which is ~0.3mW in the current experiments. By removing unnecessary exposure, the safety margin for vis-OCT for human can be significantly improved and the discomfort can be dramatically reduced. By sacrificing the spectral information from the absorption contrast, our method requires independent oxygenation measurements for calibration purpose. Fortunately, pulse oximetry utilizes the pulsatile profile and provides the reading of the overall arterial oxygen saturation, which should be equal to the values in the major retinal arterioles immediately exiting optic nerve head. This allows us to calibrate our spectral derivatives as we showed in Fig. 4(a) in every experiment. Such a calibration is easily achievable using a pulse oximetry that is widely available in the clinical settings.

Interestingly, we were also able to measure sO₂ variations in choriocapillaris right beneath the RPE under different inhalation conditions. The choroidal circulation is very different from retinal circulation, both anatomically and functionally [38]. The choroidal blood is supplied through posterior ciliary arterioles and drained through vortex vein. To our concern, the major difference between the two circulations is that the choroidal circulation has very small arteriovenous sO₂ difference, whereas the retinal circulation has a large arteriovenous sO₂ difference [39, 40]. This results in an almost uniform sO₂ level in the choriocapillaris under normal physiology conditions, which allows us to measure sO₂ over a relatively large area of choriocapillaris.

However, there are a few complications that prevent us providing absolute sO₂ measurements in choriocapillaris. First, the true sO₂ level in choriocapillaris is still unknown. The most accurate reference is the direct measurement of oxygen tension using an invasive oxygen sensitive microelectrode [38]. Although it is commonly assumed that the choroidal circulation has the systemic oxygenation as in the arterioles, the microelectrode measurements turned out to be lower than expected [38]. Without a reliable reference, we were not able to calibrate our spectral derivative as we did for the retinal vessels. Alternatively, one might consider using the least square fitting method with the full spectral range of the hemoglobin absorption contrast to characterize choriocapillaris sO₂. However, the optical properties of dispersed single blood cells in capillary network are not entirely clear, which would bring another uncertainty in the least square fitting method. Second, light passes through entire retina before choriocapillaris. The accumulated attenuation is also wavelength dependent, which would affect the spectral derivative. Nonetheless, we found that this is the first time a non-invasive measurement can be implemented to measure relative choriocapillaris sO₂ variations.

5. Conclusion

In summary, we presented a new theoretical formulation for vis-OCT angiography-based oximetry. The flow contrast was used and characterized by in vitro experiments. After proper spectral normalization, we demonstrated, for the first time, robust in vivo measurements of absolute sO₂ variation in the retinal circulation and relative sO₂ variation in the choriocapillaris. The new method requires much reduced illumination power, which greatly benefits clinical applications of vis-OCT.

Acknowledgments

We would like to acknowledge the following funding resources: National Institute of Health 1R01EY019951 and 1R24EY022883 and National Science Foundation grants CBET-1055379, CBET-1066776, and DBI-1353952. J. Yi is supported in part by a postdoctoral fellowship award from the Juvenile Diabetes Research Foundation. H. F. Zhang has financial interests in Opticent Health Inc., which, however, did not support this work.

References and links

1. J. G. Fujimoto, “Optical coherence tomography for ultrahigh resolution in vivo imaging,” Nat. Biotechnol. 21(11), 1361–1367 (2003). [CrossRef] [PubMed]

2. A. M. Zysk, F. T. Nguyen, A. L. Oldenburg, D. L. Marks, and S. A. Boppart, “Optical coherence tomography: a review of clinical development from bench to bedside,” J. Biomed. Opt. 12(5), 051403 (2007). [CrossRef] [PubMed]

3. B. White, M. Pierce, N. Nassif, B. Cense, B. Park, G. Tearney, B. Bouma, T. Chen, and J. de Boer, “In vivo dynamic human retinal blood flow imaging using ultra-high-speed spectral domain optical coherence tomography,” Opt. Express 11(25), 3490–3497 (2003). [CrossRef] [PubMed]

4. Y. Wang, B. A. Bower, J. A. Izatt, O. Tan, and D. Huang, “In vivo total retinal blood flow measurement by Fourier domain Doppler optical coherence tomography,” J. Biomed. Opt. 12(4), 041215 (2007). [CrossRef] [PubMed]

5. B. Baumann, B. Potsaid, M. F. Kraus, J. J. Liu, D. Huang, J. Hornegger, A. E. Cable, J. S. Duker, and J. G. Fujimoto, “Total retinal blood flow measurement with ultrahigh speed swept source/Fourier domain OCT,” Biomed. Opt. Express 2(6), 1539–1552 (2011). [CrossRef] [PubMed]

6. Z. Zhi, X. Yin, S. Dziennis, T. Wietecha, K. L. Hudkins, C. E. Alpers, and R. K. Wang, “Optical microangiography of retina and choroid and measurement of total retinal bloodflow in mice,” Biomed. Opt. Express 3(11), 2976–2986 (2012). [CrossRef] [PubMed]

7. Z. Zhi, W. Cepurna, E. Johnson, T. Shen, J. Morrison, and R. K. Wang, “Volumetric and quantitative imaging of retinal blood flow in rats with optical microangiography,” Biomed. Opt. Express 2(3), 579–591 (2011). [CrossRef] [PubMed]

8. R. M. Werkmeister, N. Dragostinoff, S. Palkovits, R. Told, A. Boltz, R. A. Leitgeb, M. Gröschl, G. Garhöfer, and L. Schmetterer, “Measurement of Absolute Blood Flow Velocity and Blood Flow in the Human Retina by Dual-Beam Bidirectional Doppler Fourier-Domain Optical Coherence Tomography,” Invest. Ophthalmol. Vis. Sci. 53(10), 6062–6071 (2012). [CrossRef] [PubMed]

9. C. Dai, X. Liu, H. F. Zhang, C. A. Puliafito, and S. Jiao, “Absolute Retinal Blood Flow Measurement With a Dual-Beam Doppler Optical Coherence Tomography,” Invest. Ophthalmol. Vis. Sci. 54(13), 7998–8003 (2013). [CrossRef] [PubMed]

10. L. An and R. K. Wang, “In vivo volumetric imaging of vascular perfusion within human retina and choroids with optical micro-angiography,” Opt. Express 16(15), 11438–11452 (2008). [CrossRef] [PubMed]

11. Y. Jia, O. Tan, J. Tokayer, B. Potsaid, Y. Wang, J. J. Liu, M. F. Kraus, H. Subhash, J. G. Fujimoto, J. Hornegger, and D. Huang, “Split-spectrum amplitude-decorrelation angiography with optical coherence tomography,” Opt. Express 20(4), 4710–4725 (2012). [CrossRef] [PubMed]

12. D. Y. Kim, J. Fingler, R. J. Zawadzki, S. S. Park, L. S. Morse, D. M. Schwartz, S. E. Fraser, and J. S. Werner, “Optical imaging of the chorioretinal vasculature in the living human eye,” Proc. Natl. Acad. Sci. U.S.A. 110(35), 14354–14359 (2013). [CrossRef] [PubMed]

13. Y. Wang, A. A. Fawzi, R. Varma, A. A. Sadun, X. Zhang, O. Tan, J. A. Izatt, and D. Huang, “Pilot Study of Optical Coherence Tomography Measurement of Retinal Blood Flow in Retinal and Optic Nerve Diseases,” Invest. Ophthalmol. Vis. Sci. 52(2), 840–845 (2011). [CrossRef] [PubMed]

14. Y. Jia, S. T. Bailey, D. J. Wilson, O. Tan, M. L. Klein, C. J. Flaxel, B. Potsaid, J. J. Liu, C. D. Lu, M. F. Kraus, J. G. Fujimoto, and D. Huang, “Quantitative Optical Coherence Tomography Angiography of Choroidal Neovascularization in Age-Related Macular Degeneration,” Ophthalmology 121(7), 1435–1444 (2014). [CrossRef] [PubMed]

15. S. Srinivas, O. Tan, S. Wu, M. G. Nittala, D. Huang, R. Varma, and S. R. Sadda, “Doppler Fourier-domain optical coherence tomography measured Retinal blood flow in Chinese American adults: The Chinese American Eye Study,” Invest. Ophthalmol. Vis. Sci. 55(5), 214 (2014).

16. W. Song, Q. Wei, W. Liu, T. Liu, J. Yi, N. Sheibani, A. A. Fawzi, R. A. Linsenmeier, S. Jiao, and H. F. Zhang, “A combined method to quantify the retinal metabolic rate of oxygen using photoacoustic ophthalmoscopy and optical coherence tomography,” Sci. Rep. 4, 6525 (2014). [CrossRef] [PubMed]

17. J. Yao, K. I. Maslov, Y. Zhang, Y. Xia, and L. V. Wang, “Label-free oxygen-metabolic photoacoustic microscopy in vivo,” J. Biomed. Opt. 16(7), 076003 (2011). [CrossRef] [PubMed]

18. D. J. Faber, E. G. Mik, M. C. G. Aalders, and T. G. van Leeuwen, “Light absorption of (oxy-)hemoglobin assessed by spectroscopic optical coherence tomography,” Opt. Lett. 28(16), 1436–1438 (2003). [CrossRef] [PubMed]

19. C.-W. Lu, C.-K. Lee, M.-T. Tsai, Y.-M. Wang, and C. C. Yang, “Measurement of the hemoglobin oxygen saturation level with spectroscopic spectral-domain optical coherence tomography,” Opt. Lett. 33(5), 416–418 (2008). [CrossRef] [PubMed]

20. D. J. Faber, E. G. Mik, M. C. G. Aalders, and T. G. van Leeuwen, “Toward assessment of blood oxygen saturation by spectroscopic optical coherence tomography,” Opt. Lett. 30(9), 1015–1017 (2005). [CrossRef] [PubMed]

21. L. Kagemann, G. Wollstein, M. Wojtkowski, H. Ishikawa, K. A. Townsend, M. L. Gabriele, V. J. Srinivasan, J. G. Fujimoto, and J. S. Schuman, “Spectral oximetry assessed with high-speed ultra-high-resolution optical coherence tomography,” J. Biomed. Opt. 12(4), 041212 (2007). [CrossRef] [PubMed]

22. D. J. Faber, M. C. G. Aalders, E. G. Mik, B. A. Hooper, M. J. C. van Gemert, and T. G. van Leeuwen, “Oxygen Saturation-Dependent Absorption and Scattering of Blood,” Phys. Rev. Lett. 93(2), 028102 (2004). [CrossRef] [PubMed]

23. S. Chen, J. Yi, W. Liu, V. Backman, and H. Zhang, “Monte Carlo Investigation of Optical Coherence Tomography Retinal Oximetry,” IEEE Trans. Biomed. Eng. (Accepted). [PubMed]

24. F. E. Robles, S. Chowdhury, and A. Wax, “Assessing hemoglobin concentration using spectroscopic optical coherence tomography for feasibility of tissue diagnostics,” Biomed. Opt. Express 1(1), 310–317 (2010). [CrossRef] [PubMed]

25. J. Yi and X. Li, “Estimation of oxygen saturation from erythrocytes by high-resolution spectroscopic optical coherence tomography,” Opt. Lett. 35(12), 2094–2096 (2010). [CrossRef] [PubMed]

26. F. E. Robles, C. Wilson, G. Grant, and A. Wax, “Molecular imaging true-colour spectroscopic optical coherence tomography,” Nat. Photonics 5(12), 744–747 (2011). [CrossRef] [PubMed]

27. J. Yi, Q. Wei, W. Liu, V. Backman, and H. F. Zhang, “Visible-light optical coherence tomography for retinal oximetry,” Opt. Lett. 38(11), 1796–1798 (2013). [CrossRef] [PubMed]

28. J. Yi, S. Chen, V. Backman, and H. F. Zhang, “In vivo functional microangiography by visible-light optical coherence tomography,” Biomed. Opt. Express 5(10), 3603–3612 (2014). [CrossRef] [PubMed]

29. S. P. Chong, C. W. Merkle, C. Leahy, H. Radhakrishnan, and V. J. Srinivasan, “Quantitative microvascular hemoglobin mapping using visible light spectroscopic Optical Coherence Tomography,” Biomed. Opt. Express 6(4), 1429–1450 (2015). [CrossRef] [PubMed]

30. J. Yi and V. Backman, “Imaging a full set of optical scattering properties of biological tissue by inverse spectroscopic optical coherence tomography,” Opt. Lett. 37(21), 4443–4445 (2012). [CrossRef] [PubMed]

31. J. D. Rogers, A. J. Radosevich, J. Yi, and V. Backman, “Modeling Light Scattering in Tissue as Continuous Random Media Using a Versatile Refractive Index Correlation Function,” IEEE J. Sel. Top. Quantum Electron. 20(2), 173–186 (2013). [PubMed]

32. N. Bosschaart, G. J. Edelman, M. C. Aalders, T. G. van Leeuwen, and D. J. Faber, “A literature review and novel theoretical approach on the optical properties of whole blood,” Lasers Med. Sci. 29(2), 453–479 (2014). [CrossRef] [PubMed]

33. M. Hammer, A. Roggan, D. Schweitzer, and G. Müller, “Optical properties of ocular fundus tissues-an in vitro study using the double-integrating-sphere technique and inverse Monte Carlo simulation,” Phys. Med. Biol. 40(6), 963–978 (1995). [CrossRef] [PubMed]

34. V. J. Srinivasan, J. Y. Jiang, M. A. Yaseen, H. Radhakrishnan, W. Wu, S. Barry, A. E. Cable, and D. A. Boas, “Rapid volumetric angiography of cortical microvasculature with optical coherence tomography,” Opt. Lett. 35(1), 43–45 (2010). [CrossRef] [PubMed]

35. J. Lee, V. Srinivasan, H. Radhakrishnan, and D. A. Boas, “Motion correction for phase-resolved dynamic optical coherence tomography imaging of rodent cerebral cortex,” Opt. Express 19(22), 21258–21270 (2011). [CrossRef] [PubMed]

36. O. Akça, A. G. Doufas, N. Morioka, S. Iscoe, J. Fisher, and D. I. Sessler, “Hypercapnia improves tissue oxygenation,” Anesthesiology 97(4), 801–806 (2002). [CrossRef] [PubMed]

37. G. T. Dorner, G. Garhoefer, C. Zawinka, B. Kiss, and L. Schmetterer, “Response of retinal blood flow to CO2-breathing in humans,” Eur. J. Ophthalmol. 12(6), 459–466 (2002). [PubMed]

38. N. D. Wangsa-Wirawan and R. A. Linsenmeier, “Retinal oxygen: Fundamental and clinical aspects,” Arch. Ophthalmol. 121(4), 547–557 (2003). [CrossRef] [PubMed]

39. A. Alm and A. Bill, “Blood flow and oxygen extraction in the cat uvea at normal and high intraocular pressures,” Acta Physiol. Scand. 80(1), 19–28 (1970). [CrossRef] [PubMed]

40. L. Wang, P. Törnquist, and A. Bill, “Glucose metabolism in pig outer retina in light and darkness,” Acta Physiol. Scand. 160(1), 75–81 (1997). [CrossRef] [PubMed]

	Operations procedure to extract sO₂ for any given blood vessel
1	To generate wavelength-dependent OCT angiography using STFT and motion enhanced dynamic scattering contrast according to Eq. (14)
1		Input: Raw interference spectrum, STFT window; Output: $A G (x, y, z, λ)$
2	To choose a region of interest within a vessel from en face OCT image and average the signal over x, y extents.
2		Input: $A G (x, y, z, λ)$ , ROI of a vessel (x, y); Output: $A A (z, λ)$
3	To calculate localized optical attenuation function $O D (z, λ)$ according to Eq. (4)
3		Input: $A A (z, λ)$ , source spectrum I₀(λ); Output: $O D (z, λ)$
4	To calculate partial derivative on $O D (z, λ)$ against wavelength according to Eq. (5)
4		Input: $O D (z, λ)$ , λ; Output: $\frac{\partial O D (z, λ)}{\partial λ}$ (spectral derivative)
5	To calculate average $\frac{\partial O D (z, λ)}{\partial λ}$ long the z axis and to calibrate all the major retinal arteries against systemic spO₂
5		Input: spO₂, $\frac{\partial O D (λ)}{\partial λ}$ ; Output: calibration curve
6	To convert $\frac{\partial O D (λ)}{\partial λ}$ to sO₂ according to Eq. (11)
6		Input: $\frac{\partial O D (λ)}{\partial λ}$ of selected vessel, calibration curve; Output: absolute sO₂

Measuring oxygen saturation in retinal and choroidal circulations in rats using visible light optical coherence tomography angiography

Abstract

1. Introduction

2. Theory