|
1.IntroductionMicrocirculation is the blood flow through arterioles, capillaries, and venules, which exists throughout all tissues in a living body. Microcirculation plays an important role in supplying nutrients and oxygen to cells in the tissues. Thus, insufficient supplies of oxygen in microcirculation may lead to tissue hypoxia. Tissue hypoxia is considered to be one of the factors of serious diseases, such as cystic fibrosis and chronic bronchitis.1 Moreover, it is well known that varieties of experimental shock states (e.g., sepsis and haemorrhage) are related to hypoxia.2–4 In the early stages of these shock states, the impaired microcirculation will result in the rapid onset of tissue hypoxia. However, mechanisms of these diseases due to the hypoxia condition have not been clarified. Visualizing an individual microcirculation and observing the changes in the oxygen saturation () during disease progression can be effective for understanding disease dynamics. Many imaging techniques for in microcirculation have been developed. Snapshot multispectral imaging has been used for microcirculation oximetry mainly in the retina by using the difference in optical absorption properties of oxyhemoglobin () and deoxyhemoglobin (HbR).5 Photoacoustic imaging (PAI)6 and optical coherence tomography (OCT)7 also have been introduced to measure in microcirculation. However, these techniques require complex and expensive optical hardware. Therefore, simpler and lower cost optical setups would promote easier microcirculation assessment. To visualize tissue microcirculation, orthogonal polarization spectral (OPS) imaging was developed as a real-time and noninvasive imaging technique of handheld vital microscopes (HVMs).8 This imaging technique was used to reveal alterations of the microcirculation in patients with sepsis.9 The second-generation HVM is based on sidestream dark-field (SDF) imaging.10,11 This imaging technique can acquire higher contrast microcirculation images than OPS imaging. From the SDF images of individual vessels, some physical quantities such as vessel length, vessel diameter, and velocity of red blood cells (RBCs) can be estimated.12 The first commercial SDF device (Microscan, Microvision Medical B.V., Amsterdam, The Netherlands) has been developed. To verify its effectiveness, evaluation experiments on patients with diabetes were performed.11,13 The results showed that the microcirculation density of patients with diabetes was significantly higher than that of healthy people. Compared with other imaging techniques of microcirculation such as PAI and OCT, SDF imaging is a simple and low cost technique. Hence, SDF imaging is easily applicable to clinical purposes and enables users to observe the microcirculation of human tissues. Based on SDF imaging, we have developed a trial SDF probe and spectral image-based oximetry method, which we called SDF oximetry.14,15 Our SDF oximetry utilizes two-wavelength oximetry based on the Beer–Lambert law. In another study, we applied average extinction coefficients (AECs)16 of hemoglobin to SDF oximetry and improved the accuracy of our method. However, in our SDF oximetry, there are three problems that have to be overcome. First, we have already confirmed the effectiveness of SDF oximetry in tissue-like turbid phantom experiments; however, we have not carried out the in vivo experiments. Second, since we calculate the AECs from the SDF images of the turbid phantom experiments, they are not optimized for in vivo SDF images. Third, heat from the light-emitting diodes (LEDs) used for the SDF imaging might affect hemodynamics in microcirculation since the SDF probe is in direct contact with the tissue. We need to investigate the possible effect and that is the focus of this paper. In this study, we conducted hypoxic stimulation experiments with rats by changing fraction of inspired oxygen (). As a result, we confirmed that the estimated by SDF oximetry during the hypoxic stimulation experiment corresponded to the changes in . Prior to performing the in vivo experiment, we recalculated AECs from the literature values. Then, we examined the heat influence of the LEDs. Next, we conducted the experiment with rats where we observed the time-lapsed changes in under two conditions: normoxia and hypoxia. Then, we compared these values to evaluate our method, and the results suggested that SDF oximetry could response to alterations in microcirculation correctly. 2.Materials and Methods2.1.SDF ImagingFigure 1 shows the overview of our trial SDF imaging system. Our system comprises a light source unit and SDF probe, as shown in Fig. 1(a). Figure 1(b) is a schematic illustration of an image acquisition. The SDF probe is in direct contact with the tissue and the pulsed illumination from the LEDs placed concentrically around the probe directly penetrates into the tissue. SDF imaging can only capture the light diffused in the tissues without surface reflections because the camera inside the SDF probe is optically isolated from the peripheral illumination. In SDF images, blood vessels are represented in black due to optical absorption by hemoglobin molecules of the RBCs. As an imaging device, a complementary metal oxide semiconductor (CMOS) camera (ID04MB-IP-U, iDule Corporation, Chiba, Japan) with a customized magnifier (focal length: 5.6 mm; -number: F4.0) is used. Figure 2(a) shows the two-color LEDs (SMLVN6RGB1W, Rohm Co., Ltd., Kyoto, Japan) used in the system. The spectral intensities of the LEDs are shown in Fig. 2(b) together with extinction coefficients of and HbR as a function of wavelength.17 The peak wavelengths of the two-color LEDs are 470 nm (blue) and 527 nm (green) and their full widths at half maximum are 27.3 and 37.2 nm, respectively. The wavelength band of the blue LED is sensitive to changes in hemoglobin oxygenation, whereas that of the green LED is insensitive. Therefore, we can acquire microcirculation images in different contrasts corresponding to hemoglobin oxygenation. 2.2.SDF OximetrySDF oximetry utilizes two-wavelength oximetry based on the Beer–Lambert law.14 Figure 3 shows the assumption of and in the SDF imaging. We define the backreflection light entering into the blood vessel as the incident light . When enters into the blood vessel with thickness , the Beer–Lambert law gives the transmitted light through the medium as where represents the molar extinction coefficient of the hemoglobin and represents the total hemoglobin concentration. Hemoglobin has several variants, namely methemoglobin (MetHb), carboxyhemoglobin (CoHb), and sulfhemoglobin (SulfHb). The absorption coefficients of these three hemoglobins have peaks in the range of visible light.18,19 However, concentrations of MetHb, CoHb, and SulfHb included in blood are normally lower than and HbR, thus we can ignore them in our oximetry method. Moreover, absorption and scattering of leukocytes and platelets in blood occurs.20 However, these particles account for approximately only 4% of the whole blood volume, and hence they do not contribute more to absorption or scattering than RBCs do. As a result, for blood, is approximately defined as where the range is 0 to 1. The optical density of blood is written asFinally, solving two equations of Eq. (3) for and , we obtain as follows: where for . Figure 4 shows the definition of and in SDF images. In our definition, the intensity of vascular regions is treated as and that of avascular regions treated as . Then, we can calculate by Eq. (4).The value calculated by Eq. (4) cannot ignore the influence of the bandwidth of the LEDs shown in Fig. 2(b). We have to obtain the absorbance coefficients by wavelength averaging because the spectral distribution of the LEDs is not truly monochromatic. To take into account the bandwidth of the LEDs, we apply AECs and instead of .14 When the full width at 10th maximum is , the intensity of light reaching the detector from vascular regions of and HbR can be expressed as where represents the spectral sensitivity of the camera and represents the spectral distribution of the ’th color of the LEDs. The intensity of light reaching the detector from avascular regions is expressed asHence, the AECs of and HbR can be defined as We calculated these AECs from the spectral distribution of the extinction coefficients of and HbR,14 the spectral intensity distribution of the two-color LEDs, and the spectral sensitivity characteristics of the CMOS camera, shown in Figs. 2(b) and 2(c). The calculated results were , , , and . 2.3.In Vivo Imaging PreparationFigure 5 is a photo of the experiment setup. In our experiment, we observed dermis microcirculation of male Slc:Wistar rats (260 to 280 g, 12-weeks old, ). One was used for an experiment of stability evaluation, and the other two were for a hypoxic stimulation experiment. These rats were anesthetized under 2.0% isoflurane (Escain, Pfizer, Tokyo, Japan). A heat mat (Kainuma Industrial Co., Ltd., Aichi, Japan) was laid under the rat to help maintain its body temperature. To obtain clear images, the epidermis was removed before the observation. To monitor the percutaneous arterial oxygen saturation (), a rat pulse oximeter (MouseOx Plus, STARR Life, Pennsylvania) was used. As the first experiment, to ensure the time-lapsed changes in of hypoxic stimulation are not influenced by the heat from LEDs, a rat (rat 1) was observed during an 11-min period under normoxia. The rat’s mouth and nose were covered with a mask that was connected to an inhalation tube. The SDF images were acquired under blue and green illuminations alternately every 100 ms. Next, to assess the sensitivity of the SDF oximetry in accordance with physiological changes, 2 rats (rat 2 and 3) were exposed to hypoxic stimulation by changing . Similar procedures have been employed in previous studies.21,22 They were kept to inhale the air mixed with isoflurane for the first 3 min. The oxygen concentration was then reduced to 10% oxygen by increasing nitrogen concentration together with isoflurane for 3 min to get the hypoxic stimulus state, which was followed by a switch back to the normal state for 5 min. The SDF images were acquired in the same manner as done in the first experiment. All experiments in this study were carried out in conformity with the Institutional Animal Care and Use policy of Chiba University. 3.Experiments and Discussion3.1.Time-Lapsed Stability Evaluation of SDF OximetryThe hypoxic stimulation experiment, described in Sec. 3.2, lasted 11 min. Due to the long measurement time, there are two possible influences on the estimated . The first is the heat influence caused by the LEDs. In SDF imaging, the SDF probe in which the LEDs are arranged is in direct contact with the tissue. Hence, the self-heating of LEDs caused by continuous light emission may affect microcirculation hemodynamics. The second is the influence of isoflurane. For example, Helmchen et al.23 observed anesthetic-induced changes in microcirculation hemodynamics. Therefore, we have to evaluate the stability of SDF oximetry during an 11-min period to ensure that time-lapsed changes in of hypoxic stimulation are not influenced by these two factors. Prior to performing the in vivo experiment, we measured the time-lapsed changes in the temperature of the surface of SDF probe. We used a thermometer (AM-8051E, Anritsu Meter Co., Ltd., Tokyo, Japan) with a flexible internal temperature sensor (SF-E-200-ANP, Anritsu Meter Co., Ltd., Tokyo, Japan) for temperature measurement. As shown in Fig. 6, we found that the temperature of the SDF probe surface gradually increased and reached . Then, we performed the in vivo experiment to evaluate the stability of SDF oximetry during an 11-min period. Figure 7(a) shows the SDF image obtained under the blue LED illumination at 0 min. We analyzed the time-lapsed changes in calculated by Eq. (4) for each region of interest (ROI) and calculated mean and standard deviation as shown in Fig. 7(b) and Table 1. The figure and table showed that there were no significant changes in during 11 min. The results indicated the time-lapsed stability of our SDF oximetry. It means that in an 11 min-long experiment, there is no influence of self-heating of the LEDs and isoflurane. In Fig. 7(b), large variations in the estimated value of are observed. It is because that the two SDF images were acquired under blue and green illuminations alternately every 100 ms. Since these images have a time-lag, each of their corresponding pixels used for estimation does not contain the exact same density of RBCs. If the density of RBCs within the pixel differs between two acquired SDF images, it results in an over- or underestimation of . However, adopting a moving average, we can reduce the estimation error. Table 1Mean ± standard deviation (%) of estimated SO2 (rat 1).
3.2.Observation of Time-Lapsed Changes in During Hypoxic StimulationThe hypoxic stimulation experiment in this section was repeated using two rats. We describe results for rat 2 first. As shown in Fig. 8, the response of the measured by the pulse oximeter during the hypoxic stimulation experiment corresponded to the change in . The decrease and increase in started after the change in . Figure 9 shows the time-lapsed SDF images. We could observe a slight vasoconstriction and a decrease in blood flow in accordance with the hypoxic stimulation. The vasoconstriction is caused by temporary elevation of blood pressure associated with hypoxia.24 After hypoxia, these two characteristics returned to their original state. In the analysis, we focused on two blood vessels. One had a bifurcation in the field of view (blue arrow) and the other had no bifurcation (red arrow). In particular, we predicted that the blood vessel indicated by the blue arrow was the venule because the blood flow was merged at the vessel branch. of individual vessels were calculated for each pixel from two-band SDF images. Figure 10 shows the changes in maps during hypoxic stimulation. The results clearly showed that the in microcirculation decreased during hypoxic stimulation and increased after hypoxia. Figure 11(b) shows the changes in estimated during hypoxic stimulation corresponding to the rectangular area shown in Fig. 11(a). The size of the area was set to . The gray-shaded area meant the period of hypoxic stimulation. These results mostly corresponded to the change in the , measured by the pulse oximeter shown in Fig. 8. After returning to the normoxia stimulation, the estimated increased nearly 100% almost instantaneously and then decreased gradually. This result might reflect the function of microcirculation to regulate the tissue oxygenation. Further, we focused on the difference in response to hypoxia between two blood vessels, which are marked by the red and blue arrows in Fig. 11. As shown in Fig. 11(b), the of the blue rectangular area began to change before there were observable changes in the red rectangular area. These differences between the two blood vessels may be due to the difference between arterioles and venules as mentioned in a previous study.5 In SDF imaging, arterioles and venules are distinguished by the blood flow direction at the vessel branch. However, more accurate identification requires a larger image size than that of the present SDF images (). In order to ensure the reproducibility of the hypoxic stimulation experiment, we repeated it with the other rat (rat 3). The experimental conditions were the same as the first time. Figure 12 shows that the created maps had a tendency similar to that of Fig. 10. In addition, the difference in the change of timing between different vessels was also observed as shown in Fig. 13. In the repeated experiment, the response speed for the estimated to hypoxia was faster than that of the first. We considered that to be due to individual differences between the rats or to differences of the blood vessel diameters. Figure 14 shows SDF images obtained at 0 and 7 min. There were two characteristics. First, the blood vessel indicated by the red arrow decreased in diameter compared with that indicated by the blue arrow. To compensate for the insufficient oxygen supply, the diameter of the arterioles presumably varied for the increasing blood pressure. Second, blood vessel density in the areas indicated by the yellow circles was reduced by hypoxia. We considered that the decrease in blood flow caused by hypoxia resulted in the lack of perfusion to capillaries and the reduction of blood vessel density. Our results showed that SDF oximetry could estimate of microcirculation and respond to an alteration in microcirculation in hypoxic stimulation experiments with rats. This indicates that SDF oximetry can be applied for understanding microcirculatory dysfunction in oxygen delivery induced by several diseases, as mentioned in the introduction. However, the accuracy of absolute value estimated by SDF oximetry is uncertain. Therefore, we have to compare the value estimated by SDF oximetry with the one obtained by using high accuracy oximeter (e.g., blood gas analyzer). 4.ConclusionIn this paper, we have applied our SDF oximetry to in vivo experiments and observed the time-lapsed SO2 changes in dermis microcirculation of rats under hypoxic stimulation. Before performing the in vivo experiment, we calculated the AECs from the literature values because the AECs used in our previous study were optimized only for turbid phantom experiments. Next, to confirm the time-lapsed stability of SDF oximetry, we obtained SDF images under the normoxia during an 11-min period. The results suggested that hemodynamics in microcirculation were not affected by the heat from the LEDs. Then, we conducted the in vivo experiment with rats to evaluate the sensitivity of SDF oximetry under hypoxic stimulation. This is the first in vivo experiment that observed the changes in of microcirculation based on SDF imaging. We confirmed that the estimated by SDF oximetry was in accordance with the change in . Therefore, SDF oximetry is considered to be a method that enables us to observe changes occurring in accordance with alteration of microcirculation. Moreover, we focused on not only the overall changes in microcirculation but also the difference between individual blood vessels. The speeds of these responses to hypoxic stimulation were different. These differences were thought to be due to the difference between arterioles and venules as reported previously,5 or to differences of blood vessel diameter. To raise the reliability of estimations, future work must take into account the error factors of SDF oximetry such as the influence of scattering by tissues, as mentioned in our previous study.15 As one solution to this problem, Monte Carlo simulation of photon propagation is considered. Analyzing the photon propagation in the tissue would enable researchers to understand the behavior of photons and lessen the extent of the error factors. After improving the accuracy of estimation, as future work, we will conduct in vivo experiments with a sepsis rat model by following the same procedure as described in our previous work.25 Then, we will investigate the relationship between oxygen delivery and mechanism of sepsis shock. AcknowledgmentsThis research was supported by KAKENHI (the Grants-in-aid for Scientific Research) (A) 16H01855. ReferencesR. E. Airley, J. E. Monaghan and I. J. Stratford,
“Hypoxia and disease: opportunities for novel diagnostic and therapeutic prodrug strategies,”
Pharm. J., 264
(7094), 666
–673
(2000). Google Scholar
C. G. Ellis, J. Jagger and M. Sharpe,
“The microcirculation as a functional system,”
Crit. Care, 9
(4), S3
(2005). https://doi.org/10.1186/cc3751 JBOPFO 1083-3668 Google Scholar
C. Ince and M. Sinaasappel,
“Microcirculatory oxygenation and shunting in sepsis and shock,”
Crit. Care Med., 27
(7), 1369
–1377
(1999). https://doi.org/10.1097/00003246-199907000-00031 CCMDC7 0090-3493 Google Scholar
J. W. Cannon,
“Hemorrhagic shock,”
N. Engl. J. Med., 378
(4), 370
–379
(2018). https://doi.org/10.1056/NEJMra1705649 Google Scholar
H. C. Hendargo et al.,
“Snap-shot multispectral imaging of vascular dynamics in a mouse window-chamber model,”
Opt. Lett., 40
(14), 3292
–3295
(2015). https://doi.org/10.1364/OL.40.003292 OPLEDP 0146-9592 Google Scholar
S. Hu and L. V. Wang,
“Photoacoustic imaging and characterization of the microvasculature,”
J. Biomed. Opt., 15
(1), 011101
(2010). https://doi.org/10.1117/1.3281673 Google Scholar
S. P. Chong et al.,
“Quantitative microvascular hemoglobin mapping using visible light spectroscopic optical coherence tomography,”
Biomed. Opt. Express, 6
(4), 1429
–1450
(2015). https://doi.org/10.1364/BOE.6.001429 BOEICL 2156-7085 Google Scholar
W. Groner et al.,
“Orthogonal polarization spectral imaging: a new method for study of the microcirculation,”
Nat. Med., 5
(10), 1209
–1212
(1999). https://doi.org/10.1038/13529 Google Scholar
D. De Backer et al.,
“Microvascular blood flow is altered in patients with sepsis,”
Am. J. Respir. Crit. Care Med., 166
(1), 98
–104
(2002). https://doi.org/10.1164/rccm.200109-016OC AJCMED 1073-449X Google Scholar
C. Ince,
“Sidestream dark field imaging: an improved technique to observe sublingual microcirculation,”
Crit. Care, 9
(1), P72
(2005). https://doi.org/10.1186/cc3135 Google Scholar
P. T. Goedhart et al.,
“Sidestream dark field (SDF) imaging: a novel stroboscopic LED ring-based imaging modality for clinical assessment of the microcirculation,”
Opt. Express, 15
(23), 15101
–15114
(2007). https://doi.org/10.1364/OE.15.015101 OPEXFF 1094-4087 Google Scholar
J. G. Dobbe et al.,
“Measurement of functional microcirculatory geometry and velocity distributions using automated image analysis,”
Med. Biol. Eng. Comput., 46
(7), 659
–670
(2015). https://doi.org/10.1007/s11517-008-0349-4 MBECDY 0140-0118 Google Scholar
R. Djaberi et al.,
“Non-invasive assessment of microcirculation by sidestream dark field imaging as a marker of coronary artery disease in diabetes,”
Diabetes Vasc. Dis. Res., 10
(2), 123
–134
(2013). https://doi.org/10.1177/1479164112446302 IEBEAX 0018-9294 Google Scholar
T. Kurata et al.,
“Impact of vessel diameter and bandwidth of illumination in sidestream dark-field oximetry,”
Biomed. Opt. Express, 6
(5), 1616
–1631
(2015). https://doi.org/10.1364/BOE.6.001616 BOEICL 2156-7085 Google Scholar
T. Kurata et al.,
“Correction method for influence of tissue scattering for sidestream dark-field oximetry using multicolor LEDs,”
Opt. Rev., 23
(6), 955
–967
(2016). https://doi.org/10.1007/s10043-016-0278-7 1340-6000 Google Scholar
A. J. Cohen and R. A. Laing,
“Multiple scattering analysis of retinal blood oximetry,”
IEEE Trans. Biomed. Eng., BME-23 391
–400
(1976). https://doi.org/10.1109/TBME.1976.324650 Google Scholar
D. J. Faber et al.,
“Oxygen saturation-dependent absorption and scattering of blood,”
Phys. Rev. Lett., 93
(2), 028102
(2004). https://doi.org/10.1103/PhysRevLett.93.028102 PRLTAO 0031-9007 Google Scholar
W. G. Zijlstra et al.,
“Spectrophotometry of hemoglobin: absorption spectra of bovine oxyhemoglobin, deoxyhemoglobin, carboxyhemoglobin, and methemoglobin,”
Comp. Biochem. Physiol. B: Biochem. Mol. Biol., 118
(4), 743
–749
(1997). https://doi.org/10.1016/S0305-0491(97)00230-7 CBPBB8 0305-0491 Google Scholar
A. Zwart et al.,
“A multi-wavelength spectrophotometric method for the simultaneous determination of five haemoglobin derivatives,”
Clin. Chem. Lab. Med., 19
(7), 457
–464
(1981). https://doi.org/10.1515/cclm.1981.19.7.457 Google Scholar
M. C. Meinke et al.,
“Optical properties of platelets and blood plasma and their influence on the optical behavior of whole blood in the visible to near-infrared wavelength range,”
J. Biomed. Opt., 12
(1), 014024
(2007). https://doi.org/10.1117/1.2435177 JBOPFO 1083-3668 Google Scholar
L. E. MacKenzie et al.,
“In vivo oximetry of human bulbar conjunctival and episcleral microvasculature using snapshot multispectral imaging,”
Exp. Eye Res., 149 48
–58
(2016). https://doi.org/10.1016/j.exer.2016.06.008 EXERA6 0014-4835 Google Scholar
W. Feng et al.,
“Lookup-table-based inverse model for mapping oxygen concentration of cutaneous microvessels using hyperspectral imaging,”
Opt. Express, 25
(4), 3481
–3495
(2017). https://doi.org/10.1364/OE.25.003481 OPEXFF 1094-4087 Google Scholar
F. Helmchen et al.,
“A miniature head-mounted two-photon microscope: high-resolution brain imaging in freely moving animals,”
Neuron, 31
(6), 903
–912
(2001). https://doi.org/10.1016/S0896-6273(01)00421-4 NERNET 0896-6273 Google Scholar
M. Sugimura et al.,
“Influence of acute progressive hypoxia on cardiovascular variability in conscious spontaneously hypertensive rats,”
Auton. Neurosci. Basic Clin., 141
(1), 94
–103
(2008). https://doi.org/10.1016/j.autneu.2008.05.008 Google Scholar
M. Takahashi et al.,
“Quantitative evaluation of blood flow obstruction in microcirculation with sidestream dark-field images,”
Proc. SPIE, 10068 100680A
(2017). https://doi.org/10.1117/12.2253402 PSISDG 0277-786X Google Scholar
BiographyRyohei Hashimoto is an MS student at the Graduate School of Science and Engineering, Chiba University, Japan, where he received his BS degree in engineering in 2017. Tomohiro Kurata received his MS degree from Gakushuin University in 2012 and his PhD from Chiba University in 2017. He is an R&D engineer at Takano Co., Ltd. His current research interests include image processing, diffuse reflectance spectroscopy, and multispectral imaging. Masashi Sekine received his PhD from Chiba University in 2015. He is a chief technical officer at the Center for Frontier Medical Engineering, Chiba University, and a part-time lecturer at the College of Engineering and Design, Shibaura Institute of Technology. His research interests include mechanical design, robotics, production engineering, and medical engineering. Currently, he is engaged in developing rehabilitation devices including prostheses. Kazuya Nakano received his MS and PhD degrees from Tokyo Institute of Technology in 2007 and 2013, respectively. He is a specially appointed assistant professor at the Center for Frontier Medical Engineering, Chiba University. His current research interests include diffuse reflectance spectroscopy, multispectral imaging, optical security, information, and photonics. Takashi Ohnishi received his MS and PhD degrees from Chiba University in 2010 and 2013, respectively. He is an assistant professor at the Center for Frontier Medical Engineering, Chiba University. His current research interests include image processing, image registration, and pathological images. Hideaki Haneishi received his MS degree in 1987 and his PhD in 1990 from Tokyo Institute of Technology. He joined Chiba University in 1990 as a research associate. He was a visiting research scientist at the Department of Radiology, University of Arizona, from 1995 to 1996. He has been a full professor at the Center for Frontier Medical Engineering (CFME) since 2007. Currently, he is also a director of CFME. |