|
Near-infrared spectroscopy (NIRS) is a powerful technique for noninvasive studies of biological tissue because of its relatively large penetration depth into tissue (centimeters), and its strong sensitivity to the concentration and oxygenation of hemoglobin. In particular, NIRS allows for measurements of the temporal dynamics of oxyhemoglobin ([HbO]) and deoxyhemoglobin ([Hb]) concentrations in cerebral tissue, which in human subjects feature natural oscillations at frequencies of (arterial pulsation), (respiration), [low-frequency oscillations (LFOs)],1, 2, 3 and [very low frequency oscillations (VLFOs)].2 Functional studies of brain activation that are based on a block protocol (consisting of regular repetitions of activation-rest sequences) also typically investigate hemoglobin concentration changes at the specific frequency of the stimulation blocks.2, 3 Under these common conditions, in which one is interested in hemoglobin concentration oscillations at a specific frequency or over a relatively narrow frequency band, we propose to represent such oscillations using phasor notation. Phasors are polar 2-D vectors characterized by amplitude and phase , and represent harmonic oscillations at a fixed angular frequency . In NIRS measurements of tissue hemodynamics, a specific frequency can be isolated by applying a relatively narrow bandpass filter, and the intrinsic phase and amplitude variations can be minimized (and monitored) by considering time windows corresponding to a few periods of oscillations. We define the relative phase of two measured oscillations, say and , at angular frequency , in terms of the delay time , at which the absolute value of the cross-correlation function of the two signals is maximized. Specifically, the relative phase of and is defined as [ if ]. We define the phasor amplitude as one half the peak-to-peak range of the measured oscillations, so that the amplitude of the phasor is given by . The concept of phasor representation for the oscillations of [Hb] and [HbO] is applied to the model of Fig. 1 . In this model, we consider one [Fig. 1] or two [Fig. 1] cubic regions within tissue, where [HbO] and [Hb] oscillate at a given frequency ( in this example). The hemoglobin oscillations in the time-varying tissue regions are described by the phasors and in the respective regions, whereas the hemoglobin oscillations measured at two distinct locations are described by the phasors and at two detector positions. All hemoglobin concentration phasor amplitudes in Fig. 1 are expressed in . Additional details about the model are provided in the caption of Fig. 1. The measured phasors are obtained by first translating optical intensity changes at two wavelengths, , into global tissue absorption oscillations by using the modified Beer-Lambert law4 with proper values of the differential path-length factor (DPF) at the two wavelengths considered. Then, the tissue absorption changes are expressed in terms of hemoglobin concentration oscillations by introducing the molar extinction coefficients of oxyhemoglobin and deoxyhemoglobin, and , respectively: The measured global tissue absorption changes can also be expressed in terms of the total mean photon path length and the partial mean photon path lengths in the tissue regions that feature hemoglobin oscillations [ in Fig. 1 and in Fig. 1] as follows:By combining Eqs. 1, 2 and solving the resulting linear system of two equations for and , one obtains: where , and are the ratios of partial over total mean path lengths . Equations 3, 4 show that the measured changes in oxyhemoglobin and deoxyhemoglobin concentrations are a linear combination of the hemoglobin concentration changes in the oscillating regions. The linear combinations of and in Eqs. 3, 4 correspond to weighted vector sums of the corresponding phasors, according to phasor algebra. This means that if the phase differences between the oscillations of and are all either 0 or , then it is not possible for the phase difference between and to assume any value other than 0 (in-phase oscillations) or (opposition-of-phase oscillation). This remains true in general, regardless of tissue heterogeneity, cross talk contributions in Eqs. 3, 4, or incorrect assumptions for the DPF values at the wavelengths used. This result is important, because physiological hemodynamic or oxygenation changes all induce variations in and that either have the same sign (0 phase difference) (vascular dilation/constriction, increased/decreased blood volume) or opposite signs ( phase difference) (change in oxygen consumption, speed of blood flow, or oxygenation of incoming blood).3, 5 It is only a time difference between the onset of such physiological effects (or a phase shift different than 0 or between them) that can account for a phase shift between and that is neither 0 nor .5 The two (or more) time-delayed or phase-shifted physiological effects can either occur in the same tissue region [Fig. 1] or in different tissue regions [Fig. 1]. In the first case [Fig. 1], there is a non-0, non- phase shift between and in a specific tissue region or vascular compartment. In the second case [Fig. 1], each tissue region (or vascular compartment) features in-phase or opposition-of-phase oscillations of oxyhemoglobin and deoxyhemoglobin concentrations { and }, but there is a non-0, non- phase shift between oscillations in different tissue regions. We note that the two cases can be experimentally discriminated, because in the one-region case of Fig. 1, the measured phase difference between oxyhemoglobin and deoxyhemoglobin concentrations ( and ) is independent of the locations of light source and optical detector [if the cross talk term in Eqs. 3, 4 is negligible], which is not the case of the multiple-region case [two regions in Fig. 1], because the positions of source and detector affect the values of in Eqs. 3, 4. We think that considering the interplay of time-shifted physiological processes is the key to interpreting phase shifts of [Hb] and [HbO] that are different than 0 or , which have been repeatedly observed and reported in the literature for spontaneous hemoglobin oscillations1, 2, 6 or hemodynamic responses to brain activation.3 The word “interplay” in the previous sentence (which has also been used in this context previously1) can be interpreted here in two conceptually different ways. 1. Two physiological processes within a specific tissue region actually change their relative timing, thus causing a change in the phase difference between and . 2. A redistribution of background hemoglobin concentration and/or saturation (for example, due to functional brain activation) results in a change of the relative weights in the linear combinations of tissue hemoglobin phasors [Eqs. 3, 4], so that even in the absence of any changes in the relative phases of and (for all values of ), one would still observe a change in the phase difference between and .Figure 2 further illustrates the concept of phasor representation of hemoglobin concentration oscillations by applying it to near-infrared data measured in vivo. The time traces in Fig. 2 represent bandpass filtered temporal evolutions of [Hb], [HbO], and [HbT] measured with near-infrared spectroscopy on the prefrontal cortex of a human subject during a previously reported mental workload experiment.7 Briefly, the subject was asked to keep track of the number of colored areas in a rotating cube for . The absence of colored areas is a case referred to as Workload 0 [Fig. 2], whereas the case of four different colors is referred to as Workload 4 [Fig. 2]. The periods of oscillations of the time traces of Fig. 2 ( ) and Fig. 2 ( ) are comparable. The variability in the frequency of oscillations is something to monitor because it stretches the applicability of the concept of phasors. The phasor defined from the experimental [HbT] trace is typically in excellent agreement with the phasor sum of and (see Fig. 2), confirming the validity of the proposed phasor representation. Because the cross-correlation function between [Hb] and [HbO] is a periodic function (after bandpass filtering of the hemoglobin concentration time traces), it can also be represented as a phasor whose amplitude represents the level of correlation between [Hb] and [HbO], and whose phase represents their relative phase difference. We propose to use such cross-correlation phasors to represent the temporal trends of phase relationships and correlation levels between oxyhemoglobin and deoxyhemoglobin oscillations at a specific frequency. Cross-correlation functions measured in the prefrontal cortex during Workloads 0 and 4 are presented in Fig. 2, together with their associated phasors . One can see that the level of correlation between [Hb] and [HbO], i.e., the amplitude of the cross-correlation phasors, remains about the same under Workload 0 and Workload 4 conditions (0.98 and 0.96, respectively). By contrast, the phase of the phasor, i.e., the phase difference between [Hb] and [HbO] oscillations, is greater during Workload 4 than during Workload 0 . Cross-correlation phasors can be used to generate dynamic vectorial maps of hemoglobin concentration oscillations. This can be a powerful tool to characterize networks of functional connectivity in the brain, providing more information than the commonly used correlation coefficient analysis,8 which does not take into account and can be negatively impacted by time differences (or phase shifts) between examined signals. ReferencesG. Taga, Y. Konishi, A. Maki, T. Tachibana, M. Fujiwara, and H. Koizumi,
“Spontaneous oscillation of oxy- and deoxy-hemoglobin changes with a phase difference throughout the occipital cortex of newborn infants observed using non-invasive optical topography,”
Neurosci. Lett., 282 101
–104
(2000). https://doi.org/10.1016/S0304-3940(00)00874-0 0304-3940 Google Scholar
H. Obrig, M. Neufang, R. Wenzel, M. Kohl, J. Steinbrink, K. Einhäupl, and A. Villringer,
“Spontaneous low frequency oscillations of cerebral hemodynamics and metabolism in human adults,”
Neuroimage, 12 623
–639
(2000). https://doi.org/10.1006/nimg.2000.0657 1053-8119 Google Scholar
M. Wolf, U. Wolf, V. Toronov, A. Michalos, L. A. Paunescu, J. Choi, and E. Gratton,
“Different time evolution of oxyhemoglobin and deoxyhemoglobin concentration changes in the visual and motor cortices during functional stimulation: a near-infrared spectroscopy study,”
Neuroimage, 16 704
–712
(2002). https://doi.org/10.1006/nimg.2002.1128 1053-8119 Google Scholar
D. T. Delpy, M. Cope, P. van der Zee, S. Arridge, S. Wray, and J. Wyatt,
“Estimation of optical pathlength through tissue from direct time of flight measurement,”
Phys. Med. Biol., 33 1433
–1442
(1988). https://doi.org/10.1088/0031-9155/33/12/008 0031-9155 Google Scholar
S. Fantini,
“A hemodynamic model for the physiological interpretation of in vivo measurements of the concentration and oxygen saturation of hemoglobin,”
Phys. Med. Biol., 47 N249
–N257
(2002). https://doi.org/10.1088/0031-9155/47/18/402 0031-9155 Google Scholar
M. Reinhard, E. Wehrle-Wieland, D. Grabiak, M. Roth, B. Guschlbauer, J. Timmer, C. Weiller, and A. Hetzel,
“Oscillatory cerebral hemodynamics-the macro-vs. microvascular level,”
J. Neurosurg. Sci., 250 103
–109
(2006). 0390-5616 Google Scholar
A. Sassaroli, F. Zheng, L. H. Hirshfield, A. Girouard, E. T. Solovey, R. J. K. Jacob, and S. Fantini,
“Discrimination of mental workload levels in human subjects with functional near-infrared spectroscopy,”
J. Innov. Opt. Sci., 1 227
–237
(2008). https://doi.org/10.1142/S1793545808000224 Google Scholar
B. R. White, A. Z. Snyder, A. L. Cohen, S. E. Petersen, M. E. Raichle, B. L. Schlaggar, and J. P. Culver,
“Resting-state functional connectivity in the human brain revealed with diffuse optical tomography,”
Neuroimage, 47 148
–156
(2009). https://doi.org/10.1016/j.neuroimage.2009.03.058 1053-8119 Google Scholar
|