Significance: We demonstrate the potential of probing the sO2 change under blood flow in vivo using photoacoustic (PA) imaging and sheds light on the complex relationship between RBC aggregation and oxygen delivery. Aim: To conduct in vivo assessments of the sO2 in the radial artery of healthy volunteers and simultaneously probe the relation between the sO2 and hemodynamic behavior such as red blood cell (RBC) aggregation. Approach: The effects of PA-based measurements of blood hemodynamics were studied as a function of the subjects’ age (20s, 30s, and 40s). The pulsatile blood flow in the human radial artery of 12 healthy subjects was imaged in the 700 to 900 nm optical wavelength range using a linear array-based PA system. Results: The PA power when blood velocity is minimum (Pamax) was larger than the one attained at maximum blood velocity (Pamin), consistent with predictions based on the cyclical variation of RBC aggregation during pulsatile flow. The difference between Pamin and Pamax at 800 nm (ΔPa800) increased with age (1.7, 2.2, and 2.6 dB for age group of 20s, 30s, and 40s, respectively). The sO2 computed from Pamax was larger than the one from Pamin. Conclusions: The ΔPa800 increased with participant age. The ΔPa800 metric could be a surrogate of noninvasively monitoring the age-induced changes in RBC aggregation. The sO2 change during a cycle of pulsatile blood flow also increased with age, demonstrating that RBC aggregation can affect the sO2 change. |
1.IntroductionRed blood cells (RBCs) play an important, physiologically significant role in the human body, affecting hemodynamics as well as governing an oxygen transport. The oxygen-carrying capacity of RBCs is driven by the partial pressure of oxygen in blood vessels and the tissues they supply and is systemically assessed through metrics of oxygen saturation ().1 The hemodynamic aspects of RBCs, namely their flow profiles, pressures, and velocities, across the cardiovascular system play an equally important physiological role.2 Pathophysiological changes in the hemodynamic behavior of RBCs affect a wide variety of disorders.3 A phenomenon that is affected by the hemodynamic behaviors of RBCs is their aggregability, defined as the ability of the cells to form rouleaux in the presence of plasma proteins.4 This naturally occurring phenomenon is attributed to the cyclical variation of the shear forces forming when blood velocity gradients change during flow.5 It is thus affected by the hemodynamic behavior of vessels6 as well as RBC membrane factors that impact the electrostatic and steric interactions between cells.7 The clinical significance of RBC aggregation becomes apparent in conditions where hyperaggregability has been reported. These include myocardial infarctions, bacterial infections, type 2 diabetes, and sickle cell disease.5 The study of RBC aggregation as a hemodynamic phenomenon has provided invaluable insights into the biophysical properties of this process.8,9 However, measurements of the in the presence of RBC aggregation have not been possible, impeding the assessment of the significance aggregation to the process of oxygen delivery. The is universally recognized as one of the vital signs in medicine along with pulse rate, respiratory rate, temperature, and blood pressure.10,11 This metric measures the percentage of hemoglobin binding sites that are bound with oxygen in RBCs in the blood. Oxygen metabolism is tightly regulated in the body because hypoxemia (abnormally low blood oxygen level) can lead to acute adverse effects in several organ systems.10,12–15 Recently, the decreased of COVID-19 patients measured at hospital admission was reported as one of the strongest predictors of critical illness, with silent hypoxia concurrently occurring with other symptoms such as shortness of breath.16–18 In general, clinical measurements are conducted using a pulse oximeter that attaches at the tip of the patient’s finger and is based on optical spectral differences between oxyhemoglobin (HbO) and deoxyhemoglobin (HbD). Pulse oximeters produce a measure of oxygenation which is sensitive only to the arterial vessels in the finger, not the veins. Another relevant clinical measurement is the arterial , or “,” a metric that is assessed using invasive blood draws from the radial artery followed by co-oximetry or blood gas analysis or indwelling catheter probes.19 However, in cases when peripheral perfusion is poor, as in states of hypovolemia, hypothermia, vasoconstriction, low cardiac output, and low mean arterial pressure, pulse oximeter readings become unreliable.20 Moreover, a variety of other clinical conditions such as vaso-occlusive crises in sickle cell anemia, severe anemia, methemoglobinemia, sepsis, or septic shock can influence the accuracy of pulse oximeter readings.19 Imaging techniques such as blood oxygen level-dependent magnetic resonance imaging (BOLD MRI),21 positron emission tomography (PET),22 or diffuse optical tomography or optical spectroscopy23,24 can assess the tissue oxygenation in vivo. Still, their millimeter-scale spatial resolutions limit their ability to assess the blood directly. In particular, BOLD MRI signals measure only HbD changes,25 which are likely independent of the blood while PET uses ionizing radiation. Furthermore, these techniques are not portable and cannot be used to monitor in a clinical setting or in a point of care setting. For the noninvasive measurement of in more centrally located vessels, near-infrared spectroscopy (NIRS)-based oximetry is widely used.26 However, NIRS-based tissue oximeter cannot be used in a simultaneous assessment of the relationship between the and blood dynamics. Photoacoustic (PA) imaging is an emerging technique in biomedical optics. It can provide not only anatomical structure but also functional information, complementing conventional ultrasound (US) imaging. Also, it provides greater penetration depth than conventional optical imaging methods since it relies on the detection of acoustic waves rather than ballistic photons.27,28 The high absorbance of hemoglobin within RBCs and the oxygen-dependent absorption of hemoglobin allows PA imaging to probe the microvasculature , in vitro, as well as in vivo.29–31 PA imaging is capable of noninvasive imaging of blood flow as well as simultaneous assessment of , being utilized as the state-of-the-art tool for clinical applications such as investigation on the correlation between tissue characterization and oximetry. However, no PA studies have ever examined the impact of blood’s hemodynamic behavior in the estimation of this important parameter. There are several methods to assess RBC aggregation, such as erythrocyte sedimentation rate, centrifugation methods, microscopic methods, low shear viscometry, US imaging, and analysis of light transmission (or reflection) of RBC suspension.32 The noninvasive measurement of RBC aggregation has been widely investigated using US imaging.33–35 RBC aggregation has also been used as a clinical biomarker for inflammation monitoring among other pathologies.8,9,36 However, US imaging probes the structure of RBC aggregates and does not directly measure physiological function of RBCs. As such, the assessment of both and hemodynamic behavior under blood flow needs to be independently assessed using noninvasive biomedical imaging using conventional optics37–40 and US,33–35 respectively, but not in the same setting. Our group has reported on in vitro measurements that show the correlation between the and RBC aggregation using PA spectroscopy in static and flow conditions.41–43 Through carefully controlled in vitro conditions, we demonstrated that during a pulsatile cycle in a simulated pulsatile blood flow, the PA amplitude increased with decreasing flow velocity and decreased with increasing velocity. These changes in PA amplitude were due to RBC aggregation during small velocities and disaggregation during increased velocities and consistent with the findings widely reported using conventional US methodology.33,35 In addition to RBC aggregation, the PA-based also cyclically changed in-phase with RBC aggregation. These results suggest that RBC aggregates inhibited release. In this paper, we describe the first in vivo evidence for the RBC aggregation-induced alterations in the measured in the radial artery of healthy subjects. Furthermore, a complex age-dependent correlation between the and RBC is discussed. 2.Materials and Methods2.1.Subject Recruitment and PA Imaging ProtocolThe human subject study was approved by the Research Ethic Board of Ryerson University (REB 2017-040). Healthy subjects were recruited based on the inclusion criteria listed in Table 1. Table 1Subject criteria and information.
Note: N: number of subjects, BMI: body mass index (kg/m2). The subject sat down on a chair and immersed their left arm in a degassed warm (36°C) water bath (Fig. 1). PA imaging was performed with a linear-array probe equipped PA imaging system (Vevo LAZR; LZ250-21 MHz of central frequency, 13 to 24 MHz of bandwidth and 256 elements, FUJIFILM Visualsonics, Toronto, Ontario, Canada)44 at the near-infrared wavelength (700, 750, 800, 850, and 900 nm). The pulsed-wave Doppler velocity () was measured to locate the radial artery, and the measurement system was then switched to the PA imaging mode. A PA B-mode image was acquired to measure the pulsatile blood flow in the radial artery for 10 s at each optical wavelength (). The total time that subjects’ arms were immersed in water was . The time trace was acquired for each single wavelength then time-shifted to create a combined multispectral dataset in postprocessing. The detailed methods were described in our previous study.42,45 2.2.Data Acquisition and Postprocessing: V, Pa, and sO2A representative coregistered US (gray scale) and PA (color scale) image of the radial artery was shown in Fig. 2(a). The second quarter of the full field of view was chosen as a region of interest (ROI) to maximize the PA image acquisition rate to 20 Hz. To avoid boundary-buildup or edge detection artifact from the upper and lower vessel walls, the upper and lower boundaries of the ROI were chosen at 10% margin within the vessel lumen. In the coregistered US image, the upper and lower vessel walls were tracked at each frame. For each wavelength and each subject, 200 frames of US images () were analyzed to track the vessel wall. From the upper and lower walls tracked, the upper and lower boundaries for the ROI were computed by 10% margin within the vessel lumen. For example, if the locations of two walls were at 1 and 2 mm in depth (corresponding to and , respectively), the locations of two boundaries of ROI were 1.1 and 1.9 mm in depth (corresponding to and , respectively) as shown in Fig. 2(a). The horizontal width of the ROI was fixed as a second quarter of the full field of view, whereas the vertical length (related to the vessel diameter) of the ROI was dependent on each frame. The VevoLAZR PA imaging system provides access to the prebeamformed RF data for all 256 transducer elements. These data are then beamformed postacquisition, and the amplitude of each signal is used to reconstruct the PA images shown in Fig. 2(a). The ROI is then selected from this reconstructed PA image. From beamformed radiofrequency (RF) PA signals (64 out of 256 elements) in the ROI, the PA power () was computed by taking an average of the root-mean-square of each RF signal for each , and addressed by where is the laser energy in J, which measured by Vevo LAZR, , and represent the time in the RF signal, corresponding to the upper and lower boundaries of the ROI, respectively, as shown in Fig. 2(a), P is PA amplitude in voltage, is the element number of US detector. The velocity () of the pulsatile blood flow measured for 10 s. The number of pulsatile cycles was dependent on the subject (ranging from 11 to 15 cycles during the 10-s measurement for all subjects). The was estimated using the optical absorption at two applied wavelengths, which was derived from46 where is the absorption coefficient in /cm; and are the known molar extinction coefficients of oxygenated hemoglobin (HbO) and deoxygenated hemoglobin (HbD), respectively, in /cm/; ; [HbO] and [HbD] are the molar concentrations of HbO and HbD, respectively, in ; and are two different (700 and 900 nm in this study). can hence substitute for . measured at 700 and 900 nm were synchronized to calculate the using Eq. (2), as shown in Fig. 2(b). was acquired at each of the wavelengths for 10 s, and then the wavelength was switched to the next one. The time difference between acquisition at wavelengths was 20 s, which is 10 s for wavelength switching and 10 s for acquisition at the wavelength. The assessment using our technique has been validated in previous studies47 and independently validated in pre- or clinical studies.48 Since the position and geometry of illuminated area were the same, the fluence effect could be minimized. Also, the effect of wavelength on fluence was somewhat compensated by using normalization to the energy at each wavelength. The representatives of , oxygen saturation (, red), and the Doppler velocity (, black) were shown in Fig. 2(b).Fig. 2(a) Representative US (gray scale) and PA, PA (color scale) image of the radial artery. The PA RF power () from the ROI was computed by taking an average of the root-mean-square of each RF signal for each optical wavelength (). (b) Representatives of (900-nm-thick blue and 700-nm-thin blue), oxygen saturation (, red), and the Doppler velocity (, black). The subscripts max and min represent the maximum (denoted as a circle) and minimum (denoted as a diamond) value for each parameter. ![]() The minimum ( and ) and maximum ( and ) values of and were, respectively, averaged for 10 s to compare the variation in and during the pulsatile blood flow as a function of and age. Two sample -test was conducted using built-in MALTAB function “ttest2.m,” in terms of the and changes versus age group. 3.Results3.1.PA Power versus Optical WavelengthThe PA magnitude was dependent on the optical as well as the study subjects. To compare this parameter, the and for all were normalized to the at an isosbestic point 800 nm (the point at which HbO and HbD have the same optical absorption) for each subject. The dependence on the for each subject age group is shown in Fig. 3(a). The at 800 nm was set to 0 dB as a reference [green circles in Fig. 3(a)]. Both and increased with for all groups, as shown in Fig. 3(a). In addition, the difference between and () increased with the . The average values of , , and at 700, 800, and 900 nm for all groups were tabulated (Table 2). The difference between and at 800 nm () for all subjects is shown in Fig. 3(b)-left. The error bars represent the standard deviation from the number of cycles of the pulsatile blood flow for each subject. The box and whisker plot based on the average values of of each subject for age groups is shown in Fig. 3(b)-right. The increased with age, i.e., 1.8, 2.2, and 2.7 dB for the age group of subjects in their 20s, 30s, and 40s, respectively. Fig. 3(a) PA power () as a function of optical wavelength for each group of age 20s, 30s, and 40s. The thick and thin blue lines represent and , respectively. The thick red line represents the variation from to (). The green circle indicates the isosbestic point (HbO and HbD have the same optical absorption). The error bars represent the standard deviation from the number of cycles of the pulsatile flow for all subjects, i.e., 75, 39, and 38 cycles for each age group, respectively. (b) The difference between and at an isosbestic point 800 nm () for all subjects and the corresponding box and whisker plot (based on the average value of for each subject) for each age group. The error bars represent the standard deviation from the number of cycles of the pulsatile flow for each subject. : not significant, *** : , **** : . ![]() Table 2The values of minimum (Pamin) and maximum (Pamax) of PA power and the difference between Pamin and Pamin (ΔPa) at 700, 800, and 900 nm for all groups. The error bars represent the standard deviation from the number of cycles of the pulsatile flow for all subjects.
3.2.sO2 versus Subject AgeThe estimated from () and () for each group of age 20s, 30s, and 40s were shown in Fig. 4, respectively. The average for subjects in their 20s, 30s, and 40s were 98.7%, 97.2%, and 96.7%, respectively [Fig. 5(a)]. This indicates that the difference in between the youngest and oldest subjects was 2.0% (). On the other hand, the average for subjects in their 20s, 30s, and 40s were 97.1%, 94.7%, and 93.0%, respectively [Fig. 5(a)]. The difference between the 20s and 40s groups was 4.1% (). The difference between and () during a pulsatile cycle increased with age, i.e., 1.6%, 2.5%, and 3.8% for the age group of subjects in their 20s, 30s, and 40s, respectively [Fig. 5(b)]. Fig. 4Oxygen saturation () for all subjects—subjects in their 20s (blue dots), subjects in their 30s (red diamonds), and subjects in their 40s (green square). The was computed for () and (), respectively. The error bars represent the standard deviation from the number of cycles of the pulsatile flow. ![]() Fig. 5(a) Oxygen saturation () for each group based on the average value of of each subject. The thick () and thin () symbols represent the computed using the two wavelength (700 and 900 nm) methods at and , respectively. (b) The difference between and () for each group. : not signigicant, *: , **: , ****: . ![]() 4.DiscussionSeveral studies have demonstrated the PA assessment of the .28,46,49,50 However, the complexity of hemodynamic and/or hemorheological behaviors such as RBC aggregation/disaggregation and deformation have not been considered as a possible factor that could alter the estimated . Our group has studied the feasibility of using PA spectroscopy for assessing the relation between RBC aggregation and under pulsatile blood flow in vitro using experimental and theoretical methods.41,42 This paper describes the first-in-human feasibility study on the dependence of the on the pulsatile blood flow conditions of the radial artery. The experimental results show an age dependence of the measurements in vivo. 4.1.PA Power versus Age – Aggregability of RBCsSeveral studies reported on the correlation between RBC aggregation and age.51–53 Woodward et al.51 measured hemorheological variables (blood viscosity, RBC aggregation, and fibrinogen) from subjects (25 to 74 years old) in association with cardiovascular risk factors. They found that RBC aggregation and fibrinogen increased with age. Christy et al.52 reported a significant increase in RBC aggregation with age (20 to 59 years old), establishing a significant correlation between phagocytic activity and RBC aggregability. According to Simmonds et al.53 the mechanism of the age-related increase in RBC aggregation is that aging decreased the electrostatic repulsive forces between cell surfaces, thus promoting RBC aggregation alongside with increased plasma fibrinogen concentration. During pulsatile blood flow, the dominant hemorheological characteristic is the cyclical aggregation and disaggregation phases of RBCs.5,32 During systole, the blood flow velocity is maximum, generating maximal shear rate within vessels, and as a result, the RBCs in the radial artery flow as single cells. During diastole, on the other hand, the blood flow velocity is at a minimum, resulting in a minimal shear rate, leading to the formation of RBCs rouleaux.5,32 In PA imaging, the RBC aggregates form a larger effective absorber compared to single cells, increasing the signal.41,42,54,55 In Fig. 2(b), the and Doppler velocity were out of phase during a systolic-diastolic cycle. A higher is expected in the presence of aggregates and a lower from single RBCs. This occurred at all wavelengths of illumination, and as we have shown in earlier studies, it can approximate the aggregate size.41 The PA amplitude is a function of (including the absorption cross-section) and an absorber size. According to Eq. (1), can be represented as a linear combination of , , [HbO] and [HbD], depending on . The can be derived by combining at two wavelengths as addressed in Eq. (2). Given the also varies with the absorber size affecting the , the must have a nonlinear relation with the absorber size. However, the at an isosbestic point 800 nm is not dependent on the , resulting in a linear relation between the and the absorber size, as previously demonstrated by our group.41 The increased with age [Fig. 3(b)], suggesting that the RBC aggregability also increased with age (since the PA signal is a surrogate metric of the aggregate size). 4.2.Effect of RBC Aggregation on PA Power as a Function of Absorber Size and sO2None of the abovementioned studies have examined the impact of RBC aggregation on the blood measurements in humans. It is well known that RBC aggregation and plasma viscosity play an important role in determining the overall blood viscosity.56 Realizing the fundamental function of RBCs in oxygen transport, Tateishi et al.57,58 were the first to postulate on the correlation between RBC aggregation and release from erythrocytes and its impact on the .59 In addition, the relation between RBC aggregation and the was theoretically and experimentally investigated through in vitro PA imaging by our group.41 Similar to the Tateishi’s findings and our in vitro experiments, the in vivo results in this work demonstrate the impact of RBC aggregation on the measurement. The increases due to two factors: an increase in the size of the absorber (single cells versus RBC aggregates) and the change in . The average values of as a function of for the group of subjects in their 40s are shown in Fig. 6(a). The represents measurements when nonaggregated RBCs flow in the radial artery at maximal velocity (when the shear rates inside the radial artery are highest). RBC aggregation increases the due to an increase in optical absorber size. Since the optical absorber size increases as RBCs aggregate, the increases by (2.7 dB for subjects in their 40s) for all (, depicted by the blue-black arrow). This can be schematically represented by in Fig. 6(b). The second reason for the changes in the is the change in the due to RBC aggregation. The difference between and “” can be identified as the contribution to the change in caused by RBC aggregation to the overall PA signal, (). Combining Eqs. (1) and (2), the can be expressed as a function of , Fig. 6(a) The effect of RBC aggregation on PA power () in relation to an optical absorber size and oxygen saturation (). (b) Schematic diagram of the effect of RBC aggregation on in relation to the absorber size and the . and represent the size of single RBC and the increased size due to aggregation, respectively. ![]() Since is larger than for ,60 the -induced increase in results in an increased with , as shown in Fig. 6. On the other hand, since is smaller than for , the -induced decrease in results in a decreased . This can be schematically denoted by and for and , respectively, in Fig. 6(b). The in vivo experimental results demonstrating the effect of RBC aggregation on the and how this depends on the RBC are also supported by our group’s previous in vitro work mimicking the radial artery flow conditions.41 This interpretation applies to all age groups. 4.3.Relationship between sO2 versus Agecould be derived from by combining the effect of the increase in both absorber size and due to RBC aggregation. Hence, the computed from the () were higher than that from the () for all groups [as shown in Fig. 5(a)]. In addition, both and decreased with increasing age. It has been reported that older persons have a lower than their younger counterparts.61–65 Even though the direct relation between the and age has not been fully investigated, it was reported that dysphagia (swallowing difficulties) could mediate the relation between and age.65 Older people are more likely to experience dysphagia which impairs pulmonary function and lowers . The present experimental results indicated that the decreases with age, despite the subject cohort comprising healthy individuals. The difference between and () increased with age, as shown in Fig. 5(b). This phenomenon was consistent with the relation between and age, as shown in Fig. 3(b). The represents the effects of RBC aggregability, independent of , which increased with age. Since oxygen release is inhibited by RBC aggregation,57,58 higher aggregability results in more hemoglobin molecules to be bound oxygen molecules. Despite a negative correlation between the and age, the RBC aggregation-induced increase in the () resulted in a positive correlation between and age. This is why the correlation between the and age is steeper than the correlaiton between the and age, as shown in Fig. 5(a). 4.4.Limitations of the StudyThe findings of this study suggest the feasibility of the PA assessment of both and its correlation to age, in vivo. However, there are limitations to this study that form the basis for future extension of this work. The study, whose recruitment is suspended during the global COVID-19 pandemic, enrolled a smaller than desired subject group and was mainly composed of male volunteers. The sex dependence on related to the correlation between and age must be further investigated. The laser fluctuation which sometimes occur from pulse to pulse or from wavelength to wavelength should be measured and corrected. The absolute values for the presented in this paper have not been corrected for the effects of laser fluence.66 The radial artery depth from the skin surface is . As such, the effects of fluence (spectral coloring)67 might not be as significant as it is in other PA applications. Specifically, in Eq. (3), the is calculated by a ratio of . The PA amplitude () is a function of and fluence (), so that the can be represented by the ratio of “.” Moreover, since the position and geometry were the same, variations in fluence are minimized. In addition, the normalization to the laser energy at each wavelength could also contribute to a reduction of the fluence effects. As such, the reported absolute values might not represent the true for each subject. In addition, this quantity is dependent on the subject’s skin color (leading to a larger optical path length to the radial artery, requiring fluence correction). An additional limitation is that the upper limit of age group was in 40s. Further investigation for the older age groups (50s, 60, and 70s) should be conducted to study the age dependence further. Despite these limitations, this study demonstrates the feasibility of the measurements and confirms previous in vitro findings. Future studies could be done using a portable probe, opening the potential for doing this investigation more easily. In addition, quantitative measurement of RBC aggregation by US such as the structure-factor-size-estimation68 should be further applied to this study to correlate both PA and US modalities. 5.ConclusionIn this study, the correlation between the and RBC aggregation under the pulsatile blood flow in the human radial artery, and its age dependence, was investigated using PA imaging. RBC aggregability increased with age, as observed by the age-induced increase in the metric during RBC aggregation. The change during a cycle of pulsatile blood flow also increased with age, and it was attributed to the presence of increased RBC aggregation in older subjects. This study is the first study to examine how the in vivo changes in during blood flow in human body can be assessed using PA imaging. This study demonstrates the effect of RBC aggregation on the change during a cardiac cycle in healthy volunteers and the increase in the change with age. Although these preliminary observations were conducted in only 12 healthy participants, this work demonstrates the feasibility of the measurement in vivo. Such measurements might shed light on the clinical importance of the complex relation between the blood viscosity induced by RBC aggregation and the oxygen delivery related to . As a noninvasive measurement, PA imaging of blood could be extended to the detection of blood pathologies that alter the viscosity, modifying flow behaviors and oxygen delivery. DisclosuresThe authors have no relevant financial interests in the manuscript and no other potential conflicts of interest to disclose. AcknowledgmentsThis work was funded by the Natural Sciences and Engineering Research Council of Canada/Canadian Institutes of Health Research – Collaborative Health Research Projects Grant No. 462315-2014 and an NSERC Discovery Grant (RGPIN-2017-06496). Funding to purchase the equipment was provided by the Canada Foundation for Innovation, the Ontario Ministry of Research and Innovation, and Ryerson University. E. Hysi is supported by Banting and KRESCENT Fellowships by the Canadian Institutes of Health Research and the Kidney Foundation of Canada, respectively. We thank Elizabeth Berndl at the Department of Physics at Ryerson University for providing technical support. ReferencesM. L. Ellsworth, R. N. Pittman and C. G. Ellis,
“Measurement of hemoglobin oxygen saturation in capillaries,”
Am. J. Physiol. Circ. Physiol., 252
(5), H1031
–H1040
(1987). https://doi.org/10.1152/ajpheart.1987.252.5.H1031 Google Scholar
T. W. Secomb,
“Hemodynamics,”
Compr. Physiol., 6
(2), 975
–1003
(2016). https://doi.org/10.1002/cphy.c150038 Google Scholar
R. N. Mitchell,
“Hemodynamic disorders thromboembolic disease and shock,”
Robbins and Cotran Pathologic Basis of Disease, 113
–136 9th ed.Elsevier(2014). Google Scholar
G. Barshtein, R. Ben-Ami and S. Yedgar,
“Role of red blood cell flow behavior in hemodynamics and hemostasis,”
Expert Rev. Cardiovasc. Ther., 5
(4), 743
–752
(2007). https://doi.org/10.1586/14779072.5.4.743 Google Scholar
O. K. Baskurt, B. Neu and H. J. Meiselman, Red Blood Cell Aggregation, CRC Press, Boca Raton, Florida
(2011). Google Scholar
H. H. Lipowsky,
“Microvascular rheology and hemodynamics,”
Microcirculation, 12
(1), 5
–15
(2005). https://doi.org/10.1080/10739680590894966 Google Scholar
S. Chien,
“Biophysical behavior of red cells in suspensions,”
The Red Blood Cell: Volume 2, 1031
–1133 2nd ed.Academic Press, New York
(1975). Google Scholar
F. T. H. Yu et al.,
“A local increase in red blood cell aggregation can trigger deep vein thrombosis: evidence based on quantitative cellular ultrasound imaging,”
J. Thromb. Haemost., 9
(3), 481
–488
(2011). https://doi.org/10.1111/j.1538-7836.2010.04164.x Google Scholar
J. Tripette et al.,
“In vivo venous assessment of red blood cell aggregate sizes in diabetic patients with a quantitative cellular ultrasound imaging method: proof of concept,”
PLoS One, 10
(4), e0124712
(2015). https://doi.org/10.1371/journal.pone.0124712 Google Scholar
D. W. Green and G. Kunst,
“Cerebral oximetry and its role in adult cardiac, non‐cardiac surgery and resuscitation from cardiac arrest,”
Anaesthesia, 72
(Suppl. 1), 48
–57
(2017). https://doi.org/10.1111/anae.13740 Google Scholar
R. A. L. Walton and B. D. Hansen,
“Venous oxygen saturation in critical illness,”
J. Vet. Emerg. Crit. Care, 28
(5), 387
–397
(2018). https://doi.org/10.1111/vec.12749 Google Scholar
B. D. Kent, P. D. Mitchell and W. T. Mcnicholas,
“Hypoxemia in patients with COPD: cause, effects, and disease progression,”
Int. J. Chron. Obstruct. Pulmon. Dis., 6
(1), 199
–208
(2011). https://doi.org/10.2147/COPD.S10611 Google Scholar
M. R. Hanudel et al.,
“Effects of acute kidney injury and chronic hypoxemia on fibroblast growth factor 23 levels in pediatric cardiac surgery patients,”
Pediatr. Nephrol., 31
(4), 661
–669
(2016). https://doi.org/10.1007/s00467-015-3257-5 Google Scholar
A. Puttappa et al.,
“Prolonged unexplained hypoxemia as initial presentation of cirrhosis: a case report,”
Am. J. Case Rep., 18 1
–6
(2017). https://doi.org/10.12659/AJCR.900530 Google Scholar
B. D. Kussman et al.,
“Cerebral oxygen saturation in children with congenital heart disease and chronic hypoxemia,”
Anesth. Analg., 125
(1), 234
–240
(2017). https://doi.org/10.1213/ANE.0000000000002073 Google Scholar
A. Thomas, G. Haljan and A. Mitra,
“Lung ultrasound findings in a 64-year-old woman with COVID-19,”
Can. Med. Assoc. J., 192
(15), E399
(2020). https://doi.org/10.1503/cmaj.200414 Google Scholar
R. G. Wilkerson et al.,
“Silent hypoxia: a harbinger of clinical deterioration in patients with COVID-19,”
Am. J. Emerg. Med., 38
(10), 2243.e5
–2243.e6
(2020). https://doi.org/10.1016/j.ajem.2020.05.044 Google Scholar
M. Petrilli et al.,
“Factors associated with hospital admission and critical illness among 5279 people with coronavirus disease 2019 in New York City: prospective cohort study,”
Br. Med. J., 369 m1966
(2020). https://doi.org/10.1136/bmj.m1966 Google Scholar
E. D. Chan, M. M. Chan, M. M. Chan,
“Pulse oximetry: understanding its basic principles facilitates appreciation of its limitations,”
Respir. Med., 107
(6), 789
–799
(2013). https://doi.org/10.1016/j.rmed.2013.02.004 Google Scholar
P. A. Kyriacou et al.,
“Investigation of oesophageal photoplethysmographic signals and blood oxygen saturation measurements in cardiothoracic surgery patients,”
Physiol. Meas., 23
(3), 533
–545
(2002). https://doi.org/10.1088/0967-3334/23/3/305 Google Scholar
T. Christen, D. S. Bolar and G. Zaharchuk,
“Imaging brain oxygenation with MRI using blood oxygenation approaches: methods, validation, and clinical applications,”
Am. J. Neuroradiol., 34
(6), 1113
–1123
(2013). https://doi.org/10.3174/ajnr.A3070 Google Scholar
N. Fleming et al.,
“Imaging tumour hypoxia with positron emission tomography,”
Br. J. Cancer, 112
(2), 238
–250
(2015). https://doi.org/10.1038/bjc.2014.610 Google Scholar
Y. Hoshi and Y. Yamada,
“Overview of diffuse optical tomography and its clinical applications,”
J. Biomed. Opt., 21
(9), 091312
(2016). https://doi.org/10.1117/1.JBO.21.9.091312 Google Scholar
N. C. Biswal, Y. Xu and Q. Zhu,
“Imaging tumor oxyhemoglobin and deoxyhemoglobin concentrations with ultrasound-guided diffuse optical tomography,”
Technol. Cancer Res. Treat., 10
(5), 417
–429
(2011). https://doi.org/10.7785/tcrt.2012.500219 Google Scholar
R. G. Evans et al.,
“Don’t be so BOLD: potential limitations in the use of BOLD MRI for studies of renal oxygenation,”
Kidney Int., 71
(12), 1327
–1328
(2007). https://doi.org/10.1038/sj.ki.5002321 Google Scholar
P. B. Benni et al.,
“A validation method for near-infrared spectroscopy based tissue oximeters for cerebral and somatic tissue oxygen saturation measurements,”
J. Clin. Monit. Comput., 32
(2), 269
–284
(2018). https://doi.org/10.1007/s10877-017-0015-1 Google Scholar
P. Beard,
“Biomedical photoacoustic imaging,”
Interface Focus, 1
(4), 602
–631
(2011). https://doi.org/10.1098/rsfs.2011.0028 Google Scholar
M. Xu and L. V. Wang,
“Photoacoustic imaging in biomedicine,”
Rev. Sci. Instrum., 77
(4), 041101
(2006). https://doi.org/10.1063/1.2195024 Google Scholar
F. Cao et al.,
“Photoacoustic imaging in oxygen detection,”
Appl. Sci., 7
(12), 1262
(2017). https://doi.org/10.3390/app7121262 Google Scholar
S. Hu,
“Emerging concepts in functional and molecular photoacoustic imaging,”
Curr. Opin. Chem. Biol., 33 25
–31
(2016). https://doi.org/10.1016/j.cbpa.2016.04.003 Google Scholar
S. Wang et al.,
“Recent advances in photoacoustic imaging for deep-tissue biomedical applications,”
Theranostics, 6
(13), 2394
–2413
(2016). https://doi.org/10.7150/thno.16715 Google Scholar
O. K. Baskurt and H. J. Meiselman,
“Erythrocyte aggregation: basic aspects and clinical importance,”
Clin. Hemorheol. Microcirc., 53
(1–2), 23
–37
(2013). https://doi.org/10.3233/CH-2012-1573 Google Scholar
G. Cloutier and Z. Qin,
“Ultrasound backscattering from non-aggregating and aggregating erythrocytes—a review,”
Biorheology, 34
(6), 443
–470
(1997). https://doi.org/10.3233/BIR-1997-34607 Google Scholar
T.-H. Bok et al.,
“Feasibility study of high-frequency ultrasonic blood imaging in human radial artery,”
J. Med. Biol. Eng., 35
(1), 21
–27
(2015). https://doi.org/10.1007/s40846-015-0001-3 Google Scholar
D.-G. Paeng, K.-H. Nam and K. K. Shung,
“Cyclic and radial variation of the echogenicity of blood in human carotid arteries observed by harmonic imaging,”
Ultrasound Med. Biol., 36
(7), 1118
–1124
(2010). https://doi.org/10.1016/j.ultrasmedbio.2010.03.021 Google Scholar
B. Chayer et al.,
“Pilot clinical study of quantitative ultrasound spectroscopy measurements of erythrocyte aggregation within superficial veins,”
Clin. Hemorheol. Microcirc., 74
(2), 109
–126
(2020). https://doi.org/10.3233/CH-180541 Google Scholar
M. J. R. Polson and G. L. Morris,
“Method of measuring the oxygen saturation in pulsating blood flow,”
(1995). Google Scholar
T. Aoyagi,
“Pulse oximetry: its invention, theory, and future,”
J. Anesth., 17
(4), 259
–266
(2003). https://doi.org/10.1007/s00540-003-0192-6 Google Scholar
E. Sinex,
“Pulse oximetry: principles and limitations,”
Am. J. Emerg. Med., 17
(1), 59
–66
(1999). https://doi.org/10.1016/S0735-6757(99)90019-0 Google Scholar
E. Roussakis et al.,
“Oxygen‐sensing methods in biomedicine from the macroscale to the microscale,”
Angew. Chemie Int. Ed., 54
(29), 8340
–8362
(2015). https://doi.org/10.1002/anie.201410646 Google Scholar
T.-H. Bok, E. Hysi and M. C. Kolios,
“In vitro photoacoustic spectroscopy of pulsatile blood flow: probing the interrelationship between red blood cell aggregation and oxygen saturation,”
J. Biophotonics, 11 e201700300
(2018). https://doi.org/10.1002/jbio.201700300 Google Scholar
T.-H. Bok, E. Hysi and M. C. Kolios,
“Simultaneous assessment of red blood cell aggregation and oxygen saturation under pulsatile flow using high-frequency photoacoustics,”
Biomed. Opt. Express, 7
(7), 2769
–2780
(2016). https://doi.org/10.1364/BOE.7.002769 Google Scholar
Y. X. Huang et al.,
“Human red blood cell aging: correlative changes in surface charge and cell properties,”
J. Cell. Mol. Med., 15
(12), 2634
–2642
(2011). https://doi.org/10.1111/j.1582-4934.2011.01310.x Google Scholar
A. Needles et al.,
“Development and initial application of a fully integrated photoacoustic micro-ultrasound system,”
IEEE Trans. Ultrason. Ferroelectr. Freq. Control, 60
(5), 888
–897
(2013). https://doi.org/10.1109/TUFFC.2013.2646 Google Scholar
T.-H. Bok, E. Hysi and M. C. Kolios,
“Preliminary photoacoustic imaging of the human radial artery for simultaneous assessment of red blood cell aggregation and oxygen saturation in vivo,”
in IEEE Int. Ultrason. Symp.,
(2017). https://doi.org/10.1109/ULTSYM.2017.8091879 Google Scholar
X. Wang et al.,
“Noninvasive imaging of hemoglobin concentration and oxygenation in the rat brain using high-resolution photoacoustic tomography,”
J. Biomed. Opt., 11
(2), 024015
(2006). https://doi.org/10.1117/1.2192804 Google Scholar
R. Pinto et al.,
“Feasibility of photoacoustic imaging for the non‐invasive quality management of stored blood bags,”
Vox Sang., 114 701
–710
(2019). https://doi.org/10.1111/vox.12828 Google Scholar
M. Lakshman and A. Needles,
“Screening and quantification of the tumor microenvironment with micro-ultrasound and photoacoustic imaging,”
Nat. Methods, 12
(4), iii
–v
(2015). https://doi.org/10.1038/nmeth.f.381 Google Scholar
H. F. Zhang et al.,
“Imaging of hemoglobin oxygen saturation variations in single vessels in vivo using photoacoustic microscopy,”
Appl. Phys. Lett., 90
(5), 53901
–53903
(2007). https://doi.org/10.1063/1.2435697 Google Scholar
Z. Deng and C. Li,
“Noninvasively measuring oxygen saturation of human finger-joint vessels by multi-transducer functional photoacoustic tomography,”
J. Biomed. Opt., 21
(6), 061009
(2016). https://doi.org/10.1117/1.JBO.21.6.061009 Google Scholar
M. Woodward et al.,
“Associations of blood rheology and interleukin-6 with cardiovascular risk factors and prevalent cardiovascular disease,”
Br. J. Haematol., 104
(2), 246
–257
(1999). https://doi.org/10.1046/j.1365-2141.1999.01158.x Google Scholar
R. M. Christy et al.,
“Erythrocyte aggregation and neutrophil function in an aging population,”
Gerontology, 56
(2), 175
–180
(2010). https://doi.org/10.1159/000242461 Google Scholar
M. J. Simmonds, H. J. Meiselman and O. K. Baskurt,
“Blood rheology and aging,”
J. Geriatr. Cardiol., 10
(3), 291
–301
(2013). https://doi.org/10.3969/j.issn.1671-5411.2013.03.010 Google Scholar
R. K. Saha and M. C. Kolios,
“A simulation study on photoacoustic signals from red blood cells,”
J. Acoust. Soc. Am., 129
(5), 2935
–2943
(2011). https://doi.org/10.1121/1.3570946 Google Scholar
E. Hysi, R. K. Saha and M. C. Kolios,
“Photoacoustic ultrasound spectroscopy for assessing red blood cell aggregation and oxygenation,”
J. Biomed. Opt., 17
(12), 125006
(2012). https://doi.org/10.1117/1.JBO.17.12.125006 Google Scholar
S. Chien et al.,
“Blood viscosity: influence of erythrocyte deformation,”
Science, 157
(3790), 827
–829
(1967). https://doi.org/10.1126/science.157.3790.827 Google Scholar
N. Tateishi et al.,
“ release from erythrocytes flowing in a narrow -permeable tube: effects of erythrocyte aggregation,”
Am. J. Physiol.—Hear. Circ. Physiol., 281
(1), H448
–H456
(2001). https://doi.org/10.1152/ajpheart.2001.281.1.H448 Google Scholar
N. Tateishi et al.,
“Reduced oxygen release from erythrocytes by the acceleration-induced flow shift, observed in an oxygen-permeable narrow tube,”
J. Biomech., 35
(9), 1241
–1251
(2002). https://doi.org/10.1016/S0021-9290(02)00068-4 Google Scholar
N. Maeda,
“Erythrocyte rheology and oxygen transfer in microcirculatory system,”
J. Biorheol., 29
(1), 2
–5
(2015). https://doi.org/10.17106/jbr.29.2 Google Scholar
S. Prahl,
“Tabulated molar extinction coefficient for hemoglobin in water,”
(1999) https://omlc.org/spectra/hemoglobin/summary.html Google Scholar
I. Cerveri et al.,
“Reference values of arterial oxygen tension in the middle-aged and elderly,”
Am. J. Respir. Crit. Care Med., 152
(3), 934
–941
(1995). https://doi.org/10.1164/ajrccm.152.3.7663806 Google Scholar
M. L. Vold et al.,
“Low oxygen saturation and mortality in an adult cohort: the Tromsø study,”
BMC Pulm. Med., 15
(1), 9
(2015). https://doi.org/10.1186/s12890-015-0003-5 Google Scholar
A. Rodríguez‐Molinero et al.,
“Normal respiratory rate and peripheral blood oxygen saturation in the elderly population,”
J. Am. Geriatr. Soc., 61
(12), 2238
–2240
(2013). https://doi.org/10.1111/jgs.12580 Google Scholar
R. Perez-Padilla et al.,
“Prevalance of oxygen desaturation and use of oxygen at home in adults at sea level and at moderate altitude,”
Eur. Respir. J., 27
(3), 594
–599
(2006). https://doi.org/10.1183/09031936.06.00075005 Google Scholar
N. Colodny,
“Effects of age, gender, disease, and multisystem involvement on oxygen saturation levels in dysphagic persons,”
Dysphagia, 16
(1), 48
–57
(2001). https://doi.org/10.1007/s004550000042 Google Scholar
M. N. Fadhel et al.,
“Fluence-matching technique using photoacoustic radiofrequency spectra for improving estimates of oxygen saturation,”
Photoacoustics, 19 100182
(2020). https://doi.org/10.1016/j.pacs.2020.100182 Google Scholar
S. Tzoumas et al.,
“Eigenspectra optoacoustic tomography achieves quantitative blood oxygenation imaging deep in tissues,”
Nat. Commun., 7 12121
(2016). https://doi.org/10.1038/ncomms12121 Google Scholar
F. T. H. Yu and G. Cloutier,
“Experimental ultrasound characterization of red blood cell aggregation using the structure factor size estimator,”
J. Acoust. Soc. Am., 122
(1), 645
–656
(2007). https://doi.org/10.1121/1.2735805 Google Scholar
BiographyTaehoon Bok is a postdoctoral fellow at Ryerson University. He received his BS and MS degrees in oceanography and ultrasound from Hanyang University in 2004 and 2006, respectively, and his PhD in ultrasound from Jeju National University in 2011. His current research interests include photoacoustic imaging, hemodynamics, oxygen delivery, and Monte Carlo simulation. He is a member of SPIE. Eno Hysi received his PhD from the Department of Physics at Ryerson University in May 2020 and is currently a Banting and KRESCENT Fellow at St. Michael’s Hospital, Division of Nephrology in Toronto. His research interests lie in the fields of ultrasound and photoacoustic imaging for biomedical applications include cancer treatment monitoring, the study of blood hemodynamics as well as assessing transplant organ quality. Michael C. Kolios received his BSc degree in physics (minor in computer science) from the University of Waterloo in 1991 and MSc degree and PhD degrees in medical biophysics from the University of Toronto in 1994 and 1998, respectively. He started his lab at Ryerson University in 1997 and is currently a professor in the Department of Physics and an associate dean, Research and Graduate Studies in the Faculty of Science. His research interests are in biomedical applications of ultrasound and optics. |