Spiking Laguerre Volterra networks - predicting neuronal activity from local field potentials.

OBJECTIVE: Understanding the generative mechanism between Local Field Potentials (LFP) and neuronal spiking activity is a crucial step for understanding information processing in the brain. Up to now, most approaches have relied on simply quantifying the coupling between LFP and spikes. However, very few have managed to predict the exact timing of spike occurrence based on LFP variations. APPROACH: Here, we fill this gap by proposing novel spiking Laguerre-Volterra Network (sLVN) models to describe the dynamic LFP-spike relationship. Compared to conventional artificial neural networks, the sLVNs are interpretable models that provide explainable features of the underlying dynamics. MAIN RESULTS: The proposed networks were applied on extracellular microelectrode recordings of Parkinson's Disease (PD) patients during Deep Brain Stimulation (DBS) surgery. Based on the predictability of the LFP-spike pairs, we detected three neuronal populations with unique signal characteristics and sLVN model features. SIGNIFICANCE: These clusters were indirectly associated with motor score improvement following DBS surgery, warranting further investigation into the potential of spiking activity predictability as an intraoperative biomarker for optimal DBS lead placement.

Cancer cell sedimentation in 3D cultures reveals active migration regulated by self-generated gradients and adhesion sites.

Cell sedimentation in 3D hydrogel cultures refers to the vertical migration of cells towards the bottom of the space. Understanding this poorly examined phenomenon may allow us to design better protocols to prevent it, as well as provide insights into the mechanobiology of cancer development. We conducted a multiscale experimental and mathematical examination of 3D cancer growth in triple negative breast cancer cells. Migration was examined in the presence and absence of Paclitaxel, in high and low adhesion environments and in the presence of fibroblasts. The observed behaviour was modeled by hypothesizing active migration due to self-generated chemotactic gradients. Our results did not reject this hypothesis, whereby migration was likely to be regulated by the MAPK and TGF-ß pathways. The mathematical model enabled us to describe the experimental data in absence (normalized error<40%) and presence of Paclitaxel (normalized error<10%), suggesting inhibition of random motion and advection in the latter case. Inhibition of sedimentation in low adhesion and co-culture experiments further supported the conclusion that cells actively migrated downwards due to the presence of signals produced by cells already attached to the adhesive glass surface.

On the variability of dynamic functional connectivity assessment methods.

BACKGROUND: Dynamic functional connectivity (dFC) has become an important measure for understanding brain function and as a potential biomarker. However, various methodologies have been developed for assessing dFC, and it is unclear how the choice of method affects the results. In this work, we aimed to study the results variability of commonly used dFC methods. METHODS: We implemented 7 dFC assessment methods in Python and used them to analyze the functional magnetic resonance imaging data of 395 subjects from the Human Connectome Project. We measured the similarity of dFC results yielded by different methods using several metrics to quantify overall, temporal, spatial, and intersubject similarity. RESULTS: Our results showed a range of weak to strong similarity between the results of different methods, indicating considerable overall variability. Somewhat surprisingly, the observed variability in dFC estimates was found to be comparable to the expected functional connectivity variation over time, emphasizing the impact of methodological choices on the final results. Our findings revealed 3 distinct groups of methods with significant intergroup variability, each exhibiting distinct assumptions and advantages. CONCLUSIONS: Overall, our findings shed light on the impact of dFC assessment analytical flexibility and highlight the need for multianalysis approaches and careful method selection to capture the full range of dFC variation. They also emphasize the importance of distinguishing neural-driven dFC variations from physiological confounds and developing validation frameworks under a known ground truth. To facilitate such investigations, we provide an open-source Python toolbox, PydFC, which facilitates multianalysis dFC assessment, with the goal of enhancing the reliability and interpretability of dFC studies.

Time-domain methods for quantifying dynamic cerebral blood flow autoregulation: Review and recommendations. A white paper from the Cerebrovascular Research Network (CARNet).

Cerebral Autoregulation (CA) is an important physiological mechanism stabilizing cerebral blood flow (CBF) in response to changes in cerebral perfusion pressure (CPP). By maintaining an adequate, relatively constant supply of blood flow, CA plays a critical role in brain function. Quantifying CA under different physiological and pathological states is crucial for understanding its implications. This knowledge may serve as a foundation for informed clinical decision-making, particularly in cases where CA may become impaired. The quantification of CA functionality typically involves constructing models that capture the relationship between CPP (or arterial blood pressure) and experimental measures of CBF. Besides describing normal CA function, these models provide a means to detect possible deviations from the latter. In this context, a recent white paper from the Cerebrovascular Research Network focused on Transfer Function Analysis (TFA), which obtains frequency domain estimates of dynamic CA. In the present paper, we consider the use of time-domain techniques as an alternative approach. Due to their increased flexibility, time-domain methods enable the mitigation of measurement/physiological noise and the incorporation of nonlinearities and time variations in CA dynamics. Here, we provide practical recommendations and guidelines to support researchers and clinicians in effectively utilizing these techniques to study CA.

Non-invasive estimation of arterial blood pressure fluctuations using a peripheral photoplethysmograph inside the MRI scanner.

The blood-oxygen-level-dependent (BOLD) signal measured by functional magnetic resonance imaging (fMRI) is modulated by neural activity through the neurovascular coupling effect, as well as non-neural factors of physiological origin such as heart rate, respiration, and arterial blood pressure (ABP). While the former two effects have been previously characterized, the modulation of the BOLD signal by ABP fluctuations is still poorly understood. This is largely due to the difficulty of obtaining reliable ABP measurements in the MRI environment. Here, we propose a combined experimental and mathematical modeling framework to estimate ABP fluctuations inside the MRI scanner using photoplethysmography (PPG). Specifically, we used concurrent PPG and ABP measurements obtained outside the scanner to train the mathematical model and applied it to PPG measurements obtained inside the MRI scanner. Our results suggest good agreement between the model-predicted and experimentally measured ABP fluctuations and region specific correlations with the BOLD fluctuations.

Investigation of Methodologies for Extracting Individual Brain Oscillations: Comparisons & Insights.

There is increasing evidence that the effects of non-invasive brain stimulation can be maximized when the applied intervention matches internal brain oscillations. Extracting individual brain oscillations is thus a necessary step for implementing personalized brain stimulation. In this context, different methods have been proposed for obtaining subject-specific spectral peaks from electrophysiological recordings. However, comparing the results obtained using different approaches is still lacking. Therefore, in the present work, we examined the following methodologies in terms of obtaining individual motor-related EEG spectral peaks: fast Fourier Transform analysis, power spectrum density analysis, wavelet analysis, and a principal component based time-frequency analysis. We used EEG data obtained when performing two different motor tasks - a hand grip task and a hand opening- and-closing task. Our results showed that both the motor task type and the specific method for performing the analysis had considerable impact on the extraction of subject-specific oscillation spectral peaks.Clinical Relevance-This exploratory study provides insights into the potential effects of using different methods to extract individual brain oscillations, which is important for designing personalized brain-machine-interfaces.

The default network dominates neural responses to evolving movie stories.

Neuroscientific studies exploring real-world dynamic perception often overlook the influence of continuous changes in narrative content. In our research, we utilize machine learning tools for natural language processing to examine the relationship between movie narratives and neural responses. By analyzing over 50,000 brain images of participants watching Forrest Gump from the studyforrest dataset, we find distinct brain states that capture unique semantic aspects of the unfolding story. The default network, associated with semantic information integration, is the most engaged during movie watching. Furthermore, we identify two mechanisms that underlie how the default network liaises with the amygdala and hippocampus. Our findings demonstrate effective approaches to understanding neural processes in everyday situations and their relation to conscious awareness.

Bioprinted Multicomponent Hydrogel Co-culture Tumor-Immune Model for Assessing and Simulating Tumor-Infiltrated Lymphocyte Migration and Functional Activation.

The immune response against a tumor is characterized by the interplay among components of the immune system and neoplastic cells. Here, we bioprinted a model with two distinct regions containing gastric cancer patient-derived organoids (PDOs) and tumor-infiltrated lymphocytes (TILs). The initial cellular distribution allows for the longitudinal study of TIL migratory patterns concurrently with multiplexed cytokine analysis. The chemical properties of the bioink were designed to present physical barriers that immune T-cells must breech during infiltration and migration toward a tumor with the use of an alginate, gelatin, and basal membrane mix. TIL activity, degranulation, and regulation of proteolytic activity reveal insights into the time-dependent biochemical dynamics. Regulation of the sFas and sFas-ligand present on PDOs and TILs, respectively, and the perforin and granzyme longitudinal secretion confirms TIL activation when encountering PDO formations. TIL migratory profiles were used to create a deterministic reaction-advection diffusion model. The simulation provides insights that decouple passive from active cell migration mechanisms. The mechanisms used by TILs and other adoptive cell therapeutics as they infiltrate the tumor barrier are poorly understood. This study presents a pre-screening strategy for immune cells where motility and activation across ECM environments are crucial indicators of cellular fitness.

Deep learning prediction of motor performance in stroke individuals using neuroimaging data.

The degree of motor impairment and profile of recovery after stroke are difficult to predict for each individual. Measures obtained from clinical assessments, as well as neurophysiological and neuroimaging techniques have been used as potential biomarkers of motor recovery, with limited accuracy up to date. To address this, the present study aimed to develop a deep learning model based on structural brain images obtained from stroke participants and healthy volunteers. The following inputs were used in a multi-channel 3D convolutional neural network (CNN) model: fractional anisotropy, mean diffusivity, radial diffusivity, and axial diffusivity maps obtained from Diffusion Tensor Imaging (DTI) images, white and gray matter intensity values obtained from Magnetic Resonance Imaging, as well as demographic data (e.g., age, gender). Upper limb motor function was classified into "Poor" and "Good" categories. To assess the performance of the DL model, we compared it to more standard machine learning (ML) classifiers including k-nearest neighbor, support vector machines (SVM), Decision Trees, Random Forests, Ada Boosting, and Naïve Bayes, whereby the inputs of these classifiers were the features taken from the fully connected layer of the CNN model. The highest accuracy and area under the curve values were 0.92 and 0.92 for the 3D-CNN and 0.91 and 0.91 for the SVM, respectively. The multi-channel 3D-CNN with residual blocks and SVM supported by DL was more accurate than traditional ML methods to classify upper limb motor impairment in the stroke population. These results suggest that combining volumetric DTI maps and measures of white and gray matter integrity can improve the prediction of the degree of motor impairment after stroke. Identifying the potential of recovery early on after a stroke could promote the allocation of resources to optimize the functional independence of these individuals and their quality of life.

Modeling the dynamics of cerebrovascular reactivity to carbon dioxide in fMRI under task and resting-state conditions.

Conventionally, cerebrovascular reactivity (CVR) is estimated as the amplitude of the hemodynamic response to vascular stimuli, most commonly carbon dioxide (CO2). While the CVR amplitude has established clinical utility, the temporal characteristics of CVR (dCVR) have been increasingly explored and may yield even more pathology-sensitive parameters. This work is motivated by the current need to evaluate the feasibility of dCVR modeling in various experimental conditions. In this work, we present a comparison of several recently published/utilized model-based deconvolution (response estimation) approaches for estimating the CO2 response function h(t), including maximum a posteriori likelihood (MAP), inverse logit (IL), canonical correlation analysis (CCA), and basis expansion (using Gamma and Laguerre basis sets). To aid the comparison, we devised a novel simulation framework that incorporates a wide range of SNRs, ranging from 10 to -7 dB, representative of both task and resting-state CO2 changes. In addition, we built ground-truth h(t) into our simulation framework, overcoming the conventional limitation that the true h(t) is unknown. Moreover, to best represent realistic noise found in fMRI scans, we extracted noise from in-vivo resting-state scans. Furthermore, we introduce a simple optimization of the CCA method (CCAopt) and compare its performance to these existing methods. Our findings suggest that model-based methods can accurately estimate dCVR even amidst high noise (i.e. resting-state), and in a manner that is largely independent of the underlying model assumptions for each method. We also provide a quantitative basis for making methodological choices, based on the desired dCVR parameters, the estimation accuracy and computation time. The BEL method provided the highest accuracy and robustness, followed by the CCAopt and IL methods. Of the three, the CCAopt method has the lowest computational requirements. These findings lay the foundation for wider adoption of dCVR estimation in CVR mapping.

Quantifying the Morphology and Mechanisms of Cancer Progression in 3D In-Vitro Environments: Integrating Experiments and Multiscale Models.

Mathematical models of cancer growth have become increasingly more accurate both in the space and time domains. However, the limited amount of data typically available has resulted in a larger number of qualitative rather than quantitative studies. In the present study, we provide an integrated experimental-computational framework for the quantification of the morphological characteristics and the mechanistic modelling of cancer progression in 3D environments. The proposed framework allows for the calibration of multiscale, spatiotemporal models of cancer growth using state-of-the-art 3D cell culture data, and their validation based on the resulting experimental morphological patterns using spatial point-pattern analysis techniques. We applied this framework to the study of the development of Triple Negative Breast Cancer cells cultured in Matrigel scaffolds, and validated the hypothesis of chemotactic migration using a multiscale, hybrid Keller-Segel model. The results revealed transient, non-random spatial distributions of cancer cells that consist of clustered, and dispersion patterns. The proposed model was able to describe the general characteristics of the experimental observations and suggests that chemotactic migration together with random motion was found to be a plausible mechanism leading to accumulation, during the examined time period of development. The developed framework enabled us to pursue two goals; first, the quantitative description of the morphology of cancer growth in 3D cultures using point-pattern analysis, and second, the relation of tumour morphology with underlying biophysical mechanisms that govern cancer growth and migration.

Detection and Spatiotemporal Analysis of In-vitro 3D Migratory Triple-Negative Breast Cancer Cells.

The invasion of cancer cells into the surrounding tissues is one of the hallmarks of cancer. However, a precise quantitative understanding of the spatiotemporal patterns of cancer cell migration and invasion still remains elusive. A promising approach to investigate these patterns are 3D cell cultures, which provide more realistic models of cancer growth compared to conventional 2D monolayers. Quantifying the spatial distribution of cells in these 3D cultures yields great promise for understanding the spatiotemporal progression of cancer. In the present study, we present an image processing and segmentation pipeline for the detection of 3D GFP-fluorescent triple-negative breast cancer cell nuclei, and we perform quantitative analysis of the formed spatial patterns and their temporal evolution. The performance of the proposed pipeline was evaluated using experimental 3D cell culture data, and was found to be comparable to manual segmentation, outperforming four alternative automated methods. The spatiotemporal statistical analysis of the detected distributions of nuclei revealed transient, non-random spatial distributions that consisted of clustered patterns across a wide range of neighbourhood distances, as well as dispersion for larger distances. Overall, the implementation of the proposed framework revealed the spatial organization of cellular nuclei with improved accuracy, providing insights into the 3 dimensional inter-cellular organization and its progression through time.

Transfer function analysis of dynamic cerebral autoregulation: A CARNet white paper 2022 update.

Cerebral autoregulation (CA) refers to the control of cerebral tissue blood flow (CBF) in response to changes in perfusion pressure. Due to the challenges of measuring intracranial pressure, CA is often described as the relationship between mean arterial pressure (MAP) and CBF. Dynamic CA (dCA) can be assessed using multiple techniques, with transfer function analysis (TFA) being the most common. A 2016 white paper by members of an international Cerebrovascular Research Network (CARNet) that is focused on CA strove to improve TFA standardization by way of introducing data acquisition, analysis, and reporting guidelines. Since then, additional evidence has allowed for the improvement and refinement of the original recommendations, as well as for the inclusion of new guidelines to reflect recent advances in the field. This second edition of the white paper contains more robust, evidence-based recommendations, which have been expanded to address current streams of inquiry, including optimizing MAP variability, acquiring CBF estimates from alternative methods, estimating alternative dCA metrics, and incorporating dCA quantification into clinical trials. Implementation of these new and revised recommendations is important to improve the reliability and reproducibility of dCA studies, and to facilitate inter-institutional collaboration and the comparison of results between studies.

Extracting electrophysiological correlates of functional magnetic resonance imaging data using the canonical polyadic decomposition.

The relation between electrophysiology and BOLD-fMRI requires further elucidation. One approach for studying this relation is to find time-frequency features from electrophysiology that explain the variance of BOLD time-series. Convolution of these features with a canonical hemodynamic response function (HRF) is often required to model neurovascular coupling mechanisms and thus account for time shifts between electrophysiological and BOLD-fMRI data. We propose a framework for extracting the spatial distribution of these time-frequency features while also estimating more flexible, region-specific HRFs. The core component of this method is the decomposition of a tensor containing impulse response functions using the Canonical Polyadic Decomposition. The outputs of this decomposition provide insight into the relation between electrophysiology and BOLD-fMRI and can be used to construct estimates of BOLD time-series. We demonstrated the performance of this method on simulated data while also examining the effects of simulated measurement noise and physiological confounds. Afterwards, we validated our method on publicly available task-based and resting-state EEG-fMRI data. We adjusted our method to accommodate the multisubject nature of these datasets, enabling the investigation of inter-subject variability with regards to EEG-to-BOLD neurovascular coupling mechanisms. We thus also demonstrate how EEG features for modelling the BOLD signal differ across subjects.

Modeling the Hemodynamic Response Function Using EEG-fMRI Data During Eyes-Open Resting-State Conditions and Motor Task Execution.

Being able to accurately quantify the hemodynamic response function (HRF) that links the blood oxygen level dependent functional magnetic resonance imaging (BOLD-fMRI) signal to the underlying neural activity is important both for elucidating neurovascular coupling mechanisms and improving the accuracy of fMRI-based functional connectivity analyses. In particular, HRF estimation using BOLD-fMRI is challenging particularly in the case of resting-state data, due to the absence of information about the underlying neuronal dynamics. To this end, using simultaneously recorded electroencephalography (EEG) and fMRI data is a promising approach, as EEG provides a more direct measure of neural activations. In the present work, we employ simultaneous EEG-fMRI to investigate the regional characteristics of the HRF using measurements acquired during resting conditions. We propose a novel methodological approach based on combining distributed EEG source space reconstruction, which improves the spatial resolution of HRF estimation and using block-structured linear and nonlinear models, which enables us to simultaneously obtain HRF estimates and the contribution of different EEG frequency bands. Our results suggest that the dynamics of the resting-state BOLD signal can be sufficiently described using linear models and that the contribution of each band is region specific. Specifically, it was found that sensory-motor cortices exhibit positive HRF shapes, whereas the lateral occipital cortex and areas in the parietal cortex, such as the inferior and superior parietal lobule exhibit negative HRF shapes. To validate the proposed method, we repeated the analysis using simultaneous EEG-fMRI measurements acquired during execution of a unimanual hand-grip task. Our results reveal significant associations between BOLD signal variations and electrophysiological power fluctuations in the ipsilateral primary motor cortex, particularly for the EEG beta band, in agreement with previous studies in the literature.

Removal of Transcranial Alternating Current Stimulation EEG Artifacts Using Blind Source Separation and Wavelets.

GOAL: Transcranial alternating current stimulation (tACS) is a non-invasive technology for modulating brain activity, with significant potential for improving motor and cognitive functions. To investigate the effects of tACS, many studies have used electroencephalographic (EEG) data recorded during brain stimulation. However, the large artifacts induced by tACS make the analysis of tACS-EEG recordings challenging, which in turn has prevented the implementation of closed-loop brain stimulation schemes. Here, we propose a novel combination of blind source separation (BSS) and wavelets to achieve removal of tACS-EEG artifacts with improved performance. METHODS: We examined the performance of several BSS methods both applied individually, as well as combined with the empirical wavelet transform (EWT) in terms of denoising realistic simulated and experimental tACS-EEG data. RESULTS: EWT combined with BSS yielded considerably improved performance compared to BSS alone for both simulated and experimental data. Overall, independent vector analysis (IVA) combined with EWT yielded the best performance. SIGNIFICANCE: The proposed method yields promise for quantifying the effects of tACS on simultaneously recorded EEG data, which can in turn contribute towards understanding the effects of tACS on brain activity, as well as extracting reliable biomarkers that may be used to develop closed-loop tACS strategies for modulating the underlying brain activity in real time.

A Wavelet-Based Approach for Estimating Time-Varying Connectivity in Resting-State Functional Magnetic Resonance Imaging.

Introduction: The selection of an appropriate window size, window function, and functional connectivity (FC) metric in the sliding window method is not straightforward due to the absence of ground truth. Methods: A previously proposed wavelet-based method was accordingly adjusted for estimating time-varying FC (TVFC) and was applied to a large high-quality, low-motion dataset of 400 resting-state functional magnetic resonance imaging data. Specifically, the wavelet coherence magnitude and relative phase were averaged across wavelet (frequency) scales to yield TVFC and synchronization patterns. To assess whether the observed fluctuations in TVFC were statistically significant (dynamic FC [dFC]; the distinction between TVFC and dFC is intentional), surrogate data were generated using the multivariate phase randomization (MVPR) and multivariate autoregressive randomization (MVAR) methods to define the null hypothesis of dFC absence. Results: By averaging across all frequencies, core regions of the default mode network (DMN; medial prefrontal and posterior cingulate cortices, inferior parietal lobes, hippocampal formation) were found to exhibit dFC (test-retest reproducibility of 90%) and were also synchronized in activity (-15° ≤ phase ≤15°). When averaging across distinct frequency bands, the same dynamic connections were identified, with the majority of them identified in the frequency range (0.01, 0.198) Hz, though with lower test-retest reproducibility (<66%). Additional analysis suggested that MVPR method better preserved properties (p < 10-10), including time-averaged coherence, of the original data compared with MVAR approach. Conclusions: The wavelet-based approach identified dynamic associations between the core DMN regions with fewer choices in parameters, compared with sliding window method. Impact statement We employed a wavelet-based method, previously used in the literature, and proposed modifications to assess time-varying functional connectivity in resting-state functional magnetic resonance imaging. With this approach, dynamic connections within the default mode network were identified, involving the medial prefrontal and posterior cingulate cortices, inferior parietal lobes, and hippocampal formation, which were also highly consistent in test-retest analysis (test-retest reproducibility of 90%), without the need to select window size, window function, and functional connectivity metric as with the sliding window method, whereby no consensus on the appropriate choices of hyperparameters currently exists in the literature.

A multi-measure approach for assessing the performance of fMRI preprocessing strategies in resting-state functional connectivity.

It is well established that head motion and physiological processes (e.g. cardiac and breathing activity) should be taken into consideration when analyzing and interpreting results in fMRI studies. However, even though recent studies aimed to evaluate the performance of different preprocessing pipelines there is still no consensus on the optimal strategy. This is partly due to the fact that the quality control (QC) metrics used to evaluate differences in performance across pipelines have often yielded contradictory results. Furthermore, preprocessing techniques based on physiological recordings or data decomposition techniques (e.g. aCompCor) have not been comprehensively examined. Here, to address the aforementioned issues, we propose a framework that summarizes the scores from eight previously proposed and novel QC metrics to a reduced set of two QC metrics that reflect the signal-to-noise ratio and the reduction in motion artifacts and biases in the preprocessed fMRI data. Using this framework, we evaluate the performance of three commonly used practices on the quality of data: 1) Removal of nuisance regressors from fMRI data, 2) discarding motion-contaminated volumes (i.e., scrubbing) before regression, and 3) low-pass filtering the data and the nuisance regressors before their removal. Using resting-state fMRI data from the Human Connectome Project, we show that the scores of the examined QC metrics improve the most when the global signal (GS) and about 17% of principal components from white matter (WM) are removed from the data. Finally, we observe a small further improvement with low-pass filtering at 0.20 Hz and milder variants of WM denoising, but not with scrubbing.

Identification of Beta Oscillatory Patterns During a Hand Grip Motor Task: A Comparative Analysis pre- and post-Exercise.

Electroencephalography (EEG) based Movement-Related Beta Band Desynchronization (MRBD) within the beta frequency band (13 - 30Hz) is commonly observed during motor task execution, and it has been associated with motor task performance. More recently, transient burst-like events termed beta bursts have been identified as another potential biomarker of motor function. Previous studies have reported decreased MRBD magnitude induced by exercise. However, little is known in terms of the effects of high-intensity exercise on beta burst patterns. In the present work, we investigated the modulatory effects of exercise on different beta burst features prior to, during and post motor task execution. We found that exercise mainly affected burst duration and burst rate within the left motor cortex area (M1) that is contralateral to the moving hand. Meanwhile, burst amplitude in the contralateral M1 area was affected differently by exercise, with smaller burst amplitude values observed during the movement preparation phase and larger magnitude during as well as post motor task execution. Since MRBD and beta burst patterns are closely associated with motor task performance, results from the present study can promote understanding about the association between exercise induced neural plasticity changes and motor performance, which can be further used for designing a subject-specific training therapy for improving motor function.