Search

BACKGROUND: This study examined whether chronic neuropathic pain, modulated by a local anesthetic block, is associated with cortical magnetic field changes. METHODS: In a group of 20 patients with pain caused by unilateral traumatic peripheral nerve injury, a local block with lidocaine 1% was administered and the cortical effects were measured and compared with a control group. The global field power (GFP), describing distribution of cortical activation after median and ulnar nerve stimulation, was plotted and calculated. The effects on the affected hemisphere and the unaffected hemisphere (UH) before and after a block of the injured nerve were statistically evaluated. RESULTS: Major differences based on the GFP curves, at a component between 50 ms - 90 ms (M70), were found in patients: in the affected hemisphere the M70 GFP peak values were statistically significantly larger in comparison with the UH, and the GFP curves differed morphologically. Interestingly, the mean UH responses were reduced in comparison with the control group, a finding suggesting that the UH is also part of the cortical changes. At M70, the GFP curves and values in the affected hemisphere were modulated by a local block of the median or the ulnar nerve. The most likely location of cortical adaptation is in the primary somatosensory cortex. CONCLUSIONS: Cortical activation is enhanced in the affected hemisphere compared with the UH and is modulated by a local block. The UH in neuropathic pain changes as well. Evoked fields may offer an opportunity to monitor the effectiveness of treatments of neuropathic pain in humans. [ABSTRACT FROM AUTHOR]

Objective: The aim of glioblastoma surgery is to maximize the extent of resection while preserving functional integrity. Standards are lacking for surgical decision-making, and previous studies indicate treatment variations. These shortcomings reflect the need to evaluate larger populations from different care teams. In this study, the authors used probability maps to quantify and compare surgical decision-making throughout the brain by 12 neurosurgical teams for patients with glioblastoma., Methods: The study included all adult patients who underwent first-time glioblastoma surgery in 2012-2013 and were treated by 1 of the 12 participating neurosurgical teams. Voxel-wise probability maps of tumor location, biopsy, and resection were constructed for each team to identify and compare patient treatment variations. Brain regions with different biopsy and resection results between teams were identified and analyzed for patient functional outcome and survival., Results: The study cohort consisted of 1087 patients, of whom 363 underwent a biopsy and 724 a resection. Biopsy and resection decisions were generally comparable between teams, providing benchmarks for probability maps of resections and biopsies for glioblastoma. Differences in biopsy rates were identified for the right superior frontal gyrus and indicated variation in biopsy decisions. Differences in resection rates were identified for the left superior parietal lobule, indicating variations in resection decisions., Conclusions: Probability maps of glioblastoma surgery enabled capture of clinical practice decisions and indicated that teams generally agreed on which region to biopsy or to resect. However, treatment variations reflecting clinical dilemmas were observed and pinpointed by using the probability maps, which could therefore be useful for quality-of-care discussions between surgical teams for patients with glioblastoma.

Purpose: To improve the robustness of deep learning-based glioblastoma segmentation in a clinical setting with sparsified datasets., Materials and Methods: In this retrospective study, preoperative T1-weighted, T2-weighted, T2-weighted fluid-attenuated inversion recovery, and postcontrast T1-weighted MRI from 117 patients (median age, 64 years; interquartile range [IQR], 55-73 years; 76 men) included within the Multimodal Brain Tumor Image Segmentation (BraTS) dataset plus a clinical dataset (2012-2013) with similar imaging modalities of 634 patients (median age, 59 years; IQR, 49-69 years; 382 men) with glioblastoma from six hospitals were used. Expert tumor delineations on the postcontrast images were available, but for various clinical datasets, one or more sequences were missing. The convolutional neural network, DeepMedic, was trained on combinations of complete and incomplete data with and without site-specific data. Sparsified training was introduced, which randomly simulated missing sequences during training. The effects of sparsified training and center-specific training were tested using Wilcoxon signed rank tests for paired measurements., Results: A model trained exclusively on BraTS data reached a median Dice score of 0.81 for segmentation on BraTS test data but only 0.49 on the clinical data. Sparsified training improved performance (adjusted P < .05), even when excluding test data with missing sequences, to median Dice score of 0.67. Inclusion of site-specific data during sparsified training led to higher model performance Dice scores greater than 0.8, on par with a model based on all complete and incomplete data. For the model using BraTS and clinical training data, inclusion of site-specific data or sparsified training was of no consequence., Conclusion: Accurate and automatic segmentation of glioblastoma on clinical scans is feasible using a model based on large, heterogeneous, and partially incomplete datasets. Sparsified training may boost the performance of a smaller model based on public and site-specific data. Supplemental material is available for this article. Published under a CC BY 4.0 license., Competing Interests: Disclosures of Conflicts of Interest: R.S.E. disclosed no relevant relationships. M.V. Activities related to the present article: institution receives grant from Netherlands Organization for Scientific Research (NOW) (project number 10-10400-96-14003); institution receives grant from Dutch Cancer Society (VU2014-7113). Activities not related to the present article: disclosed no relevant relationships. Other relationships: disclosed no relevant relationships. D.M.J.M. disclosed no relevant relationships. F.B. Activities related to the present article: institution receives grant from Netherlands Organization for Scientific Research (NOW) for PICTURE project. Activities not related to the present article: author paid as board member of Roche, Bayer, and Merck (DSMB and Steering Committees); author is consultant for IXICO; institution receives grants from EU-H2020, UKMSS, NWO, MRC, HIHR-BRCUCLH; author receives royalties from Springer books; institution paid for educational presentations by Biogen (PML educational website). Other relationships: disclosed no relevant relationships. H.V. Activities related to the present article: disclosed no relevant relationships. Activities not related to the present article: institution is consultant for Merck (multiple sclerosis brain imaging consulting); institution receives grants from Teva, Novartis, Merck (research grants for multiple sclerosis studies). Other relationships: disclosed no relevant relationships. M.v.H. Activities related to the present article: institution receives grant from Dutch Cancer Society. Activities not related to the present article: institution receives grants from MRC Proximity to Discovery, MRC Studentship (PhD student); author paid for lectures by ESTRO (Honorarium -Sept 2018, ESTRO [ATP Meeting] €500 Honorarium -Feb 2019, ESTRO [IGRT Meeting] €500); author receives travel accommodations from IPEM, UKIO, AMC, AIFM, IGT Network, ICR, ESTRO, Elekta, NKI (conference and meeting expenses such as travel, accommodations, etc. IPEM July 2019, UKIO June 2019, July 2018, AMC May 2019, Jan 2019, BIR April 2019, March 2019, AIFM March 2019, IGT Network March 2019, Nov 2018, ICR March 2019, ESTRO Feb 2019, Sept 2018 Elekta June 2018, NKI July 2018). Other relationships: disclosed no relevant relationships. L.B. disclosed no relevant relationships. M.C.N. disclosed no relevant relationships. M.R. disclosed no relevant relationships. T.S. disclosed no relevant relationships. M.S.B. disclosed no relevant relationships. S.H.J. disclosed no relevant relationships. B.K. disclosed no relevant relationships. G.W. disclosed no relevant relationships. J.F. disclosed no relevant relationships. P.A.J.T.R. disclosed no relevant relationships. E.M. disclosed no relevant relationships. P.C.D.W.H. Activities related to the present article: institution receives grant from Dutch Cancer Society (VU2014-7113); BrainLab provided SmartBrush software; ZonMW (This research is part of the program Innovative Medical Devices Initiative with project number 10-10400-96-14003, which is financed by the Netherlands Organization for Scientific Research). Activities not related to the present article: disclosed no relevant relationships. Other relationships: disclosed no relevant relationships. J.C.d.M. Activities related to the present article: institution received grant from ZonMW (ZonMW paid a research grant to VUmc to cover salaries of PhD students involved. Activities not related to the present article: disclosed no relevant relationships. Other relationships: disclosed no relevant relationships. M.G.W. disclosed no relevant relationships., (2020 by the Radiological Society of North America, Inc.)

To summarize the distribution of glioma location within a patient population, registration of individual MR images to anatomical reference space is required. In this study, we quantified the accuracy of MR image registration to anatomical reference space with linear and non-linear transformations using estimated tumor targets of glioblastoma and lower-grade glioma, and anatomical landmarks at pre- and post-operative time-points using six commonly used registration packages (FSL, SPM5, DARTEL, ANTs, Elastix, and NiftyReg). Routine clinical pre- and post-operative, post-contrast T1-weighted images of 20 patients with glioblastoma and 20 with lower-grade glioma were collected. The 2009a Montreal Neurological Institute brain template was used as anatomical reference space. Tumors were manually segmented in the patient space and corresponding healthy tissue was delineated as a target volume in the anatomical reference space. Accuracy of the tumor alignment was quantified using the Dice score and the Hausdorff distance. To measure the accuracy of general brain alignment, anatomical landmarks were placed in patient and in anatomical reference space, and the landmark distance after registration was quantified. Lower-grade gliomas were registered more accurately than glioblastoma. Registration accuracy for pre- and post-operative MR images did not differ. SPM5 and DARTEL registered tumors most accurate, and FSL least accurate. Non-linear transformations resulted in more accurate general brain alignment than linear transformations, but tumor alignment was similar between linear and non-linear transformation. We conclude that linear transformation suffices to summarize glioma locations in anatomical reference space., (Copyright © 2020 Visser, Petr, Müller, Eijgelaar, Hendriks, Witte, Barkhof, van Herk, Mutsaerts, Vrenken, de Munck and De Witt Hamer.)

Objective: Decisions in glioblastoma surgery are often guided by presumed eloquence of the tumor location. The authors introduce the "expected residual tumor volume" (eRV) and the "expected resectability index" (eRI) based on previous decisions aggregated in resection probability maps. The diagnostic accuracy of eRV and eRI to predict biopsy decisions, resectability, functional outcome, and survival was determined., Methods: Consecutive patients with first-time glioblastoma surgery in 2012-2013 were included from 12 hospitals. The eRV was calculated from the preoperative MR images of each patient using a resection probability map, and the eRI was derived from the tumor volume. As reference, Sawaya's tumor location eloquence grades (EGs) were classified. Resectability was measured as observed extent of resection (EOR) and residual volume, and functional outcome as change in Karnofsky Performance Scale score. Receiver operating characteristic curves and multivariable logistic regression were applied., Results: Of 915 patients, 674 (74%) underwent a resection with a median EOR of 97%, functional improvement in 71 (8%), functional decline in 78 (9%), and median survival of 12.8 months. The eRI and eRV identified biopsies and EORs of at least 80%, 90%, or 98% better than EG. The eRV and eRI predicted observed residual volumes under 10, 5, and 1 ml better than EG. The eRV, eRI, and EG had low diagnostic accuracy for functional outcome changes. Higher eRV and lower eRI were strongly associated with shorter survival, independent of known prognostic factors., Conclusions: The eRV and eRI predict biopsy decisions, resectability, and survival better than eloquence grading and may be useful preoperative indices to support surgical decisions.

Objective: We investigated the potential added value of high-density resting-state EEG by addressing differences with healthy individuals and associations with Fugl-Meyer motor assessment of the upper extremity (FM-UE) scores in chronic stroke., Methods: Twenty-one chronic stroke survivors with initial upper limb paresis and eleven matched controls were included. Group differences regarding resting-state EEG parameters (Delta Alpha ratio (DAR) and pairwise-derived Brain Symmetry Index (BSI)) and associations with FM-UE were investigated, as well as lateralization of BSI and the value of different frequency bands., Results: Chronic stroke survivors showed higher BSI compared to controls (p < 0.001), most pronounced in delta and theta frequency bands (p < 0.0001; p < 0.001). In the delta and theta band, BSI was significantly negatively associated with FM-UE (both p = 0.008) corrected for confounding factors. DAR showed no differences between groups nor association with FM-UE. Directional BSI showed increased power in the affected versus the unaffected hemisphere., Conclusions: Asymmetry in spectral power between hemispheres was present in chronic stroke, most pronounced in low frequencies and related to upper extremity motor function deficit., Significance: BSI is related to motor impairment and higher in chronic stroke patients compared to healthy controls, suggesting that BSI may be a marker of selective motor control., (Copyright © 2019 International Federation of Clinical Neurophysiology. Published by Elsevier B.V. All rights reserved.)

Distinguishing tumors from normal brain cells is important but challenging in glioma surgery due to the lack of clear interfaces between the two. The ability of label-free third harmonic generation (THG) microscopy in combination with automated image analysis to quantitatively detect glioma infiltration in fresh, unprocessed tissue in real time is assessed. The THG images reveal increased cellularity in grades II-IV glioma samples from 23 patients, as confirmed by subsequent hematoxylin and eosin histology. An automated image quantification workflow is presented for quantitative assessment of the imaged cellularity as a reflection of the degree of glioma invasion. The cellularity is validated in three ways: 1) Quantitative comparison of THG imaging with fluorescence microscopy of nucleus-stained samples demonstrates that THG reflects the true tissue cellularity. 2) Thresholding of THG cellularity differentiates normal brain from glioma infiltration, with 96.6% sensitivity and 95.5% specificity, in nearly perfect (93%) agreement with pathologists. 3) In one patient, a good correlation between THG cellularity and preoperative magnetic resonance and positron emission tomography imaging is demonstrated. In conclusion, quantitative real-time THG microscopy accurately assesses glioma infiltration in ex vivo human brain samples, and therefore holds strong potential for improving the accuracy of surgical resection., Competing Interests: The authors declare no conflict of interest.

Background: Hippocampal avoidance prophylactic cranial irradiation (HA-PCI) techniques have been developed to reduce radiation damage to the hippocampus. An inter-observer hippocampus delineation analysis was performed and the influence of the delineation variability on dose to the hippocampus was studied., Materials and Methods: For five patients, seven observers delineated both hippocampi on brain MRI. The intra-class correlation (ICC) with absolute agreement and the generalized conformity index (CI gen ) were computed. Median surfaces over all observers' delineations were created for each patient and regional outlining differences were analysed. HA-PCI dose plans were made from the median surfaces and we investigated whether dose constraints in the hippocampus could be met for all delineations., Results: The ICC for the left and right hippocampus was 0.56 and 0.69, respectively, while the CI gen ranged from 0.55 to 0.70. The posterior and anterior-medial hippocampal regions had most variation with SDs ranging from approximately 1 to 2.5 mm. The mean dose (D mean ) constraint was met for all delineations, but for the dose received by 1% of the hippocampal volume (D 1% ) violations were observed., Conclusion: The relatively low ICC and CI gen indicate that delineation variability among observers for both left and right hippocampus was large. The posterior and anterior-medial border have the largest delineation inaccuracy. The hippocampus D mean constraint was not violated.

Objective: The objective is to present a proof-of-concept of a semi-automatic method to reduce hippocampus segmentation time on magnetic resonance images (MRI)., Materials and Methods: FAst Segmentation Through SURface Fairing (FASTSURF) is based on a surface fairing technique which reconstructs the hippocampus from sparse delineations. To validate FASTSURF, simulations were performed in which sparse delineations extracted from full manual segmentations served as input. On three different datasets with different diagnostic groups, FASTSURF hippocampi were compared to the original segmentations using Jaccard overlap indices and percentage volume differences (PVD). In one data set for which back-to-back scans were available, unbiased estimates of overlap and PVD were obtained. Using longitudinal scans, we compared hippocampal atrophy rates measured by manual, FASTSURF and two automatic segmentations (FreeSurfer and FSL-FIRST)., Results: With only seven input contours, FASTSURF yielded mean Jaccard indices ranging from 72(±4.3)% to 83(±2.6)% and PVDs ranging from 0.02(±2.40)% to 3.2(±3.40)% across the three datasets. Slightly poorer results were obtained for the unbiased analysis, but the performance was still considerably better than both tested automatic methods with only five contours., Conclusions: FASTSURF segmentations have high accuracy and require only a fraction of the delineation effort of fully manual segmentation. Atrophy rate quantification based on completely manual segmentation is well reproduced by FASTSURF. Therefore, FASTSURF is a promising tool to be implemented in clinical workflow, provided a future prospective validation confirms our findings., Competing Interests: HV has received research grants from Novartis, Teva, MerckSerono and Pfizer, and a speaker honorarium from Novartis, but these funders did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. All funds were paid directly to his institution. There are no patents, products in development or marketed products to declare. This does not alter our adherence to all the PLOS ONE policies on sharing data and materials.

Brain volume loss, or atrophy, has been proven to be an important characteristic of neurological diseases such as Alzheimer's disease and multiple sclerosis. To use atrophy rate as a reliable clinical biomarker and to increase statistical power in clinical treatment trials, measurement variability needs to be minimized. Among other sources, systematic differences between different MR scanners are suspected to contribute to this variability. In this study we developed and performed initial validation tests of an MR-compatible phantom and analysis software for robust and reliable evaluation of the brain volume loss. The phantom contained three inflatable models of brain structures, i.e. cerebral hemisphere, putamen, and caudate nucleus. Software to reliably quantify volumes form the phantom images was also developed. To validate the method, the phantom was imaged using 3D T1-weighted protocols at three clinical 3T MR scanners from different vendors. Calculated volume change from MRI was compared with the known applied volume change using ICC and mean absolute difference. As assessed by the ICC, the agreement between our developed software and the applied volume change for different structures ranged from 0.999-1 for hemisphere, 0.976-0.998 for putamen, and 0.985-0.999 for caudate nucleus. The mean absolute differences between measured and applied volume change were 109-332 μL for hemisphere, 2.9-11.9 μL for putamen, and 2.2-10.1 μL for caudate nucleus. This method offers a reliable and robust measurement of volume change using MR images and could potentially be used to standardize clinical measurement of atrophy rates., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: Tumor segmentation of glioma on MRI is a technique to monitor, quantify and report disease progression. Manual MRI segmentation is the gold standard but very labor intensive. At present the quality of this gold standard is not known for different stages of the disease, and prior work has mainly focused on treatment-naive glioblastoma. In this paper we studied the inter-rater agreement of manual MRI segmentation of glioblastoma and WHO grade II-III glioma for novices and experts at three stages of disease. We also studied the impact of inter-observer variation on extent of resection and growth rate., Methods: In 20 patients with WHO grade IV glioblastoma and 20 patients with WHO grade II-III glioma (defined as non-glioblastoma) both the enhancing and non-enhancing tumor elements were segmented on MRI, using specialized software, by four novices and four experts before surgery, after surgery and at time of tumor progression. We used the generalized conformity index (GCI) and the intra-class correlation coefficient (ICC) of tumor volume as main outcome measures for inter-rater agreement., Results: For glioblastoma, segmentations by experts and novices were comparable. The inter-rater agreement of enhancing tumor elements was excellent before surgery (GCI 0.79, ICC 0.99) poor after surgery (GCI 0.32, ICC 0.92), and good at progression (GCI 0.65, ICC 0.91). For non-glioblastoma, the inter-rater agreement was generally higher between experts than between novices. The inter-rater agreement was excellent between experts before surgery (GCI 0.77, ICC 0.92), was reasonable after surgery (GCI 0.48, ICC 0.84), and good at progression (GCI 0.60, ICC 0.80). The inter-rater agreement was good between novices before surgery (GCI 0.66, ICC 0.73), was poor after surgery (GCI 0.33, ICC 0.55), and poor at progression (GCI 0.36, ICC 0.73). Further analysis showed that the lower inter-rater agreement of segmentation on postoperative MRI could only partly be explained by the smaller volumes and fragmentation of residual tumor. The median interquartile range of extent of resection between raters was 8.3% and of growth rate was 0.22 mm/year., Conclusion: Manual tumor segmentations on MRI have reasonable agreement for use in spatial and volumetric analysis. Agreement in spatial overlap is of concern with segmentation after surgery for glioblastoma and with segmentation of non-glioblastoma by non-experts., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

Objective: To compare the performance of different methods for determining hippocampal atrophy rates using longitudinal MRI scans in aging and Alzheimer's disease (AD)., Background: Quantifying hippocampal atrophy caused by neurodegenerative diseases is important to follow the course of the disease. In dementia, the efficacy of new therapies can be partially assessed by measuring their effect on hippocampal atrophy. In radiotherapy, the quantification of radiation-induced hippocampal volume loss is of interest to quantify radiation damage. We evaluated plausibility, reproducibility and sensitivity of eight commonly used methods to determine hippocampal atrophy rates using test-retest scans., Materials and Methods: Manual, FSL-FIRST, FreeSurfer, multi-atlas segmentation (MALF) and non-linear registration methods (Elastix, NiftyReg, ANTs and MIRTK) were used to determine hippocampal atrophy rates on longitudinal T1-weighted MRI from the ADNI database. Appropriate parameters for the non-linear registration methods were determined using a small training dataset (N = 16) in which two-year hippocampal atrophy was measured using test-retest scans of 8 subjects with low and 8 subjects with high atrophy rates. On a larger dataset of 20 controls, 40 mild cognitive impairment (MCI) and 20  AD patients, one-year hippocampal atrophy rates were measured. A repeated measures ANOVA analysis was performed to determine differences between controls, MCI and AD patients. For each method we calculated effect sizes and the required sample sizes to detect one-year volume change between controls and MCI (N CTRL&#95;MCI ) and between controls and AD (N CTRL&#95;AD ). Finally, reproducibility of hippocampal atrophy rates was assessed using within-session rescans and expressed as an average distance measure D Ave , which expresses the difference in atrophy rate, averaged over all subjects. The same D Ave was used to determine the agreement between different methods., Results: Except for MALF, all methods detected a significant group difference between CTRL and AD, but none could find a significant difference between the CTRL and MCI. FreeSurfer and MIRTK required the lowest sample sizes (FreeSurfer: N CTRL&#95;MCI  = 115, N CTRL&#95;AD  = 17 with D Ave  = 3.26%; MIRTK: N CTRL&#95;MCI  = 97, N CTRL&#95;AD  = 11 with D Ave  = 3.76%), while ANTs was most reproducible (N CTRL&#95;MCI  = 162, N CTRL&#95;AD  = 37 with D Ave  = 1.06%), followed by Elastix (N CTRL&#95;MCI  = 226, N CTRL&#95;AD  = 15 with D Ave  = 1.78%) and NiftyReg (N CTRL&#95;MCI  = 193, N CTRL&#95;AD  = 14 with D Ave  = 2.11%). Manually measured hippocampal atrophy rates required largest sample sizes to detect volume change and were poorly reproduced (N CTRL&#95;MCI  = 452, N CTRL&#95;AD  = 87 with D Ave  = 12.39%). Atrophy rates of non-linear registration methods also agreed best with each other., Discussion and Conclusion: Non-linear registration methods were most consistent in determining hippocampal atrophy and because of their better reproducibility, methods, such as ANTs, Elastix and NiftyReg, are preferred for determining hippocampal atrophy rates on longitudinal MRI. Since performances of non-linear registration methods are well comparable, the preferred method would mostly depend on computational efficiency., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

Purpose: The aim of glioblastoma surgery is to maximize the extent of resection while preserving functional integrity, which depends on the location within the brain. A standard to compare these decisions is lacking. We present a volumetric voxel-wise method for direct comparison between two multidisciplinary teams of glioblastoma surgery decisions throughout the brain., Methods: Adults undergoing first-time glioblastoma surgery from 2012 to 2013 performed by two neuro-oncologic teams were included. Patients had had a diagnostic biopsy or resection. Preoperative tumors and postoperative residues were segmented on magnetic resonance imaging in three dimensions and registered to standard brain space. Voxel-wise probability maps of tumor location, biopsy, and resection were constructed for each team to compare patient referral bias, indication variation, and treatment variation. To evaluate the quality of care, subgroups of differentially resected brain regions were analyzed for survival and functional outcome., Results: One team included 101 patients, and the other included 174; 91 tumors were biopsied, and 181 were resected. Patient characteristics were largely comparable between teams. Distributions of tumor locations were dissimilar, suggesting referral bias. Distributions of biopsies were similar, suggesting absence of indication variation. Differentially resected regions were identified in the anterior limb of the right internal capsule and the right caudate nucleus, indicating treatment variation. Patients with (n = 12) and without (n = 6) surgical removal in these regions had similar overall survival and similar permanent neurologic deficits., Conclusion: Probability maps of tumor location, biopsy, and resection provide additional information that can inform surgical decision making across multidisciplinary teams for patients with glioblastoma.

Third harmonic generation (THG) microscopy shows great potential for instant pathology of brain tissue during surgery. However, the rich morphologies contained and the noise associated makes image restoration, necessary for quantification of the THG images, challenging. Anisotropic diffusion filtering (ADF) has been recently applied to restore THG images of normal brain, but ADF is hard-to-code, time-consuming and only reconstructs salient edges. This work overcomes these drawbacks by expressing ADF as a tensor regularized total variation model, which uses the Huber penalty and the L 1 norm for tensor regularization and fidelity measurement, respectively. The diffusion tensor is constructed from the structure tensor of ADF yet the tensor decomposition is performed only in the non-flat areas. The resulting model is solved by an efficient and easy-to-code primal-dual algorithm. Tests on THG brain tumor images show that the proposed model has comparable denoising performance as ADF while it much better restores weak edges and it is up to 60% more time efficient., (© 2018 The Authors. Journal of Biophotonics published by WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Introduction: Amyloid plaque deposition in the brain is an early pathological change in Alzheimer's disease (AD), causing disrupted synaptic connections. Brain network disruptions in AD have been demonstrated with eigenvector centrality (EC), a measure that identifies central regions within networks. Carrying an apolipoprotein (APOE)-ε4 allele is a genetic risk for AD, associated with increased amyloid deposition. We studied whether APOE-ε4 carriership is associated with EC disruptions in cognitively normal individuals., Methods: A total of 261 healthy middle-aged to older adults (mean age 56.6 years) were divided into high-risk (APOE-ε4 carriers) and low-risk (noncarriers) groups. EC was computed from resting-state functional MRI data. Clusters of between-group differences were assessed with a permutation-based method. Correlations between cluster mean EC with brain volume, CSF biomarkers, and psychological test scores were assessed., Results: Decreased EC in the visual cortex was associated with APOE-ε4 carriership, a genetic risk factor for AD. EC differences were correlated with age, CSF amyloid levels, and scores on the trail-making and 15-object recognition tests., Conclusion: Our findings suggest that the APOE-ε4 genotype affects brain connectivity in regions previously found to be abnormal in AD as a sign of very early disease-related pathology. These differences were too subtle in healthy elderly to use EC for single-subject prediction of APOE genotype., (© 2018 The Authors. Brain and Behavior published by Wiley Periodicals, Inc.)

In searching for clinical biomarkers of the somatosensory function, we studied reproducibility of somatosensory potentials (SEP) evoked by finger stimulation in healthy subjects. SEPs induced by electrical stimulation and especially after median nerve stimulation is a method widely used in the literature. It is unclear, however, if the EEG recordings after finger stimulation are reproducible within the same subject. We tested in five healthy subjects the consistency and reproducibility of responses through bootstrapping as well as test-retest recordings. We further evaluated the possibility to discriminate activity of different fingers both at electrode and at source level. The lack of consistency and reproducibility suggest responses to finger stimulation to be unreliable, even with reasonably high signal-to-noise ratio and adequate number of trials. At sources level, somatotopic arrangement of the fingers representation was only found in one of the subjects. Although finding distinct locations of the different fingers activation was possible, our protocol did not allow for non-overlapping dipole representations of the fingers. We conclude that despite its theoretical advantages, we cannot recommend the use of somatosensory potentials evoked by finger stimulation to extract clinical biomarkers.

Background: The accuracy of source reconstruction depends on the spatial configuration of the neural sources underlying encephalographic signals, the temporal distance of the source activity, the level and structure of noise in the recordings, and - of course - on the employed inverse method. This plenitude of factors renders a definition of 'spatial resolution' of the electro-encephalogram (EEG) a challenge., New Method: A proper definition of spatial resolution requires a ground truth. We used data from numerical simulations of two dipoles changed with waveforms resembling somatosensory evoked potentials peaking at 20, 30, 50, 100 ms. We varied inter-dipole distances and added noise to the simulated scalp recordings with distinct signal-to-noise ratios (SNRs). Prior to inverse modeling we pre-whitened the simulated data and the leadfield. We tested a two-dipole fit, sc-MUSIC, and sc-eLORETA and assessed their accuracy via the distance between the simulated and estimated sources., Results: To quantify the spatial resolution of EEG, we introduced the notion of separability, i.e. the separation of two dipolar sources with a certain inter-dipole distance. Our results indicate separability of two sources in the presence of realistic noise with SNR up to 3 if they are 11 mm or further apart., Comparison With Existing Methods: In the presence of realistic noise, spatial pre-whitening appears mandatory preprocessing step irrespective of the inverse method employed., Conclusions: Separability is a legitimate measure to quantify EEG's spatial resolution. An optimal resolution in source reconstruction requires spatial pre-whitening as a crucial pre-processing step., (Copyright © 2018 Elsevier B.V. All rights reserved.)

Third harmonic generation (THG) microscopy is a label-free imaging technique that shows great potential for rapid pathology of brain tissue during brain tumor surgery. However, the interpretation of THG brain images should be quantitatively linked to images of more standard imaging techniques, which so far has been done qualitatively only. We establish here such a quantitative link between THG images of mouse brain tissue and all-nuclei-highlighted fluorescence images, acquired simultaneously from the same tissue area. For quantitative comparison of a substantial pair of images, we present here a segmentation workflow that is applicable for both THG and fluorescence images, with a precision of 91.3 % and 95.8 % achieved respectively. We find that the correspondence between the main features of the two imaging modalities amounts to 88.9 %, providing quantitative evidence of the interpretation of dark holes as brain cells. Moreover, 80 % bright objects in THG images overlap with nuclei highlighted in the fluorescence images, and they are 2 times smaller than the dark holes, showing that cells of different morphologies can be recognized in THG images. We expect that the described quantitative comparison is applicable to other types of brain tissue and with more specific staining experiments for cell type identification., (© 2017 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Maladaptive changes in the involvement of striatal and frontal cortical regions in drug use are thought to underlie the progression to habitual drug use and loss of cognitive control over drug intake that occur with accumulating drug experience. The present experiments focus on changes in neuronal activity in these regions associated with short-term (10 days) and long-term (60 days) self-administration of cocaine. Quantitative in situ hybridization for the immediate early gene Mkp1 was combined with statistical parametric mapping to assess the distribution of neuronal activity. We hypothesized that neuronal activity in striatum would increase in its dorsal part and that activity in frontal cortex would decrease with prolonged cocaine self-administration experience. Expression of Mkp1 was profoundly increased after cocaine self-administration, and the magnitude of this effect was greater after short-term compared to long-term self-administration. Increased neuronal activity was seen in both dorsal and ventral sectors of the striatum after 10 days exposure to cocaine. However, enhanced activity was restricted to dorsomedial and dorsocentral striatum after 60 days cocaine self-administration. In virtually all medial prefrontal and most orbitofrontal areas, increased expression of Mkp1 was observed after 10 days of cocaine taking, whereas after 60 days, enhanced expression was restricted to caudal parts of medial prefrontal and caudomedial parts of orbitofrontal cortex. Our data reveal functional changes in cellular activity in striatum and frontal cortex with increasing cocaine self-administration experience. These changes might reflect the neural processes that underlie the descent from recreational drug taking to compulsive cocaine use.

Motivation: The morphologies contained in 3D third harmonic generation (THG) images of human brain tissue can report on the pathological state of the tissue. However, the complexity of THG brain images makes the usage of modern image processing tools, especially those of image filtering, segmentation and validation, to extract this information challenging., Results: We developed a salient edge-enhancing model of anisotropic diffusion for image filtering, based on higher order statistics. We split the intrinsic 3-phase segmentation problem into two 2-phase segmentation problems, each of which we solved with a dedicated model, active contour weighted by prior extreme. We applied the novel proposed algorithms to THG images of structurally normal ex-vivo human brain tissue, revealing key tissue components-brain cells, microvessels and neuropil, enabling statistical characterization of these components. Comprehensive comparison to manually delineated ground truth validated the proposed algorithms. Quantitative comparison to second harmonic generation/auto-fluorescence images, acquired simultaneously from the same tissue area, confirmed the correctness of the main THG features detected., Availability and Implementation: The software and test datasets are available from the authors., Contact: z.zhang@vu.nl., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author 2017. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com)

Purpose: Precise and reproducible hippocampus outlining is important to quantify hippocampal atrophy caused by neurodegenerative diseases and to spare the hippocampus in whole brain radiation therapy when performing prophylactic cranial irradiation or treating brain metastases. This study aimed to quantify systematic differences between methods by comparing regional volume and outline reproducibility of manual, FSL-FIRST and FreeSurfer hippocampus segmentations., Materials and Methods: This study used a dataset from ADNI (Alzheimer's Disease Neuroimaging Initiative), including 20 healthy controls, 40 patients with mild cognitive impairment (MCI), and 20 patients with Alzheimer's disease (AD). For each subject back-to-back (BTB) T1-weighted 3D MPRAGE images were acquired at time-point baseline (BL) and 12 months later (M12). Hippocampi segmentations of all methods were converted into triangulated meshes, regional volumes were extracted and regional Jaccard indices were computed between the hippocampi meshes of paired BTB scans to evaluate reproducibility. Regional volumes and Jaccard indices were modelled as a function of group (G), method (M), hemisphere (H), time-point (T), region (R) and interactions., Results: For the volume data the model selection procedure yielded the following significant main effects G, M, H, T and R and interaction effects G-R and M-R. The same model was found for the BTB scans. For all methods volumes reduces with the severity of disease. Significant fixed effects for the regional Jaccard index data were M, R and the interaction M-R. For all methods the middle region was most reproducible, independent of diagnostic group. FSL-FIRST was most and FreeSurfer least reproducible., Discussion/conclusion: A novel method to perform detailed analysis of subtle differences in hippocampus segmentation is proposed. The method showed that hippocampal segmentation reproducibility was best for FSL-FIRST and worst for Freesurfer. We also found systematic regional differences in hippocampal segmentation between different methods reinforcing the need of adopting harmonized protocols., Competing Interests: HV has received research grants from Novartis, Teva, MerckSerono and Pfizer, commercial companies, for this study. There are no patents, products in development or marketed products to declare. This does not alter our adherence to all the PLOS ONE policies on sharing data and materials.

Objective: Subtle motion of an epileptic patient examined with co-registered EEG and functional MRI (EEG-fMRI) may often lead to spurious fMRI activation patterns when true epileptic spikes are contaminated with motion artefacts. In recent years, methods relying on reference signals for correcting these subtle movements in the EEG have emerged. In this study, the performance of two reference-based devices are compared to the template-based method with regard to their ability to remove movement-related artifacts in EEG measured during scanning., Methods: Measurements were performed with a novel double layer cap consisting of 29 EEG and 29 reference electrodes, and with a current loop cap consisting of 60 electrodes and three current loop wires attached to the cap. EEG was acquired inside the scanner during resting state, as well as when the subject was performing a cued movement task. For the double layer cap recordings, newly developed artifact removal algorithms are introduced and both reference signal-based methods are compared to a template-based correction method., Results: The BCG artifacts occurring at resting state could be removed successfully by both the reference signal-based methods as well as by the template-based method. However, the reference signal-based methods were also capable of removing EEG artifacts induced by subtle movements, whereas the template-based method failed to remove these artifacts., Conclusion: Reference signal-based methods enable to correct for artifacts due to subtle movements, which are not removed by commonly used template-based removal algorithms., Significance: Sensitivity of EEG-fMRI analysis in patients with focal epilepsy is improved by avoiding erroneous detections of subtle movements as epileptic spikes in the EEG.

Purpose To investigate the diagnostic accuracy of an image-based classifier to distinguish between Alzheimer disease (AD) and behavioral variant frontotemporal dementia (bvFTD) in individual patients by using gray matter (GM) density maps computed from standard T1-weighted structural images obtained with multiple imagers and with independent training and prediction data. Materials and Methods The local institutional review board approved the study. Eighty-four patients with AD, 51 patients with bvFTD, and 94 control subjects were divided into independent training (n = 115) and prediction (n = 114) sets with identical diagnosis and imager type distributions. Training of a support vector machine (SVM) classifier used diagnostic status and GM density maps and produced voxelwise discrimination maps. Discriminant function analysis was used to estimate suitability of the extracted weights for single-subject classification in the prediction set. Receiver operating characteristic (ROC) curves and area under the ROC curve (AUC) were calculated for image-based classifiers and neuropsychological z scores. Results Training accuracy of the SVM was 85% for patients with AD versus control subjects, 72% for patients with bvFTD versus control subjects, and 79% for patients with AD versus patients with bvFTD (P ≤ .029). Single-subject diagnosis in the prediction set when using the discrimination maps yielded accuracies of 88% for patients with AD versus control subjects, 85% for patients with bvFTD versus control subjects, and 82% for patients with AD versus patients with bvFTD, with a good to excellent AUC (range, 0.81-0.95; P ≤ .001). Machine learning-based categorization of AD versus bvFTD based on GM density maps outperforms classification based on neuropsychological test results. Conclusion The SVM can be used in single-subject discrimination and can help the clinician arrive at a diagnosis. The SVM can be used to distinguish disease-specific GM patterns in patients with AD and those with bvFTD as compared with normal aging by using common T1-weighted structural MR imaging. (©) RSNA, 2015.

Electroencephalography (EEG) benefits from accurate head models. Dipole source modelling errors can be reduced from over 1cm to a few millimetres by replacing generic head geometry and conductivity with tailored ones. When adequate head geometry is available, electrical impedance tomography (EIT) can be used to infer the conductivities of head tissues. In this study, the boundary element method (BEM) is applied with three-compartment (scalp, skull and brain) subject-specific head models. The optimal injection of small currents to the head with a modular EIT current injector, and voltage measurement by an EEG amplifier is first sought by simulations. The measurement with a 64-electrode EEG layout is studied with respect to three noise sources affecting EIT: background EEG, deviations from the fitting assumption of equal scalp and brain conductivities, and smooth model geometry deviations from the true head geometry. The noise source effects were investigated depending on the positioning of the injection and extraction electrode and the number of their combinations used sequentially. The deviation from equal scalp and brain conductivities produces rather deterministic errors in the three conductivities irrespective of the current injection locations. With a realistic measurement of around 2 min and around 8 distant distinct current injection pairs, the error from the other noise sources is reduced to around 10% or less in the skull conductivity. The analysis of subsequent real measurements, however, suggests that there could be subject-specific local thinnings in the skull, which could amplify the conductivity fitting errors. With proper analysis of multiplexed sinusoidal EIT current injections, the measurements on average yielded conductivities of 340 mS/m (scalp and brain) and 6.6 mS/m (skull) at 2 Hz. From 11 to 127 Hz, the conductivities increased by 1.6% (scalp and brain) and 6.7% (skull) on the average. The proper analysis was ensured by using recombination of the current injections into virtual ones, avoiding problems in location-specific skull morphology variations. The observed large intersubject variations support the need for in vivo measurement of skull conductivity, resulting in calibrated subject-specific head models., (Copyright © 2015 Elsevier Inc. All rights reserved.)

In this paper we introduce a covariance framework for the analysis of single subject EEG and MEG data that takes into account observed temporal stationarity on small time scales and trial-to-trial variations. We formulate a model for the covariance matrix, which is a Kronecker product of three components that correspond to space, time and epochs/trials, and consider maximum likelihood estimation of the unknown parameter values. An iterative algorithm that finds approximations of the maximum likelihood estimates is proposed. Our covariance model is applicable in a variety of cases where spontaneous EEG or MEG acts as source of noise and realistic noise covariance estimates are needed, such as in evoked activity studies, or where the properties of spontaneous EEG or MEG are themselves the topic of interest, like in combined EEG-fMRI experiments in which the correlation between EEG and fMRI signals is investigated. We use a simulation study to assess the performance of the estimator and investigate the influence of different assumptions about the covariance factors on the estimated covariance matrix and on its components. We apply our method to real EEG and MEG data sets., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Spatial independent component analysis (ICA) is increasingly being used to extract resting-state networks from fMRI data. Previous studies showed that ICA also reveals independent components (ICs) related to the seizure onset zone. However, it is currently unknown how these epileptic ICs depend on the presence of interictal epileptic discharges (IEDs) in the EEG. The goal of this study was to explore the relation between ICs obtained from fMRI epochs during the occurrence of IEDs in the EEG and those without IEDs. fMRI data sets with co-registered EEG were retrospectively selected of patients from whom the location of the epileptogenic zone was confirmed by outcome of surgery (n = 8). The fMRI data were split into two epochs: one with IEDs visible in scalp EEG and one without. Spatial ICA was applied to the fMRI data of each part separately. The maps of all resulting components were compared to the resection area and the EEG-fMRI correlation pattern by computing a spatial correlation coefficient to detect the epilepsy-related component. For all patients, except one, there was a remarkable resemblance between the epilepsy-related components selected during epochs with IEDs and those without IEDs. These findings suggest that epilepsy-related ICs are not dependent on the presence of IEDs in scalp EEG. Since these epileptic ICs showed partial overlap with resting-state networks of healthy volunteers (n = 10), our study supports the need for new ways to classify epileptic ICs.

Anti-epileptic drugs (AEDs) have a global effect on the neurophysiology of the brain which is most likely reflected in functional brain activity recorded with EEG and fMRI. These effects may cause substantial inter-subject variability in studies where EEG correlated functional MRI (EEG-fMRI) is used to determine the epileptogenic zone in patients who are candidate for epilepsy surgery. In the present study the effects on resting state fMRI are quantified in conditions with AED administration and after withdrawal of AEDs. EEG-fMRI data were obtained from 10 patients in the condition that the patient was on the steady-state maintenance doses of AEDs as prescribed (condition A) and after withdrawal of AEDs (condition B), at the end of a clinically standard pre-surgical long term video-EEG monitoring session. Resting state networks (RSN) were extracted from fMRI. The epileptic component (ICE) was identified by selecting the RSN component with the largest overlap with the EEG-fMRI correlation pattern. Changes in RSN functional connectivity between conditions A and B were quantified. EEG-fMRI correlation analysis was successful in 30% and 100% of the cases in conditions A and B, respectively. Spatial patterns of ICEs are comparable in conditions A and B, except for one patient for whom it was not possible to identify the ICE in condition A. However, the resting state functional connectivity is significantly increased in the condition after withdrawal of AEDs (condition B), which makes resting state fMRI potentially a new tool to study AED effects. The difference in sensitivity of EEG-fMRI in conditions A and B, which is not related to the number of epileptic EEG events occurring during scanning, could be related to the increased functional connectivity in condition B.

Background: Cognitive dysfunction in multiple sclerosis (MS) has a large impact on the quality of life and is poorly understood., Objective: The aim of this study was to investigate functional network integrity in MS, and relate this to cognitive dysfunction and physical disability., Methods: Resting state fMRI scans were included of 128 MS patients and 50 controls. Eigenvector centrality mapping (ECM) was applied, a graph analysis technique that ranks the importance of brain regions based on their connectivity patterns. Significant ECM changes were related to physical disability and cognitive dysfunction., Results: In MS patients, ECM values were increased in bilateral thalamus and posterior cingulate (PCC) areas, and decreased in sensorimotor and ventral stream areas. Sensorimotor ECM decreases were related to higher EDSS (rho = -0.24, p = 0.007), while ventral stream decreases were related to poorer average cognition (rho = 0.23, p = 0.009). The thalamus displayed increased connectivity to sensorimotor and ventral stream areas., Conclusion: In MS, areas in the ventral stream and sensorimotor cortex appear to become less central in the entire functional network of the brain, which is associated with clinico-cognitive dysfunction. The thalamus, however, displays increased connectivity with these areas. These findings may aid in further elucidating the function of functional reorganization processes in MS., (© The Author(s) 2013.)

fMRI signals during rest are strongly correlated with heart rate variations. These heart rate/fMRI associations may influence the results of brain activation studies, particularly if heart rate is affected by the task. To assess the contribution of task-related heart rate changes on fMRI brain activation related to executive processing, we co-registered the electrocardiogram with fMRI in 91 subjects during an interference task (color-word Stroop) and during a planning task (Tower of London; ToL). We found that both Stroop interference and ToL planning significantly increased heart rate in the scanner and confirmed significant main effects of heart rate regressors on the fMRI signals. Nevertheless, statistical contrasts that test for increased fMRI during Stroop interference and ToL planning were not significantly influenced by inclusion of heart rate regressors. We conclude therefore that fMRI changes associated with heart rate changes do not impact strongly on higher-order fMRI effects in these commonly used executive function tasks, but routinely adding a correction seems prudent., (Copyright © 2014 Elsevier B.V. All rights reserved.)

Recent imaging studies have demonstrated functional brain network changes in patients with Alzheimer's disease (AD). Eigenvector centrality (EC) is a graph analytical measure that identifies prominent regions in the brain network hierarchy and detects localized differences between patient populations. This study used voxel-wise EC mapping (ECM) to analyze individual whole-brain resting-state functional magnetic resonance imaging (MRI) scans in 39 AD patients (age 67 ± 8) and 43 healthy controls (age 69 ± 7). Between-group differences were assessed by a permutation-based method. Associations of EC with biomarkers for AD pathology in cerebrospinal fluid (CSF) and Mini Mental State Examination (MMSE) scores were assessed using Spearman correlation analysis. Decreased EC was found bilaterally in the occipital cortex in AD patients compared to controls. Regions of increased EC were identified in the anterior cingulate and paracingulate gyrus. Across groups, frontal and occipital EC changes were associated with pathological concentrations of CSF biomarkers and with cognition. In controls, decreased EC values in the occipital regions were related to lower MMSE scores. Our main finding is that ECM, a hypothesis-free and computationally efficient analysis method of functional MRI (fMRI) data, identifies changes in brain network organization in AD patients that are related to cognition and underlying AD pathology. The relation between AD-like EC changes and cognitive performance suggests that resting-state fMRI measured EC is a potential marker of disease severity for AD., (Copyright © 2013 Wiley Periodicals, Inc.)

EEG-correlated functional MRI (EEG-fMRI) visualizes brain regions associated with interictal epileptiform discharges (IEDs). This technique images the epileptiform network, including multifocal, superficial and deeply situated cortical areas. To understand the role of EEG-fMRI in presurgical evaluation, its results should be validated relative to a gold standard. For that purpose, EEG-fMRI data were acquired for a heterogeneous group of surgical candidates (n=16) who were later implanted with subdural grids and strips (ECoG). The EEG-fMRI correlation patterns were systematically compared with brain areas involved in IEDs ECoG, using a semi-automatic analysis method, as well as to the seizure onset zone, resected area, and degree of seizure freedom. In each patient at least one of the EEG-fMRI areas was concordant with an interictally active ECoG area, always including the early onset area of IEDs in the ECoG data. This confirms that EEG-fMRI reflects a pattern of onset and propagation of epileptic activity. At group level, 76% of the BOLD regions that were covered with subdural grids, were concordant with interictally active ECoG electrodes. Due to limited spatial sampling, 51% of the BOLD regions were not covered with electrodes and could, therefore, not be validated. From an ECoG perspective it appeared that 29% of the interictally active ECoG regions were missed by EEG-fMRI and that 68% of the brain regions were correctly identified as inactive with EEG-fMRI. Furthermore, EEG-fMRI areas included the complete seizure onset zone in 83% and resected area in 93% of the data sets. No clear distinction was found between patients with a good or poor surgical outcome: in both patient groups, EEG-fMRI correlation patterns were found that were either focal or widespread. In conclusion, by comparison of EEG-fMRI with interictal invasive EEG over a relatively large patient population we were able to show that the EEG-fMRI correlation patterns are spatially accurate at the level of neurosurgical units (i.e. anatomical brain regions) and reflect the underlying network of IEDs. Therefore, we expect that EEG-fMRI can play an important role for the determination of the implantation strategy., (Copyright © 2013 Elsevier Inc. All rights reserved.)

Objective: To evaluate the viability of MEG source reconstruction in the presence of large interference due to orthodontic material., Methods: We recorded the magnetic fields following a simple hand movement and following electrical stimulation of the median nerve (somatosensory evoked field -SEF). These two tasks were performed twice, once with and once without artificial dental artefacts. Temporal Signal Space Separation (tSSS) was applied to spatially filter the data and source reconstruction was performed according to standard procedures for pre-surgical mapping of eloquent cortex, applying dipole fitting to the SEF data and beamforming to the hand movement data., Results: Comparing the data with braces to the data without braces, the observed distances between the activations following hand movement in the two conditions were on average 6.4 and 4.5 mm for the left and right hand, respectively, whereas the dipole localisation errors for the SEF were 4.1 and 5.4 mm, respectively. Without tSSS it was generally not possible to obtain reliable dipole fit or beamforming results when wearing braces., Conclusion: We confirm that tSSS is a required and effective pre-processing step for data recorded with the Elekta-MEG system. Moreover, we have shown that even the presence of large interference from orthodontic material does not significantly alter the results from dipole localisation or beamformer analysis, provided the data are spatially filtered by tSSS., Significance: State-of-the-art signal processing techniques enable the use of MEG for pre-surgical evaluation in a much larger clinical population than previously thought possible., (Copyright © 2012 International Federation of Clinical Neurophysiology. Published by Elsevier Ireland Ltd. All rights reserved.)

Co-registered EEG and functional MRI (EEG/fMRI) is a potential clinical tool for planning invasive EEG in patients with epilepsy. In addition, the analysis of EEG/fMRI data provides a fundamental insight into the precise physiological meaning of both fMRI and EEG data. Routine application of EEG/fMRI for localization of epileptic sources is hampered by large artefacts in the EEG, caused by switching of scanner gradients and heartbeat effects. Residuals of the ballistocardiogram (BCG) artefacts are similarly shaped as epileptic spikes, and may therefore cause false identification of spikes. In this study, new ideas and methods are presented to remove gradient artefacts and to reduce BCG artefacts of different shapes that mutually overlap in time. Gradient artefacts can be removed efficiently by subtracting an average artefact template when the EEG sampling frequency and EEG low-pass filtering are sufficient in relation to MR gradient switching (Gonçalves et al., 2007). When this is not the case, the gradient artefacts repeat themselves at time intervals that depend on the remainder between the fMRI repetition time and the closest multiple of the EEG acquisition time. These repetitions are deterministic, but difficult to predict due to the limited precision by which these timings are known. Therefore, we propose to estimate gradient artefact repetitions using a clustering algorithm, combined with selective averaging. Clustering of the gradient artefacts yields cleaner EEG for data recorded during scanning of a 3T scanner when using a sampling frequency of 2048 Hz. It even gives clean EEG when the EEG is sampled with only 256 Hz. Current BCG artefacts-reduction algorithms based on average template subtraction have the intrinsic limitation that they fail to deal properly with artefacts that overlap in time. To eliminate this constraint, the precise timings of artefact overlaps were modelled and represented in a sparse matrix. Next, the artefacts were disentangled with a least squares procedure. The relevance of this approach is illustrated by determining the BCG artefacts in a data set consisting of 29 healthy subjects recorded in a 1.5 T scanner and 15 patients with epilepsy recorded in a 3 T scanner. Analysis of the relationship between artefact amplitude, duration and heartbeat interval shows that in 22% (1.5T data) to 30% (3T data) of the cases BCG artefacts show an overlap. The BCG artefacts of the EEG/fMRI data recorded on the 1.5T scanner show a small negative correlation between HBI and BCG amplitude. In conclusion, the proposed methodology provides a substantial improvement of the quality of the EEG signal without excessive computer power or additional hardware than standard EEG-compatible equipment., (Copyright © 2012 Elsevier Inc. All rights reserved.)

EEG-correlated functional MRI (EEG-fMRI) has been used to indicate brain regions associated with interictal epileptiform discharges (IEDs). This technique enables the delineation of the complete epileptiform network, including multifocal and deeply situated cortical areas. Before EEG-fMRI can be used as an additional diagnostic tool in the preoperative work-up, its added value should be assessed in relation to intracranial EEG recorded from depth electrodes (SEEG) or from the cortex (ECoG), currently the clinical standard. In this study, we propose a framework for the analysis of the SEEG data to investigate in a quantitative way whether EEG-fMRI reflects the same cortical areas as identified by the IEDs present in SEEG recordings. For that purpose, the data of both modalities were analyzed with a general linear model at the same time scale and within the same spatial domain. The IEDs were used as predictors in the model, yielding for EEG-fMRI the brain voxels that were related to the IEDs and, similarly for SEEG, the electrodes that were involved. Finally, the results of the regression analysis were projected on the anatomical MRI of the patients. To explore the usefulness of this quantitative approach, a sample of five patients was studied who both underwent EEG-fMRI and SEEG recordings. For clinical validation, the results of the SEEG analysis were compared to the standard visual review of IEDs in SEEG and to the identified seizure onset zone, the resected area, and outcome of surgery. SEEG analysis revealed a spatial pattern for the most frequent and dominant IEDs present in the data of all patients. The electrodes with the highest correlation values were in good concordance with the electrodes that showed maximal amplitude during those events in the SEEG recordings. These results indicate that the analysis of SEEG data at the time scale of EEG-fMRI, using the same type of regression model, is a promising way to validate EEG-fMRI data. In fact, the BOLD areas with a positive hemodynamic response function were closely related to the spatial pattern of IEDs in the SEEG recordings in four of the five patients. The areas of significant BOLD that were not located in the vicinity of depth electrodes, were mainly characterized by negative hemodynamic responses. Furthermore, the area with a positive hemodynamic response function overlapped with the resected area in three patients, while it was located at the edge of the resection area for one. To conclude, the results of this study encourage the application of EEG-fMRI to guide the implantation of depth electrodes as prerequisite for successful epilepsy surgery., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Motor dominance is well established, but sensory dominance is much less clear. We therefore studied the cortical evoked magnetic fields using magnetoencephalography (MEG) in a group of 20 healthy right handed subjects in order to examine whether standard electrical stimulation of the median and ulnar nerve demonstrated sensory lateralization. The global field power (GFP) curves, as an indication of cortical activation, did not depict sensory lateralization to the dominant left hemisphere. Comparison of the M20, M30, and M70 peak latencies and GFP values exhibited no statistical differences between the hemispheres, indicating no sensory hemispherical dominance at these latencies for each nerve. Field maps at these latencies presented a first and second polarity reversal for both median and ulnar stimulation. Spatial dipole position parameters did not reveal statistical left-right differences at the M20, M30 and M70 peaks for both nerves. Neither did the dipolar strengths at M20, M30 and M70 show a statistical left-right difference for both nerves. Finally, the Laterality Indices of the M20, M30 and M70 strengths did not indicate complete lateralization to one of the hemispheres. After electrical median and ulnar nerve stimulation no evidence was found for sensory hand dominance in brain responses of either hand, as measured by MEG. The results can provide a new assessment of patients with sensory dysfunctions or perceptual distortion when sensory dominance occurs way beyond the estimated norm.

The analysis of simultaneous EEG and fMRI data is generally based on the extraction of regressors of interest from the EEG, which are correlated to the fMRI data in a general linear model setting. In more advanced approaches, the spatial information of EEG is also exploited by assuming underlying dipole models. In this study, we present a semi automatic and efficient method to determine electrode positions from electrode gel artifacts, facilitating the integration of EEG and fMRI in future EEG/fMRI data models. In order to visualize all electrode artifacts simultaneously in a single view, a surface rendering of the structural MRI is made using a skin triangular mesh model as reference surface, which is expanded to a "pancake view". Then the electrodes are determined with a simple mouse click for each electrode. Using the geometry of the skin surface and its transformation to the pancake view, the 3D coordinates of the electrodes are reconstructed in the MRI coordinate frame. The electrode labels are attached to the electrode positions by fitting a template grid of the electrode cap in which the labels are known. The correspondence problem between template and sample electrodes is solved by minimizing a cost function over rotations, shifts and scalings of the template grid. The crucial step here is to use the solution of the so-called "Hungarian algorithm" as a cost function, which makes it possible to identify the electrode artifacts in arbitrary order. The template electrode grid has to be constructed only once for each cap configuration. In our implementation of this method, the whole procedure can be performed within 15 min including import of MRI, surface reconstruction and transformation, electrode identification and fitting to template. The method is robust in the sense that an electrode template created for one subject can be used without identification errors for another subject for whom the same EEG cap was used. Furthermore, the method appears to be robust against spurious or missing artifacts. We therefore consider the proposed method as a useful and reliable tool within the larger toolbox required for the analysis of co-registered EEG/fMRI data., (Copyright © 2011 Elsevier Inc. All rights reserved.)

Eigenvector centrality mapping (ECM) has recently emerged as a measure to spatially characterize connectivity in functional brain imaging by attributing network properties to voxels. The main obstacle for widespread use of ECM in functional magnetic resonance imaging (fMRI) is the cost of computing and storing the connectivity matrix. This article presents fast ECM (fECM), an efficient algorithm to estimate voxel-wise eigenvector centralities from fMRI time series. Instead of explicitly storing the connectivity matrix, fECM computes matrix-vector products directly from the data, achieving high accelerations for computing voxel-wise centralities in fMRI at standard resolutions for multivariate analyses, and enabling high-resolution analyses performed on standard hardware. We demonstrate the validity of fECM at cluster and voxel levels, using synthetic and in vivo data. Results from synthetic data are compared to the theoretical gold standard, and local centrality changes in fMRI data are measured after experimental intervention. A simple scheme is presented to generate time series with prescribed covariances that represent a connectivity matrix. These time series are used to construct a 4D dataset whose volumes consist of separate regions with known intra- and inter-regional connectivities. The fECM method is tested and validated on these synthetic data. Resting-state fMRI data acquired after real-versus-sham repetitive transcranial magnetic stimulation show fECM connectivity changes in resting-state network regions. A comparison of analyses with and without accounting for motion parameters demonstrates a moderate effect of these parameters on the centrality estimates. Its computational speed and statistical sensitivity make fECM a good candidate for connectivity analyses of multimodality and high-resolution functional neuroimaging data.

Objective: Synchronization between distributed rhythms in the brain is commonly assessed by estimating the synchronization strength from simultaneous measurements. This approach, however, does not elucidate the phase dynamics that underlies synchronization. For this, an explicit dynamical model is required. Based on the assumption that the recorded rhythms can be described as weakly coupled oscillators, we propose a method for characterizing their phase-interaction dynamics., Methods: We propose to model ongoing magnetoencephalographic (MEG) oscillations as weakly coupled oscillators. Based on this model, the phase interactions between simultaneously recorded signals are characterized by estimating the modulation in instantaneous frequency as a function of their phase difference. Furthermore, we mathematically derive the effect of volume conduction on the model and show how indices for strength and direction of coupling can be derived., Results: The methodology is tested using simulations and is applied to ongoing occipital-frontal MEG oscillations of healthy subjects in the alpha and beta bands during rest. The simulations show that the model is robust against the presence of noise, short observation times, and model violations. The application to MEG data shows that the model can reconstruct the observed occipital-frontal phase difference distributions. Furthermore, it suggests that phase locking in the alpha and beta band is established by qualitatively different mechanisms., Conclusion: When the recorded rhythms are assumed to be weakly coupled oscillators, a dynamical model for the phase interactions can be fitted to data. The model is able to reconstruct the observed phase difference distribution, and hence, provides a dynamical explanation for observed phase locking., (Copyright © 2011 Wiley-Liss, Inc.)