Background and Hypothesis We used the uniquely high combined spatial and temporal resolution of magnetoencephalography to characterize working memory (WM)-related modulation of beta band activity in neuroleptic-free patients with schizophrenia in comparison to a large sample of performance-matched healthy controls. We also tested for effects of antipsychotic medication on identified differences in these same patients. Study Design Inpatients with schizophrenia (n = 21) or psychotic disorder not otherwise specified (n = 4) completed N-back and control tasks during magnetoencephalography while on placebo and during antipsychotic medication treatment, in a blinded, randomized, counterbalanced manner. Healthy, performance-matched controls (N = 100) completed the same tasks. WM-related neural activation was estimated as beta band (14-30 Hz) desynchronization throughout the brain in successive 400 ms time windows. Voxel-wise statistical comparisons were performed between controls and patients while off-medication at each time window. Significant clusters resulting from this between-groups analysis were then used as regions-of-interest, the activations of which were compared between on- and off-medication conditions in patients. Study Results Controls showed beta-band desynchronization (activation) of a fronto-parietal network immediately preceding correct button press responses-the time associated with WM updating and task execution. Altered activation in medication-free patients occurred largely during this time, in prefrontal, parietal, and visual cortices. Medication altered patients' neural responses such that the activation time courses in these regions-of-interest more closely resembled those of controls. Conclusions These findings demonstrate that WM-related beta band alterations in schizophrenia are time-specific and associated with neural systems targeted by antipsychotic medications. Future studies may investigate this association by examining its potential neurochemical basis.

A single-nucleotide polymorphism in the promoter region of the Matrix Metalloproteinase-9 (MMP9) gene, rs3918242, has been shown to affect MMP9 expression in macrophages and was associated with schizophrenia by two independent groups. However, rs3918242's effects on MMP9 expression were not replicable in cell lines or brain tissue. Additionally, publically available data indicate that rs3918242 genotype is related not to MMP9 expression, but rather to expression of SLC12A5, a nearby gene coding for a K+/Cl- cotransporter, whose expression has also been related to schizophrenia. Here, we studied brain structure and function in healthy participants stratified by rs3918242 genotype using structural MRI (N = 298), functional MRI during an N-back working memory task (N = 554), and magnetoencephalography (MEG) during the same task (N = 190). We found rs3918242 was associated with gray matter volume (GMV) in the insula and dorsolateral prefrontal cortex bilaterally, closely replicated in discovery and replication samples; and with inferior parietal lobule (IPL) GMV when the samples were meta-analytically combined. Additionally, using both fMRI and MEG, rs3918242 was associated with right IPL working memory-related activation, replicated in two cohorts and across imaging modalities. These convergent results provide further impetus for examinations of the relationship of SLC12A5 with brain structure and function in neuropsychiatric disease.

In human electrophysiology research, the high gamma part of the power spectrum (similar to>60 Hz) is a relatively new area of investigation. Despite a low signal-to-noise ratio, evidence exists that it contains significant information about activity in local cortical networks. Here, using magnetoencephalography (MEG), we found high gamma activity when comparing data from an n-back working memory task to resting data in a large sample of normal volunteers. Initial analysis of power spectra from 0-back, 2-back, and rest trials showed three frequency bands exhibiting task-related differences: alpha, beta, and high gamma. Unlike alpha and beta, the high gamma spectrum was broad, without a peak at a single frequency. In addition, power in high gamma was highest for the 2-back and lowest during rest, while the opposite pattern occurred in the other bands. Beamformer source localization of each of the three frequency bands revealed a distinct set of sources for high gamma. These included several regions of prefrontal cortex that exhibited greater power when both n-back conditions were compared to rest. A subset of these regions had more power when the 2-back was compared to 0-back, which indicates a role in working memory performance. Our results show that high gamma will be important for understanding cortical processing during cognitive and other tasks. Furthermore, data from human intracortical recordings suggest that high gamma is the aggregate of spiking in local cortical networks, which implies that MEG could serve to bridge experimental modalities by noninvasively observing task-related modulation of spiking rates.

Top-down modulation of sensory processing is a critical neural mechanism subserving numerous important cognitive roles, one of which may be to inform lower-order sensory systems of the current 'task at hand' by conveying behavioral context to these systems. Accumulating evidence indicates that top-down cortical influences are carried by directed interareal synchronization of oscillatory neuronal populations, with recent results pointing to beta-frequency oscillations as particularly important for top-down processing. However, it remains to be determined if top-down beta-frequency oscillations indeed convey behavioral context. We measured spectral Granger Causality (sGC) using local field potentials recorded from microelectrodes chronically implanted in visual areas V1/V2, V4, and TEO of two rhesus macaque monkeys, and applied multivariate pattern analysis to the spatial patterns of top-down sGC. We decoded behavioral context by discriminating patterns of top-down (V4/TEO-to-V1/V2) beta-peak sGC for two different task rules governing correct responses to identical visual stimuli. The results indicate that top-down directed influences are carried to visual cortex by beta oscillations, and differentiate task demands even before visual stimulus processing. They suggest that top-down beta-frequency oscillatory processes coordinate processing of sensory information by conveying global knowledge states to early levels of the sensory cortical hierarchy independently of bottom-up stimulus-driven processing.

Temporal latencies between constituent parts of any complex dynamical system are an important factor in determining its function and stability. This is the case in neural systems where precise timing of the arrival of neural spikes and signals that propagate along axons is of fundamental importance. While there is a great effort today in reconstructing brain connectivity maps, there is no reliable framework to determine latencies in such networks. To address this, we use a multi-modality approach, combining structural and functional networks, each with their own latency estimates. Using structural information obtained with Mean Apparent Propagator (MAP) MRI data, we estimate the conduction velocity distribution (CVD) from the measured axon diameter distribution (AAD) in white matter. Using functional information derived from transcranial magnetic stimulation (TMS) evoked potentials (EP) via electroencephalography (EEG), as well as time-series analysis of magnetoencephalography (MEG) and electroencephalography (EEG) recordings, we have an independent assessment of the functional brain network connectivity. In deriving the latency from such data it is important to distinguish between temporal delays between neural events and true latencies that we aim to reconstruct. Latencies obtained from time-series analysis can only establish Wiener-Granger causality while incorporating TMS measurements and structural information about axons morphology and white matter path lengths provide additional information that can aid in deducing the true causality. In this abstract, we present preliminary results, with a focus on theoretical aspects and simulations. Our goal here is two-fold: 1) to improve the estimation of latencies from neural recordings using both stationary and non-stationary time-series approaches, and 2) to develop approaches that combine results from different imaging modalities, e.g., diffusion MRI and fMRI. In 1), we compare

Recently, independent components analysis (ICA) of resting state magnetoencephalography (MEG) recordings has revealed resting state networks (RSNs) that exhibit fluctuations of band-limited power envelopes. Most of the work in this area has concentrated on networks derived from the power envelope of beta bandpass-filtered data. Although research has demonstrated that most networks show maximal correlation in the beta band, little is known about how spatial patterns of correlations may differ across frequencies. This study analyzed MEG data from 18 healthy subjects to determine if the spatial patterns of RSNs differed between delta, theta, alpha, beta, gamma, and high gamma frequency bands. To validate our method, we focused on the sensorimotor network, which is well-characterized and robust in both MEG and functional magnetic resonance imaging (fMRI) resting state data. Synthetic aperture magnetometry (SAM) was used to project signals into anatomical source space separately in each band before a group temporal ICA was performed over all subjects and bands. This method preserved the inherent correlation structure of the data and reflected connectivity derived from single-band ICA, but also allowed identification of spatial spectral modes that are consistent across subjects. The implications of these results on our understanding of sensorimotor function are discussed, as are the potential applications of this technique. Hum Brain Mapp 38:779-791, 2017. © 2016 Wiley Periodicals, Inc.

Anxiety disorders affect approximately 1 in 5 (18%) Americans within a given 1 year period, placing a substantial burden on the national health care system. Therefore, there is a critical need to understand the neural mechanisms mediating anxiety symptoms. We used unbiased, multimodal, data-driven, whole-brain measures of neural activity (magnetoencephalography) and connectivity (fMRI) to identify the regions of the brain that contribute most prominently to sustained anxiety. We report that a single brain region, the intraparietal sulcus (IPS), shows both elevated neural activity and global brain connectivity during threat. The IPS plays a key role in attention orienting and may contribute to the hypervigilance that is a common symptom of pathological anxiety. Hyperactivation of this region during elevated state anxiety may account for the paradoxical facilitation of performance on tasks that require an external focus of attention, and impairment of performance on tasks that require an internal focus of attention.

Functional neuroimaging techniques including magnetoencephalography (MEG) have demonstrated that the brain is organized into networks displaying correlated activity. Group connectivity differences between healthy controls and participants with major depressive disorder (MDD) can be detected using temporal independent components analysis (ICA) on beta-bandpass filtered Hilbert envelope MEG data. However, the response of these networks to treatment is unknown. Ketamine, an N-methyl-D-aspartate (NMDA) receptor antagonist, exerts rapid antidepressant effects. We obtained MEG recordings before and after open-label infusion of 0.5mg/kg ketamine in MDD subjects (N=13) and examined networks previously shown to differ between healthy individuals and those with MDD. Connectivity between the amygdala and an insulo-temporal component decreased post-ketamine in MDD subjects towards that observed in control subjects at baseline. Decreased baseline connectivity of the subgenual anterior cingulate cortex (sgACC) with a bilateral precentral network had previously been observed in MDD compared to healthy controls, and the change in connectivity post-ketamine was proportional to the change in sgACC glucose metabolism in a subset (N=8) of subjects receiving [11F]FDG-PET imaging. Ketamine appeared to reduce connectivity, regardless of whether connectivity was abnormally high or low compared to controls at baseline. These preliminary findings suggest that sgACC connectivity may be directly related to glutamate levels.

Purpose: To describe and optimize an automated beamforming technique followed by identification of locations with excess kurtosis (g2) for efficient detection and localization of interictal spikes in patients with medically refractory epilepsy.Methods: Synthetic aperture magnetometry with g2 averaged over a sliding time window (SAMepi) was performed in seven patients with focal epilepsy and five healthy volunteers. The effect of varied window lengths on detection of spiking activity was evaluated.Results: Sliding window lengths of 0.5 to 10 seconds performed similarly, with 0.5- and 1-second windows detecting spiking activity in 1 of the 3 virtual sensor locations with highest kurtosis. These locations were concordant with the region of eventual surgical resection in these seven patients who remained seizure-free at 1 year. Average g2 values increased with increasing sliding window length in all subjects. In healthy volunteers, kurtosis values stabilized in data sets longer than 2 minutes.Conclusions: SAMepi using g2 averaged over 1-second sliding time windows in data sets of at least 2 minutes of duration reliably identified interictal spiking and the presumed seizure focus in these seven patients. Screening the five locations with highest kurtosis values for spiking activity is an efficient and accurate technique for localizing interictal activity using magnetoencephalography.Significance: SAMepi should be applied using the parameter values and procedure described for optimal detection and localization of interictal spikes. Use of this screening procedure could significantly improve the efficiency of magnetoencephalography analysis if clinically validated.

In typical magnetoencephalography and/or electroencephalography functional connectivity analysis, researchers select one of several methods that measure a relationship between regions to determine connectivity, such as coherence, power correlations, and others. However, it is largely unknown if some are more suited than others for various types of investigations. In this study, the authors investigate seven connectivity metrics to evaluate which, if any, are sensitive to audiovisual integration by contrasting connectivity when tracking an audiovisual object versus connectivity when tracking a visual object uncorrelated with the auditory stimulus. The authors are able to assess the metrics' performances at detecting audiovisual integration by investigating connectivity between auditory and visual areas. Critically, the authors perform their investigation on a whole-cortex all-to-all mapping, avoiding confounds introduced in seed selection. The authors find that amplitude-based connectivity measures in the beta band detect strong connections between visual and auditory areas during audiovisual integration, specifically between V4/V5 and auditory cortices in the right hemisphere. Conversely, phase-based connectivity measures in the beta band as well as phase and power measures in alpha, gamma, and theta do not show connectivity between audiovisual areas. The authors postulate that while beta power correlations detect audiovisual integration in the current experimental context, it may not always be the best measure to detect connectivity. Instead, it is likely that the brain utilizes a variety of mechanisms in neuronal communication that may produce differential types of temporal relationships.

Functional magnetic resonance imaging (fMRI) studies have revealed the existence of robust, interconnected brain networks exhibiting correlated low frequency fluctuations during rest, which can be derived by examining inherent spatio-temporal patterns in functional scans independent of any a priori model. In order to explore the electrophysiological underpinnings of these networks, analogous techniques have recently been applied to magnetoencephalography (MEG) data, revealing similar networks that exhibit correlated low frequency fluctuations in the power envelope of beta band (14-30 Hz) power. However, studies to date using this technique have concentrated on healthy subjects, and no method has yet been presented for group comparisons. We extended the ICA resting state MEG method to enable group comparisons, and demonstrate the technique in a sample of subjects with major depressive disorder (MDD). We found that the intrinsic resting state networks evident in fMRI appeared to be disrupted in individuals with MDD compared to healthy participants, particularly in the subgenual cingulate, although the electrophysiological correlates of this are unknown. Networks extracted from a combined group of healthy and MDD participants were examined for differences between groups. Individuals with MDD showed reduced correlations between the subgenual anterior cingulate (sgACC) and hippocampus in a network with primary nodes in the precentral and middle frontal gyri. Individuals with MDD also showed increased correlations between insulo-temporal nodes and amygdala compared to healthy controls. To further support our methods and findings, we present test/re-test reliability on independent recordings acquired within the same session. Our results demonstrate that group analyses are possible with the resting state MEG-independent component analysis (ICA) technique, highlighting a new pathway for analysis and discovery. This study also provides the first evidence of altered sgACC connectivity with a motor network. This finding, reliable across multiple sessions, suggests that the sgACC may partially mediate the psychomotor symptoms of MDD via synchronized changes in beta-band power, and expands the idea of the sgACC as a hub region mediating cognitive and emotional symptomatic domains in MDD. Findings of increased connectivity between the amygdala and cortical nodes further support the role of amygdalar networks in mediated depressive symptomatology. Published by Elsevier Inc.

This study investigated whether catechol-O-methyltransferase (COMT) Val/Met polymorphism was associated with variation in event-related desynchronization/synchronization (ERD/ERS) of responses during working memory (WM). 11 Val/Val and 11 Met/Met homozygous participants underwent magnetoencephalography (MEG) while performing a WM task. In contrast to small effects behaviourally, during the delay period Val/Val individuals showed lower ERS in the gamma band (Hz 30-50) in frontal regions, increased ERS in the alpha band (Hz 8-12) in the right frontal and parietal regions and increased ERD in the beta band (Hz 14-30) in the left fronto-temporal regions as compared with Met/Met homozygous individuals. During the response period Val/Val participants showed greater beta ERD in the prefrontal and parietotemporal regions. These results demonstrate that COMT genotype has a strong impact on brain responses (oscillatory activity) during WM performance likely a consequence of compensatory activity during the delay and response periods.

The processing of social information in the human brain is widely distributed neuroanatomically and finely orchestrated over time. However, a detailed account of the spatiotemporal organization of these key neural underpinnings of human social cognition remains to be elucidated. Here, we applied functional magnetic resonance imaging (fMRI) and magnetoencephalography (MEG) in the same participants to investigate spatial and temporal neural patterns evoked by viewing videos of facial muscle configurations. We show that observing the emergence of expressions elicits sustained blood oxygenation level-dependent responses in the superior temporal sulcus (STS), a region implicated in processing meaningful biological motion. We also found corresponding event-related changes in sustained MEG betaband (14-30 Hz) oscillatory activity in the STS, consistent with the possible role of beta-band activity in visual perception. Dynamically evolving fearful and happy expressions elicited early (0-400 ms) transient beta-band activity in sensorimotor cortex that persisted beyond 400 ms, at which time it became accompanied by a frontolimbic spread (400-1000 ms). In addition, individual differences in sustained STS beta-band activity correlated with speed of emotion recognition, substantiating the behavioral relevance of these signals. This STS beta-band activity showed valence-specific coupling with the time courses of facial movements as they emerged into full-blown fearful and happy expressions (negative and positive coupling, respectively). These data offer new insights into the perceptual relevance and orchestrated function of the STS and interconnected pathways in social-emotion cognition.

This chapter contains sections titled: Introduction Autonomous, Intermittent, Hierarchical Motions, from Brain Proteins Fluctuations to Emergent Magnetic Fields Magnetic Field Induction and Turbulence; Its Maintenance, Decay, and Modulation Localizing a Time-Varying Entropy Measure of Turbulence, Rank Vector Entropy (RVE) [35, 107], Using a Linearly Constrained Minimum Variance (LCMV) Beamformer Such as Synthetic Aperture Magnetometry (SAM) [25, 34], Yields State and Function- Related Localized Increases and Decreases in the RVE Estimate Potential Implications of the MHD Approach to MEG Magnetic Fields for Understanding the Mechanisms of Action and Clinical Applications of the Family of TMS (Transcranial Magnetic Stimulation) Human Brain Therapies

Complex networks have been observed to comprise small-world properties, believed to represent an optimal organization of local specialization and global integration of information processing at reduced wiring cost. Here, we applied magnitude squared coherence to resting magnetoencephalographic time series in reconstructed source space, acquired from controls and patients with schizophrenia, and generated frequency-dependent adjacency matrices modeling functional connectivity between virtual channels. After configuring undirected binary and weighted graphs, we found that all human networks demonstrated highly localized clustering and short characteristic path lengths. The most conservatively thresholded networks showed efficient wiring, with topographical distance between connected vertices amounting to one-third as observed in surrogate randomized topologies. Nodal degrees of the human networks conformed to a heavy-tailed exponentially truncated power-law, compatible with the existence of hubs, which included theta and alpha bilateral cerebellar tonsil, beta and gamma bilateral posterior cingulate, and bilateral thalamus across all frequencies. We conclude that all networks showed small-worldness, minimal physical connection distance, and skewed degree distributions characteristic of physically-embedded networks, and that these calculations derived from graph theoretical mathematics did not quantifiably distinguish between subject populations, independent of bandwidth. However, post-hoc measurements of edge computations at the scale of the individual vertex revealed trends of reduced gamma connectivity across the posterior medial parietal cortex in patients, an observation consistent with our prior resting activation study that found significant reduction of synthetic aperture magnetometry gamma power across similar regions. The basis of these small differences remains unclear.

What are the functional neuroimaging measurements required for more fully characterizing the events and locations of neocortical activity? A prime assumption has been that modulation of cortical activity will inevitably be reflected in changes in energy utilization (for the most part) changes of glucose and oxygen consumption. Are such a measures complete and sufficient? More direct measures of cortical electrophysiological activity show event or task-related modulation of amplitude or band-limited oscillatory power. Using magnetoencephalography (MEG), these measures have been shown to correlate well with energy utilization sensitive BOLD fMRI. In this paper, we explore the existence of state changes in electrophysiological cortical activity that can occur independently of changes in averaged amplitude, source power or indices of metabolic rates. In addition, we demonstrate that such state changes can be described by applying a new measure of complexity, rank vector entropy (RVE), to source waveform estimates from beamformer-processed MEG. RVE is a non-parametric symbolic dynamic informational entropy measure that accommodates the wide dynamic range of measured brain signals while resolving its temporal variations. By representing the measurements by their rank values, RVE overcomes the problem of defining embedding space partitions without resorting to signal compression. This renders RVE-independent of absolute signal amplitude. In addition, this approach is robust, being relatively free of tunable parameters. We present examples of task-free and task-dependent MEG demonstrating that RVE provides new information by uncovering hidden dynamical structure in the apparent turbulent (or chaotic) dynamics of spontaneous cortical activity.

Empirical studies over the past two decades have provided support for the hypothesis that schizophrenia is characterized by altered connectivity patterns in functional brain networks. These alterations have been proposed as genetically mediated diagnostic biomarkers and are thought to underlie altered cognitive functions such as working memory. However, the nature of this dysconnectivity remains far from understood. In this study, we perform an extensive analysis of functional connectivity patterns extracted from MEG data in 14 subjects with schizophrenia and 14 healthy controls during a 2-back working memory task. We investigate uni-, bi- and multivariate properties of sensor time series by computing wavelet entropy of and correlation between time series, and by constructing binary networks of functional connectivity both within and between classical frequency bands ([Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text]). Networks are based on the mutual information between wavelet time series, and estimated for each trial window separately, enabling us to consider both network topology and network dynamics. We observed significant decreases in time series entropy and significant increases in functional connectivity in the schizophrenia group in comparison to the healthy controls and identified an inverse relationship between these measures across both subjects and sensors that varied over frequency bands and was more pronounced in controls than in patients. The topological organization of connectivity was altered in schizophrenia specifically in high frequency [Formula: see text] and [Formula: see text] band networks as well as in the [Formula: see text]-[Formula: see text] cross-frequency networks. Network topology varied over trials to a greater extent in patients than in controls, suggesting disease-associated alterations in dynamic network properties of brain function. Our results identify signatures of aberrant neurophysiological behavior in schizophrenia across uni-, bi- and multivariate scales and lay the groundwork for further clinical studies that might lead to the discovery of new intermediate phenotypes.