In the cortex, the interactions among neurons give rise to transient coherent activity patterns that underlie perception, cognition, and action. Recently, it was actively debated whether the most basic interactions, i.e., the pairwise correlations between neurons or groups of neurons, suffice to explain those observed activity patterns. So far, the evidence reported is controversial. Importantly, the overall organization of neuronal interactions and the mechanisms underlying their generation, especially those of high-order interactions, have remained elusive. Here we show that higher-order interactions are required to properly account for cortical dynamics such as ongoing neuronal avalanches in the alert monkey and evoked visual responses in the anesthetized cat. A Gaussian interaction model that utilizes the observed pairwise correlations and event rates and that applies intrinsic thresholding identifies those higher-order interactions correctly, both in cortical local field potentials and spiking activities. This allows for accurate prediction of large neuronal population activities as required, e.g., in brain-machine interface paradigms. Our results demonstrate that higher-order interactions are inherent properties of cortical dynamics and suggest a simple solution to overcome the apparent formidable complexity previously thought to be intrinsic to those interactions.

Imaging of neuronal depolarization in the brain is a major goal in neuroscience, but no technique currently exists that could image neural activity over milliseconds throughout the whole brain. Electrical impedance tomography (EIT) is an emerging medical imaging technique which can produce tomographic images of impedance changes with non-invasive surface electrodes. We report EIT imaging of impedance changes in rat somatosensory cerebral cortex with a resolution of 2 ms and < 200 μm during evoked potentials using epicortical arrays with 30 electrodes. Images were validated with local field potential recordings and current source-sink density analysis. Our results demonstrate that EIT can image neural activity in a volume 7 × 5 × 2 mm in somatosensory cerebral cortex with reduced invasiveness, greater resolution and imaging volume than other methods. Modeling indicates similar resolutions are feasible throughout the entire brain so this technique, uniquely, has the potential to image functional connectivity of cortical and subcortical structures.

Multifrequency electrical impedance tomography (MFEIT) exploits the dependence of tissue impedance on frequency to recover an image of conductivity. MFEIT could provide emergency diagnosis of pathologies such as acute stroke, brain injury and breast cancer. We present a method for performing MFEIT using spectral constraints. Boundary voltage data is employed directly to reconstruct the volume fraction distribution of component tissues using a nonlinear method. Given that the reconstructed parameter is frequency independent, this approach allows for the simultaneous use of all multifrequency data, thus reducing the degrees of freedom of the reconstruction problem. Furthermore, this method allows for the use of frequency difference data in a nonlinear reconstruction algorithm. Results from empirical phantom measurements suggest that our fraction reconstruction method points to a new direction for the development of multifrequency EIT algorithms in the case that the spectral constraints are known, and may provide a unifying framework for static EIT imaging.

Epilepsy affects approximately 50 million people worldwide, and 20–30% of these cases are refractory to antiepileptic drugs. Many patients with intractable epilepsy can benefit from surgical resection of the tissue generating the seizures; however, difficulty in precisely localising seizure foci has limited the number of patients undergoing surgery as well as potentially lowered its effectiveness. Here we demonstrate a novel imaging method for monitoring rapid changes in cerebral tissue impedance occurring during interictal and ictal activity, and show that it can reveal the propagation of pathological activity in the cortex. Cortical impedance was recorded simultaneously to ECoG using a 30-contact electrode mat placed on the exposed cortex of anaesthetised rats, in which interictal spikes (IISs) and seizures were induced by cortical injection of 4-aminopyridine (4-AP), picrotoxin or penicillin. We characterised the tissue impedance responses during IISs and seizures, and imaged these responses in the cortex using Electrical Impedance Tomography (EIT). We found a fast, transient drop in impedance occurring as early as 12 ms prior to the IISs, followed by a steep rise in impedance within ~ 120 ms of the IIS. EIT images of these impedance changes showed that they were co-localised and centred at a depth of 1 mm in the cortex, and that they closely followed the activity propagation observed in the surface ECoG signals. The fast, pre-IIS impedance drop most likely reflects synchronised depolarisation in a localised network of neurons, and the post-IIS impedance increase reflects the subsequent shrinkage of extracellular space caused by the intense activity. EIT could also be used to picture a steady rise in tissue impedance during seizure activity, which has been previously described. Thus, our results demonstrate that EIT can detect and localise different physiological changes during interictal and ictal activity and, in conjunction with ECoG, may in future improve the localisation of seizure foci in the clinical setting.

A method is presented for reconstructing images of fast neural evoked activity in rat cerebral cortex recorded with electrical impedance tomography (EIT) and a 6 × 5 mm(2) epicortical planar 30 electrode array. A finite element model of the rat brain and inverse solution with Tikhonov regularization were optimized in order to improve spatial resolution and accuracy. The optimized FEM mesh had 7 M tetrahedral elements, with finer resolution (0.05 mm) near the electrodes. A novel noise-based image processing technique based on t-test significance improved depth localization accuracy from 0.5 to 0.1 mm. With the improvements, a simulated perturbation 0.5 mm in diameter could be localized in a region 4 × 5 mm(2) under the centre of the array to a depth of 1.4 mm, thus covering all six layers of the cerebral cortex with an accuracy of <0.1 mm. Simulated deep brain hippocampal or thalamic activity could be localized with an accuracy of 0.5 mm with a 256 electrode array covering the brain. Parallel studies have achieved a temporal resolution of 2 ms for imaging fast neural activity by EIT during evoked activity; this encourages the view that fast neural EIT can now resolve the propagation of depolarization-related fast impedance changes in cerebral cortex and deeper in the brain with a resolution equal or greater to the dimension of a cortical column.