Laura Marzetti

Temporal dynamics of spontaneous MEG activity in brain networks

Functional MRI (fMRI) studies have shown that low-frequency (<0.1 Hz) spontaneous fluctuations of the blood oxygenation level dependent (BOLD) signal during restful wakefulness are coherent within distributed large-scale cortical and subcortical networks (resting state networks, RSNs). The neuronal mechanisms underlying RSNs remain poorly understood. Here, we describe magnetoencephalographic correspondents of two well-characterized RSNs: the dorsal attention and the default mode networks. Seed-based correlation mapping was performed using time-dependent MEG power reconstructed at each voxel within the brain. The topography of RSNs computed on the basis of extended (5 min) epochs was similar to that observed with fMRI but confined to the same hemisphere as the seed region. Analyses taking into account the nonstationarity of MEG activity showed transient formation of more complete RSNs, including nodes in the contralateral hemisphere. Spectral analysis indicated that RSNs manifest in MEG as synchronous modulation of band-limited power primarily within the theta, alpha, and beta bands—that is, in frequencies slower than those associated with the local electrophysiological correlates of event-related BOLD responses.

A Cortical Core for Dynamic Integration of Functional Networks in the Resting Human Brain

We used magneto-encephalography to study the temporal dynamics of band-limited power correlation at rest within and across six brain networks previously defined by prior functional magnetic resonance imaging (fMRI) studies. Epochs of transiently high within-network band limited power (BLP) correlation were identified and correlation of BLP time-series across networks was assessed in these epochs. These analyses demonstrate that functional networks are not equivalent with respect to cross-network interactions. The default-mode network and the posterior cingulate cortex, in particular, exhibit the highest degree of transient BLP correlation with other networks especially in the 14-25 Hz (β band) frequency range. Our results indicate that the previously demonstrated neuroanatomical centrality of the PCC and DMN has a physiological counterpart in the temporal dynamics of network interaction at behaviorally relevant timescales. This interaction involved subsets of nodes from other networks during periods in which their internal correlation was low.

Frequency specific interactions of MEG resting state activity within and across brain networks as revealed by the multivariate interaction measure.

Resting state networks (RSNs) are sets of brain regions exhibiting temporally coherent activity fluctuations in the absence of imposed task structure. RSNs have been extensively studied with fMRI in the infra-slow frequency range (nominally <10(-1)Hz). The topography of fMRI RSNs reflects stationary temporal correlation over minutes. However, neuronal communication occurs on a much faster time scale, at frequencies nominally in the range of 10(0)-10(2)Hz. We examined phase-shifted interactions in the delta (2-3.5 Hz), theta (4-7 Hz), alpha (8-12 Hz) and beta (13-30 Hz) frequency bands of resting-state source space MEG signals. These analyses were conducted between nodes of the dorsal attention network (DAN), one of the most robust RSNs, and between the DAN and other networks. Phase shifted interactions were mapped by the multivariate interaction measure (MIM), a measure of true interaction constructed from the maximization of imaginary coherency in the virtual channels comprised of voxel signals in source space. Non-zero-phase interactions occurred between homologous left and right hemisphere regions of the DAN in the delta and alpha frequency bands. Even stronger non-zero-phase interactions were detected between networks. Visual regions bilaterally showed phase-shifted interactions in the alpha band with regions of the DAN. Bilateral somatomotor regions interacted with DAN nodes in the beta band. These results demonstrate the existence of consistent, frequency specific phase-shifted interactions on a millisecond time scale between cortical regions within RSN as well as across RSNs.

The use of standardized infinity reference in EEG coherency studies

The study of large scale interactions in the brain from EEG signals is a promising method for the identification of functional networks. However, the validity of a large scale parameter is limited by two factors: the use of a non-neutral reference and the artifactual self-interactions between the measured EEG signals introduced by volume conduction. In this paper, we propose an approach to study large scale EEG coherency in which these factors are eliminated. Artifactual self-interaction by volume conduction is eliminated by using the imaginary part of the complex coherency as a measure of interaction and the Reference Electrode Standardization Technique (REST) is used for the approximate standardization of the reference of scalp EEG recordings to a point at infinity that, being far from all possible neural sources, acts like a neutral virtual reference. The application of our approach to simulated and real EEG data shows that the detection of interaction, as opposed to artifacts due to reference and volume conduction, is a goal that can be achieved from the study of a large scale parameter.

Adding dynamics to the Human Connectome Project with MEG

The Human Connectome Project (HCP) seeks to map the structural and functional connections between network elements in the human brain. Magnetoencephalography (MEG) provides a temporally rich source of information on brain network dynamics and represents one source of functional connectivity data to be provided by the HCP. High quality MEG data will be collected from 50 twin pairs both in the resting state and during performance of motor, working memory and language tasks. These data will be available to the general community. Additionally, using the cortical parcellation scheme common to all imaging modalities, the HCP will provide processing pipelines for calculating connection matrices as a function of time and frequency. Together with structural and functional data generated using magnetic resonance imaging methods, these data represent a unique opportunity to investigate brain network connectivity in a large cohort of normal adult human subjects. The analysis pipeline software and the dynamic connectivity matrices that it generates will all be made freely available to the research community.

Estimating true brain connectivity from EEG/MEG data invariant to linear and static transformations in sensor space

The imaginary part of coherency is a measure to investigate the synchronization of brain sources on the EEG/MEG sensor level, robust to artifacts of volume conduction meaning that independent sources cannot generate a significant result. It does not mean, however, that volume conduction is irrelevant when true interactions are present. Here, we analyze in detail the possibilities to construct measures of true brain interactions which are strictly invariant to linear spatial transformations of the sensor data. Specifically, such measures can be constructed from maximization of imaginary coherency in virtual channels, bivariate measures as a corrected variate of imaginary coherence, and global measures indicating the total interaction contained within a space or between two spaces. A complete theoretic framework on this question is provided for second order statistical moments. Relations to existing linear and nonlinear approaches are presented. We applied the methods to resting state EEG data, showing clear interactions at all bands, and to a combined measurement of EEG and MEG during rest condition and a finger tapping task. We found that MEG was capable of observing brain interactions which were not observable in the EEG data.

Understanding brain connectivity from EEG data by identifying systems composed of interacting sources.

In understanding and modeling brain functioning by EEG/MEG, it is not only important to be able to identify active areas but also to understand interference among different areas. The EEG/MEG signals result from the superimposition of underlying brain source activities volume conducted through the head. The effects of volume conduction produce spurious interactions in the measured signals. It is fundamental to separate true source interactions from noise and to unmix the contribution of different systems composed by interacting sources in order to understand interference mechanisms. As a prerequisite, we consider the problem of unmixing the contribution of uncorrelated sources to a measured field. This problem is equivalent to the problem of unmixing the contribution of different uncorrelated compound systems composed by interacting sources. To this end, we develop a principal component analysis-based method, namely, the source principal component analysis (sPCA), which exploits the underlying assumption of orthogonality for sources, estimated from linear inverse methods, for the extraction of essential features in signal space. We then consider the problem of demixing the contribution of correlated sources that comprise each of the compound systems identified by using sPCA. While the sPCA orthogonality assumption is sufficient to separate uncorrelated systems, it cannot separate the individual components within each system. To address that problem, we introduce the Minimum Overlap Component Analysis (MOCA), employing a pure spatial criterion to unmix pairs of correlates (or coherent) sources. The proposed methods are tested in simulations and applied to EEG data from human micro and alpha rhythms.

Local and remote effects of transcranial direct current stimulation on the electrical activity of the motor cortical network

We systematically investigated the effects of cathodal and anodal Transcranial Direct Current Stimulation (CtDCS, AtDCS) on the electric activity of primary motor cortex during a motor task. High‐density electroencephalography was used to define the spatial diffusion of tDCS after effects. Ten healthy subjects performed a finger tapping task with the right hand before and after three separate sessions of 20 minutes of Sham, AtDCS or CtDCS over left primary motor cortex (M1). During movement, we found an increment of low alpha band Event‐Related Desynchronization (ERD) in bilateral central, frontal areas and in the left inferior parietal region, as well as an increment of beta ERD in fronto‐central and parieto‐occipital regions, after AtDCs compared to Sham and CtDCS. In the rest pre‐movement period, after Sham as well as AtDCS, we documented an increment of low alpha band power over the course of pre‐ and post‐stimulation recording sessions, localized in the sensorimotor and parieto‐occipital regions. On the contrary, after CtDCS no increment of low alpha power was found. Finally beta band coherence among signals from left sensorimotor cortex and activity of bilateral parietal, occipital and right frontal regions was higher after AtDCS compared with Sham condition. Similarly, theta coherence with parietal and frontal regions was enhanced after AtDCS. We hypothesize that the local modulation of membrane polarization, as well as long‐lasting synaptic modification induced by tDCS over M1, could result in changes of both local band power and functional architecture of the motor network. Hum Brain Mapp 35:2220–2232, 2014.

A signal-processing pipeline for magnetoencephalography resting-state networks.

To study functional connectivity using magnetoencephalographic (MEG) data, the high-quality source-level reconstruction of brain activity constitutes a critical element. MEG resting-state networks (RSNs) have been documented by means of a dedicated processing pipeline: MEG recordings are decomposed by independent component analysis (ICA) into artifact and brain components (ICs); next, the channel maps associated with the latter ones are projected into the source space and the resulting voxel-wise weights are used to linearly combine the IC time courses. An extensive description of the proposed pipeline is provided here, along with an assessment of its performances with respect to alternative approaches. The following investigations were carried out: (1) ICA decomposition algorithm. Synthetic data are used to assess the sensitivity of the ICA results to the decomposition algorithm, by testing FastICA, INFOMAX, and SOBI. FastICA with deflation approach, a standard solution, provides the best decomposition. (2) Recombination of brain ICs versus subtraction of artifactual ICs (at the channel level). Both the recombination of the brain ICs in the sensor space and the classical procedure of subtracting the artifactual ICs from the recordings provide a suitable reconstruction, with a lower distortion using the latter approach. (3) Recombination of brain ICs after localization versus localization of artifact-corrected recordings. The brain IC recombination after source localization, as implemented in the proposed pipeline, provides a lower source-level signal distortion. (4) Detection of RSNs. The accuracy in source-level reconstruction by the proposed pipeline is confirmed by an improved specificity in the retrieval of RSNs from experimental data.

Fractal Dimension of EEG Activity Senses Neuronal Impairment in Acute Stroke

The brain is a self-organizing system which displays self-similarities at different spatial and temporal scales. Thus, the complexity of its dynamics, associated to efficient processing and functional advantages, is expected to be captured by a measure of its scale-free (fractal) properties. Under the hypothesis that the fractal dimension (FD) of the electroencephalographic signal (EEG) is optimally sensitive to the neuronal dysfunction secondary to a brain lesion, we tested the FD’s ability in assessing two key processes in acute stroke: the clinical impairment and the recovery prognosis. Resting EEG was collected in 36 patients 4–10 days after a unilateral ischemic stroke in the middle cerebral artery territory and 19 healthy controls. National Health Institute Stroke Scale (NIHss) was collected at T0 and 6 months later. Highuchi FD, its inter-hemispheric asymmetry (FDasy) and spectral band powers were calculated for EEG signals. FD was smaller in patients than in controls (1.447±0.092 vs 1.525±0.105) and its reduction was paired to a worse acute clinical status. FD decrease was associated to alpha increase and beta decrease of oscillatory activity power. Larger FDasy in acute phase was paired to a worse clinical recovery at six months. FD in our patients captured the loss of complexity reflecting the global system dysfunction resulting from the structural damage. This decrease seems to reveal the intimate nature of structure-function unity, where the regional neural multi-scale self-similar activity is impaired by the anatomical lesion. This picture is coherent with neuronal activity complexity decrease paired to a reduced repertoire of functional abilities. FDasy result highlights the functional relevance of the balance between homologous brain structures’ activities in stroke recovery.