PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2017; 12(3): e0173448.
Published online 2017 March 20. doi:  10.1371/journal.pone.0173448
PMCID: PMC5358752

The impact of high grade glial neoplasms on human cortical electrophysiology

Nader Pouratian, Editor

Abstract

Objective

The brain’s functional architecture of interconnected network-related oscillatory patterns in discrete cortical regions has been well established with functional magnetic resonance imaging (fMRI) studies or direct cortical electrophysiology from electrodes placed on the surface of the brain, or electrocorticography (ECoG). These resting state networks exhibit a robust functional architecture that persists through all stages of sleep and under anesthesia. While the stability of these networks provides a fundamental understanding of the organization of the brain, understanding how these regions can be perturbed is also critical in defining the brain’s ability to adapt while learning and recovering from injury.

Methods

Patients undergoing an awake craniotomy for resection of a tumor were studied as a unique model of an evolving injury to help define how the cortical physiology and the associated networks were altered by the presence of an invasive brain tumor.

Results

This study demonstrates that there is a distinct pattern of alteration of cortical physiology in the setting of a malignant glioma. These changes lead to a physiologic sequestration and progressive synaptic homogeneity suggesting that a de-learning phenomenon occurs within the tumoral tissue compared to its surroundings.

Significance

These findings provide insight into how the brain accommodates a region of “defunctionalized” cortex. Additionally, these findings may have important implications for emerging techniques in brain mapping using endogenous cortical physiology.

Introduction

The brain’s functional architecture of interconnected network-related oscillatory patterns in discrete cortical regions has been well established [1, 2]. The signal modalities identifying these networked regions have largely been either resting-state metabolic activity studied with functional magnetic resonance imaging (fMRI) studies [18] or direct cortical electrophysiology from electrodes placed on the surface of the brain, or electrocorticography (ECoG) [2, 912]. One of the earliest and most robust links between cortical electrophysiological signals and the radiographic blood oxygenation level dependent (BOLD) signal are those physiologic signals within the infraslow frequency band (<0.5 Hz) or the slow cortical potential (SCP) [2]. These resting state SCP networks exhibit a robust functional architecture that persists through all stages of sleep and under anesthesia and represent the physiologic correlate of the BOLD imaging signal [2, 13]. While the stability of these networks provides a fundamental understanding of the organization of the brain, understanding how these regions can be perturbed or altered is also critical in defining the brain’s ability to adapt while learning and recovering from injury.

Deliberate training and learning appears to strengthen resting state connectivity between brain regions that are mutually engaged in a learning task [14]. FMRI studies in animal models of stroke and in human patients recovering from cerebrovascular insults also have shown evidence to suggest that rehabilitation can alter resting state connectivity after injury [1517]. These findings highlight that these networks, while quite durable through multiple states of consciousness, are not immutable. Namely, that with intensive stimuli they can be subject to augmentation (in the case of learning) or restoration (in the context of injury). However, how the brain responds to a defunctionalized region is incompletely understood. The manner in which resting state functional architecture is altered in the context of an invasive brain tumor has been explored in a small number of studies using magnetoencephalography[18, 19] and more recently using resting state functional MRI[20], however this question has never been investigated using invasively acquired ECoG. The limited experience investigating the impact of an invasive glioma on resting state networks is referable in part to technical barriers with functional imaging in the context of an invasive lesion, specifically the decoupling effect these lesions exert on the neurovascular relationship which can make interpretation of their imaging more challenging [2023]. While invasive cortical physiology has been used to identify and study resting state networks associated with functional networks such as motor and language [9, 24, 25], this modality has not been used in the context of a pathologic lesion.

To address this gap, patients undergoing awake craniotomies for resection of a tumor were studied as a unique model to define how the cortical physiology and the associated networks were altered by the presence of an invasive brain tumor. Distinct from acute injuries such as a stroke, which immediately and radically reduce the function of a region of brain, tumors that invade the brain can provide insights into how the cortex accommodates to a subregion becoming non-functional. In this subacute injury model, we hypothesized that in the setting of an infiltrative glioblastoma, resting state connectivity would be perturbed in the cortex. This study demonstrates that there is a distinct pattern of alteration of cortical physiology in the setting of a malignant glioma. These changes lead to a physiologic sequestration and progressive synaptic homogeneity suggesting that a de-learning phenomenon occurs within the tumoral tissue compared to its surroundings.

Materials and methods

Subjects

Six patients undergoing surgical treatment for glioblastoma were included in this study. Study protocol was approved by the Human Research and Protection Organization at Washington University School of Medicine. Before inclusion, all patients provided written informed consent for study participation. Exclusion criteria included the absence of a glioblastoma on pathologic analysis of neoplastic tissue. See Table 1 for demographic and clinical information. Each patient underwent craniotomy for the resection of their tumor. See Figs Figs11 and and22 for intraoperative photo and preoperative imaging respectively. Following acquisition of a baseline ECoG signal, surgical resection of the patient’s neoplastic tissue proceeded in the usual fashion.

Fig 1
Intraoperative Photograph.
Fig 2
Preoperative MR Imaging.
Table 1
Patient Demographics & Clinical Information.

Data collection

With the patient sedated with remifentanil and midazolam for surgical exposure prior to planned awake craniotomy, the cortical surface electrode array was connected to amplifiers and an additional 1x4 cortical electrode strip was placed to provide the reference and ground. Five minutes of baseline cortical signal was acquired with the patient sedated prior to wake-up for planned awake tumor resection. Electrode grid placement was selected based on the anatomical relationship between the tumor and local functional cortex.

ECoG signals were acquired from the electrode array (Ad-Tech Medical Instrument Corporation, Racine WI, USA) at a sampling frequency of 1,200 Hz using g.tech biosignal amplifiers (Guger Technologies, Graz, Austria). The data were stored on a Dell PC running BCI2000 software [26]. The ECoG electrode array contained either 48 or 64 2.3mm diameter platinum electrodes spaced 10mm apart. Patient data is available in Supporting Information, S1S12 Files.

Electrode selection

Electrodes were defined as overlying tumor, peritumor or distant cortex based on analysis of both intraoperative photography and the patient’s preoperative MRI. In 4 patients, grossly neoplastic tissue presented to the cortical surface. In these cases, intraoperative photography was used exclusively to define electrodes overlying tumoral tissue. Peritumoral tissue was defined as cortex underlying perimeter electrodes immediately adjacent to tumor electrodes located within 1 cm of tumoral tissue. Distant electrodes were defined as all remaining electrodes. In the case of the 2 patients whose tumor did not present to the cortical surface, intraoperative photography was used to define the orientation of the grid relative to the exposed cortical surface. This was then related to the patient’s preoperative MRI to define the margins of the tumor relative to the exposed cortex. Electrodes were again defined as overlying tumoral tissue, peritumoral cortex or distant cortex in the same way. Electrodes determined to be overlying peritumor cortex were excluded from further analysis due to their position as a relative transition zone between the tumor and distant cortex.

Signal analysis

Signals from every electrode were visually inspected and those electrodes identified as having poor signal to noise characteristics (amplitude greater than 10x that of the majority of electrodes in the grid) were excluded from further analysis (33 of 304 electrodes, 11%). Remaining signals were re-referenced to the amplifier-specific common mean for further analysis. The 5-minute acquisition was analyzed as a single time epoch.

Power analysis was performed by Welch’s Method [27]. The time-varying power was computed for frequencies ranging from 0.5 to 150 Hz. Known noise bands at 60 and 120 Hz were avoided by band-pass filtering 2 Hz bands centered on each known noise peak. Analysis was cut off at 150 Hz as determined by the signal’s noise floor. Analysis of trends in the absolute band power and the variance in absolute band power followed for the subdivided “tumor” and “distant” signals. Findings between these two source signals were then compared. Power spectral curves were then fit to a power-law function by converting frequency and power to logarithmic space and fitting a linear regression curve. ECoG data characteristically has a knee in the low-frequency alpha/beta range and therefore the power-law structure was fit after this observed knee [28]. The error was determined by summation of squares of difference between the fit and observed curves. The ratio of these sum-of-square errors with observed variance was calculated to produce R2 values, which represent the percentage of variability explained by the power-law spectral model regression fit (Table 2). The R2 values quantify a goodness-of-fit supporting maintenance of the power-law structure in both tumor and distant electrodes.

Table 2
Power-law spectral fit.

Networked interactions were also evaluated. Signals from every electrode were first visually inspected and those electrodes identified as noisy (as above) were excluded from further analysis. Signals from the remaining electrodes were re-referenced to the global mean and low-pass filtered using a 3rd order lowpass digital Butterworth filter (cutoff = 0.5Hz) to isolate the SCP. The temporal correlation of the SCP over the five-minute epoch was then calculated for all electrode pairs by first computing the covariance (Eq 1) between each pair.

Covariance(m,n)=iN(electrodeimelectrodem¯)(electrodeinelectroden¯)N
(1)

Here, electrodemi and electrodeni represent the ith of N values of the SCP signal from the mth and nth electrodes, respectively. The correlation coefficient between each electrode pair was then obtained by Eq (2).

CorrelationCoefficient(m,n)=Covariance(m,n)Covariance(m,m)Covariance(n,n)
(2)

Electrical potentials measured from the surface of the brain have a distance-dependent inherent correlation strength[29, 30]. These local spatial correlations tend to be stronger in lower frequencies than higher frequencies[31]. Here, we investigate low frequency correlations within relatively confined regions defined by the presence or absence of a tumor. To control for this effect of local spatial correlation, we employed a method of distance regression of the correlation values between electrode pairs. Prior literature suggests the relationship is not entirely linear[31], therefore, we fit a spline regression curve to the distance and correlation coefficients of electrode pairs. The fitted regression was then subtracted from each correlation value and this corrected value, detrended for distance, was used for further analysis.

Next, ECoG signal was subdivided into those signals recorded from “tumor” electrodes and signals recorded from “distant” electrodes. For each electrode grouping, tumor and distant, a central, non-noisy electrode channel was selected as a representative seed electrode. First, we investigated mean ECoG correlation strength among tumor electrodes compared to distant electrodes to evaluate for alteration of signaling capacity within tumor-infiltrated brain. In this comparison, positive and negative correlations were considered to have equal meaning with respect to the strength of shared information, therefore, the absolute value of the correlation coefficient was used. The correlation between the seed electrode and all other electrodes within each group was calculated for each patient. These correlation values were then averaged across all electrode pairs within each group and compared using a two-sample t-test for each patient. The Bonferroni correction for number of patients was applied to correct for multiple comparisons.

Next, we tested the distribution of correlations for each region. Pairs of tumor electrodes and distant electrodes were evaluated as a group to determine if they represented samples of significantly different distributions. The Kolmogorov-Smirnov test of difference in sample distributions was used for this purpose. To correct for fewer electrodes overlying tumor compared to distant cortical regions in each subject, we employed a permutation resampling method to achieve equal sample size within each group. Correlations where binned by distance between electrode pairs for each group. Observed samples where then resampled for each group to achieve equal numbers of channel pairs with respect to binned distance. The distribution of correlations was normalized for each group and reported by frequency count area-under-the-curve (AUC).

Finally, topographical correlation maps relative to the seed electrode were computed for each patient and each electrode group. The Pearson correlation coefficient of each electrode with the seed electrode was computed and plotted in a grid configuration matching that of the ECoG electrode array. This was repeated for all pairs of electrodes with the tumor seed and distant seed, respectively. The resulting correlation maps provide a measure of topography of connectivity within the tumor and distant regions. To quantify the relationship between tumor and distant topographies, we calculated spatial correlation coefficients between tumor and distant topography maps for each patient. The spatial correlation is performed similar to the temporal correlation described in Eq 1, iterating over electrodes at successive locations instead of ECoG samples at successive time points. Positive correlations represent strength of similarity in topography, while negative correlations represent strength of reciprocal topography, such as elevations in one map corresponding to depression in the other.

Results

Alterations in power spectra

Global power was analyzed from 0.5–150 Hz. A consistent reduction in overall power was identified for tumor electrodes when compared against distant electrodes. Statistical significance across subjects as indicated by 95% confidence intervals was reached in 5 of 6 subjects (Fig 3). The remaining subject’s power spectra appeared consistent with those reaching statistical significance but noise inherent within the intraoperative setting limited the ability to achieve statistical significance. This power reduction was seen across all frequency bands suggesting maintenance of power-law scaling in these signals despite an overall reduction in amplitude. Thus, while tumors lead to a decrease in global cortical activity, they do not result in a fundamental change in the nature of the scaling relationship.

Fig 3
Power spectra comparing tumor (red) and distant (blue) spectral densities from each individual subject.

Maintenance of connectivity

Cortical connectivity, defined by correlation of SCP oscillations, was maintained when comparing time series data from electrodes overlying neoplastic versus distant cortical tissue. Correlation strength among tumor electrodes did not demonstrate a decrease among pairs of tumor electrodes relative to distant electrodes as would be expected if SCP connectivity was disrupted by the tumor (Fig 4A). Likewise, comparison between the distributions of correlation values across electrode pairs for both tumor and distant regions did not demonstrate a statistically significant difference as determined by a two-sample Kolmogorov-Smirnov test, p>0.05. This supports the notion that SCP connectivity is not reduced by the presence of a tumor (Fig 4B). It should be noted that the effect of regressing distance from correlations may artificially elevate correlation coefficients above 1.0 or below -1.0 due to intrinsically strong correlations despite a large distance. This corrected correlation value maintains the same properties of a standard correlation coefficient, namely being a dimensionless quantity that provides a relative measure of correlation strength and can be interpreted as such despite slight deviations from the standard, expected range of correlation values. The tumor and distant group correlation distributions are presented in the corrected form as correlation values which have been detrended for distance (Fig 4B).

Fig 4
A. Bar plot demonstrating comparison between average correlation values between all tumor electrode pairs (red) and all distant electrode pairs (blue). This suggests that connectivity is maintained within cortical regions invaded by glioma. B. Bar histogram ...

Disruption of local topography

Despite maintenance of connectivity overall, there is a perturbation in the cortical surface topography of correlation values relative to neoplastic tissue compared to non-neoplastic tissue. Said another way, the topographic structure of tumor and distant cortical sites are similar in distribution but separate in morphology. Spatial correlation comparisons between surface topographies for both tumor and distant seed electrodes demonstrate an internally consistent, positively correlated, topographic structure within the tumor compared to its surroundings that is absent when selecting for a distant electrode (Fig 5). Five of the six patients were found to have statistically significant reciprocal topographies between tumor and distant seed maps by spatial correlation analysis. This suggests that areas infiltrated by tumor show correlated activity that is separable from surrounding cortical activity. For alpha = 0.05 / 6 = 0.00833 (Bonferroni corrected for N = 6):

Fig 5
Patient specific correlation matrices demonstrating representative correlation maps with seeds placed centrally within each patient’s tumor (left column) and remotely within distant cortex (right column) for each individual subject.

Patients A, B and C significant at p < 0.001

Patients D and E significant at p < 0.005

Patient F did not show any statistically significant relationship between tumor and distant seed maps.

Discussion

Using invasively acquired electrophysiologic signal from human subjects, this study characterizes the impact of a high-grade primary glial neoplasm on cortical physiologic phenomena. The physiologic behavior of electrodes overlying neoplastic tissue was compared to that of electrodes overlying distant, normal tissue. Tumor infiltrated tissue appeared to have a reduced global power in the majority of patients, while maintaining an infra-slow scale network structure. Interestingly, this maintained network structure appeared predominantly limited to the region involved with the tumor. Taken together, these findings support the possibility of a “de-learning” phenomenon in which poorly functioning tissue alters its connectivity patterns to which it is ascribed prior to development of the neoplastic process and this becomes progressively sequestered.

Physiologic implications

The preservation of connectivity within the tumor as demonstrated by an unchanged mean correlation value between electrodes related to the underlying neoplasm and those related to distant cortical tissue indicates that the capability for connectivity is maintained at a cellular level despite the aggressive and infiltrative behavior of the underlying high grade glial neoplastic process. The negative-shift of these SCPs likely reflects the result of slow rhythmic depolarization of apical dendrites in superficial cortical layers [32, 33]. Thus, despite the tumor infiltrating more deeply in the white matter and even transgressing the gray matter the more superficial mechanisms of synaptic homeostasis seem to be preserved. This is in keeping with findings in peritumoral regions in a murine glioma model in which neurons are found to demonstrate functional changes including increased spontaneous activity which was reversible with administration of cytotoxic necrotizing factor 1 suggesting reversibility to these changes in the peritumoral area compared to the tumor core[34]. That said, on the topographic level it appears that there is an alteration in network structure that appears to be constrained to the region of the tumor as defined by covariance analysis of the ECoG signal. Namely, that the resting state network of internal intratumoral connectivity, is limited to the general area of the tumor itself.

On the spectral scale, there also are perturbations. There is the global power reduction within the tumor when compared to its surroundings. This reduction in overall power across frequency bands suggests that both thalamo-cortical and cortico-cortical processing is limited within the neoplastic process when compared to its surroundings [3537]. In addition to this global power reduction there is also a subtle flattening of the μβ peak in patients D & F. Mu (8-12Hz) and beta (18-26Hz) rhythms are oscillations that are thought to represent post-synaptic potentials associated with thalamocortical modulation of cortex [38]. These rhythms tend to be spectrally confined and anatomically broad. Their absence in power in electrodes overlying tumor-infiltrated regions suggests that impairment in the thalamic interactions with more superficial cortical regions has occurred. Given that gliomas are often found primarily in the white matter and likely to impact cortical connections with deeper structures, the loss in narrow band power peaks also provides “lesional” support for these mu and beta rhythms reflecting a disconnection between the thalamus and cortex. In the absence of the narrow band spectral oscillator peaks the residual signal is a lower power broad band signal that follows a power law distribution. This likely reflects residual asynchronous neuronal firing in cortex [28].

Taken together, the networked topography and the spectral analysis show the region of brain invaded by the tumor to be relatively homogeneous. While it is important to note that the physiologic measurement was taken at a single point in time (immediately prior to tumor resection), that the infiltrative nature of the glial neoplastic process may lead to progressive regional functional decline which manifests as intratumoral physiologic isolation is intriguing to consider. This is in keeping with the graded changes seen in resting state fMRI data in the presence of a glioblastoma where the tumor and contralateral internal control timeseries were significantly different from each other with the timeseries of the non-enhancing peritumoral tissue more closely resembling that of the tumor than that of the contralateral internal control[20]. This physiologic sequestration notion supports a de-learning phenomenon occurring within the tumor. Previous reports indicate that a degree of plasticity exists within SCP network connectivity [14, 39]. The corollary of that phenomenon appears to be functioning here in which the tumor creates a de-learning process. In the same way that repetitive task exposure and execution can strengthen SCP connectivity between two regions, the neoplastic invasion appears to reduce overall activity (e.g. global power reduction), reduce thalamocortical interaction (loss of mu-beta peaks), and reduce SCP connectivity between the tumor and distant regions. Cumulatively, this leads to maintained intratumoral connectivity and functional sequestration.

Clinical implications

Approximately 20,000 new primary glial neoplasms are diagnosed each year [40]. Of these, roughly half are aggressive, grade IV glioblastomas [41]. Overwhelmingly, the most common locations for these tumors to occur are within the frontal, temporal and parietal lobes [42]. These cortical regions largely comprise functional cortex involved in both motor and language systems, among others. Many of these tumors are in close proximity to these functional cortical regions. This situation often requires employing mapping techniques during their resection in order to maintain safe distance from adjacent functional cortex. Burgeoning techniques for mapping functional areas are currently emerging that use endogenous electrophysiology either during a task or at rest to complement more classic methodologies of electro-cortical stimulation (ECS)[4346]. All the early demonstrations of these techniques applied during the resting state, however, have used invasively monitored epilepsy subjects without distorting brain lesions.

For task-based applications, amplitude modulation in higher gamma frequency bands has been used to identify focal areas associated with cortical activation. Activity in the high gamma frequency range (> 40 Hz) are thought to be produced by local cortical circuits [47]. Because of the anatomically constrained nature of high frequency amplitude modulation, these gamma rhythms have been focused on as the optimal signal for use in task-based clinical brain mapping [48, 49]. In the context of a brain tumor and the demonstration of reduced power across all frequencies it will be important to test if the sensitivity and specificity determined in non-lesional human cortex will still be maintained. Task independent mapping techniques have also been employed [50]. As mentioned, the cortical physiology that most correlates with fMRI-defined resting state networks are covariant regions defined by very low frequency oscillations (< 0.5Hz), known as slow cortical potentials. [51] Given that these networks are durable through sleep wake cycles, and anesthesia [5254], they were evaluated as a possible tool for localizing eloquent cortex. Breshears et al., demonstrated that these SCP defined networks identified sensorimotor cortex using a data-driven approach in patients under anesthesia and awake with a high sensitivity, suggesting that resting-state networks may be useful for tailoring subsequent stimulation mapping. This must be guided however by the knowledge that BOLD activations are reduced in the setting of an infiltrating glioma [55, 56]. If the neoplastic tissue is anatomically sequestered from surrounding, non-neoplastic cortex, it remains an open question as to whether it would be safe to resect.

Limitations

The limitations of this study include a relatively small sample size of six patients. A second limitation was limited electrode grid coverage. Each patient’s grid coverage was dictated exclusively by clinical circumstances alone, therefore our findings should only be interpreted in the context of limited electrode grid coverage. Finally, our findings are applicable exclusively to high-grade glial neoplasms. As such, the consistency of cortical physiologic phenomena in lower grade glial neoplasms and metastatic disease still remains unknown.

Conclusions

We report a unique assessment of invasive human cortical physiology related to glioblastoma. Observations from our analysis include the maintenance of the capability for connectivity within the tumor but a physiologic sequestration of the tumor from the surrounding cortex coupled with an overall reduction in global power of these regions. These findings suggest a possible de-learning phenomenon that may result in a functional isolation of the region of brain invaded by the tumor. In addition to providing insight into how the brain accommodates a region of “defunctionalized” cortex, these findings may have important clinical implications including applications for emerging techniques in brain mapping using endogenous cortical physiology.

Supporting information

S1 File

PatientAdata.mat.

Patient A data file.

(MAT)

S2 File

PatientAparams.mat.

Patient A acquisition parameters.

(MAT)

S3 File

PatientBdata.mat.

Patient B data file.

(MAT)

S4 File

PatientBparams.mat.

Patient B acquisition parameters.

(MAT)

S5 File

PatientCdata.mat.

Patient C data file.

(MAT)

S6 File

PatientCparams.mat.

Patient C acquisition parameters.

(MAT)

S7 File

PatientDdata.mat.

Patient D data file.

(MAT)

S8 File

PatientDparams.mat.

Patient D acquisition parameters.

(MAT)

S9 File

PatientEdata.mat.

Patient E data file.

(MAT)

S10 File

PatientEparams.mat.

Patient E acquisition parameters.

(MAT)

S11 File

PatientFdata.mat.

Patient F data file.

(MAT)

S12 File

PatientFparams.mat.

Patient F acquisition parameters.

(MAT)

Abbreviations

ECoG
electrocorticography
MRI
magnetic resonance imaging
SCP
slow cortical potential
Hz
hertz
mm
millimeters
BOLD
blood oxygen level dependent
ROI
region of interest

Funding Statement

The authors received no specific funding for this work.

Data Availability

Data Availability

Data files have been uploaded as Supporting Information.

References

1. Biswal B, Yetkin FZ, Haughton VM, Hyde JS. Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn Reson Med. 1995;34(4):537–41. Epub 1995/10/01. [PubMed]
2. He BJ, Snyder AZ, Zempel JM, Smyth MD, Raichle ME. Electrophysiological correlates of the brain's intrinsic large-scale functional architecture. Proc Natl Acad Sci U S A. 2008;105(41):16039–44. Epub 2008/10/10. PubMed Central PMCID: PMC2564983. doi: 10.1073/pnas.0807010105 [PubMed]
3. Cordes D, Haughton VM, Arfanakis K, Wendt GJ, Turski PA, Moritz CH, et al. Mapping functionally related regions of brain with functional connectivity MR imaging. AJNR Am J Neuroradiol. 2000;21(9):1636–44. [PubMed]
4. Arfanakis K, Cordes D, Haughton VM, Moritz CH, Quigley MA, Meyerand ME. Combining independent component analysis and correlation analysis to probe interregional connectivity in fMRI task activation datasets. Magn Reson Imaging. 2000;18(8):921–30. [PubMed]
5. Hampson M, Peterson BS, Skudlarski P, Gatenby JC, Gore JC. Detection of functional connectivity using temporal correlations in MR images. Hum Brain Mapp. 2002;15(4):247–62. [PubMed]
6. Greicius MD, Krasnow B, Reiss AL, Menon V. Functional connectivity in the resting brain: a network analysis of the default mode hypothesis. Proc Natl Acad Sci U S A. 2003;100(1):253–8. PubMed Central PMCID: PMCPMC140943. doi: 10.1073/pnas.0135058100 [PubMed]
7. Di X, Biswal BB. Dynamic brain functional connectivity modulated by resting-state networks. Brain Struct Funct. 2013. Epub 2013/10/01. [PMC free article] [PubMed]
8. Tertel K, Tandon N, Ellmore TM. Probing brain connectivity by combined analysis of diffusion MRI tractography and electrocorticography. Comput Biol Med. 2011;41(12):1092–9. Epub 2010/12/07. PubMed Central PMCID: PMC3085041. doi: 10.1016/j.compbiomed.2010.11.004 [PMC free article] [PubMed]
9. Greenlee JD, Behroozmand R, Larson CR, Jackson AW, Chen F, Hansen DR, et al. Sensory-motor interactions for vocal pitch monitoring in non-primary human auditory cortex. PloS one. 2013;8(4):e60783 PubMed Central PMCID: PMC3620048. doi: 10.1371/journal.pone.0060783 [PMC free article] [PubMed]
10. Ko AL, Weaver KE, Hakimian S, Ojemann JG. Identifying functional networks using endogenous connectivity in gamma band electrocorticography. Brain Connect. 2013;3(5):491–502. Epub 2013/07/25. PubMed Central PMCID: PMC3796331. doi: 10.1089/brain.2013.0157 [PMC free article] [PubMed]
11. Gunduz A, Brunner P, Daitch A, Leuthardt EC, Ritaccio AL, Pesaran B, et al. Neural correlates of visual-spatial attention in electrocorticographic signals in humans. Front Hum Neurosci. 2011;5:89 Epub 2011/11/03. PubMed Central PMCID: PMC3202224. doi: 10.3389/fnhum.2011.00089 [PMC free article] [PubMed]
12. Crone NE, Sinai A, Korzeniewska A. High-frequency gamma oscillations and human brain mapping with electrocorticography. Prog Brain Res. 2006;159:275–95. Epub 2006/10/31. doi: 10.1016/S0079-6123(06)59019-3 [PubMed]
13. Breshears JD, Roland JL, Sharma M, Gaona CM, Freudenburg ZV, Tempelhoff R, et al. Stable and dynamic cortical electrophysiology of induction and emergence with propofol anesthesia. Proc Natl Acad Sci U S A. 2010;107(49):21170–5. Epub 2010/11/17. PubMed Central PMCID: PMC3000270. doi: 10.1073/pnas.1011949107 [PubMed]
14. Lewis CM, Baldassarre A, Committeri G, Romani GL, Corbetta M. Learning sculpts the spontaneous activity of the resting human brain. Proc Natl Acad Sci U S A. 2009;106(41):17558–63. Epub 2009/10/07. PubMed Central PMCID: PMC2762683. doi: 10.1073/pnas.0902455106 [PubMed]
15. Borstad AL, Bird T, Choi S, Goodman L, Schmalbrock P, Nichols-Larsen DS. Sensorimotor training and neural reorganization after stroke: a case series. J Neurol Phys Ther. 2013;37(1):27–36. Epub 2013/02/13. doi: 10.1097/NPT.0b013e318283de0d [PMC free article] [PubMed]
16. Jiang L, Xu H, Yu C. Brain connectivity plasticity in the motor network after ischemic stroke. Neural Plast. 2013;2013:924192 Epub 2013/06/06. PubMed Central PMCID: PMC3655657. doi: 10.1155/2013/924192 [PMC free article] [PubMed]
17. Grefkes C, Ward NS. Cortical Reorganization After Stroke: How Much and How Functional? Neuroscientist. 2013. Epub 2013/06/19. [PubMed]
18. van Dellen E, Douw L, Hillebrand A, Ris-Hilgersom IH, Schoonheim MM, Baayen JC, et al. MEG network differences between low- and high-grade glioma related to epilepsy and cognition. PloS one. 2012;7(11):e50122 PubMed Central PMCID: PMCPMC3498183. doi: 10.1371/journal.pone.0050122 [PMC free article] [PubMed]
19. Bosma I, Stam CJ, Douw L, Bartolomei F, Heimans JJ, van Dijk BW, et al. The influence of low-grade glioma on resting state oscillatory brain activity: a magnetoencephalography study. J Neurooncol. 2008;88(1):77–85. PubMed Central PMCID: PMCPMC2295256. doi: 10.1007/s11060-008-9535-3 [PMC free article] [PubMed]
20. Chow DS, Horenstein CI, Canoll P, Lignelli A, Hillman EM, Filippi CG, et al. Glioblastoma Induces Vascular Dysregulation in Nonenhancing Peritumoral Regions in Humans. AJR Am J Roentgenol. 2016:1–9. [PMC free article] [PubMed]
21. Aubert A, Costalat R, Duffau H, Benali H. Modeling of pathophysiological coupling between brain electrical activation, energy metabolism and hemodynamics: insights for the interpretation of intracerebral tumor imaging. Acta biotheoretica. 2002;50(4):281–95. [PubMed]
22. Aubert A, Costalat R. A model of the coupling between brain electrical activity, metabolism, and hemodynamics: application to the interpretation of functional neuroimaging. NeuroImage. 2002;17(3):1162–81. [PubMed]
23. Sakatani K, Murata Y, Fujiwara N, Hoshino T, Nakamura S, Kano T, et al. Comparison of blood-oxygen-level-dependent functional magnetic resonance imaging and near-infrared spectroscopy recording during functional brain activation in patients with stroke and brain tumors. Journal of biomedical optics. 2007;12(6):062110 doi: 10.1117/1.2823036 [PubMed]
24. Ko AL, Weaver KE, Hakimian S, Ojemann JG. Identifying functional networks using endogenous connectivity in gamma band electrocorticography. Brain Connect. 2013;3(5):491–502. PubMed Central PMCID: PMC3796331. doi: 10.1089/brain.2013.0157 [PMC free article] [PubMed]
25. Gunduz A, Brunner P, Daitch A, Leuthardt EC, Ritaccio AL, Pesaran B, et al. Neural correlates of visual-spatial attention in electrocorticographic signals in humans. Front Hum Neurosci. 2011;5:89 PubMed Central PMCID: PMC3202224. doi: 10.3389/fnhum.2011.00089 [PMC free article] [PubMed]
26. Schalk G, McFarland DJ, Hinterberger T, Birbaumer N, Wolpaw JR. BCI2000: a general-purpose brain-computer interface (BCI) system. IEEE Trans Biomed Eng. 2004;51(6):1034–43. Epub 2004/06/11. doi: 10.1109/TBME.2004.827072 [PubMed]
27. Welch PD. The Use of Fast Fourier Transform for the Estimation of Power Spectra: A Method Based on Time Averaging Over Short, Modified Periodograms. IEEE Transactions on Audio Electroacoustics. 1967;AU-15:70–3.
28. Miller KJ, Sorensen LB, Ojemann JG, den Nijs M. Power-law scaling in the brain surface electric potential. PLoS computational biology. 2009;5(12):e1000609 PubMed Central PMCID: PMC2787015. doi: 10.1371/journal.pcbi.1000609 [PMC free article] [PubMed]
29. Leopold DA, Murayama Y, Logothetis NK. Very slow activity fluctuations in monkey visual cortex: implications for functional brain imaging. Cereb Cortex. 2003;13(4):422–33. [PubMed]
30. Liu X, Yanagawa T, Leopold DA, Fujii N, Duyn JH. Robust Long-Range Coordination of Spontaneous Neural Activity in Waking, Sleep and Anesthesia. Cereb Cortex. 2015;25(9):2929–38. PubMed Central PMCID: PMCPMC4537437. doi: 10.1093/cercor/bhu089 [PMC free article] [PubMed]
31. Leopold DA, Maier A. Ongoing physiological processes in the cerebral cortex. Neuroimage. 2012;62(4):2190–200. PubMed Central PMCID: PMCPMC3288739. doi: 10.1016/j.neuroimage.2011.10.059 [PMC free article] [PubMed]
32. Birbaumer N, Elbert T, Canavan AG, Rockstroh B. Slow potentials of the cerebral cortex and behavior. Physiological reviews. 1990;70(1):1–41. [PubMed]
33. Mitzdorf U. Current source-density method and application in cat cerebral cortex: investigation of evoked potentials and EEG phenomena. Physiological reviews. 1985;65(1):37–100. [PubMed]
34. Vannini E, Olimpico F, Middei S, Ammassari-Teule M, de Graaf EL, McDonnell L, et al. Electrophysiology of glioma: a Rho GTPase-activating protein reduces tumor growth and spares neuron structure and function. Neuro Oncol. 2016. [PubMed]
35. Felleman DJ, Van Essen DC. Distributed hierarchical processing in the primate cerebral cortex. Cereb Cortex. 1991;1(1):1–47. Epub 1991/01/01. [PubMed]
36. Rosanova M, Casali A, Bellina V, Resta F, Mariotti M, Massimini M. Natural frequencies of human corticothalamic circuits. J Neurosci. 2009;29(24):7679–85. Epub 2009/06/19. doi: 10.1523/JNEUROSCI.0445-09.2009 [PubMed]
37. White NS, Alkire MT. Impaired thalamocortical connectivity in humans during general-anesthetic-induced unconsciousness. NeuroImage. 2003;19(2 Pt 1):402–11. Epub 2003/06/20. [PubMed]
38. Niedermeyer E, Lopes da Silva FH. Electroencephalography, basic principles, clinical applications, and related fields Baltimore: Urban & Schwarzenberg; 1982. x, 752 p. p.
39. Albert NB, Robertson EM, Miall RC. The resting human brain and motor learning. Curr Biol. 2009;19(12):1023–7. Epub 2009/05/12. PubMed Central PMCID: PMC2701987. doi: 10.1016/j.cub.2009.04.028 [PMC free article] [PubMed]
40. Eheman C, Henley SJ, Ballard-Barbash R, Jacobs EJ, Schymura MJ, Noone AM, et al. Annual Report to the Nation on the status of cancer, 1975–2008, featuring cancers associated with excess weight and lack of sufficient physical activity. Cancer. 2012;118(9):2338–66. Epub 2012/03/31. doi: 10.1002/cncr.27514 [PMC free article] [PubMed]
41. Davis FG, McCarthy BJ, Berger MS. Centralized databases available for describing primary brain tumor incidence, survival, and treatment: Central Brain Tumor Registry of the United States; Surveillance, Epidemiology, and End Results; and National Cancer Data Base. Neuro Oncol. 1999;1(3):205–11. Epub 2001/09/14. PubMed Central PMCID: PMC1920744. [PMC free article] [PubMed]
42. Larjavaara S, Mantyla R, Salminen T, Haapasalo H, Raitanen J, Jaaskelainen J, et al. Incidence of gliomas by anatomic location. Neuro Oncol. 2007;9(3):319–25. Epub 2007/05/25. PubMed Central PMCID: PMC1907421. doi: 10.1215/15228517-2007-016 [PMC free article] [PubMed]
43. Bello L, Acerbi F, Giussani C, Baratta P, Taccone P, Songa V, et al. Intraoperative language localization in multilingual patients with gliomas. Neurosurgery. 2006;59(1):115–25; discussion -25. [PubMed]
44. Leuthardt EC, Miller K, Anderson NR, Schalk G, Dowling J, Miller J, et al. Electrocorticographic frequency alteration mapping: a clinical technique for mapping the motor cortex. Neurosurgery. 2007;60(4 Suppl 2):260–70; discussion 70–1. [PubMed]
45. Roland JL, Hacker CD, Breshears JD, Gaona CM, Hogan RE, Burton H, et al. Brain mapping in a patient with congenital blindness—a case for multimodal approaches. Front Hum Neurosci. 2013;7:431 PubMed Central PMCID: PMCPMC3728570. doi: 10.3389/fnhum.2013.00431 [PMC free article] [PubMed]
46. Wang Y, Fifer MS, Flinker A, Korzeniewska A, Cervenka MC, Anderson WS, et al. Spatial-temporal functional mapping of language at the bedside with electrocorticography. Neurology. 2016;86(13):1181–9. PubMed Central PMCID: PMCPMC4818563. doi: 10.1212/WNL.0000000000002525 [PMC free article] [PubMed]
47. Lopes da Silva FH, Pfurtscheller G, editors. Event-Related Desynchronization. Handbook of Electroencephalography and Clinical Neurophysiology. Elsevier, Amsterdam: Elsevier; 1999.
48. Wu M, Wisneski K, Schalk G, Sharma M, Roland J, Breshears J, et al. Electrocorticographic frequency alteration mapping for extraoperative localization of speech cortex. Neurosurgery. 2010;66(2):E407–9. doi: 10.1227/01.NEU.0000345352.13696.6F [PubMed]
49. Brunner P, Ritaccio AL, Lynch TM, Emrich JF, Wilson JA, Williams JC, et al. A practical procedure for real-time functional mapping of eloquent cortex using electrocorticographic signals in humans. Epilepsy & behavior: E&B. 2009;15(3):278–86. PubMed Central PMCID: PMC2754703. [PMC free article] [PubMed]
50. Breshears JD, Gaona CM, Roland JL, Sharma M, Bundy DT, Shimony JS, et al. Mapping sensorimotor cortex with slow cortical potential resting-state networks while awake and under anesthesia. Neurosurgery. 2012;71(2):305–16; discussion 16. doi: 10.1227/NEU.0b013e318258e5d1 [PMC free article] [PubMed]
51. He B, Snyder A, Zempel J, Smyth M, Raichle M. Electrophysiological correlates of the brain's intrinsic large-scale functional architecture. Proceedings of the National Academy of Sciences. 2008;105(41):16039. [PubMed]
52. Larson-Prior LJ, Zempel JM, Nolan TS, Prior FW, Snyder AZ, Raichle ME. Cortical network functional connectivity in the descent to sleep. Proc Natl Acad Sci U S A. 2009;106(11):4489–94. Epub 2009/03/04. PubMed Central PMCID: PMC2657465. doi: 10.1073/pnas.0900924106 [PubMed]
53. Breshears JD, Roland JL, Sharma M, Gaona CM, Freudenburg ZV, Tempelhoff R, et al. Stable and dynamic cortical electrophysiology of induction and emergence with propofol anesthesia. Proc Natl Acad Sci U S A. 2010;107(49):21170–5. Epub 2010/11/17. PubMed Central PMCID: PMC3000270. doi: 10.1073/pnas.1011949107 [PubMed]
54. Vincent JL, Patel GH, Fox MD, Snyder AZ, Baker JT, Van Essen DC, et al. Intrinsic functional architecture in the anaesthetized monkey brain. Nature. 2007;447(7140):83–6. Epub 2007/05/04. doi: 10.1038/nature05758 [PubMed]
55. Jiang Z, Krainik A, David O, Salon C, Tropres I, Hoffmann D, et al. Impaired fMRI activation in patients with primary brain tumors. NeuroImage. 2010;52(2):538–48. doi: 10.1016/j.neuroimage.2010.04.194 [PubMed]
56. Holodny AI, Schulder M, Liu WC, Maldjian JA, Kalnin AJ. Decreased BOLD functional MR activation of the motor and sensory cortices adjacent to a glioblastoma multiforme: implications for image-guided neurosurgery. AJNR Am J Neuroradiol. 1999;20(4):609–12. [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science