Diffusion-based spatial priors for imaging

Multiple Source Detection based on Spatial Clustering and Its Applications on Wearable OPM-MEG

OBJECTIVE

Magnetoencephalography (MEG) is a non-invasive technique that measures the magnetic fields of brain activity. In particular, a new type of optically pumped magnetometer (OPM)-based wearable MEG system has been developed in recent years. Source localization in MEG can provide signals and locations of brain activity. However, conventional source localization methods face the difficulty of accurately estimating multiple sources. The present study presented a new parametric method to estimate the number of sources and localize multiple sources. In addition, we applied the proposed method to a constructed wearable OPM-MEG system.

METHODS

We used spatial clustering of the dipole spatial distribution to detect sources. The spatial distribution of dipoles was obtained by segmenting the MEG data temporally into slices and then estimating the parameters of the dipoles on each data slice using the particle swarm optimization algorithm. Spatial clustering was performed using the spatial-temporal density-based spatial clustering of applications with a noise algorithm. The performance of our approach for detecting multiple sources was compared with that of four typical benchmark algorithms using the OPM-MEG sensor configuration.

RESULTS

The simulation results showed that the proposed method had the best performance for detecting multiple sources. Moreover, the effectiveness of the method was verified by a multimodel sensory stimuli experiment on a real constructed 31-channel OPM-MEG.

CONCLUSION

Our study provides an effective method for the detection of multiple sources.

SIGNIFICANCE

With the improvement of the source localization methods, MEG may have a wider range of applications in neuroscience and clinical research.

Imaging somatosensory cortex responses measured by OPM-MEG: Variational free energy-based spatial smoothing estimation approach

In recent years, optically pumped magnetometer (OPM)-based magnetoencephalography (MEG) has shown potential for analyzing brain activity. It has a flexible sensor configuration and comparable sensitivity to conventional SQUID-MEG. We constructed a 32-channel OPM-MEG system and used it to measure cortical responses to median and ulnar nerve stimulations. Traditional magnetic source imaging methods tend to blur the spatial extent of sources. Accurate estimation of the spatial size of the source is important for studying the organization of brain somatotopy and for pre-surgical functional mapping. We proposed a new method called variational free energy-based spatial smoothing estimation (FESSE) to enhance the accuracy of mapping somatosensory cortex responses. A series of computer simulations based on the OPM-MEG showed better performance than the three types of competing methods under different levels of signal-to-noise ratios, source patch sizes, and co-registration errors. FESSE was then applied to the source imaging of the OPM-MEG experimental data.

Optimizing EEG Source Reconstruction with Concurrent fMRI-Derived Spatial Priors

Reconstructing EEG sources involves a complex pipeline, with the inverse problem being the most challenging. Multiple inversion algorithms are being continuously developed, aiming to tackle the non-uniqueness of this problem, which has been shown to be partially circumvented by including prior information in the inverse models. Despite a few efforts, there are still current and persistent controversies regarding the inversion algorithm of choice and the optimal set of spatial priors to be included in the inversion models. The use of simultaneous EEG-fMRI data is one approach to tackle this problem. The spatial resolution of fMRI makes fMRI derived spatial priors very convenient for EEG reconstruction, however, only task activation maps and resting-state networks (RSNs) have been explored so far, overlooking the recent, but already accepted, notion that brain networks exhibit dynamic functional connectivity fluctuations. The lack of a systematic comparison between different source reconstruction algorithms, considering potentially more brain-informative priors such as fMRI, motivates the search for better reconstruction models. Using simultaneous EEG-fMRI data, here we compared four different inversion algorithms (minimum norm, MN; low resolution electromagnetic tomography, LORETA; empirical Bayes beamformer, EBB; and multiple sparse priors, MSP) under a Bayesian framework (as implemented in SPM), each with three different sets of priors consisting of: (1) those specific to the algorithm; (2) those specific to the algorithm plus fMRI task activation maps and RSNs; and (3) those specific to the algorithm plus fMRI task activation maps and RSNs and network modules of task-related dFC states estimated from the dFC fluctuations. The quality of the reconstructed EEG sources was quantified in terms of model-based metrics, namely the expectation of the posterior probability P(model|data) and variance explained of the inversion models, and the overlap/proportion of brain regions known to be involved in the visual perception tasks that the participants were submitted to, and RSN templates, with/within EEG source components. Model-based metrics suggested that model parsimony is preferred, with the combination MSP and priors specific to this algorithm exhibiting the best performance. However, optimal overlap/proportion values were found using EBB and priors specific to this algorithm and fMRI task activation maps and RSNs or MSP and considering all the priors (algorithm priors, fMRI task activation maps and RSNs and dFC state modules), respectively, indicating that fMRI spatial priors, including dFC state modules, might contain useful information to recover EEG source components reflecting neuronal activity of interest. Our main results show that providing fMRI spatial derived priors that reflect the dynamics of the brain might be useful to map neuronal activity more accurately from EEG-fMRI. Furthermore, this work paves the way towards a more informative selection of the optimal EEG source reconstruction approach, which may be critical in future studies.

Anterior cingulate gyrus acts as a moderator of the relationship between problematic mobile phone use and depressive symptoms in college students

This study aimed to investigate the brain grey matter volume (GMV) related to problematic mobile phone use (PMPU), and whether these regions of GMV play a potential moderating role in the relationship between PMPU and depressive symptoms. We recruited 266 students who underwent magnetic resonance imaging (MRI) scanning. PMPU and depressive symptoms were assessed by a self-rating questionnaire for adolescent PMPU and patient health questionnaire-9, respectively. A multiple regression model was performed to detect GMV and white matter (WM) integrity associated with PMPU by voxel-based morphometry (VBM) and tract-based spatial statistics (TBSS) methods, and the moderating analysis was conducted by PROCESS using SPSS software. VBM analysis found an inverse correlation between the GMV of the anterior cingulate gyrus (ACC) and right fusiform gyrus (FFG) with PMPU (P_FDR < 0.05), and TBSS analysis revealed that fractional anisotropy (FA) in the body of the corpus callosum was negatively correlated with PMPU. The correlation between PMPU and depressive symptoms was moderated by the GMV of the ACC. These results suggest that the GMV of the ACC and right FFG, as well as FA in the body of the corpus callosum, was related to PMPU, and we further found that increased GMV of the ACC could reduce the relationship between PMPU and depressive symptoms in college students.

Comparison of Machine Learning Techniques Based Brain Source Localization Using EEG Signals

BACKGROUND

The brain is the most complex organ of the human body with millions of connections and activations. The electromagnetic signals are generated inside the brain due to a mental or physical task performed. These signals excite a bunch of neurons within a particular lobe depending upon the nature of the task performed. To localize this activity, certain machine learning (ML) techniques in conjunction with a neuroimaging technique (M/EEG, fMRI, PET) are developed. Different ML techniques are provided in the literature for brain source localization. Among them, the most common are: minimum norm estimation (MNE), low resolution brain electromagnetic tomography (LORETA) and Bayesian framework based multiple sparse priors (MSP).

AIMS

In this research work, EEG is used as a neuroimaging technique.

METHODS

EEG data is synthetically generated at SNR=5dB. Afterwards, ML techniques are applied to estimate the active sources. Each dataset is run for multiple trials (>40). The performance is analyzed using free energy and localization error as performance indicators. Furthermore, MSP is applied with a variant number of patches to observe the impact of patches on source localization.

RESULTS

It is observed that with an increased number of patches, the sources are localized with more precision and accuracy as expressed in terms of free energy and localization error, respectively.

CONCLUSION

The patches optimization within the Bayesian Framework produces improved results in terms of free energy and localization error.

Gaussian process linking functions for mind, brain, and behavior

The link between mind, brain, and behavior has mystified philosophers and scientists for millennia. Recent progress has been made by forming statistical associations between manifest variables of the brain (e.g., electroencephalogram [EEG], functional MRI [fMRI]) and manifest variables of behavior (e.g., response times, accuracy) through hierarchical latent variable models. Within this framework, one can make inferences about the mind in a statistically principled way, such that complex patterns of brain-behavior associations drive the inference procedure. However, previous approaches were limited in the flexibility of the linking function, which has proved prohibitive for understanding the complex dynamics exhibited by the brain. In this article, we propose a data-driven, nonparametric approach that allows complex linking functions to emerge from fitting a hierarchical latent representation of the mind to multivariate, multimodal data. Furthermore, to enforce biological plausibility, we impose both spatial and temporal structure so that the types of realizable system dynamics are constrained. To illustrate the benefits of our approach, we investigate the model's performance in a simulation study and apply it to experimental data. In the simulation study, we verify that the model can be accurately fitted to simulated data, and latent dynamics can be well recovered. In an experimental application, we simultaneously fit the model to fMRI and behavioral data from a continuous motion tracking task. We show that the model accurately recovers both neural and behavioral data and reveals interesting latent cognitive dynamics, the topology of which can be contrasted with several aspects of the experiment.

Cortical atrophy mediates the accumulating effects of vascular risk factors on cognitive decline in the Alzheimer's disease spectrum

There are increasing concerns regarding the association of vascular risk factors (VRFs) and cognitive decline in the Alzheimer's disease (AD) spectrum. Currently, we investigated whether the accumulating effects of VRFs influenced gray matter volumes and subsequently led to cognitive decline in the AD spectrum. Mediation analysis was used to explore the association among VRFs, cortical atrophy, and cognition in the AD spectrum. 123 AD spectrum were recruited and VRF scores were constructed. Multivariate linear regression analysis revealed that higher VRF scores were correlated with lower Mini-Mental State Examination scores and higher Alzheimer's Disease Assessment Scale-Cognitive Subscale scores, indicating higher VRF scores lead to severer cognitive decline in the AD spectrum. In addition, subjects with higher VRF scores suffered severe cortical atrophy, especially in medial prefrontal cortex and medial temporal lobe. More importantly, common circuits of VRFs- and cognitive decline associated with gray matter atrophy were identified. Further, using mediation analysis, we demonstrated that cortical atrophy regions significantly mediated the relationship between VRF scores and cognitive decline in the AD spectrum. These findings highlight the importance of accumulating risk in the vascular contribution to AD spectrum, and targeting VRFs may provide new strategies for the therapeutic and prevention of AD.

Altered brain structural and cognitive impairment in end-stage renal disease patients with secondary hyperparathyroidism

BACKGROUND

Cognitive impairment has received attention as an important problem in patients with end-stage renal disease, although end-stage renal disease patients with secondary hyperparathyroidism have not been studied.

PURPOSE

To assess the pattern of brain volume changes in end-stage renal disease patients with secondary hyperparathyroidism by using voxel-based morphometry and correlating these measures with clinical markers and the Montreal Cognitive Assessment scores.

MATERIAL AND METHODS

Fifty end-stage renal disease patients with no anatomical abnormalities in conventional MRI (25 patients with secondary hyperparathyroidism, 14 men, mean age 42.20 ± 7.53 years; 25 patients without secondary hyperparathyroidism, 15 men, mean age 41.96 ± 6.17 years) were selected in this study. All patients underwent laboratory tests, neuropsychological tests, and MRI. Voxel-based morphometry analysis was performed to detect regional gray matter volume differences between the two groups. The relationships between abnormal gray matter volume and clinical markers and Montreal Cognitive Assessment scores were investigated.

RESULTS

Voxel-based morphometry revealed increased gray matter volume in end-stage renal disease patients with secondary hyperparathyroidism in the bilateral caudate and bilateral thalamus compared with non- secondary hyperparathyroidism end-stage renal disease patients (P < 0.05, FWE corrected). Regarding the laboratory and neuropsychological tests, we found significant correlations between volume in these brain regions and intact parathyroid hormone levels and negative correlations with the Montreal Cognitive Assessment scores. There were no significant associations between brain volume changes and other clinical data (disease duration, urea, creatinine, and uric acid levels).

CONCLUSION

Our results showed significantly increased gray matter volume in end-stage renal disease patients with secondary hyperparathyroidism, which was associated with intact parathyroid hormone levels and cognitive impairment. Serum intact parathyroid hormone levels may be a risk factor for cognitive impairment in end-stage renal disease patients with secondary hyperparathyroidism.

Influence of Patient-Specific Head Modeling on EEG Source Imaging

Electromagnetic source imaging (ESI) techniques have become one of the most common alternatives for understanding cognitive processes in the human brain and for guiding possible therapies for neurological diseases. However, ESI accuracy strongly depends on the forward model capabilities to accurately describe the subject's head anatomy from the available structural data. Attempting to improve the ESI performance, we enhance the brain structure model within the individual-defined forward problem formulation, combining the head geometry complexity of the modeled tissue compartments and the prior knowledge of the brain tissue morphology. We validate the proposed methodology using 25 subjects, from which a set of magnetic-resonance imaging scans is acquired, extracting the anatomical priors and an electroencephalography signal set needed for validating the ESI scenarios. Obtained results confirm that incorporating patient-specific head models enhances the performed accuracy and improves the localization of focal and deep sources.

Gray Matter Structural Network Disruptions in Survivors of Acute Lymphoblastic Leukemia with Chemotherapy Treatment

OBJECTIVES

Neuroimaging studies of acute lymphoblastic leukemia (ALL) during chemotherapy treatment have shown alterations in structure, function, and connectivity in several brain regions, suggesting neurobiological impairment that might influence the large-scale brain network. This study aimed to detect the alterations in the topological organization of structural covariance networks of ALL patients.

METHODS

This study included 28 ALL patients undergoing chemotherapy and 20 matched healthy controls. We calculated the gray matter volume of 90 brain regions based on an automated anatomical labeling template and applied graph theoretical analysis to compare the topological parameters of the gray matter structural networks between the two groups.

RESULTS

The results demonstrated that both the ALL and healthy control groups exhibited a small-world topology across the range of densities. Compared to healthy controls, ALL patients had less highly interactive nodes and a reduced degree/betweenness in temporal regions, which may contribute to impaired memory and executive functions in these patients.

CONCLUSION

These results reveal that ALL patients undergoing chemotherapy treatment may have decreased regional connectivity and reduced efficiency of their structural covariance network. This is the first report of anomalous large-scale gray matter structural networks in ALL patients undergoing chemotherapy treatment and provides new insights regarding the neurobiological mechanisms underlying the chemo-brain network.

Bayesian Electromagnetic Spatio-Temporal Imaging of Extended Sources Based on Matrix Factorization

Accurate estimation of the locations and extents of neural sources from electroencephalography and magnetoencephalography (E/MEG) is challenging, especially for deep and highly correlated neural activities. In this study, we proposed a new fully data-driven source imaging method, source imaging based on spatio-temporal basis function (SI-STBF), which is built upon a Bayesian framework, to address this issue. The SI-STBF is based on the factorization of a source matrix as a product of a sparse coding matrix and a temporal basis function (TBF) matrix, which includes a few TBFs. The prior of the TBF is set in the empirical Bayesian manner. Similarly, for the spatial constraint, the SI-STBF assumes the prior covariance of the coding matrix as a weighted sum of several spatial covariance components. Both the TBFs and the coding matrix are learned from E/MEG simultaneously through variational Bayesian inference. To enable inference on high-resolution source space, we derived a scalable algorithm using convex analysis. The performance of the SI-STBF was assessed using both simulated and experimental E/MEG recordings. Compared with L₂-norm constrained methods, the SI-STBF is superior in reconstructing extended sources with less spatial diffusion and less localization error. By virtue of the spatio-temporal factorization of source matrix, the SI-STBF also produces more accurate estimations than spatial-only constraint method for high correlated and deep sources.

Relationship between extraversion personality and gray matter volume and functional connectivity density in healthy young adults: an fMRI study

Extraversion and neuroticism are two main dimensions of Eysenck's personality. We assessed the relationship between extraversion and neuroticism with brain structure and function by voxel-based morphometry (VBM) and functional connectivity density (FCD). The resting state functional magnetic resonance image and high resolution structural T1 weighted images of 100 young healthy subjects were used in analysis. Our results showed that extraversion was negatively correlated with gray matter volume (GMV) of the bilateral putamen, and it was negatively correlated with FCD in the precuneus. No associations between neuroticism and brain structure and function changes. Overall, our results suggested that several brain regions involved in shaping of extraversion traits among young individuals, which may provide a neurobiological basis of extraversion.

Quantifying the performance of MEG source reconstruction using resting state data

In magnetoencephalography (MEG) research there are a variety of inversion methods to transform sensor data into estimates of brain activity. Each new inversion scheme is generally justified against a specific simulated or task scenario. The choice of this scenario will however have a large impact on how well the scheme performs. We describe a method with minimal selection bias to quantify algorithm performance using human resting state data. These recordings provide a generic, heterogeneous, and plentiful functional substrate against which to test different MEG recording and reconstruction approaches. We used a Hidden Markov model to spatio-temporally partition data into self-similar dynamic states. To test the anatomical precision that could be achieved, we then inverted these data onto libraries of systematically distorted subject-specific cortical meshes and compared the quality of the fit using cross validation and a Free energy metric. This revealed which inversion scheme was able to identify the least distorted (most accurate) anatomical models, and allowed us to quantify an upper bound on the mean anatomical distortion accordingly. We used two resting state datasets, one recorded with head-casts and one without. In the head-cast data, the Empirical Bayesian Beamformer (EBB) algorithm showed the best mean anatomical discrimination (3.7 mm) compared with Minimum Norm/LORETA (6.0 mm) and Multiple Sparse Priors (9.4 mm). This pattern was replicated in the second (conventional dataset) although with a marginally poorer (non-significant) prediction of the missing (cross-validated) data. Our findings suggest that the abundant resting state data now commonly available could be used to refine and validate MEG source reconstruction methods and/or recording paradigms.

Reproducibility of EEG-MEG fusion source analysis of interictal spikes: Relevance in presurgical evaluation of epilepsy

Fusion of electroencephalography (EEG) and magnetoencephalography (MEG) data using maximum entropy on the mean method (MEM-fusion) takes advantage of the complementarities between EEG and MEG to improve localization accuracy. Simulation studies demonstrated MEM-fusion to be robust especially in noisy conditions such as single spike source localizations (SSSL). Our objective was to assess the reliability of SSSL using MEM-fusion on clinical data. We proposed to cluster SSSL results to find the most reliable and consistent source map from the reconstructed sources, the so-called consensus map. Thirty-four types of interictal epileptic discharges (IEDs) were analyzed from 26 patients with well-defined epileptogenic focus. SSSLs were performed on EEG, MEG, and fusion data and consensus maps were estimated using hierarchical clustering. Qualitative (spike-to-spike reproducibility rate, SSR) and quantitative (localization error and spatial dispersion) assessments were performed using the epileptogenic focus as clinical reference. Fusion SSSL provided significantly better results than EEG or MEG alone. Fusion found at least one cluster concordant with the clinical reference in all cases. This concordant cluster was always the one involving the highest number of spikes. Fusion yielded highest reproducibility (SSR EEG = 55%, MEG = 71%, fusion = 90%) and lowest localization error. Also, using only few channels from either modality (21EEG + 272MEG or 54EEG + 25MEG) was sufficient to reach accurate fusion. MEM-fusion with consensus map approach provides an objective way of finding the most reliable and concordant generators of IEDs. We, therefore, suggest the pertinence of SSSL using MEM-fusion as a valuable clinical tool for presurgical evaluation of epilepsy.

Serial correlations in single-subject fMRI with sub-second TR

When performing statistical analysis of single-subject fMRI data, serial correlations need to be taken into account to allow for valid inference. Otherwise, the variability in the parameter estimates might be under-estimated resulting in increased false-positive rates. Serial correlations in fMRI data are commonly characterized in terms of a first-order autoregressive (AR) process and then removed via pre-whitening. The required noise model for the pre-whitening depends on a number of parameters, particularly the repetition time (TR). Here we investigate how the sub-second temporal resolution provided by simultaneous multislice (SMS) imaging changes the noise structure in fMRI time series. We fit a higher-order AR model and then estimate the optimal AR model order for a sequence with a TR of less than 600 ms providing whole brain coverage. We show that physiological noise modelling successfully reduces the required AR model order, but remaining serial correlations necessitate an advanced noise model. We conclude that commonly used noise models, such as the AR(1) model, are inadequate for modelling serial correlations in fMRI using sub-second TRs. Rather, physiological noise modelling in combination with advanced pre-whitening schemes enable valid inference in single-subject analysis using fast fMRI sequences.

No structural cerebral MRI changes related to fatigue in patients with primary Sjögren's syndrome

Objective

Whether or not chronic fatigue is reflected in structural changes in the brain is a matter of debate. Primary SS (pSS) is characterized by dryness of the mouth and eyes, migrating muscle and joint pain and prominent fatigue. We aimed to investigate whether the severity of fatigue in pSS was associated with cerebral MRI findings.

Methods

Fatigue was measured with the fatigue visual analog scale in 65 patients with pSS. Global grey matter (GM) and white matter volumes were estimated from magnetic resonance T1 images, and associations between fatigue and brain volumes were assessed in regression models. Voxel-based morphometric analyses of GM were performed to investigate possible associations between fatigue and GM volume changes in particular brain regions.

Results

The fatigue scores in the patient group were spread across a wide range. Global volume analyses showed no significant effect of GM volumes and white matter volumes on fatigue. Voxel-wise analyses of GM did not identify any particular brain region associated with fatigue.

Conclusion

Fatigue is a dominant phenomenon in pSS patients but is not reflected in structural abnormalities in the brain as visualized by conventional MRI. Our findings support the hypothesis of fatigue as a physiological phenomenon that does not lead to vascular changes or neuronal or glial death or damage.

Adaptive smoothing based on Gaussian processes regression increases the sensitivity and specificity of fMRI data

Temporal and spatial filtering of fMRI data is often used to improve statistical power. However, conventional methods, such as smoothing with fixed-width Gaussian filters, remove fine-scale structure in the data, necessitating a tradeoff between sensitivity and specificity. Specifically, smoothing may increase sensitivity (reduce noise and increase statistical power) but at the cost loss of specificity in that fine-scale structure in neural activity patterns is lost. Here, we propose an alternative smoothing method based on Gaussian processes (GP) regression for single subjects fMRI experiments. This method adapts the level of smoothing on a voxel by voxel basis according to the characteristics of the local neural activity patterns. GP-based fMRI analysis has been heretofore impractical owing to computational demands. Here, we demonstrate a new implementation of GP that makes it possible to handle the massive data dimensionality of the typical fMRI experiment. We demonstrate how GP can be used as a drop-in replacement to conventional preprocessing steps for temporal and spatial smoothing in a standard fMRI pipeline. We present simulated and experimental results that show the increased sensitivity and specificity compared to conventional smoothing strategies. Hum Brain Mapp 38:1438-1459, 2017. © 2016 Wiley Periodicals, Inc.

Emotion Discrimination Using Spatially Compact Regions of Interest Extracted from Imaging EEG Activity

Lately, research on computational models of emotion had been getting much attention due to their potential for understanding the mechanisms of emotions and their promising broad range of applications that potentially bridge the gap between human and machine interactions. We propose a new method for emotion classification that relies on features extracted from those active brain areas that are most likely related to emotions. To this end, we carry out the selection of spatially compact regions of interest that are computed using the brain neural activity reconstructed from Electroencephalography data. Throughout this study, we consider three representative feature extraction methods widely applied to emotion detection tasks, including Power spectral density, Wavelet, and Hjorth parameters. Further feature selection is carried out using principal component analysis. For validation purpose, these features are used to feed a support vector machine classifier that is trained under the leave-one-out cross-validation strategy. Obtained results on real affective data show that incorporation of the proposed training method in combination with the enhanced spatial resolution provided by the source estimation allows improving the performed accuracy of discrimination in most of the considered emotions, namely: dominance, valence, and liking.

Intracranial EEG potentials estimated from MEG sources: A new approach to correlate MEG and iEEG data in epilepsy

Detection of epileptic spikes in MagnetoEncephaloGraphy (MEG) requires synchronized neuronal activity over a minimum of 4cm2. We previously validated the Maximum Entropy on the Mean (MEM) as a source localization able to recover the spatial extent of the epileptic spike generators. The purpose of this study was to evaluate quantitatively, using intracranial EEG (iEEG), the spatial extent recovered from MEG sources by estimating iEEG potentials generated by these MEG sources. We evaluated five patients with focal epilepsy who had a pre-operative MEG acquisition and iEEG with MRI-compatible electrodes. Individual MEG epileptic spikes were localized along the cortical surface segmented from a pre-operative MRI, which was co-registered with the MRI obtained with iEEG electrodes in place for identification of iEEG contacts. An iEEG forward model estimated the influence of every dipolar source of the cortical surface on each iEEG contact. This iEEG forward model was applied to MEG sources to estimate iEEG potentials that would have been generated by these sources. MEG-estimated iEEG potentials were compared with measured iEEG potentials using four source localization methods: two variants of MEM and two standard methods equivalent to minimum norm and LORETA estimates. Our results demonstrated an excellent MEG/iEEG correspondence in the presumed focus for four out of five patients. In one patient, the deep generator identified in iEEG could not be localized in MEG. MEG-estimated iEEG potentials is a promising method to evaluate which MEG sources could be retrieved and validated with iEEG data, providing accurate results especially when applied to MEM localizations. Hum Brain Mapp 37:1661-1683, 2016. © 2016 Wiley Periodicals, Inc.

Objective Bayesian fMRI analysis-a pilot study in different clinical environments

Functional MRI (fMRI) used for neurosurgical planning delineates functionally eloquent brain areas by time-series analysis of task-induced BOLD signal changes. Commonly used frequentist statistics protect against false positive results based on a p-value threshold. In surgical planning, false negative results are equally if not more harmful, potentially masking true brain activity leading to erroneous resection of eloquent regions. Bayesian statistics provides an alternative framework, categorizing areas as activated, deactivated, non-activated or with low statistical confidence. This approach has not yet found wide clinical application partly due to the lack of a method to objectively define an effect size threshold. We implemented a Bayesian analysis framework for neurosurgical planning fMRI. It entails an automated effect-size threshold selection method for posterior probability maps accounting for inter-individual BOLD response differences, which was calibrated based on the frequentist results maps thresholded by two clinical experts. We compared Bayesian and frequentist analysis of passive-motor fMRI data from 10 healthy volunteers measured on a pre-operative 3T and an intra-operative 1.5T MRI scanner. As a clinical case study, we tested passive motor task activation in a brain tumor patient at 3T under clinical conditions. With our novel effect size threshold method, the Bayesian analysis revealed regions of all four categories in the 3T data. Activated region foci and extent were consistent with the frequentist analysis results. In the lower signal-to-noise ratio 1.5T intra-operative scanner data, Bayesian analysis provided improved brain-activation detection sensitivity compared with the frequentist analysis, albeit the spatial extents of the activations were smaller than at 3T. Bayesian analysis of fMRI data using operator-independent effect size threshold selection may improve the sensitivity and certainty of information available to guide neurosurgery.

Discrimination of cortical laminae using MEG

Typically MEG source reconstruction is used to estimate the distribution of current flow on a single anatomically derived cortical surface model. In this study we use two such models representing superficial and deep cortical laminae. We establish how well we can discriminate between these two different cortical layer models based on the same MEG data in the presence of different levels of co-registration noise, Signal-to-Noise Ratio (SNR) and cortical patch size. We demonstrate that it is possible to make a distinction between superficial and deep cortical laminae for levels of co-registration noise of less than 2 mm translation and 2° rotation at SNR > 11 dB. We also show that an incorrect estimate of cortical patch size will tend to bias layer estimates. We then use a 3D printed head-cast (Troebinger et al., 2014) to achieve comparable levels of co-registration noise, in an auditory evoked response paradigm, and show that it is possible to discriminate between these cortical layer models in real data.

•

We evaluate necessary recording precision to distinguish superficial/deep laminae.

•

For coregistration error of < 2 mm/2° we can distinguish between these laminar models.

•

Incorrect assumptions about cortical patch size bias these layer estimates.

•

Initial results suggest that the auditory M100 derives from deep cortical layers.

Algorithmic procedures for Bayesian MEG/EEG source reconstruction in SPM

The MEG/EEG inverse problem is ill-posed, giving different source reconstructions depending on the initial assumption sets. Parametric Empirical Bayes allows one to implement most popular MEG/EEG inversion schemes (Minimum Norm, LORETA, etc.) within the same generic Bayesian framework. It also provides a cost-function in terms of the variational Free energy-an approximation to the marginal likelihood or evidence of the solution. In this manuscript, we revisit the algorithm for MEG/EEG source reconstruction with a view to providing a didactic and practical guide. The aim is to promote and help standardise the development and consolidation of other schemes within the same framework. We describe the implementation in the Statistical Parametric Mapping (SPM) software package, carefully explaining each of its stages with the help of a simple simulated data example. We focus on the Multiple Sparse Priors (MSP) model, which we compare with the well-known Minimum Norm and LORETA models, using the negative variational Free energy for model comparison. The manuscript is accompanied by Matlab scripts to allow the reader to test and explore the underlying algorithm.

Random location of multiple sparse priors for solving the MEG/EEG inverse problem

MEG/EEG brain imaging has become an important tool in neuroimaging. Current techniques based in Bayesian approaches require an a-priori definition of patch locations on the cortical manifold. Too many patches results in a complex optimisation problem, too few an under sampling of the solution space. In this work random locations of the possible active regions of the brain are proposed to iteratively arrive at a solution. We use Bayesian model averaging to combine different possible solutions. The proposed methodology was tested with synthetic MEG datasets reducing the localisation error of the approaches based on fixed locations. Real data from a visual attention study was used for validation.

Fully Bayesian inference for structural MRI: application to segmentation and statistical analysis of T2-hypointensities

Aiming at iron-related T2-hypointensity, which is related to normal aging and neurodegenerative processes, we here present two practicable approaches, based on Bayesian inference, for preprocessing and statistical analysis of a complex set of structural MRI data. In particular, Markov Chain Monte Carlo methods were used to simulate posterior distributions. First, we rendered a segmentation algorithm that uses outlier detection based on model checking techniques within a Bayesian mixture model. Second, we rendered an analytical tool comprising a Bayesian regression model with smoothness priors (in the form of Gaussian Markov random fields) mitigating the necessity to smooth data prior to statistical analysis. For validation, we used simulated data and MRI data of 27 healthy controls (age: ; range, ). We first observed robust segmentation of both simulated T2-hypointensities and gray-matter regions known to be T2-hypointense. Second, simulated data and images of segmented T2-hypointensity were analyzed. We found not only robust identification of simulated effects but also a biologically plausible age-related increase of T2-hypointensity primarily within the dentate nucleus but also within the globus pallidus, substantia nigra, and red nucleus. Our results indicate that fully Bayesian inference can successfully be applied for preprocessing and statistical analysis of structural MRI data.

MEG source localization of spatially extended generators of epileptic activity: comparing entropic and hierarchical bayesian approaches

Localizing the generators of epileptic activity in the brain using Electro-EncephaloGraphy (EEG) or Magneto-EncephaloGraphy (MEG) signals is of particular interest during the pre-surgical investigation of epilepsy. Epileptic discharges can be detectable from background brain activity, provided they are associated with spatially extended generators. Using realistic simulations of epileptic activity, this study evaluates the ability of distributed source localization methods to accurately estimate the location of the generators and their sensitivity to the spatial extent of such generators when using MEG data. Source localization methods based on two types of realistic models have been investigated: (i) brain activity may be modeled using cortical parcels and (ii) brain activity is assumed to be locally smooth within each parcel. A Data Driven Parcellization (DDP) method was used to segment the cortical surface into non-overlapping parcels and diffusion-based spatial priors were used to model local spatial smoothness within parcels. These models were implemented within the Maximum Entropy on the Mean (MEM) and the Hierarchical Bayesian (HB) source localization frameworks. We proposed new methods in this context and compared them with other standard ones using Monte Carlo simulations of realistic MEG data involving sources of several spatial extents and depths. Detection accuracy of each method was quantified using Receiver Operating Characteristic (ROC) analysis and localization error metrics. Our results showed that methods implemented within the MEM framework were sensitive to all spatial extents of the sources ranging from 3 cm² to 30 cm², whatever were the number and size of the parcels defining the model. To reach a similar level of accuracy within the HB framework, a model using parcels larger than the size of the sources should be considered.

EEG/fMRI fusion based on independent component analysis: integration of data-driven and model-driven methods

Simultaneous electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) provide complementary noninvasive information of brain activity, and EEG/fMRI fusion can achieve higher spatiotemporal resolution than each modality separately. This focuses on independent component analysis (ICA)-based EEG/fMRI fusion. In order to appreciate the issues, we first describe the potential and limitations of the developed fusion approaches: fMRI-constrained EEG imaging, EEG-informed fMRI analysis, and symmetric fusion. We then outline some newly developed hybrid fusion techniques using ICA and the combination of data-/model-driven methods, with special mention of the spatiotemporal EEG/fMRI fusion (STEFF). Finally, we discuss the current trend in methodological development and the existing limitations for extrapolating neural dynamics.

Modeling brain activation in fMRI using group MRF

Noise confounds present serious complications to functional magnetic resonance imaging (fMRI) analysis. The amount of discernible signals within a single dataset of a subject is often inadequate to obtain satisfactory intra-subject activation detection. To remedy this limitation, we propose a novel group Markov random field (GMRF) that extends each subject's neighborhood system to other subjects to enable information coalescing. A distinct advantage of GMRF over standard fMRI group analysis is that no stringent one-to-one voxel correspondence is required. Instead, intra- and inter-subject neighboring voxels are jointly regularized to encourage spatially proximal voxels to be assigned similar labels across subjects. Our proposed group-extended graph structure thus provides an effective means for handling inter-subject variability. Also, adopting a group-wise approach by integrating group information into intra-subject activation, as opposed to estimating a single average group map, permits inter-subject differences to be characterized and studied. GMRF can be elegantly implemented as a single MRF, thus enabling all subjects' activation maps to be simultaneously and collaboratively segmented with global optimality guaranteed in the case of binary labeling. We validate our technique on synthetic and real fMRI data and demonstrate GMRF's superior performance over standard fMRI analysis.

A general Bayesian treatment for MEG source reconstruction incorporating lead field uncertainty

There is uncertainty introduced when a cortical surface based model derived from an anatomical MRI is used to reconstruct neural activity with MEG data. This is a specific case of a problem with uncertainty in parameters on which M/EEG lead fields depend non-linearly. Here we present a general mathematical treatment of any such problem with a particular focus on co-registration. We use a Metropolis search followed by Bayesian Model Averaging over multiple sparse prior source inversions with different headlocation/orientation parameters. Based on MEG data alone we can locate the cortex to within 4mm at empirically realistic signal to noise ratios. We also show that this process gives improved posterior distributions on the estimated current distributions, and can be extended to make inference on the locations of local maxima by providing confidence intervals for each source.

Connectivity-informed fMRI activation detection

A growing interest has emerged in studying the correlation structure of spontaneous and task-induced brain activity to elucidate the functional architecture of the brain. In particular, functional networks estimated from resting state (RS) data were shown to exhibit high resemblance to those evoked by stimuli. Motivated by these findings, we propose a novel generative model that integrates RS-connectivity and stimulus-evoked responses under a unified analytical framework. Our model permits exact closed-form solutions for both the posterior activation effect estimates and the model evidence. To learn RS networks, graphical LASSO and the oracle approximating shrinkage technique are deployed. On a cohort of 65 subjects, we demonstrate increased sensitivity in fMRI activation detection using our connectivity-informed model over the standard univariate approach. Our results thus provide further evidence for the presence of an intrinsic relationship between brain activity during rest and task, the exploitation of which enables higher detection power in task-driven studies.

Optimizing the performance of local canonical correlation analysis in fMRI using spatial constraints

The benefits of locally adaptive statistical methods for fMRI research have been shown in recent years, as these methods are more proficient in detecting brain activations in a noisy environment. One such method is local canonical correlation analysis (CCA), which investigates a group of neighboring voxels instead of looking at the single voxel time course. The value of a suitable test statistic is used as a measure of activation. It is customary to assign the value to the center voxel for convenience. The method without constraints is prone to artifacts, especially in a region of localized strong activation. To compensate for these deficiencies, the impact of different spatial constraints in CCA on sensitivity and specificity are investigated. The ability of constrained CCA (cCCA) to detect activation patterns in an episodic memory task has been studied. This research shows how any arbitrary contrast of interest can be analyzed by cCCA and how accurate P-values optimized for the contrast of interest can be computed using nonparametric methods. Results indicate an increase of up to 20% in detecting activation patterns for some of the advanced cCCA methods, as measured by ROC curves derived from simulated and real fMRI data.

fMRI functional networks for EEG source imaging

The brain exhibits temporally coherent networks (TCNs) involving numerous cortical and sub-cortical regions both during the rest state and during the performance of cognitive tasks. TCNs represent the interactions between different brain areas, and understanding such networks may facilitate electroencephalography (EEG) source estimation. We propose a new method for examining TCNs using scalp EEG in conjunction with data obtained by functional magnetic resonance imaging (fMRI). In this approach, termed NEtwork based SOurce Imaging (NESOI), multiple TCNs derived from fMRI with independent component analysis (ICA) are used as the covariance priors of the EEG source reconstruction using Parametric Empirical Bayesian (PEB). In contrast to previous applications of PEB in EEG source imaging with smoothness or sparseness priors, TCNs play a fundamental role among the priors used by NESOI. NESOI achieves an efficient integration of the high temporal resolution EEG and TCN derived from the high spatial resolution fMRI. Using synthetic and real data, we directly compared the performance of NESOI with other distributed source inversion methods, with and without the use of fMRI priors. Our results indicated that NESOI is a potentially useful approach for EEG source imaging.

Incorporating FMRI functional networks in EEG source imaging: a Bayesian model comparison approach

Brain functional networks extracted from fMRI can improve the accuracy of EEG source localization. However, the coupling between EEG and fMRI remains poorly understood, i.e., whether fMRI networks provide information about the magnitude of neural activity, and whether neural sources demonstrate temporal correlations within each network. In this paper, we present an improved version of the NEtwork-based SOurce Imaging method (iNESOI) through Bayesian model comparison. Different models correspond to various matching between EEG and fMRI, and the appropriate one is selected by data with the model evidence. Synthetic and real data tests show that iNESOI has potential to select the appropriate fMRI priors to reach a better source reconstruction than some other typical approaches.

A topographic latent source model for fMRI data

We describe and evaluate a new statistical generative model of functional magnetic resonance imaging (fMRI) data. The model, topographic latent source analysis (TLSA), assumes that fMRI images are generated by a covariate-dependent superposition of latent sources. These sources are defined in terms of basis functions over space. The number of parameters in the model does not depend on the number of voxels, enabling a parsimonious description of activity patterns that avoids many of the pitfalls of traditional voxel-based approaches. We develop a multi-subject extension where latent sources at the subject-level are perturbations of a group-level template. We evaluate TLSA according to prediction, reconstruction and reproducibility. We show that it compares favorably to a Naive Bayes model while using fewer parameters. We also describe a hypothesis testing framework that can be used to identify significant latent sources.

Detection and characterization of single-trial fMRI bold responses: paradigm free mapping

This work presents a novel method of mapping the brain's response to single stimuli in space and time without prior knowledge of the paradigm timing: paradigm free mapping (PFM). This method is based on deconvolution of the hemodynamic response from the voxel time series assuming a linear response and using a ridge-regression algorithm. Statistical inference is performed by defining a spatio-temporal t-statistic and by controlling for multiple comparisons using the false discovery rate procedure. The methodology was validated on five subjects who performed self-paced and visually cued finger tapping at 7 Tesla, with moderate (TR = 2 s) and high (TR = 0.4 s) temporal resolution. The results demonstrate that detection of single-trial BOLD events is feasible without a priori information on the stimulus paradigm. The proposed method opens up the possibility of designing temporally unconstrained paradigms to study the cortical response to unpredictable mental events.

A parallel framework for simultaneous EEG/fMRI analysis: methodology and simulation

Concurrent EEG/fMRI recordings represent multiple, simultaneously active, regionally overlapping neuronal mass responses. To address the problems caused by the overlapping nature of these responses, we propose a parallel framework for Spatial-Temporal EEG/fMRI Fusion (STEFF). This technique adopts Independent Component Analysis (ICA) to recover the time-course and spatial mapping components from EEG and fMRI separately. These components are then linked concurrently in the spatial and temporal domain using an Empirical Bayesian (EB) model. This approach enables information one modality to be utilized as priors for the other and hence improves the spatial (for EEG) or temporal (for fMRI) resolution of the other modality. Consequently, STEFF achieves flexible and sparse matching among EEG and fMRI components with common neuronal substrates. Simulations under realistic noise conditions indicated that STEFF is a feasible and physiologically reasonable hybrid approach for spatiotemporal mapping of cognitive processing in the human brain.

Bayesian spatiotemporal model of fMRI data

This research describes a new Bayesian spatiotemporal model to analyse block-design BOLD fMRI studies. In the temporal dimension, we parameterise the hemodynamic response function's (HRF) shape with a potential increase of signal and a subsequent exponential decay. In the spatial dimension, we use Gaussian Markov random fields (GMRF) priors on activation characteristics parameters (location and magnitude) that embody our prior knowledge that evoked responses are spatially contiguous and locally homogeneous. The result is a spatiotemporal model with a small number of parameters, all of them interpretable. Simulations from the model are performed in order to ascertain the performance of the sampling scheme and the ability of the posterior to estimate model parameters, as well as to check the model sensitivity to signal to noise ratio. Results are shown on synthetic data and on real data from a block-design fMRI experiment.

Combined spatial and non-spatial prior for inference on MRI time-series

When modelling FMRI and other MRI time-series data, a Bayesian approach based on adaptive spatial smoothness priors is a compelling alternative to using a standard generalized linear model (GLM) on presmoothed data. Another benefit of the Bayesian approach is that biophysical prior information can be incorporated in a principled manner; however, this requirement for a fixed non-spatial prior on a parameter would normally preclude using spatial regularization on that same parameter. We have developed a Gaussian-process-based prior to apply adaptive spatial regularization while still ensuring that the fixed biophysical prior is correctly applied on each voxel. A parameterized covariance matrix provides separate control over the variance (the diagonal elements) and the between-voxel correlation (due to off-diagonal elements). Analysis proceeds using evidence optimization (EO), with variational Bayes (VB) updates used for some parameters. The method can also be applied to non-linear forward models by using a linear Taylor expansion centred on the latest parameter estimates. Applying the method to FMRI with a constrained haemodynamic response function (HRF) shape model shows improved fits in simulations, compared to using either the non-spatial or spatial-smoothness prior alone. We also analyse multi-inversion arterial spin labelling data using a non-linear perfusion model to estimate cerebral blood flow and bolus arrival time. By combining both types of prior information, this new prior performs consistently well across a wider range of situations than either prior alone, and provides better estimates when both types of prior information are relevant.

High-resolution fMRI: overcoming the signal-to-noise problem

Increasing the spatial resolution in functional Magnetic Resonance Imaging (fMRI) inherently lowers the signal-to-noise ratio (SNR). In order to still detect functionally significant activations in high-resolution images, spatial smoothing of the data is required. However, conventional non-adaptive smoothing comes with a reduced effective resolution, foiling the benefit of the higher acquisition resolution. We show how our recently proposed structural adaptive smoothing procedure for functional MRI data can improve signal detection of high-resolution fMRI experiments regardless of the lower SNR. The procedure is evaluated on human visual and sensory-motor mapping experiments. In these applications, the higher resolution could be fully utilized and high-resolution experiments were outperforming normal resolution experiments by means of both statistical significance and information content.

Inferring functional connections between neurons

A central question in neuroscience is how interactions between neurons give rise to behavior. In many electrophysiological experiments, the activity of a set of neurons is recorded while sensory stimuli or movement tasks are varied. Tools that aim to reveal underlying interactions between neurons from such data can be extremely useful. Traditionally, neuroscientists have studied these interactions using purely descriptive statistics (cross-correlograms or joint peri-stimulus time histograms). However, the interpretation of such data is often difficult, particularly as the number of recorded neurons grows. Recent research suggests that model-based, maximum likelihood methods can improve these analyses. In addition to estimating neural interactions, application of these techniques has improved decoding of external variables, created novel interpretations of existing electrophysiological data, and may provide new insight into how the brain represents information.

Diffusion-based spatial priors for functional magnetic resonance images

We recently outlined a Bayesian scheme for analyzing fMRI data using diffusion-based spatial priors [Harrison, L.M., Penny, W., Ashburner, J., Trujillo-Barreto, N., Friston, K.J., 2007. Diffusion-based spatial priors for imaging. NeuroImage 38, 677–695]. The current paper continues this theme, applying it to a single-subject functional magnetic resonance imaging (fMRI) study of the auditory system. We show that spatial priors on functional activations, based on diffusion, can be formulated in terms of the eigenmodes of a graph Laplacian. This allows one to discard eigenmodes with small eigenvalues, to provide a computationally efficient scheme. Furthermore, this formulation shows that diffusion-based priors are a generalization of conventional Laplacian priors [Penny, W.D., Trujillo-Barreto, N.J., Friston, K.J., 2005. Bayesian fMRI time series analysis with spatial priors. NeuroImage 24, 350–362]. Finally, we show how diffusion-based priors are a special case of Gaussian process models that can be inverted using classical covariance component estimation techniques like restricted maximum likelihood [Patterson, H.D., Thompson, R., 1974. Maximum likelihood estimation of components of variance. Paper presented at: 8th International Biometrics Conference (Constanta, Romania)]. The convention in SPM is to smooth data with a fixed isotropic Gaussian kernel before inverting a mass-univariate statistical model. This entails the strong assumption that data are generated smoothly throughout the brain. However, there is no way to determine if this assumption is supported by the data, because data are smoothed before statistical modeling. In contrast, if a spatial prior is used, smoothness is estimated given non-smoothed data. Explicit spatial priors enable formal model comparison of different prior assumptions, e.g., that data are generated from a stationary (i.e., fixed throughout the brain) or non-stationary spatial process. Indeed, for the auditory data we provide strong evidence for a non-stationary process, which concurs with a qualitative comparison of predicted activations at the boundary of functionally selective regions.