Modeling brain dynamics in brain tumor patients using the virtual brain

Presurgical planning for brain tumor resection aims at delineating eloquent tissue in the vicinity of the lesion to spare during surgery. To this end, noninvasive neuroimaging techniques such as functional MRI and diffusion-weighted imaging fiber tracking are currently employed. However, taking into account this information is often still insufficient, as the complex nonlinear dynamics of the brain impede straightforward prediction of functional outcome after surgical intervention. Large-scale brain network modeling carries the potential to bridge this gap by integrating neuroimaging data with biophysically based models to predict collective brain dynamics. As a first step in this direction, an appropriate computational model has to be selected, after which suitable model parameter values have to be determined. To this end, we simulated large-scale brain dynamics in 25 human brain tumor patients and 11 human control participants using The Virtual Brain, an open-source neuroinformatics platform. Local and global model parameters of the Reduced Wong-Wang model were individually optimized and compared between brain tumor patients and control subjects. In addition, the relationship between model parameters and structural network topology and cognitive performance was assessed. Results showed (1) significantly improved prediction accuracy of individual functional connectivity when using individually optimized model parameters; (2) local model parameters that can differentiate between regions directly affected by a tumor, regions distant from a tumor, and regions in a healthy brain; and (3) interesting associations between individually optimized model parameters and structural network topology and cognitive performance

An automated method for identifying an independent component analysis-based language-related resting-state network in brain tumor subjects for surgical planning

As a noninvasive and "task-free" technique, resting-state functional magnetic resonance imaging (rs-fMRI) has been gradually applied to pre-surgical functional mapping. Independent component analysis (ICA)-based mapping has shown advantage, as no a priori information is required. We developed an automated method for identifying language network in brain tumor subjects using ICA on rs-fMRI. In addition to standard processing strategies, we applied a discriminability-index-based component identification algorithm to identify language networks in three different groups. The results from the training group were validated in an independent group of healthy human subjects. For the testing group, ICA and seed-based correlation were separately computed and the detected language networks were assessed by intra-operative stimulation mapping to verify reliability of application in the clinical setting. Individualized language network mapping could be automatically achieved for all subjects from the two healthy groups except one (19/20, success rate = 95.0%). In the testing group (brain tumor patients), the sensitivity of the language mapping result was 60.9%, which increased to 87.0% (superior to that of conventional seed-based correlation [47.8%]) after extending to a radius of 1 cm. We established an automatic and practical component identification method for rs-fMRI-based pre-surgical mapping and successfully applied it to brain tumor patients

Tumor Tissue Detection using Blood-Oxygen-Level-Dependent Functional MRI based on Independent Component Analysis

Accurate delineation of gliomas from the surrounding normal brain areas helps maximize tumor resection and improves outcome. Blood-oxygen-level-dependent (BOLD) functional MRI (fMRI) has been routinely adopted for presurgical mapping of the surrounding functional areas. For completely utilizing such imaging data, here we show the feasibility of using presurgical fMRI for tumor delineation. In particular, we introduce a novel method dedicated to tumor detection based on independent component analysis (ICA) of resting-state fMRI (rs-fMRI) with automatic tumor component identification. Multi-center rs-fMRI data of 32 glioma patients from three centers, plus the additional proof-of-concept data of 28 patients from the fourth center with non-brain musculoskeletal tumors, are fed into individual ICA with different total number of components (TNCs). The best-fitted tumor-related components derived from the optimized TNCs setting are automatically determined based on a new template-matching algorithm. The success rates are 100%, 100% and 93.75% for glioma tissue detection for the three centers, respectively, and 85.19% for musculoskeletal tumor detection. We propose that the high success rate could come from the previously overlooked ability of BOLD rs-fMRI in characterizing the abnormal vascularization, vasomotion and perfusion caused by tumors. Our findings suggest an additional usage of the rs-fMRI for comprehensive presurgical assessment

BOLD Coupling between Lesioned and Healthy Brain Is Associated with Glioma Patients’ Recovery

Predicting functional outcomes after surgery and early adjuvant treatment is difficult due to the complex, extended, interlocking brain networks that underpin cognition. The aim of this study was to test glioma functional interactions with the rest of the brain, thereby identifying the risk factors of cognitive recovery or deterioration. Seventeen patients with diffuse non-enhancing glioma (aged 22–56 years) were longitudinally MRI scanned and cognitively assessed before and after surgery and during a 12-month recovery period (55 MRI scans in total after exclusions). We initially found, and then replicated in an independent dataset, that the spatial correlation pattern between regional and global BOLD signals (also known as global signal topography) was associated with tumour occurrence. We then estimated the coupling between the BOLD signal from within the tumour and the signal extracted from different brain tissues. We observed that the normative global signal topography is reorganised in glioma patients during the recovery period. Moreover, we found that the BOLD signal within the tumour and lesioned brain was coupled with the global signal and that this coupling was associated with cognitive recovery. Nevertheless, patients did not show any apparent disruption of functional connectivity within canonical functional networks. Understanding how tumour infiltration and coupling are related to patients’ recovery represents a major step forward in prognostic development.Consejeria de Economia, Innovacion, Ciencia y Empleo.Junta de Andalucia CV20-45250; A-TIC-080-UGR18; B-TIC-586-UGR20; P20-0052

BOLD Coupling between Lesioned and Healthy Brain Is Associated with Glioma Patients’ Recovery

This article belongs to the Special Issue Perioperative Imaging and Mapping Methods in Glioma Patients.[Simple Summary] Glioma, a type of brain tumour, affects not only the function of immediately adjacent brain tissue but also that in more distant areas, potentially impacting cognitive function after its surgical removal. Here, 17 patients with glioma had brain scans and tests of cognitive function during treatment and recovery. We investigated the effects of glioma on the brain, and what happens during recovery, using the brain’s “global signal” detected with magnetic resonance imaging (MRI). We found that the signal from gliomas was synchronised with the global signal in all patients and that this synchronisation was associated with the recovery of cognition after surgery. Specifically, patients with a greater reduction in glioma–global signal synchronisation following surgery were more likely to have a larger number of newly acquired cognitive difficulties. Together, these results suggest that the interaction between gliomas and the brain can predict how patients recover their cognitive abilities, which is important for their quality of life.[Abstract] Predicting functional outcomes after surgery and early adjuvant treatment is difficult due to the complex, extended, interlocking brain networks that underpin cognition. The aim of this study was to test glioma functional interactions with the rest of the brain, thereby identifying the risk factors of cognitive recovery or deterioration. Seventeen patients with diffuse non-enhancing glioma (aged 22–56 years) were longitudinally MRI scanned and cognitively assessed before and after surgery and during a 12-month recovery period (55 MRI scans in total after exclusions). We initially found, and then replicated in an independent dataset, that the spatial correlation pattern between regional and global BOLD signals (also known as global signal topography) was associated with tumour occurrence. We then estimated the coupling between the BOLD signal from within the tumour and the signal extracted from different brain tissues. We observed that the normative global signal topography is reorganised in glioma patients during the recovery period. Moreover, we found that the BOLD signal within the tumour and lesioned brain was coupled with the global signal and that this coupling was associated with cognitive recovery. Nevertheless, patients did not show any apparent disruption of functional connectivity within canonical functional networks. Understanding how tumour infiltration and coupling are related to patients’ recovery represents a major step forward in prognostic development.This research was supported by the Guarantors of Brain, Cancer Research UK Cambridge Centre, The Brain Tumour Charity and the EMERGIA Junta de Andalucia program. Y.E. is funded by a Royal Society Dorothy Hodgkin Research Fellowship (DHF130100). JMG is funded by the Ministerio de Ciencia e Innovación (España)/FEDER under the RTI2018-098913-B100 project, by the Consejería de Economía, Innovación, Ciencia y Empleo (Junta de Andalucía) and FEDER under CV20-45250, A-TIC-080-UGR18, B-TIC-586-UGR20 and P20-00525 projects. MA was funded by a Cambridge Trust—Yousef Jameel Scholarship. This research was also supported by the NIHR Cambridge Biomedical Research Centre (BRC-1215-20014). SJP (NIHR Career Development Fellowship, CDF-2018-11-ST2-003) is funded by the National Institute for Health Research (NIHR) for this research project

Bold coupling between lesioned and healthy brain is associated with glioma patients’ recovery

Predicting functional outcomes after surgery and early adjuvant treatment is difficult due to the complex, extended, interlocking brain networks that underpin cognition. The aim of this study was to test glioma functional interactions with the rest of the brain, thereby identifying the risk factors of cognitive recovery or deterioration. Seventeen patients with diffuse non-enhancing glioma (aged 22–56 years) were longitudinally MRI scanned and cognitively assessed before and after surgery and during a 12-month recovery period (55 MRI scans in total after exclusions). We initially found, and then replicated in an independent dataset, that the spatial correlation pattern between regional and global BOLD signals (also known as global signal topography) was associated with tumour occurrence. We then estimated the coupling between the BOLD signal from within the tumour and the signal extracted from different brain tissues. We observed that the normative global signal topography is reorganised in glioma patients during the recovery period. Moreover, we found that the BOLD signal within the tumour and lesioned brain was coupled with the global signal and that this coupling was associated with cognitive recovery. Nevertheless, patients did not show any apparent disruption of functional connectivity within canonical functional networks. Understanding how tumour infiltration and coupling are related to patients’ recovery represents a major step forward in prognostic development.</p

Correlation between brain functional connectivity and neurocognitive function in patients with left frontal glioma

The association between neurocognitive function (NCF) impairment and brain cortical functional connectivity in glioma patients remains unclear. The correlations between brain oscillatory activity or functional connectivity and NCF measured by the Wechsler Adult Intelligence Scale full-scale intelligence quotient scores (WAIS FSIQ), the Wechsler Memory Scale-revised general memory scores (WMS-R GM), and the Western aphasia battery aphasia quotient scores (WAB AQ) were evaluated in 18 patients with left frontal glioma using resting-state electroencephalography (EEG). Current source density (CSD) and lagged phase synchronization (LPS) were analyzed using exact low-resolution electromagnetic tomography (eLORETA). Although 2 and 2 patients scored in the borderline range of WAIS FSIQ and WMS-R GM, respectively, the mean WAIS FSIQ, WMS-R GM, and WAB AQ values of all patients were within normal limits, and none had aphasia. In the correlation analysis, lower WMS-R GM was associated with a higher LPS value between the right anterior prefrontal cortex and the left superior parietal lobule in the beta1 band (13-20 Hz, R = - 0.802, P = 0.012). These findings suggest that LPS evaluated by scalp EEG is associated with memory function in patients with left frontal glioma and mild NCF disorders

Global Effects of Focal Brain Tumors on Functional Complexity and Network Robustness: A Prospective Cohort Study.

BACKGROUND: Neurosurgical management of brain tumors has entered a paradigm of supramarginal resections that demands thorough understanding of peritumoral functional effects. Historically, the effects of tumors have been believed to be local, and long-range effects have not been considered. OBJECTIVE: To test the hypothesis that tumors affect the brain globally, producing long-range gradients in cortical function. METHODS: Resting-state functional magnetic resonance imaging (fMRI) data were acquired from 11 participants with glioblastoma and split into discovery and validation datasets in a single-center prospective cohort study. Fractal complexity was computed with a wavelet-based estimator of the Hurst exponent. Distance-related effects of the tumors were tested with a tumor mask-dilation technique and parcellation of the underlying Hurst maps. RESULTS: Fractal complexity demonstrates a penumbra of suppression in the peritumoral region. At a global level, as distance from the tumor increases, this initial suppression is balanced by a subsequent overactivity before finally normalizing. These effects were best fit by a quadratic model and were consistent across different network construction pipelines. The Hurst exponent was correlated with graph theory measures of centrality including network robustness, but graph theory measures did not demonstrate distance-dependent effects. CONCLUSION: This work provides evidence supporting the theory that focal brain tumors produce long-range gradients in function. Consequently, the effects of focal lesions need to be interpreted in terms of the global changes on functional complexity and network architecture rather than purely in terms of functional localization. Determining whether peritumoral changes represent potential plasticity may facilitate extended resection of tumors without functional cost.MGH is funded by the Wellcome Trust Neuroscience in Psychiatry Network with additional support from the National Institute for Health Research Cambridge Biomedical Research Centre. The imaging studies were funded by an NIHR Clinician Scientist Fellowship for SJP (NIHR/CS/009/011)

Interrogating spatiotemporal patterns of resting state neuronal and hemodynamic activity in the awake mouse model

Since the advent of functional magnetic resonance imaging (fMRI) and the rise in popularity of its use for resting state functional connectivity mapping (rs-FCM) to non-invasively detect correlated networks of brain activity in human and animal models, many resting state FCM studies have reported differences in these networks under pathologies such as Alzheimer’s or schizophrenia, highlighting the potential for the method’s diagnostic relevance. A common underlying assumption of this analysis, however, is that the blood oxygen level dependent (BOLD) signal of fMRI is a direct measurement of local neural activity. The BOLD signal is in fact a measurement of the local changes in concentration of deoxy-hemoglobin (HbR). Thus, it is imperative that neurovascular coupling—the relationship between neuronal activity and subsequent hemodynamic activity—be better characterized to enable accurate interpretation of resting state fMRI in the context of clinical usage. This dissertation first describes the development and utility of WFOM paradigm for the robust and easily adaptable imaging of simultaneous neuronal and hemodynamic activity in awake mouse models of health or disease in strains with genetically encoded fluorescent calcium reporters. Subsequent exploration of resting state WFOM data collected in Thy1-GCaMP3 and Thy1-GCaMP6f mouse strains is then presented, namely the characterization of spatiotemporal patterns of neuronal and hemodynamic activity and different modulatory depths of neuronal activity via a toolbox of unsupervised blind source separation (e.g. k-means clustering) and supervised (e.g. non-negative least squares, Pearson correlation) analysis tools. The presence of these different modulatory depths of neuronal activity were then confirmed in another Thy1-jRGECO1a mouse strain using the same imaging scheme. Finally, the dissertation documents the application of the WFOM paradigm and select analysis tools to a novel mouse model of diffusely infiltrating glioma, through which neuronal and hemodynamic activity changes during diffusely infiltrating glioma development which impact temporal coherence of the tumor region activity relative to non-tumor regions activity were recorded and analyzed. The paradigm also allowed for recording of numerous spontaneous occurrences of interictal neuronal activity during which neurovascular coupling is modified in the tumor, as well as occurrences of non-convulsive generalized seizure activity (during which neurovascular is non-linear and cortex eventually suffers hypoxia). The detection of spatiotemporal patterns and different modulatory depths of activity in the awake mouse cortex, as well as observation of changes in functional activity in the context of diffusely infiltrating glioma, provide us with new insights into the possible mechanisms underlying variations in resting state connectivity networks found in resting state fMRI studies comparing health and disease states

Assessment of hemodynamic and connectivity alterations in brain gliomas through sparse DCM

The Master thesis is focused on the assessment of the alterations in the hemodynamic response function in 12 patients with brain gliomas through sparse Dynamic Causal Model, after it has been explored the brain connectivity through the Effective connectivity matrices for each subject.The Master thesis is focused on the assessment of the alterations in the hemodynamic response function in 12 patients with brain gliomas through sparse Dynamic Causal Model, after it has been explored the brain connectivity through the Effective connectivity matrices for each subject