Retinal oxygen extraction in humans.

Adequate function of the retina is dependent on proper oxygen supply. In humans, the inner retina is oxygenated via the retinal circulation. We present a method to calculate total retinal oxygen extraction based on measurement of total retinal blood flow using dual-beam bidirectional Doppler optical coherence tomography and measurement of oxygen saturation by spectrophotometry. These measurements were done on 8 healthy subjects while breathing ambient room air and 100% oxygen. Total retinal blood flow was 44.3 ± 9.0 μl/min during baseline and decreased to 18.7 ± 4.2 μl/min during 100% oxygen breathing (P < 0.001) resulting in a pronounced decrease in retinal oxygen extraction from 2.33 ± 0.51 μl(O2)/min to 0.88 ± 0.14 μl(O2)/min during breathing of 100% oxygen. The method presented in this paper may have significant potential to study oxygen metabolism in hypoxic retinal diseases such as diabetic retinopathy.

In 2010, a total of 185 million people were legally blind worldwide1. In industrialized countries, the three major causes of blindness are age-related macular degeneration, glaucoma and diabetic retinopathy. Alterations in ocular blood flow and oxygenation have been implicated in all of these conditions23.

The human retina has a complex system of vascular supply to ensure adequate oxygenation. The inner retina including the retinal ganglion cells (RGC) is nourished by the retinal circulation whereas the outer retina is supplied by the choroidal vasculature4. As in the brain, the retina has a continuous demand for oxygen and, as such, a lack of oxygen delivery will be associated with hypoxia leading to retinal disease56.

Most of the knowledge about retinal oxygenation arises from animal studies using microelectrodes789. This technique allows for the mapping of oxygen gradients with high spatial and temporal resolution and retinal oxygen profiles in pigs10, cats11, rats12 and monkeys1314 were published based on this method. Except for the monkey fovea, oxygen profiles in these species share a lot of similarities. Under light-adapted conditions, they are characterized by a steep fall between the choriocapillaris and the photoreceptor inner segments5. Subsequently, oxygen tension decreases more gradually through the outer retina and has a value of approximately 20 mmHg in the inner retina5.

In humans, measurement of oxygen saturation with microelectrodes is not applicable due to the invasive nature of this technique. Here we describe a new method for the measurement of total retinal oxygen extraction in healthy subjects. This was done by combining measurements of total retinal blood flow using dual-beam bidirectional Doppler Fourier domain optical coherence tomography15 (FD-OCT) with measurements of oxygen saturation in retinal vessels using fundus reflectometry1617. A mathematical model was developed to translate these measurements into total retinal oxygen extraction. Measurements were done during both breathing of ambient room air and during 100% oxygen breathing. During systemic hyperoxia, there is evidence that the amount of oxygen extracted from the retinal circulation shows a pronounced decrease, and we believe that this occurs because oxygen diffuses from the choroid to the inner retina181920.

Research Design and Methods

Subjects

The study protocol was approved by the Ethics Committee of the Medical University of Vienna and followed the guidelines set forth in the Declaration of Helsinki. A total of eight healthy subjects (six male and two female) between 18 and 35 years were included in the study after signing written informed consent. All subjects passed a screening examination before the study day that included physical examination, blood draw to assess hematological status and chemistry, 12-lead electro cardiogram, assessment of visual acuity, slit lamp biomicroscopy, funduscopy and measurement of intraocular pressure (IOP). Exclusion criteria were ametropia ≥3 diopter, anisometropia ≥3 diopter, other ocular abnormalities and any clinically relevant illness as judged by the investigators, blood donation and/or intake of any medication in the three weeks prior to the study. Participants had to abstain from beverages containing alcohol or caffeine in the 12 hours before the study day.

Protocol

After instillation of one drop of tropicamide (Mydriatikum AGEPHA, Vienna, Austria) into the study eye, a resting period of at least 20 minutes was scheduled. Thereafter, retinal blood velocities were measured using dual-beam bidirectional Doppler FD-OCT. Fundus photographs were taken using the Retinal Vessel Analyser (RVA, Imedos Systems UG, Jena, Germany) in order to quantify vessel diameters and oxygen saturation. A capillary blood sample was obtained from the earlobe for blood gas analysis. Blood pressure and pulse rate were measured non-invasively. Oxygen (gases for human use, Messer; Vienna, Austria) was provided via a partially expanded reservoir bag at atmospheric pressure. For gas delivery to the subject, a face mask which covered mouth and nose connected to a two-valve system, preventing the subject from rebreathing was used. After a resting period of 20 minutes, the phase of inhaling 100% oxygen that lasted 30 minutes began. Retinal hemodynamic measurements were performed 15 minutes after the start of inhalation.

Measurement of retinal blood flow and oxygen saturation

The measurements were performed with a dual-beam Doppler FD-OCT system, described previously152122. Two orthogonally polarized probe beams in the interferometer’s sample arm are focused onto the retina under a known angle Δα between them. The system records the spectra of the two channels as a function of frequency using two identical spectrometers. When the probe beams fall onto moving red blood cells (RBC), the back reflected light in each channel is Doppler-shifted. These Doppler shifts Φ1 and Φ2 can be obtained by calculating the difference of phase of the Fourier transform of subsequent A-line recordings at the spectrometer’s CCD cameras which are operated at a line scan rate of τ =27 μs. Then, the phase difference ΔΦ=Φ1 − Φ2 between the two beams and system channels is calculated and the absolute flow velocity can be obtained as

In equation (1), λ0=838.8 nm is the light source’s central wavelength, n = 1.37 is the refractive index of blood and β is the angle of the vessel’s blood flow velocity vector with respect to the plane spanned by the two probe beams. The angle Δα was calculated for each subject’s eye individually based on the separation of the two probe beams at the pupil plane, the eye length and the ametropia value. All phase images were corrected for bulk-motion and phase wrapping23.

The OCT-setup was integrated into the RVA, which consists of a modified fundus camera (FF450plus; Carl Zeiss Meditec AG, Jena, Germany) and a CCD camera allowing for the measurement of retinal vessel diameters and oxygen saturation. The OCT beams are coupled into the RVA using a dichroic mirror, which ensures independent operation of the OCT system and the fundus camera. Using a rotatable beam displacer, the detection plane spanned by the two probe beams can either be aligned vertically or horizontally relative to the eye. This is required because without rotation of the detection plane, vessels with an orientation perpendicular to this plane would not be assessable, since the angle β in equation (1) would become zero. To determine the total retinal blood flow, OCT measurements are performed in a rectangular scanning pattern around the optic nerve head (ONH) (Fig. 1), which ensures that the mean velocities (vA,i, vV,j) in all the blood vessels leaving or entering the ONH can be measured. The measurement time per single location is approximately 5 s, which allows for averaging over several pulse periods. The angle β was determined from the fundus photographs recorded with the RVA’s camera.

Figure 1

Sample measurement in a healthy subject.

(a) Fundus image and OCT scanning patterns used to measure all retinal vessels entering the optic nerve head to measure total retinal blood flow. Black bars indicate the measurement locations. (b,c) Magnification of the fundus image during breathing room air and 100% oxygen, respectively. Vasoconstriction of retinal artery and vein can clearly be seen. (d,e) Time course of retinal arterial and venous blood velocity during breathing room air and breathing 100% oxygen respectively over approximately 4 cardiac cycles.

Concomitantly with the OCT measurements, the diameters of the measured vessels (dA,i, dV,j) were determined with the RVA using static vessel analysis. The analyzing software was modified to allow for calculation of absolute values in μm based on the individual axial eye length, the axial refraction for each subject and the optical characteristics of the fundus camera system. Since the RVA measurements are less reliable for very small vessels, vessel diameters below 65 μm were measured in the OCT images and determined by comparison with the RVA results from larger vessels as described in detail previously15.

Using the mean velocity v from the OCT recordings and the vessel diameter d from the RVA measurements, blood flow Q in each individual vessel entering or leaving the ONH can be calculated as

A is the cross sectional area of the vessel and is calculated assuming a circular cross section:

Due to the equation of continuity and because the retina is an end organ, the blood flow in the central vessels can be calculated based on these measurements. By summing up the blood flow values QA,i in all measured arteries, one can determine the blood flow QA,tot in the central retinal artery, and by summing up the blood flow values QV,j in all measured veins, it is possible to calculate the blood flow QV,tot in the central retinal vein. In formula (4) and formula (5), #A is the total number of measured arteries and #V is the total number of measured veins.

Obviously, QA,tot needs to be equal to QV,tot, which provides a validity check for each measurement. Total retinal blood flow was calculated as the arithmetic mean of QA,tot and QV,tot:

The fundus photographs were used to extract the oxygen saturation in retinal arteries and veins1624. This is based on spectral analysis of light at selected wavelengths, which is reflected at the fundus. In the present system, two fundus pictures at wavelengths of 610 and 545 nm, respectively, are taken. Extraction of the retinal oxygen saturation is based on the fact that oxygenated hemoglobin has different light absorption characteristics as compared to deoxygenated hemoglobin. The isosbestic wavelength is 548 nm; it is defined as the point in the light spectrum where oxygenated and deoxygenated hemoglobin show identical absorption. On the other hand, oxygenated hemoglobin is nearly transparent if it is illuminated with light at a wavelength of 610 nm. Quantifying the contrast at these wavelengths enables the determination of the relation between oxygenated and total hemoglobin and the calculation of the oxygen saturation. The oxygen saturation was measured in all retinal arteries (SaO2,A) and all retinal veins (SaO).

Model to calculate total retinal oxygen extraction

Measurements as described in the previous section were used to calculate the total retinal oxygen extraction. The visible arterial branch arteries at the posterior pole of the eye arise from the central retinal artery, which bifurcates at the center of the optic disk. The visible retinal branch veins merge into the central retinal vein. When measuring at some distance from the ONH, it is therefore necessary to correct the oxygen saturation readings for the oxygen loss across the arterial or venous wall between the point of measurement and the point where the vessels merge. The oxygen loss is described in formula (7) and depends on the diameter and the length of the vessel25:

where Q is blood flow, Hb is the hemoglobin concentration, Sin − Sout is the difference in oxygen saturation in the specific vessel between the bifurcation point at the disk and the point where saturation is measured, R is the radius of the vessel, L is the distance between the point of measurement and the center of the ONH and is the oxygen loss through the vascular wall. The value 1.35 is the amount of oxygen in [ml] which can bind to 1 g hemoglobin (Fig. 2).

Figure 2

% change in total retinal blood flow (Q), arterial oxygen content (cO2,CRA), venous oxygen content (cO2,CRV), arteriovenous difference in oxygen content (cO2,DIFF) and retinal oxygen extraction (extO2) in response to 100% oxygen breathing (n = 8).

Data are presented as means ± SD.

Obviously, the oxygen loss across the wall depends on how many oxygen molecules diffuse into the vitreous or the retinal tissue. This value is described as the oxygen flux . Measurements in the illuminated cat retina have shown that is . In the present model, this value was used for the calculation. To calculate the corrected oxygen saturation (cSaO2) when a vein leaves the retina, the estimated difference (Sin − Sout) in the oxygen saturation has to be subtracted from the measured value (SaO2). In the arteries, on the other hand, the saturation at a point L from the ONH will be lower than the actual saturation of the central retinal artery, so the corrected value when the vessel joins the central retinal artery is higher than the measured value. Thus,

The value cSaO2 is used for the further calculations. The corrected oxygen saturations for all arteries in the posterior pole of the eye are used to calculate the mean value of the oxygen saturation in the central retinal artery (SaO2,CRA):

where #A is the total number of measured arteries.

The calculation of the oxygen saturation in the central retinal vein is more difficult, as the blood in the central retinal vein is a mixture of all veins in the back of the eye, which have different oxygen saturations. The corrected venous values in each vein also need to be weighted by the blood flow in each vein in order to calculate the oxygen saturation in the central retinal vein. The corrected oxygen saturation of each vein (SaO2,V) must be multiplied with the blood flow QV,j in the specific vessel and divided by the total flow in all veins QV,tot to obtain the central retinal vein’s oxygen saturation:

where #V is the number of measured retinal veins.

To calculate the retinal oxygen extraction of the human retina, we need to evaluate the amount of oxygen in the retinal vessels. The amount of oxygen includes two components: the amount bound to hemoglobin, and the amount dissolved in plasma. The amount of oxygen which is bound to hemoglobin can be easily calculated using the oxygen saturation, the hemoglobin concentration in the blood and the factor 1.35, as 1.35 ml oxygen is bound to 1 g hemoglobin. The fraction of dissolved oxygen can only be estimated, since it is not possible to measure the oxygen partial pressure in retinal vessels in humans. The oxygen partial pressure has to be multiplied by 0.003 ml O2/mmHg26. Under physiological conditions, the amount of oxygen that is dissolved is low. During 100% oxygen breathing, dissolved oxygen becomes, however, a significant proportion of the bloods oxygen content27.

To evaluate the arterial oxygen partial pressure, we used an arterialized blood sample from the earlobe. Venous PO2 was estimated from the saturation values by using the oxygen-binding curve28 (P50: 26.8 mmHg) at a PCO2 of 37 mmHg and a temperature of 37 degrees. Using these corrections, the oxygen content in the retinal vessels can finally be calculated as follows:

In this equation, Hb is the hemoglobin concentration in g per ml, is the oxygen saturation and is the oxygen partial pressure; ves stands for the type of vessel (artery or vein).

The oxygen content of a vessel indicates how much oxygen is transported in a certain period of time. The oxygen extraction is the amount of oxygen which is consumed by the retinal tissues; it is calculated by subtracting the venous from the arterial oxygen content and multiplying by the blood flow, as given in formula (13).

This formula was used to estimate total oxygen extraction from the retinal circulation. To obtain data per 100 g tissue, a value of 326 mg was chosen as total weight of the human retina29.

Measurement of IOP and systemic hemodynamics

Systolic, diastolic and mean arterial blood pressures (SBP, DBP, MAP) were measured on the upper arm by an automated oscillometric device (Infinity Delta, Dräger, Vienna, Austria). The same device was used to record pulse rate and systemic oxygen saturation by a finger pulse oximeter.

Blood gas analysis

Arterialized capillary blood from the earlobe was collected from a lancet incision into a thin glass capillary tube. Arterial pH, PCO2, and PO2 were determined using an automatic blood gas analysis system (AVL 912; CO-Oxilite, Graz, Austriax). Hemoglobin concentration (Hb) was measured photometrically (Sysmex XE 500, Kobe, Japan).

Statistical analysis

Data are presented as means ±SD. Paired t-tests were used to compare baseline values to values obtained during 100% oxygen breathing. A p-value < 0.05 was considered the level of significance.

Results

Baseline IOP in healthy subjects was 16.8 ± 1.4 mmHg. Breathing 100% oxygen did not induce a significant effect on either blood pressure or HR (Table 1). Retinal blood flow under room air conditions was 44.3 ± 9.0 μl/min. On average, 11.0 ± 2.2 arteries and 8.1 ± 1.6 veins were measured. There was no significant difference between QA,tot (44.2 ± 10.6 μl/min) and QV,tot (44.4 ± 7.6 μl/min), the absolute value of the mean difference between the 2 values was 8.0 ± 4.5%. During air breathing, the mean value of was 95.3 ± 1.9%, and the mean value of was 68.0 ± 3.9%. With corrected oxygen content values of  = 0.181 ± 0.008  and  = 0.128 ± 0.011 , this resulted in a total retinal oxygen extraction () of 2.33 ± 0.51 μl(O2)/min. Analysis of both histological30 and OCT31 data for total and inner retinal thickness and numerical integration of these gave values of 317 mg and 137 mg for the total retinal mass and the mass of the inner retina, respectively (unpublished data). The weight of the retina was found to be in good agreement with the literature value of 326 mg29. This gives an inner retinal fraction of 0.43. However, the retinal circulation also supplies about 11 percent of the outer retina in macaque14, yielding that it provides blood to about half of the retinal weight. Taking these values into account, the result for retinal oxygen extraction can be written as 1.42 ± 0.31 ml(O2)/min/100g tissue.

Table 1

Blood pressure and pulse rate at baseline and during 100% oxygen breathing (n = 8). Data are presented as means ± SD.

	Baseline	100% oxygen breathing	p-value
Systolic blood pressure (mmHg)	119.5 ± 6.3	121.6 ± 7.1	0.48
Diastolic blood pressure (mmHg)	69.8 ± 11.9	72.8 ± 7.7	0.41
Mean blood pressure (mmHg)	87.9 ± 7.4	91.8 ± 5.0	0.19
Heart rate (beats/min)	69.1 ± 16.4	66.1 ± 9.8	0.33

The percent changes of retinal parameters during oxygen breathing are depicted in Fig. 2. During 100% oxygen breathing, we observed a pronounced decrease in retinal blood flow to a value of 18.7 ± 4.2 μl/min (p < 0.001). Again, there was no difference between QA,tot (19.0 ± 4.1 μl/min) and QV,tot (18.4 ± 4.5 μl/min), and the average absolute value of the difference was small (9.2 ± 3.8%). Both and increased to values of 99.4 ± 0.3% (p = 0.004) and 76.6 ± 5.2% (p = 0.015), respectively. Hence, corrected oxygen saturation values during 100% oxygen breathing also increased ( = 0.194 ± 0.010 , p < 0.001 and  = 0.137 ± 0.013  , p = 0.012). Since both retinal blood flow and the arterio-venous difference in corrected oxygen content decreased, we also observed a pronounced decrease in retinal oxygen extraction during 100% oxygen breathing ( = 0.88 ± 0.14 μl(O2)/min or 0.54 ± 0.28 ml(O2)/min/100g tissue).

Table 2 shows data obtained in one healthy subject in all measured blood vessels. Figure 1 shows a sample measurement during room air breathing (b,d) and breathing 100% oxygen (c,e).

Table 2

Values as obtained in one healthy subject in all retinal vessels. (I = inferior, N = nasal, S = superior, T = temporal).

Vessel	Artery(A)/ Vein(V)	Velocity (mm/s)	Diameter (μm)	Flow (μl/min)	Oxygen saturation (%)	Oxygen content (ml(O2)/ml)
I2	A	19.48	111.6	10.23	93.7	0.941
I3	A	7.53	70.2	1.56	93.2	0.950
I5	A	10.26	79.1	2.71	93.5	0.946
N1	A	7.94	80.8	2.19	92.4	0.938
N3	A	13.45	79.7	3.60	91.9	0.928
S2	A	7.78	65.6	1.41	77.8	0.796
S4	A	10.70	91.1	3.74	91.2	0.922
T2	A	10.22	82.1	3.14	94.3	0.953
T5	A	9.38	83.3	2.74	92.6	0.938
T7	A	15.79	109.1	7.93	91.5	0.920
I1	V	12.83	135.7	9.96	62.7	0.622
I4	V	5.34	93.5	1.97	67.2	0.653
I6	V	5.69	107.2	2.76	55.5	0.540
N2	V	6.31	108.6	3.14	59	0.576
N4	V	7.02	109.5	3.55	71.44	0.702
S1	V	8.40	165.5	9.71	67.8	0.671
S3	V	4.64	103.1	2.08	65.6	0.637
T1	V	2.72	80.7	0.75	57.1	0.529
T3	V	1.96	55	0.88	63.4	0.610
T4	V	7.87	71.1	1.68	62.4	0.607
T6	V	7.53	95.9	2.92	59.2	0.579

Discussion

We present a method that is capable of non-invasively measuring total retinal oxygen extraction in humans. The value of 2.33 μl(O2)/min or 1.42 ml(O2)/min/100 g tissue blood is considerably lower than the value of 8 ml(O2)/min/100g tissue estimated previously32. The main difference is that the authors of this previous study estimated a flow rate of retinal blood as high as 170 ml/100 g tissue/min, which is much higher than the values obtained in the present study. Our values for retinal blood flow are, however, in the same range as those obtained by other authors using a variety of different methods1533343536373839.

Our data are also lower than those obtained from animal experiments, but the difference is less pronounced. Under photopic conditions, several authors reported oxygen consumption values for the inner retina using different techniques. Employing Fick’s principle, a value of 3.9 ml(O2)/min/100 g tissue blood was reported in the pig retina40. Lower values of 2.3 ml(O2)/min/100 g tissue blood were reported in the rat by using the Warburg technique, an in vitro application of Fick’s law41. Values were also reported for inner retinal oxygen consumption in the cat42. By modeling data from microelectrode experiments during occlusion of the central retinal artery, the authors reported 3.9 ml(O2)/min/100 g tissue blood. It needs, however, to be considered that data from experiments using microelectrodes represent local values from the center of the retina. By contrast, our data as well as the data obtained via Fick’s method represent total retinal oxygen extraction including the peripheral retina, which contains less RGCs and less microvessels443.

Another study reported values for inner retinal oxygen extraction in the rat of 516 nl(O2)/min44. In this study, oxygen tension was measured by phosphorescence lifetime imaging and total retinal blood flow was measured by red-free and fluorescent microsphere imaging. In another recent study in rats, the inner retinal oxygen extraction was reported to be 300 nl(O2)/min45. In this study, total retinal blood flow was measured via Doppler FD-OCT, which is feasible as in anesthetized animals it is, due to the absence of sample motion, not required to use a dual-beam system that is insensitive to changes in the Doppler angle. Photoacoustic ophthalmoscopy was used to measure oxygen saturation. Using literature data for the weight of the rat’s retina of 14 mg46 and considering that 50% of oxygen consumption takes place in the inner retina this can be converted to 7.4 ml(O2)/min/100 g tissue (39) and 4.3 ml(O2)/min/100 g tissue (40) for the two studies, respectively.

Our data indicate that during 100% oxygen breathing there is a pronounced decrease in oxygen extraction from the retinal circulation. This is caused by a pronounced decrease in retinal blood flow as well as a decrease in arterio-venous oxygen difference. Our results are in good agreement with our previous data in which oxygen extraction was estimated from measurements on single arteries and veins20. In animal studies, it has been shown that inner retinal oxygen tension is well regulated during 100% oxygen breathing18194748. In these experiments, delivery of oxygen by the retinal circulation shows a pronounced decrease due to a pronounced decrease in retinal blood flow. On the other hand, choroidal blood flow shows almost no reaction to systemic hyperoxia49, which leads to significant diffusion of oxygen from the choroid to the inner retina.

The results of the present study are critically dependent on the validity of the underlying techniques. Measurement of retinal vessel diameters based on fundus photography is the gold standard technique providing excellent reproducibility and sensitivity50. The validity of dual-beam bidirectional Doppler FD-OCT is shown by several previous experiments. There is a high degree of association with blood velocities as measured with laser Doppler velocimetry, a technique widely used in clinical research, the law of mass conservation is well fulfilled at retinal bifurcations, and Murray’s law is fulfilled in our measurements152223. In addition, values of total retinal arterial flow and total retinal venous flow are in good agreement in this study as it was the case in our previous experiments15, as can be expected for an end organ.

Measurements of oxygen saturation in retinal vessels were based on a two wavelength approach of fundus reflectometry. This technique, which is based on the different absorption characteristics of oxy- and deoxy-hemoglobin has several limitations. When photons enter the camera, the optical pathway is generally unknown. Hence, not all photons arising from the locations of the blood vessels have necessarily travelled through the entire vessel due to scattering. In addition, scattered photons from the areas adjacent to the vessel may have undergone absorption in the melanin of the retinal pigment epithelium. As such it does not come as a surprise that Monte Carlo simulation indicates an influence of vessel diameter and melanin on the accuracy of oxygen saturation measurements51. In the present system, this is overcome by compensating for these parameters based on a large dataset of normal values16. Validity of data as obtained with two wavelength systems as employed in the present study is strongly dependent on this calibration process. With the device produced by Imedos this calibration is based on retinal vessel reflectance spectra with a 2-nm resolution52. Nevertheless other investigators reported lower values for SaO2,V in healthy volunteers with a mean value of 60%53. This has an effect on our calculation of oxygen extraction based on equation (13) and could lead to a slight underestimation of values with our technology. On the other hand, experimental evidence indicates that reduced arterial oxygen values in patients with Eisenmenger’s syndrome54, chronic obstructive pulmonary disease55 and breathing gas mixtures with reduced fraction of inspired oxygen56 are well reflected in the retinal oxygen saturation measurements. Hence, the difference between values during normoxia and hyperoxia are not affected by this limitation.

In conclusion, we present a system that is capable of measuring retinal oxygen extraction in humans. This technique may have considerable potential for diagnosis, and risk stratification treatment monitoring in patients with retinal vascular disease.

Additional Information

How to cite this article: Werkmeister, R. M. et al. Retinal oxygen extraction in humans. Sci. Rep. 5, 15763; doi: 10.1038/srep15763 (2015).