Modeling BOLD-fMRI Hemodynamics via Multidimensional Decomposition of Electrophysiology Data: A Simulation Study

We propose a framework for studying the electrophysiological correlates of BOLD-fMRI. This framework relies on structured coupled matrix-tensor factorization (sCMTF), a joint multidimensional decomposition which reveals dynamical interactions between LFP/EEG oscillatory features and BOLD-fMRI data. We test whether LFP/EEG-BOLD co-fluctuations and regional hemodynamic response functions can be estimated by sCMTF using whole-brain modelling of restingstate activity. We produce permuted datasets to show that our framework extracts EEG/LFP temporal patterns that correlate significantly with BOLD signal fluctuations. Our framework is also capable of estimating HRFs that accurately embodies simulated hemodynamics, with a word of caution regarding initialization of the sCMTF algorithm.

Variability in the analysis of a single neuroimaging dataset by many teams

Data analysis workflows in many scientific domains have become increasingly complex and flexible. Here we assess the effect of this flexibility on the results of functional magnetic resonance imaging by asking 70 independent teams to analyse the same dataset, testing the same 9 ex-ante hypotheses1. The flexibility of analytical approaches is exemplified by the fact that no two teams chose identical workflows to analyse the data. This flexibility resulted in sizeable variation in the results of hypothesis tests, even for teams whose statistical maps were highly correlated at intermediate stages of the analysis pipeline. Variation in reported results was related to several aspects of analysis methodology. Notably, a meta-analytical approach that aggregated information across teams yielded a significant consensus in activated regions. Furthermore, prediction markets of researchers in the field revealed an overestimation of the likelihood of significant findings, even by researchers with direct knowledge of the dataset2-5. Our findings show that analytical flexibility can have substantial effects on scientific conclusions, and identify factors that may be related to variability in the analysis of functional magnetic resonance imaging. The results emphasize the importance of validating and sharing complex analysis workflows, and demonstrate the need for performing and reporting multiple analyses of the same data. Potential approaches that could be used to mitigate issues related to analytical variability are discussed.

Individualized growth prediction of mice skin tumors with maximum likelihood estimators

BACKGROUND & OBJECTIVE

In this work, we focus on estimating the parameters of the Gompertz model in order to predict tumor growth. The estimation is based on measurements from mice skin tumors of de novo carcinogenesis. The main objective is to compare the Maximum Likelihood estimator with the best performance from our previous work with the Non-linear Least Squares estimator which is commonly used in the literature to estimate the growth parameters of the Gompertz model.

METHODS

To describe tumor growth, we propose a stochastic model which is based on the Gompertz growth function. The principle of Maximum Likelihood is used to estimate both the growth rate and the carrying capacity of the Gompertz function, along with the characteristics of the additive Gaussian process and measurement noise. Moreover, we examine whether a Maximum A Posteriori estimator is able to utilize any available prior knowledge in order to improve the predictions.

RESULTS

Experimental data from a total of 24 tumors in 8 mice (3 tumors each) were used to study the performance of the proposed methods with respect to prediction accuracy. Our results show that the Maximum Likelihood estimator is able to provide, in most cases, more accurate predictions. Moreover, the Maximum A Posteriori estimator has the potential to correct potentially non-realistic estimates for the carrying capacity at early growth stages.

CONCLUSION

In most cases, the Maximum Likelihood estimator is able to provide more reliable predictions for the tumor's growth on individual test subjects. The Maximum A Posteriori estimator, it has the potential to improve the prediction when the available experimental data do not provide adequate information by utilizing prior knowledge about the unknown parameters.

Functional brain networks of patients with epilepsy exhibit pronounced multiscale periodicities, which correlate with seizure onset

Epileptic seizure detection and prediction by using noninvasive measurements such as scalp EEG signals or invasive, intracranial recordings, has been at the heart of epilepsy studies for at least three decades. To this end, the most common approach has been to consider short‐length recordings (several seconds to a few minutes) around a seizure, aiming to identify significant changes that occur before or during seizures. An inherent assumption in this approach is the presence of a relatively constant EEG activity in the interictal period, which is interrupted by seizure occurrence. Here, we examine this assumption by using long‐duration scalp EEG data (21–94 hr) in nine patients with epilepsy, based on which we construct functional brain networks. Our results reveal that these networks vary over time in a periodic fashion, exhibiting multiple peaks at periods ranging between 1 and 24 hr. The effects of seizure onset on the functional brain network properties were found to be considerably smaller in magnitude compared to the changes due to these inherent periodic cycles. Importantly, the properties of the identified network periodic components (instantaneous phase) were found to be strongly correlated to seizure onset, especially for the periodicities around 3 and 5 hr. These correlations were found to be largely absent between EEG signal periodicities and seizure onset, suggesting that higher specificity may be achieved by using network‐based metrics. In turn, this implies that more robust seizure detection and prediction can be achieved if longer term underlying functional brain network periodic variations are taken into account.

Assessment of dynamic cerebral autoregulation in humans: Is reproducibility dependent on blood pressure variability?

We tested the influence of blood pressure variability on the reproducibility of dynamic cerebral autoregulation (DCA) estimates. Data were analyzed from the 2^nd CARNet bootstrap initiative, where mean arterial blood pressure (MABP), cerebral blood flow velocity (CBFV) and end tidal CO2 were measured twice in 75 healthy subjects. DCA was analyzed by 14 different centers with a variety of different analysis methods. Intraclass Correlation (ICC) values increased significantly when subjects with low power spectral density MABP (PSD-MABP) values were removed from the analysis for all gain, phase and autoregulation index (ARI) parameters. Gain in the low frequency band (LF) had the highest ICC, followed by phase LF and gain in the very low frequency band. No significant differences were found between analysis methods for gain parameters, but for phase and ARI parameters, significant differences between the analysis methods were found. Alternatively, the Spearman-Brown prediction formula indicated that prolongation of the measurement duration up to 35 minutes may be needed to achieve good reproducibility for some DCA parameters. We conclude that poor DCA reproducibility (ICC<0.4) can improve to good (ICC > 0.6) values when cases with low PSD-MABP are removed, and probably also when measurement duration is increased.

Cerebrovascular Pulsatility During Rest and Exercise Reflects Hemodynamic Impairment in Stroke and Cerebral Small Vessel Disease

Although aerobic exercise is recommended as a core component of stroke rehabilitation, knowledge of acute cerebrovascular responses in patients is limited. This study tested the hypothesis that older adults with chronic stroke or cerebral small vessel disease (SVD) exhibit a greater increase in pulsatile hemodynamics during exercise compared with young and age-matched healthy adults. Middle cerebral artery blood flow velocity was acquired during 20 min of moderate intensity cycling in 51 participants from four groups (young, old, SVD and stroke). During rest, only the stroke group had a higher pulsatility index (PI) compared with the young group (1.02 ± 0.17 vs 0.83 ± 0.13; p = 0.038). During exercise, however, the SVD group exhibited a larger increase in PI (68 ± 20% relative to rest) than the young (47 ± 19%), old (45 ± 17%) and stroke (40 ± 25%) groups (p < 0.05, for each). The stress of aerobic exercise may reveal arterial dysfunction associated with latent and overt cerebrovascular disease.

Unfolding the Effects of Acute Cardiovascular Exercise on Neural Correlates of Motor Learning Using Convolutional Neural Networks

Cardiovascular exercise is known to promote the consolidation of newly acquired motor skills. Previous studies seeking to understand the neural correlates underlying motor memory consolidation that is modulated by exercise, have relied so far on using traditional statistical approaches for a priori selected features from neuroimaging data, including EEG. With recent advances in machine learning, data-driven techniques such as deep learning have shown great potential for EEG data decoding for brain-computer interfaces, but have not been explored in the context of exercise. Here, we present a novel Convolutional Neural Network (CNN)-based pipeline for analysis of EEG data to study the brain areas and spectral EEG measures modulated by exercise. To the best of our knowledge, this work is the first one to demonstrate the ability of CNNs to be trained in a limited sample size setting. Our approach revealed discriminative spectral features within a refined frequency band (27–29 Hz) as compared to the wider beta bandwidth (15–30 Hz), which is commonly used in data analyses, as well as corresponding brain regions that were modulated by exercise. These results indicate the presence of finer EEG spectral features that could have been overlooked using conventional hypothesis-driven statistical approaches. Our study thus demonstrates the feasibility of using deep network architectures for neuroimaging analysis, even in small-scale studies, to identify robust brain biomarkers and investigate neuroscience-based questions.

Effects of motion related outliers in dynamic functional connectivity using the sliding window method

BACKGROUND

It has been suggested that the use of window functions, other than the rectangular, in the sliding window method, may be beneficial for reducing the effects of motion-related outliers in the time-series, when assessing dynamic functional connectivity (dFC) in resting-state fMRI (rs-fMRI).

METHODOLOGY

Ten window functions for a wide range of window lengths (20-150 s) combined with Pearson and Kendall correlation metrics, were investigated. One hundred high quality rs-fMRI datasets from healthy controls, were used to systematically assess the effect of varying the window function and length on dFC assessment. To this end, two approaches were implemented: a) simulated outliers were added to the experimental data and b) the experimental data were divided into low and high motion subgroups.

RESULTS

The presence of experimental motion-noise tended to inflate the number of dynamic connections for longer (≥100 s) wide-shaped windows, while shorter (20-30 s) narrow-shaped windows exhibited increased sensitivity in the presence of simulated outliers. Moreover, window sizes from 60 s to 90 s were mildly affected by motion-related effects. In most cases, the number of dynamic connections increased, and gradually lower frequencies were captured, with an increasing window size.

CONCLUSIONS

Subject motion considerably affects the obtained dFC patterns; thus, it is preferable to perform motion artefact removal in the pre-processing stage rather than using alternative window functions to mitigate their effects. Provided that motion-noise is not excessive, the choice of a rectangular window is adequate. Finally, low frequency oscillations in functional connectivity seem to play an important role in the context of dFC assessment.

Identification of physiological response functions to correct for fluctuations in resting-state fMRI related to heart rate and respiration

Functional magnetic resonance imaging (fMRI) is widely viewed as the gold standard for studying brain function due to its high spatial resolution and non-invasive nature. However, it is well established that changes in breathing patterns and heart rate strongly influence the blood oxygen-level dependent (BOLD) fMRI signal and this, in turn, can have considerable effects on fMRI studies, particularly resting-state studies. The dynamic effects of physiological processes are often quantified by using convolution models along with simultaneously recorded physiological data. In this context, physiological response function (PRF) curves (cardiac and respiratory response functions), which are convolved with the corresponding physiological fluctuations, are commonly employed. While it has often been suggested that the PRF curves may be region- or subject-specific, it is still an open question whether this is the case. In the present study, we propose a novel framework for the robust estimation of PRF curves and use this framework to rigorously examine the implications of using population-, subject-, session- and scan-specific PRF curves. The proposed framework was tested on resting-state fMRI and physiological data from the Human Connectome Project. Our results suggest that PRF curves vary significantly across subjects and, to a lesser extent, across sessions from the same subject. These differences can be partly attributed to physiological variables such as the mean and variance of the heart rate during the scan. The proposed methodological framework can be used to obtain robust scan-specific PRF curves from data records with duration longer than 5 min, exhibiting significantly improved performance compared to previously defined canonical cardiac and respiration response functions. Besides removing physiological confounds from the BOLD signal, accurate modeling of subject- (or session-/scan-) specific PRF curves is of importance in studies that involve populations with altered vascular responses, such as aging subjects.

Older adults exhibit a more pronounced modulation of beta oscillations when performing sustained and dynamic handgrips

Muscle contractions are associated with a decrease in beta oscillatory activity, known as movement-related beta desynchronization (MRBD). Older adults exhibit a MRBD of greater amplitude compared to their younger counterparts, even though their beta power remains higher both at rest and during muscle contractions. Further, a modulation in MRBD has been observed during sustained and dynamic pinch contractions, whereby beta activity during periods of steady contraction following a dynamic contraction is elevated. However, how the modulation of MRBD is affected by aging has remained an open question. In the present work, we investigated the effect of aging on the modulation of beta oscillations and their putative link with motor performance. We collected magnetoencephalography (MEG) data from younger and older adults during a resting-state period and motor handgrip paradigms, which included sustained and dynamic contractions, to quantify spontaneous and motor-related beta oscillatory activity. Beta power at rest was found to be significantly increased in the motor cortex of older adults. During dynamic hand contractions, MRBD was more pronounced in older participants in frontal, premotor and motor brain regions. These brain areas also exhibited age-related decreases in cortical thickness; however, the magnitude of MRBD and cortical thickness were not found to be associated after controlling for age. During sustained hand contractions, MRBD exhibited a decrease in magnitude compared to dynamic contraction periods in both groups and did not show age-related differences. This suggests that the amplitude change in MRBD between dynamic and sustained contractions is larger in older compared to younger adults. We further probed for a relationship between beta oscillations and motor behaviour and found that greater MRBD in primary motor cortices was related to degraded motor performance beyond age, but our results suggested that age-related differences in beta oscillations were not predictive of motor performance.

Dynamic Cerebral Autoregulation Reproducibility Is Affected by Physiological Variability

Parameters describing dynamic cerebral autoregulation (DCA) have limited reproducibility. In an international, multi-center study, we evaluated the influence of multiple analytical methods on the reproducibility of DCA. Fourteen participating centers analyzed repeated measurements from 75 healthy subjects, consisting of 5 min of spontaneous fluctuations in blood pressure and cerebral blood flow velocity signals, based on their usual methods of analysis. DCA methods were grouped into three broad categories, depending on output types: (1) transfer function analysis (TFA); (2) autoregulation index (ARI); and (3) correlation coefficient. Only TFA gain in the low frequency (LF) band showed good reproducibility in approximately half of the estimates of gain, defined as an intraclass correlation coefficient (ICC) of >0.6. None of the other DCA metrics had good reproducibility. For TFA-like and ARI-like methods, ICCs were lower than values obtained with surrogate data (p < 0.05). For TFA-like methods, ICCs were lower for the very LF band (gain 0.38 ± 0.057, phase 0.17 ± 0.13) than for LF band (gain 0.59 ± 0.078, phase 0.39 ± 0.11, p ≤ 0.001 for both gain and phase). For ARI-like methods, the mean ICC was 0.30 ± 0.12 and for the correlation methods 0.24 ± 0.23. Based on comparisons with ICC estimates obtained from surrogate data, we conclude that physiological variability or non-stationarity is likely to be the main reason for the poor reproducibility of DCA parameters.

Modeling of the BOLD signal using event-related simultaneous EEG-fMRI and convolutional sparse coding analysis

In this work, we employ simultaneous EEG-fMRI data acquired during a visually-guided attention task along with convolutional sparse coding (CSC) analysis to extract transient events from the EEG. Subsequently, we use these events in a standard voxel-wise fMRI analysis and compare the resultant activation maps with maps obtained using the subjects' response time (RT) in detection of visual target stimuli. We also employ FIR models to obtain HRF estimates using the detected CSC events. Our results show concordance between the resultant activation maps and consistent HRF shapes for most of the subjects, suggesting that CSC can be used as a tool for the detection of reliable events in the EEG.

Assessment of dynamic functional connectivity in resting-state fMRI using the sliding window technique

INTRODUCTION

Recent studies related to assessing functional connectivity (FC) in resting-state functional magnetic resonance imaging have revealed that the resulting connectivity patterns exhibit considerable fluctuations (dynamic FC [dFC]). A widely applied method for quantifying dFC is the sliding window technique. According to this method, the data are divided into segments with the same length (window size) and a correlation metric is employed to assess the connectivity within these segments, whereby the window size is often empirically chosen.

METHODS

In this study, we rigorously investigate the assessment of dFC using the sliding window approach. Specifically, we perform a detailed comparison between different correlation metrics, including Pearson, Spearman and Kendall correlation, Pearson and Spearman partial correlation, Mutual Information (MI), Variation of Information (VI), Kullback-Leibler divergence, Multiplication of Temporal Derivatives and Inverse Covariance.

RESULTS

Using test-retest datasets, we show that MI and VI yielded the most consistent results by achieving high reliability with respect to dFC estimates for different window sizes. Subsequent hypothesis testing, based on multivariate phase randomization surrogate data generation, allowed the identification of dynamic connections between the posterior cingulate cortex and regions in the frontal lobe and inferior parietal lobes, which were overall in agreement with previous studies.

CONCLUSIONS

In the case of MI and VI, a window size of at least 120 s was found to be necessary for detecting dFC for some of the previously identified dynamically connected regions.

Graph Theoretical Characteristics of EEG-Based Functional Brain Networks in Patients With Epilepsy: The Effect of Reference Choice and Volume Conduction

It is well-established that both volume conduction and the choice of recording reference (montage) affect the correlation measures obtained from scalp EEG, both in the time and frequency domains. As a result, a number of correlation measures have been proposed aiming to reduce these effects. In our previous work, we have showed that scalp-EEG based functional brain networks in patients with epilepsy exhibit clear periodic patterns at different time scales and that these patterns are strongly correlated to seizure onset, particularly at shorter time scales (around 3 and 5 h), which has important clinical implications. In the present work, we use the same long-duration clinical scalp EEG data (multiple days) to investigate the extent to which the aforementioned results are affected by the choice of reference choice and correlation measure, by considering several widely used montages as well as correlation metrics that are differentially sensitive to the effects of volume conduction. Specifically, we compare two standard and commonly used linear correlation measures, cross-correlation in the time domain, and coherence in the frequency domain, with measures that account for zero-lag correlations: corrected cross-correlation, imaginary coherence, phase lag index, and weighted phase lag index. We show that the graphs constructed with corrected cross-correlation and WPLI are more stable across different choices of reference. Also, we demonstrate that all the examined correlation measures revealed similar periodic patterns in the obtained graph measures when the bipolar and common reference (Cz) montage were used. This includes circadian-related periodicities (e.g., a clear increase in connectivity during sleep periods as compared to awake periods), as well as periodicities at shorter time scales (around 3 and 5 h). On the other hand, these results were affected to a large degree when the average reference montage was used in combination with standard cross-correlation, coherence, imaginary coherence, and PLI, which is likely due to the low number of electrodes and inadequate electrode coverage of the scalp. Finally, we demonstrate that the correlation between seizure onset and the brain network periodicities is preserved when corrected cross-correlation and WPLI were used for all the examined montages. This suggests that, even in the standard clinical setting of EEG recording in epilepsy where only a limited number of scalp EEG measurements are available, graph-theoretic quantification of periodic patterns using appropriate montage, and correlation measures corrected for volume conduction provides useful insights into seizure onset.

Functional and effective reorganization of the aging brain during unimanual and bimanual hand movements

Motor performance decline observed during aging is linked to changes in brain structure and function, however, the precise neural reorganization associated with these changes remains largely unknown. We investigated the neurophysiological correlates of this reorganization by quantifying functional and effective brain network connectivity in elderly individuals (n = 11; mean age = 67.5 years), compared to young adults (n = 12; mean age = 23.7 years), while they performed visually-guided unimanual and bimanual handgrips inside the magnetoencephalography (MEG) scanner. Through a combination of principal component analysis and Granger causality, we observed age-related increases in functional and effective connectivity in whole-brain, task-related motor networks. Specifically, elderly individuals demonstrated (i) greater information flow from contralateral parietal and ipsilateral secondary motor regions to the left primary motor cortex during the unimanual task and (ii) decreased interhemispheric temporo-frontal communication during the bimanual task. Maintenance of motor performance and task accuracy in elderly was achieved by hyperactivation of the task-specific motor networks, reflecting a possible mechanism by which the aging brain recruits additional resources to counteract known myelo- and cytoarchitectural changes. Furthermore, resting-state sessions acquired before and after each motor task revealed that both older and younger adults maintain the capacity to adapt to task demands via network-wide increases in functional connectivity. Collectively, our study consolidates functional connectivity and directionality of information flow in systems-level cortical networks during aging and furthers our understanding of neuronal flexibility in motor processes.

A Novel Framework for Estimating Time-Varying Multivariate Autoregressive Models and Application to Cardiovascular Responses to Acute Exercise

OBJECTIVE

We present a novel modeling framework for identifying time-varying (TV) couplings between time-series of biomedical relevance.

METHODS

The proposed methodology is based on multivariate autoregressive (MVAR) models, which have been extensively used to study couplings between biosignals. Contrary to the standard estimation methods that assume time-invariant relationships, we propose a modified recursive Kalman filter (KF) to track changes in the model parameters. We perform model order selection and hyperparameter optimization simultaneously using Genetic Algorithms, greatly improving accuracy and computation time. In addition, we address the effect of residual heteroscedasticity, possibly associated with event-related changes or phase transitions during a given experimental protocol, on the TV-MVAR coupling measures by using Generalized Autoregressive Conditional Heteroskedasticity (GARCH) models to fit the TV-MVAR residuals.

RESULTS

Using simulated data, we show that the proposed framework yields more accurate parameter estimates compared to the conventional KF, particularly when the true system parameters exhibit different rate of variations over time. Furthermore, by accounting for heteroskedasticity, we obtain more accurate estimates of the strength and directionality of the underlying couplings. We also use our approach to investigate TV hemodynamic interactions during exercise in young and old healthy adults, as well as individuals with chronic stroke. We extract TV coupling patterns that reflect well known exercise-induced effects on the underlying regulatory mechanisms with excellent time resolution.

CONCLUSION AND SIGNIFICANCE

The proposed modeling framework can be used to efficiently quantify TV couplings between biosignals. It is fully automated and does not require prior knowledge of the system TV characteristics.

Reproducibility of dynamic cerebral autoregulation parameters: a multi-centre, multi-method study

OBJECTIVE

Different methods to calculate dynamic cerebral autoregulation (dCA) parameters are available. However, most of these methods demonstrate poor reproducibility that limit their reliability for clinical use. Inter-centre differences in study protocols, modelling approaches and default parameter settings have all led to a lack of standardisation and comparability between studies. We evaluated reproducibility of dCA parameters by assessing systematic errors in surrogate data resulting from different modelling techniques.

APPROACH

Fourteen centres analysed 22 datasets consisting of two repeated physiological blood pressure measurements with surrogate cerebral blood flow velocity signals, generated using Tiecks curves (autoregulation index, ARI 0-9) and added noise. For reproducibility, dCA methods were grouped in three broad categories: 1. Transfer function analysis (TFA)-like output; 2. ARI-like output; 3. Correlation coefficient-like output. For all methods, reproducibility was determined by one-way intraclass correlation coefficient analysis (ICC).

MAIN RESULTS

For TFA-like methods the mean (SD; [range]) ICC gain was 0.71 (0.10; [0.49-0.86]) and 0.80 (0.17; [0.36-0.94]) for VLF and LF (p  =  0.003) respectively. For phase, ICC values were 0.53 (0.21; [0.09-0.80]) for VLF, and 0.92 (0.13; [0.44-1.00]) for LF (p  <  0.001). Finally, ICC for ARI-like methods was equal to 0.84 (0.19; [0.41-0.94]), and for correlation-like methods, ICC was 0.21 (0.21; [0.056-0.35]).

SIGNIFICANCE

When applied to realistic surrogate data, free from the additional exogenous influences of physiological variability on cerebral blood flow, most methods of dCA modelling showed ICC values considerably higher than what has been reported for physiological data. This finding suggests that the poor reproducibility reported by previous studies may be mainly due to the inherent physiological variability of cerebral blood flow regulatory mechanisms rather than related to (stationary) random noise and the signal analysis methods.

Modeling of dynamic cerebrovascular reactivity to spontaneous and externally induced CO2 fluctuations in the human brain using BOLD-fMRI

In this work, we investigate the regional characteristics of the dynamic interactions between arterial CO₂ and BOLD (dynamic cerebrovascular reactivity - dCVR) during normal breathing and hypercapnic, externally induced step CO₂ challenges. To obtain dCVR curves at each voxel, we use a custom set of basis functions based on the Laguerre and gamma basis sets. This allows us to obtain robust dCVR estimates both in larger regions of interest (ROIs), as well as in individual voxels. We also implement classification schemes to identify brain regions with similar dCVR characteristics. Our results reveal considerable variability of dCVR across different brain regions, as well as during different experimental conditions (normal breathing and hypercapnic challenges), suggesting a differential response of cerebral vasculature to spontaneous CO₂ fluctuations and larger, externally induced CO₂ changes that are possibly associated with the underlying differences in mean arterial CO₂ levels. The clustering results suggest that anatomically distinct brain regions are characterized by different dCVR curves that in some cases do not exhibit the standard, positive valued curves that have been previously reported. They also reveal a consistent set of dCVR cluster shapes for resting and forcing conditions, which exhibit different distribution patterns across brain voxels.

Prediction of the Time to Syncope Occurrence in Patients Diagnosed with Vasovagal Syncope

OBJECTIVE

In this study we aimed to predict the time to syncope occurrence (TSO) in patients with vasovagal syncope (VVS), solely based on measurements recorded during the supine position of the head-up tilt (HUT) testing protocol.

METHODS

We extracted various time and frequency domain features related to morphological aspects of arterial blood pressure (ABP) and the electrocardiogram (ECG) raw signals as well as to dynamic interactions between beat-to-beat ABP, heart rate, and cerebral blood flow velocity. From these we identified the most predictive features related to TSO.

RESULTS

Specifically, when no orthostatic stress is involved, TSO in VVS patients can be predicted with high accuracy from a set of only five ECG features.

Tumor growth modeling: Parameter estimation with Maximum Likelihood methods

BACKGROUND & OBJECTIVE

In this work, we focus on estimating the parameters of the widely used Gompertz tumor growth model, based on measurements of the tumor's volume. Being able to accurately describe the dynamics of tumor growth on an individual basis is very important both for growth prediction and designing personalized, optimal therapy schemes (e.g. when using model predictive control).

METHODS

Our analysis aims to compute both the growth rate and the carrying capacity of the Gompertz function, along with the characteristics of the additive Gaussian process and measurement noise of the system. Three methods based on Maximum Likelihood estimation are proposed. The first utilizes an assumption regarding the measurement noise that simplifies the problem, the second combines the Extended Kalman Filter and Maximum Likelihood estimation, and the third is a nonstandard exact form of Maximum Likelihood estimation, where numerical integration is used to approximate the likelihood of the measurements, along with a novel way to reduce the required computations.

RESULTS

Synthetic data were used in order to perform extensive simulations aiming to compare the methods' effectiveness, with respect to the accuracy of the estimation. The proposed methods are able to estimate the growth dynamics, even when the noise characteristics are not estimated accurately. Another very important finding is that the methods perform best in the case that corresponds to the problem needed to be solved when dealing with experimental data.

CONCLUSION

Using nonstandard, problem specific techniques can improve the estimation accuracy and best exploit the available data.

Identification of Time-Varying Cortico-cortical and Cortico-Muscular Coherence during Motor Tasks with Multivariate Autoregressive Models

Neural populations coordinate at fast subsecond time-scales during rest and task execution. As a result, functional brain connectivity assessed with different neuroimaging modalities (EEG, MEG, fMRI) may also change over different time scales. In addition to the more commonly used sliding window techniques, the General Linear Kalman Filter (GLFK) approach has been proposed to estimate time-varying brain connectivity. In the present work, we propose a modification of the GLFK approach to model timevarying connectivity. We also propose a systematic method to select the hyper-parameters of the model. We evaluate the performance of the method using MEG and EMG data collected from 12 young subjects performing two motor tasks (unimanual and bimanual hand grips), by quantifying time-varying cortico-cortical and corticomuscular coherence (CCC and CMC). The CMC results revealed patterns in accordance with earlier findings, as well as an improvement in both time and frequency resolution compared to sliding window approaches. These results suggest that the proposed methodology is able to unveil accurate time-varying connectivity patterns with an excellent time resolution.

Investigation Of Interaction Between Physiological Signals And fMRI Dynamic Functional Connectivity Using Independent Component Analysis

The blood oxygen level dependent (BOLD) fMRI signal is influenced not only by neuronal activity but also by fluctuations in physiological signals, including respiration, arterial CO2 and heart rate/ heart rate variability (HR/HRV). Even spontaneous physiological signal fluctuations have been shown to influence the BOLD fMRI signal in a regionally specific manner. Consequently, estimates of functional connectivity between different brain regions, performed when the subject is at rest, may be confounded by the effects of physiological signal fluctuations. In addition, resting functional connectivity has been shown to vary with respect to time (dynamic functional connectivity - DFC), with the sources of this variation not fully elucidated. The effect of physiological factors on dynamic (time-varying) resting-state functional connectivity has not been studied extensively, to our knowledge. In our previous study, we investigated the effect of heart rate (HR) and end-tidal CO2 (PETCO2) on the time-varying network degree of three well-described RSNs (DMN, SMN and Visual Network) using mask-based and seed-based analysis, and we identified brain-heart interactions which were more pronounced in specific frequency bands. Here, we extend this work, by estimating DFC and its corresponding network degree for the RSNs, employing a data-driven approach to extract the RSNs (low-and high-dimensional Independent Component Analysis (ICA)), which we subsequently correlate with the characteristics of simultaneously collected physiological signals. The results confirm that physiological signals have a modulatory effect on resting-state, fMRI-based DFC.

Effects of continuous positive airway pressure and isocapnic-hypoxia on cerebral autoregulation in patients with obstructive sleep apnoea

KEY POINTS

Altered cerebral autoregulation (CA) in obstructive sleep apnoea (OSA) patients may contribute to increased stroke risk in this population; the gold standard treatment for OSA is continuous positive airway pressure, which improves cerebrovascular regulation and may decrease the risk of stroke. Isocapnic-hypoxia impairs CA in healthy subjects, but it remains unknown in OSA whether impaired CA is further exacerbated by isocapnic-hypoxia and whether it is improved by treatment with continuous positive airway pressure. During normoxia, CA was altered in the more severe but not in the less severe OSA patients, while, in contrast, during isocapnic-hypoxia, CA was similar between groups and tended to improve in patients with more severe OSA compared to normoxia. From a clinical perspective, one month of continuous positive airway pressure treatment does not improve CA. From a physiological perspective, this study suggests that sympathetic overactivity may be responsible for altered CA in the more severe OSA patients.

ABSTRACT

Cerebral autoregulation (CA) impairment may contribute to the increased risk of stroke associated with obstructive sleep apnoea (OSA). It is unknown if impaired CA is further exacerbated by isocapnic-hypoxia and whether it is improved by treatment of OSA with continuous positive airway pressure (CPAP). CA was assessed during wakefulness in 53 OSA patients (50.3 ± 9.3 years) and 21 controls (49.8 ± 8.6 years) at baseline and following a minimum of 1 month of effective CPAP therapy (OSA patients, n = 40). Control participants (n = 21) performed a follow-up visit to control for time effects within OSA patients between baseline and the post-CPAP visit. Beat-by-beat middle cerebral artery blood flow velocity and mean arterial blood pressure (MBP), and breath-by-breath end-tidal partial pressure of CO₂ (P ET ,CO2) were monitored. CA was determined during normoxia and isocapnic-hypoxia using transfer function (phase and gain) and coherence analysis (including multiple and partial coherence (using MBP and P ET ,CO2 as inputs)) in the very low frequency range (0.03-0.07 Hz). OSA patients were divided into two subgroups (less severe and more severe) based upon the median respiratory disturbance index (RDI). During normoxia, the more severe OSA patients (RDI 45.9 ± 10.3) exhibited altered CA compared to controls and the less severe OSA patients (RDI 24.5 ± 5.9). In contrast, during isocapnic-hypoxia, CA was similar between groups. CPAP had no effect on CA. In conclusion, CA is altered in the more severe OSA patients during normoxia but not during isocapnic-hypoxia and CPAP treatment does not impact CA.

Estimation of voxel-wise dynamic cerebrovascular reactivity curves from resting-state fMRI data

In this work, we investigate the linear dynamic interactions between fluctuations in arterial CO2 that occur during normal breathing, and the BOLD fMRI signal. We cast this problem within a systems-theoretic framework, where we employ functional expansions for the estimation of the impulse responses in large regions of interest, as well as in individual voxels. We also implement classification schemes in order to identify different brain regions with similar cerebrovascular reactivity characteristics. Our results reveal that it is feasible to obtain reliable estimates of cerebrovascular reactivity curves from resting-state data and that these curves exhibit considerable variability across different brain regions that may be related to the underlying anatomy.

Unsupervised detection and removal of muscle artifacts from scalp EEG recordings using canonical correlation analysis, wavelets and random forests

OBJECTIVE

This paper proposes supervised and unsupervised algorithms for automatic muscle artifact detection and removal from long-term EEG recordings, which combine canonical correlation analysis (CCA) and wavelets with random forests (RF).

METHODS

The proposed algorithms first perform CCA and continuous wavelet transform of the canonical components to generate a number of features which include component autocorrelation values and wavelet coefficient magnitude values. A subset of the most important features is subsequently selected using RF and labelled observations (supervised case) or synthetic data constructed from the original observations (unsupervised case). The proposed algorithms are evaluated using realistic simulation data as well as 30min epochs of non-invasive EEG recordings obtained from ten patients with epilepsy.

RESULTS

We assessed the performance of the proposed algorithms using classification performance and goodness-of-fit values for noisy and noise-free signal windows. In the simulation study, where the ground truth was known, the proposed algorithms yielded almost perfect performance. In the case of experimental data, where expert marking was performed, the results suggest that both the supervised and unsupervised algorithm versions were able to remove artifacts without affecting noise-free channels considerably, outperforming standard CCA, independent component analysis (ICA) and Lagged Auto-Mutual Information Clustering (LAMIC).

CONCLUSION

The proposed algorithms achieved excellent performance for both simulation and experimental data. Importantly, for the first time to our knowledge, we were able to perform entirely unsupervised artifact removal, i.e. without using already marked noisy data segments, achieving performance that is comparable to the supervised case.

SIGNIFICANCE

Overall, the results suggest that the proposed algorithms yield significant future potential for improving EEG signal quality in research or clinical settings without the need for marking by expert neurophysiologists, EMG signal recording and user visual inspection.