Brainstorm-DUNEuro: An integrated and user-friendly Finite Element Method for modeling electromagnetic brain activity.

Human brain activity generates scalp potentials (electroencephalography - EEG), intracranial potentials (iEEG), and external magnetic fields (magnetoencephalography - MEG). These electrophysiology (e-phys) signals can often be measured simultaneously for research and clinical applications. The forward problem involves modeling these signals at their sensors for a given equivalent current dipole configuration within the brain. While earlier researchers modeled the head as a simple set of isotropic spheres, today's magnetic resonance imaging (MRI) data allow for a detailed anatomic description of brain structures and anisotropic characterization of tissue conductivities. We present a complete pipeline, integrated into the Brainstorm software, that allows users to automatically generate an individual and accurate head model based on the subject's MRI and calculate the electromagnetic forward solution using the finite element method (FEM). The head model generation is performed by integrating the latest tools for MRI segmentation and FEM mesh generation. The final head model comprises the five main compartments: white-matter, gray-matter, CSF, skull, and scalp. The anisotropic brain conductivity model is based on the effective medium approach (EMA), which estimates anisotropic conductivity tensors from diffusion-weighted imaging (DWI) data. The FEM electromagnetic forward solution is obtained through the DUNEuro library, integrated into Brainstorm, and accessible with either a user-friendly graphical interface or scripting. With tutorials and example data sets available in an open-source format on the Brainstorm website, this integrated pipeline provides access to advanced FEM tools for electromagnetic modeling to a broader neuroscience community.

Variability of Single Pulse Electrical Stimulation Responses Recorded with Intracranial Electroencephalography in Epileptic Patients.

Cohort studies of brain stimulations performed with stereo-electroencephalographic (SEEG) electrodes in epileptic patients allow to derive large scale functional connectivity. It is known, however, that brain responses to electrical or magnetic stimulation techniques are not always reproducible. Here, we study variability of responses to single pulse SEEG electrical stimulation. We introduce a second-order probability analysis, i.e. we extend estimation of connection probabilities, defined as the proportion of responses trespassing a statistical threshold (determined in terms of Z-score with respect to spontaneous neuronal activity before stimulation) over all responses and derived from a number of individual measurements, to an analysis of pairs of measurements.Data from 445 patients were processed. We found that variability between two equivalent measurements is substantial in particular conditions. For long ( > ~ 90 mm) distances between stimulating and recording sites, and threshold value Z = 3, correlation between measurements drops almost to zero. In general, it remains below 0.5 when the threshold is smaller than Z = 4 or the stimulating current intensity is 1 mA. It grows with an increase of either of these factors. Variability is independent of interictal spiking rates in the stimulating and recording sites.We conclude that responses to SEEG stimulation in the human brain are variable, i.e. in a subject at rest, two stimulation trains performed at the same electrode contacts and with the same protocol can give discrepant results. Our findings highlight an advantage of probabilistic interpretation of such results even in the context of a single individual.

Advances in human intracranial electroencephalography research, guidelines and good practices.

Since the second-half of the twentieth century, intracranial electroencephalography (iEEG), including both electrocorticography (ECoG) and stereo-electroencephalography (sEEG), has provided an intimate view into the human brain. At the interface between fundamental research and the clinic, iEEG provides both high temporal resolution and high spatial specificity but comes with constraints, such as the individual's tailored sparsity of electrode sampling. Over the years, researchers in neuroscience developed their practices to make the most of the iEEG approach. Here we offer a critical review of iEEG research practices in a didactic framework for newcomers, as well addressing issues encountered by proficient researchers. The scope is threefold: (i) review common practices in iEEG research, (ii) suggest potential guidelines for working with iEEG data and answer frequently asked questions based on the most widespread practices, and (iii) based on current neurophysiological knowledge and methodologies, pave the way to good practice standards in iEEG research. The organization of this paper follows the steps of iEEG data processing. The first section contextualizes iEEG data collection. The second section focuses on localization of intracranial electrodes. The third section highlights the main pre-processing steps. The fourth section presents iEEG signal analysis methods. The fifth section discusses statistical approaches. The sixth section draws some unique perspectives on iEEG research. Finally, to ensure a consistent nomenclature throughout the manuscript and to align with other guidelines, e.g., Brain Imaging Data Structure (BIDS) and the OHBM Committee on Best Practices in Data Analysis and Sharing (COBIDAS), we provide a glossary to disambiguate terms related to iEEG research.

A brain atlas of axonal and synaptic delays based on modelling of cortico-cortical evoked potentials.

Epilepsy presurgical investigation may include focal intracortical single-pulse electrical stimulations with depth electrodes, which induce cortico-cortical evoked potentials at distant sites because of white matter connectivity. Cortico-cortical evoked potentials provide a unique window on functional brain networks because they contain sufficient information to infer dynamical properties of large-scale brain connectivity, such as preferred directionality and propagation latencies. Here, we developed a biologically informed modelling approach to estimate the neural physiological parameters of brain functional networks from the cortico-cortical evoked potentials recorded in a large multicentric database. Specifically, we considered each cortico-cortical evoked potential as the output of a transient stimulus entering the stimulated region, which directly propagated to the recording region. Both regions were modelled as coupled neural mass models, the parameters of which were estimated from the first cortico-cortical evoked potential component, occurring before 80 ms, using dynamic causal modelling and Bayesian model inversion. This methodology was applied to the data of 780 patients with epilepsy from the F-TRACT database, providing a total of 34 354 bipolar stimulations and 774 445 cortico-cortical evoked potentials. The cortical mapping of the local excitatory and inhibitory synaptic time constants and of the axonal conduction delays between cortical regions was obtained at the population level using anatomy-based averaging procedures, based on the Lausanne2008 and the HCP-MMP1 parcellation schemes, containing 130 and 360 parcels, respectively. To rule out brain maturation effects, a separate analysis was performed for older (>15 years) and younger patients (<15 years). In the group of older subjects, we found that the cortico-cortical axonal conduction delays between parcels were globally short (median = 10.2 ms) and only 16% were larger than 20 ms. This was associated to a median velocity of 3.9 m/s. Although a general lengthening of these delays with the distance between the stimulating and recording contacts was observed across the cortex, some regions were less affected by this rule, such as the insula for which almost all efferent and afferent connections were faster than 10 ms. Synaptic time constants were found to be shorter in the sensorimotor, medial occipital and latero-temporal regions, than in other cortical areas. Finally, we found that axonal conduction delays were significantly larger in the group of subjects younger than 15 years, which corroborates that brain maturation increases the speed of brain dynamics. To our knowledge, this study is the first to provide a local estimation of axonal conduction delays and synaptic time constants across the whole human cortex in vivo, based on intracerebral electrophysiological recordings.

The Sensitivity of Ear-EEG: Evaluating the Source-Sensor Relationship Using Forward Modeling.

Ear-EEG allows to record brain activity in every-day life, for example to study natural behaviour or unhindered social interactions. Compared to conventional scalp-EEG, ear-EEG uses fewer electrodes and covers only a small part of the head. Consequently, ear-EEG will be less sensitive to some cortical sources. Here, we perform realistic electromagnetic simulations to compare cEEGrid ear-EEG with 128-channel cap-EEG. We compute the sensitivity of ear-EEG for different cortical sources, and quantify the expected signal loss of ear-EEG relative to cap-EEG. Our results show that ear-EEG is most sensitive to sources in the temporal cortex. Furthermore, we show how ear-EEG benefits from a multi-channel configuration (i.e. cEEGrid). The pipelines presented here can be adapted to any arrangement of electrodes and can therefore provide an estimate of sensitivity to cortical regions, thereby increasing the chance of successful experiments using ear-EEG.

Mapping the Insula with Stereo-Electroencephalography: The Emergence of Semiology in Insula Lobe Seizures.

OBJECTIVE: Insula epilepsy is rare and can be evaluated effectively by Stereotactic intracerebral EEG (SEEG). Many previous studies of insulo-opercular seizures have been unable to separate insular and opercular onset. With adequate sampling of the insula, this study shows this is possible. METHODS: We analyzed intrainsular dynamics and extrainsular propagation in 12 patients with "pure" insula epilepsy (n = 9) or insular and only deepest opercular involvement (n = 3) at seizure onset. Review of semiology defined clinical groups, agglomerative cluster, and principal component analysis of semiological features were performed. Quantitative epileptogenicity, and intrainsular and extrainsular propagation were computed via time frequency analysis and epileptogenicity mapping. RESULTS: Seizure onset patterns were heterogeneous; the seizure onset zone was focal. Seizure onset and first ictal change within insula functional subdivision correlated with aura and reflex component. No paninsular spread occurred; contralateral insular spread was very early. While the discharge was intrainsular, clinical signs related to aura or vegetative signs. Extrainsular propagation was early and related to the emergence of the majority of clinical signs. Cluster analysis found an anterior, intermediate, and posterior insula seizure onset group. The largest principal component separated anterior insula manifestations, including early hypermotor signs, early recovery, and no aura from posterior insula features of early dystonia, early tonic motor features, and sensorimotor aura. INTERPRETATION: Aura is vital to identifying seizure onset and relates to insula functional subdivision. Seizures are heterogenous; extrainsular propagation occurs early, accounting for most of the semiology. With adequate sampling, "pure" insula epilepsy can be identified and focal curative resection is possible. ANN NEUROL 2020;88:477-488.

Epileptogenicity Mapping: A Quantitative Approach to Identify the Seizure Onset.

During the presurgical evaluation of patients with focal refractory epilepsies, the spatial mapping of the seizure onset zone (SOZ) and seizure propagation networks critically depends on the use of different features extracted from the intracranial electroencephalogram (IEEG). The identification of the SOZ is usually based on visual inspection by highly qualified neurophysiologists. However, quantitative IEEG analyses have recently been developed by exploiting signal and image characteristics in order to improve and expedite the SOZ detection. Here, the authors briefly review some of the latest methods proposed by different research groups and then present the recent implementation in Brainstorm software.

Integrated open-source software for multiscale electrophysiology.

The methods for electrophysiology in neuroscience have evolved tremendously over the recent years with a growing emphasis on dense-array signal recordings. Such increased complexity and augmented wealth in the volume of data recorded, have not been accompanied by efforts to streamline and facilitate access to processing methods, which too are susceptible to grow in sophistication. Moreover, unsuccessful attempts to reproduce peer-reviewed publications indicate a problem of transparency in science. This growing problem could be tackled by unrestricted access to methods that promote research transparency and data sharing, ensuring the reproducibility of published results. Here, we provide a free, extensive, open-source software that provides data-analysis, data-management and multi-modality integration solutions for invasive neurophysiology. Users can perform their entire analysis through a user-friendly environment without the need of programming skills, in a tractable (logged) way. This work contributes to open-science, analysis standardization, transparency and reproducibility in invasive neurophysiology.

Brainstorm Pipeline Analysis of Resting-State Data From the Open MEG Archive.

We present a simple, reproducible analysis pipeline applied to resting-state magnetoencephalography (MEG) data from the Open MEG Archive (OMEGA). The data workflow was implemented with Brainstorm, which like OMEGA is free and openly accessible. The proposed pipeline produces group maps of ongoing brain activity decomposed in the typical frequency bands of electrophysiology. The procedure is presented as a technical proof of concept for streamlining a broader range and more sophisticated studies of resting-state electrophysiological data. It also features the recently introduced extension of the brain imaging data structure (BIDS) to MEG data, highlighting the scalability and generalizability of Brainstorm analytical pipelines to other, and potentially larger data volumes.

MEG/EEG Group Analysis With Brainstorm.

Brainstorm is a free, open-source Matlab and Java application for multimodal electrophysiology data analytics and source imaging [primarily MEG, EEG and depth recordings, and integration with MRI and functional near infrared spectroscopy (fNIRS)]. We also provide a free, platform-independent executable version to users without a commercial Matlab license. Brainstorm has a rich and intuitive graphical user interface, which facilitates learning and augments productivity for a wider range of neuroscience users with little or no knowledge of scientific coding and scripting. Yet, it can also be used as a powerful scripting tool for reproducible and shareable batch processing of (large) data volumes. This article describes these Brainstorm interactive and scripted features via illustration through the complete analysis of group data from 16 participants in a MEG vision study.

Probabilistic functional tractography of the human cortex revisited.

In patients with pharmaco-resistant focal epilepsies investigated with intracranial electroencephalography (iEEG), direct electrical stimulations of a cortical region induce cortico-cortical evoked potentials (CCEP) in distant cerebral cortex, which properties can be used to infer large scale brain connectivity. In 2013, we proposed a new probabilistic functional tractography methodology to study human brain connectivity. We have now been revisiting this method in the F-TRACT project (f-tract.eu) by developing a large multicenter CCEP database of several thousand stimulation runs performed in several hundred patients, and associated processing tools to create a probabilistic atlas of human cortico-cortical connections. Here, we wish to present a snapshot of the methods and data of F-TRACT using a pool of 213 epilepsy patients, all studied by stereo-encephalography with intracerebral depth electrodes. The CCEPs were processed using an automated pipeline with the following consecutive steps: detection of each stimulation run from stimulation artifacts in raw intracranial EEG (iEEG) files, bad channels detection with a machine learning approach, model-based stimulation artifact correction, robust averaging over stimulation pulses. Effective connectivity between the stimulated and recording areas is then inferred from the properties of the first CCEP component, i.e. onset and peak latency, amplitude, duration and integral of the significant part. Finally, group statistics of CCEP features are implemented for each brain parcel explored by iEEG electrodes. The localization (coordinates, white/gray matter relative positioning) of electrode contacts were obtained from imaging data (anatomical MRI or CT scans before and after electrodes implantation). The iEEG contacts were repositioned in different brain parcellations from the segmentation of patients' anatomical MRI or from templates in the MNI coordinate system. The F-TRACT database using the first pool of 213 patients provided connectivity probability values for 95% of possible intrahemispheric and 56% of interhemispheric connections and CCEP features for 78% of intrahemisheric and 14% of interhemispheric connections. In this report, we show some examples of anatomo-functional connectivity matrices, and associated directional maps. We also indicate how CCEP features, especially latencies, are related to spatial distances, and allow estimating the velocity distribution of neuronal signals at a large scale. Finally, we describe the impact on the estimated connectivity of the stimulation charge and of the contact localization according to the white or gray matter. The most relevant maps for the scientific community are available for download on f-tract. eu (David et al., 2017) and will be regularly updated during the following months with the addition of more data in the F-TRACT database. This will provide an unprecedented knowledge on the dynamical properties of large fiber tracts in human.

IntrAnat Electrodes: A Free Database and Visualization Software for Intracranial Electroencephalographic Data Processed for Case and Group Studies.

In some cases of pharmaco-resistant and focal epilepsies, intracranial recordings performed epidurally (electrocorticography, ECoG) and/or in depth (stereoelectroencephalography, SEEG) can be required to locate the seizure onset zone and the eloquent cortex before surgical resection. In SEEG, each electrode contact records brain's electrical activity in a spherical volume of 3 mm diameter approximately. The spatial coverage is around 1% of the brain and differs between patients because the implantation of electrodes is tailored for each case. Group studies thus need a large number of patients to reach a large spatial sampling, which can be achieved more easily using a multicentric approach such as implemented in our F-TRACT project (f-tract.eu). To facilitate group studies, we developed a software-IntrAnat Electrodes-that allows to perform virtual electrode implantation in patients' neuroanatomy and to overlay results of epileptic and functional mapping, as well as resection masks from the surgery. IntrAnat Electrodes is based on a patient database providing multiple search criteria to highlight various group features. For each patient, the anatomical processing is based on a series of software publicly available. Imaging modalities (Positron Emission Tomography (PET), anatomical MRI pre-implantation, post-implantation and post-resection, functional MRI, diffusion MRI, Computed Tomography (CT) with electrodes) are coregistered. The 3D T1 pre-implantation MRI gray/white matter is segmented and spatially normalized to obtain a series of cortical parcels using different neuroanatomical atlases. On post-implantation images, the user can position 3D models of electrodes defined by their geometry. Each electrode contact is then labeled according to its position in the anatomical atlases, to the class of tissue (gray or white matter, cerebro-spinal fluid) and to its presence inside or outside the resection mask. Users can add more functionally informed labels on contact, such as clinical responses after electrical stimulation, cortico-cortical evoked potentials, gamma band activity during cognitive tasks or epileptogenicity. IntrAnat Electrodes software thus provides a means to visualize multimodal data. The contact labels allow to search for patients in the database according to multiple criteria representing almost all available data, which is to our knowledge unique in current SEEG software. IntrAnat Electrodes will be available in the forthcoming release of BrainVisa software and tutorials can be found on the F-TRACT webpage.

MEG-BIDS, the brain imaging data structure extended to magnetoencephalography.

We present a significant extension of the Brain Imaging Data Structure (BIDS) to support the specific aspects of magnetoencephalography (MEG) data. MEG measures brain activity with millisecond temporal resolution and unique source imaging capabilities. So far, BIDS was a solution to organise magnetic resonance imaging (MRI) data. The nature and acquisition parameters of MRI and MEG data are strongly dissimilar. Although there is no standard data format for MEG, we propose MEG-BIDS as a principled solution to store, organise, process and share the multidimensional data volumes produced by the modality. The standard also includes well-defined metadata, to facilitate future data harmonisation and sharing efforts. This responds to unmet needs from the multimodal neuroimaging community and paves the way to further integration of other techniques in electrophysiology. MEG-BIDS builds on MRI-BIDS, extending BIDS to a multimodal data structure. We feature several data-analytics software that have adopted MEG-BIDS, and a diverse sample of open MEG-BIDS data resources available to everyone.

OMEGA: The Open MEG Archive.

In contrast with other imaging modalities, there is presently a scarcity of fully open resources in magnetoencephalography (MEG) available to the neuroimaging community. Here we present a collaborative effort led by the McConnell Brain Imaging Centre of the Montreal Neurological Institute, and the Université de Montréal to build and share a centralised repository to curate MEG data in raw and processed form for open dissemination. The Open MEG Archive (OMEGA, omega.bic.mni.mcgill.ca) is bound to become a continuously expanding repository of multimodal data with a primary focus on MEG, in addition to storing anatomical MRI volumes, demographic participant data and questionnaires, and other forms of electrophysiological data such as EEG. The OMEGA initiative offers both the technological framework for multi-site MEG data aggregation, and serves as one of the largest freely available resting-state and eventually task-related MEG datasets presently available.

Brainstorm: a user-friendly application for MEG/EEG analysis.

Brainstorm is a collaborative open-source application dedicated to magnetoencephalography (MEG) and electroencephalography (EEG) data visualization and processing, with an emphasis on cortical source estimation techniques and their integration with anatomical magnetic resonance imaging (MRI) data. The primary objective of the software is to connect MEG/EEG neuroscience investigators with both the best-established and cutting-edge methods through a simple and intuitive graphical user interface (GUI).

Neural signature of the conscious processing of auditory regularities.

Can conscious processing be inferred from neurophysiological measurements? Some models stipulate that the active maintenance of perceptual representations across time requires consciousness. Capitalizing on this assumption, we designed an auditory paradigm that evaluates cerebral responses to violations of temporal regularities that are either local in time or global across several seconds. Local violations led to an early response in auditory cortex, independent of attention or the presence of a concurrent visual task, whereas global violations led to a late and spatially distributed response that was only present when subjects were attentive and aware of the violations. We could detect the global effect in individual subjects using functional MRI and both scalp and intracerebral event-related potentials. Recordings from 8 noncommunicating patients with disorders of consciousness confirmed that only conscious individuals presented a global effect. Taken together these observations suggest that the presence of the global effect is a signature of conscious processing, although it can be absent in conscious subjects who are not aware of the global auditory regularities. This simple electrophysiological marker could thus serve as a useful clinical tool.