Glioblastoma remodelling of human neural circuits decreases survival.

Gliomas synaptically integrate into neural circuits1,2. Previous research has demonstrated bidirectional interactions between neurons and glioma cells, with neuronal activity driving glioma growth1-4 and gliomas increasing neuronal excitability2,5-8. Here we sought to determine how glioma-induced neuronal changes influence neural circuits underlying cognition and whether these interactions influence patient survival. Using intracranial brain recordings during lexical retrieval language tasks in awake humans together with site-specific tumour tissue biopsies and cell biology experiments, we find that gliomas remodel functional neural circuitry such that task-relevant neural responses activate tumour-infiltrated cortex well beyond the cortical regions that are normally recruited in the healthy brain. Site-directed biopsies from regions within the tumour that exhibit high functional connectivity between the tumour and the rest of the brain are enriched for a glioblastoma subpopulation that exhibits a distinct synaptogenic and neuronotrophic phenotype. Tumour cells from functionally connected regions secrete the synaptogenic factor thrombospondin-1, which contributes to the differential neuron-glioma interactions observed in functionally connected tumour regions compared with tumour regions with less functional connectivity. Pharmacological inhibition of thrombospondin-1 using the FDA-approved drug gabapentin decreases glioblastoma proliferation. The degree of functional connectivity between glioblastoma and the normal brain negatively affects both patient survival and performance in language tasks. These data demonstrate that high-grade gliomas functionally remodel neural circuits in the human brain, which both promotes tumour progression and impairs cognition.

Phase Alignment of Low-Frequency Neural Activity to the Amplitude Envelope of Speech Reflects Evoked Responses to Acoustic Edges, Not Oscillatory Entrainment.

The amplitude envelope of speech is crucial for accurate comprehension. Considered a key stage in speech processing, the phase of neural activity in the theta-delta bands (1-10 Hz) tracks the phase of the speech amplitude envelope during listening. However, the mechanisms underlying this envelope representation have been heavily debated. A dominant model posits that envelope tracking reflects entrainment of endogenous low-frequency oscillations to the speech envelope. Alternatively, envelope tracking reflects a series of evoked responses to acoustic landmarks within the envelope. It has proven challenging to distinguish these two mechanisms. To address this, we recorded MEG while participants (n = 12, 6 female) listened to natural speech, and compared the neural phase patterns to the predictions of two computational models: an oscillatory entrainment model and a model of evoked responses to peaks in the rate of envelope change. Critically, we also presented speech at slowed rates, where the spectro-temporal predictions of the two models diverge. Our analyses revealed transient theta phase-locking in regular speech, as predicted by both models. However, for slow speech, we found transient theta and delta phase-locking, a pattern that was fully compatible with the evoked response model but could not be explained by the oscillatory entrainment model. Furthermore, encoding of acoustic edge magnitudes was invariant to contextual speech rate, demonstrating speech rate normalization of acoustic edge representations. Together, our results suggest that neural phase-locking to the speech envelope is more likely to reflect discrete representation of transient information rather than oscillatory entrainment.SIGNIFICANCE STATEMENT This study probes a highly debated topic in speech perception: the neural mechanisms underlying the cortical representation of the temporal envelope of speech. It is well established that the slow intensity profile of the speech signal, its envelope, elicits a robust brain response that "tracks" these envelope fluctuations. The oscillatory entrainment model posits that envelope tracking reflects phase alignment of endogenous neural oscillations. Here the authors provide evidence for a distinct mechanism. They show that neural speech envelope tracking arises from transient evoked neural responses to rapid increases in the speech envelope. Explicit computational modeling provides direct and compelling evidence that evoked responses are the primary mechanism underlying cortical speech envelope representations, with no evidence for oscillatory entrainment.

Cortical Synchrony and Information Flow during Transition from Wakefulness to Light Non-Rapid Eye Movement Sleep.

Sleep is a highly stereotyped phenomenon, requiring robust spatiotemporal coordination of neural activity. Understanding how the brain coordinates neural activity with sleep onset can provide insights into the physiological functions subserved by sleep and the pathologic phenomena associated with sleep onset. We quantified whole-brain network changes in synchrony and information flow during the transition from wakefulness to light non-rapid eye movement (NREM) sleep, using MEG imaging in a convenient sample of 14 healthy human participants (11 female; mean 63.4 years [SD 11.8 years]). We furthermore performed computational modeling to infer excitatory and inhibitory properties of local neural activity. The transition from wakefulness to light NREM was identified to be encoded in spatially and temporally specific patterns of long-range synchrony. Within the delta band, there was a global increase in connectivity from wakefulness to light NREM, which was highest in frontoparietal regions. Within the theta band, there was an increase in connectivity in fronto-parieto-occipital regions and a decrease in temporal regions from wakefulness to Stage 1 sleep. Patterns of information flow revealed that mesial frontal regions receive hierarchically organized inputs from broad cortical regions upon sleep onset, including direct inflow from occipital regions and indirect inflow via parieto-temporal regions within the delta frequency band. Finally, biophysical neural mass modeling demonstrated changes in the anterior-to-posterior distribution of cortical excitation-to-inhibition with increased excitation-to-inhibition model parameters in anterior regions in light NREM compared with wakefulness. Together, these findings uncover whole-brain corticocortical structure and the orchestration of local and long-range, frequency-specific cortical interactions in the sleep-wake transition.SIGNIFICANCE STATEMENT Our work uncovers spatiotemporal cortical structure of neural synchrony and information flow upon the transition from wakefulness to light non-rapid eye movement sleep. Mesial frontal regions were identified to receive hierarchically organized inputs from broad cortical regions, including both direct inputs from occipital regions and indirect inputs via the parieto-temporal regions within the delta frequency range. Biophysical neural mass modeling revealed a spatially heterogeneous, anterior-posterior distribution of cortical excitation-to-inhibition. Our findings shed light on the orchestration of local and long-range cortical neural structure that is fundamental to sleep onset, and support an emerging view of cortically driven regulation of sleep homeostasis.

Neuronal synchrony abnormalities associated with subclinical epileptiform activity in early-onset Alzheimer's disease.

Since the first demonstrations of network hyperexcitability in scientific models of Alzheimer's disease, a growing body of clinical studies have identified subclinical epileptiform activity and associated cognitive decline in patients with Alzheimer's disease. An obvious problem presented in these studies is lack of sensitive measures to detect and quantify network hyperexcitability in human subjects. In this study we examined whether altered neuronal synchrony can be a surrogate marker to quantify network hyperexcitability in patients with Alzheimer's disease. Using magnetoencephalography (MEG) at rest, we studied 30 Alzheimer's disease patients without subclinical epileptiform activity, 20 Alzheimer's disease patients with subclinical epileptiform activity and 35 age-matched controls. Presence of subclinical epileptiform activity was assessed in patients with Alzheimer's disease by long-term video-EEG and a 1-h resting MEG with simultaneous EEG. Using the resting-state source-space reconstructed MEG signal, in patients and controls we computed the global imaginary coherence in alpha (8-12 Hz) and delta-theta (2-8 Hz) oscillatory frequencies. We found that Alzheimer's disease patients with subclinical epileptiform activity have greater reductions in alpha imaginary coherence and greater enhancements in delta-theta imaginary coherence than Alzheimer's disease patients without subclinical epileptiform activity, and that these changes can distinguish between Alzheimer's disease patients with subclinical epileptiform activity and Alzheimer's disease patients without subclinical epileptiform activity with high accuracy. Finally, a principal component regression analysis showed that the variance of frequency-specific neuronal synchrony predicts longitudinal changes in Mini-Mental State Examination in patients and controls. Our results demonstrate that quantitative neurophysiological measures are sensitive biomarkers of network hyperexcitability and can be used to improve diagnosis and to select appropriate patients for the right therapy in the next-generation clinical trials. The current results provide an integrative framework for investigating network hyperexcitability and network dysfunction together with cognitive and clinical correlates in patients with Alzheimer's disease.

Clinical Validation of the Champagne Algorithm for Evoked Response Source Localization in Magnetoencephalography.

Magnetoencephalography (MEG) is a robust method for non-invasive functional brain mapping of sensory cortices due to its exceptional spatial and temporal resolution. The clinical standard for MEG source localization of functional landmarks from sensory evoked responses is the equivalent current dipole (ECD) localization algorithm, known to be sensitive to initialization, noise, and manual choice of the number of dipoles. Recently many automated and robust algorithms have been developed, including the Champagne algorithm, an empirical Bayesian algorithm, with powerful abilities for MEG source reconstruction and time course estimation (Wipf et al. 2010; Owen et al. 2012). Here, we evaluate automated Champagne performance in a clinical population of tumor patients where there was minimal failure in localizing sensory evoked responses using the clinical standard, ECD localization algorithm. MEG data of auditory evoked potentials and somatosensory evoked potentials from 21 brain tumor patients were analyzed using Champagne, and these results were compared with equivalent current dipole (ECD) fit. Across both somatosensory and auditory evoked field localization, we found there was a strong agreement between Champagne and ECD localizations in all cases. Given resolution of 8mm voxel size, peak source localizations from Champagne were below 10mm of ECD peak source localization. The Champagne algorithm provides a robust and automated alternative to manual ECD fits for clinical localization of sensory evoked potentials and can contribute to improved clinical MEG data processing workflows.

Interhemispheric Auditory Cortical Synchronization in Asymmetric Hearing Loss.

OBJECTIVES: Auditory cortical activation of the two hemispheres to monaurally presented tonal stimuli has been shown to be asynchronous in normal hearing (NH) but synchronous in the extreme case of adult-onset asymmetric hearing loss (AHL) with single-sided deafness. We addressed the wide knowledge gap between these two anchoring states of interhemispheric temporal organization. The objectives of this study were as follows: (1) to map the trajectory of interhemispheric temporal reorganization from asynchrony to synchrony using magnitude of interaural threshold difference as the independent variable in a cross-sectional study and (2) to evaluate reversibility of interhemispheric synchrony in association with hearing in noise performance by amplifying the aidable poorer ear in a repeated measures, longitudinal study. DESIGN: The cross-sectional and longitudinal cohorts were comprised of 49 subjects (AHL; N = 21; 11 male, 10 female; mean age = 48 years) and NH (N = 28; 16 male, 12 female; mean age = 45 years). The maximum interaural threshold difference of the two cohorts spanned from 0 to 65 dB. Magnetoencephalography analyses focused on latency of the M100 peak response from auditory cortex in both hemispheres between 50 msec and 150 msec following monaural tonal stimulation at the frequency (0.5, 1, 2, 3, or 4 kHz) corresponding to the maximum and minimum interaural threshold difference for better and poorer ears separately. The longitudinal AHL cohort was drawn from three subjects in the cross-sectional AHL cohort (all male; ages 49 to 60 years; varied AHL etiologies; no amplification for at least 2 years). All longitudinal study subjects were treated by monaural amplification of the poorer ear and underwent repeated measures examination of the M100 response latency and quick speech in noise hearing in noise performance at baseline, and postamplification months 3, 6, and 12. RESULTS: The M100 response peak latency values in the ipsilateral hemisphere lagged those in the contralateral hemisphere for all stimulation conditions. The mean (SD) interhemispheric latency difference values (ipsilateral less contralateral) to better ear stimulation for three categories of maximum interaural threshold difference were as follows: NH (≤ 10 dB)-8.6 (3.0) msec; AHL (15 to 40 dB)-3.0 (1.2) msec; AHL (≥ 45 dB)-1.4 (1.3) msec. In turn, the magnitude of difference values were used to define interhemispheric temporal organization states of asynchrony, mixed asynchrony and synchrony, and synchrony, respectively. Amplification of the poorer ear in longitudinal subjects drove interhemispheric organization change from baseline synchrony to postamplification asynchrony and hearing in noise performance improvement in those with baseline impairment over a 12-month period. CONCLUSIONS: Interhemispheric temporal organization in AHL was anchored between states of asynchrony in NH and synchrony in single-sided deafness. For asymmetry magnitudes between 15 and 40 dB, the intermediate mixed state of asynchrony and synchrony was continuous and reversible. Amplification of the poorer ear in AHL improved hearing in noise performance and restored normal temporal organization of auditory cortices in the two hemispheres. The return to normal interhemispheric asynchrony from baseline synchrony and improvement in hearing following monoaural amplification of the poorer ear evolved progressively over a 12-month period.

Sensorimotor Cortical Oscillations during Movement Preparation in 16p11.2 Deletion Carriers.

Sensorimotor deficits are prevalent in many neurodevelopmental disorders like autism, including one of its common genetic etiologies, a 600 kb reciprocal deletion/duplication at 16p11.2. We have previously shown that copy number variations of 16p11.2 impact regional brain volume, white matter integrity, and early sensory responses in auditory cortex. Here, we test the hypothesis that abnormal cortical neurophysiology is present when genes in the 16p11.2 region are haploinsufficient, and in humans that this in turn may account for behavioral deficits specific to deletion carriers. We examine sensorimotor cortical network activity in males and females with 16p11.2 deletions compared with both typically developing individuals, and those with duplications of 16p11.2, using magnetoencephalographic imaging during preparation of overt speech or hand movements in tasks designed to be easy for all participants. In deletion carriers, modulation of beta oscillations (12-30 Hz) were increased during both movement types over effector-specific regions of motor cortices compared with typically developing individuals or duplication carriers, with no task-related performance differences between cohorts, even when corrected for their own cognitive and sensorimotor deficits. Reduced left hemispheric language specialization was observed in deletion carriers but not in duplication carriers. Neural activity over sensorimotor cortices in deletion carriers was linearly related to clinical measures of speech and motor impairment. These findings link insufficient copy number repeats at 16p11.2 to excessive neural activity (e.g., increased beta oscillations) in motor cortical networks for speech and hand motor control. These results have significant implications for understanding the neural basis of autism and related neurodevelopmental disorders.SIGNIFICANCE STATEMENT The recurrent ∼600 kb deletion at 16p11.2 (BP4-BP5) is one of the most common genetic etiologies of ASD and, more generally, of neurodevelopmental disorders. Here, we use high-resolution magnetoencephalographic imaging (MEG-I) to define with millisecond precision the underlying neurophysiological signature of motor impairments for individuals with 16p11.2 deletions. We identify significant increases in beta (12-30 Hz) suppression in sensorimotor cortices related to performance during speech and hand movement tasks. These findings not only provide a neurophysiological phenotype for the clinical presentation of this genetic deletion, but also guide our understanding of how genetic variation encodes for neural oscillatory dynamics.

Global resting-state functional connectivity of neural oscillations in tinnitus with and without hearing loss.

This study examined global resting-state functional connectivity of neural oscillations in individuals with chronic tinnitus and normal and impaired hearing. We tested the hypothesis that distinct neural oscillatory networks are engaged in tinnitus with and without hearing loss. In both tinnitus groups, with and without hearing loss, we identified multiple frequency band-dependent regions of increased and decreased global functional connectivity. We also found that the auditory domain of tinnitus severity, assayed by the Tinnitus Functional Index, was associated with global functional connectivity in both auditory and nonauditory regions. These findings provide candidate biomarkers to target and monitor treatments for tinnitus with and without hearing loss.

MEG imaging of recurrent gliomas reveals functional plasticity of hemispheric language specialization.

In patients with gliomas, changes in hemispheric specialization for language determined by magnetoencephalography (MEG) were analyzed to elucidate the impact of treatment and tumor recurrence on language networks. Demonstration of reorganization of language networks in these patients has significant implications on the prevention of postoperative functional loss and recovery. Whole-brain activity during an auditory verb generation task was estimated from MEG recordings in a group of 73 patients with recurrent gliomas. Hemisphere of language dominance was estimated using the language laterality index (LI), a measure derived from the task. The initial scan was performed prior to resection; patients subsequently underwent surgery and adjuvant treatment. A second scan was performed upon recurrence prior to repeat resection. The relationship between the shift in LI between scans and demographics, anatomic location, pathology, and adjuvant treatment was analyzed. Laterality shifts were observed between scans; the median percent change was 29.1% across all patients. Laterality shift magnitude and relative direction were associated with the initial position of language dominance; patients with increased lateralization experienced greater shifts than those presenting more bilateral representation. A change in LI from left or right to bilateral (or vice versa) occurred in 23.3% of patients; complete switch occurred in 5.5% of patients. Patients with tumors within the language-dominant hemisphere experienced significantly greater shifts than those with contralateral tumors. The majority of patients with glioma experience shifts in language network organization over time which correlate with the relative position of language lateralization and tumor location.

The Contribution of the Corpus Callosum to Language Lateralization.

The development of hemispheric lateralization for language is poorly understood. In one hypothesis, early asymmetric gene expression assigns language to the left hemisphere. In an alternate view, language is represented a priori in both hemispheres and lateralization emerges via cross-hemispheric communication through the corpus callosum. To address this second hypothesis, we capitalized on the high temporal and spatial resolution of magnetoencephalographic imaging to measure cortical activity during language processing, speech preparation, and speech execution in 25 participants with agenesis of the corpus callosum (AgCC) and 21 matched neurotypical individuals. In contrast to strongly lateralized left hemisphere activations for language in neurotypical controls, participants with complete or partial AgCC exhibited bilateral hemispheric activations in both auditory or visually driven language tasks, with complete AgCC participants showing significantly more right hemisphere activations than controls or than individuals with partial AgCC. In AgCC individuals, language laterality positively correlated with verbal IQ. These findings suggest that the corpus callosum helps to drive language lateralization. SIGNIFICANCE STATEMENT: The role that corpus callosum development has on the hemispheric specialization of language is poorly understood. Here, we used magnetoencephalographic imaging during linguistic tests (verb generation, picture naming) to test for hemispheric dominance in patients with agenesis of the corpus callosum (AgCC) and found reduced laterality (i.e., greater likelihood of bilaterality or right hemisphere dominance) in this cohort compared with controls, especially in patients with complete agenesis. Laterality was positively correlated with behavioral measures of verbal intelligence. These findings provide support for the hypothesis that the callosum aids in functional specialization throughout neural development and that the loss of this mechanism correlates with impairments in verbal performance.

Global and regional functional connectivity maps of neural oscillations in focal epilepsy.

Intractable focal epilepsy is a devastating disorder with profound effects on cognition and quality of life. Epilepsy surgery can lead to seizure freedom in patients with focal epilepsy; however, sometimes it fails due to an incomplete delineation of the epileptogenic zone. Brain networks in epilepsy can be studied with resting-state functional connectivity analysis, yet previous investigations using functional magnetic resonance imaging or electrocorticography have produced inconsistent results. Magnetoencephalography allows non-invasive whole-brain recordings, and can be used to study both long-range network disturbances in focal epilepsy and regional connectivity at the epileptogenic zone. In magnetoencephalography recordings from presurgical epilepsy patients, we examined: (i) global functional connectivity maps in patients versus controls; and (ii) regional functional connectivity maps at the region of resection, compared to the homotopic non-epileptogenic region in the contralateral hemisphere. Sixty-one patients were studied, including 30 with mesial temporal lobe epilepsy and 31 with focal neocortical epilepsy. Compared with a group of 31 controls, patients with epilepsy had decreased resting-state functional connectivity in widespread regions, including perisylvian, posterior temporo-parietal, and orbitofrontal cortices (P < 0.01, t-test). Decreased mean global connectivity was related to longer duration of epilepsy and higher frequency of consciousness-impairing seizures (P < 0.01, linear regression). Furthermore, patients with increased regional connectivity within the resection site (n = 24) were more likely to achieve seizure postoperative seizure freedom (87.5% with Engel I outcome) than those with neutral (n = 15, 64.3% seizure free) or decreased (n = 23, 47.8% seizure free) regional connectivity (P < 0.02, chi-square). Widespread global decreases in functional connectivity are observed in patients with focal epilepsy, and may reflect deleterious long-term effects of recurrent seizures. Furthermore, enhanced regional functional connectivity at the area of resection may help predict seizure outcome and aid surgical planning.

Global MEG Resting State Functional Connectivity in Children with Autism and Sensory Processing Dysfunction.

Sensory processing dysfunction not only affects most individuals with autism spectrum disorder (ASD), but at least 5% of children without ASD also experience dysfunctional sensory processing. Our understanding of the relationship between sensory dysfunction and resting state brain activity is still emerging. This study compared long-range resting state functional connectivity of neural oscillatory behavior in children aged 8-12 years with autism spectrum disorder (ASD; N=18), those with sensory processing dysfunction (SPD; N=18) who do not meet ASD criteria, and typically developing control participants (TDC; N=24) using magnetoencephalography (MEG). Functional connectivity analyses were performed in the alpha and beta frequency bands, which are known to be implicated in sensory information processing. Group differences in functional connectivity and associations between sensory abilities and functional connectivity were examined. Distinct patterns of functional connectivity differences between ASD and SPD groups were found only in the beta band, but not in the alpha band. In both alpha and beta bands, ASD and SPD cohorts differed from the TDC cohort. Somatosensory cortical beta-band functional connectivity was associated with tactile processing abilities, while higher-order auditory cortical alpha-band functional connectivity was associated with auditory processing abilities. These findings demonstrate distinct long-range neural synchrony alterations in SPD and ASD that are associated with sensory processing abilities. Neural synchrony measures could serve as potential sensitive biomarkers for ASD and SPD.

Abnormal gamma phase-amplitude coupling in the parahippocampal cortex is associated with network hyperexcitability in Alzheimer's disease.

While animal models of Alzheimer's disease (AD) have shown altered gamma oscillations (∼40 Hz) in local neural circuits, the low signal-to-noise ratio of gamma in the resting human brain precludes its quantification via conventional spectral estimates. Phase-amplitude coupling (PAC) indicating the dynamic integration between the gamma amplitude and the phase of low-frequency (4-12 Hz) oscillations is a useful alternative to capture local gamma activity. In addition, PAC is also an index of neuronal excitability as the phase of low-frequency oscillations that modulate gamma amplitude, effectively regulates the excitability of local neuronal firing. In this study, we sought to examine the local neuronal activity and excitability using gamma PAC, within brain regions vulnerable to early AD pathophysiology-entorhinal cortex and parahippocampus, in a clinical population of patients with AD and age-matched controls. Our clinical cohorts consisted of a well-characterized cohort of AD patients (n = 50; age, 60 ± 8 years) with positive AD biomarkers, and age-matched, cognitively unimpaired controls (n = 35; age, 63 ± 5.8 years). We identified the presence or the absence of epileptiform activity in AD patients (AD patients with epileptiform activity, AD-EPI+, n = 20; AD patients without epileptiform activity, AD-EPI-, n = 30) using long-term electroencephalography (LTM-EEG) and 1-hour long magnetoencephalography (MEG) with simultaneous EEG. Using the source reconstructed MEG data, we computed gamma PAC as the coupling between amplitude of the gamma frequency (30-40 Hz) with phase of the theta (4-8 Hz) and alpha (8-12 Hz) frequency oscillations, within entorhinal and parahippocampal cortices. We found that patients with AD have reduced gamma PAC in the left parahippocampal cortex, compared to age-matched controls. Furthermore, AD-EPI+ patients showed greater reductions in gamma PAC than AD-EPI- in bilateral parahippocampal cortices. In contrast, entorhinal cortices did not show gamma PAC abnormalities in patients with AD. Our findings demonstrate the spatial patterns of altered gamma oscillations indicating possible region-specific manifestations of network hyperexcitability within medial temporal lobe regions vulnerable to AD pathophysiology. Greater deficits in AD-EPI+ suggests that reduced gamma PAC is a sensitive index of network hyperexcitability in AD patients. Collectively, the current results emphasize the importance of investigating the role of neural circuit hyperexcitability in early AD pathophysiology and explore its potential as a modifiable contributor to AD pathobiology.

Language mapping with navigated repetitive TMS: proof of technique and validation.

OBJECTIVE: Lesion-based mapping of speech pathways has been possible only during invasive neurosurgical procedures using direct cortical stimulation (DCS). However, navigated transcranial magnetic stimulation (nTMS) may allow for lesion-based interrogation of language pathways noninvasively. Although not lesion-based, magnetoencephalographic imaging (MEGI) is another noninvasive modality for language mapping. In this study, we compare the accuracy of nTMS and MEGI with DCS. METHODS: Subjects with lesions around cortical language areas underwent preoperative nTMS and MEGI for language mapping. nTMS maps were generated using a repetitive TMS protocol to deliver trains of stimulations during a picture naming task. MEGI activation maps were derived from adaptive spatial filtering of beta-band power decreases prior to overt speech during picture naming and verb generation tasks. The subjects subsequently underwent awake language mapping via intraoperative DCS. The language maps obtained from each of the 3 modalities were recorded and compared. RESULTS: nTMS and MEGI were performed on 12 subjects. nTMS yielded 21 positive language disruption sites (11 speech arrest, 5 anomia, and 5 other) while DCS yielded 10 positive sites (2 speech arrest, 5 anomia, and 3 other). MEGI isolated 32 sites of peak activation with language tasks. Positive language sites were most commonly found in the pars opercularis for all three modalities. In 9 instances the positive DCS site corresponded to a positive nTMS site, while in 1 instance it did not. In 4 instances, a positive nTMS site corresponded to a negative DCS site, while 169 instances of negative nTMS and DCS were recorded. The sensitivity of nTMS was therefore 90%, specificity was 98%, the positive predictive value was 69% and the negative predictive value was 99% as compared with intraoperative DCS. MEGI language sites for verb generation and object naming correlated with nTMS sites in 5 subjects, and with DCS sites in 2 subjects. CONCLUSION: Maps of language function generated with nTMS correlate well with those generated by DCS. Negative nTMS mapping also correlates with negative DCS mapping. In our study, MEGI lacks the same level of correlation with intraoperative mapping; nevertheless it provides useful adjunct information in some cases. nTMS may offer a lesion-based method for noninvasively interrogating language pathways and be valuable in managing patients with peri-eloquent lesions.

Dynamics of hemispheric dominance for language assessed by magnetoencephalographic imaging.

OBJECTIVE: The goal of the current study was to examine the dynamics of language lateralization using magnetoencephalographic (MEG) imaging, to determine the sensitivity and specificity of MEG imaging, and to determine whether MEG imaging can become a viable alternative to the intracarotid amobarbital procedure (IAP), the current gold standard for preoperative language lateralization in neurosurgical candidates. METHODS: MEG was recorded during an auditory verb generation task and imaging analysis of oscillatory activity was initially performed in 21 subjects with epilepsy, brain tumor, or arteriovenous malformation who had undergone IAP and MEG. Time windows and brain regions of interest that best discriminated between IAP-determined left or right dominance for language were identified. Parameters derived in the retrospective analysis were applied to a prospective cohort of 14 patients and healthy controls. RESULTS: Power decreases in the beta frequency band were consistently observed following auditory stimulation in inferior frontal, superior temporal, and parietal cortices; similar power decreases were also seen in inferior frontal cortex prior to and during overt verb generation. Language lateralization was clearly observed to be a dynamic process that is bilateral for several hundred milliseconds during periods of auditory perception and overt speech production. Correlation with the IAP was seen in 13 of 14 (93%) prospective patients, with the test demonstrating a sensitivity of 100% and specificity of 92%. INTERPRETATION: Our results demonstrate excellent correlation between MEG imaging findings and the IAP for language lateralization, and provide new insights into the spatiotemporal dynamics of cortical speech processing.

Resting functional connectivity in patients with brain tumors in eloquent areas.

OBJECTIVE: Resection of brain tumors adjacent to eloquent areas represents a challenge in neurosurgery. If maximal resection is desired without inducing postoperative neurological deficits, a detailed knowledge of the functional topography in and around the tumor is crucial. The aim of the present work is to evaluate the value of preoperative magnetoencephalography (MEG) imaging of functional connectivity to predict the results of intraoperative electrical stimulation (IES) mapping, the clinical gold standard for neurosurgical localization of functional areas. METHODS: Resting-state whole-cortex MEG recordings were obtained from 57 consecutive subjects with focal brain tumors near or within motor, sensory, or language areas. Neural activity was estimated using adaptive spatial filtering algorithms, and the mean imaginary coherence between the rest of the brain and voxels in and around brain tumors were compared to the mean imaginary coherence between the rest of the brain and contralesional voxels as an index of functional connectivity. IES mapping was performed in all subjects. The cortical connectivity pattern near the tumor was compared to the IES results. RESULTS: Maps with decreased resting-state functional connectivity in the entire tumor area had a negative predictive value of 100% for absence of eloquent cortex during IES. Maps showing increased resting-state functional connectivity within the tumor area had a positive predictive value of 64% for finding language, motor, or sensory cortical sites during IES mapping. INTERPRETATION: Preoperative resting state MEG connectivity analysis is a useful noninvasive tool to evaluate the functionality of the tissue surrounding tumors within eloquent areas, and could potentially contribute to surgical planning and patient counseling.

Magnetoencephalography Imaging Reveals Abnormal Information Flow in Temporal Lobe Epilepsy.

Background/Introduction: Widespread network disruption has been hypothesized to be an important predictor of outcomes in patients with refractory temporal lobe epilepsy (TLE). Most studies examining functional network disruption in epilepsy have largely focused on the symmetric bidirectional metrics of the strength of network connections. However, a more complete description of network dysfunction impacts in epilepsy requires an investigation of the potentially more sensitive directional metrics of information flow. Methods: This study describes a whole-brain magnetoencephalography-imaging approach to examine resting-state directional information flow networks, quantified by phase-transfer entropy (PTE), in patients with TLE compared with healthy controls (HCs). Associations between PTE and clinical characteristics of epilepsy syndrome are also investigated. Results: Deficits of information flow were specific to alpha-band frequencies. In alpha band, while HCs exhibit a clear posterior-to-anterior directionality of information flow, in patients with TLE, this pattern of regional information outflow and inflow was significantly altered in the frontal and occipital regions. The changes in information flow within the alpha band in selected brain regions were correlated with interictal spike frequency and duration of epilepsy. Conclusions: Impaired information flow is an important dimension of network dysfunction associated with the pathophysiological mechanisms of TLE.

Network connectivity predicts effectiveness of responsive neurostimulation in focal epilepsy.

Responsive neurostimulation is a promising treatment for drug-resistant focal epilepsy; however, clinical outcomes are highly variable across individuals. The therapeutic mechanism of responsive neurostimulation likely involves modulatory effects on brain networks; however, with no known biomarkers that predict clinical response, patient selection remains empiric. This study aimed to determine whether functional brain connectivity measured non-invasively prior to device implantation predicts clinical response to responsive neurostimulation therapy. Resting-state magnetoencephalography was obtained in 31 participants with subsequent responsive neurostimulation device implantation between 15 August 2014 and 1 October 2020. Functional connectivity was computed across multiple spatial scales (global, hemispheric, and lobar) using pre-implantation magnetoencephalography and normalized to maps of healthy controls. Normalized functional connectivity was investigated as a predictor of clinical response, defined as percent change in self-reported seizure frequency in the most recent year of clinic visits relative to pre-responsive neurostimulation baseline. Area under the receiver operating characteristic curve quantified the performance of functional connectivity in predicting responders (≥50% reduction in seizure frequency) and non-responders (<50%). Leave-one-out cross-validation was furthermore performed to characterize model performance. The relationship between seizure frequency reduction and frequency-specific functional connectivity was further assessed as a continuous measure. Across participants, stimulation was enabled for a median duration of 52.2 (interquartile range, 27.0-62.3) months. Demographics, seizure characteristics, and responsive neurostimulation lead configurations were matched across 22 responders and 9 non-responders. Global functional connectivity in the alpha and beta bands were lower in non-responders as compared with responders (alpha, pfdr < 0.001; beta, pfdr < 0.001). The classification of responsive neurostimulation outcome was improved by combining feature inputs; the best model incorporated four features (i.e. mean and dispersion of alpha and beta bands) and yielded an area under the receiver operating characteristic curve of 0.970 (0.919-1.00). The leave-one-out cross-validation analysis of this four-feature model yielded a sensitivity of 86.3%, specificity of 77.8%, positive predictive value of 90.5%, and negative predictive value of 70%. Global functional connectivity in alpha band correlated with seizure frequency reduction (alpha, P = 0.010). Global functional connectivity predicted responder status more strongly, as compared with hemispheric predictors. Lobar functional connectivity was not a predictor. These findings suggest that non-invasive functional connectivity may be a candidate personalized biomarker that has the potential to predict responsive neurostimulation effectiveness and to identify patients most likely to benefit from responsive neurostimulation therapy. Follow-up large-cohort, prospective studies are required to validate this biomarker. These findings furthermore support an emerging view that the therapeutic mechanism of responsive neurostimulation involves network-level effects in the brain.

Clinical Validation of the Champagne Algorithm for Epilepsy Spike Localization.

Magnetoencephalography (MEG) is increasingly used for presurgical planning in people with medically refractory focal epilepsy. Localization of interictal epileptiform activity, a surrogate for the seizure onset zone whose removal may prevent seizures, is challenging and depends on the use of multiple complementary techniques. Accurate and reliable localization of epileptiform activity from spontaneous MEG data has been an elusive goal. One approach toward this goal is to use a novel Bayesian inference algorithm-the Champagne algorithm with noise learning-which has shown tremendous success in source reconstruction, especially for focal brain sources. In this study, we localized sources of manually identified MEG spikes using the Champagne algorithm in a cohort of 16 patients with medically refractory epilepsy collected in two consecutive series. To evaluate the reliability of this approach, we compared the performance to equivalent current dipole (ECD) modeling, a conventional source localization technique that is commonly used in clinical practice. Results suggest that Champagne may be a robust, automated, alternative to manual parametric dipole fitting methods for localization of interictal MEG spikes, in addition to its previously described clinical and research applications.

NUTMEG: Open Source Software for M/EEG Source Reconstruction.

Neurodynamic Utility Toolbox for Magnetoencephalo- and Electroencephalography (NUTMEG) is an open-source MATLAB-based toolbox for the analysis and reconstruction of magnetoencephalography/electroencephalography data in source space. NUTMEG includes a variety of options for the user in data import, preprocessing, source reconstruction, and functional connectivity. A group analysis toolbox allows the user to run a variety of inferential statistics on their data in an easy-to-use GUI-driven format. Importantly, NUTMEG features an interactive five-dimensional data visualization platform. A key feature of NUTMEG is the availability of a large menu of interference cancelation and source reconstruction algorithms. Each NUTMEG operation acts as a stand-alone MATLAB function, allowing the package to be easily adaptable and scripted for the more advanced user for interoperability with other software toolboxes. Therefore, NUTMEG enables a wide range of users access to a complete "sensor-to- source-statistics" analysis pipeline.