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.

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.

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.

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.

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.

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.

Relationship of fast- and slow-timescale neuronal dynamics in human MEG and SEEG.

A growing body of evidence suggests that the neuronal dynamics are poised at criticality. Neuronal avalanches and long-range temporal correlations (LRTCs) are hallmarks of such critical dynamics in neuronal activity and occur at fast (subsecond) and slow (seconds to hours) timescales, respectively. The critical dynamics at different timescales can be characterized by their power-law scaling exponents. However, insight into the avalanche dynamics and LRTCs in the human brain has been largely obtained with sensor-level MEG and EEG recordings, which yield only limited anatomical insight and results confounded by signal mixing. We investigated here the relationship between the human neuronal dynamics at fast and slow timescales using both source-reconstructed MEG and intracranial stereotactical electroencephalography (SEEG). Both MEG and SEEG revealed avalanche dynamics that were characterized parameter-dependently by power-law or truncated-power-law size distributions. Both methods also revealed robust LRTCs throughout the neocortex with distinct scaling exponents in different functional brain systems and frequency bands. The exponents of power-law regimen neuronal avalanches and LRTCs were strongly correlated across subjects. Qualitatively similar power-law correlations were also observed in surrogate data without spatial correlations but with scaling exponents distinct from those of original data. Furthermore, we found that LRTCs in the autonomous nervous system, as indexed by heart-rate variability, were correlated in a complex manner with cortical neuronal avalanches and LRTCs in MEG but not SEEG. These scalp and intracranial data hence show that power-law scaling behavior is a pervasive but neuroanatomically inhomogeneous property of neuronal dynamics in central and autonomous nervous systems.

Automatic segmentation of deep intracerebral electrodes in computed tomography scans.

BACKGROUND: Invasive monitoring of brain activity by means of intracerebral electrodes is widely practiced to improve pre-surgical seizure onset zone localization in patients with medically refractory seizures. Stereo-Electroencephalography (SEEG) is mainly used to localize the epileptogenic zone and a precise knowledge of the location of the electrodes is expected to facilitate the recordings interpretation and the planning of resective surgery. However, the localization of intracerebral electrodes on post-implant acquisitions is usually time-consuming (i.e., manual segmentation), it requires advanced 3D visualization tools, and it needs the supervision of trained medical doctors in order to minimize the errors. In this paper we propose an automated segmentation algorithm specifically designed to segment SEEG contacts from a thresholded post-implant Cone-Beam CT volume (0.4 mm, 0.4 mm, 0.8 mm). The algorithm relies on the planned position of target and entry points for each electrode as a first estimation of electrode axis. We implemented the proposed algorithm into DEETO, an open source C++ prototype based on ITK library. RESULTS: We tested our implementation on a cohort of 28 subjects in total. The experimental analysis, carried out over a subset of 12 subjects (35 multilead electrodes; 200 contacts) manually segmented by experts, show that the algorithm: (i) is faster than manual segmentation (i.e., less than 1s/subject versus a few hours) (ii) is reliable, with an error of 0.5 mm ± 0.06 mm, and (iii) it accurately maps SEEG implants to their anatomical regions improving the interpretability of electrophysiological traces for both clinical and research studies. Moreover, using the 28-subject cohort we show here that the algorithm is also robust (error < 0.005 mm) against deep-brain displacements (< 12 mm) of the implanted electrode shaft from those planned before surgery. CONCLUSIONS: Our method represents, to the best of our knowledge, the first automatic algorithm for the segmentation of SEEG electrodes. The method can be used to accurately identify the neuroanatomical loci of SEEG electrode contacts by a non-expert in a fast and reliable manner.

Phase and amplitude correlations in resting-state activity in human stereotactical EEG recordings.

Inter-areal interactions of neuronal oscillations may be a key mechanism in the coordination of anatomically distributed neuronal processing. In humans, invasive stereo-electroencephalography (SEEG) is emerging as a reference method for electrophysiological recordings because of its excellent spatial and temporal resolution. It could thus be also considered an optimal method for mapping neuronal inter-areal interactions. However, the common bipolar (BP) referencing of SEEG data may both confuse signals from distinct sources and suppress true neuronal interactions whereas the alternative monopolar (MP) reference yields data contaminated by volume conduction. We advance here a novel referencing scheme for SEEG data where electrodes in grey matter are referenced to closest white-matter (CW) electrodes. Using a 22 subject cohort and these three referencing schemes, we observed that both inter-areal phase and amplitude correlations decayed as function of distance and frequency but remained significant and stable across distances up to 10cm. Furthermore, we found that deep and superficial cortical laminae exhibit distinct spectral profiles of oscillation power as well as distinct patterns of inter-areal phase and amplitude interactions. These effects were qualitatively similar in MP and CW but distorted with BP referencing. Importantly CW was not influenced by the apparent large-scale volume conduction inherent to MP. We thus demonstrate here that with CW referencing, the superior anatomical accuracy of SEEG can be leveraged to yield accurate quantification and qualitatively novel insight into phase and amplitude interactions in human brain activity.

Bistability breaks-off deterministic responses to intracortical stimulation during non-REM sleep.

During non-rapid eye movement (NREM) sleep (stage N3), when consciousness fades, cortico-cortical interactions are impaired while neurons are still active and reactive. Why is this? We compared cortico-cortical evoked-potentials recorded during wakefulness and NREM by means of time-frequency analysis and phase-locking measures in 8 epileptic patients undergoing intra-cerebral stimulations/recordings for clinical evaluation. We observed that, while during wakefulness electrical stimulation triggers a chain of deterministic phase-locked activations in its cortical targets, during NREM the same input induces a slow wave associated with an OFF-period (suppression of power>20Hz), possibly reflecting a neuronal down-state. Crucially, after the OFF-period, cortical activity resumes to wakefulness-like levels, but the deterministic effects of the initial input are lost, as indicated by a sharp drop of phase-locked activity. These findings suggest that the intrinsic tendency of cortical neurons to fall into a down-state after a transient activation (i.e. bistability) prevents the emergence of stable patterns of causal interactions among cortical areas during NREM. Besides sleep, the same basic neurophysiological dynamics may play a role in pathological conditions in which thalamo-cortical information integration and consciousness are impaired in spite of preserved neuronal activity.

A digital repository with an extensible data model for biobanking and genomic analysis management.

MOTIVATION: Molecular biology laboratories require extensive metadata to improve data collection and analysis. The heterogeneity of the collected metadata grows as research is evolving in to international multi-disciplinary collaborations and increasing data sharing among institutions. Single standardization is not feasible and it becomes crucial to develop digital repositories with flexible and extensible data models, as in the case of modern integrated biobanks management. RESULTS: We developed a novel data model in JSON format to describe heterogeneous data in a generic biomedical science scenario. The model is built on two hierarchical entities: processes and events, roughly corresponding to research studies and analysis steps within a single study. A number of sequential events can be grouped in a process building up a hierarchical structure to track patient and sample history. Each event can produce new data. Data is described by a set of user-defined metadata, and may have one or more associated files. We integrated the model in a web based digital repository with a data grid storage to manage large data sets located in geographically distinct areas. We built a graphical interface that allows authorized users to define new data types dynamically, according to their requirements. Operators compose queries on metadata fields using a flexible search interface and run them on the database and on the grid. We applied the digital repository to the integrated management of samples, patients and medical history in the BIT-Gaslini biobank. The platform currently manages 1800 samples of over 900 patients. Microarray data from 150 analyses are stored on the grid storage and replicated on two physical resources for preservation. The system is equipped with data integration capabilities with other biobanks for worldwide information sharing. CONCLUSIONS: Our data model enables users to continuously define flexible, ad hoc, and loosely structured metadata, for information sharing in specific research projects and purposes. This approach can improve sensitively interdisciplinary research collaboration and allows to track patients' clinical records, sample management information, and genomic data. The web interface allows the operators to easily manage, query, and annotate the files, without dealing with the technicalities of the data grid.

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.

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.

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.

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.

A repository based on a dynamically extensible data model supporting multidisciplinary research in neuroscience.

BACKGROUND: Robust, extensible and distributed databases integrating clinical, imaging and molecular data represent a substantial challenge for modern neuroscience. It is even more difficult to provide extensible software environments able to effectively target the rapidly changing data requirements and structures of research experiments. There is an increasing request from the neuroscience community for software tools addressing technical challenges about: (i) supporting researchers in the medical field to carry out data analysis using integrated bioinformatics services and tools; (ii) handling multimodal/multiscale data and metadata, enabling the injection of several different data types according to structured schemas; (iii) providing high extensibility, in order to address different requirements deriving from a large variety of applications simply through a user runtime configuration. METHODS: A dynamically extensible data structure supporting collaborative multidisciplinary research projects in neuroscience has been defined and implemented. We have considered extensibility issues from two different points of view. First, the improvement of data flexibility has been taken into account. This has been done through the development of a methodology for the dynamic creation and use of data types and related metadata, based on the definition of "meta" data model. This way, users are not constrainted to a set of predefined data and the model can be easily extensible and applicable to different contexts. Second, users have been enabled to easily customize and extend the experimental procedures in order to track each step of acquisition or analysis. This has been achieved through a process-event data structure, a multipurpose taxonomic schema composed by two generic main objects: events and processes. Then, a repository has been built based on such data model and structure, and deployed on distributed resources thanks to a Grid-based approach. Finally, data integration aspects have been addressed by providing the repository application with an efficient dynamic interface designed to enable the user to both easily query the data depending on defined datatypes and view all the data of every patient in an integrated and simple way. RESULTS: The results of our work have been twofold. First, a dynamically extensible data model has been implemented and tested based on a "meta" data-model enabling users to define their own data types independently from the application context. This data model has allowed users to dynamically include additional data types without the need of rebuilding the underlying database. Then a complex process-event data structure has been built, based on this data model, describing patient-centered diagnostic processes and merging information from data and metadata. Second, a repository implementing such a data structure has been deployed on a distributed Data Grid in order to provide scalability both in terms of data input and data storage and to exploit distributed data and computational approaches in order to share resources more efficiently. Moreover, data managing has been made possible through a friendly web interface. The driving principle of not being forced to preconfigured data types has been satisfied. It is up to users to dynamically configure the data model for the given experiment or data acquisition program, thus making it potentially suitable for customized applications. CONCLUSIONS: Based on such repository, data managing has been made possible through a friendly web interface. The driving principle of not being forced to preconfigured data types has been satisfied. It is up to users to dynamically configure the data model for the given experiment or data acquisition program, thus making it potentially suitable for customized applications.

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.

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.

XTENS - an eXTensible environment for neuroscience.

The XTENS (eXTensible Environment for NeuroScience) platform consists in an highly extensible environment for collaborative work that improve repeatability of experiment and provides data storage and analysis capabilities. The platform is divided in repository and application domains, branched in services with different purpose. The first domain is the central component of the platform and consists in a multimodal repository with a client-server architecture. The second one provides remote tools for image and signal visualization and analysis. The main issue for such a platform is not only to provide an extensible collaborative environment, but also to build a development platform for testing models and algorithms in neuroscience. For these reasons a Grid approach has been considered. Both computational and data Grids infrastructures can be exploited to analyze and share large datasets of distributed data. The architecture has been deployed to support surgical planning for patients affected by drug resistant epilepsy. In that scenario, a complex analysis for a fully multimodal dataset including different image modalities, EEG and video is required to localize the origin of the ictal discharge and critical brain areas. As first results, prototype versions of both repository and application domain components are presented.