Dipole Source Localization of Mouse Electroencephalogram Using the Fieldtrip Toolbox

The mouse model is an important research tool in neurosciences to examine brain function and diseases with genetic perturbation in different brain regions. However, the limited techniques to map activated brain regions under specific experimental manipulations has been a drawback of the mouse model compared to human functional brain mapping. Here, we present a functional brain mapping method for fast and robust in vivo brain mapping of the mouse brain. The method is based on the acquisition of high density electroencephalography (EEG) with a microarray and EEG source estimation to localize the electrophysiological origins. We adapted the Fieldtrip toolbox for the source estimation, taking advantage of its software openness and flexibility in modeling the EEG volume conduction. Three source estimation techniques were compared: Distribution source modeling with minimum-norm estimation (MNE), scanning with multiple signal classification (MUSIC), and single-dipole fitting. Known sources to evaluate the performance of the localization methods were provided using optogenetic tools. The accuracy was quantified based on the receiver operating characteristic (ROC) analysis. The mean detection accuracy was high, with a false positive rate less than 1.3% and 7% at the sensitivity of 90% plotted with the MNE and MUSIC algorithms, respectively. The mean center-to-center distance was less than 1.2 mm in single dipole fitting algorithm. Mouse microarray EEG source localization using microarray allows a reliable method for functional brain mapping in awake mouse opening an access to cross-species study with human brain.     

EEGNET: An Open Source Tool for Analyzing and Visualizing M/EEG Connectome

The brain is a large-scale complex network often referred to as the “connectome”. Exploring the dynamic behavior of the connectome is a challenging issue as both excellent time and space resolution is required. In this context Magneto/Electroencephalography (M/EEG) are effective neuroimaging techniques allowing for analysis of the dynamics of functional brain networks at scalp level and/or at reconstructed sources. However, a tool that can cover all the processing steps of identifying brain networks from M/EEG data is still missing. In this paper, we report a novel software package, called EEGNET, running under MATLAB (Math works, inc), and allowing for analysis and visualization of functional brain networks from M/EEG recordings. EEGNET is developed to analyze networks either at the level of scalp electrodes or at the level of reconstructed cortical sources. It includes i) Basic steps in preprocessing M/EEG signals, ii) the solution of the inverse problem to localize / reconstruct the cortical sources, iii) the computation of functional connectivity among signals collected at surface electrodes or/and time courses of reconstructed sources and iv) the computation of the network measures based on graph theory analysis. EEGNET is the unique tool that combines the M/EEG functional connectivity analysis and the computation of network measures derived from the graph theory. The first version of EEGNET is easy to use, flexible and user friendly. EEGNET is an open source tool and can be freely downloaded from this webpage: https://sites.google.com/site/eegnetworks/.     

Modeling the effect of sleep regulation on a neural mass model

In mammals, sleep is categorized by two main sleep stages, rapid eye movement (REM) and non-REM (NREM) sleep that are known to fulfill different functional roles, the most notable being the consolidation of memory. While REM sleep is characterized by brain activity similar to wakefulness, the EEG activity changes drastically with the emergence of K-complexes, sleep spindles and slow oscillations during NREM sleep. These changes are regulated by circadian and ultradian rhythms, which emerge from an intricate interplay between multiple neuronal populations in the brainstem, forebrain and hypothalamus and the resulting varying levels of neuromodulators. Recently, there has been progress in the understanding of those rhythms both from a physiological as well as theoretical perspective. However, how these neuromodulators affect the generation of the different EEG patterns and their temporal dynamics is poorly understood. Here, we build upon previous work on a neural mass model of the sleeping cortex and investigate the effect of those neuromodulators on the dynamics of the cortex and the corresponding transition between wakefulness and the different sleep stages. We show that our simplified model is sufficient to generate the essential features of human EEG over a full day. This approach builds a bridge between sleep regulatory networks and EEG generating neural mass models and provides a valuable tool for model validation.     

Possible microwave mechanisms of the mammalian nervous system

The hypothesis is examined that the living mammal generates and uses electromagnetic waves in the lower microwave frequency region as an integral part of the functioning of central and peripheral nervous systems. Analysis of the potential energy of a protein integral to the neural membrane compared to that of an extracellular positive ion yields major known features of action potential generation, and identification of the integral protein as a microwave emitter and absorber by changes in rotational energy state. Prolate spheroidal analysis of the adult human brain/skull cavity shows capability to support modes in the range 200 MHz to 3 GHz; spatial mode properties correspond to gross anatomy and neuromorphology of the brain. Phase-lock loop interaction between lower microwave modes and action potential conduction results in pulse microwave/action potential generation observable by EEG instrumentation as brain waves; alpha waves occur if the corpus callosum is the major neural tract involved. Spatially consistent Lorentz forces of standing microwaves provide physical basis for correspondence of spatial properties of microwave modes with brain anatomy, and for the formation of myelin sheath and the nodes of Ranvier on larger neurons.     

A latent state space model for estimating brain dynamics from electroencephalogram (EEG) data

Modern neuroimaging technologies have substantially advanced the measurement of brain activity. Electroencephalogram (EEG) as a noninvasive neuroimaging technique measures changes in electrical voltage on the scalp induced by brain cortical activity. With its high temporal resolution, EEG has emerged as an increasingly useful tool to study brain connectivity. Challenges with modeling EEG signals of complex brain activity include interactions among unknown sources, low signal-to-noise ratio, and substantial between-subject heterogeneity. In this work, we propose a state space model that jointly analyzes multichannel EEG signals and learns dynamics of different sources corresponding to brain cortical activity. Our model borrows strength from spatially correlated measurements and uses low-dimensional latent states to explain all observed channels. The model can account for patient heterogeneity and quantify the effect of a subject's covariates on the latent space. The EM algorithm, Kalman filtering, and bootstrap resampling are used to fit the state space model and provide comparisons between patient diagnostic groups. We apply the developed approach to a case-control study of alcoholism and reveal significant attenuation of brain activity in response to visual stimuli in alcoholic subjects compared to healthy controls.     

Des sources neuronales aux signaux EEG : modèle et interprétation des activités épileptiques

This work deals with the interpretation of electrophysiological patients recorded in epileptic patients candidate to surgery. This issue is addressed through a physiologically relevant model for the generation of scalp and intracerebral electroencephalographic (EEG) signals. The proposed model is based on a spatiotemporal representation of the sources of brain activity, which combines a distributed dipole source model and a model of coupled neuronal populations. Signals recorded by sensors (scalp and intracerebral) are then computed by solving the forward problem in the head volume conductor. In this paper, the EEG generation model is used to study the influence of some source-related parameters (spatial extent, position, synchronization) on simulated signals, during epileptic transient activity (interictal spikes). Results show that the model allows for studying, on the one hand, the relationship between the spatiotemporal organization of neuronal sources and the properties of the observed signals and, on the other hand, the relationship between surface and depth EEG signals.     

Dynamics of spontaneous activity in random networks with multiple neuron subtypes and synaptic noise

During slow-wave sleep, brain electrical activity is dominated by the slow (< 1 Hz) electroencephalogram (EEG) oscillations characterized by the periodic transitions between active (or Up) and silent (or Down) states in the membrane voltage of the cortical and thalamic neurons. Sleep slow oscillation is believed to play critical role in consolidation of recent memories. Past computational studies, based on the Hodgkin-Huxley type neuronal models, revealed possible intracellular and network mechanisms of the neuronal activity during sleep, however, they failed to explore the large-scale cortical network dynamics depending on collective behavior in the large populations of neurons. In this new study, we developed a novel class of reduced discrete time spiking neuron models for large-scale network simulations of wake and sleep dynamics. In addition to the spiking mechanism, the new model implemented nonlinearities capturing effects of the leak current, the Ca²⁺ dependent K⁺ current and the persistent Na⁺ current that were found to be critical for transitions between Up and Down states of the slow oscillation. We applied the new model to study large-scale two-dimensional cortical network activity during slow-wave sleep. Our study explained traveling wave dynamics and characteristic synchronization properties of transitions between Up and Down states of the slow oscillation as observed in vivo in recordings from cats. We further predict a critical role of synaptic noise and slow adaptive currents for spike sequence replay as found during sleep related memory consolidation.     

Electroencephalographic imaging of higher brain function.

High temporal resolution is necessary to resolve the rapidly changing patterns of brain activity that underlie mental function. Electroencephalography (EEG) provides temporal resolution in the millisecond range. However, traditional EEG technology and practice provide insufficient spatial detail to identify relationships between brain electrical events and structures and functions visualized by magnetic resonance imaging or positron emission tomography. Recent advances help to overcome this problem by recording EEGs from more electrodes, by registering EEG data with anatomical images, and by correcting the distortion caused by volume conduction of EEG signals through the skull and scalp. In addition, statistical measurements of sub-second interdependences between EEG time-series recorded from different locations can help to generate hypotheses about the instantaneous functional networks that form between different cortical regions during perception, thought and action. Example applications are presented from studies of language, attention and working memory. Along with its unique ability to monitor brain function as people perform everyday activities in the real world, these advances make modern EEG an invaluable complement to other functional neuroimaging modalities.     

Dynamic Imaging of Coherent Sources Reveals Different Network Connectivity Underlying the Generation and Perpetuation of Epileptic Seizures

The concept of focal epilepsies includes a seizure origin in brain regions with hyper synchronous activity (epileptogenic zone and seizure onset zone) and a complex epileptic network of different brain areas involved in the generation, propagation, and modulation of seizures. The purpose of this work was to study functional and effective connectivity between regions involved in networks of epileptic seizures. The beginning and middle part of focal seizures from ictal surface EEG data were analyzed using dynamic imaging of coherent sources (DICS), an inverse solution in the frequency domain which describes neuronal networks and coherences of oscillatory brain activities. The information flow (effective connectivity) between coherent sources was investigated using the renormalized partial directed coherence (RPDC) method. In 8/11 patients, the first and second source of epileptic activity as found by DICS were concordant with the operative resection site; these patients became seizure free after epilepsy surgery. In the remaining 3 patients, the results of DICS / RPDC calculations and the resection site were discordant; these patients had a poorer post-operative outcome. The first sources as found by DICS were located predominantly in cortical structures; subsequent sources included some subcortical structures: thalamus, Nucl. Subthalamicus and cerebellum. DICS seems to be a powerful tool to define the seizure onset zone and the epileptic networks involved. Seizure generation seems to be related to the propagation of epileptic activity from the primary source in the seizure onset zone, and maintenance of seizures is attributed to the perpetuation of epileptic activity between nodes in the epileptic network. Despite of these promising results, this proof of principle study needs further confirmation prior to the use of the described methods in the clinical praxis.     

Comparative power spectral analysis of simultaneous elecroencephalographic and magnetoencephalographic recordings in humans suggests non-resistive extracellular media

The resistive or non-resistive nature of the extracellular space in the brain is still debated, and is an important issue for correctly modeling extracellular potentials. Here, we first show theoretically that if the medium is resistive, the frequency scaling should be the same for electroencephalogram (EEG) and magnetoencephalogram (MEG) signals at low frequencies (<10 Hz). To test this prediction, we analyzed the spectrum of simultaneous EEG and MEG measurements in four human subjects. The frequency scaling of EEG displays coherent variations across the brain, in general between 1/f and 1/f ², and tends to be smaller in parietal/temporal regions. In a given region, although the variability of the frequency scaling exponent was higher for MEG compared to EEG, both signals consistently scale with a different exponent. In some cases, the scaling was similar, but only when the signal-to-noise ratio of the MEG was low. Several methods of noise correction for environmental and instrumental noise were tested, and they all increased the difference between EEG and MEG scaling. In conclusion, there is a significant difference in frequency scaling between EEG and MEG, which can be explained if the extracellular medium (including other layers such as dura matter and skull) is globally non-resistive.     

Effect of high dose cytosine arabinoside on quantitative EEG in patients with acute myeloid leukemia

Background EEG activity is considered an index of functional state of brain. Chemotherapy (CT), used for non-central nervous system (CNS) cancer, can cross the blood brain barrier and contribute to changes in the functional state of brain that can alter background EEG activity. Quantitative EEG (qEEG) is superior to conventional EEG in the detection of subtle alterations of EEG background activity and for this reason, the use of qEEG might assist the clinician in evaluating the possible effect of CT on the CNS. The nucleoside analog cytosine arabinoside (Ara-C) is one of the milestone chemotherapeutic agents used for treatment of acute myeloid leukemia (AML). Our observational study evaluates the possible effect of Ara-C on the qEEG of patients with AML, without CNS involvement. We conducted an observational study on newly diagnosed AML patients without CNS involvement, undergoing treatment with Ara-C to analyze the possible effect of Ara-C high doses on EEG background activity using qEEG analyses. A total of nine AML patients, 5 with Ara-C i.v. high dose (≥3 g/m² die), 4 with standard dose (100 mg/m² die) underwent qEEG (at rest, during hyperpnoea, mental arithmetic task and blocking reaction). We compared the EEG background activity of the two groups at baseline and after 6 months. Statistical analysis showed no significant differences between the two groups in mean relative power for all frequency bands, at rest and during hyperpnoea, mental arithmetic task and blocking reaction. Our data indicate that high dose Ara-C i.v. did not induce significant changes on EEG background activity in our patients. Future research in this area could include prospective studies that would combine qEEG and neuropsychological testing to assess the impact of CT on brain functions.     

A large-scale brain network mechanism for increased seizure propensity in Alzheimer’s disease

People with Alzheimer’s disease (AD) are 6-10 times more likely to develop seizures than the healthy aging population. Leading hypotheses largely consider hyperexcitability of local cortical tissue as primarily responsible for increased seizure prevalence in AD. However, in the general population of people with epilepsy, large-scale brain network organization additionally plays a role in determining seizure likelihood and phenotype. Here, we propose that alterations to large-scale brain network organization seen in AD may contribute to increased seizure likelihood. To test this hypothesis, we combine computational modelling with electrophysiological data using an approach that has proved informative in clinical epilepsy cohorts without AD. EEG was recorded from 21 people with probable AD and 26 healthy controls. At the time of EEG acquisition, all participants were free from seizures. Whole brain functional connectivity derived from source-reconstructed EEG recordings was used to build subject-specific brain network models of seizure transitions. As cortical tissue excitability was increased in the simulations, AD simulations were more likely to transition into seizures than simulations from healthy controls, suggesting an increased group-level probability of developing seizures at a future time for AD participants. We subsequently used the model to assess seizure propensity of different regions across the cortex. We found the most important regions for seizure generation were those typically burdened by amyloid-beta at the early stages of AD, as previously reported by in-vivo and post-mortem staging of amyloid plaques. Analysis of these spatial distributions also give potential insight into mechanisms of increased susceptibility to generalized (as opposed to focal) seizures in AD vs controls. This research suggests avenues for future studies testing patients with seizures, e.g. co-morbid AD/epilepsy patients, and comparisons with PET and MRI scans to relate regional seizure propensity with AD pathologies.     

Fully Automated Whole-Head Segmentation with Improved Smoothness and Continuity,with Theory Reviewed

Individualized current-flow models are needed for precise targeting of brain structures using transcranial electrical or magnetic stimulation (TES/TMS). The same is true for current-source reconstruction in electroencephalography and magnetoencephalography (EEG/MEG). The first step in generating such models is to obtain an accurate segmentation of individual head anatomy, including not only brain but also cerebrospinal fluid (CSF), skull and soft tissues, with a field of view (FOV) that covers the whole head. Currently available automated segmentation tools only provide results for brain tissues, have a limited FOV, and do not guarantee continuity and smoothness of tissues, which is crucially important for accurate current-flow estimates. Here we present a tool that addresses these needs. It is based on a rigorous Bayesian inference framework that combines image intensity model, anatomical prior (atlas) and morphological constraints using Markov random fields (MRF). The method is evaluated on 20 simulated and 8 real head volumes acquired with magnetic resonance imaging (MRI) at 1 mm3 resolution. We find improved surface smoothness and continuity as compared to the segmentation algorithms currently implemented in Statistical Parametric Mapping (SPM). With this tool, accurate and morphologically correct modeling of the whole-head anatomy for individual subjects may now be feasible on a routine basis. Code and data are fully integrated into SPM software tool and are made publicly available. In addition, a review on the MRI segmentation using atlas and the MRF over the last 20 years is also provided, with the general mathematical framework clearly derived.     

Semi-automatic classification of skeletal morphology in genetically altered mice using flat-panel volume computed tomography

Rapid progress in exploring the human and mouse genome has resulted in the generation of a multitude of mouse models to study gene functions in their biological context. However, effective screening methods that allow rapid noninvasive phenotyping of transgenic and knockout mice are still lacking. To identify murine models with bone alterations in vivo, we used flat-panel volume computed tomography (fpVCT) for high-resolution 3-D imaging and developed an algorithm with a computational intelligence system. First, we tested the accuracy and reliability of this approach by imaging discoidin domain receptor 2- (DDR2-) deficient mice, which display distinct skull abnormalities as shown by comparative landmark-based analysis. High-contrast fpVCT data of the skull with 200 microm isotropic resolution and 8-s scan time allowed segmentation and computation of significant shape features as well as visualization of morphological differences. The application of a trained artificial neuronal network to these datasets permitted a semi-automatic and highly accurate phenotype classification of DDR2-deficient compared to C57BL/6 wild-type mice. Even heterozygous DDR2 mice with only subtle phenotypic alterations were correctly determined by fpVCT imaging and identified as a new class. In addition, we successfully applied the algorithm to classify knockout mice lacking the DDR1 gene with no apparent skull deformities. Thus, this new method seems to be a potential tool to identify novel mouse phenotypes with skull changes from transgenic and knockout mice on the basis of random mutagenesis as well as from genetic models. However for this purpose, new neuronal networks have to be created and trained. In summary, the combination of fpVCT images with artificial neuronal networks provides a reliable, novel method for rapid, cost-effective, and noninvasive primary screening tool to detect skeletal phenotypes in mice.     

Multiphoton Microscopy of Cleared Mouse Brain Expressing YFP

Multiphoton microscopy of intrinsic fluorescence and second harmonic generation (SHG) of whole mouse organs is made possible by optically clearing the organ before imaging.^1,2 However, for organs that contain fluorescent proteins such as GFP and YFP, optical clearing protocols that use methanol dehydration and clear using benzyl alcohol:benzyl benzoate (BABB) while unprotected from light³ do not preserve the fluorescent signal. The protocol presented here is a novel way in which to perform whole organ optical clearing on mouse brain while preserving the fluorescence signal of YFP expressed in neurons. Altering the optical clearing protocol such that the organ is dehydrated using an ethanol graded series has been found to reduce the damage to the fluorescent proteins and preserve their fluorescent signal for multiphoton imaging.⁴ Using an optimized method of optical clearing with ethanol-based dehydration and clearing by BABB while shielded from light, we show high-resolution multiphoton images of yellow fluorescent protein (YFP) expression in the neurons of a mouse brain more than 2 mm beneath the tissue surface.     

Photoinduced resonance phenomena in the human electroencephalogram as a function of frequency, intensity, and duration of stimulation

The features of resonance phenomena in high-resolution EEG structure were analyzed for two intensities and three values of duration of exposure to 20 constant frequencies of intermittent photic stimulation in a range of 1-20 Hz with 1 Hz steps. It was shown that with a 6 s step duration, an irregular activation of multiple spectral EEG components for both light intensities occurs. With longer durations (12 and 18 s) of fixed-frequency stimulation, the EEG reactions are of resonance nature. Low-intensity flashes cause only the resonance activation of the intrinsic oscillator in the range of dominant alpha-EEG frequency. During a more intensive stimulation, the resonance EEG phenomena are observed for the whole range of stimulation frequencies. The interval of 6-12 s is supposed to be the relaxation period for a system of brain electrical activity generation. After this time, the low-intensity stimuli cause the adaptation of the system to light, whereas more intensive flashes cause more pronounced resonance EEG phenomena and physiological effects.     

Binding of 3H-naloxone in the mouse brain: effect of ions and tolerance development

The binding of ³H-naloxone (spec. act. 5.2 Ci/mmol) in a crude mitochondrial fraction of the whole mouse brain was examined. Binding was reversed by the narcotic agonists levorphanol, morphine and 1-methadone but not by dextrorphan. Levorphanol sensitive (specific) ³H-naloxone binding was blocked by Na+, Li+, Ca++, Mg++ and Mn++ but not by K+. When the crude mitochondrial fraction was separated on a discontinuous sucrose gradient, the highest activity of specific binding was found in the nerve ending particle fraction. Animals made physically dependent by 3 day morphine pellet implantation did not show an increased binding affinity for ³H-nalovxone. The implantation of a 10 mg naloxone pellet increased the apparent total number of binding sites on the second and third day of implantation.     

Resected Brain Tissue,Seizure Onset Zone and Quantitative EEG Measures: Towards Prediction of Post-Surgical Seizure Control

BackgroundEpilepsy surgery is a potentially curative treatment option for pharmacoresistent patients. If non-invasive methods alone do not allow to delineate the epileptogenic brain areas the surgical candidates undergo long-term monitoring with intracranial EEG. Visual EEG analysis is then used to identify the seizure onset zone for targeted resection as a standard procedure.MethodsDespite of its great potential to assess the epileptogenicty of brain tissue, quantitative EEG analysis has not yet found its way into routine clinical practice. To demonstrate that quantitative EEG may yield clinically highly relevant information we retrospectively investigated how post-operative seizure control is associated with four selected EEG measures evaluated in the resected brain tissue and the seizure onset zone. Importantly, the exact spatial location of the intracranial electrodes was determined by coregistration of pre-operative MRI and post-implantation CT and coregistration with post-resection MRI was used to delineate the extent of tissue resection. Using data-driven thresholding, quantitative EEG results were separated into normally contributing and salient channels.ResultsIn patients with favorable post-surgical seizure control a significantly larger fraction of salient channels in three of the four quantitative EEG measures was resected than in patients with unfavorable outcome in terms of seizure control (median over the whole peri-ictal recordings). The same statistics revealed no association with post-operative seizure control when EEG channels contributing to the seizure onset zone were studied.ConclusionsWe conclude that quantitative EEG measures provide clinically relevant and objective markers of target tissue, which may be used to optimize epilepsy surgery. The finding that differentiation between favorable and unfavorable outcome was better for the fraction of salient values in the resected brain tissue than in the seizure onset zone is consistent with growing evidence that spatially extended networks might be more relevant for seizure generation, evolution and termination than a single highly localized brain region (i.e. a “focus”) where seizures start.     

New perspectives in brain information processing

Brain cortex activity, as variously recorded by scalp or cortical electrodes in the electroencephalography (EEG) frequency range, probably reflects the basic strategy of brain information processing. Various hypotheses have been advanced to interpret this phenomenon, the most popular of which is that suitable combinations of excitatory and inhibitory neurons behave as assemblies of oscillators susceptible to synchronization and desynchronization. Implicit in this view is the assumption that EEG potentials are epiphenomena of action potentials, which is consistent with the argument that voltage variations in dendritic membranes reproduce the postsynaptic effects of targeting neurons. However, this classic argument does not really fit the discovery that firing synchronization over extended brain areas often appears to be established in about 1 ms, which is a small fraction of any EEG frequency component period. This is in contrast with the fact that all computational models of dynamic systems formed by more or less weakly interacting oscillators of near frequencies take more than one period to reach synchronization. The discovery that the somatodendritic membranes of specialized populations of neurons exhibit intrinsic subthreshold oscillations (ISOs) in the EEG frequency range, together with experimental evidence that short inhibitory stimuli are capable of resetting ISO phases, radically changes the scheme described above and paves the way to a novel view. This paper aims to elucidate the nature of ISO generation mechanisms, to explain the reasons for their reliability in starting and stopping synchronized firing, and to indicate their potential in brain information processing. The need for a repertoire of extraneuronal regulation mechanisms, putatively mediated by astrocytes, is also inferred. Lastly, the importance of ISOs for the brain as a parallel recursive machine is briefly discussed.     

A data-driven model of the generation of human EEG based on a spatially distributed stochastic wave equation

We discuss a model for the dynamics of the primary current density vector field within the grey matter of human brain. The model is based on a linear damped wave equation, driven by a stochastic term. By employing a realistically shaped average brain model and an estimate of the matrix which maps the primary currents distributed over grey matter to the electric potentials at the surface of the head, the model can be put into relation with recordings of the electroencephalogram (EEG). Through this step it becomes possible to employ EEG recordings for the purpose of estimating the primary current density vector field, i.e. finding a solution of the inverse problem of EEG generation. As a technique for inferring the unobserved high-dimensional primary current density field from EEG data of much lower dimension, a linear state space modelling approach is suggested, based on a generalisation of Kalman filtering, in combination with maximum-likelihood parameter estimation. The resulting algorithm for estimating dynamical solutions of the EEG inverse problem is applied to the task of localising the source of an epileptic spike from a clinical EEG data set; for comparison, we apply to the same task also a non-dynamical standard algorithm.