White Matter Information Flow Mapping from Diffusion MRI and EEG

International audienceThe human brain can be described as a network of specialized and spatially distributed regions. The activity of individual regions can be estimated using electroencephalography and the structure of the network can be measured using diffusion magnetic resonance imaging. However, the communication between the different cortical regions occurring through the white matter, coined information flow, cannot be observed by either modalities independently. Here, we present a new method to infer information flow in the white matter of the brain from joint diffusion MRI and EEG measurements. This is made possible by the millisecond resolution of EEG which makes the transfer of information from one region to another observable. A subject specific Bayesian network is built which captures the possible interactions between brain regions at different times. This network encodes the connections between brain regions detected using diffusion MRI tractography derived white matter bundles and their associated delays. By injecting the EEG measurements as evidence into this model, we are able to estimate the directed dynamical functional connectivity whose delays are supported by the diffusion MRI derived structural connectivity. We present our results in the form of information flow diagrams that trace transient communication between cortical regions over a functional data window. The performance of our algorithm under different noise levels is assessed using receiver operating characteristic curves on simulated data. In addition, using the well-characterized visual motor network as grounds to test our model, we present the information flow obtained during a reaching task following left or right visual stimuli. These promising results present the transfer of information from the eyes to the primary motor cortex. The information flow obtained using our technique can also be projected back to the anatomy and animated to produce videos of the information path through the white matter, opening a new window into multi-modal dynamic brain connectivity

Tracking dynamic interactions between structural and functional connectivity : a TMS/EEG-dMRI study

Transcranial magnetic stimulation (TMS) in combination with neuroimaging techniques allows to measure the effects of a direct perturbation of the brain. When coupled with high-density electroencephalography (TMS/hd-EEG), TMS pulses revealed electrophysiological signatures of different cortical modules in health and disease. However, the neural underpinnings of these signatures remain unclear. Here, by applying multimodal analyses of cortical response to TMS recordings and diffusion magnetic resonance imaging (dMRI) tractography, we investigated the relationship between functional and structural features of different cortical modules in a cohort of awake healthy volunteers. For each subject, we computed directed functional connectivity interactions between cortical areas from the source-reconstructed TMS/hd-EEG recordings and correlated them with the correspondent structural connectivity matrix extracted from dMRI tractography, in three different frequency bands (alpha, beta, gamma) and two sites of stimulation (left precuneus and left premotor). Each stimulated area appeared to mainly respond to TMS by being functionally elicited in specific frequency bands, that is, beta for precuneus and gamma for premotor. We also observed a temporary decrease in the whole-brain correlation between directed functional connectivity and structural connectivity after TMS in all frequency bands. Notably, when focusing on the stimulated areas only, we found that the structure-function correlation significantly increases over time in the premotor area controlateral to TMS. Our study points out the importance of taking into account the major role played by different cortical oscillations when investigating the mechanisms for integration and segregation of information in the human brain

Genetic and Neuroanatomical Support for Functional Brain Network Dynamics in Epilepsy

Focal epilepsy is a devastating neurological disorder that affects an overwhelming number of patients worldwide, many of whom prove resistant to medication. The efficacy of current innovative technologies for the treatment of these patients has been stalled by the lack of accurate and effective methods to fuse multimodal neuroimaging data to map anatomical targets driving seizure dynamics. Here we propose a parsimonious model that explains how large-scale anatomical networks and shared genetic constraints shape inter-regional communication in focal epilepsy. In extensive ECoG recordings acquired from a group of patients with medically refractory focal-onset epilepsy, we find that ictal and preictal functional brain network dynamics can be accurately predicted from features of brain anatomy and geometry, patterns of white matter connectivity, and constraints complicit in patterns of gene coexpression, all of which are conserved across healthy adult populations. Moreover, we uncover evidence that markers of non-conserved architecture, potentially driven by idiosyncratic pathology of single subjects, are most prevalent in high frequency ictal dynamics and low frequency preictal dynamics. Finally, we find that ictal dynamics are better predicted by white matter features and more poorly predicted by geometry and genetic constraints than preictal dynamics, suggesting that the functional brain network dynamics manifest in seizures rely on - and may directly propagate along - underlying white matter structure that is largely conserved across humans. Broadly, our work offers insights into the generic architectural principles of the human brain that impact seizure dynamics, and could be extended to further our understanding, models, and predictions of subject-level pathology and response to intervention

Getting the best outcomes from epilepsy surgery.

Neurosurgery is an underutilized treatment that can potentially cure drug-refractory epilepsy. Careful, multidisciplinary presurgical evaluation is vital for selecting patients and to ensure optimal outcomes. Advances in neuroimaging have improved diagnosis and guided surgical intervention. Invasive electroencephalography allows the evaluation of complex patients who would otherwise not be candidates for neurosurgery. We review the current state of the assessment and selection of patients and consider established and novel surgical procedures and associated outcome data. We aim to dispel myths that may inhibit physicians from referring and patients from considering neurosurgical intervention for drug-refractory focal epilepsies. Ann Neurol 2018;83:676-690

Brain enhancement through cognitive training: A new insight from brain connectome

Owing to the recent advances in neurotechnology and the progress in understanding of brain cognitive functions, improvements of cognitive performance or acceleration of learning process with brain enhancement systems is not out of our reach anymore, on the contrary, it is a tangible target of contemporary research. Although a variety of approaches have been proposed, we will mainly focus on cognitive training interventions, in which learners repeatedly perform cognitive tasks to improve their cognitive abilities. In this review article, we propose that the learning process during the cognitive training can be facilitated by an assistive system monitoring cognitive workloads using electroencephalography (EEG) biomarkers, and the brain connectome approach can provide additional valuable biomarkers for facilitating leaners' learning processes. For the purpose, we will introduce studies on the cognitive training interventions, EEG biomarkers for cognitive workload, and human brain connectome. As cognitive overload and mental fatigue would reduce or even eliminate gains of cognitive training interventions, a real-time monitoring of cognitive workload can facilitate the learning process by flexibly adjusting difficulty levels of the training task. Moreover, cognitive training interventions should have effects on brain sub-networks, not on a single brain region, and graph theoretical network metrics quantifying topological architecture of the brain network can differentiate with respect to individual cognitive states as well as to different individuals' cognitive abilities, suggesting that the connectome is a valuable approach for tracking the learning progress. Although only a few studies have exploited the connectome approach for studying alterations of the brain network induced by cognitive training interventions so far, we believe that it would be a useful technique for capturing improvements of cognitive function

The ENIGMA Stroke Recovery Working Group: Big data neuroimaging to study brain–behavior relationships after stroke

The goal of the Enhancing Neuroimaging Genetics through Meta‐Analysis (ENIGMA) Stroke Recovery working group is to understand brain and behavior relationships using well‐powered meta‐ and mega‐analytic approaches. ENIGMA Stroke Recovery has data from over 2,100 stroke patients collected across 39 research studies and 10 countries around the world, comprising the largest multisite retrospective stroke data collaboration to date. This article outlines the efforts taken by the ENIGMA Stroke Recovery working group to develop neuroinformatics protocols and methods to manage multisite stroke brain magnetic resonance imaging, behavioral and demographics data. Specifically, the processes for scalable data intake and preprocessing, multisite data harmonization, and large‐scale stroke lesion analysis are described, and challenges unique to this type of big data collaboration in stroke research are discussed. Finally, future directions and limitations, as well as recommendations for improved data harmonization through prospective data collection and data management, are provided

Incorporation of anisotropic conductivities in EEG source analysis

The electroencephalogram (EEG) is a measurement of brain activity over a period of time by placing electrodes at the scalp (surface EEG) or in the brain (depth EEG) and is used extensively in the clinical practice. In the past 20 years, EEG source analysis has been increasingly used as a tool in the diagnosis of neurological disorders (like epilepsy) and in the research of brain functionality. EEG source analysis estimates the origin of brain activity given the electrode potentials measured at the scalp. This involves solving an inverse problem where a forward solution, which depends on the source parameters, is fitted to the given set of electrode potentials. The forward solution are the electrode potentials caused by a source in a given head model. The head model is dependent on the geometry and the conductivity. Often an isotropic conductivity (i.e. the conductivity is equal in all directions) is used, although the skull and white matter have an anisotropic conductivity (i.e. the conductivity can differ depending on the direction the current flows). In this dissertation a way to incorporate the anisotropic conductivities is presented and the effect of not incorporating these anisotropic conductivities is investigated. Spherical head models are simple head models where an analytical solution to the forward problem exists. A small simulation study in a 5 shell spherical head model was performed to investigate the estimation error due to neglecting the anisotropic properties of skull and white matter. The results show that the errors in the dipole location can be larger than 15 mm, which is unacceptable for an accurate dipole estimation in the clinical practice. Therefore, anisotropic conductivities have to be included in the head model. However, these spherical head models are not representative for the human head. Realistic head models are usually made from magnetic resonance scans through segmentation and are a better approximation to the geometry of the human head. To solve the forward problem in these head models numerical methods are needed. In this dissertation we proposed a finite difference technique that can incorporate anisotropic conductivities. Moreover, by using the reciprocity theorem the forward calculation time during an dipole source estimation procedure can be significantly reduced. By comparing the analytical solution for the dipole estimation problem with the one using the numerical method, the anisotropic finite difference with reciprocity method (AFDRM) is validated. Therefore, a cubic grid is made on the 5 shell spherical head model. The electrode potentials are obtained in the spherical head model with anisotropic conductivities by solving the forward problem using the analytical solution. Using these electrode potentials the inverse problem was solved in the spherical head model using the AFDRM. In this way we can determine the location error due to using the numerical technique. We found that the incorporation of anisotropic conductivities results in a larger location error when the head models are fully isotropical conducting. Furthermore, the location error due to the numerical technique is smaller if the cubic grid is made finer. To minimize the errors due to the numerical technique, the cubic grid should be smaller than or equal to 1 mm. Once the numerical technique is validated, a realistic head model can now be constructed. As a cubic grid should be used of at most 1 mm, the use of segmented T1 magnetic resonance images is best suited the construction. The anisotropic conductivities of skull and white matter are added as follows: The anisotropic conductivity of the skull is derived by calculating the normal and tangential direction to the skull at each voxel. The conductivity in the tangential direction was set 10 times larger than the normal direction. The conductivity of the white matter was derived using diffusion weighted magnetic resonance imaging (DW-MRI), a technique that measures the diffusion of water in several directions. As diffusion is larger along the nerve fibers, it is assumed that the conductivity along the nerve fibers is larger than the perpendicular directions to the nerve bundle. From the diffusion along each direction, the conductivity can be derived using two approaches. A simplified approach takes the direction with the largest diffusion and sets the conductivity along that direction 9 times larger than the orthogonal direction. However, by calculating the fractional anisotropy, a well-known measure indicating the degree of anisotropy, we can appreciate that a fractional anisotropy of 0.8715 is an overestimation. In reality, the fractional anisotorpy is mostly smaller and variable throughout the white matter. A realistic approach was therefore presented, which states that the conductivity tensor is a scaling of the diffusion tensor. The volume constraint is used to determine the scaling factor. A comparison between the realistic approach and the simplified approach was made. The results showed that the location error was on average 4.0 mm with a maximum of 10 mm. The orientation error was found that the orientation could range up to 60 degrees. The large orientation error was located at regions where the anisotropic ratio was low using the realistic approach but was 9 using the simplified approach. Furthermore, as the DW-MRI can also be used to measure the anisotropic diffusion in a gray matter voxel, we can derive a conductivity tensor. After investigating the errors due to neglecting these anisotropic conductivities of the gray matter, we found that the location error was very small (average dipole location error: 2.8 mm). The orientation error was ranged up to 40 degrees, although the mean was 5.0 degrees. The large errors were mostly found at the regions that had a high anisotropic ratio in the anisotropic conducting gray matter. Mostly these effects were due to missegmentation or to partial volume effects near the boundary interfaces of the gray and white matter compartment. After the incorporation of the anisotropic conductivities in the realistic head model, simulation studies can be performed to investigate the dipole estimation errors when these anisotropic conductivities of the skull and brain tissues are not taken into account. This can be done by comparing the solution to the dipole estimation problem in a head model with anisotropic conductivities with the one in a head model, where all compartments are isotropic conducting. This way we determine the error when a simplified head model is used instead of a more realistic one. When the anisotropic conductivity of both the skull and white matter or the skull only was neglected, it was found that the location error between the original and the estimated dipole was on average, 10 mm (maximum: 25 mm). When the anisotropic conductivity of the brain tissue was neglected, the location error was much smaller (an average location error of 1.1 mm). It was found that the anisotropy of the skull acts as an extra shielding of the electrical activity as opposed to an isotropic skull. Moreover, we saw that if the dipole is close to a highly anisotropic region, the potential field is changed reasonable in the near vicinity of the location of the dipole. In reality EEG contains noise contributions. These noise contribution will interact with the systematical error by neglecting anisotropic conductivities. The question we wanted to solve was “Is it worthwhile to incorporate anisotropic conductivities, even if the EEG contains noise?” and “How much noise should the EEG contain so that incorporating anisotropic conductivities improves the accuracy of EEG source analysis?”. When considering the anisotropic conductivities of the skull and brain tissues and the skull only, the location error due to the noise and neglecting the anisotropic conductivities is larger then the location error due to noise only. When only neglecting the anisotropic conductivities of the brain tissues only, the location error due to noise is similar to the location error due to noise and neglecting the anisotropic conductivities. When more advanced MR techniques can be used a better model to construct the anisotropic conductivities of the soft brain tissues can be used, which could result in larger errors even in the presence of noise. However, this is subject to further investigation. This suggests that the anisotropic conductivities of the skull should be incorporated. The technique presented in the dissertation can be used to epileptic patients in the presurgical evaluation. In this procedure patients are evaluated by means of medical investigations to determine the cause of the epileptic seizures. Afterwards, a surgical procedure can be performed to render the patient seizure free. A data set from a patiënt was obtained from a database of the Reference Center of Refractory Epilepsy of the Department of Neurology and the Department of Radiology of the Ghent University Hospital (Ghent, Belgium). The patient was monitored with a video/EEG monitoring with scalp and with implanted depth electrodes. An MR image was taken from the patient with the implanted depth electrodes, therefore, we could pinpoint the hippocampus as the onset zone of the epileptic seizures. The patient underwent a resective surgery removing the hippocampus, which rendered the patient seizure free. As DW-MRI images were not available, the head model constructed in chapter 4 and 5 was used. A neuroradiologist aligned the hippocampus in the MR image from which the head model was constructed. A spike was picked from a dataset and was used to estimate the source in a head model where all compartments were isotropic conducting, on one hand, and where the skull and brain tissues were anisotropic conducting, on the other. It was found that using the anisotropic head model, the source was estimated closer to the segmented hippocampus than the isotropic head model. This example shows the possibilities of this technique and allows us to apply it in the clinical practice. Moreover, a thorough validation of the technique has yet to be performed. There is a lot of discussion in the clinical community whether the spikes and epileptical seizures originate from the same origin in the brain. This question can be solved by applying our technique in patient studies

Loss of brain inter-frequency hubs in Alzheimer's disease

Alzheimer's disease (AD) causes alterations of brain network structure and function. The latter consists of connectivity changes between oscillatory processes at different frequency channels. We proposed a multi-layer network approach to analyze multiple-frequency brain networks inferred from magnetoencephalographic recordings during resting-states in AD subjects and age-matched controls. Main results showed that brain networks tend to facilitate information propagation across different frequencies, as measured by the multi-participation coefficient (MPC). However, regional connectivity in AD subjects was abnormally distributed across frequency bands as compared to controls, causing significant decreases of MPC. This effect was mainly localized in association areas and in the cingulate cortex, which acted, in the healthy group, as a true inter-frequency hub. MPC values significantly correlated with memory impairment of AD subjects, as measured by the total recall score. Most predictive regions belonged to components of the default-mode network that are typically affected by atrophy, metabolism disruption and amyloid-beta deposition. We evaluated the diagnostic power of the MPC and we showed that it led to increased classification accuracy (78.39%) and sensitivity (91.11%). These findings shed new light on the brain functional alterations underlying AD and provide analytical tools for identifying multi-frequency neural mechanisms of brain diseases.Comment: 27 pages, 6 figures, 3 tables, 3 supplementary figure

Evaluation of rigid registration methods for whole head imaging in diffuse optical tomography

Functional brain imaging has become an important neuroimaging technique for the study of brain organization and development. Compared to other imaging techniques, diffuse optical tomography (DOT) is a portable and low-cost technique that can be applied to infants and hospitalized patients using an atlas-based light model. For DOT imaging, the accuracy of the forward model has a direct effect on the resulting recovered brain function within a field of view and so the accuracy of the spatially normalized atlas-based forward models must be evaluated. Herein, the accuracy of atlas-based DOT is evaluated on models that are spatially normalized via a number of different rigid registration methods on 24 subjects. A multileveled approach is developed to evaluate the correlation of the geometrical and sensitivity accuracies across the full field of view as well as within specific functional subregions. Results demonstrate that different registration methods are optimal for recovery of different sets of functional brain regions. However, the “nearest point to point” registration method, based on the EEG 19 landmark system, is shown to be the most appropriate registration method for image quality throughout the field of view of the high-density cap that covers the whole of the optically accessible cortex