Brain Decoding Using fMRI Images for Multiple Subjects through Deep Learning

Substantial information related to human cerebral conditions can be decoded through various noninvasive evaluating techniques like fMRI. Exploration of the neuronal activity of the human brain can divulge the thoughts of a person like what the subject is perceiving, thinking, or visualizing. Furthermore, deep learning techniques can be used to decode the multifaceted patterns of the brain in response to external stimuli. Existing techniques are capable of exploring and classifying the thoughts of the human subject acquired by the fMRI imaging data. fMRI images are the volumetric imaging scans which are highly dimensional as well as require a lot of time for training when fed as an input in the deep learning network. However, the hassle for more efficient learning of highly dimensional high-level features in less training time and accurate interpretation of the brain voxels with less misclassification error is needed. In this research, we propose an improved CNN technique where features will be functionally aligned. The optimal features will be selected after dimensionality reduction. The highly dimensional feature vector will be transformed into low dimensional space for dimensionality reduction through autoadjusted weights and combination of best activation functions. Furthermore, we solve the problem of increased training time by using Swish activation function, making it denser and increasing efficiency of the model in less training time. Finally, the experimental results are evaluated and compared with other classifiers which demonstrated the supremacy of the proposed model in terms of accuracy.

Nonlinear visuoauditory integration in the mouse superior colliculus

Sensory information from different modalities is processed in parallel, and then integrated in associative brain areas to improve object identification and the interpretation of sensory experiences. The Superior Colliculus (SC) is a midbrain structure that plays a critical role in integrating visual, auditory, and somatosensory input to assess saliency and promote action. Although the response properties of the individual SC neurons to visuoauditory stimuli have been characterized, little is known about the spatial and temporal dynamics of the integration at the population level. Here we recorded the response properties of SC neurons to spatially restricted visual and auditory stimuli using large-scale electrophysiology. We then created a general, population-level model that explains the spatial, temporal, and intensity requirements of stimuli needed for sensory integration. We found that the mouse SC contains topographically organized visual and auditory neurons that exhibit nonlinear multisensory integration. We show that nonlinear integration depends on properties of auditory but not visual stimuli. We also find that a heuristically derived nonlinear modulation function reveals conditions required for sensory integration that are consistent with previously proposed models of sensory integration such as spatial matching and the principle of inverse effectiveness.

Estimating the dimensionality of the manifold underlying multi-electrode neural recordings

It is generally accepted that the number of neurons in a given brain area far exceeds the number of neurons needed to carry any specific function controlled by that area. For example, motor areas of the human brain contain tens of millions of neurons that control the activation of tens or at most hundreds of muscles. This massive redundancy implies the covariation of many neurons, which constrains the population activity to a low-dimensional manifold within the space of all possible patterns of neural activity. To gain a conceptual understanding of the complexity of the neural activity within a manifold, it is useful to estimate its dimensionality, which quantifies the number of degrees of freedom required to describe the observed population activity without significant information loss. While there are many algorithms for dimensionality estimation, we do not know which are well suited for analyzing neural activity. The objective of this study was to evaluate the efficacy of several representative algorithms for estimating the dimensionality of linearly and nonlinearly embedded data. We generated synthetic neural recordings with known intrinsic dimensionality and used them to test the algorithms' accuracy and robustness. We emulated some of the important challenges associated with experimental data by adding noise, altering the nature of the embedding of the low-dimensional manifold within the high-dimensional recordings, varying the dimensionality of the manifold, and limiting the amount of available data. We demonstrated that linear algorithms overestimate the dimensionality of nonlinear, noise-free data. In cases of high noise, most algorithms overestimated the dimensionality. We thus developed a denoising algorithm based on deep learning, the "Joint Autoencoder", which significantly improved subsequent dimensionality estimation. Critically, we found that all algorithms failed when the intrinsic dimensionality was high (above 20) or when the amount of data used for estimation was low. Based on the challenges we observed, we formulated a pipeline for estimating the dimensionality of experimental neural data.

Dynamic Electrode-to-Image (DETI) mapping reveals the human brain's spatiotemporal code of visual information

A number of neuroimaging techniques have been employed to understand how visual information is transformed along the visual pathway. Although each technique has spatial and temporal limitations, they can each provide important insights into the visual code. While the BOLD signal of fMRI can be quite informative, the visual code is not static and this can be obscured by fMRI’s poor temporal resolution. In this study, we leveraged the high temporal resolution of EEG to develop an encoding technique based on the distribution of responses generated by a population of real-world scenes. This approach maps neural signals to each pixel within a given image and reveals location-specific transformations of the visual code, providing a spatiotemporal signature for the image at each electrode. Our analyses of the mapping results revealed that scenes undergo a series of nonuniform transformations that prioritize different spatial frequencies at different regions of scenes over time. This mapping technique offers a potential avenue for future studies to explore how dynamic feedforward and recurrent processes inform and refine high-level representations of our visual world.

The visual information that we sample from our environment undergoes a series of neural modifications, with each modification state (or visual code) consisting of a unique distribution of responses across neurons along the visual pathway. However, current noninvasive neuroimaging techniques provide an account of that code that is coarse with respect to time or space. Here, we present dynamic electrode-to-image (DETI) mapping, an analysis technique that capitalizes on the high temporal resolution of EEG to map neural signals to each pixel within a given image to reveal location-specific modifications of the visual code. The DETI technique reveals maps of features that are associated with the neural signal at each pixel and at each time point. DETI mapping shows that real-world scenes undergo a series of nonuniform modifications over both space and time. Specifically, we find that the visual code varies in a location-specific manner, likely reflecting that neural processing prioritizes different features at different image locations over time. DETI mapping therefore offers a potential avenue for future studies to explore how each modification state informs and refines the conceptual meaning of our visual world.

A strategy for mapping biophysical to abstract neuronal network models applied to primary visual cortex

A fundamental challenge for the theoretical study of neuronal networks is to make the link between complex biophysical models based directly on experimental data, to progressively simpler mathematical models that allow the derivation of general operating principles. We present a strategy that successively maps a relatively detailed biophysical population model, comprising conductance-based Hodgkin-Huxley type neuron models with connectivity rules derived from anatomical data, to various representations with fewer parameters, finishing with a firing rate network model that permits analysis. We apply this methodology to primary visual cortex of higher mammals, focusing on the functional property of stimulus orientation selectivity of receptive fields of individual neurons. The mapping produces compact expressions for the parameters of the abstract model that clearly identify the impact of specific electrophysiological and anatomical parameters on the analytical results, in particular as manifested by specific functional signatures of visual cortex, including input-output sharpening, conductance invariance, virtual rotation and the tilt after effect. Importantly, qualitative differences between model behaviours point out consequences of various simplifications. The strategy may be applied to other neuronal systems with appropriate modifications.

A hierarchy of theoretical approaches to study a neuronal network depends on a tradeoff between biological fidelity and mathematical tractibility. Biophysically-detailed models consider cellular mechanisms and anatomically defined synaptic circuits, but are often too complex to reveal insights into fundamental principles. In contrast, increasingly abstract reduced models facilitate analytical insights. To better ground the latter to the underlying biology, we describe a systematic procedure to move across the model hierarchy that allows understanding how changes in biological parameters—physiological, pathophysiological, or because of new data—impact the behaviour of the network. We apply this approach to mammalian primary visual cortex, and examine how the different models in the hierarchy reproduce functional signatures of this area, in particular the tuning of neurons to the orientation of a visual stimulus. Our work provides a navigation of the complex parameter space of neural network models faithful to biology, as well as highlighting how simplifications made for mathematical convenience can fundamentally change their behaviour.

Unveiling functions of the visual cortex using task-specific deep neural networks

The human visual cortex enables visual perception through a cascade of hierarchical computations in cortical regions with distinct functionalities. Here, we introduce an AI-driven approach to discover the functional mapping of the visual cortex. We related human brain responses to scene images measured with functional MRI (fMRI) systematically to a diverse set of deep neural networks (DNNs) optimized to perform different scene perception tasks. We found a structured mapping between DNN tasks and brain regions along the ventral and dorsal visual streams. Low-level visual tasks mapped onto early brain regions, 3-dimensional scene perception tasks mapped onto the dorsal stream, and semantic tasks mapped onto the ventral stream. This mapping was of high fidelity, with more than 60% of the explainable variance in nine key regions being explained. Together, our results provide a novel functional mapping of the human visual cortex and demonstrate the power of the computational approach.

Topology-preserving smoothing of retinotopic maps

Retinotopic mapping, i.e., the mapping between visual inputs on the retina and neuronal activations in cortical visual areas, is one of the central topics in visual neuroscience. For human observers, the mapping is obtained by analyzing functional magnetic resonance imaging (fMRI) signals of cortical responses to slowly moving visual stimuli on the retina. Although it is well known from neurophysiology that the mapping is topological (i.e., the topology of neighborhood connectivity is preserved) within each visual area, retinotopic maps derived from the state-of-the-art methods are often not topological because of the low signal-to-noise ratio and spatial resolution of fMRI. The violation of topological condition is most severe in cortical regions corresponding to the neighborhood of the fovea (e.g., < 1 degree eccentricity in the Human Connectome Project (HCP) dataset), significantly impeding accurate analysis of retinotopic maps. This study aims to directly model the topological condition and generate topology-preserving and smooth retinotopic maps. Specifically, we adopted the Beltrami coefficient, a metric of quasiconformal mapping, to define the topological condition, developed a mathematical model to quantify topological smoothing as a constrained optimization problem, and elaborated an efficient numerical method to solve the problem. The method was then applied to V1, V2, and V3 simultaneously in the HCP dataset. Experiments with both simulated and real retinotopy data demonstrated that the proposed method could generate topological and smooth retinotopic maps.

Behavioral correlates of cortical semantic representations modeled by word vectors

The quantitative modeling of semantic representations in the brain plays a key role in understanding the neural basis of semantic processing. Previous studies have demonstrated that word vectors, which were originally developed for use in the field of natural language processing, provide a powerful tool for such quantitative modeling. However, whether semantic representations in the brain revealed by the word vector-based models actually capture our perception of semantic information remains unclear, as there has been no study explicitly examining the behavioral correlates of the modeled brain semantic representations. To address this issue, we compared the semantic structure of nouns and adjectives in the brain estimated from word vector-based brain models with that evaluated from human behavior. The brain models were constructed using voxelwise modeling to predict the functional magnetic resonance imaging (fMRI) response to natural movies from semantic contents in each movie scene through a word vector space. The semantic dissimilarity of brain word representations was then evaluated using the brain models. Meanwhile, data on human behavior reflecting the perception of semantic dissimilarity between words were collected in psychological experiments. We found a significant correlation between brain model- and behavior-derived semantic dissimilarities of words. This finding suggests that semantic representations in the brain modeled via word vectors appropriately capture our perception of word meanings.

Statistical power: Implications for planning MEG studies

Statistical power is key for robust, replicable science. Here, we systematically explored how numbers of trials and subjects affect statistical power in MEG sensor-level data. More specifically, we simulated "experiments" using the MEG resting-state dataset of the Human Connectome Project (HCP). We divided the data in two conditions, injected a dipolar source at a known anatomical location in the "signal condition", but not in the "noise condition", and detected significant differences at sensor level with classical paired t-tests across subjects, using amplitude, squared amplitude, and global field power (GFP) measures. Group-level detectability of these simulated effects varied drastically with anatomical origin. We thus examined in detail which spatial properties of the sources affected detectability, looking specifically at the distance from closest sensor and orientation of the source, and at the variability of these parameters across subjects. In line with previous single-subject studies, we found that the most detectable effects originate from source locations that are closest to the sensors and oriented tangentially with respect to the head surface. In addition, cross-subject variability in orientation also affected group-level detectability, boosting detection in regions where this variability was small and hindering detection in regions where it was large. Incidentally, we observed a considerable covariation of source position, orientation, and their cross-subject variability in individual brain anatomical space, making it difficult to assess the impact of each of these variables independently of one another. We thus also performed simulations where we controlled spatial properties independently of individual anatomy. These additional simulations confirmed the strong impact of distance and orientation and further showed that orientation variability across subjects affects detectability, whereas position variability does not. Importantly, our study indicates that strict unequivocal recommendations as to the ideal number of trials and subjects for any experiment cannot be realistically provided for neurophysiological studies and should be adapted according to the brain regions under study.

Aberrant Brain Functional Connectivity Strength and Effective Connectivity in Patients with Type 2 Diabetes Mellitus

Alterations of brain functional connectivity in patients with type 2 diabetes mellitus (T2DM) have been reported by resting-state functional magnetic resonance imaging studies, but the underlying precise neuropathological mechanism remains unclear. This study is aimed at investigating the implicit alterations of functional connections in T2DM by integrating functional connectivity strength (FCS) and Granger causality analysis (GCA) and further exploring their associations with clinical characteristics. Sixty T2DM patients and thirty-three sex-, age-, and education-matched healthy controls (HC) were recruited. Global FCS analysis of resting-state functional magnetic resonance imaging was performed to explore seed regions with significant differences between the two groups; then, GCA was applied to detect directional effective connectivity (EC) between the seeds and other brain regions. Correlations of EC with clinical variables were further explored in T2DM patients. Compared with HC, T2DM patients showed lower FCS in the bilateral fusiform gyrus, right superior frontal gyrus (SFG), and right postcentral gyrus, but higher FCS in the right supplementary motor area (SMA). Moreover, altered directional EC was found between the left fusiform gyrus and bilateral lingual gyrus and right medial frontal gyrus (MFG), as well as between the right SFG and bilateral frontal regions. In addition, triglyceride, insulin, and plasma glucose levels were correlated with the abnormal EC of the left fusiform, while disease duration and cognitive function were associated with the abnormal EC of the right SFG in T2DM patients. These results suggest that T2DM patients show aberrant brain function connectivity strength and effective connectivity which is associated with the diabetes-related metabolic characteristics, disease duration, and cognitive function, providing further insights into the complex neural basis of diabetes.

A meta-analysis of functional magnetic resonance imaging studies of divergent thinking using activation likelihood estimation

There are conflicting findings regarding brain regions and networks underpinning creativity, with divergent thinking tasks commonly used to study this. A handful of meta-analyses have attempted to synthesise findings on neural mechanisms of divergent thinking. With the rapid proliferation of research and recent developments in fMRI meta-analysis approaches, it is timely to reassess the regions activated during divergent thinking creativity tasks. Of particular interest is examining the evidence regarding large-scale brain networks proposed to be key in divergent thinking and extending this work to consider the role of the semantic control network. Studies utilising fMRI with healthy participants completing divergent thinking tasks were systematically identified, with 20 studies meeting the criteria. Activation Likelihood Estimation was then used to integrate the neuroimaging results across studies. This revealed four clusters: the left inferior parietal lobe; the left inferior frontal and precentral gyrus; the superior and medial frontal gyrus and the right cerebellum. These regions are key in the semantic network, important for flexible retrieval of stored knowledge, highlighting the role of this network in divergent thinking.

Prefiltering based on experimental paradigm for analysis of fMRI complex brain networks

Brain networks offers a new insight about connections between function and anatomical regions of human brain. We present results from brain networks built from functional magnetic resonance images during finger tapping paradigm. Pearson voxel-voxel correlation in time and frequency domains were performed for all subjects. Besides this standard framework we have implemented a new approach consisting in filtering the data with respect to the fMRI paradigm (finger tapping) in order to obtain a better understanding of the network involved in the execution of the task. The main topological graph measures have been compared in both cases: voxel-voxel correlation and voxel-paradigm filtering plus voxel-voxel correlation. With the standard voxel-voxel correlation a clearly free-scale network was obtained. On the other hand, when we prefiltered the paradigm we obtained two different kind of networks: 1) free-scale; 2) random-like. To our best knowledge, this behaviour is reported here for first time for brain networks. We suggest that paradigm signal prefiltering can provide more infomation about the brain networks.

Evaluating the Acute Effect of Stereoscopic Recovery by Dichoptic Stimulation Using Electroencephalogram

Amblyopia is a common developmental disorder in adolescents and children. Stereoscopic loss is a symptom of amblyopia that can seriously affect the quality of patient's life. Recent studies have shown that the push-pull perceptual learning protocol had a positive effect on stereoscopic recovery. In this study, we developed a stereoscopic training method using a polarized visualization system according to the push-pull protocol. Dichoptic stimulation for 36 anisometropic and amblyopic subjects and 33 children with normal visual acuity (VA) has been conducted. Electroencephalogram (EEG) was used to evaluate the neurophysiological changes before, during, and after stimulation. For the anisometropic and amblyopic subjects, the statistical analysis demonstrated significant differences (p < 0.01) in the beta rhythm at the middle temporal and occipital lobes, while the EEG from the normal VA subjects indicated no significant changes when comparing the results before and after training. We concluded that the dichoptic training in our study can activate the middle temporal visual area and visual cortex. The EEG changes can be used to evaluate the training effects. This study also found that the beta band EEG acquired during visual stimulation at the dorsal visual stream can be potentially used for predicting acute training effect. The results facilitated the optimization of the individual training plan.

Glucagon-Like Peptide-1 Receptor Agonist Differentially Affects Brain Activation in Response to Visual Food Cues in Lean and Obese Individuals with Type 2 Diabetes Mellitus

BACKGROUND

To investigate the effects of a glucagon-like peptide-1 receptor agonist on functional brain activation in lean and obese individuals with type 2 diabetes mellitus (T2DM) in response to visual food cues.

METHODS

In a randomized, single-blinded, crossover study, 15 lean and 14 obese individuals with T2DM were administered lixisenatide or normal saline subcutaneously with a 1-week washout period. We evaluated brain activation in response to pictures of high-calorie food, low-calorie food, and nonfood using functional magnetic resonance imaging and measured appetite and caloric intake in participants who were given access to an ad libitum buffet.

RESULTS

Obese individuals with T2DM showed significantly greater activation of the hypothalamus, pineal gland, parietal cortex (high-calorie food vs. low-calorie food, P<0.05), orbitofrontal cortex (high-calorie food vs. nonfood, P<0.05), and visual cortex (food vs. nonfood, P<0.05) than lean individuals with T2DM. Lixisenatide injection significantly reduced the functional activation of the fusiform gyrus and lateral ventricle in obese individuals with T2DM compared with that in lean individuals with T2DM (nonfood vs. high-calorie food, P<0.05). In addition, in individuals who decreased their caloric intake after lixisenatide injection, there were significant interaction effects between group and treatment in the posterior cingulate, medial frontal cortex (high-calorie food vs. low-calorie food, P<0.05), hypothalamus, orbitofrontal cortex, and temporal lobe (food vs. nonfood, P<0.05).

CONCLUSION

Brain responses to visual food cues were different in lean and obese individuals with T2DM. In addition, acute administration of lixisenatide differentially affected functional brain activation in these individuals, especially in those who decreased their caloric intake after lixisenatide injection.

Non-Invasive Functional-Brain-Imaging with an OPM-based Magnetoencephalography System

A non-invasive functional-brain-imaging system based on optically-pumped-magnetometers (OPM) is presented. The OPM-based magnetoencephalography (MEG) system features 20 OPM channels conforming to the subject's scalp. We have conducted two MEG experiments on three subjects: assessment of somatosensory evoked magnetic field (SEF) and auditory evoked magnetic field (AEF) using our OPM-based MEG system and a commercial MEG system based on superconducting quantum interference devices (SQUIDs). We cross validated the robustness of our system by calculating the distance between the location of the equivalent current dipole (ECD) yielded by our OPM-based MEG system and the ECD location calculated by the commercial SQUID-based MEG system. We achieved sub-centimeter accuracy for both SEF and AEF responses in all three subjects. Due to the proximity (12 mm) of the OPM channels to the scalp, it is anticipated that future OPM-based MEG systems will offer enhanced spatial resolution as they will capture finer spatial features compared to traditional MEG systems employing SQUIDs.

Quantitative analysis of intracranial electrocorticography signals using the concept of statistical parametric mapping

Statistical parametric mapping (SPM) is a technique with which one can delineate brain activity statistically deviated from the normative mean, and has been commonly employed in noninvasive neuroimaging and EEG studies. Using the concept of SPM, we developed a novel technique for quantification of the statistical deviation of an intracranial electrocorticography (ECoG) measure from the nonepileptic mean. We validated this technique using data previously collected from 123 patients with drug-resistant epilepsy who underwent resective epilepsy surgery. We determined how the measurement of statistical deviation of modulation index (MI) from the non-epileptic mean (rated by z-score) improved the performance of seizure outcome classification model solely based on conventional clinical, seizure onset zone (SOZ), and neuroimaging variables. Here, MI is a summary measure quantifying the strength of in-situ coupling between high-frequency activity at >150 Hz and slow wave at 3-4 Hz. We initially generated a normative MI atlas showing the mean and standard deviation of slow-wave sleep MI of neighboring non-epileptic channels of 47 patients, whose ECoG sampling involved all four lobes. We then calculated 'MI z-score' at each electrode site. SOZ had a greater 'MI z-score' compared to non-SOZ in the remaining 76 patients. Subsequent multivariate logistic regression analysis and receiver operating characteristic analysis to the combined data of all patients revealed that the full regression model incorporating all predictor variables, including SOZ and 'MI z-score', best classified the seizure outcome with sensitivity/specificity of 0.86/0.76. The model excluding 'MI z-score' worsened its sensitivity/specificity to 0.86/0.48. Furthermore, the leave-one-out analysis successfully cross-validated the full regression model. Measurement of statistical deviation of MI from the non-epileptic mean on invasive recording is technically feasible. Our analytical technique can be used to evaluate the utility of ECoG biomarkers in epilepsy presurgical evaluation.

Using centrality measures to extract core pattern of brain dynamics during the resting state

The patterns of brain dynamics were studied during resting state on a macroscopic scale for control subjects and multiple sclerosis patients. Macroscopic brain dynamics is defined after successive coarse-grainings and selection of significant patterns and transitions based on Markov representation of brain activity. The resulting networks show that control dynamics is merely organized according to a single principal pattern whereas patients dynamics depict more variable patterns. Centrality measures are used to extract core dynamical pattern in brain dynamics and classification technique allow to define MS dynamics with relevant error rate.

Spatial neuronal synchronization and the waveform of oscillations: Implications for EEG and MEG

Neuronal oscillations are ubiquitous in the human brain and are implicated in virtually all brain functions. Although they can be described by a prominent peak in the power spectrum, their waveform is not necessarily sinusoidal and shows rather complex morphology. Both frequency and temporal descriptions of such non-sinusoidal neuronal oscillations can be utilized. However, in non-invasive EEG/MEG recordings the waveform of oscillations often takes a sinusoidal shape which in turn leads to a rather oversimplified view on oscillatory processes. In this study, we show in simulations how spatial synchronization can mask non-sinusoidal features of the underlying rhythmic neuronal processes. Consequently, the degree of non-sinusoidality can serve as a measure of spatial synchronization. To confirm this empirically, we show that a mixture of EEG components is indeed associated with more sinusoidal oscillations compared to the waveform of oscillations in each constituent component. Using simulations, we also show that the spatial mixing of the non-sinusoidal neuronal signals strongly affects the amplitude ratio of the spectral harmonics constituting the waveform. Finally, our simulations show how spatial mixing can affect the strength and even the direction of the amplitude coupling between constituent neuronal harmonics at different frequencies. Validating these simulations, we also demonstrate these effects in real EEG recordings. Our findings have far reaching implications for the neurophysiological interpretation of spectral profiles, cross-frequency interactions, as well as for the unequivocal determination of oscillatory phase.

Two-step paretial least square regression classifiers in brain-state decoding using functional magnetic resonance imaging

Multivariate analysis methods have been widely applied to decode brain states from functional magnetic resonance imaging (fMRI) data. Among various multivariate analysis methods, partial least squares regression (PLSR) is often used to select relevant features for decoding brain states. However, PLSR is seldom directly used as a classifier to decode brain states from fMRI data. It is unclear how PLSR classifiers perform in brain-state decoding using fMRI. In this study, we propose two types of two-step PLSR classifiers that use PLSR/sparse PLSR (SPLSR) to select features and PLSR for classification to improve the performance of the PLSR classifier. The results of simulated and real fMRI data demonstrated that the PLSR classifier using PLSR/SPLSR to select features outperformed both the PLSR classifier using a general linear model (GLM) and the support vector machine (SVM) using PLSR/SPLSR/GLM in most cases. Moreover, PLSR using SPLSR to select features showed the best performance among all of the methods. Compared to GLM, PLSR is more sensitive in selecting the voxels that are specific to each task. The results suggest that the performance of the PLSR classifier can be largely improved when the PLSR classifier is combined with the feature selection methods of SPLSR and PLSR.

Detection and analysis of spatiotemporal patterns in brain activity

There is growing evidence that population-level brain activity is often organized into propagating waves that are structured in both space and time. Such spatiotemporal patterns have been linked to brain function and observed across multiple recording methodologies and scales. The ability to detect and analyze these patterns is thus essential for understanding the working mechanisms of neural circuits. Here we present a mathematical and computational framework for the identification and analysis of multiple classes of wave patterns in neural population-level recordings. By drawing a conceptual link between spatiotemporal patterns found in the brain and coherent structures such as vortices found in turbulent flows, we introduce velocity vector fields to characterize neural population activity. These vector fields are calculated for both phase and amplitude of oscillatory neural signals by adapting optical flow estimation methods from the field of computer vision. Based on these velocity vector fields, we then introduce order parameters and critical point analysis to detect and characterize a diverse range of propagating wave patterns, including planar waves, sources, sinks, spiral waves, and saddle patterns. We also introduce a novel vector field decomposition method that extracts the dominant spatiotemporal structures in a recording. This enables neural data to be represented by the activity of a small number of independent spatiotemporal modes, providing an alternative to existing dimensionality reduction techniques which separate space and time components. We demonstrate the capabilities of the framework and toolbox with simulated data, local field potentials from marmoset visual cortex and optical voltage recordings from whole mouse cortex, and we show that pattern dynamics are non-random and are modulated by the presence of visual stimuli. These methods are implemented in a MATLAB toolbox, which is freely available under an open-source licensing agreement.

Structured activity such as propagating wave patterns at the level of neural circuits can arise from highly variable firing activity of individual neurons. This property makes the brain, a quintessential example of a complex system, analogous to other complex physical systems such as turbulent fluids, in which structured patterns like vortices similarly emerge from molecules that behave irregularly. In this study, by uniquely adapting techniques for the identification of coherent structures in fluid turbulence, we develop new analytical and computational methods for the reliable detection of a diverse range of propagating wave patterns in large-scale neural recordings, for comprehensive analysis and visualization of these patterns, and for analysis of their dominant spatiotemporal modes. We demonstrate that these methods can be used to uncover the essential spatiotemporal properties of neural population activity recorded by different modalities, thus offering new insights into understanding the working mechanisms of neural systems.

Assessing robustness against potential publication bias in Activation Likelihood Estimation (ALE) meta-analyses for fMRI

The importance of integrating research findings is incontrovertible and procedures for coordinate-based meta-analysis (CBMA) such as Activation Likelihood Estimation (ALE) have become a popular approach to combine results of fMRI studies when only peaks of activation are reported. As meta-analytical findings help building cumulative knowledge and guide future research, not only the quality of such analyses but also the way conclusions are drawn is extremely important. Like classical meta-analyses, coordinate-based meta-analyses can be subject to different forms of publication bias which may impact results and invalidate findings. The file drawer problem refers to the problem where studies fail to get published because they do not obtain anticipated results (e.g. due to lack of statistical significance). To enable assessing the stability of meta-analytical results and determine their robustness against the potential presence of the file drawer problem, we present an algorithm to determine the number of noise studies that can be added to an existing ALE fMRI meta-analysis before spatial convergence of reported activation peaks over studies in specific regions is no longer statistically significant. While methods to gain insight into the validity and limitations of results exist for other coordinate-based meta-analysis toolboxes, such as Galbraith plots for Multilevel Kernel Density Analysis (MKDA) and funnel plots and egger tests for seed-based d mapping, this procedure is the first to assess robustness against potential publication bias for the ALE algorithm. The method assists in interpreting meta-analytical results with the appropriate caution by looking how stable results remain in the presence of unreported information that may differ systematically from the information that is included. At the same time, the procedure provides further insight into the number of studies that drive the meta-analytical results. We illustrate the procedure through an example and test the effect of several parameters through extensive simulations. Code to generate noise studies is made freely available which enables users to easily use the algorithm when interpreting their results.

Brain functional BOLD perturbation modelling for forward fMRI and inverse mapping

Purpose

To computationally separate dynamic brain functional BOLD responses from static background in a brain functional activity for forward fMRI signal analysis and inverse mapping.

Methods

A brain functional activity is represented in terms of magnetic source by a perturbation model: χ = χ₀ +δχ, with δχ for BOLD magnetic perturbations and χ₀ for background. A brain fMRI experiment produces a timeseries of complex-valued images (T2* images), whereby we extract the BOLD phase signals (denoted by δP) by a complex division. By solving an inverse problem, we reconstruct the BOLD δχ dataset from the δP dataset, and the brain χ distribution from a (unwrapped) T2* phase image. Given a 4D dataset of task BOLD fMRI, we implement brain functional mapping by temporal correlation analysis.

Results

Through a high-field (7T) and high-resolution (0.5mm in plane) task fMRI experiment, we demonstrated in detail the BOLD perturbation model for fMRI phase signal separation (P + δP) and reconstructing intrinsic brain magnetic source (χ and δχ). We also provided to a low-field (3T) and low-resolution (2mm) task fMRI experiment in support of single-subject fMRI study. Our experiments show that the δχ-depicted functional map reveals bidirectional BOLD χ perturbations during the task performance.

Conclusions

The BOLD perturbation model allows us to separate fMRI phase signal (by complex division) and to perform inverse mapping for pure BOLD δχ reconstruction for intrinsic functional χ mapping. The full brain χ reconstruction (from unwrapped fMRI phase) provides a new brain tissue image that allows to scrutinize the brain tissue idiosyncrasy for the pure BOLD δχ response through an automatic function/structure co-localization.

High-dimensional therapeutic inference in the focally damaged human brain

See Thiebaut de Schotten and Foulon (doi:10.1093/brain/awx332) for a scientific commentary on this article.Though consistency across the population renders the extraordinarily complex functional anatomy of the human brain surveyable, the inverse inference-from common functional maps to individual behaviour-is constrained by marked individual deviation from the population mean. Such inference is fundamental to the evaluation of therapeutic interventions in focal brain injury, where the impact of an induced structural change in the brain is quantified by its behavioural consequences, inevitably refracted through the lens of lesion-outcome relations. Current therapeutic evaluations do not incorporate inferences to the individual outcome derived from a detailed specification of the lesion anatomy, relying only on reductive parameters such as lesion volume and crudely discretised location. Examining 1172 patients with anatomically registered focal brain lesions, here we show that such low-dimensional models are highly insensitive to therapeutic effects. In contrast, high-dimensional models supported by machine learning dramatically improve sensitivity by leveraging complex individuating patterns in the functional architecture of the brain. The failure to replicate in humans positive interventional effects in experimental animals is thus revealed to have a remediable inferential cause, forcing a radical re-evaluation of therapeutic inference in the human brain.

Sparse feature learning for multi-class Parkinson's disease classification

This paper solves the multi-class classification problem for Parkinson's disease (PD) analysis by a sparse discriminative feature selection framework. Specifically, we propose a framework to construct a least square regression model based on the Fisher's linear discriminant analysis (LDA) and locality preserving projection (LPP). This framework utilizes the global and local information to select the most relevant and discriminative features to boost classification performance. Differing in previous methods for binary classification, we perform a multi-class classification for PD diagnosis. Our proposed method is evaluated on the public available Parkinson's progression markers initiative (PPMI) datasets. Extensive experimental results indicate that our proposed method identifies highly suitable regions for further PD analysis and diagnosis and outperforms state-of-the-art methods.

Sharing brain mapping statistical results with the neuroimaging data model

Only a tiny fraction of the data and metadata produced by an fMRI study is finally conveyed to the community. This lack of transparency not only hinders the reproducibility of neuroimaging results but also impairs future meta-analyses. In this work we introduce NIDM-Results, a format specification providing a machine-readable description of neuroimaging statistical results along with key image data summarising the experiment. NIDM-Results provides a unified representation of mass univariate analyses including a level of detail consistent with available best practices. This standardized representation allows authors to relay methods and results in a platform-independent regularized format that is not tied to a particular neuroimaging software package. Tools are available to export NIDM-Result graphs and associated files from the widely used SPM and FSL software packages, and the NeuroVault repository can import NIDM-Results archives. The specification is publically available at: http://nidm.nidash.org/specs/nidm-results.html.

Similar Spectral Power Densities Within the Schumann Resonance and a Large Population of Quantitative Electroencephalographic Profiles: Supportive Evidence for Koenig and Pobachenko

In 1954 and 1960 Koenig and his colleagues described the remarkable similarities of spectral power density profiles and patterns between the earth-ionosphere resonance and human brain activity which also share magnitudes for both electric field (mV/m) and magnetic field (pT) components. In 2006 Pobachenko and colleagues reported real time coherence between variations in the Schumann and brain activity spectra within the 6-16 Hz band for a small sample. We examined the ratios of the average potential differences (~3 μV) obtained by whole brain quantitative electroencephalography (QEEG) between rostral-caudal and left-right (hemispheric) comparisons of 238 measurements from 184 individuals over a 3.5 year period. Spectral densities for the rostral-caudal axis revealed a powerful peak at 10.25 Hz while the left-right peak was 1.95 Hz with beat-differences of ~7.5 to 8 Hz. When global cerebral measures were employed, the first (7-8 Hz), second (13-14 Hz) and third (19-20 Hz) harmonics of the Schumann resonances were discernable in averaged QEEG profiles in some but not all participants. The intensity of the endogenous Schumann resonance was related to the 'best-of-fitness' of the traditional 4-class microstate model. Additional measurements demonstrated real-time coherence for durations approximating microstates in spectral power density variations between Schumann frequencies measured in Sudbury, Canada and Cumiana, Italy with the QEEGs of local subjects. Our results confirm the measurements reported by earlier researchers that demonstrated unexpected similarities in the spectral patterns and strengths of electromagnetic fields generated by the human brain and the earth-ionospheric cavity.

Graph Matching: Relax at Your Own Risk

Graph matching-aligning a pair of graphs to minimize their edge disagreements-has received wide-spread attention from both theoretical and applied communities over the past several decades, including combinatorics, computer vision, and connectomics. Its attention can be partially attributed to its computational difficulty. Although many heuristics have previously been proposed in the literature to approximately solve graph matching, very few have any theoretical support for their performance. A common technique is to relax the discrete problem to a continuous problem, therefore enabling practitioners to bring gradient-descent-type algorithms to bear. We prove that an indefinite relaxation (when solved exactly) almost always discovers the optimal permutation, while a common convex relaxation almost always fails to discover the optimal permutation. These theoretical results suggest that initializing the indefinite algorithm with the convex optimum might yield improved practical performance. Indeed, experimental results illuminate and corroborate these theoretical findings, demonstrating that excellent results are achieved in both benchmark and real data problems by amalgamating the two approaches.

Structural and Functional Alterations in Right Dorsomedial Prefrontal and Left Insular Cortex Co-Localize in Adolescents with Aggressive Behaviour: An ALE Meta-Analysis

Recent neuroimaging work has suggested that aggressive behaviour (AB) is associated with structural and functional brain abnormalities in processes subserving emotion processing and regulation. However, most neuroimaging studies on AB to date only contain relatively small sample sizes. To objectively investigate the consistency of previous structural and functional research in adolescent AB, we performed a systematic literature review and two coordinate-based activation likelihood estimation meta-analyses on eight VBM and nine functional neuroimaging studies in a total of 783 participants (408 [224AB/184 controls] and 375 [215 AB/160 controls] for structural and functional analysis respectively). We found 19 structural and eight functional foci of significant alterations in adolescents with AB, mainly located within the emotion processing and regulation network (including orbitofrontal, dorsomedial prefrontal and limbic cortex). A subsequent conjunction analysis revealed that functional and structural alterations co-localize in right dorsomedial prefrontal cortex and left insula. Our results are in line with meta-analytic work as well as structural, functional and connectivity findings to date, all of which make a strong point for the involvement of a network of brain areas responsible for emotion processing and regulation, which is disrupted in AB. Increased knowledge about the behavioural and neuronal underpinnings of AB is crucial for the development of novel and implementation of existing treatment strategies. Longitudinal research studies will have to show whether the observed alterations are a result or primary cause of the phenotypic characteristics in AB.

Persistent vegetative state and minimally conscious state: a systematic review and meta-analysis of diagnostic procedures

BACKGROUND

Acute brain damage can cause major disturbances of consciousness, ranging all the way to the persistent vegetative state (PVS), which is also known as "unresponsive wakefulness syndrome". PVS can be hard to distinguish from a state of minimal preserved consciousness ("minimally conscious state," MCS); the rate of misdiagnosis is high and has been estimated at 37-43%. In contrast, PVS is easily distinguished from brain death. We discuss the various diagnostic techniques that can be used to determine whether a patient is minimally conscious or in a persistent vegetative state.

METHODS

This article is based on a systematic review of pertinent literature and on a quantitative meta-analysis of the sensitivity and specificity of new diagnostic methods for the minimally conscious state.

RESULTS

We identified and evaluated 20 clinical studies involving a total of 906 patients with either PVS or MCS. The reported sensitivities and specificities of the various techniques used to diagnose MCS vary widely. The sensitivity and specificity of functional MRI-based techniques are 44% and 67%, respectively (with corresponding 95% confidence intervals of 19%-72% and 55%-77%); those of quantitative EEG are 90% and 80%, respectively (95% CI, 69%-97% and 66%-90%). EEG, event-related potentials, and imaging studies can also aid in prognostication. Contrary to prior assumptions, 10% to 24% of patients in PVS can regain consciousness, sometimes years after the event, but only with marked functional impairment.

CONCLUSION

The basic diagnostic evaluation for differentiating PVS from MCS consists of a standardized clinical examination. In the future, modern diagnostic techniques may help identify patients who are in a subclinical minimally conscious state.

An antidote to the imager's fallacy, or how to identify brain areas that are in limbo

Traditionally, fMRI data are analyzed using statistical parametric mapping approaches. Regardless of the precise thresholding procedure, these approaches ultimately divide the brain in regions that do or do not differ significantly across experimental conditions. This binary classification scheme fosters the so-called imager's fallacy, where researchers prematurely conclude that region A is selectively involved in a certain cognitive task because activity in that region reaches statistical significance and activity in region B does not. For such a conclusion to be statistically valid, however, a test on the differences in activation across these two regions is required. Here we propose a simple GLM-based method that defines an "in-between" category of brain regions that are neither significantly active nor inactive, but rather "in limbo". For regions that are in limbo, the activation pattern is inconclusive: it does not differ significantly from baseline, but neither does it differ significantly from regions that do show significant changes from baseline. This pattern indicates that measurement was insufficiently precise. By directly testing differences in activation, our procedure helps reduce the impact of the imager's fallacy. The method is illustrated using concrete examples.

EEG source connectivity analysis: from dense array recordings to brain networks

The recent past years have seen a noticeable increase of interest for electroencephalography (EEG) to analyze functional connectivity through brain sources reconstructed from scalp signals. Although considerable advances have been done both on the recording and analysis of EEG signals, a number of methodological questions are still open regarding the optimal way to process the data in order to identify brain networks. In this paper, we analyze the impact of three factors that intervene in this processing: i) the number of scalp electrodes, ii) the combination between the algorithm used to solve the EEG inverse problem and the algorithm used to measure the functional connectivity and iii) the frequency bands retained to estimate the functional connectivity among neocortical sources. Using High-Resolution (hr) EEG recordings in healthy volunteers, we evaluated these factors on evoked responses during picture recognition and naming task. The main reason for selection this task is that a solid literature background is available about involved brain networks (ground truth). From this a priori information, we propose a performance criterion based on the number of connections identified in the regions of interest (ROI) that belong to potentially activated networks. Our results show that the three studied factors have a dramatic impact on the final result (the identified network in the source space) as strong discrepancies were evidenced depending on the methods used. They also suggest that the combination of weighted Minimum Norm Estimator (wMNE) and the Phase Synchronization (PS) methods applied on High-Resolution EEG in beta/gamma bands provides the best performance in term of topological distance between the identified network and the expected network in the above-mentioned cognitive task.

Coupled Intrinsic Connectivity Distribution analysis: a method for exploratory connectivity analysis of paired FMRI data

We present a novel voxel-based connectivity approach for paired functional magnetic resonance imaging (fMRI) data collected under two different conditions labeled the Coupled Intrinsic Connectivity Distribution (coupled-ICD). Our proposed method jointly models both conditions to incorporate additional paired information into the connectivity metric. Voxel-based connectivity holds promise as a clinical tool to characterize a wide range of neurological and psychiatric diseases, and monitor their treatment. As such, examining paired connectivity data such as scans acquired pre- and post-intervention is an important application for connectivity methodologically. When presented with data from paired conditions, conventional voxel-based methods analyze each condition separately. However, summarizing each connection separately can misrepresent patterns of changes in connectivity. We show that commonly used methods can underestimate functional changes and subsequently introduce and evaluate our solution to this problem, the coupled-ICD metric, using two studies: 1) healthy controls scanned awake and under anesthesia, and 2) cocaine-dependent subjects and healthy controls scanned while being presented with relaxing or drug-related imagery cues. The coupled-ICD approach detected differences between paired conditions in similar brain regions as the conventional approaches while also revealing additional changes in regions not identified using conventional voxel-based connectivity analyses. Follow-up seed-based analyses on data independent from the voxel-based results also showed connectivity differences between conditions in regions detected by coupled-ICD. This approach of jointly analyzing paired resting-state scans provides a new and important tool with many applications for clinical and basic neuroscience research.

Knowledge-guided robust MRI brain extraction for diverse large-scale neuroimaging studies on humans and non-human primates

Accurate and robust brain extraction is a critical step in most neuroimaging analysis pipelines. In particular, for the large-scale multi-site neuroimaging studies involving a significant number of subjects with diverse age and diagnostic groups, accurate and robust extraction of the brain automatically and consistently is highly desirable. In this paper, we introduce population-specific probability maps to guide the brain extraction of diverse subject groups, including both healthy and diseased adult human populations, both developing and aging human populations, as well as non-human primates. Specifically, the proposed method combines an atlas-based approach, for coarse skull-stripping, with a deformable-surface-based approach that is guided by local intensity information and population-specific prior information learned from a set of real brain images for more localized refinement. Comprehensive quantitative evaluations were performed on the diverse large-scale populations of ADNI dataset with over 800 subjects (55∼90 years of age, multi-site, various diagnosis groups), OASIS dataset with over 400 subjects (18∼96 years of age, wide age range, various diagnosis groups), and NIH pediatrics dataset with 150 subjects (5∼18 years of age, multi-site, wide age range as a complementary age group to the adult dataset). The results demonstrate that our method consistently yields the best overall results across almost the entire human life span, with only a single set of parameters. To demonstrate its capability to work on non-human primates, the proposed method is further evaluated using a rhesus macaque dataset with 20 subjects. Quantitative comparisons with popularly used state-of-the-art methods, including BET, Two-pass BET, BET-B, BSE, HWA, ROBEX and AFNI, demonstrate that the proposed method performs favorably with superior performance on all testing datasets, indicating its robustness and effectiveness.

Principal feature analysis: a multivariate feature selection method for fMRI data

Brain decoding with functional magnetic resonance imaging (fMRI) requires analysis of complex, multivariate data. Multivoxel pattern analysis (MVPA) has been widely used in recent years. MVPA treats the activation of multiple voxels from fMRI data as a pattern and decodes brain states using pattern classification methods. Feature selection is a critical procedure of MVPA because it decides which features will be included in the classification analysis of fMRI data, thereby improving the performance of the classifier. Features can be selected by limiting the analysis to specific anatomical regions or by computing univariate (voxel-wise) or multivariate statistics. However, these methods either discard some informative features or select features with redundant information. This paper introduces the principal feature analysis as a novel multivariate feature selection method for fMRI data processing. This multivariate approach aims to remove features with redundant information, thereby selecting fewer features, while retaining the most information.

Population level inference for multivariate MEG analysis

Multivariate analysis is a very general and powerful technique for analysing Magnetoencephalography (MEG) data. An outstanding problem however is how to make inferences that are consistent over a group of subjects as to whether there are condition-specific differences in data features, and what are those features that maximise these differences. Here we propose a solution based on Canonical Variates Analysis (CVA) model scoring at the subject level and random effects Bayesian model selection at the group level. We apply this approach to beamformer reconstructed MEG data in source space. CVA estimates those multivariate patterns of activation that correlate most highly with the experimental design; the order of a CVA model is then determined by the number of significant canonical vectors. Random effects Bayesian model comparison then provides machinery for inferring the optimal order over the group of subjects. Absence of a multivariate dependence is indicated by the null model being the most likely. This approach can also be applied to CVA models with a fixed number of canonical vectors but supplied with different feature sets. We illustrate the method by identifying feature sets based on variable-dimension MEG power spectra in the primary visual cortex and fusiform gyrus that are maximally discriminative of data epochs before versus after visual stimulation.

Identifying preseizure state in intracranial EEG data using diffusion kernels

The goal of this study is to identify preseizure changes in intracranial EEG (icEEG). A novel approach based on the recently developed diffusion map framework, which is considered to be one of the leading manifold learning methods, is proposed. Diffusion mapping provides dimensionality reduction of the data as well as pattern recognition that can be used to distinguish different states of the patient, for example, interictal and preseizure. A new algorithm, which is an extension of diffusion maps, is developed to construct coordinates that generate efficient geometric representations of the complex structures in the icEEG data. In addition, this method is adapted to the icEEG data and enables the extraction of the underlying brain activity. The algorithm is tested on icEEG data recorded from several electrode contacts from a patient being evaluated for possible epilepsy surgery at the Yale-New Haven Hospital. Numerical results show that the proposed approach provides a distinction between interictal and preseizure states.

A mixed L2 norm regularized HRF estimation method for rapid event-related fMRI experiments

Brain state decoding or "mind reading" via multivoxel pattern analysis (MVPA) has become a popular focus of functional magnetic resonance imaging (fMRI) studies. In brain decoding, stimulus presentation rate is increased as fast as possible to collect many training samples and obtain an effective and reliable classifier or computational model. However, for extremely rapid event-related experiments, the blood-oxygen-level-dependent (BOLD) signals evoked by adjacent trials are heavily overlapped in the time domain. Thus, identifying trial-specific BOLD responses is difficult. In addition, voxel-specific hemodynamic response function (HRF), which is useful in MVPA, should be used in estimation to decrease the loss of weak information across voxels and obtain fine-grained spatial information. Regularization methods have been widely used to increase the efficiency of HRF estimates. In this study, we propose a regularization framework called mixed L2 norm regularization. This framework involves Tikhonov regularization and an additional L2 norm regularization term to calculate reliable HRF estimates. This technique improves the accuracy of HRF estimates and significantly increases the classification accuracy of the brain decoding task when applied to a rapid event-related four-category object classification experiment. At last, some essential issues such as the impact of low-frequency fluctuation (LFF) and the influence of smoothing are discussed for rapid event-related experiments.

Iterative cross-correlation analysis of resting state functional magnetic resonance imaging data

Seed-based cross-correlation analysis (sCCA) and independent component analysis have been widely employed to extract functional networks from the resting state functional magnetic resonance imaging data. However, the results of sCCA, in terms of both connectivity strength and network topology, can be sensitive to seed selection variations. ICA avoids the potential problems due to seed selection, but choosing which component(s) to represent the network of interest could be subjective and problematic. In this study, we proposed a seed-based iterative cross-correlation analysis (siCCA) method for resting state brain network analysis. The method was applied to extract default mode network (DMN) and stable task control network (STCN) in two independent datasets acquired from normal adults. Compared with the networks obtained by traditional sCCA and ICA, the resting state networks produced by siCCA were found to be highly stable and independent on seed selection. siCCA was used to analyze DMN in first-episode major depressive disorder (MDD) patients. It was found that, in the MDD patients, the volume of DMN negatively correlated with the patients' social disability screening schedule scores.

Evidence of diagnostic specificity in the neural correlates of facial affect processing in bipolar disorder and schizophrenia: a meta-analysis of functional imaging studies

BACKGROUND

Schizophrenia (SZ) and bipolar disorder (BD) may overlap in etiology and phenomenology but differ with regard to emotional processing. We used facial affect as a probe for emotional processing to determine whether there are diagnosis-related differences between SZ and BD in the function of the underlying neural circuitry.

METHOD

Functional magnetic resonance imaging (fMRI) studies published up to 30 April 2012 investigating facial affect processing in patients with SZ or BD were identified through computerized and manual literature searches. Activation foci from 29 studies encompassing 483 healthy individuals, 268 patients with SZ and 267 patients with BD were subjected to voxel-based quantitative meta-analysis using activation likelihood estimation (ALE).

RESULTS

Compared to healthy individuals, when emotional facial stimuli were contrasted to neutral stimuli, patients with BD showed overactivation within the parahippocampus/amygdala and thalamus and reduced engagement within the ventrolateral prefrontal cortex (PFC) whereas patients with SZ showed underactivation throughout the entire facial affect processing network and increased activation in visual processing regions within the cuneus. Patients with BD showed greater thalamic engagement compared to patients with SZ; in the reverse comparison, patients with SZ showed greater engagement in posterior associative visual cortices.

CONCLUSIONS

During facial affect processing, patients with BD show overactivation in subcortical regions and underactivation in prefrontal regions of the facial affect processing network, consistent with the notion of reduced emotional regulation. By contrast, overactivation within visual processing regions coupled with reduced engagement of facial affect processing regions points to abnormal visual integration as the core underlying deficit in SZ.

Decoding the anatomical network of spatial attention

The study of stroke patients with modern lesion-symptom analysis techniques has yielded valuable insights into the representation of spatial attention in the human brain. Here we introduce an approach--multivariate pattern analysis--that no longer assumes independent contributions of brain regions but rather quantifies the joint contribution of multiple brain regions in determining behavior. In a large sample of stroke patients, we found patterns of damage more predictive of spatial neglect than the best-performing single voxel. In addition, modeling multiple brain regions--those that are frequently damaged and, importantly, spared--provided more predictive information than modeling single regions. Interestingly, we also found that the superior temporal gyrus demonstrated a consistent ability to improve classifier performance when added to other regions, implying uniquely predictive information. In sharp contrast, classifier performance for both the angular gyrus and insular cortex was reliably enhanced by the addition of other brain regions, suggesting these regions lack independent predictive information for spatial neglect. Our findings highlight the utility of multivariate pattern analysis in lesion mapping, furnishing neuroscience with a modern approach for using lesion data to study human brain function.

Reliable identification of deep sulcal pits: the effects of scan session, scanner, and surface extraction tool

Sulcal pit analysis has been providing novel insights into brain function and development. The purpose of this study was to evaluate the reliability of sulcal pit extraction with respect to the effects of scan session, scanner, and surface extraction tool. Five subjects were scanned 4 times at 3 MRI centers and other 5 subjects were scanned 3 times at 2 MRI centers, including 1 test-retest session. Sulcal pits were extracted on the white matter surfaces reconstructed with both Montreal Neurological Institute and Freesurfer pipelines. We estimated similarity of the presence of sulcal pits having a maximum value of 1 and their spatial difference within the same subject. The tests showed high similarity of the sulcal pit presence and low spatial difference. The similarity was more than 0.90 and the spatial difference was less than 1.7 mm in most cases according to different scan sessions or scanners, and more than 0.85 and about 2.0 mm across surface extraction tools. The reliability of sulcal pit extraction was more affected by the image processing-related factors than the scan session or scanner factors. Moreover, the similarity of sulcal pit distribution appeared to be largely influenced by the presence or absence of the sulcal pits on the shallow and small folds. We suggest that our sulcal pit extraction from MRI is highly reliable and could be useful for clinical applications as an imaging biomarker.

Publication bias in neuroimaging research: implications for meta-analyses

Neuroimaging and the neurosciences have made notable advances in sharing activation results through detailed databases, making meta-analysis of the published research faster and easier. However, the effect of publication bias in these fields has not been previously addressed or accounted for in the developed meta-analytic methods. In this article, we examine publication bias in functional magnetic resonance imaging (fMRI) for tasks involving working memory in the frontal lobes (Brodmann Areas 4, 6, 8, 9, 10, 37, 45, 46, and 47). Seventy-four studies were selected from the literature and the effect of publication bias was examined using a number of regression-based techniques. Pearson's r correlation coefficient and Cohen's d effect size estimates were computed for the activation in each study and compared to the study sample size using Egger's regression, Macaskill's regression, and the 'Trim and Fill' method. Evidence for publication bias was identified in this body of literature (p < 0.01 for each test), generally, though was neither task- nor sub-region-dependent. While we focused our analysis on this subgroup of brain mapping studies, we believe our findings generalize to the brain imaging literature as a whole and databases seeking to curate their collective results. While neuroimaging databases of summary effects are of enormous value to the community, the potential publication bias should be considered when performing meta-analyses based on database contents.

Age differences in symbolic representations of motor sequence learning

We recently reported that young adults (YA) preferentially recruit cerebellar lobule HVI for symbolic motor sequence learning [3]. Learning magnitude in the symbolic condition was correlated with activation level in lobule HVI. Here, we evaluated age differences in the symbolic representation of motor sequence learning. Fourteen YA and 14 older adults (OA) performed the alternating serial reaction time task (ASRT) under conditions in which the spatial processing component was selectively eliminated from stimulus presentation (spatial versus symbolic), response execution (manual versus vocal), or both. Results showed that OA had reduced learning magnitudes relative to YA. Using the cerebellum lobule HVI as a region-of-interest, we found that OA had significantly lower activation in this region than YA during the symbolic learning conditions (FWE, P<0.05). Similar to YA, OA also showed a significant correlation between learning magnitude and cerebellar activation in the symbolic conditions. These results suggest that although YA and OA recruit similar neural networks during implicit learning, OA under-recruit relevant brain areas which may partially explain their implicit sequence learning deficits.

Special surgical considerations for functional brain mapping

The development of functional mapping techniques gives neurosurgeons many options for preoperative planning. Integrating functional and anatomic data can inform patient selection and surgical planning and makes functional mapping more accessible than when only invasive studies were available. However, the applications of functional mapping to neurosurgical patients are still evolving. Functional imaging remains complex and requires an understanding of the underlying physiologic and imaging characteristics. Neurosurgeons must be accustomed to interpreting highly processed data. Successful implementation of functional image-guided procedures requires efficient interactions between neurosurgeon, neurologist, radiologist, neuropsychologist, and others, but promises to enhance the care of patients.

Comparative study of SVM methods combined with voxel selection for object category classification on fMRI data

BACKGROUND

Support vector machine (SVM) has been widely used as accurate and reliable method to decipher brain patterns from functional MRI (fMRI) data. Previous studies have not found a clear benefit for non-linear (polynomial kernel) SVM versus linear one. Here, a more effective non-linear SVM using radial basis function (RBF) kernel is compared with linear SVM. Different from traditional studies which focused either merely on the evaluation of different types of SVM or the voxel selection methods, we aimed to investigate the overall performance of linear and RBF SVM for fMRI classification together with voxel selection schemes on classification accuracy and time-consuming.

METHODOLOGY/PRINCIPAL FINDINGS

Six different voxel selection methods were employed to decide which voxels of fMRI data would be included in SVM classifiers with linear and RBF kernels in classifying 4-category objects. Then the overall performances of voxel selection and classification methods were compared. Results showed that: (1) Voxel selection had an important impact on the classification accuracy of the classifiers: in a relative low dimensional feature space, RBF SVM outperformed linear SVM significantly; in a relative high dimensional space, linear SVM performed better than its counterpart; (2) Considering the classification accuracy and time-consuming holistically, linear SVM with relative more voxels as features and RBF SVM with small set of voxels (after PCA) could achieve the better accuracy and cost shorter time.

CONCLUSIONS/SIGNIFICANCE

The present work provides the first empirical result of linear and RBF SVM in classification of fMRI data, combined with voxel selection methods. Based on the findings, if only classification accuracy was concerned, RBF SVM with appropriate small voxels and linear SVM with relative more voxels were two suggested solutions; if users concerned more about the computational time, RBF SVM with relative small set of voxels when part of the principal components were kept as features was a better choice.

A simple model can unify a broad range of phenomena in retinotectal map development

A paradigm model system for studying the development of patterned connections in the nervous system is the topographic map formed by retinal axons in the optic tectum/superior colliculus. Starting in the 1970s, a series of computational models have been proposed to explain map development in both normal conditions, and perturbed conditions where the retina and/or tectum/superior colliculus are altered. This stands in contrast to more recent models that have often been simpler than older ones, and tend to address more limited data sets, but include more recent genetic manipulations. The original exploration of many of the early models was one-dimensional and limited by the computational resources of the time. This leaves open the ability of these early models to explain both map development in two dimensions, and the genetic manipulation data that have only appeared more recently. In this article, we show that a two-dimensional and updated version of the XBAM model (eXtended Branch Arrow Model), first proposed in 1982, reproduces a range of surgical map manipulations not yet demonstrated by more modern models. A systematic exploration of the parameter space of this model in two dimensions also reveals richer behavior than that apparent from the original one-dimensional versions. Furthermore, we show that including a specific type of axon-axon interaction can account for the map collapse recently observed when particular receptor levels are genetically manipulated in a subset of retinal ganglion cells. Together these results demonstrate that balancing multiple influences on map development seems to be necessary to explain many biological phenomena in retinotectal map formation, and suggest important constraints on the underlying biological variables.

An integrated multimodality MR brain imaging study: gray matter tissue loss mediates the association between cerebral hypoperfusion and Alzheimer's disease

Multimodality MR image processing and analysis aims to determine to what extent information from various imaging modalities is redundant or complementary and how changes in various regions of the brain, detected by various modalities, interact with each other to produce cognitive and functional changes. Here, we presented a multimodality image processing framework to integrate unique and complementary observations on anatomical and physiological brain changes measured by structural and perfusion MR imaging modalities. A unique aspect to this study is that we performed a test on mediation hypothesis that requires a model based on measurements from both MR imaging modalities. Our findings from these integrated multi-modality analysis were congruent with previously published results on neuropathology of AD patients and supplied a comprehensive explanation on the integration between structural atrophy and physiological deterioration due to AD.