Effect of tumor resection on the characteristics of functional brain networks

Network topology in brain tumor patients with and without structural epilepsy: a prospective MEG study

Background

It was proposed that network topology is altered in brain tumor patients. However, there is no consensus on the pattern of these changes and evidence on potential drivers is lacking.

Objectives

We aimed to characterize neurooncological patients' network topology by analyzing glial brain tumors (GBTs) and brain metastases (BMs) with respect to the presence of structural epilepsy.

Methods

Network topology derived from resting state magnetoencephalography was compared between (1) patients and controls, (2) GBTs and BMs, and (3) patients with (PSEs) and without structural epilepsy (PNSEs). Eligible patients were investigated from February 2019 to March 2021. We calculated whole brain (WB) connectivity in six frequency bands, network topological parameters (node degree, average shortest path length, local clustering coefficient) and performed a stratification, where differences in power were identified. For data analysis, we used Fieldtrip, Brain Connectivity MATLAB toolboxes, and in-house built scripts.

Results

We included 41 patients (21 men), with a mean age of 60.1 years (range 23-82), of those were: GBTs (n = 23), BMs (n = 14), and other histologies (n = 4). Statistical analysis revealed a significantly decreased WB node degree in patients versus controls in every frequency range at the corrected level (p1-30Hz = 0.002, pγ = 0.002, pβ = 0.002, pα = 0.002, pθ = 0.024, and pδ = 0.002). At the descriptive level, we found a significant augmentation for WB local clustering coefficient (p1-30Hz = 0.031, pδ = 0.013) in patients compared to controls, which did not persist the false discovery rate correction. No differences regarding networks of GBTs compared to BMs were identified. However, we found a significant increase in WB local clustering coefficient (pθ = 0.048) and decrease in WB node degree (pα = 0.039) in PSEs versus PNSEs at the uncorrected level.

Conclusion

Our data suggest that network topology is altered in brain tumor patients. Histology per se might not, however, tumor-related epilepsy seems to influence the brain's functional network. Longitudinal studies and analysis of possible confounders are required to substantiate these findings.

Graph Theoretical Analysis of Brain Network Characteristics in Brain Tumor Patients: A Systematic Review

Graph theory is a branch of mathematics that allows for the characterization of complex networks, and has rapidly grown in popularity in network neuroscience in recent years. Researchers have begun to use graph theory to describe the brain networks of individuals with brain tumors to shed light on disrupted networks. This systematic review summarizes the current literature on graph theoretical analysis of magnetic resonance imaging data in the brain tumor population with particular attention paid to treatment effects and other clinical factors. Included papers were published through June 24^th, 2020. Searches were conducted on Pubmed, PsycInfo, and Web of Science using the search terms (graph theory OR graph analysis) AND (brain tumor OR brain tumour OR brain neoplasm) AND (MRI OR EEG OR MEG). Studies were eligible for inclusion if they: evaluated participants with a primary brain tumor, used graph theoretical analyses on structural or functional MRI data, MEG, or EEG, were in English, and were an empirical research study. Seventeen papers met criteria for inclusion. Results suggest alterations in network properties are often found in people with brain tumors, although the directions of differences are inconsistent and few studies reported effect sizes. The most consistent finding suggests increased network segregation. Changes are most prominent with more intense treatment, in hub regions, and with factors such as faster tumor growth. The use of graph theory to study brain tumor patients is in its infancy, though some conclusions can be drawn. Future studies should focus on treatment factors, changes over time, and correlations with functional outcomes to better identify those in need of early intervention.

Overall survival time prediction for high-grade glioma patients based on large-scale brain functional networks

High-grade glioma (HGG) is a lethal cancer with poor outcome. Accurate preoperative overall survival (OS) time prediction for HGG patients is crucial for treatment planning. Traditional presurgical and noninvasive OS prediction studies have used radiomics features at the local lesion area based on the magnetic resonance images (MRI). However, the highly complex lesion MRI appearance may have large individual variability, which could impede accurate individualized OS prediction. In this paper, we propose a novel concept, namely brain connectomics-based OS prediction. It is based on presurgical resting-state functional MRI (rs-fMRI) and the non-local, large-scale brain functional networks where the global and systemic prognostic features rather than the local lesion appearance are used to predict OS. We propose that the connectomics features could capture tumor-induced network-level alterations that are associated with prognosis. We construct both low-order (by means of sparse representation with regional rs-fMRI signals) and high-order functional connectivity (FC) networks (characterizing more complex multi-regional relationship by synchronized dynamics FC time courses). Then, we conduct a graph-theoretic analysis on both networks for a jointly, machine-learning-based individualized OS prediction. Based on a preliminary dataset (N = 34 with bad OS, mean OS, ~400 days; N = 34 with good OS, mean OS, ~1030 days), we achieve a promising OS prediction accuracy (86.8%) on separating the individuals with bad OS from those with good OS. However, if using only conventionally derived descriptive features (e.g., age and tumor characteristics), the accuracy is low (63.2%). Our study highlights the importance of the rs-fMRI and brain functional connectomics for treatment planning.

Information diffusion backbones in temporal networks

Progress has been made in understanding how temporal network features affect the percentage of nodes reached by an information diffusion process. In this work, we explore further: which node pairs are likely to contribute to the actual diffusion of information, i.e., appear in a diffusion trajectory? How is this likelihood related to the local temporal connection features of the node pair? Such deep understanding of the role of node pairs is crucial to tackle challenging optimization problems such as which kind of node pairs or temporal contacts should be stimulated in order to maximize the prevalence of information spreading. We start by using Susceptible-Infected (SI) model, in which an infected (information possessing) node could spread the information to a susceptible node with a given infection probability β whenever a contact happens between the two nodes, as the information diffusion process. We consider a large number of real-world temporal networks. First, we propose the construction of an information diffusion backbone GB(β) for a SI spreading process with an infection probability β on a temporal network. The backbone is a weighted network where the weight of each node pair indicates how likely the node pair appears in a diffusion trajectory starting from an arbitrary node. Second, we investigate the relation between the backbones with different infection probabilities on a temporal network. We find that the backbone topology obtained for low and high infection probabilities approach the backbone GB(β → 0) and GB(β = 1), respectively. The backbone GB(β → 0) equals the integrated weighted network, where the weight of a node pair counts the total number of contacts in between. Finally, we explore node pairs with what local connection features tend to appear in GB(β = 1), thus actually contribute to the global information diffusion. We discover that a local connection feature among many other features we proposed, could well identify the (high-weight) links in GB(β = 1). This local feature encodes the time that each contact occurs, pointing out the importance of temporal features in determining the role of node pairs in a dynamic process.

The Direction of Tumour Growth in Glioblastoma Patients

Generating MR-derived growth pattern models for glioblastoma multiforme (GBM) has been an attractive approach in neuro-oncology, suggesting a distinct pattern of lesion spread with a tendency in growing along the white matter (WM) fibre direction for the invasive component. However, the direction of growth is not much studied in vivo. In this study, we sought to study the dominant directions of tumour expansion/shrinkage pre-treatment. We examined fifty-six GBMs at two time-points: at radiological diagnosis and as part of the pre-operative planning, both with contrast-enhanced T1-weighted MRIs. The tumour volumes were semi-automatically segmented. A non-linear registration resulting in a deformation field characterizing the changes between the two time points was used together with the segmented tumours to determine the dominant directions of tumour change. To compute the degree of alignment between tumour growth vectors and WM fibres, an angle map was calculated. Our results demonstrate that tumours tend to grow predominantly along the WM, as evidenced by the dominant vector population with the maximum alignments. Our findings represent a step forward in investigating the hypothesis that tumour cells tend to migrate preferentially along the WM.

Altered brain structural networks in attention deficit/hyperactivity disorder children revealed by cortical thickness

This study investigated the cortical thickness and topological features of human brain anatomical networks related to attention deficit/hyperactivity disorder. Data were collected from 40 attention deficit/hyperactivity disorder children and 40 normal control children. Interregional correlation matrices were established by calculating the correlations of cortical thickness between all pairs of cortical regions (68 regions) of the whole brain. Further thresholds were applied to create binary matrices to construct a series of undirected and unweighted graphs, and global, local, and nodal efficiencies were computed as a function of the network cost. These experimental results revealed abnormal cortical thickness and correlations in attention deficit/hyperactivity disorder, and showed that the brain structural networks of attention deficit/hyperactivity disorder subjects had inefficient small-world topological features. Furthermore, their topological properties were altered abnormally. In particular, decreased global efficiency combined with increased local efficiency in attention deficit/hyperactivity disorder children led to a disorder-related shift of the network topological structure toward regular networks. In addition, nodal efficiency, cortical thickness, and correlation analyses revealed that several brain regions were altered in attention deficit/hyperactivity disorder patients. These findings are in accordance with a hypothesis of dysfunctional integration and segregation of the brain in patients with attention deficit/hyperactivity disorder and provide further evidence of brain dysfunction in attention deficit/hyperactivity disorder patients by observing cortical thickness on magnetic resonance imaging.

Brain networks under attack: robustness properties and the impact of lesions

A growing number of studies approach the brain as a complex network, the so-called 'connectome'. Adopting this framework, we examine what types or extent of damage the brain can withstand-referred to as network 'robustness'-and conversely, which kind of distortions can be expected after brain lesions. To this end, we review computational lesion studies and empirical studies investigating network alterations in brain tumour, stroke and traumatic brain injury patients. Common to these three types of focal injury is that there is no unequivocal relationship between the anatomical lesion site and its topological characteristics within the brain network. Furthermore, large-scale network effects of these focal lesions are compared to those of a widely studied multifocal neurodegenerative disorder, Alzheimer's disease, in which central parts of the connectome are preferentially affected. Results indicate that human brain networks are remarkably resilient to different types of lesions, compared to other types of complex networks such as random or scale-free networks. However, lesion effects have been found to depend critically on the topological position of the lesion. In particular, damage to network hub regions-and especially those connecting different subnetworks-was found to cause the largest disturbances in network organization. Regardless of lesion location, evidence from empirical and computational lesion studies shows that lesions cause significant alterations in global network topology. The direction of these changes though remains to be elucidated. Encouragingly, both empirical and modelling studies have indicated that after focal damage, the connectome carries the potential to recover at least to some extent, with normalization of graph metrics being related to improved behavioural and cognitive functioning. To conclude, we highlight possible clinical implications of these findings, point out several methodological limitations that pertain to the study of brain diseases adopting a network approach, and provide suggestions for future research.

Emerging hubs in phantom perception connectomics

Brain networks are small-world networks typically characterized by the presence of hubs, i.e. nodes that have significantly greater number of links in comparison to other nodes in the network. These hubs act as short cuts in the network and promote long-distance connectivity. Long-distance connections increase the efficiency of information transfer but also increase the cost of the network. Brain disorders are associated with an altered brain connectome which reflects either as a complete change in the network topology, as in, the replacement of hubs or as an alteration in the connectivity between the hubs while retaining network structure. The current study compares the network topology of binary and weighted networks in tinnitus patients and healthy controls by studying the hubs of the two networks in different oscillatory bands. The EEG of 311 tinnitus patients and 256 control subjects are recorded, pre-processed and source-localized using sLORETA. The hubs of the different binary and weighted networks are identified using different measures of network centrality. The results suggest that the tinnitus and control networks are distinct in all the frequency bands but substantially overlap in the gamma frequency band. The differences in network topology in the tinnitus and control groups in the delta, theta and the higher beta bands are driven by a change in hubs as well as network connectivity; in the alpha band by changes in hubs alone and in the gamma band by changes in network connectivity. Thus the brain seems to employ different frequency band-dependent adaptive mechanisms trying to compensate for auditory deafferentation.

Time-dependent degree-degree correlations in epileptic brain networks: from assortative to dissortative mixing

We investigate the long-term evolution of degree-degree correlations (assortativity) in functional brain networks from epilepsy patients. Functional networks are derived from continuous multi-day, multi-channel electroencephalographic data, which capture a wide range of physiological and pathophysiological activities. In contrast to previous studies which all reported functional brain networks to be assortative on average, even in case of various neurological and neurodegenerative disorders, we observe large fluctuations in time-resolved degree-degree correlations ranging from assortative to dissortative mixing. Moreover, in some patients these fluctuations exhibit some periodic temporal structure which can be attributed, to a large extent, to daily rhythms. Relevant aspects of the epileptic process, particularly possible pre-seizure alterations, contribute marginally to the observed long-term fluctuations. Our findings suggest that physiological and pathophysiological activity may modify functional brain networks in a different and process-specific way. We evaluate factors that possibly influence the long-term evolution of degree-degree correlations.

The Union of Shortest Path Trees of Functional Brain Networks

Communication between brain regions is still insufficiently understood. Applying concepts from network science has shown to be successful in gaining insight in the functioning of the brain. Recent work has implicated that especially shortest paths in the structural brain network seem to play a major role in the communication within the brain. So far, for the functional brain network, only the average length of the shortest paths has been analyzed. In this article, we propose to construct the union of shortest path trees (USPT) as a new topology for the functional brain network. The minimum spanning tree, which has been successful in a lot of recent studies to comprise important features of the functional brain network, is always included in the USPT. After interpreting the link weights of the functional brain network as communication probabilities, the USPT of this network can be uniquely defined. Using data from magnetoencephalography, we applied the USPT as a method to find differences in the network topology of multiple sclerosis patients and healthy controls. The new concept of the USPT of the functional brain network also allows interesting interpretations and may represent the highways of the brain.

Neural network alterations underlie cognitive deficits in brain tumor patients

PURPOSE OF REVIEW

Brain tumor patients suffer from cognitive deficits, regardless of tumor grade or location. Deficits have a general character, falling in the domains of attention, working memory, information processing speed, and executive functioning. This review explores a new, brain network-based view of these deficits in brain tumor patients.

RECENT FINDINGS

Network theory has evolved within the fields of mathematics and sociology and has resulted in its application to many complex systems, such as social networks, traffic flow networks, and biological protein networks. In the brain, a network can be constructed by assessing either functional or anatomical connections between brain areas, and subsequently extracting their overarching network patterns. Important brain network features are local specialization, operationalized by local clustering, and global integration or path length. Widespread disturbances in network topology are found in brain tumor patients, which relate to their cognitive problems. Furthermore, changes in network topology in response to oncological interventions, particularly tumor resection, go hand in hand with cognitive outcome.

SUMMARY

Cognitive deficits in brain tumor patients are reflected in whole-brain network disturbances. Possible future clinical use of these findings mostly concerns prognostics and tailoring treatment strategies.

Small-World Brain Functional Networks in Children With Attention-Deficit/Hyperactivity Disorder Revealed by EEG Synchrony

We investigated the topologic properties of human brain attention-related functional networks associated with Multi-Source Interference Task (MSIT) performance using electroencephalography (EEG). Data were obtained from 13 children diagnosed with attention-deficit/hyperactivity disorder (ADHD) and 13 normal control children. Functional connectivity between all pairwise combinations of EEG channels was established by calculating synchronization likelihood (SL). The cluster coefficients and path lengths were computed as a function of degree K. The results showed that brain attention functional networks of normal control subjects had efficient small-world topologic properties, whereas these topologic properties were altered in ADHD. In particular, increased local characteristics combined with decreased global characteristics in ADHD led to a disorder-related shift of the network topologic structure toward ordered networks. These findings are consistent with a hypothesis of dysfunctional segregation and integration of the brain in ADHD, and enhance our understanding of the underlying pathophysiologic mechanism of this illness.

Incidental and intentional learning of verbal episodic material differentially modifies functional brain networks

Learning- and memory-related processes are thought to result from dynamic interactions in large-scale brain networks that include lateral and mesial structures of the temporal lobes. We investigate the impact of incidental and intentional learning of verbal episodic material on functional brain networks that we derive from scalp-EEG recorded continuously from 33 subjects during a neuropsychological test schedule. Analyzing the networks' global statistical properties we observe that intentional but not incidental learning leads to a significantly increased clustering coefficient, and the average shortest path length remains unaffected. Moreover, network modifications correlate with subsequent recall performance: the more pronounced the modifications of the clustering coefficient, the higher the recall performance. Our findings provide novel insights into the relationship between topological aspects of functional brain networks and higher cognitive functions.

Assortative mixing in functional brain networks during epileptic seizures

We investigate assortativity of functional brain networks before, during, and after one-hundred epileptic seizures with different anatomical onset locations. We construct binary functional networks from multi-channel electroencephalographic data recorded from 60 epilepsy patients; and from time-resolved estimates of the assortativity coefficient, we conclude that positive degree-degree correlations are inherent to seizure dynamics. While seizures evolve, an increasing assortativity indicates a segregation of the underlying functional network into groups of brain regions that are only sparsely interconnected, if at all. Interestingly, assortativity decreases already prior to seizure end. Together with previous observations of characteristic temporal evolutions of global statistical properties and synchronizability of epileptic brain networks, our findings may help to gain deeper insights into the complicated dynamics underlying generation, propagation, and termination of seizures.

Cognitive dysfunction in early multiple sclerosis: altered centrality derived from resting-state functional connectivity using magneto-encephalography

BACKGROUND

Cognitive dysfunction in multiple sclerosis (MS) is frequent. Insight into underlying mechanisms would help to develop therapeutic strategies.

OBJECTIVE

To explore the relationship of cognitive performance to patterns of nodal centrality derived from magneto-encephalography (MEG).

METHODS

34 early relapsing-remitting MS patients (median EDSS 2.0) and 28 age- and gender-matched healthy controls (HC) had a MEG, a neuropsychological assessment and structural MRI. Resting-state functional connectivity was determined by the synchronization likelihood. Eigenvector Centrality (EC) was used to quantify for each sensor its connectivity and importance within the network. A cognition-score was calculated, and normalized grey and white matter volumes were determined. EC was compared per sensor and frequency band between groups using permutation testing, and related to cognition.

RESULTS

Patients had lower grey and white matter volumes than HC, male patients lower cognitive performance than female patients. In HC, EC distribution showed highest nodal centrality over bi-parietal sensors ("hubs"). In patients, nodal centrality was even higher bi-parietally (theta-band) but markedly lower left temporally (upper alpha- and beta-band). Lower cognitive performance correlated to decreased nodal centrality over left temporal (lower alpha-band) and right temporal (beta-band) sensors, and to increased nodal centrality over right parieto-temporal sensors (beta-band). Network changes were most pronounced in male patients.

CONCLUSIONS

Partial functional disconnection of the temporal regions was associated with cognitive dysfunction in MS; increased centrality in parietal hubs may reflect a shift from temporal to possibly less efficient parietal processing. To better understand patterns and dynamics of these network changes, longitudinal studies are warranted, also addressing the influence of gender.

Altered small-world brain networks in schizophrenia patients during working memory performance

Impairment of working memory (WM) performance in schizophrenia patients (SZ) is well-established. Compared to healthy controls (HC), SZ patients show aberrant blood oxygen level dependent (BOLD) activations and disrupted functional connectivity during WM performance. In this study, we examined the small-world network metrics computed from functional magnetic resonance imaging (fMRI) data collected as 35 HC and 35 SZ performed a Sternberg Item Recognition Paradigm (SIRP) at three WM load levels. Functional connectivity networks were built by calculating the partial correlation on preprocessed time courses of BOLD signal between task-related brain regions of interest (ROIs) defined by group independent component analysis (ICA). The networks were then thresholded within the small-world regime, resulting in undirected binarized small-world networks at different working memory loads. Our results showed: 1) at the medium WM load level, the networks in SZ showed a lower clustering coefficient and less local efficiency compared with HC; 2) in SZ, most network measures altered significantly as the WM load level increased from low to medium and from medium to high, while the network metrics were relatively stable in HC at different WM loads; and 3) the altered structure at medium WM load in SZ was related to their performance during the task, with longer reaction time related to lower clustering coefficient and lower local efficiency. These findings suggest brain connectivity in patients with SZ was more diffuse and less strongly linked locally in functional network at intermediate level of WM when compared to HC. SZ show distinctly inefficient and variable network structures in response to WM load increase, comparing to stable highly clustered network topologies in HC.

Altered small-world brain networks in schizophrenia patients during working memory performance

Impairment of working memory (WM) performance in schizophrenia patients (SZ) is well-established. Compared to healthy controls (HC), SZ patients show aberrant blood oxygen level dependent (BOLD) activations and disrupted functional connectivity during WM performance. In this study, we examined the small-world network metrics computed from functional magnetic resonance imaging (fMRI) data collected as 35 HC and 35 SZ performed a Sternberg Item Recognition Paradigm (SIRP) at three WM load levels. Functional connectivity networks were built by calculating the partial correlation on preprocessed time courses of BOLD signal between task-related brain regions of interest (ROIs) defined by group independent component analysis (ICA). The networks were then thresholded within the small-world regime, resulting in undirected binarized small-world networks at different working memory loads. Our results showed: 1) at the medium WM load level, the networks in SZ showed a lower clustering coefficient and less local efficiency compared with HC; 2) in SZ, most network measures altered significantly as the WM load level increased from low to medium and from medium to high, while the network metrics were relatively stable in HC at different WM loads; and 3) the altered structure at medium WM load in SZ was related to their performance during the task, with longer reaction time related to lower clustering coefficient and lower local efficiency. These findings suggest brain connectivity in patients with SZ was more diffuse and less strongly linked locally in functional network at intermediate level of WM when compared to HC. SZ show distinctly inefficient and variable network structures in response to WM load increase, comparing to stable highly clustered network topologies in HC.

Identifying important nodes in weighted functional brain networks: a comparison of different centrality approaches

We compare different centrality metrics which aim at an identification of important nodes in complex networks. We investigate weighted functional brain networks derived from multichannel electroencephalograms recorded from 23 healthy subject under resting-state eyes-open or eyes-closed conditions. Although we observe the metrics strength, closeness, and betweenness centrality to be related to each other, they capture different spatial and temporal aspects of important nodes in these networks associated with behavioral changes. Identifying and characterizing of these nodes thus benefits from the application of several centrality metrics.

Altered small-world brain networks in temporal lobe in patients with schizophrenia performing an auditory oddball task

The functional architecture of the human brain has been extensively described in terms of complex networks characterized by efficient small-world features. Recent functional magnetic resonance imaging (fMRI) studies have found altered small-world topological properties of brain functional networks in patients with schizophrenia (SZ) during the resting state. However, little is known about the small-world properties of brain networks in the context of a task. In this study, we investigated the topological properties of human brain functional networks derived from fMRI during an auditory oddball (AOD) task. Data were obtained from 20 healthy controls and 20 SZ; A left and a right task-related network which consisted of the top activated voxels in temporal lobe of each hemisphere were analyzed separately. All voxels were detected by group independent component analysis. Connectivity of the left and right task-related networks were estimated by partial correlation analysis and thresholded to construct a set of undirected graphs. The small-worldness values were decreased in both hemispheres in SZ. In addition, SZ showed longer shortest path length and lower global efficiency only in the left task-related networks. These results suggested small-world attributes are altered during the AOD task-related networks in SZ which provided further evidences for brain dysfunction of connectivity in SZ.