|Home | About | Journals | Submit | Contact Us | Français|
Low frequency oscillations (LFO) of cerebral vessels are believed to reflect cerebral autoregulation. We investigated day-to-day and hemispheric variations in 0·1 Hz LFO with near infrared spectroscopy (NIRS) and transcranial Doppler (TCD) to determine phase shift and gain of oxygenated haemoglobin (oxy-Hb) and the velocity of the middle cerebral artery (Vmca) to the arterial blood pressure (ABP). The direct left–right phase shifts of oxyHb and Vmca were also assessed. We examined 44 healthy volunteers by simultaneous recordings of ABP, oxyHb and Vmca during spontaneous and paced breathing at 6 breaths per minute on two separate days.
The variation between hemispheres had a prediction interval (PI) of ±39° for ABP–oxyHb phase shift and ±69% for gain. ABP–Vmca showed ±57° PI phase shift and ±158% PI for gain. The variation from day to day showed ±61° PI for ABP–oxyHb phase shift and ±297% PI for gain. ABP–Vmca showed ±45° PI phase shift and ±166% PI for gain. We found a linear relation between phase shift of oxyHb and Vmca at paced breathing (P = 0·0005), but not at rest (P = 0·235).
Our results show that LFO phase shift ABP–oxyHb may be used as a robust measurement of differences in autoregulation between hemispheres and over time. In addition, we found a strong relation between oxyHb and Vmca during paced breathing. Gain showed too large variation for clinical use, as the SD was up to 100-fold of mean values.
Near infrared spectroscopy (NIRS) is a noninvasive optical method to investigate regional changes in oxygenated (oxyHb) and deoxygenated haemoglobin (deoxyHb) in the outermost layers of the cerebral cortex . Among the advantages of NIRS are the high temporal resolution and the lack of discomfort during recording, which ensure high compliance for both patients and healthy volunteers. Few previous studies have used NIRS to investigate spontaneous 0·1 Hz low frequency oscillations (LFO) of cortical cerebral vessels [2–4], whereas most studies used transcranial Doppler (TCD) to examine LFO of middle cerebral artery velocity (Vmca) [5–10]. These studies analysed LFO and examined the time relation (phase shift) and the amplitude ratio (gain) between corresponding oscillations in arterial blood pressure (ABP) and cerebral vessels. The phase shift at 0·1 Hz is suggested to reflect the cerebral blood flow autoregulation  and is known to be altered in vascular disorders such as ischaemic stroke [11–13] and carotid artery disease [6,14]. We aimed to describe the inter-hemisphere and day-to-day variability of 0·1 Hz LFO phase shift and gain on oxyHb and Vmca measurements to the ABP in healthy volunteers. In addition, we examined the linear relation between oxyHb and Vmca as well as the possible variation caused by age and gender. We also investigated the direct phase shift and gain variation between left–right hemispheres for oxyHb and Vmca values.
We recruited 44 healthy volunteers with no history of cardio-vascular or neurological disease. Subjects using daily medication of any kind were excluded, except oral contraceptives, and two subjects using a cholesterol-lowering drug (simvastatine). None of the subjects had a history of alcohol or substance abuse, and none of the subjects had a body mass index >30. We recruited five men and five women in three age groups of 27–34, 37–44 and 47–54 years and a fourth group of six men and eight women aged >57 years (range, 57–69; mean, 63 years). Four subjects were left-handed (9%).
The Biomedical Research Ethics committee in the Capital Region of Denmark approved the study (H-B-2008-039). All subjects gave informed consent to participate in the study, which was undertaken in accordance with the Helsinki Declaration of 1964, as revised in Edinburgh in 2000.
Near infrared spectroscopy measurements were performed using continuous wave NIRS (NIRS2; TechEn Inc, Milford, MA, USA). The NIRS optodes were placed bilateral on the forehead with one light source at two wavelengths (690 and 830 nm) and two detectors on each side, avoiding the midline sinus. The distance between source and detectors was 3 cm, with the detectors lateral to the source. Thus, the detectors measured the frontal cortex region in the territory supplied by the middle cerebral artery (MCA). NIRS recordings were acquired at a sampling rate of 200 Hz.
The Vmca was recorded on both sides using a TCD monitor (Multidop-X4, DWL). The ultrasound probes were fixed using a head band. Identification of the middle cerebral artery was performed as previously described . All recordings were performed by the same trained physician (DP).
Continuous noninvasive ABP recording was achieved via a finger plethysmograph (CNAP500) using the subjects’ left hand, positioned at heart level. End-tidal CO2 partial pressure (PetCO2) was measured continuously using an open mask that caused no respiratory resistance (ProPac Encore). Data from TCD and finger plethysmograph were recorded synchronously with the NIRS data for subsequent offline analysis.
Each subject was examined on two different days within 1 week (range, 1–7 days; mean, 3·0 days), and both examinations were performed at the same time of the day. Before the experiments, each subject underwent a general physical and neurological examination. Coffee, tea or tobacco was not allowed 8 h before the investigations. All experiments were performed with the subjects placed in a supine position in a silent room with a constant temperature and the light dimmed. After 15 min of rest in the supine position, data acquisition was started with a 10-min trial of spontaneous breathing. This was followed by a 5-min trial with paced breathing at a rate of six respiration cycles per minute (0·1 Hz).
All signals were analysed in successive 50% overlapping time segments of 100 s. The NIRS light intensity at 690 and 830 nm were first converted to time series of variations in oxyHb concentrations using the modified Beer–Lambert law. On each time segment, we then computed the power spectra of all signals (ABP, oxyHb, Vmca), the coherence spectra for ABP–oxyHb and ABP–Vmca and the complex transfer function for ABP–oxyHb and ABP–Vmca, using Matlab (The MathWorks, Inc., Natick, MA, USA) functions pwelch, mscohere and tfestimate, respectively. All three functions are based on the Welch's averaged periodogram method, and we used 50% overlap 50-s Hamming windowing of the data in each segment. For two signals x and y, the Matlab function tfestimate computes the complex transfer function TFx-y as the quotient of the cross-power spectral density of the two signals Pxy and the power spectral density of the input signal Pxx. The magnitude-squared coherence Cxy between the signal x and y is obtained with mscohere as Cxy = (abs(Pxy))2/(Pxx . Pyy). The phase shift and gain between ABP oscillations and oxyHb and Vmca were obtained as the phase and absolute values of the complex transfer function, respectively. The LFO frequency was selected as the frequency with maximal ABP oscillation power in the 0·09–0·11 Hz range (Fig. 1). In this manner, for each run, we obtained 5 (5-min run) or 11 (10-min run) overlapping segments each characterized by a coherence value, a phase shift, and a gain for ABP–oxyHb and ABP–Vmca. We selected the segments for which the coherence was above an arbitrary threshold of 0·7 and averaged the gain (linear average) and the phase shift (circular average) over these selected segments. If the total number of segments with a coherence above 0·7 was smaller than 3 for a specific run, we discarded that run. Each NIRS channel was investigated independently and averaged subsequently.
Data are presented as mean values with standard deviation (SD), confidence interval (CI) or prediction interval (PI) which is ± 2 SD.
The primary end-points were as follows:
Secondary end-points were as follows:
For calculation of PI for ABP–oxyHb and ABP–Vmca and the relation between ABP–oxyHb and ABP–Vmca, we used a mixed model, which takes into account the correlations between sides and days of each subject. PI was chosen to establish a normal range that may be used in future studies examining cerebrovascular diseases.
Gain was measured in arbitrary unit (a.u.). Owing to nonnormal distribution of gain data and increasing variability with increasing level of gain for both ABP–oxyHb and ABP–Vmca, we applied a logarithmic transformation to the gain data. As a subtraction of logarithmic values results in a ratio after back-transformation, all gain data are presented in per cent for each recording, except mean values which are presented as a geometric mean that can be interpreted much like a median. All analyses were performed with SAS 9.2 (SAS Institute Inc., Cary, NC, USA). Five per cent (P = 0·05) was chosen as the level of significance.
The final data analysis comprised 44 healthy subjects. In seven subjects, the Vmca recordings on one or both sides were missing because of an insufficient TCD signal. The recordings at rest were missing in three subjects on 1 day because of technical problems. The percentage of recordings excluded because of low coherence in three or more windows were phase shift rest oxyHb 24% and Vmca 10%; phase shift 6 bpm oxyHb 16% and Vmca 15%; gain rest oxyHb 24% and Vmca 17%; gain 6 bpm oxyHb 17% and Vmca 15%. DeoxyHb showed very low coherence and was not analysed further. Mean arterial blood pressure (MAP) in the different age groups were 27–34 years: 81 mmHg, 37–44 years: 82 mmHg, 47–54 years: 90 mmHg and >57 years: 91 mmHg. There were no significant differences between groups (P = 0.072).
At rest, Vmca oscillations were 57° (PI – 137; 24), 1·6 s, ahead of ABP. The oxyHb oscillations were 15° (PI – 45; 75), 0·4 s, after ABP, thus 72° and 2 s after Vmca (Table 1). The direct relation for oxyHb left was 4° (PI – 29; 21), 0·1 s, after right, and for Vmca left was 5° (PI – 55; 65), 0·1 s, ahead of right (Table 2).
At rest, gain was 0·242 a.u. (PI 0·003–2·146) for ABP–oxyHb and 1·346 a.u. (PI 0·309–5·859) for ABP–Vmca. Gain was 0·950 a.u. (PI 0·544–1·356) for oxyHb left–right and 0·947 a.u. (PI 0·090–1·804) for Vmca left–right. PetCO2 did not change (32·6–32·2 mmHg, P = 0·197) at rest, whereas during paced breathing PetCO2 decreased (33·6–31·8 mmHg, P < 0·001).We found no significant differences between age groups and gender in any variables (Table 3).
There were no significant differences between hemispheres in phase shift and gain for ABP–oxyHb or ABP–Vmca. The mean differences in phase shift between the two hemispheres at rest ranged from 3 to 6° (PI 39) for ABP–oxyHb and from 0 to 2° (PI 57) for ABP–Vmca. Differences in gain at rest were 6% (PI 69%) in ABP–oxyHb and ranged from 2% to 8% (PI 158%) for ABP–Vmca (Table 4).
We found no significant difference between day 1 and day 2 in phase shift and gain for ABP–oxyHb and ABP–Vmca. The mean differences in phase shifts from day to day at rest ranged from 2 to 5° (PI 61) for ABP–oxyHb and from 2 to 10° (PI 45) for ABP–Vmca (Table 4). Differences in gain at rest ranged from 2% to 11% (PI 297%) in ABP–oxyHb and from 4% to 8% (PI 166%) for ABP–Vmca. For the direct left–right phase shift and gain, we found no significant differences for either oxyHb or Vmca (Table 2).
We investigated the linear relation between ABP–oxyHb and ABP–Vmca and found a positive significant linear relation during paced breathing for both phase shift (P = 0·0005, slope = 0·35, Fig. 2) and gain (P = 0·0025, slope = 0·36). We found no linear relation at rest for phase shift (P = 0·235), see Fig. 3, and gain (P = 0·131).
For phase shift at 6 bpm, Vmca was ahead of ABP with 54° (PI – 147; 38), 1·5 s, and oxyHb was 7° (PI – 60; 74), 0·2 s after ABP, thus a difference of 60°, 1·7 s. Gain during 6 bpm was 0·294 a.u. (PI 0·039–2·208) for ABP–oxyHb and 0·942 a.u. (PI 0·207–4·278) for ABP–Vmca (Table 5).
In the present study, we have showed the variation of LFO phase shift and gain for ABP–oxyHb and ABP–Vmca between hemispheres and days. In addition, we found a positive linear relation for phase shift and gain between oxyHb and ABP–Vmca during paced breathing, but not at rest. Gain showed too large variation for clinical use, as the SD was up to 100-fold of the mean values.
Obrig et al.  reported an ABP–oxyHb phase shift of zero in healthy volunteers (n = 3), and Payne et al.  using a larger source–detector separation of 5 cm and subjects in a sitting position, reported 7° ± 4° (n = 9), which is in agreement with the phase shift value of 15 found in the present study (n = 44). Whether the differences in examination protocol of the two studies have any influence on the results remain unknown. In the present study, the Vmca phase shift value at rest was 57°, which is similar to previous studies [6,9,16].
The variation between hemispheres is of great interest as cerebrovascular diseases, such as stroke and carotid artery disease, often present with unilateral changes. Previous TCD studies have shown that there is a difference in ABP–Vmca phase shift and gain between the affected and the contralateral side in minor stroke [5,11] and decreasing phase shift over time on the contralateral side after major stroke. Reinhard et al.  reported a phase shift difference of ABP–Vmca between hemispheres in healthy controls to be ± 1·5° (SD 7·7), which is in agreement with the present study, even though we found a higher variance (PI 57). Both studies used identical TCD set-up, except that Reinhard et al.  investigated the subjects with the upper body slightly elevated, which could cause changes in sympathetic tone affecting the LFO. For gain data, Reinhard et al.  reported the difference between hemispheres to be approximately ± 40 (SD 140), but we presented gain data in arbitrary units, and therefore, the results between the two studies cannot be directly compared. Furthermore, we found that variability increased as gain increased, which made it necessary to log-transform the gain data. This suggests that LFO data should be carefully examined before postprocessing, as transformation of data might become necessary.
Few TCD studies have investigated variation over time in different groups of patients to detect changes in the disease [5,7,18,19]. While a day-to-day variation of 26% on absolute Vmca values previously has been reported , we found no NIRS or TCD studies on LFO variation from day to day in healthy subjects at rest. Reproducibility between days of Vmca were examined by Reinhard et al.  in patients with unilateral carotid artery stenosis, which for both phase shift and gain showed significant correlation between days. However, one would expect that both gain and phase shift values would show significant day to day correlations. Thus, when reporting reproducibility, we believe that it is more useful to report the exact day-to-day variation than the correlation. Although we did not statistically compare inter-hemispheric and day-to-day variations between ABP–oxyHb and ABP–Vmca, the variations appear to be of the same magnitude, which suggests that NIRS recordings could also be used as a robust tool in assessing cerebrovascular autoregulation of cortical vessels between hemispheres and over time.
We investigated the phase shift relation between spontaneous oscillations of ABP–oxyHb and ABP–Vmca at rest and found no linear relation. The clear lack of linear relation might be due to local autoregulatory mechanisms, which could be an indicator of intact autoregulation in healthy volunteers. A previous study has reported a significant positive correlation between ABP–oxyHb and ABP–Vmca LFO at paced breathing . This is in agreement with the present study showing a significant positive linear relation only at paced breathing with a slope of 0·35.
Payne et al.  examined the effects of arterial blood gas levels on phase shift and found that the phase shift is strongly affected by the level of CO2. Fluctuations in CO2 also drive changes in the phase shift during paced breathing  and this explains the strong linear relationship found here between ABP–oxyHb and ABP–Vmca during paced breathing. The slope of the linear relation is <1 as the oxyHb and Vmca phase shift have different sensitivities to CO2 . Furthermore, Peng et al.  showed that, at rest, the power contained in the CO2 trace at 0·1 Hz is negligible, which explains why we did not find a linear relation between ABP–oxyHb and ABP–Vmca at rest.
Paced breathing induces oscillations in the vascular system by changes in the intrathoracic pressure, and paced breathing at 6 bpm is used to obtain strong regular oscillations for ABP, NIRS variables and Vmca [2,3]. Our phase shift values for ABP–oxyHb and ABP–Vmca using paced breathing at 6 bpm correspond with previous findings of ABP–oxyHb: +20°  and for ABP–Vmca: –65° . The phase shift of ABP–Vmca between hemispheres during paced breathing were previously examined in healthy controls in a carotid artery stenosis study  and showed a difference of 5·45° (SD 5·71), which is of the same magnitude as our results. Reinhard et al.  also investigated differences between spontaneous and paced breathing. The authors concluded that paced breathing was more reproducible, as paced breathing phase shift in Vmca showed a higher correlation between days than spontaneous breathing, which is in contrast to our findings.
Breathing amplitudes and phase probably vary across study days and subjects, which could explain the high variation for paced breathing. The CO2 drive controlling the phase shift during paced breathing should be compensated for . Furthermore, we found a significant decrease in PetCO2 during paced breathing despite attempting to keep the PetCO2 stable. Several healthy volunteers were unable to perform paced breathing for 5 min, which would likely be an even larger problem in patients with cerebrovascular or other diseases.
A previous study investigating age and LFOs in oxyHb showed that LFO amplitude declined with age , whereas another group investigating both oxyHb and Vmca found no differences for phase shift related to age . Thomsen et al.  reported that Vmca values were higher in women compared to men, and velocity was shown to decrease with age. Based on these finding, it could be suspected that there might be difference in oscillations between both gender and across age groups. However, we found no statistical differences between male and female volunteers or between the four age groups. As our healthy subjects did not undergo further imaging evaluation, we cannot rule out the existence of underlying subclinical cerebrovascular disease, which could have affected our findings. It has been suggested that the age changes may be related to a decline in spontaneous activity in microvascular smooth muscle cells and stiffening of the vessels . However, we did not examine the amplitude directly, but used gain, which is the amplitude ratio between oscillations in oxyHb and Vmca and corresponding oscillations in ABP. ABP oscillations could be affected by similar changes and therefore could leave the gain unaffected with no change with increasing age.
In conclusion, the present study has established the variation between hemispheres and over time for ABP–oxyHb and ABP–Vmca LFO phase shift, and the results indicate that phase shift for ABP–oxyHb in LFOs, as well as ABP–Vmca, can be used as a robust measurement of autoregulation status in patients with localized unilateral cerebrovascular disease. Thus, using bilateral NIRS and TCD technique may be a noninvasive tool to assess which vascular territories are affected by cerebrovascular disease and whether there might be a change over time in cerebral autoregulation. We also found a linear relation between LFO ABP–oxyHb and ABP-Vmca during paced breathings, but not at rest. NIRS and TCD measure different vascular compartments and therefore may provide complimentary information of cerebral autoregulation.
This study was supported by the Lundbeck Foundation via the Lundbeck Foundation Center for Neurovascular Signaling (LUCENS).
Dorte Phillip contributed to design of the study, acquisition of data, analysis of data and writing the draft of the paper. Henrik Schytz contributed to design of the study, analysis/interpretation of data and writing the paper. Juliette Selb contributed to analysis and interpretation of data and writing the paper. Stephen Payne contributed to interpretation of data and writing the paper. Helle K Iversen contributed to design of the study, interpretation of data and writing the paper. Lene T Skovgaard contributed to statistical analysis of data and writing the statistics. David A Boas contributed to analysis and interpretation of data and writing the paper. Messoud Ashina, the study promoter, contributed to design of the study, interpretation of data and writing the paper.
Conflict of interest
There was no conflict of interest in relation to this study