 Research
 Open Access
 Published:
Cognitive tasks during walking affect cerebral blood flow signal features in middle cerebral arteries and their correlation to gait characteristics
Behavioral and Brain Functions volume 11, Article number: 29 (2015)
Abstract
Gait is a complex process involving both cognitive and sensory ability and is strongly impacted by the environment. In this paper, we propose to study of the impact of a cognitive task during gait on the cerebral blood flow velocity, the blood flow signal features and the correlation of gait and blood flow features through a dual task methodology. Both cerebral blood flow velocity and gait characteristics of eleven participants with no history of brain or gait conditions were recorded using transcranial Doppler on midcerebral artery while on a treadmill. The cognitive task was induced by a backward counting starting from 10,000 with decrement of 7. Central blood flow velocity raw and envelope features were extracted in both time, frequency and timescale domain; informationtheoretic metrics were also extracted and statistical significances were inspected. A similar feature extraction was performed on the stride interval signal. Statistical differences between the cognitive and baseline trials, between the left and right midcerebral arteries signals and the impact of the antropometric variables where studied using linear mixed models. No statistical differences were found between the left and right midcerebral arteries flows or the baseline and cognitive state gait features, while statistical differences for specific features were measured between cognitive and baseline states. These statistical differences found between the baseline and cognitive states show that cognitive process has an impact on the cerebral activity during walking. The state was found to have an impact on the correlation between the gait and blood flow features.
Background
Walking is a complex sensorycognitive interaction which has various demands depending on the environment [1–3]. Initially, cognition and motor control of gait have been believed to be two completely autonomous processes, with walking regarded as an automatic motor function, independent of any cognitive tasks [4, 5]. Walking was thought to be automatic as it is generated by spinal cord oscillating circuits, or at best, locomotor centers in the brain stem but without cortical input under usual conditions. If that were true, there should be little if any interference of cognitive functioning with walking—but recent research shows that a cognitive load has an effect on gait [6–8]. In fact among people with pathological conditions, the dualtask methodology was mainly used to observe the effect of concurrent stimuli while checking the gait state [9]. This technique may underline a cognitivemotor interference that indicates a conflict between concurrent tasks (i.e. a motor and a cognitive challenges), as there may be a deterioration of one or both of the tasks [9, 10].
In general, a dual task (i.e. a cognitive stimulus/walking) causes concurrent demands for attention/cognitive resources. Hence, it generates a cognitivemotor interference that implies tripping, falling, physical instability and/or a decline in performance [2, 4, 11]. These falls occur mainly in elderly subjects and have a strong impact on the health of elderly subjects [6]. These events also have an impact on the health costs [12] and the risk of falling is associated with gait disruptions [13]. Prior studies have illustrated that individuals with neurodegenerative disorders, those recovering from stroke or elderly adults, are more prone to the effects of cognitive load on gait [7, 9, 14]. On the other hand, just a few investigations have investigated the cerebral repercussions of a dualtask in the case of ablebodied young participants [15]. Moreover, there has been limited research on the effects of cognitive load on the cerebrovascular system during gait.
Thus, there is a growing interest in clarifying the correlation between motor control and cognition. Brain imaging methods have revealed activation of various cerebral regions associated with higher cognitive functions during walking (i.e. the dorsolateral prefrontal cortex and anterior cingulate cortex [16–19]).
However, because of the lack of portability of these devices and the fact that they require the subject to lie down during acquisition, the gait is usually modeled as a series of feet tapping in the fMRI case [20, 21] and it is studied by the injection of a radiotracer, followed by the subject’s walk, and a scan for the PET/SPECT case [22, 23]. These devices, because of their lack of portability and the strong constraints associated with their use, are not optimal for performing a real time study of the brain during walking. These drawbacks led us to chose CBFV recordings using transcranial Doppler to assess brain function during walking.
Since it has been proven that CBFV modifications and neural activity can be correlated [24–28], we use transcranial Doppler recordings to monitor the hemodynamic activities of the main cerebral arteries [29] in order to study a motorcognition interaction. This noninvasive ultrasound diagnostic tool first introduced by Aaslid et al. measures the cerebral blood flow velocities (CBFV) [30]. Readings are taken with ultrasonic transducers, placed bilaterally in the transtemporal window of the skull of one participant, which allows the monitoring of the left and right side of the circle of Willis’ cerebral arteries [31–34]. Most studies focused on activities of the middle cerebral artery (MCA) given that the MCA carries more than 80 % of blood to the brain [35]. Moreover, previous papers highlighted brain perfusion changes during neural cognitive challenges [36–38]. In concern to physical performance, global cerebral blood flow (CBF) appears to be increased, unchanged or decreased during stimuli [39–41]. Nevertheless, a regional increase of CBF was noticed during physical exercise [42].
The current study was focused on studying the associations between the cerebrovascular system and gait and understanding the repercussions of cognitive load during walking on the system. Specially, we examined CBFV signals in MCA. Our major contributions are the statistical differences found in the central blood flow velocity signals between the baseline and cognitive states during walking. We examined both raw signals and the envelope signals, which are maximum peak velocity outcomes. In addition, baseline is actually walking—and on a treadmill which is a dedicated stepping pattern facilitator, walking could be even more automatic than usual—yet the blood flow changed with added cognitive processing; and the changes in CBF consistency mirrored in part the changes in gait consistency—at least for the performance of some—illustrating the potential interference of cognitive functioning and walking.
Methods
Data acquisition
For this preliminary data collection, eleven participants were recruited (4 females and 7 males, ages ranging from 19 to 23 years). None of these subjects had history of concussion, heart murmurs, migraines, strokes or other brain and gait conditions. Participants were asked to signed the University of Pittsburgh Institutional Review Board approved consent form and the procedure was explained prior to the beginning of the experiment and data collection.
The participants were briefed to walk at a pace of 2 mph on the Noraxon MR3 treadmill and to remain thoughtfree during this 6 min baseline testing period. Upon the completion of this trial, the participants were asked to count backwards from 10,000 in decrements of 7 while walking for 6 min for another trial referred as cognitive. Approximately 1 week later, each participant repeated the experiment.
To observe stride time, Noraxon SciFit treadmill using capacitive sensor technology to analyze individual foot pressure (ranging from 1 to \(120~\,N /cm ^2\)) was used at a sampling frequency of 100 Hz. The treadmill data was recorded and extracted using the manufacturer’s software.
A SONARA TCD System (Carefusion, San Diego, CA, USA) was used to measure bloodflow velocity in the midcerebral artery. Two 2 MHz transducers were used to gather simultaneous bilateral CBFV acquisitions from the leftMCA and rightMCA. The transducer were planted on transtemporal windows [43] to reach the MCA bloodflows. The position, angle and insonation depth of transducers were adjusted in order to get the correct MCA signal [34, 44]. Once these adjustments were performed, the transducers were fixed with a headset on both sides of the subject’s head. The data was extracted as audio files sampled at \(44.1\,~kHz\), representing the cerebral blood flow from the RMCA and LMCA. These signals were resampled at 8820 Hz (factor 5) to increase the feature extraction’s speed. The downsampled signals, referred as raw signals, are composed of multiple sinusoidal components due to the parabolic speed CBFV distribution [37]. In this study, these raw CBFV along with the envelope of the CBFV were collected. The raw signals are composed of the various velocities of blood particles in cerebral arteries and envelope signals constitute the maximal Doppler shift [37, 45].
Furthermore, the endtidal carbon dioxide \(ETCO _2\) (BCI Capnocheck Sleep Capnograph, Smiths Medical, Waukesha, Wisconsin, USA) was monitored along with respiration rate, electrocardiogram, head movement and skin conductance via a multisystem physiological data monitoring system (NexusX, Mindmedia, The Netherlands).
Analysis of stride interval time series
The mean, standard deviation, coefficient of variation (which is the ratio of the standard deviation to the mean) and the spectral exponent were estimated for stride interval time series extracted from both feet.
To determine the spectral exponents of the stride interval signals, the average wavelet coefficient (AWC) method was used. This method can be summarized as follows [46, 47]:

1.
The wavelet transform of the centred version of the signal, \(WV _x(\tau ,s)\), is computed using the Daubechies 12 wavelet.

2.
The mean average with respect to translation coefficient \(\tau\) of the magnitude of the wavelet transform is evaluated for each scale. This quantity is called:
$$\begin{aligned} WV _x(s) = \left\langle WV _x(\tau ,s) \right\rangle _\tau \end{aligned}$$(1) 
3.
The loglog plot of \(WV _x(s)\) as a function of the scale s is plotted. The linear regression of the resulting graph is computed. The slope of this line is \(H_{fBm} + 1/2\) [46].

4.
The spectral exponent \(\beta\) is evaluated using the expression:
$$\begin{aligned} \beta = 2 H_{fBm} + 1 \end{aligned}$$(2)
The number of scales were chosen with respect of the signal length, using the following expression:
Analysis of MCA signals
Statistical features
In this study, four distinct statistical features where evaluated. The standard deviation, the skewness, the kurtosis along with the crosscorrelation coefficient of right and left MCA signals were extracted and compared. The first three parameters characterize the shape of the signal’s distribution [48], while the last one characterizes the similarity between two signals.
The kurtosis of a distribution can be expressed as [48]:
where \(\mu\) is the mean of the signal and n is its length.
The expression of the skewness is [48]:
In this study, similarity between right MCA and left MCA signals was calculated as follows [49]:
where X and Y represent signals from the right and the left side of the MCA.
Informationtheoretic features
With regard to the informationtheoretic feature space, the Lempel–Ziv complexity (LZC), which represents the amount of new pattern formation in finite time sequences built upon the original signal (which can be interpreted as the randomness, the predictability and the regularity of a given discretetime signal) [50] and the entropy rate, which represents the statistic degree of recurrence of patterns in a stochastic process [51].
In order to obtain a finite time sequence to compute the LZC, the signal is divided into \(\alpha\) finite equal spaces using \(\alpha 1\) threshold values given as \(T_h  h \in \{1, \ldots , \alpha  1\}\) [52].
Then, portions of the quantized signal \(X_1^n=\{x_1,x_2,\ldots ,x_n\}\) are assembled to shape blocks such that [53]:
where the length of each block is given as \(L = j  l + 1\), and the length of the signal is n. The blocks are a time series of successive data.
Then, for each length L, an analyse of each block is performed the following way: a counter c(j) is defined to illustrate the amount of new pattern formation. If the sequence represented in a block has not appeared in a previous analysis, this counter is incremented by one.
Eventually, the LZC is computed as given herein bellow:
where c(n) is the value that the counter takes when the entire analysis is performed and \(\alpha\) represents the total number of quantized levels in the signal, chosen such as \(\alpha = 100\) in that study.
To compute the entropy rate \(\rho\) of a given stochastic process, the pattern distribution is first normalized to feature zeros mean and unit variance. This distribution is then quantized using 10 equally divided discrete amplitude levels, taking integer values between 0 and 9. The quantized distribution \(\tilde{X} = \{\tilde{x}_1, \tilde{x}_2,\ldots ,\tilde{x}_n\}\) is then decomposed into blocks comprised of successive samples, with length L such that \(10 \le L \le 30, \; L \in \mathbb {Z}^+\). Subsequently, the distribution made up of the different blocks was encoded into \(\Omega _L=\{ \omega _1, \omega _2,\ldots ,\omega _{nL+1}\}\) such as [54]:
with \(w_i\) varying from 0 to \(9 \times (1  10^L)/(110) = 10^L  1\).
The Shannon entropy \(S (L)\), which represents the degree of complexity of \(\Omega _L\), is defined by [54]:
where \(p_{\Omega _L}(j)\) is the probability of the value j in \(\Omega _L\), which is approximated by its sample frequency in this study.
The normalized conditional entropy is then given as [55]:
where \({\rm S} (1) \cdot {\rm pe} _\%(L)\) is a correction term defined by the multiplication of the percentage of patterns with length L arising only once in \(\Omega _L\), \({\rm pe} _\%(L)\), with \({\rm S} (1)\) the conditional entropy estimation of the process with unit length L, which is the Shannon entropy of white Gaussian noise process. This term corrects the underestimation of \({\rm S} (L){\rm S} (L1)\) for large lengths L [56]. Considering the opposite nature of the variation of the terms in the denominator (the first term decreases with L while the second term increases), the function \({\rm N} (L)\) exhibits a minimum \(\min _{L}\left[ N(L)\right]\). This minimum is the best estimation of the normalized conditional entropy, and it can be seen as an indicator of complexity of the process. Conversely, the complement to this indicator, given as \(\rho = 1  \min _{L}\left[ N(L)\right]\), is an index of the regularity of the stochastic process, ranging between 0 and 1 [55].
For comparison between two probability density functions purposes, one can use the crossentropy rate. This index quantifies the amount of mutual information between two given distributions, and aims to predict the data of a considered signal using the previous and current information found in another signal. The two distributions X and Y were normalized and quantized in a similar fashion than the one used to compute the entropy rate, giving as a result the two quantized distribution \(\tilde{X} = \{\tilde{x}_1, \tilde{x}_2,\ldots ,\tilde{x}_n\}\) and \(\tilde{Y} = \{\tilde{y}_1, \tilde{y}_2,\ldots ,\tilde{y}_n\}\).
Eventually, the crossentropy rate \(\Omega _L^{XY}\), which is the information amount that can be found in one of the samples of the quantized process \(\tilde{X}\) when a pattern of \(L1\) samples of the quantized signal \(\tilde{Y}\) is assumed, was encoded using the following code [55]:
where \({\rm S} _X(L)\), \({\rm S} _Y(L)\) and \({\rm S} _{XY}\) are the Shannon entropies of the distributions X, Y and \(\Omega _L^{XY}\).
The normalized crossentropy was then computed as:
where \({\rm pe} _{XY}(L)\) is the percentage of arrangements of length L that were present only once in \(\Omega _L^{XY}\) and \({\rm S} _X(1) \cdot {\rm pe} _{XY}(L)\) is a corrective term added for the same reasons as the one stated herein above. Similarly to the previously given algorithm, \({\rm S} _X(1)\) is the conditional entropy estimation of the stochastic process X for a unit length. The function \({\rm NC} _{XY}(L)\) features a minimum \(\min _L\left[ {\rm NC} _{XY}(L), {\rm NC} _{YX}(L) \right]\). The synchronization index is then defined by \(\Lambda _{XY}=1 \min _L\left[ {\rm NC} _{XY}(L), {\rm NC} _{YX}(L) \right]\), and it ranges from 0 (when X and Y are independent) to 1 (when X and Y are synchronized stochastic processes).
Frequency features
Features were also extracted with regard to the frequency domain: the peak frequency, spectral centroid and bandwidth were computed as characterization indexes for the extracted from the TCD signals [45].
The peak frequency, given as \(f_p\), is the frequency location where the largest spectral power can be found, and it is thus computed as:
where \(F_X(f)\) is the Fourier transform of the signal X and \(f_{\max }\) the spectrum’s maximum frequency, equal to the sampling frequency divided by two.
The spectral centroid, \(f_c\), can be interpreted as the spectrum’s center of mass and is calculated as follows [57]:
The bandwidth BW, representing the dispersion of the spectrum, was computed as [45]:
Timefrequency features
To extract features related to the timescale domain, the signal was decomposed into 10 levels \(W = \left\{ a_{10}, d_{10}, d_{9}, \ldots , d_{1} \right\}\) (where \(a_{10}\) is the approximation coefficient and \(d_{k}\) represents detail coefficient at the kth level [58]) using a discrete wavelet transform approach and the Meyer wavelet as a mother wavelet. The relative wavelet energy from both the approximation coefficient and detail coefficients are then computed as [59]:
where \(\Vert \cdot \Vert\) is the Euclidian norm, defined as \(\Vert \varvec{x} = [x_1,\ldots , x_n] \Vert = \left( {\sum x_i^2}\right) ^{1/2}\).
These metrics describe the relative energies repartition within distinct frequency bands based upon the ratio of the kth level of decomposition to the total energy of the signal.
The wavelet entropy \(\Omega\) is an indicator of the order of the decomposed signal [59]. It is a representation of the clustering of the wavelet energies comprised in the different decomposition levels.
where \(\Xi _{\cdot }\) are the relative energies defined herein above.
Statistical test
SAS version 9.3 (SAS Institute, Inc., Cary, North Carolina) was used for all statistical analyses. Our overall strategy was to keep the two trials separate and fit a sufficiently complex model in order to account for betweentrial correlation to maximize the use of information and statistical power. For simpler methods such as correlation analyses, we averaged participant measures across the two trials. To examine significance of systematic differences in measures between the two trials we fitted a series of linear mixed models using the SAS MIXED procedure with each gait/blood flow measure as the dependent variable; trial as the factor of interest; and a subject random effect to account for multiple correlated measurements from the same participant. We computed baseline condition vs cognitive task difference for each measure, and fitted a similar mixed model with the said differences as dependent variables but in an interceptonly model to obtain statistical significance for the betweencondition difference under two correlated trials. To examine whether participant characteristics (age, gender, height, weight, BMI), we fit another series of models with each gait/blood flow measure under each task condition as the dependent variable; each participant characteristic as the fixed effect of interest; and a participant random effect to account for multiple trials. To examine correlation coefficients across measures, we first averaged all measures across the two trials. Next, we used Pearson correlation coefficients between gait measures and blood flow measures; as well as between left and right sides for the blood flow measures. Additionally, we obtained partial correlations adjusted for participant characteristics. We are aware of the large number of measures and statistical tests performed in our analysis, and thus use false discovery rate (FDR) methodology for adjusting pvalues for multiplicity [60].
Results
The endtidal carbon dioxide level does not influence the mean diameter of the middle cerebral arteries [61]. Consequently, it is not taken into consideration. The results are presented in tables in the form of \(({\rm mean} ~\pm ~ {\rm standard deviation} )\) where the baseline period is pointed out by a “B” and the cognitive task is indicated by a “C”. RMCA indicates the right MCA, while LMCA indicates the left MCA. The testretest were not found to introduce statistical differences if the FDR adjusted pvalues are considered.
Stride interval features
Feature extracted from the subjects’ stride intervals are provided in Table 1. No statistical differences between the baseline (denoted as “B”) and cognitive states (denoted as “C”) were found if the FDR adjusted p values are considered.
MCA signal features
Time features
Table 2 summarizes timedomain feature values for the raw and the envelope signals. For every metrics considered in that table, no statistical differences between the left side and the right side were observed. However, statistical differences between cognitive state and baseline state were observed for the standard deviation of the envelope signal (\(p \approx 0.018\)) for both RMCA and LMCA, and for both the envelope and raw signals’ RMCA and LMCA crosscorrelation (\(p \approx 0.018\)).
Moreover, the inspection of the mean peak blood flow velocity revealed an increase during the cognitive task (\(57.8~\pm ~9.0 {\rm cm} ~s ^{1}\) in the baseline trial vs \(59.5 ~\pm ~ 10.3 {\rm cm} ~s ^{1}\) in the cognitive trial for the RMCA and \(59.1~\pm ~8.4 {\rm cm} ~s ^{1}\) in the baseline trial vs \(59.3~\pm ~10.1 {\rm cm} ~s ^{1}\) in the cognitive trial for the LMCA).
Additionally, a similar increase in the standard deviation of the envelope signal (the increase if of 1.25 \(cm ~s ^{1}\) for the RMCA and of 1.41 \(cm ~s ^{1}\) for the LMCA) can be observed. This increase in the standard deviation of the blood flow velocity is characteristic of an increase in the blood flow: because the higher standard deviation characterizes a wider blood flow velocity dispersion, this means that there is an increase in the range covered by the minimal and maximal blood flow velocity, which is characteristic of a global increase of the blood flow velocity.
Informationtheoretic features
Informationtheoretic features from both CBFV raw and envelope signals are presented in Table 3. The LZC and synchronization index of the envelope signals exhibited significant statistical differences between the two states (\(p < 0.02\)).
Frequency features
Table 4 summarizes the frequency characteristics of raw and envelope signals. These metrics do not display any significant statistical differences between the left and right MCA raw and envelope signals. The spectral centroid of the RMCA raw signal was found to feature statistical differences (\(p = 0.038\)), while the bandwidth of the raw CBFV shows statistical differences for the LMCA, with \(p < 0.01\).
Timefrequency features
Table 5 shows the wavelet entropy values for raw and envelope signals, while Figure 1 displays the wavelet decomposition of the raw CBFV signal. The envelope signal’s wavelet decomposition is not displayed because almost the entire relative wavelet energy can be found in \(a_{10}\) (99 %). When the wavelet entropies are considered, the RMCA and LMCA CBFV wavelet entropies feature statistical differences between baseline and cognitive states, with respectively \(p = 5.4~\times ~10^{3}\) and \(p = 1.1~\times ~10^{2}\).
Figure 1 exhibits that the majority of wavelet energy is massed in high order detail coefficients (mainly for \(7 \le k \le 10,~k \in \mathbb {Z}^+\)). The coefficient \(d_{10}\) accounts for the highest percentage of wavelet energy (between 70 % for the RMCA in the baseline state and 58 % for the LMCA in the cognitive state). The energy contained is the \(d_{10}\) band is higher in the baseline state than in the cognitive state. A similar observation can be drawn with energy repartition in the \(d_{10}\) of the RMCA being higher than the energy of the LMCA wavelet decomposition. The \(d_9\) band behaves differently, and no particular trend emerges in this detail band. The behavior of the \(d_8\) band is the inverse of the \(d_{10}\) band: indeed, the energy proportion of the \(d_8\) band increases from the baseline state to the cognitive state and from the RMCA to the LMCA. Statistical differences between the baseline and cognitive state can be found in the bands \(d_{10}\) and \(d_8\) for the raw RMCA and LMCA signals (\(p~\le ~0.03\)).
Correlation between left and right MCA blood flow
Statistical significance for the correlation between the left and right MCA CBFV are only found for the wavelet decomposition of the raw CBFV signals: the levels \(d_5\) to \(d_7\) exhibit a correlation coefficient greater than 0.8 with \(p \le 0.005\). The levels \(d_2\) and \(d_1\) feature a correlation coefficient greater than 0.9 with \(p \ll 0.001\).
Impact of anthropometric variables on features
The anthropometric variables were not found to have an impact on any of the features.
Correlation between gait and MCA features
The Pearson correlation coefficient \(\rho\) were thresholded using the following rule:
Figure 2 displays the correlation of the gait features and the CBFV features. A darkgray rectangle represents an absolute correlation of at least 0.7 in the baseline state, while a lightgray rectangle represents a correlation in the cognitive state. Some features exhibit correlations in both states (e.g. LMCA raw signal synchronization index is always correlated to the right foot standard deviation or RMCA raw signal kurtosis is correlated to both the right foot standard deviation and coefficient of variation), while some features’ correlation changes from one state to another (e.g. the correlation between the RMCA Lempel–Zif complexity and the right foot mean stride interval observed in the baseline state is not present in the cognitive state, or the correlation between the left foot standard deviation and the entropy rate of the LMCA raw signal observed in the cognitive state is not present in the baseline state). From this figure, it is clear that with the cognitive processing added task there was a greater number of correlated gait and blood flow velocity features: in the baseline test, 9.1 % of the potential relations between gait and cerebral blood flow were correlated, while in the cognitive processing added task condition, 14.5 % of the potential gait and blood flow measures were related.
More specifically, all the striderelated features except the right foot spectral exponent coefficient are found to be correlated with at least one feature extracted from the cerebral blood flow signals. Indeed, four blood flow related features are found to be correlated to both the left and right feet mean stride interval (peak frequency of the RMCA raw signal, standard deviation of the LMCA raw signal, both the skewness and kurtosis of the RMCA envelope signals). The standard deviation and coefficient of variation of the right foot stride interval both are correlated with 10 of the features of the CBFV (with 9 of the correlated being the same for the left and right foot). The same stride related features when the left foot is considered are correlated with 11 and 10 features of the CBFV, respectively. Lastly, the spectral exponent coefficient for the left foot stride interval is correlated to two of the central blood flow velocity signals features: the RMCA raw signal peak frequency and the RMCA envelope signal wavelet entropy.
Discussion
Our major finding is that there are statistical differences in cerebral blood flow velocities in MCA between the baseline and cognitive states during walking. These differences were observed both in the raw signals and in the envelope signals (i.e., the maximal velocities) associated with cerebral blood flow velocities.
The statistical differences revealed in this study between blood flow during walking and blood flow during walking with an added cognitive processing task illustrate the suggested performance pattern seen in older adults with greater difficulty walking under challenging conditions such as walking and talking or walking and thinking. Under the challenge of the added cognitive task, it was more difficult to maintain usual gait performance. The difficulty manifested in gait as largely greater inconsistency in the pattern of strides than in the mean of the stride interval time series. This difficulty was also mirrored in the cerebral blood flow. Thus, under the cognitive processing task condition, the range of gait performance was expanded for some variables, despite no change in the group mean for the gait characteristics. The expanded range of gait performance may underlie the different and greater number of correlations between gait and blood flow in the walking plus cognitive processing condition. A cognitively challenging task during walking may alter the range of gait performance, which is also associated with greater cerebral blood flow or ‘work’ of the brain.
Among the stride measures, the standard deviation of the stride interval exhibited the strongest correlation with blood flow velocity features. This finding may indicate that while the cognitive processing did not disrupt gait described by the mean of the stride interval, the added cognition may have influenced the consistency of gait performance (eg, standard deviation). While all participants accomplished the task of walking and thinking on the treadmill with no change in the mean stride characteristics, the difference in the level of cognitive challenge for participants was apparent in the inconsistency of gait. Greater cognitive challenge for some was associated with greater blood flow. Moreover, the increase in the bandwidth of the raw signals denotes a change in individual behavior of the erythrocytes, as a higher bandwidth value denotes higher spectral spread, which is caused by more erratic flow of red blood cells. The lack of statistical differences between the first and second trials denotes that the extracted features tend to be trialindependent.
Limitations of the study
A potential limitation of this study comes from the fact that the participants were instructed to walk on a treadmill: indeed, changes in the gait characteristics are very likely concealed by the walking rhythm induced by the treadmill.
Conclusion
In this paper, the dual task methodology exhibited that the cerebral blood flow velocity signals in the baseline and cognitive states during walking are statistically different. These differences can be found in extracted features for both the raw and envelope of the LMCA and RMCA central blood flow velocity signals. While the changes in gait features were subtle, the different and more relations between gait and blood flow with the added cognitive processing during walking suggests cerebral blood flow velocity may represent the work of the brain when thinking and walking.
References
 1.
Sheridan P, Solomont J, Kowall N, Hausdorff J. Influence of executive function on locomotor function: divided attention increases gait variability in Alzheimer’s disease. J Am Geriatr Soc. 2003;51(11):1633–7.
 2.
AlYahyaa E, Dawesa H, Smithc L, Dennisa A, Howellsa K, Cockburna J. Cognitive motor interference while walking: a systematic review and metaanalysis. Neurosci Biobehav Rev. 2011;35(3):715–28.
 3.
Lajoie Y, Teasdale N, Bard C, Fleury M. Attentional demands for static and dynamic equilibrium. Exp Brain Res. 1993;97(1):139–44.
 4.
YogevSeligmann G, Hausdorff J, Giladi N. The role of executive function and attention in gait. Mov Disord. 2008;23(3):329–42.
 5.
Fraizer E, Mitra S. Methodological and interpretive issues in posturecognition dualtasking in upright stance. Gait Posture. 2008;27(2):271–9.
 6.
Giladi N. Gait and mental function: the interplay between walking, behavior and cognition. J Neural Transm. 2007;114(10):1241–2.
 7.
Eladio Martin E, Bajcsy R. Analysis of the effect of cognitive load on gait with offtheshelf accelerometers. In: COGNITIVE 2011: The Third International Conference on Advanced Cognitive Technologies and Applications; 2011. p. 1–6.
 8.
Allali G, Van Der Meulen M, Beauchet O, Rieger SW, Vuilleumier P, Assal F. The neural basis of agerelated changes in motor imagery of gait: an fMRI study. J Gerontol A Biol Sci Med Sci. 2013; 207:1389–98.
 9.
PlummerD’Amatoa P, Altmannb L, Saracinoc D, Foxd E, Behrmand A, Marsiskef M. Interactions between cognitive tasks and gait after stroke: a dual task study. Gait Posture. 2008;27(4):683–8.
 10.
Rao A, Jasim U, Gillman A, Louis E. Cognitive motor interference during dualtask gait in essential tremor. Gait Posture. 2013;38(3):403–9.
 11.
Hollman J, Kovash F, Kubik J, Linbo R. Agerelated differences in spatiotemporal markers of gait stability during dual task walking. Gait Posture. 2007;26(1):113–9.
 12.
Beauchamp MK, Sibley KM, Lakhani B, Romano J, Mathur S, Goldstein RS, Brooks D. Impairments in systems underlying control of balance in copdbalance in COPD. CHEST J. 2012;141(6):1496–503.
 13.
Axer H, Axer M, Sauer H, Witte OW, Hagemann G. Falls and gait disorders in geriatric neurology. Clin Neurol Neurosurg. 2010;112(4):265–74.
 14.
Ijmker T, Lamoth C. Gait and cognition: the relationship between gait stability and variability with executive function in persons with and without dementia. Gait Posture. 2012;35(1):126–31.
 15.
Woollacott M, ShumwayCook A. Attention and the control of posture and gait: a review of an emerging area of research. Gait Posture. 2002;16(1):1–14.
 16.
Yogev G, Giladi N, Peretz C, Springer S, Simon E, Hausdorff J. Dual tasking, gait rhythmicity, and Parkinson’s disease: which aspects of gait are attention demanding? Eur J Neurosci. 2005;22(5):1248–56.
 17.
Harada T, Miyai I, Suzuki M, Kubota K. Gait capacity affects cortical activation patterns related to speed control in the elderly. Exp Brain Res. 2009;193(3):445–54.
 18.
Suzuki M, Miyai I, Ono T, Oda I, Konishi I, Kochiyama T, Kubota K. Prefrontal and premotor cortices are involved in adapting walking and running speed on the treadmill: an optical imaging study. Neuroimage. 2004;23(3):1020–6.
 19.
Iseki K, Hanakawa T, Shinozaki J, Nankaku M, Fukuyama H. Neural mechanisms involved in mental imagery and observation of gait. Neuroimage. 2008;41(3):1021–31.
 20.
Mehta JP, Verber MD, Wieser JA, Schmit BD, SchindlerIvens SM. A novel technique for examining human brain activity associated with pedaling using fMRI. J Neurosci Methods. 2009;179(2):230–9.
 21.
Godde B, VoelckerRehage C. More automation and less cognitive control of imagined walking movements in highversus lowfit older adults. Front Aging Neurosci. 2010;2.
 22.
Fukuyama H, Ouchi Y, Matsuzaki S, Nagahama Y, Yamauchi H, Ogawa M, Kimura J, Shibasaki H. Brain functional activity during gait in normal subjects: a SPECT study. Neurosci Lett. 1997;228(3):183–6.
 23.
Higuchi M, Okamura N, Ishii K, Bereczki D, Sasaki H. 18FFDG PET mapping of regional brain activity in runners. J Sports Med Phys Fitness. 2001;41:11–7.
 24.
Hartje W, Ringelstein EB, Kistinger B, Fabianek D, Willmes K. Transcranial Doppler ultrasonic assessment of middle cerebral artery blood flow velocity changes during verbal and visuospatial cognitive tasks. Neuropsychologia. 1994;32(12):1443–52.
 25.
Droste DW, Harders AG, Rastogi E. A transcranial Doppler study of blood flow velocity in the middle cerebral arteries performed at rest and during mental activities. Stroke. 1989;20(8):1005–11.
 26.
BullaHellwig M, Vollmer J, Götzen A, Skreczek W, Hartje W. Hemispheric asymmetry of arterial blood flow velocity changes during verbal and visuospatial tasks. Neuropsychologia. 1996;34(10):987–91.
 27.
Kelley RE, Chang JY, Scheinman NJ, Levin BE, Duncan RC, Lee SC. Transcranial Doppler assessment of cerebral flow velocity during cognitive tasks. Stroke. 1992;23(1):9–14.
 28.
Vingerhoets G, Stroobant N. Lateralization of cerebral blood flow velocity changes during cognitive task: a simultaneous bilateral transcranial Doppler study. Stroke. 1999;30(10):2152–8.
 29.
Brass L, Pavlakis S, DeVivo D. Transcranial Doppler measurements of the middle cerebral artery. Effect of hematocrit. Stroke. 1988;19(12):1466–9.
 30.
Aaslid A, Markwalder TM, Nornes H. Noninvasive transcranial Doppler ultrasound recording of flow velocity in basal cerebral arteries. J Neurosurg. 1982;57(6):769–74.
 31.
Reid JM, Spencer MP. Ultrasonic Doppler technique for imaging blood vessels. Science. 1972;176(4040):1235–6.
 32.
Bishop CC, Powell S, Rutt D, Browse NL. Transcranial Doppler measurement of middle cerebral artery blood flow velocity: a validation study. Stroke. 1986;17(5):913–5.
 33.
Larsen F, Olsen K, Hansen B, Paulson O, Knudsen G. Transcranial Doppler is valid for determination of the lower limit of cerebral blood flow autoregulation. Stroke. 1994;25(10):1985–8.
 34.
White H, Venkatesh B. Applications of transcranial Doppler in the ICU: a review. Intensive Care Med. 2006;32(7):981–94.
 35.
Lindegaard K, Lundar T, Wiberg J, Sjøberg D, Aaslid R, Nornes H. Variations in middle cerebral artery blood flow investigated with noninvasive transcranial blood velocity measurements. Stroke. 1987;18(6):1025–30.
 36.
Fox PT, Raichle ME. Focal physiological uncoupling of cerebral blood flow and oxidative metabolism during somatosensory stimulation in human subjects. Proc Natl Acad Sci USA. 1986;83(4):1140–4.
 37.
Deppe M, Ringelstein EB, Knecht S. The investigation of functional brain lateralization by transcranial Doppler sonography. Neuroimage. 2004;21(3):1124–46.
 38.
Phillips A, Ainslie P, Krassioukov A, Warburton D. Regulation of cerebral blood flow after spinal cord injury. J Neurotrauma. 2013;30:1551–63.
 39.
Globus M, Melamed E, Keren A, Tzivoni D, Granot C, Lavy S, Stern S. Effect of exercise on cerebral circulation. J Cereb Blood Flow Metab. 1983;3(3):287–90.
 40.
Herholz K, Buskies W, Rist M, Pawlik G, Hollmann W, Heiss W. Regional cerebral blood flow in man at rest and during exercise. J Neurol. 1987;234(1):9–13.
 41.
Madsen P, Sperling B, Warming T, Schmidt J, Secher N, Wildschiodtz G, Holm S, Lassen N. Middle cerebral artery blood velocity and cerebral blood flow and \(o_2\) uptake during dynamic exercise. J Appl Physiol. 1993;74(1):245–50.
 42.
Jorgensen L, Perko G, Secher N. Regional cerebral artery mean flow velocity and blood flow during dynamic exercise in humans. J Appl Physiol. 1992;73(5):1825–30.
 43.
Schneider P, Rossman M, Bernstein E, Torem S, Ringelstein E, Otis S. Effect of internal carotid artery occlusion on intracranial hemodynamics. Transcranial Doppler evaluation and clinical correlation. Stroke. 1988;19(5):589–93.
 44.
Alexandrov A, Sloan M, Wong L, Douville C, Razumovsky A, Koroshetz W, Kaps M, Tegeler C. Practice standards for transcranial Doppler ultrasound: part I—test performance. J Neuroimaging. 2007;17(1):11–8.
 45.
Sejdić E, Kalika D, Czarnek N. An analysis of restingstate functional transcranial Doppler recordings from middle cerebral arteries. PLoS One. 2013;8(2):e554051–9.
 46.
Schaefer A, Brach JS, Perera S, Sejdić E. A comparative analysis of spectral exponent estimation techniques for \(1/f^\beta\) processes with applications to the analysis of stride interval time series. J Neurosci Methods. 2014;222:118–30.
 47.
Simonsen I, Hansen A, Nes OM. Determination of the Hurst exponent by use of wavelet transforms. Phys Rev E. 1998;58:2779–87.
 48.
Papoulis A. Probability, Random Variables, and Stochastic Processes. New York: WCB/McGrawHill; 1991.
 49.
Tiao G, Box G. Modeling multiple time series with applications. J Am Stat Assoc. 1981;76(376):802–16.
 50.
Ahmed S, Shahjahan M, Murase K. A Lempel–Ziv Complexitybased neural network pruning algorithm. Int J Neural Syst. 2011;21(5):427–41.
 51.
Pincus S, Gladstone I, Ehrenkranz R. A regularity statistic fore medical data analysis. J Clin Monit. 1991;7(4):335–45.
 52.
Hu J, Gao J, Principe J. Analysis of biomedical signals by the Lempel–Ziv Complexity: the effect of finite data size. IEEE Trans Biomed Eng. 2006;53(12):2606–9.
 53.
Lempel A, Ziv J. On the complexity of finite sequences. IEEE Trans Inf Theory. 1976;22(1):75–81.
 54.
Porta A, Guzzetti S, Montano N, Furlan R, Pagani M, Somers V. Entropy, entropy rate, and pattern classification as tools to typify complexity in short heart period variability series. IEEE Trans Inf Theory. 2011;48(11):1282–91.
 55.
Porta A, Guzzetti S, Montano N, Pagani M, Somers V, Malliani A, Baselli G, Cerutti S. Information domain analysis of cardiovascular variability signals: evaluation of regularity, synchronisation and coordination. Med Biol Eng Comput. 2000;38(2):180–8.
 56.
Porta A, Baselli G, Liberati D, Montano N, Cogliati C, GnecchiRuscone T, Malliani A, Cerutti S. Measuring regularity by means of a corrected conditional entropy in sympathetic outflow. Biol Cybern. 1998;78(1):71–8.
 57.
Vergara L, Gosalbéz J, Fuente J, Miralles R, Bosch I. Measurement of cement porosity by centroid frequency profiles of ultrasonic grain noise. Biomed Eng Online. 2004;84(12):2315–24.
 58.
Hilton M. Wavelet and wavelet packet compression of electrocardiograms. IEEE Trans Biomed Eng. 1997;44(5):394–402.
 59.
Rosso O, Blanco S, Yordanova J, Kolev V, Figliola A, Schürmann M, Basar E. Wavelet entropy: a new tool for analysis of short duration brain electrical signals. J Neurosci Methods. 2001;105(1):65–75.
 60.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B (Methodological). 1995;57:289–300.
 61.
Giller C, Bowman G, Dyer H, Mootz L, Krippner W. Cerebral arterial diameters during changes in blood pressure and carbon dioxide during craniotomy. Neurosurgery. 1993;32(5):737–42.
Authors’ contributions
AG wrote the main manuscript, participated in the feature extraction and created all the presented figures. HB participated in the feature extraction. JV helped with the interpretation of the results and wrote part of the discussion section. SP helped with the statistical analysis and wrote the related part in the methods section. ST and TS collected the data. ES helped with the writing of the manuscript and the correction of the multiple drafts. All authors read and approved the final manuscript.
Acknowledgements
The work was supported by the Pittsburgh Older Americans Independence Center (NIA P30 AG024827).
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.
Author information
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Gatouillat, A., Bleton, H., VanSwearingen, J. et al. Cognitive tasks during walking affect cerebral blood flow signal features in middle cerebral arteries and their correlation to gait characteristics. Behav Brain Funct 11, 29 (2015). https://doi.org/10.1186/s1299301500739
Received:
Accepted:
Published:
Keywords
 Transcranial Doppler
 Gait accelerometry
 Cerebral blood flow velocity
 Walking
 Signal features