Neuronal synchrony and critical bistability: Mechanistic biomarkers for localizing the epileptogenic network.

OBJECTIVE: Postsurgical seizure freedom in drug-resistant epilepsy (DRE) patients varies from 30% to 80%, implying that in many cases the current approaches fail to fully map the epileptogenic zone (EZ). We aimed to advance a novel approach to better characterize epileptogenicity and investigate whether the EZ encompasses a broader epileptogenic network (EpiNet) beyond the seizure zone (SZ) that exhibits seizure activity. METHODS: We first used computational modeling to test putative complex systems-driven and systems neuroscience-driven mechanistic biomarkers for epileptogenicity. We then used these biomarkers to extract features from resting-state stereoelectroencephalograms recorded from DRE patients and trained supervised classifiers to localize the SZ against gold standard clinical localization. To further explore the prevalence of pathological features in an extended brain network outside of the clinically identified SZ, we also used unsupervised classification. RESULTS: Supervised SZ classification trained on individual features achieved accuracies of .6-.7 area under the receiver operating characteristic curve (AUC). Combining all criticality and synchrony features further improved the AUC to .85. Unsupervised classification discovered an EpiNet-like cluster of brain regions, in which 51% of brain regions were outside of the SZ. Brain regions in the EpiNet-like cluster engaged in interareal hypersynchrony and locally exhibited high-amplitude bistability and excessive inhibition, which was strikingly similar to the high seizure risk regime revealed by our computational modeling. SIGNIFICANCE: The finding that combining biomarkers improves SZ localization accuracy indicates that the novel mechanistic biomarkers for epileptogenicity employed here yield synergistic information. On the other hand, the discovery of SZ-like brain dynamics outside of the clinically defined SZ provides empirical evidence of an extended pathophysiological EpiNet.

Rhythmicity of neuronal oscillations delineates their cortical and spectral architecture.

Neuronal oscillations are commonly analyzed with power spectral methods that quantify signal amplitude, but not rhythmicity or 'oscillatoriness' per se. Here we introduce a new approach, the phase-autocorrelation function (pACF), for the direct quantification of rhythmicity. We applied pACF to human intracerebral stereoelectroencephalography (SEEG) and magnetoencephalography (MEG) data and uncovered a spectrally and anatomically fine-grained cortical architecture in the rhythmicity of single- and multi-frequency neuronal oscillations. Evidencing the functional significance of rhythmicity, we found it to be a prerequisite for long-range synchronization in resting-state networks and to be dynamically modulated during event-related processing. We also extended the pACF approach to measure 'burstiness' of oscillatory processes and characterized regions with stable and bursty oscillations. These findings show that rhythmicity is double-dissociable from amplitude and constitutes a functionally relevant and dynamic characteristic of neuronal oscillations.

EEG-based machine learning models for the prediction of phenoconversion time and subtype in isolated rapid eye movement sleep behavior disorder.

STUDY OBJECTIVES: Isolated rapid eye movement sleep behavior disorder (iRBD) is a prodromal stage of α-synucleinopathies and eventually phenoconverts to overt neurodegenerative diseases including Parkinson's disease (PD), dementia with Lewy bodies (DLB), and multiple system atrophy (MSA). Associations of baseline resting-state electroencephalography (EEG) with phenoconversion have been reported. In this study, we aimed to develop machine learning models to predict phenoconversion time and subtype using baseline EEG features in patients with iRBD. METHODS: At baseline, resting-state EEG and neurological assessments were performed on patients with iRBD. Calculated EEG features included spectral power, weighted phase lag index, and Shannon entropy. Three models were used for survival prediction, and four models were used for α-synucleinopathy subtype prediction. The models were externally validated using data from a different institution. RESULTS: A total of 236 iRBD patients were followed up for up to 8 years (mean 3.5 years), and 31 patients converted to α-synucleinopathies (16 PD, 9 DLB, 6 MSA). The best model for survival prediction was the random survival forest model with an integrated Brier score of 0.114 and a concordance index of 0.775. The K-nearest neighbor model was the best model for subtype prediction with an area under the receiver operating characteristic curve of 0.901. Slowing of the EEG was an important feature for both models. CONCLUSIONS: Machine learning models using baseline EEG features can be used to predict phenoconversion time and its subtype in patients with iRBD. Further research including large sample data from many countries is needed to make a more robust model.

NLP-based tools for localization of the epileptogenic zone in patients with drug-resistant focal epilepsy.

Epilepsy surgery is an option for people with focal onset drug-resistant (DR) seizures but a delayed or incorrect diagnosis of epileptogenic zone (EZ) location limits its efficacy. Seizure semiological manifestations and their chronological appearance contain valuable information on the putative EZ location but their interpretation relies on extensive experience. The aim of our work is to support the localization of EZ in DR patients automatically analyzing the semiological description of seizures contained in video-EEG reports. Our sample is composed of 536 descriptions of seizures extracted from Electronic Medical Records of 122 patients. We devised numerical representations of anamnestic records and seizures descriptions, exploiting Natural Language Processing (NLP) techniques, and used them to feed Machine Learning (ML) models. We performed three binary classification tasks: localizing the EZ in the right or left hemisphere, temporal or extra-temporal, and frontal or posterior regions. Our computational pipeline reached performances above 70% in all tasks. These results show that NLP-based numerical representation combined with ML-based classification models may help in localizing the origin of the seizures relying only on seizures-related semiological text data alone. Accurate early recognition of EZ could enable a more appropriate patient management and a faster access to epilepsy surgery to potential candidates.

Critical-like Brain Dynamics in a Continuum from Second- to First-Order Phase Transition.

The classic brain criticality hypothesis postulates that the brain benefits from operating near a continuous second-order phase transition. Slow feedback regulation of neuronal activity could, however, lead to a discontinuous first-order transition and thereby bistable activity. Observations of bistability in awake brain activity have nonetheless remained scarce and its functional significance unclear. Moreover, there is no empirical evidence to support the hypothesis that the human brain could flexibly operate near either a first- or second-order phase transition despite such a continuum being common in models. Here, using computational modeling, we found bistable synchronization dynamics to emerge through elevated positive feedback and occur exclusively in a regimen of critical-like dynamics. We then assessed bistability in vivo with resting-state MEG in healthy adults (7 females, 11 males) and stereo-electroencephalography in epilepsy patients (28 females, 36 males). This analysis revealed that a large fraction of the neocortices exhibited varying degrees of bistability in neuronal oscillations from 3 to 200 Hz. In line with our modeling results, the neuronal bistability was positively correlated with classic assessment of brain criticality across narrow-band frequencies. Excessive bistability was predictive of epileptic pathophysiology in the patients, whereas moderate bistability was positively correlated with task performance in the healthy subjects. These empirical findings thus reveal the human brain as a one-of-a-kind complex system that exhibits critical-like dynamics in a continuum between continuous and discontinuous phase transitions.SIGNIFICANCE STATEMENT In the model, while synchrony per se was controlled by connectivity, increasing positive local feedback led to gradually emerging bistable synchrony with scale-free dynamics, suggesting a continuum between second- and first-order phase transitions in synchrony dynamics inside a critical-like regimen. In resting-state MEG and SEEG, bistability of ongoing neuronal oscillations was pervasive across brain areas and frequency bands and was observed only with concurring critical-like dynamics as the modeling predicted. As evidence for functional relevance, moderate bistability was positively correlated with executive functioning in the healthy subjects, and excessive bistability was associated with epileptic pathophysiology. These findings show that critical-like neuronal dynamics in vivo involves both continuous and discontinuous phase transitions in a frequency-, neuroanatomy-, and state-dependent manner.

Brain criticality predicts individual levels of inter-areal synchronization in human electrophysiological data.

Neuronal oscillations and their synchronization between brain areas are fundamental for healthy brain function. Yet, synchronization levels exhibit large inter-individual variability that is associated with behavioral variability. We test whether individual synchronization levels are predicted by individual brain states along an extended regime of critical-like dynamics - the Griffiths phase (GP). We use computational modelling to assess how synchronization is dependent on brain criticality indexed by long-range temporal correlations (LRTCs). We analyze LRTCs and synchronization of oscillations from resting-state magnetoencephalography and stereo-electroencephalography data. Synchronization and LRTCs are both positively linearly and quadratically correlated among healthy subjects, while in epileptogenic areas they are negatively linearly correlated. These results show that variability in synchronization levels is explained by the individual position along the GP with healthy brain areas operating in its subcritical and epileptogenic areas in its supercritical side. We suggest that the GP is fundamental for brain function allowing individual variability while retaining functional advantages of criticality.

Influence of adaptive denoising on Diffusion Kurtosis Imaging at 3T and 7T.

BACKGROUND AND OBJECTIVE: Choosing the most appropriate denoising method to improve the quality of diagnostic images maximally is key in pre-processing of diffusion MRI images. Recent advancements in acquisition and reconstruction techniques have questioned traditional noise estimation methods favoring adaptive denoising frameworks, circumventing the need to know a priori information that is hardly available in a clinical setting. In this observational study, we compared two innovative adaptive techniques sharing some features, Patch2Self and Nlsam, through application on reference adult data at 3T and 7T. The primary aim was identifying the most effective method in case of Diffusion Kurtosis Imaging (DKI) data - particularly susceptible to noise and signal fluctuations - at 3T and 7T fields. A side goal consisted of investigating the dependence of kurtosis metrics' variability with respect to the magnetic field on the adopted denoising methodology. METHODS: For comparison purposes, we focused on qualitative and quantitative analysis of DKI data and related microstructural maps before and after applying the two denoising approaches. Specifically, we assessed computational efficiency, preservation of anatomical details via perceptual metrics, consistency of microstructure model fitting, alleviation of degeneracies in model estimation, and joint variability with varying field strength and denoising method. RESULTS: Accounting for all these factors, Patch2Self framework has turned out to be specifically suitable for DKI data, with improving performance at 7T. Nlsam method is more robust in alleviating degenerate black voxels while introducing some blurring, which in turn is reflected in an overall loss of image sharpness. Regarding the impact of denoising on field-dependent variability, both methods have been shown to make variations from standard to Ultra-High Field more concordant with theoretical evidence, claiming that kurtosis metrics are sensitive to susceptibility-induced background gradients, directly proportional to the magnetic field strength and sensitive to the microscopic distribution of iron and myelin. CONCLUSIONS: This study serves as a proof-of-concept stressing the need for an accurate choice of a denoising methodology, specifically tailored for the data under analysis and allowing higher spatial resolution acquisition within clinically compatible timings, with all the potential benefits that improving suboptimal quality of diagnostic images entails.

Phase and amplitude electroencephalography correlations change with disease progression in people with idiopathic rapid eye-movement sleep behavior disorder.

STUDY OBJECTIVES: Increased phase synchronization in electroencephalography (EEG) bands might reflect the activation of compensatory mechanisms of cognitive decline in people with neurodegenerative diseases. Here, we investigated whether altered large-scale couplings of brain oscillations could be linked to the balancing of cognitive decline in a longitudinal cohort of people with idiopathic rapid eye-movement sleep behavior disorder (iRBD). METHODS: We analyzed 18 patients (17 males, 69.7 ± 7.5 years) with iRBD undergoing high-density EEG (HD-EEG), presynaptic dopaminergic imaging, and clinical and neuropsychological (NPS) assessments at two time points (time interval 24.2 ± 5.9 months). We thus quantified the HD-EEG power distribution, orthogonalized amplitude correlation, and weighted phase-lag index at both time points and correlated them with clinical, NPS, and imaging data. RESULTS: Four patients phenoconverted at follow-up (three cases of parkinsonism and one of dementia). At the group level, NPS scores decreased over time, without reaching statistical significance. However, alpha phase synchronization increased and delta amplitude correlations decreased significantly at follow-up compared to baseline. Both large-scale network connectivity metrics were significantly correlated with NPS scores but not with sleep quality indices or presynaptic dopaminergic imaging data. CONCLUSIONS: These results suggest that increased alpha phase synchronization and reduced delta amplitude correlation may be considered electrophysiological signs of an active compensatory mechanism of cognitive impairment in people with iRBD. Large-scale functional modifications may be helpful biomarkers in the characterization of prodromal stages of alpha-synucleinopathies.

Diffusion Kurtosis Imaging of Neonatal Spinal Cord in Clinical Routine.

Diffusion kurtosis imaging (DKI) has undisputed advantages over the more classical diffusion magnetic resonance imaging (dMRI) as witnessed by the fast-increasing number of clinical applications and software packages widely adopted in brain imaging. However, in the neonatal setting, DKI is still largely underutilized, in particular in spinal cord (SC) imaging, because of its inherently demanding technological requirements. Due to its extreme sensitivity to non-Gaussian diffusion, DKI proves particularly suitable for detecting complex, subtle, fast microstructural changes occurring in this area at this early and critical stage of development, which are not identifiable with only DTI. Given the multiplicity of congenital anomalies of the spinal canal, their crucial effect on later developmental outcome, and the close interconnection between the SC region and the brain above, managing to apply such a method to the neonatal cohort becomes of utmost importance. This study will (i) mention current methodological challenges associated with the application of advanced dMRI methods, like DKI, in early infancy, (ii) illustrate the first semi-automated pipeline built on Spinal Cord Toolbox for handling the DKI data of neonatal SC, from acquisition setting to estimation of diffusion measures, through accurate adjustment of processing algorithms customized for adult SC, and (iii) present results of its application in a pilot clinical case study. With the proposed pipeline, we preliminarily show that DKI is more sensitive than DTI-related measures to alterations caused by brain white matter injuries in the underlying cervical SC.

Genuine cross-frequency coupling networks in human resting-state electrophysiological recordings.

Phase synchronization of neuronal oscillations in specific frequency bands coordinates anatomically distributed neuronal processing and communication. Typically, oscillations and synchronization take place concurrently in many distinct frequencies, which serve separate computational roles in cognitive functions. While within-frequency phase synchronization has been studied extensively, less is known about the mechanisms that govern neuronal processing distributed across frequencies and brain regions. Such integration of processing between frequencies could be achieved via cross-frequency coupling (CFC), either by phase-amplitude coupling (PAC) or by n:m-cross-frequency phase synchrony (CFS). So far, studies have mostly focused on local CFC in individual brain regions, whereas the presence and functional organization of CFC between brain areas have remained largely unknown. We posit that interareal CFC may be essential for large-scale coordination of neuronal activity and investigate here whether genuine CFC networks are present in human resting-state (RS) brain activity. To assess the functional organization of CFC networks, we identified brain-wide CFC networks at mesoscale resolution from stereoelectroencephalography (SEEG) and at macroscale resolution from source-reconstructed magnetoencephalography (MEG) data. We developed a novel, to our knowledge, graph-theoretical method to distinguish genuine CFC from spurious CFC that may arise from nonsinusoidal signals ubiquitous in neuronal activity. We show that genuine interareal CFC is present in human RS activity in both SEEG and MEG data. Both CFS and PAC networks coupled theta and alpha oscillations with higher frequencies in large-scale networks connecting anterior and posterior brain regions. CFS and PAC networks had distinct spectral patterns and opposing distribution of low- and high-frequency network hubs, implying that they constitute distinct CFC mechanisms. The strength of CFS networks was also predictive of cognitive performance in a separate neuropsychological assessment. In conclusion, these results provide evidence for interareal CFS and PAC being 2 distinct mechanisms for coupling oscillations across frequencies in large-scale brain networks.

Comparison of Methods to Identify Modules in Noisy or Incomplete Brain Networks.

Community structure, or "modularity," is a fundamentally important aspect in the organization of structural and functional brain networks, but their identification with community detection methods is confounded by noisy or missing connections. Although several methods have been used to account for missing data, the performance of these methods has not been compared quantitatively so far. In this study, we compared four different approaches to account for missing connections when identifying modules in binary and weighted networks using both Louvain and Infomap community detection algorithms. The four methods are "zeros," "row-column mean," "common neighbors," and "consensus clustering." Using Lancichinetti-Fortunato-Radicchi benchmark-simulated binary and weighted networks, we find that "zeros," "row-column mean," and "common neighbors" approaches perform well with both Louvain and Infomap, whereas "consensus clustering" performs well with Louvain but not Infomap. A similar pattern of results was observed with empirical networks from stereotactical electroencephalography data, except that "consensus clustering" outperforms other approaches on weighted networks with Louvain. Based on these results, we recommend any of the four methods when using Louvain on binary networks, whereas "consensus clustering" is superior with Louvain clustering of weighted networks. When using Infomap, "zeros" or "common neighbors" should be used for both binary and weighted networks. These findings provide a basis to accounting for noisy or missing connections when identifying modules in brain networks.

Cortical response to levodopa in Parkinson's disease patients with dyskinesias.

Levodopa-induced dyskinesias are a common and disabling side effect of dopaminergic therapy in Parkinson's disease, but their neural mechanisms in vivo are still poorly understood. Besides striatal pathology, the importance of cortical dysfunction has been increasingly recognized. The supplementary motor area in particular, may have a relevant role in dyskinesias onset given its involvement in endogenously generated actions. The aim of the present study was to investigate the levodopa-related cortical excitability changes along with the emergence of levodopa-induced peak-of-dose dyskinesias in subjects with Parkinson's disease. Thirteen patients without dyskinesias and ten with dyskinesias received 200/50 mg fast-acting oral levodopa/benserazide following overnight withdrawal (12 hr) from their dopaminergic medication. We targeted transcranial magnetic stimulation to the supplementary motor area, ipsilateral to the most dopamine-depleted striatum defined with single-photon emission computed tomography with [123 I]N-ω-fluoropropyl-2ß-carbomethoxy-3ß-(4-iodophenyl)nortropane, and recorded transcranial magnetic stimulation-evoked potentials with high-density electroencephalography before and at 30, 60, and 180 min after levodopa/benserazide intake. Clinical improvement from levodopa/benserazide paralleled the increase in cortical excitability in both groups. Subjects with dyskinesias showed higher fluctuation of cortical excitability in comparison to non-dyskinetic patients, possibly reflecting dyskinetic movements. Together with endogenous brain oscillation, levodopa-related dynamics of brain state could influence the therapeutic response of neuromodulatory interventions.

Phase matters: A role for the subthalamic network during gait.

The role of the subthalamic nucleus in human locomotion is unclear although relevant, given the troublesome management of gait disturbances with subthalamic deep brain stimulation in patients with Parkinson's disease. We investigated the subthalamic activity and inter-hemispheric connectivity during walking in eight freely-moving subjects with Parkinson's disease and bilateral deep brain stimulation. In particular, we compared the subthalamic power spectral densities and coherence, amplitude cross-correlation and phase locking value between resting state, upright standing, and steady forward walking. We observed a phase locking value drop in the ß-frequency band (≈13-35Hz) during walking with respect to resting and standing. This modulation was not accompanied by specific changes in subthalamic power spectral densities, which was not related to gait phases or to striatal dopamine loss measured with [123I]N-ω-fluoropropyl-2ß-carbomethoxy-3ß-(4-iodophenyl)nortropane and single-photon computed tomography. We speculate that the subthalamic inter-hemispheric desynchronization in the ß-frequency band reflects the information processing of each body side separately, which may support linear walking. This study also suggests that in some cases (i.e. gait) the brain signal, which could allow feedback-controlled stimulation, might derive from network activity.

Improvement in White Matter Tract Reconstruction with Constrained Spherical Deconvolution and Track Density Mapping in Low Angular Resolution Data: A Pediatric Study and Literature Review.

INTRODUCTION: Diffusion-weighted magnetic resonance imaging (DW-MRI) allows noninvasive investigation of brain structure in vivo. Diffusion tensor imaging (DTI) is a frequently used application of DW-MRI that assumes a single main diffusion direction per voxel, and is therefore not well suited for reconstructing crossing fiber tracts. Among the solutions developed to overcome this problem, constrained spherical deconvolution with probabilistic tractography (CSD-PT) has provided superior quality results in clinical settings on adult subjects; however, it requires particular acquisition parameters and long sequences, which may limit clinical usage in the pediatric age group. The aim of this work was to compare the results of DTI with those of track density imaging (TDI) maps and CSD-PT on data from neonates and children, acquired with low angular resolution and low b-value diffusion sequences commonly used in pediatric clinical MRI examinations. MATERIALS AND METHODS: We analyzed DW-MRI studies of 50 children (eight neonates aged 3-28 days, 20 infants aged 1-8 months, and 22 children aged 2-17 years) acquired on a 1.5 T Philips scanner using 34 gradient directions and a b-value of 1,000 s/mm2. Other sequence parameters included 60 axial slices; acquisition matrix, 128 × 128; average scan time, 5:34 min; voxel size, 1.75 mm × 1.75 mm × 2 mm; one b = 0 image. For each subject, we computed principal eigenvector (EV) maps and directionally encoded color TDI maps (DEC-TDI maps) from whole-brain tractograms obtained with CSD-PT; the cerebellar-thalamic, corticopontocerebellar, and corticospinal tracts were reconstructed using both CSD-PT and DTI. Results were compared by two neuroradiologists using a 5-point qualitative score. RESULTS: The DEC-TDI maps obtained presented higher anatomical detail than EV maps, as assessed by visual inspection. In all subjects, white matter (WM) tracts were successfully reconstructed using both tractography methodologies. The mean qualitative scores of all tracts obtained with CSD-PT were significantly higher than those obtained with DTI (p-value < 0.05 for all comparisons). CONCLUSION: CSD-PT can be successfully applied to DW-MRI studies acquired at 1.5 T with acquisition parameters adapted for pediatric subjects, thus providing TDI maps with greater anatomical detail. This methodology yields satisfactory results for clinical purposes in the pediatric age group.

Retrospective evaluation and SEEG trajectory analysis for interactive multi-trajectory planner assistant.

PURPOSE: Focal epilepsy is a neurological disease that can be surgically treated by removing area of the brain generating the seizures. The stereotactic electroencephalography (SEEG) procedure allows patient brain activity to be recorded in order to localize the onset of seizures through the placement of intracranial electrodes. The planning phase can be cumbersome and very time consuming, and no quantitative information is provided to neurosurgeons regarding the safety and efficacy of their trajectories. In this work, we present a novel architecture specifically designed to ease the SEEG trajectory planning using the 3D Slicer platform as a basis. METHODS: Trajectories are automatically optimized following criteria like vessel distance and insertion angle. Multi-trajectory optimization and conflict resolution are optimized through a selective brute force approach based on a conflict graph construction. Additionally, electrode-specific optimization constraints can be defined, and an advanced verification module allows neurosurgeons to evaluate the feasibility of the trajectory. RESULTS: A retrospective evaluation was performed using manually planned trajectories on 20 patients: the planning algorithm optimized and improved trajectories in 98% of cases. We were able to resolve and optimize the remaining 2% by applying electrode-specific constraints based on manual planning values. In addition, we found that the global parameters used discards 68% of the manual planned trajectories, even when they represent a safe clinical choice. CONCLUSIONS: Our approach improved manual planned trajectories in 98% of cases in terms of quantitative indexes, even when applying more conservative criteria with respect to actual clinical practice. The improved multi-trajectory strategy overcomes the previous work limitations and allows electrode optimization within a tolerable time span.

A new tool for touch-free patient registration for robot-assisted intracranial surgery: application accuracy from a phantom study and a retrospective surgical series.

OBJECTIVE The purpose of this study was to compare the accuracy of Neurolocate frameless registration system and frame-based registration for robotic stereoelectroencephalography (SEEG). METHODS The authors performed a 40-trajectory phantom laboratory study and a 127-trajectory retrospective analysis of a surgical series. The laboratory study was aimed at testing the noninferiority of the Neurolocate system. The analysis of the surgical series compared Neurolocate-based SEEG implantations with a frame-based historical control group. RESULTS The mean localization errors (LE) ± standard deviations (SD) for Neurolocate-based and frame-based trajectories were 0.67 ± 0.29 mm and 0.76 ± 0.34 mm, respectively, in the phantom study (p = 0.35). The median entry point LE was 0.59 mm (interquartile range [IQR] 0.25-0.88 mm) for Neurolocate-registration-based trajectories and 0.78 mm (IQR 0.49-1.08 mm) for frame-registration-based trajectories (p = 0.00002) in the clinical study. The median target point LE was 1.49 mm (IQR 1.06-2.4 mm) for Neurolocate-registration-based trajectories and 1.77 mm (IQR 1.25-2.5 mm) for frame-registration-based trajectories in the clinical study. All the surgical procedures were successful and uneventful. CONCLUSIONS The results of the phantom study demonstrate the noninferiority of Neurolocate frameless registration. The results of the retrospective surgical series analysis suggest that Neurolocate-based procedures can be more accurate than the frame-based ones. The safety profile of Neurolocate-based registration should be similar to that of frame-based registration. The Neurolocate system is comfortable, noninvasive, easy to use, and potentially faster than other registration devices.

SEEG assistant: a 3DSlicer extension to support epilepsy surgery.

BACKGROUND: In the evaluation of Stereo-Electroencephalography (SEEG) signals, the physicist's workflow involves several operations, including determining the position of individual electrode contacts in terms of both relationship to grey or white matter and location in specific brain regions. These operations are (i) generally carried out manually by experts with limited computer support, (ii) hugely time consuming, and (iii) often inaccurate, incomplete, and prone to errors. RESULTS: In this paper we present SEEG Assistant, a set of tools integrated in a single 3DSlicer extension, which aims to assist neurosurgeons in the analysis of post-implant structural data and hence aid the neurophysiologist in the interpretation of SEEG data. SEEG Assistant consists of (i) a module to localize the electrode contact positions using imaging data from a thresholded post-implant CT, (ii) a module to determine the most probable cerebral location of the recorded activity, and (iii) a module to compute the Grey Matter Proximity Index, i.e. the distance of each contact from the cerebral cortex, in order to discriminate between white and grey matter location of contacts. Finally, exploiting 3DSlicer capabilities, SEEG Assistant offers a Graphical User Interface that simplifies the interaction between the user and the tools. SEEG Assistant has been tested on 40 patients segmenting 555 electrodes, and it has been used to identify the neuroanatomical loci and to compute the distance to the nearest cerebral cortex for 9626 contacts. We also performed manual segmentation and compared the results between the proposed tool and gold-standard clinical practice. As a result, the use of SEEG Assistant decreases the post implant processing time by more than 2 orders of magnitude, improves the quality of results and decreases, if not eliminates, errors in post implant processing. CONCLUSIONS: The SEEG Assistant Framework for the first time supports physicists by providing a set of open-source tools for post-implant processing of SEEG data. Furthermore, SEEG Assistant has been integrated into 3D Slicer, a software platform for the analysis and visualization of medical images, overcoming limitations of command-line tools.

Modular co-organization of functional connectivity and scale-free dynamics in the human brain.

Scale-free neuronal dynamics and interareal correlations are emergent characteristics of spontaneous brain activity. How such dynamics and the anatomical patterns of neuronal connectivity are mutually related in brain networks has, however, remained unclear. We addressed this relationship by quantifying the network colocalization of scale-free neuronal activity-both neuronal avalanches and long-range temporal correlations (LRTCs)-and functional connectivity (FC) by means of intracranial and noninvasive human resting-state electrophysiological recordings. We found frequency-specific colocalization of scale-free dynamics and FC so that the interareal couplings of LRTCs and the propagation of neuronal avalanches were most pronounced in the predominant pathways of FC. Several control analyses and the frequency specificity of network colocalization showed that the results were not trivial by-products of either brain dynamics or our analysis approach. Crucially, scale-free neuronal dynamics and connectivity also had colocalized modular structures at multiple levels of network organization, suggesting that modules of FC would be endowed with partially independent dynamic states. These findings thus suggest that FC and scale-free dynamics-hence, putatively, neuronal criticality as well-coemerge in a hierarchically modular structure in which the modules are characterized by dense connectivity, avalanche propagation, and shared dynamic states.

Global and local complexity of intracranial EEG decreases during NREM sleep.

Key to understanding the neuronal basis of consciousness is the characterization of the neural signatures of changes in level of consciousness during sleep. Here we analysed three measures of dynamical complexity on spontaneous depth electrode recordings from 10 epilepsy patients during wakeful rest (WR) and different stages of sleep: (i) Lempel-Ziv complexity, which is derived from how compressible the data are; (ii) amplitude coalition entropy, which measures the variability over time of the set of channels active above a threshold; (iii) synchrony coalition entropy, which measures the variability over time of the set of synchronous channels. When computed across sets of channels that are broadly distributed across multiple brain regions, all three measures decreased substantially in all participants during early-night non-rapid eye movement (NREM) sleep. This decrease was partially reversed during late-night NREM sleep, while the measures scored similar to WR during rapid eye movement (REM) sleep. This global pattern was in almost all cases mirrored at the local level by groups of channels located in a single region. In testing for differences between regions, we found elevated signal complexity in the frontal lobe. These differences could not be attributed solely to changes in spectral power between conditions. Our results provide further evidence that the level of consciousness correlates with neural dynamical complexity.