28 March 2016 Blood-pressure-induced oscillations of deoxy- and oxyhemoglobin concentrations are in-phase in the healthy breast and out-of-phase in the healthy brain
Author Affiliations +
J. of Biomedical Optics, 21(10), 101410 (2016). doi:10.1117/1.JBO.21.10.101410
We present a near-infrared spectroscopy (NIRS) study of local hemodynamics in the breast and the brain (prefrontal cortex) of healthy volunteers in a protocol involving periodic perturbations to the systemic arterial blood pressure. These periodic perturbations were achieved by cyclic inflation (to a pressure of 200 mmHg) and deflation (at frequencies of 0.046, 0.056, 0.063, 0.071, and 0.083 Hz) of two pneumatic cuffs wrapped around the subject’s thighs. As a result of these systemic perturbations, the concentrations of deoxy- and oxyhemoglobin in tissue (D and O, respectively) oscillate at the set frequency. We found that the oscillations of D and O in breast tissue are in-phase at all frequencies considered, a result that we attribute to dominant contributions from blood volume oscillations. In contrast, D and O oscillations in brain tissue feature a frequency-dependent phase difference, which we attribute to significant contributions from cerebral blood flow oscillations. Frequency-resolved measurements of D and O oscillations are exploited by the technique of coherent hemodynamics spectroscopy for the assessment of cerebrovascular parameters and cerebral autoregulation. We show the relevant physiological information content of NIRS measurements of oscillatory hemodynamics, which have qualitatively distinct features in the healthy breast and healthy brain.
Tgavalekos, Kainerstorfer, Sassaroli, and Fantini: Blood-pressure-induced oscillations of deoxy- and oxyhemoglobin concentrations are in-phase in the healthy breast and out-of-phase in the healthy brain



Near-infrared spectroscopy (NIRS) is a noninvasive optical method for local measurements of deoxyhemoglobin (D) and oxyhemoglobin (O) concentrations in tissue and has been used for a variety of applications such as functional brain imaging1 and breast cancer imaging.2,3 Recently, there has been interest in spontaneous as well as induced oscillations of D and O, because their amplitudes and relative phases can be related to physiological quantities such as blood volume, blood flow, oxygen consumption, and cerebral autoregulation (AR). Oscillations of D and O occur naturally within the human body at the heart rate (1  Hz), the respiratory rate (0.25  Hz), and at the frequency of spontaneous low-frequency oscillations (0.1  Hz). Oscillations can also be induced at specific frequencies with various protocols, including paced breathing,4 thigh cuff occlusions,5 and squat-stand maneuvers.6

The phase difference between D and O oscillations is relevant to the study of both the diseased and healthy human brains. In a brain study of patients with unilateral carotid obstruction, it was found that, in a paced breathing protocol, the phase lag between D and O oscillations was significantly larger on the side with the obstruction than on the contralateral healthy side.7 These findings were attributed to impaired versus normal cerebral AR, the mechanism that maintains cerebral blood flow in response to changes in the mean arterial pressure (MAP). In a sleep study of healthy subjects, NIRS measurements from the brain revealed altered phase differences between D and O low-frequency (0.1  Hz) spontaneous oscillations when comparing non-rapid eye movement(REM) with REM sleep states,8 a finding that was interpreted in terms of oscillations in blood flow velocity and blood volume.

Coherent hemodynamics spectroscopy (CHS) is a method for studying hemodynamic physiological processes and tissue oxygen consumption in terms of the phases and amplitudes of deoxyhemoglobin, oxyhemoglobin, and total hemoglobin (T) concentrations at multiple frequencies of oscillations.9 The amplitude and phase of D, O, and T oscillations at angular frequency ω can be described with phasor notation: D(ω), O(ω), and T(ω). CHS spectra are given by the frequency dependence of the magnitude and phase of the phasor ratios D(ω)/O(ω) and O(ω)/T(ω). In phasor algebra, this corresponds to the following amplitude ratios and phase differences as a function of frequency: |D|/|O|, |O|/|T|, Arg(D)Arg(O), and Arg(O)Arg(T), where “Arg” is the argument of the phasor, i.e., its phase angle with respect to the x-axis. By means of a hemodynamic model,9 the measured CHS spectra can be related to fundamental physiological parameters like the capillary and venous blood transit times and the cutoff frequency of the cerebral AR process.

CHS can be applied to any tissue, not just the brain. In this study, we have also chosen to study the breast tissue, because CHS may yield valuable information for cancer assessment as a result of the hemodynamic and metabolic perturbations associated with breast cancer. The utility of assessing hemodynamic changes with NIRS in breast tissue has been investigated in the literature. In a study looking at breast hemodynamics, a dual-breast near-infrared tomographic system measured different responses to a Valsalva maneuver in a comparison between the healthy and tumor-bearing breasts on the same subject.10 The return to baseline of deoxyhemoglobin concentration in the tumor-bearing breast lagged behind the return to baseline in the healthy breast. Carp et al.11 have studied the breast hemodynamic response to compression in order to investigate physiological differences of healthy and diseased tissues. Breast hemodynamics during periodic oscillations have also been studied with a protocol in which subjects used a respiratory device to alternate between breathing oxygen and carbogen.12 Total hemoglobin concentration response was measured with NIRS, and the reported results support that the total hemoglobin responses of cancerous breasts were less correlated to the changes in inspired oxygen concentrations than the total hemoglobin responses of healthy breasts. Flexman et al.13 have studied the response of breast tissue to a breath-holding protocol and found differences between malignant and healthy tissues during breath holding and recovery. The spatial and temporal study of breast hemodynamics may result in a powerful tool for the characterization and detection of abnormalities in the tissue vasculature and metabolism. CHS may therefore be a useful addition to this area of study.

Oscillations in blood volume and blood flow have individual effects on the oscillations of D and O that are measured with NIRS. Blood volume oscillations, by themselves, result in synchronous (or in-phase) oscillations of D and O. Blood flow oscillations, by themselves, result in oscillations of D and O that are in opposition of phase (π-phase difference). Such individual volume and flow effects are, in general, out-of-phase with each other even in the case in which the blood volume and blood flow oscillations are synchronous. This result follows from the delayed D and O responses to blood flow changes (this delay is caused by the finite blood transit time in the microvasculature), which we modeled by treating the microvasculature as a low-pass filter for the driving blood flow dynamics (input) and the resulting D and O dynamics (output).9 According to this analysis, the phase difference between the overall, measured D and O oscillations reflects the relative contributions to them from blood volume and blood flow changes. In-phase D and O oscillations suggest volume-dominated effects, and opposition-of-phase D and O oscillations suggest flow-dominated effects. The relative contributions to D and O oscillations from oscillations in blood volume and blood flow (and, when applicable, oxygen consumption) may depend on the frequency of oscillations, and this fact is exploited by CHS for a quantitative study of tissue hemodynamics and oxidative metabolism.9

Even though blood volume and blood flow changes may be interconnected (e.g., a change in blood volume in a compartment may reflect an imbalance between flow in and flow out of that compartment), we note that NIRS is not sensitive to blood flow dynamics occurring outside the capillary compartment. One of the assumptions of our model is that blood volume changes occur synchronously in the three compartments (arterial, capillary, and venous), and for this reason, we use the general term “blood volume change” without any further clarification. On the contrary, blood flow changes in the context of this work refer only to capillary flow changes.

In this study, we have used NIRS to concurrently measure blood-pressure–induced hemoglobin oscillations in the breast and the brain (prefrontal cortex) of healthy volunteers. The measured CHS spectra in the breast and the brain show that D and O oscillations are in-phase (i.e., synchronous) in the healthy breast and out-of-phase (i.e., asynchronous) in the healthy brain. We attribute these results to different relative contributions of blood volume and blood flow to the measured oscillations of D and O in the two tissues examined. While we have focused on breast and brain tissues, a similar hemodynamic study in other tissues can yield similar information about the contribution of blood volume and blood flow to hemoglobin concentration oscillations measured with NIRS.




Human Subjects and Data Acquisition

Eleven healthy female subjects (age range: 24 to 32 years old) participated in the study. All subjects were premenopausal, had no history of vascular disorders, and had no known risks for breast cancer. The Tufts University Institutional Review Board approved the experimental protocol, and the subjects provided written informed consent prior to the experiment. Figure 1 shows the experimental setup. The NIRS measurements were performed with a frequency-domain commercial NIRS instrument (OxiplexTS, ISS Inc., Champaign, Illinois). Optical probes connected to the spectrometer delivered light at two wavelengths, 690 and 830 nm, at a source–detector distance of 35 mm. The brain probe was placed against the left side of the subject’s forehead, to access tissue in the prefrontal cortex, and secured with a flexible headband. The breast probe was placed 20  mm above the areola of the left breast and secured with adhesive medical tape. Continuous arterial blood pressure (ABP) was recorded with a beat-to-beat blood pressure monitoring system (NIBP100D, BIOPAC Systems, Inc., Goleta, California). Pneumatic thigh cuffs were wrapped around both of the subject’s thighs and connected to an automated cuff inflation system (E-20 Rapid Cuff Inflation System, D.E. Hokanson, Inc., Bellevue, Washington). The air pressure in the thigh cuffs was continuously monitored with a digital manometer (Series 626 Pressure Transmitter, Dwyer Instruments, Inc., Michigan City, Indiana). Analog outputs of the ABP monitor and the thigh cuff pressure monitor were fed to auxiliary inputs of the NIRS instrument for concurrent recordings with the NIRS data. All signals were sampled synchronously at an acquisition rate of 12.5 Hz.

Fig. 1

Experimental setup. Signals from the NIRS instrument, ABP monitor, and thigh cuff manometer were recorded synchronously. Optical data from the spectrometer were used to compute the changes in O and D (ΔO and ΔD) over time (t).


Baseline measurements, during which the subjects were resting and the thigh cuffs were deflated, were recorded for 2 min for each participant. Then the thigh cuffs were periodically inflated (to a pressure of 200  mmHg) and deflated for six periods at five different frequencies, fi(fi=ωi/2π): 0.046, 0.056, 0.063, 0.071, and 0.083 Hz. The cyclic inflations of the thigh cuffs induced periodic changes in the ABP.5 The associated changes in O and D in the brain and breast were measured with NIRS. Each set of cyclic cuff inflations was followed by 20 s of rest during which the baseline conditions were re-established.


Data Analysis

Data analysis and processing were performed with MATLAB® (Mathworks Inc., Natick, Massachusetts). Optical intensity changes from the brain and the breast were translated to changes in the concentrations of oxyhemoglobin [ΔO(t)], deoxyhemoglobin [ΔD(t)], and total hemoglobin [ΔT(t)] using the modified Beer–Lambert law. The signals were analyzed with the following process for each of the five frequencies of thigh cuff oscillations, fi (i=1,,5). The time traces were filtered with a linear phase, finite impulse response, bandpass filter based on the Parks–McClellan algorithm.14 The filter had a width of 0.01 Hz that was centered at fi. The Hilbert transform was applied to the bandpass-filtered signals in order to obtain the instantaneous amplitude and phase at that frequency.15 Experiments on tissue-like phantoms with static optical properties were used to determine the noise floor of the intensity measurements, which were translated to the detection threshold for oscillations in hemoglobin concentrations in tissue. From these phantom experiments, a threshold of 0.015  μM was considered the noise level for the amplitude of oscillations and applied to the instantaneous amplitudes of deoxy- and oxyhemoglobin to remove data points with lower amplitudes.16 The thigh cuff manometer signal was then used to determine the time ranges at which the six periods of induced oscillations of frequency fi occurred. The magnitudes and phases of D(ωi)/O(ωi) and O(ωi)/T(ωi) were computed and averaged within the time ranges of the measurements at each fi to obtain CHS spectra [namely |D|/|O|, |O|/|T|, Arg(D)Arg(O), and Arg(O)Arg(T)].

In the case of brain measurements, the CHS spectra showed typical features previously reported by our group17 and allowed for the assessment of six hemodynamic parameters by fitting the CHS spectra with a hemodynamic model.9 Briefly, the model describes how O, D, and T are related to microcirculation blood volume, blood flow, and metabolic rate of oxygen. The key equations for the model are




where S and its superscripts describe the oxygen saturation of hemoglobin in the arterial (a), capillary (c), or venous (v) compartment, and S(c) denotes a spatial average of the capillary saturation. The phasor that describes the oscillations of cerebral blood volume cbv(ω) and baseline cerebral blood volume CBV0 may also have superscripts to indicate their values for a specific vascular compartment. The Fåhraeus factor F(c) accounts for reduced hematocrit in the capillaries. Three transfer functions are notated as with subscript and superscript specifiers of the nature of the filter and vascular compartment (if applicable): RC low-pass (RC-LP) in the capillary compartment, Gaussian low-pass (G-LP) in the venous compartment, and RC high-pass (RC-HP) for AR.9 We used a nonlinear fitting procedure in MATLAB (function “lsqcurvefit”) with a trust region reflective algorithm to fit our measurements with the model. The six hemodynamic parameters fitted for in our analysis are the capillary blood transit time (t(c)), the venous blood transit time (t(v)), the ratio of capillary-to-venous baseline blood volumes corrected for the small-to-large vessel hematocrit ratio (F(c)CBV0(c)/CBV0(v)), the ratio of amplitudes of arterial-to-venous blood volume oscillations (ΔCBV(a)/ΔCBV(v)), the AR cutoff frequency (fc(AR)), and the product of the inverse of the Grubb’s exponent (k) times the ratio of venous-to-total baseline blood volumes (kCBV0(v)/CBV0). The search regions for the fitting parameters are defined by ranges found in the literature.18 The analytical expressions derived from the hemodynamic model and the simplifying assumptions relating these six hemodynamic parameters and the CHS spectra are described in detail in Kainerstorfer et al.19

In the case of breast measurements, because the phasors D(ωi), O(ωi), and T(ωi) were found to be in-phase with each other at all measured frequencies, the amplitude ratio |O(ωi)|/|T(ωi)| provides a measure of the oxygen saturation of hemoglobin in the volume-oscillating vascular compartments,9 which we indicate here with SV. This measure of hemoglobin saturation is a weighted average of the saturations of the volume-oscillating compartments, where each vascular compartment is weighted according to its relative contributions to the overall blood volume oscillations:





Figure 2 shows the bandpass-filtered time traces of changes in ABP, O, and D in the brain and the breast for subject no. 8. By inspection of Fig. 2, one can clearly see the oscillations in systemic ABP and local concentrations of oxy- and deoxyhemoglobin that are elicited by the cyclic thigh cuff occlusions (identified by the vertical gray bars in Fig. 2).

Fig. 2

Time traces for subject no. 8. The shaded regions indicate the times at which the cuffs were inflated to a pressure of 200 mmHg. The cuffs were timed such that each frequency has six periods for analysis. The order of frequencies in this figure is: 0.063, 0.083, 0.046, 0.071, and 0.056 Hz. (a) Changes in mean ABP; (b) changes in oxyhemoglobin concentration measured in brain tissue; (c) changes in deoxyhemoglobin concentration measured in brain tissue; (d) changes in oxyhemoglobin concentration measured in breast tissue; and (e) changes in deoxyhemoglobin concentration measured in breast tissue. All the temporal traces shown in this figure were bandpass filtered (0.01 to 0.11 Hz, width: 0.1 Hz), whereas data analysis at each frequency was performed by applying a narrower bandpass filter (width: 0.01 Hz) centered at the frequency of interest.


The group averages of CHS spectra were computed for the 11 subjects. A variance threshold was used to select the induced frequencies in each subject that successfully maintained a consistent phase difference: the standard deviation of Arg(D)Arg(O) and the standard deviation of Arg(O)Arg(T) both needed to be <45  deg. The number of subjects for whom the measured data passed the threshold was as follows for brain data at each frequency: 0.046 Hz, n=10; 0.056 Hz, n=11; 0.063 Hz, n=10; 0.071 Hz, n=10; and 0.083 Hz, n=11. For breast data, the corresponding numbers of subjects were as follows: 0.046 Hz, n=11; 0.056 Hz, n=11; 0.063 Hz, n=11; 0.071 Hz, n=11; and 0.083 Hz, n=10. On almost all the subjects and all the frequencies, we could measure reliable and coherent hemodynamics on both brain and breast tissues. For the cases that passed the threshold requirement, the mean phase differences and amplitude ratios (representing the CHS spectra) were computed and are shown in Fig. 3. To compute the mean values and standard deviations of Arg(D)Arg(O) and Arg(O)Arg(T), we used the standard methods in circular statistics.20 For |D|/|O| and |O|/|T|, we computed the group means, and the errors bars in Fig. 3 represent the standard error of the mean.

Fig. 3

Group-averaged CHS spectra measured on the brain (circles) and breast (triangles) of 11 human subjects at five frequencies of induced hemodynamic oscillations. The lines through the brain spectra (solid lines) are the best fits of the data with the hemodynamic model. The horizontal lines through the breast spectra (dashed lines) are at the average values over the five frequencies. (a) Phase difference Arg(D)Arg(O); (b) amplitude ratio |D|/|O|; (c) phase difference Arg(O)Arg(T); and (d) amplitude ratio |O|/|T|.


For the brain spectra, the best fits with the hemodynamic model are also plotted (solid lines in Fig. 3). The brain CHS spectra show quantitative features similar to those found in previous measurements in healthy subjects.17 Over the frequency range measured here (0.046 to 0.083 Hz), the phase spectrum for Arg(D)Arg(O) shows a negative slope of 0.8  deg/mHz, whereas the amplitude spectra have values of 0.25 for |D|/|O| and 1.2 for |O|/|T| (we observe that |O|/|T| is a phasor amplitude ratio, so that its value >1 is not incompatible with the phasor relationship T=O+D). The oscillations of D lag behind the oscillations of O by 210 deg to 235 deg, whereas the oscillations of O lag behind the oscillations of T by 5 deg to 10 deg.

The CHS spectra measured in the breast did not show frequency dependence for any of the subjects, and a horizontal straight line at the value of the mean for each spectrum is shown in Fig. 3 (dashed lines). The mean of Arg(D)Arg(O) is 2  deg±14  deg and the mean of Arg(O)Arg(T) is 0  deg±2  deg, both of which indicate a nonsignificant difference from zero, so that the oscillations of O, D, and T in the breast can be described to be in-phase. The mean value of |D|/|O| is 0.21±0.01, whereas the mean value of |O|/|T| is 0.83±0.01 (in this case of in-phase D, O, and T phasors, the ratio |O|/|T| must be 1).

The phase difference between the oscillations of deoxy- and oxyhemoglobin concentrations [Arg(D)Arg(O)] was the feature of greatest interest for this study, because it provides dynamic information about the relative contributions of volume and flow changes in the probed tissue. As an illustrative example of the phase difference of D and O oscillations, Fig. 4(a) shows the time traces of D and O changes (ΔD and ΔO) and Fig. 4(b) shows the corresponding phasor diagrams of D and O. Figure 4 refers to subject no. 1 and a frequency of cyclic cuff inflation of 0.056 Hz [the time periods of cuff inflation are shaded in Fig. 4(a)]. The time traces in Fig. 4(a) show the out-of-phase behavior of deoxy- and oxyhemoglobin oscillations in brain tissue (top panels of Fig. 4) and the in-phase behavior of the oscillations in breast tissue (bottom panels of Fig. 4). In this example, Arg(D)Arg(O) is 185  deg±2  deg in the brain and 4  deg±1  deg in the breast. In Fig. 4(b), O is the phase reference (set at 0 deg), and the phasor amplitudes are drawn to scale with respect to each other.

Fig. 4

(a) Time traces of changes in oxy- and deoxyhemoglobin concentrations (ΔO and ΔD, respectively) measured during the cyclic cuff occlusion protocol at 0.056 Hz for subject no. 1. The periods of time when the thigh cuffs were inflated are shaded. The top panel shows the out-of-phase behavior of O and D oscillations in the healthy brain. The bottom panel reports the breast data and shows the in-phase behavior of O and D oscillations in the healthy breast. (b) Phasor diagrams corresponding to the oscillations reported in (a). The phasor diagram in the top panel shows the out-of-phase O and D phasors in the healthy brain, whereas the one in the bottom panel shows the approximately in-phase O and D phasors in the healthy breast.


The fits to the brain CHS spectra with the hemodynamic model (solid lines in Fig. 3) yielded the set of six hemodynamic parameters defined in Sec. 2.2. The mean values of these parameters (± the standard error of the mean) over the 11 subjects were as follows: t(c)=0.9±0.2  s, t(v)=5.4±1.1  s, (F(c)CBV0(c)/CBV0(v))=1.3±0.3, (ΔCBV(a)/ΔCBV(v))=5.2±1.7, fc(AR)  =0.03±0.02  Hz, (kCBV0(v)/CBV0)=0.8±0.2. Most parameters are in good agreement with values previously reported for the prefrontal cortex of 11 healthy human subjects during a paced breathing protocol.19 Specifically, the previously reported values of the six parameters were 0.9±0.2  s, 1.3±0.3  s, 1.1±0.3, 2.9±0.9, 0.035±0.002  Hz, and 0.6±0.1, respectively.19 While some of the parameters deviate from our previous results, the reproducibility of t(c) and fc(AR) confirms the robustness of quantitative CHS measurements of cerebral blood flow (which is inversely related to t(c)) and cerebral AR (whose effectiveness scales with fc(AR)).

The case of breast CHS spectra is qualitatively different. Because Arg(O)Arg(T) and Arg(D)Arg(O) are essentially zero, we can conclude that the blood volume oscillations are the dominant source of hemoglobin concentration oscillations. Therefore, oxy-, deoxy-, and total hemoglobin concentrations all oscillate in-phase. Under these conditions, |O|/|T| specifies the average hemoglobin saturation of the volume-oscillating vascular compartments SV, for which we found a value of 83%±1% [see Fig. 3(d)]. We hypothesize that the highly compliant venous compartment provides the greatest contributions to these volume oscillations driven by blood pressure perturbations. This is similar to the case of spiroximetry, which was proposed to measure venous saturation from the blood volume oscillations associated with respiration (which also modulates ABP). In studies where spiroximetry was applied to skeletal muscles, the average venous saturation was found within the range 70% to 80%.17,21 In an NIRS study using venous occlusion to measure venous saturation in the human forearm, S(v) ranged from 50% to 80%.22 The value found by us for SV in the breast is somewhat greater than the reported values of venous saturation in skeletal muscle. This result may be due to a different balance of blood flow and oxygen consumption in breast tissue compared with skeletal muscle or to greater contributions from the arterial compartment to   SV.



The model described in Refs. 9 and 23 indicates that a pure blood flow oscillation induces a phase shift of π between D and O, whereas a pure blood volume oscillation induces synchronous oscillations of D and O. Equations 1 and 2, introduced in Sec. 2.2, are the analytical equations in the hemodynamic model that we used in order to solve for the six unknown parameters. By solving D(ω)/O(ω) and O(ω)/T(ω) instead of solving separately for the phasors D(ω), O(ω), and T(ω), we can cancel out unknown or frequency-dependent scaling factors for the amplitudes and phases of oscillations. However, to best explain the in-phase versus out-of-phase contributions from blood volume and blood flow oscillations, the equations for D(ω) and O(ω) are shown here as a sum of contributions from blood volume (subscript V), blood flow (subscript F), and oxygen consumption (subscript O)


O(ω)=  OV(ω)+OF(ω)+OO(ω)


D(ω)=  DV(ω)+DF(ω)+DO(ω)

Table 1 contains the terms and signs associated with each of the contributions.9,23 In the results reported here, we have assumed that there was no modulation of the cerebral metabolic rate (CMRO2) or tissue oxygen consumption during the experiments. Table 1 (first and second rows) shows that an oscillation in any component (arterial, capillary, or venous) of the blood volume generates in-phase oscillations OV(ω) and DV(ω). Table 1 also shows that a given blood flow phasor cbf(ω) (i.e., blood flow oscillations) generates phase-lagged oscillations OF(ω) [because of the negative phase of the complex factor in square brackets in the third row of Table 1 (for details, see Ref. 9)] and oscillations DF(ω) that are in opposition of phase with OF(ω) (because of the opposite sign of the terms in the third and fourth rows of Table 1). More precisely, even if blood volume and blood flow oscillations are in-phase, OF(ω) and DF(ω) are lagging behind OV(ω) and DV(ω). Note also that OF(ω) and DF(ω) are frequency-dependent, reflecting the nature of the low-pass filters associated with the capillary and venous compartments. A combination of both blood flow and blood volume oscillations induces frequency-dependent oscillations O(ω) and D(ω), which feature a relative phase that depends on the relative contributions of flow and volume oscillations, the blood transit time in the microvasculature, and the presence of any phenomena (such as AR mechanisms) that link flow and volume changes.

Table 1

Contributions to O(ω) and D(ω) from oscillations in blood volume (subscript V) or blood flow (subscript F) at frequency ω. The oscillations in blood volume and blood flow are represented by the phasors cbv(ω) and cbf(ω), respectively. This table also shows the effect of oscillations in the metabolic rate of oxygen [phasor cmro2(ω)] that are neglected in this work. ctHb is the concentration of hemoglobin in blood. Other symbols are defined in the text.


Figures 3(a) and 4 show that deoxy- and oxyhemoglobin oscillations are in-phase in the breast tissue and out-of-phase in the brain. This is an indication that the hemodynamics measured with NIRS are dominated by blood volume dynamics in healthy breast tissue, whereas they result from both blood volume and blood flow dynamics in brain tissue. We hypothesize that this contrasting behavior is due, at least in part, to a difference in the elastic properties of breast and brain tissues; in particular, the presence of the rigid skull limits the vascular expansion in brain tissue that, instead, is not constrained in breast tissue. Therefore, one may expect that the contributions of blood volume and blood flow to D and O are different in breast and brain tissues.

We observe that while the relative phase of D and O oscillations is related to the level of flow AR, one may not associate a lack of phase shift between D and O oscillations with a lack of AR. In fact, even in the absence of AR, when blood flow oscillations passively follow MAP oscillations so that the cerebral blood flow and MAP (and, to a first approximation, cerebral blood volume) are in-phase with each other, D and O oscillations are still, in general, out-of-phase with each other. This is because, as described in Sec. 1 and shown in Table 1, OV(ω) is in-phase with cbv(ω), whereas OF(ω) lags behind cbf(ω) as a result of the finite blood transit time in the microvasculature.

The phase relationship between D and O oscillations in breast tissue may be affected by the presence of breast cancer. Breast tumors have a larger elastic modulus than healthy breast tissue,24 and they are also known to have an abnormal vasculature.25 For both of these reasons, D and O oscillations in breast tumors may no longer be dominated by blood volume effects, as in healthy breast tissue, since the stiffer tissue and abnormal blood vessels may result in a less compliant tumor vasculature. Based on the model and results described here, D and O oscillations would become increasingly out-of-phase with each other as the contribution of blood volume oscillations decreases.

The in-phase oscillations of deoxy- and oxyhemoglobin in healthy breast tissue permitted us to compute the oxygen saturation of hemoglobin in the volume-oscillating compartments, SV=|OV|/|T|, simply as the ratio |O|/|T|. To measure the saturation of the volume-oscillating compartments in the brain, one must take into account the blood flow contributions to the measured hemoglobin concentration. Kainerstorfer et al.26 have described in detail how to account for contributions from blood flow in order to measure SV in a general case where deoxy- and oxyhemoglobin oscillations are not in-phase with each other.

This work underlines the value of quantitative measurements of hemodynamic oscillations, be they induced by controlled perturbations in the ABP (as done here) or spontaneously occurring. However, in order to interpret such quantitative measurements of oscillatory hemodynamics, it is important to understand the dynamic relationship between the measured quantities (deoxy- and oxyhemoglobin concentrations in the case of NIRS) and the underlying physiological processes. In particular, this work highlights the qualitatively distinct behavior of different tissues (brain and breast) in response to the same systemic perturbation in ABP.

For noninvasive optical studies of the human brain, it is important to consider the fact that, in addition to brain tissue, extracerebral tissue (scalp, skull, and so on) also contributes to the measured optical signals. In relation to the study reported here, this means that extracerebral hemodynamics may act as a confounding factor for the cerebral hemodynamics that we intend to investigate and compare with breast tissue hemodynamics. However, while potentially confounding contributions from extracerebral hemodynamics may affect our measurements of the relative phase shift of O and D oscillations, they would not affect the main finding of this work that cerebral O and D oscillations are out-of-phase, whereas O and D oscillations in the healthy breast are in-phase. Nevertheless, in an effort to quantify the effect of extracerebral hemodynamics on our dynamic NIRS measurements, we are currently exploring the application of a two-layer diffusion model to discriminate superficial and deep-tissue hemodynamics.27



This study has shown that blood-pressure–induced hemodynamic oscillations behave qualitatively differently in the healthy breast and in the healthy brain. Specifically, the different phase delays between D and O oscillations observed in the breast (0  deg in the range 0.04 to 0.08 Hz) and the brain (200  deg at 0.06 Hz) may be attributed to different relative contributions of blood flow and blood volume oscillations to the measured oxy- and deoxyhemoglobin concentrations. The situation may be different in pathologic conditions, which may affect the vascular compliance, the microvascular architecture, normal physiological mechanisms, and the elastic properties of tissue. In this case, a quantitative assessment of hemodynamic oscillations, possibly as a function of frequency as done in CHS, may offer diagnostically relevant information for pathologic conditions. The broad implications of this work are that CHS may be sensitive to perturbations to the integrity of the vasculature and hemodynamics in the brain (i.e., subarachnoid hemorrhage, traumatic brain injury, stroke, and so on) and in the breast (i.e., breast cancer). A thorough characterization of hemodynamic oscillations in tissue and their accurate interpretation can have a significant impact in the study of microvascular circulation and the assessment of its integrity.


We thank Anaïs Leproux, Michael Ghijsen, and Bruce Tromberg for useful discussions. This research was supported by the National Institutes of Health (Grant No. R01-CA154774).



M. Ferrari and V. Quaresima, “A brief review on the history of human functional near-infrared spectroscopy (fNIRS) development and fields of application,” Neuroimage 63(2), 921–935 (2012).NEIMEF1053-8119http://dx.doi.org/10.1016/j.neuroimage.2012.03.049Google Scholar


B. J. Tromberg et al., “Assessing the future of diffuse optical imaging technologies for breast cancer management,” Med. Phys. 35(6), 2443–2451 (2008).MPHYA60094-2405http://dx.doi.org/10.1118/1.2919078Google Scholar


S. Fantini and A. Sassaroli, “Near-infrared optical mammography for breast cancer detection with intrinsic contrast,” Ann. Biomed. Eng. 40(2), 398–407 (2012).ABMECF0090-6964http://dx.doi.org/10.1007/s10439-011-0404-4Google Scholar


M. Reinhard et al., “Spatial mapping of dynamic cerebral autoregulation by multichannel near-infrared spectroscopy in high-grade carotid artery disease,” J. Biomed. Opt. 19(9), 097005 (2014).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.19.9.097005Google Scholar


R. Aaslid et al., “Asymmetric dynamic cerebral autoregulatory response to cyclic stimuli,” Stroke 38(5), 1465–1469 (2007).SJCCA70039-2499http://dx.doi.org/10.1161/STROKEAHA.106.473462Google Scholar


J. A. H. R. Claassen, B. D. Levine and R. Zhang, “Dynamic cerebral autoregulation during repeated squat-stand maneuvers,” J. Appl. Physiol. 106, 153–160 (2009).http://dx.doi.org/10.1152/japplphysiol.90822.2008Google Scholar


M. Reinhard et al., “Oscillatory cerebral hemodynamics—the macro- vs. microvascular level,” J. Neurol. Sci. 250, 103–109 (2006).JNSCAG0022-510Xhttp://dx.doi.org/10.1016/j.jns.2006.07.011Google Scholar


M. L. Pierro et al., “Phase-amplitude investigation of spontaneous low-frequency oscillations of cerebral hemodynamics with near-infrared spectroscopy: a sleep study in human subjects,” Neuroimage 63(3), 1571–1584 (2012).NEIMEF1053-8119http://dx.doi.org/10.1016/j.neuroimage.2012.07.015Google Scholar


S. Fantini, “Dynamic model for the tissue concentration and oxygen saturation of hemoglobin in relation to blood volume, flow velocity, and oxygen consumption: implications for functional neuroimaging and coherent hemodynamics spectroscopy (CHS),” Neuroimage 85, 202–221 (2014).NEIMEF1053-8119http://dx.doi.org/10.1016/j.neuroimage.2013.03.065Google Scholar


C. H. Schmitz et al., “Design and implementation of dynamic near-infrared optical tomographic imaging instrumentation for simultaneous dual-breast measurements,” Appl. Opt. 44(11), 2140–2153 (2005).APOPAI0003-6935http://dx.doi.org/10.1364/AO.44.002140Google Scholar


S. A. Carp et al., “Dynamic functional and mechanical response of breast tissue to compression,” Opt. Express 16(20), 16064–16078 (2008).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.16.016064Google Scholar


C. M. Carpenter et al., “Inspired gas-induced vascular change in tumors with magnetic-resonance-guided near-infrared imaging: human breast pilot study,” J. Biomed. Opt. 15(3), 036026 (2010).JBOPFO1083-3668http://dx.doi.org/10.1117/1.3430729Google Scholar


M. L. Flexman et al., “Optical biomarkers for breast cancer derived from dynamic diffuse optical tomography,” J. Biomed. Opt. 18(9), 096012 (2013).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.18.9.096012Google Scholar


T. Parks and J. McClellan, “Chebyshev approximation for nonrecursive digital filters with linear phase,” IEEE Trans. Circuit Theory 19(2), 189–194 (1972).IECTAF0018-9324http://dx.doi.org/10.1109/TCT.1972.1083419Google Scholar


B. Boashash, “Estimating and interpreting the instantaneous frequency of a signal—part 1: fundamentals,” Proc. IEEE 80(4), 520–538 (1992).IEEPAD0018-9219http://dx.doi.org/10.1109/5.135376Google Scholar


M. L. Pierro et al., “Reduced speed of microvascular blood flow in hemodialysis patients versus healthy controls: a coherent hemodynamics spectroscopy study,” J. Biomed. Opt. 19(2), 026005 (2014).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.19.2.026005Google Scholar


M. L. Pierro et al., “Validation of a novel hemodynamic model for coherent hemodynamics spectroscopy (CHS) and functional brain studies with fNIRS and fMRI,” Neuroimage 85, 222–233 (2014).NEIMEF1053-8119http://dx.doi.org/10.1016/j.neuroimage.2013.03.037Google Scholar


J. M. Kainerstorfer, A. Sassaroli and S. Fantini, “Coherent hemodynamics spectroscopy in a single step,” Biomed. Opt. Express 5(10), 3403–3616 (2014).BOEICL2156-7085http://dx.doi.org/10.1364/BOE.5.003403Google Scholar


J. M. Kainerstorfer et al., “Practical steps for applying a new dynamic model to near-infrared spectroscopy measurements of hemodynamic oscillations and transient changes: implications for cerebrovascular and functional brain studies,” Acad. Radiol. 21(2), 185–196 (2014).http://dx.doi.org/10.1016/j.acra.2013.10.012Google Scholar


J. H. Zar, “Circular distributions,” in Biostatistical Analysis, 5th ed., pp. 612–623, Pearson Prentice Hall, Upper Saddle River, New Jersey (2010).Google Scholar


M. A. Franceschini et al., “Near-infrared spiroximetry: noninvasive measurements of venous saturation in piglets and human subjects,” J. Appl. Physiol. 92, 372–384 (2002).Google Scholar


C. W. Yoxail and A. M. Weindling, “Measurement of venous oxyhaemoglobin saturation in the adult human forearm by near infrared spectroscopy with venous occlusion,” Med. Biol. Eng. Comput. 35, 331–336 (1997).MBECDY0140-0118http://dx.doi.org/10.1007/BF02534086Google Scholar


S. Fantini, “A new hemodynamic model shows that temporal perturbations of cerebral blood flow and metabolic rate of oxygen cannot be measured individually using functional near-infrared spectroscopy,” Physiol. Meas. 35(1), N1–N9 (2014).PMEAE30967-3334http://dx.doi.org/10.1088/0967-3334/35/1/N1Google Scholar


A. Samani, J. Zubovits and D. Plewes, “Elastic moduli of normal and pathological human breast tissues: an inversion-technique-based investigation of 169 samples,” Phys. Med. Biol. 52(6), 1565–1576 (2007).PHMBA70031-9155http://dx.doi.org/10.1088/0031-9155/52/6/002Google Scholar


P. Vaupel, F. Kallinowski and P. Okunieff, “Blood flow, oxygen and nutrient supply, and metabolic microenvironment of human tumors: a review,” Cancer Res. 49(23), 6449–6465 (1989).CNREA80008-5472Google Scholar


J. M. Kainerstorfer, A. Sassaroli and S. Fantini, “Optical oximetry of volume-oscillating vascular compartments: contributions from oscillatory blood flow,” J. Biomed. Opt. 21(10), 101408 (2016).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.21.10.101408Google Scholar


S. Fantini et al., “Non-invasive assessment of cerebral microcirculation with diffuse optics and coherent hemodynamics spectroscopy,” Proc. SPIE 9690, 96900S (2016).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.2213139Google Scholar

Biographies for the authors are not available.

Kristen T. Tgavalekos, Jana M. Kainerstorfer, Angelo Sassaroli, Sergio Fantini, "Blood-pressure-induced oscillations of deoxy- and oxyhemoglobin concentrations are in-phase in the healthy breast and out-of-phase in the healthy brain," Journal of Biomedical Optics 21(10), 101410 (28 March 2016). http://dx.doi.org/10.1117/1.JBO.21.10.101410
Submission: Received ; Accepted






Blood circulation

Near infrared spectroscopy

Back to Top