Beamforming is a spatial filtering based source reconstruction method for EEG and MEG that allows the estimation of neuronal activity at a particular location within the brain. The computation of the location specific filter depends solely on an estimate of the data covariance matrix and on the forward model. Increasing the number of M/EEG sensors, increases the quantity of data required for accurate covariance matrix estimation. Often however we have a prior hypothesis about the site of, or the signal of interest. Here we show how this prior specification, in combination with optimal estimations of data dimensionality, can give enhanced beamformer performance for relatively short data segments. Specifically we show how temporal (Bayesian Principal Component Analysis) and spatial (lead field projection) methods can be combined to produce improvements in source estimation over and above employing the approaches individually.
•This paper concerns optimising beamformer analysis for anatomical ROIs.•Channel reduction is performed using an ROI projection and Bayesian PCA.•This improves covariance matrix estimation for a given data length.•The proposed approach results in improvements in source estimation.
Beamforming; Regions of interest; Bayesian PCA
This technical paper offers a critical re-evaluation of (spectral) Granger causality measures in the analysis of biological timeseries. Using realistic (neural mass) models of coupled neuronal dynamics, we evaluate the robustness of parametric and nonparametric Granger causality. Starting from a broad class of generative (state-space) models of neuronal dynamics, we show how their Volterra kernels prescribe the second-order statistics of their response to random fluctuations; characterised in terms of cross-spectral density, cross-covariance, autoregressive coefficients and directed transfer functions. These quantities in turn specify Granger causality — providing a direct (analytic) link between the parameters of a generative model and the expected Granger causality. We use this link to show that Granger causality measures based upon autoregressive models can become unreliable when the underlying dynamics is dominated by slow (unstable) modes — as quantified by the principal Lyapunov exponent. However, nonparametric measures based on causal spectral factors are robust to dynamical instability. We then demonstrate how both parametric and nonparametric spectral causality measures can become unreliable in the presence of measurement noise. Finally, we show that this problem can be finessed by deriving spectral causality measures from Volterra kernels, estimated using dynamic causal modelling.
•This paper describes the evaluation of expected Granger causality measures.•It uses these measures to quantify problems with dynamical instability and noise.•These problems are resolved by basing Granger measures on DCM estimates.
Granger causality; Dynamic causal modelling; Effective connectivity; Functional connectivity; Dynamics; Cross spectra; Neurophysiology
This paper tests the hypothesis that patients with schizophrenia have a deficit in selectively attending to predictable events. We used dynamic causal modeling (DCM) of electrophysiological responses – to predictable and unpredictable visual targets – to quantify the effective connectivity within and between cortical sources in the visual hierarchy in 25 schizophrenia patients and 25 age-matched controls. We found evidence for marked differences between normal subjects and schizophrenia patients in the strength of extrinsic backward connections from higher hierarchical levels to lower levels within the visual system. In addition, we show that not only do schizophrenia subjects have abnormal connectivity but also that they fail to adjust or optimize this connectivity when events can be predicted. Thus, the differential intrinsic recurrent connectivity observed during processing of predictable versus unpredictable targets was markedly attenuated in schizophrenia patients compared with controls, suggesting a failure to modulate the sensitivity of neurons responsible for passing sensory information of prediction errors up the visual cortical hierarchy. The findings support the proposed role of abnormal connectivity in the neuropathology and pathophysiology of schizophrenia.
Dynamic causal modeling; Schizophrenia; EEG; Prediction; Connectivity
The type I interferon (IFN) response protects cells from viral infection by inducing hundreds of interferon-stimulated genes (ISGs), some of which encode direct antiviral effectors1–3. Recent screening studies have begun to catalogue ISGs with antiviral activity against several RNA and DNA viruses4–13. However, antiviral ISG specificity across multiple distinct classes of viruses remains largely unexplored. Here we used an ectopic expression assay to screen a library of more than 350 human ISGs for effects on 14 viruses representing 7 families and 11 genera. We show that 47 genes inhibit one or more viruses, and 25 genes enhance virus infectivity. Comparative analysis reveals that the screened ISGs target positive-sense single-stranded RNA viruses more effectively than negative-sense single-stranded RNA viruses. Gene clustering highlights the cytosolic DNA sensor cyclic GMP-AMP synthase (cGAS, also known as MB21D1) as a gene whose expression also broadly inhibits several RNA viruses. In vitro, lentiviral delivery of enzymatically active cGAS triggers a STING-dependent, IRF3-mediated antiviral program that functions independently of canonical IFN/STAT1 signalling. In vivo, genetic ablation of murine cGAS reveals its requirement in the antiviral response to two DNA viruses, and an unappreciated contribution to the innate control of an RNA virus. These studies uncover new paradigms for the preferential specificity of IFN-mediated antiviral pathways spanning several virus families.
The integration of auditory feedback with vocal motor output is important for the control of voice fundamental frequency (F0). We used a pitch-shift paradigm where subjects respond to an alteration, or shift, of voice pitch auditory feedback with a reflexive change in F0. We presented varying magnitudes of pitch shifted auditory feedback to subjects during vocalization and passive listening and measured event related potentials (ERP’s) to the feedback shifts. Shifts were delivered at +100 and +400 cents (200 ms duration). The ERP data were modeled with Dynamic Causal Modeling (DCM) techniques where the effective connectivity between the superior temporal gyrus (STG), inferior frontal gyrus and premotor areas were tested. We compared three main factors; the effect of intrinsic STG connectivity, STG modulation across hemispheres and the specific effect of hemisphere. A Bayesian model selection procedure was used to make inference about model families. Results suggest that both intrinsic STG and left to right STG connections are important in the identification of self-voice error and sensory motor integration. We identified differences in left to right STG connections between 100 cent and 400 cent shift conditions suggesting that self and non-self voice error are processed differently in the left and right hemisphere. These results also highlight the potential of DCM modeling of ERP responses to characterize specific network properties of forward models of voice control.
Vocalization; Auditory feedback; ERP; DCM; Audio-vocal integration; Pitch shift
Movement is accompanied by changes in the degree to which neurons in corticobasal ganglia loops synchronize their activity within discrete frequency ranges. Although two principal frequency bands—beta (15–30 Hz) and gamma (60–90 Hz)—have been implicated in motor control, the precise functional correlates of their activities remain unclear. Local field potential (LFP) recordings in humans with Parkinson's disease undergoing surgery for deep brain stimulation to the subthalamic nucleus (STN) indicate that spectral changes both anticipate movement and occur perimovement. The extent to which such changes are modulated by cognitive factors involved in making a correct response seems critical in characterizing the functional associations of these oscillations. Accordingly, by recording LFP activity from the STN in parkinsonian patients, we demonstrate that perimovement beta and gamma reactivity is modulated by task complexity in a dopamine-dependent manner, despite the dynamics of the movement remaining unchanged. In contrast, spectral changes occurring in anticipation of future movement were limited to the beta band and, although modulated by dopaminergic therapy, were not modulated by task complexity. Our findings suggest two dopamine-dependent processes indexed by spectral changes in the STN: (1) an anticipatory activity reflected in the beta band that signals the likelihood of future action but does not proactively change with the cognitive demands of the potential response, and (2) perimovement activity that involves reciprocal beta and gamma band changes and is not exclusively related to explicit motor processing. Rather perimovement activity can also vary with, and may reflect, the cognitive complexity of the task.
Antiviral responses must be tightly regulated to rapidly defend against infection while minimizing inflammatory damage. Type 1 interferons (IFN-I) are crucial mediators of antiviral responses1 and their transcription is regulated by a variety of transcription factors2; principal amongst these is the family of interferon regulatory factors (IRFs)3. The IRF gene regulatory networks are complex and contain multiple feedback loops. The tools of systems biology are well suited to elucidate the complex interactions that give rise to precise coordination of the interferon response. Here we have used an unbiased systems approach to predict that a member of the forkhead family of transcription factors, FOXO3, is a negative regulator of a subset of antiviral genes. This prediction was validated using macrophages isolated from Foxo3-null mice. Genome-wide location analysis combined with gene deletion studies identified the Irf7 gene as a critical target of FOXO3. FOXO3 was identified as a negative regulator of Irf7 transcription and we have further demonstrated that FOXO3, IRF7 and IFN-I form a coherent feed-forward regulatory circuit. Our data suggest that the FOXO3-IRF7 regulatory circuit represents a novel mechanism for establishing the requisite set points in the interferon pathway that balances the beneficial effects and deleterious sequelae of the antiviral response.
Functional neurosurgical techniques provide a unique opportunity to explore patterns of interaction between the cerebral cortex and basal ganglia in patients with Parkinson's disease (PD). Previous work using simultaneous magnetoencephalographic (MEG) and local field potential (LFP) recordings from the region of the subthalamic nucleus (STNr) has characterised resting patterns of connectivity in the alpha and beta frequency bands and their modulation by dopaminergic medication. Recently we have also characterised the effect of movement on patterns of gamma band coherence between the STNr and cortical sites. Here we specifically investigate how the prominent coherence between the STNr and temporal cortex in the alpha band is modulated by movement both on and off dopaminergic medication in patients following the insertion of Deep Brain Stimulation (DBS) electrodes. We show that movement is associated with a suppression of local alpha power in the temporal cortex and STNr that begins about 2 s prior to a self-paced movement and is independent of dopaminergic status. In contrast, the peak reduction in coherence between these sites occurs after movement onset and is more marked in the on than in the off dopaminergic medication state. The difference in alpha band coherence on and off medication was found to correlate with the drug related improvement in clinical parameters. Overall, the movement-related behaviour of activities in the alpha band in patients with PD serves to highlight the role of dopamine in modulating large-scale, interregional synchronisation.
► We studied subthalamo-cortical coherence in the alpha band in Parkinson's patients. ► A decrease in coherence occurred with movement facilitated by dopamine. ► This effect correlated with clinical improvement.
Oscillations; Human; Intracranial recordings
Functional neurosurgery has afforded the opportunity to assess interactions between populations of neurons in the human cerebral cortex and basal ganglia in patients with Parkinson’s disease (PD). Interactions occur over a wide range of frequencies, and the functional significance of those above 30 Hz is particularly unclear. Do they improve movement and, if so, in what way? We acquired simultaneously magnetoencephalography (MEG) and direct recordings from the subthalamic nucleus (STN) in 17 PD patients. We examined the effect of synchronous and sequential finger movements and of the dopamine prodrug levodopa on induced power in the contralateral primary motor cortex (M1) and STN and on the coherence between the two structures. We observed discrete peaks in M1 and STN power over 60-90 Hz and 300-400 Hz. All these power peaks increased with movement and levodopa treatment. Only STN activity over 60-90 Hz was coherent with activity in M1. Directionality analysis showed that STN gamma activity at 60-90 Hz tended to drive gamma activity in M1. The effects of levodopa on both local and distant synchronisation over 60-90 Hz correlated with the degree of improvement in bradykinesia-rigidity, as did local STN activity at 300-400 Hz. Despite this, there were no effects of movement type, nor interactions between movement type and levodopa in the STN, nor in the coherence between STN and M1. We conclude that synchronisation over 60-90 Hz in the basal ganglia cortical network is prokinetic, but likely through a modulatory effect rather than any involvement in explicit motor processing.
Oscillatory activity in the beta frequency band has been shown to be modulated during the preparation and execution of voluntary movements at both cortical and subcortical levels. The exaggeration of beta activity in the basal ganglia of patients with Parkinson's disease has heightened interest in this phenomenon. However, the precise function, if any, subserved by modulations in beta activity remains unclear. Here we test the hypothesis that beta reactivity can be dissociated from processing of specific actions and can index the salience of cues with respect to future behavior in a way that might help prospectively prioritize resources. To this end we used an experimental paradigm designed to dissociate salient warning cues from processing of specific motor or cognitive actions. We recorded local field potential activity from the subthalamic nucleus of humans undergoing functional neurosurgery for the treatment of Parkinson's disease, while the same patients were on or off the dopamine prodrug levodopa. In this way we demonstrate that beta reactivity is indeed dependent on the salience of cues with respect to future motor and cognitive action and is promoted by dopamine. The loss of normal beta encoding of saliency may underlie some of the motor and cognitive features of basal ganglia disorders such as Parkinson's disease.
Magnetoencephalographic (MEG) recordings are a rich source of information about the neural dynamics underlying cognitive processes in the brain, with excellent temporal and good spatial resolution. In recent years there have been considerable advances in MEG hardware developments and methods. Sophisticated analysis techniques are now routinely applied and continuously improved, leading to fascinating insights into the intricate dynamics of neural processes. However, the rapidly increasing level of complexity of the different steps in a MEG study make it difficult for novices, and sometimes even for experts, to stay aware of possible limitations and caveats. Furthermore, the complexity of MEG data acquisition and data analysis requires special attention when describing MEG studies in publications, in order to facilitate interpretation and reproduction of the results. This manuscript aims at making recommendations for a number of important data acquisition and data analysis steps and suggests details that should be specified in manuscripts reporting MEG studies. These recommendations will hopefully serve as guidelines that help to strengthen the position of the MEG research community within the field of neuroscience, and may foster discussion in order to further enhance the quality and impact of MEG research.
Magnetoencephalography; MEG; Acquisition; Analysis; Connectivity; Source localization; Guidelines; Recommendations; Reproducible research; Spectral analysis
In Kilner et al. [Kilner, J.M., Kiebel, S.J., Friston, K.J., 2005. Applications of random field theory to electrophysiology. Neurosci. Lett. 374, 174–178.] we described a fairly general analysis of induced responses—in electromagnetic brain signals—using the summary statistic approach and statistical parametric mapping. This involves localising induced responses—in peristimulus time and frequency—by testing for effects in time–frequency images that summarise the response of each subject to each trial type. Conventionally, these time–frequency summaries are estimated using post‐hoc averaging of epoched data. However, post‐hoc averaging of this sort fails when the induced responses overlap or when there are multiple response components that have variable timing within each trial (for example stimulus and response components associated with different reaction times). In these situations, it is advantageous to estimate response components using a convolution model of the sort that is standard in the analysis of fMRI time series. In this paper, we describe one such approach, based upon ordinary least squares deconvolution of induced responses to input functions encoding the onset of different components within each trial. There are a number of fundamental advantages to this approach: for example; (i) one can disambiguate induced responses to stimulus onsets and variably timed responses; (ii) one can test for the modulation of induced responses—over peristimulus time and frequency—by parametric experimental factors and (iii) one can gracefully handle confounds—such as slow drifts in power—by including them in the model. In what follows, we consider optimal forms for convolution models of induced responses, in terms of impulse response basis function sets and illustrate the utility of deconvolution estimators using simulated and real MEG data.
► We propose a new approach to analysis of induced responses in M/EEG. ► The General Linear Model is used to model continuous power as in fMRI 1st-level. ► The results can be presented as conventional time–frequency images. ► Our method is better for experiments with variable timing and overlapping events.
EEG; MEG; ERSP; Induced responses; Time–frequency analysis; General linear model; Convolution; Statistical parametric mapping
Stroke results in reorganization of residual brain networks. The functional role of brain regions within these networks remains unclear, particularly those in the contralesional hemisphere. We studied 25 stroke patients with a range of motor impairment and 23 healthy age-matched controls using magnetoencephalography (MEG) and electromyography (EMG) to measure oscillatory signals from the brain and affected muscles simultaneously during a simple isometric hand grip, from which cortico-muscular coherence (CMC) was calculated. Peaks of cortico-muscular coherence in both the beta and gamma bands were found in the contralateral sensorimotor cortex in all healthy controls, but were more widespread in stroke patients, including some peaks found in the contralesional hemisphere (7 patients for beta coherence and 5 for gamma coherence). Neither the coherence value nor the distance of the coherence peak from the mean of controls correlated with impairment. Peak CMC in the contralesional hemisphere was found not only in some highly impaired patients, but also in some patients with good functional recovery. Our results provide evidence that a wide range of cortical brain regions, including some in the contralesional hemisphere, may have influence over EMG activity in the affected muscles after stroke thereby supporting functional recovery.
► We examined cortico-muscular coherence location in stroke patients and controls. ► The location of peak coherence was more widely distributed in stroke patients. ► In some patients, peak coherence was found in the contralesional hemisphere. ► The location of coherence in patients did not correlate with impairment. ► Contralesional hemisphere can support functional motor recovery after stroke.
CMC, cortico-muscular coherence; MEG, magnetoencephalography; EMG, electromyography; fMRI, functional magnetic resonance imaging; TMS, transcranial magnetic stimulation; M1, primary motor cortex; PCA, principal component analysis; MVC, maximum voluntary contraction; DICS, dynamic imaging of coherent sources; Magnetoencephalography; Cortico-muscular coherence; Stroke recovery; Motor; Brain
A multimodal neuroimaging study of virtual spatial navigation extends the role of the hippocampal theta rhythm to human memory and self-directed learning.
The hippocampus is crucial for episodic or declarative memory and the theta rhythm has been implicated in mnemonic processing, but the functional contribution of theta to memory remains the subject of intense speculation. Recent evidence suggests that the hippocampus might function as a network hub for volitional learning. In contrast to human experiments, electrophysiological recordings in the hippocampus of behaving rodents are dominated by theta oscillations reflecting volitional movement, which has been linked to spatial exploration and encoding. This literature makes the surprising cross-species prediction that the human hippocampal theta rhythm supports memory by coordinating exploratory movements in the service of self-directed learning. We examined the links between theta, spatial exploration, and memory encoding by designing an interactive human spatial navigation paradigm combined with multimodal neuroimaging. We used both non-invasive whole-head Magnetoencephalography (MEG) to look at theta oscillations and Functional Magnetic Resonance Imaging (fMRI) to look at brain regions associated with volitional movement and learning. We found that theta power increases during the self-initiation of virtual movement, additionally correlating with subsequent memory performance and environmental familiarity. Performance-related hippocampal theta increases were observed during a static pre-navigation retrieval phase, where planning for subsequent navigation occurred. Furthermore, periods of the task showing movement-related theta increases showed decreased fMRI activity in the parahippocampus and increased activity in the hippocampus and other brain regions that strikingly overlap with the previously observed volitional learning network (the reverse pattern was seen for stationary periods). These fMRI changes also correlated with participant's performance. Our findings suggest that the human hippocampal theta rhythm supports memory by coordinating exploratory movements in the service of self-directed learning. These findings directly extend the role of the hippocampus in spatial exploration in rodents to human memory and self-directed learning.
Neural activity both within and across brain regions can oscillate in different frequency ranges (such as alpha, gamma, and theta frequencies), and these different ranges are associated with distinct functions. In behaving rodents, for example, theta rhythms (4–12 Hz) in the hippocampus are prominent during the initiation of movement and have been linked to spatial exploration. Recent evidence in humans, however, suggests that the human hippocampus is involved in guiding self-directed learning. This suggests that the human hippocampal theta rhythm supports memory by coordinating exploratory movements in the service of self-directed learning. In this study, we tested whether there is a human analogue for the movement-initiation-related theta rhythm found in the rodent hippocampus by using a virtual navigation paradigm, combined with non-invasive recordings and functional imaging techniques. Our recordings showed that, indeed, theta power increases are linked to movement initiation. We also examined the relationship to memory encoding, and we found that hippocampal theta oscillations related to pre-retrieval planning predicted memory performance. Imaging results revealed that periods of the task showing movement-related theta also showed increased activity in the hippocampus, as well as other brain regions associated with self-directed learning. These findings directly extend the role of the hippocampal theta rhythm in rodent spatial exploration to human memory and self-directed learning.
Statistical parametric mapping (SPM) locates significant clusters based on a ratio of signal to noise (a ‘contrast’ of the parameters divided by its standard error) meaning that very low noise regions, for example outside the brain, can attain artefactually high statistical values. Similarly, the commonly applied preprocessing step of Gaussian spatial smoothing can shift the peak statistical significance away from the peak of the contrast and towards regions of lower variance. These problems have previously been identified in positron emission tomography (PET) (Reimold et al., 2006) and voxel-based morphometry (VBM) (Acosta-Cabronero et al., 2008), but can also appear in functional magnetic resonance imaging (fMRI) studies. Additionally, for source-reconstructed magneto- and electro-encephalography (M/EEG), the problems are particularly severe because sparsity-favouring priors constrain meaningfully large signal and variance to a small set of compactly supported regions within the brain. (Acosta-Cabronero et al., 2008) suggested adding noise to background voxels (the ‘haircut’), effectively increasing their noise variance, but at the cost of contaminating neighbouring regions with the added noise once smoothed. Following theory and simulations, we propose to modify – directly and solely – the noise variance estimate, and investigate this solution on real imaging data from a range of modalities.
► Statistical parametric mapping judges significance with a signal-to-noise ratio. ► Low noise, e.g. outside the brain, can yield artefactually high statistical values. ► Spatial smoothing can shift peaks substantially towards regions of low variance. ► Source-reconstructed M/EEG data exhibits the problem particularly severely. ► The problem can be addressed by modifying the noise variance estimate.
EEG, electroencephalography; fMRI, functional magnetic resonance imaging; FWHM, full-width at half-maximum; GM, grey matter; MEG, magnetoencephalography; MIP, maximum intensity projection; MNI, Montreal Neurological Institute; ResMS, residual mean squares; SPM, statistical parametric mapping; PET, positron emission tomography; VBM, voxel-based morphometry; SPM; Low variance; VBM; MEG; EEG; Source reconstruction
Optimal decision-making is guided by evaluating the outcomes of previous decisions. Prediction errors are theoretical teaching signals which integrate two features of an outcome: its inherent value and prior expectation of its occurrence. To uncover the magnetic signature of prediction errors in the human brain we acquired magnetoencephalographic (MEG) data while participants performed a gambling task. Our primary objective was to use formal criteria, based upon an axiomatic model (Caplin and Dean, 2008a), to determine the presence and timing profile of MEG signals that express prediction errors. We report analyses at the sensor level, implemented in SPM8, time locked to outcome onset. We identified, for the first time, a MEG signature of prediction error, which emerged approximately 320 ms after an outcome and expressed as an interaction between outcome valence and probability. This signal followed earlier, separate signals for outcome valence and probability, which emerged approximately 200 ms after an outcome. Strikingly, the time course of the prediction error signal, as well as the early valence signal, resembled the Feedback-Related Negativity (FRN). In simultaneously acquired EEG data we obtained a robust FRN, but the win and loss signals that comprised this difference wave did not comply with the axiomatic model. Our findings motivate an explicit examination of the critical issue of timing embodied in computational models of prediction errors as seen in human electrophysiological data.
► We identified, for the first time, an MEG signature of a human prediction error. ► The waveform resembled the Feedback-Related Negativity (FRN) signal in EEG. ► MEG effects of probability and valence were emerged before the prediction error signals, 200 ms after the outcome. ► The EEG data revealed classic FRN which was modulated by probability.
Decision-making; Prediction error; Reward; MEG; Feedback-related negativity; Error-related negativity
Cortico-basal ganglia-thalamocortical circuits are severely disrupted by the dopamine depletion of Parkinson's disease (PD), leading to pathologically exaggerated beta oscillations. Abnormal rhythms, found in several circuit nodes are correlated with movement impairments but their neural basis remains unclear. Here, we used dynamic causal modelling (DCM) and the 6-hydroxydopamine-lesioned rat model of PD to examine the effective connectivity underlying these spectral abnormalities. We acquired auto-spectral and cross-spectral measures of beta oscillations (10–35 Hz) from local field potential recordings made simultaneously in the frontal cortex, striatum, external globus pallidus (GPe) and subthalamic nucleus (STN), and used these data to optimise neurobiologically plausible models. Chronic dopamine depletion reorganised the cortico-basal ganglia-thalamocortical circuit, with increased effective connectivity in the pathway from cortex to STN and decreased connectivity from STN to GPe. Moreover, a contribution analysis of the Parkinsonian circuit distinguished between pathogenic and compensatory processes and revealed how effective connectivity along the indirect pathway acquired a strategic importance that underpins beta oscillations. In modelling excessive beta synchrony in PD, these findings provide a novel perspective on how altered connectivity in basal ganglia-thalamocortical circuits reflects a balance between pathogenesis and compensation, and predicts potential new therapeutic targets to overcome dysfunctional oscillations.
Parkinson's disease is a progressive age-related neurodegenerative disorder that severely disrupts movement. The major pathology in Parkinson's disease is the degeneration of a group of neurons that contain a chemical known as dopamine. Treatment of Parkinsonism includes pharmacological interventions that aim to replace dopamine and more recently, implanted devices that aim to restore movement through electrical stimulation of the brain's movement circuits. Understanding the electrical properties that emerge as a result of depleted dopamine may reveal new avenues for developing these technologies. By combining a novel model-based approach with multi-site electrophysiological recordings from an animal model of Parkinson's disease we provide empirical evidence for a link between abnormal electrical activity in the Parkinsonian brain and its physiological basis. We have examined the connections along the brain's motor circuits, and found an abnormality in inter-area connections in a particular neural pathway, a pathway critically dependent on dopamine. The scheme makes strong and testable predictions about which neural pathways are significantly altered in the pathological state and so represent empirically motivated therapeutic targets.
We address the problem of controlling false positive rates in mass-multivariate tests for electromagnetic responses in compact regions of source space. We show that mass-univariate thresholds based on sensor level multivariate thresholds (approximated using Roy's union–intersection principle) are unduly conservative. We then consider a Bonferroni correction for source level tests based on the number of unique lead-field extrema. For a given source space, the sensor indices corresponding to the maxima and minima (for each dipolar lead field) are listed, and the number of unique extrema is given by the number of unique pairs in this list. Using a multivariate beamformer formulation, we validate this heuristic against empirical permutation thresholds for mass-univariate and mass-multivariate tests (of induced and evoked responses) for a variety of source spaces, using simulated and real data. We also show that the same approximations hold when dealing with a cortical manifold (rather than a volume) and for mass-multivariate minimum norm solutions. We demonstrate that the mass-multivariate framework is not restricted to tests on a single contrast of effects (cf, Roy's maximum root) but also accommodates multivariate effects (cf, Wilk's lambda).
► We aim to estimate the number of independent tests made in MEG source space. ► We trial a heuristic based on the number of unique lead field extrema. ► We compare Bonferroni corrected tests using this heuristic to permutation methods. ► The heuristic performs well for both mass-univariate and mass-multivariate tests.
SPM is a free and open source software written in MATLAB (The MathWorks, Inc.). In addition to standard M/EEG preprocessing, we presently offer three main analysis tools: (i) statistical analysis of scalp-maps, time-frequency images, and volumetric 3D source reconstruction images based on the general linear model, with correction for multiple comparisons using random field theory; (ii) Bayesian M/EEG source reconstruction, including support for group studies, simultaneous EEG and MEG, and fMRI priors; (iii) dynamic causal modelling (DCM), an approach combining neural modelling with data analysis for which there are several variants dealing with evoked responses, steady state responses (power spectra and cross-spectra), induced responses, and phase coupling. SPM8 is integrated with the FieldTrip toolbox , making it possible for users to combine a variety of standard analysis methods with new schemes implemented in SPM and build custom analysis tools using powerful graphical user interface (GUI) and batching tools.
We review recent methodological developments within a parametric empirical Bayesian (PEB) framework for reconstructing intracranial sources of extracranial electroencephalographic (EEG) and magnetoencephalographic (MEG) data under linear Gaussian assumptions. The PEB framework offers a natural way to integrate multiple constraints (spatial priors) on this inverse problem, such as those derived from different modalities (e.g., from functional magnetic resonance imaging, fMRI) or from multiple replications (e.g., subjects). Using variations of the same basic generative model, we illustrate the application of PEB to three cases: (1) symmetric integration (fusion) of MEG and EEG; (2) asymmetric integration of MEG or EEG with fMRI, and (3) group-optimization of spatial priors across subjects. We evaluate these applications on multi-modal data acquired from 18 subjects, focusing on energy induced by face perception within a time–frequency window of 100–220 ms, 8–18 Hz. We show the benefits of multi-modal, multi-subject integration in terms of the model evidence and the reproducibility (over subjects) of cortical responses to faces.
source reconstruction; bioelectromagnetic signals; data fusion; neuroimaging
Motivation: Histone acetylation (HAc) is associated with open chromatin, and HAc has been shown to facilitate transcription factor (TF) binding in mammalian cells. In the innate immune system context, epigenetic studies strongly implicate HAc in the transcriptional response of activated macrophages. We hypothesized that using data from large-scale sequencing of a HAc chromatin immunoprecipitation assay (ChIP-Seq) would improve the performance of computational prediction of binding locations of TFs mediating the response to a signaling event, namely, macrophage activation.
Results: We tested this hypothesis using a multi-evidence approach for predicting binding sites. As a training/test dataset, we used ChIP-Seq-derived TF binding site locations for five TFs in activated murine macrophages. Our model combined TF binding site motif scanning with evidence from sequence-based sources and from HAc ChIP-Seq data, using a weighted sum of thresholded scores. We find that using HAc data significantly improves the performance of motif-based TF binding site prediction. Furthermore, we find that within regions of high HAc, local minima of the HAc ChIP-Seq signal are particularly strongly correlated with TF binding locations. Our model, using motif scanning and HAc local minima, improves the sensitivity for TF binding site prediction by ∼50% over a model based on motif scanning alone, at a false positive rate cutoff of 0.01.
Availability: The data and software source code for model training and validation are freely available online at http://magnet.systemsbiology.net/hac.
Contact: firstname.lastname@example.org; email@example.com
Supplementary information: Supplementary data are available at Bioinformatics online.
Insight into how brain structures interact is critical for understanding the principles of functional brain architectures and may lead to better diagnosis and therapy for neuropsychiatric disorders. We recorded, simultaneously, magnetoencephalographic (MEG) signals and subcortical local field potentials (LFP) in a Parkinson's disease (PD) patient with bilateral deep brain stimulation (DBS) electrodes in the subthalamic nucleus (STN). These recordings offer a unique opportunity to characterize interactions between the subcortical structures and the neocortex. However, high-amplitude artefacts appeared in the MEG. These artefacts originated from the percutaneous extension wire, rather than from the actual DBS electrode and were locked to the heart beat. In this work, we show that MEG beamforming is capable of suppressing these artefacts and quantify the optimal regularization required. We demonstrate how beamforming makes it possible to localize cortical regions whose activity is coherent with the STN-LFP, extract artefact-free virtual electrode time-series from regions of interest and localize cortical areas exhibiting specific task-related power changes. This furnishes results that are consistent with previously reported results using artefact-free MEG data. Our findings demonstrate that physiologically meaningful information can be extracted from heavily contaminated MEG signals and pave the way for further analysis of combined MEG-LFP recordings in DBS patients.
The innate immune system is a two-edged sword; it is absolutely required for host defense against infection but, uncontrolled, can trigger a plethora of inflammatory diseases. Here we used systems biology approaches to predict and validate a gene regulatory network involving a dynamic interplay between the transcription factors NF-κB, C/EBPδ, and ATF3 that controls inflammatory responses. We mathematically modeled transcriptional regulation of Il6 and Cebpd genes and experimentally validated the prediction that the combination of an initiator (NF-κB), an amplifier (C/EBPδ) and an attenuator (ATF3) forms a regulatory circuit that discriminates between transient and persistent Toll-like receptor 4-induced signals. Our results suggest a mechanism that enables the innate immune system to detect the duration of infection and to respond appropriately.
The aim of this paper is to describe a simple procedure for
electromagnetic (EEG or MEG) source reconstruction, in the context of group
studies. This entails a simple extension of existing source reconstiruction
techniques based upon the inversion of hierarchical models. The extension
ensures that evoked or induced responses are reconstructed in the same subset of
sources, over subjects. Effectively, the procedure aligns the deployment of
reconstructed activity over subjects and increases, substantially, the detection
of differences between evoked or induced responses at the group or
Hierarchical Bayes; Source reconstruction; EEG; MEG; Restricted maximum likelihood; Automatic relevance determination