Threat and Reward Imminence Processing in the Human Brain.

In the human brain, aversive and appetitive processing have been studied with controlled stimuli in rather static settings. In addition, the extent to which aversive-related and appetitive-related processing engage distinct or overlapping circuits remains poorly understood. Here, we sought to investigate the dynamics of aversive and appetitive processing while male and female participants engaged in comparable trials involving threat avoidance or reward seeking. A central goal was to characterize the temporal evolution of responses during periods of threat or reward imminence. For example, in the aversive domain, we predicted that the bed nucleus of the stria terminalis (BST), but not the amygdala, would exhibit anticipatory responses given the role of the former in anxious apprehension. We also predicted that the periaqueductal gray (PAG) would exhibit threat-proximity responses based on its involvement in proximal-threat processes, and that the ventral striatum would exhibit threat-imminence responses given its role in threat escape in rodents. Overall, we uncovered imminence-related temporally increasing ("ramping") responses in multiple brain regions, including the BST, PAG, and ventral striatum, subcortically, and dorsal anterior insula and anterior midcingulate, cortically. Whereas the ventral striatum generated anticipatory responses in the proximity of reward as expected, it also exhibited threat-related imminence responses. In fact, across multiple brain regions, we observed a main effect of arousal. In other words, we uncovered extensive temporally evolving, imminence-related processing in both the aversive and appetitive domain, suggesting that distributed brain circuits are dynamically engaged during the processing of biologically relevant information regardless of valence, findings further supported by network analysis.SIGNIFICANCE STATEMENT In the human brain, aversive and appetitive processing have been studied with controlled stimuli in rather static settings. Here, we sought to investigate the dynamics of aversive/appetitive processing while participants engaged in trials involving threat avoidance or reward seeking. A central goal was to characterize the temporal evolution of responses during periods of threat or reward imminence. We uncovered imminence-related temporally increasing ("ramping") responses in multiple brain regions, including the bed nucleus of the stria terminalis, periaqueductal gray, and ventral striatum, subcortically, and dorsal anterior insula and anterior midcingulate, cortically. Overall, we uncovered extensive temporally evolving, imminence-related processing in both the aversive and appetitive domain, suggesting that distributed brain circuits are dynamically engaged during the processing of biologically relevant information regardless of valence.

The Spiraling Cognitive-Emotional Brain: Combinatorial, Reciprocal, and Reentrant Macro-organization.

This article proposes a framework for understanding the macro-scale organization of anatomical pathways in the mammalian brain. The architecture supports flexible behavioral decisions across a spectrum of spatio-temporal scales. The proposal emphasizes the combinatorial, reciprocal, and reentrant connectivity-called CRR neuroarchitecture-between cortical, BG, thalamic, amygdala, hypothalamic, and brainstem circuits. Thalamic nuclei, especially midline/intralaminar nuclei, are proposed to act as hubs routing the flow of signals between noncortical areas and pFC. The hypothalamus also participates in multiregion circuits via its connections with cortex and thalamus. At slower timescales, long-range behaviors integrate signals across levels of the neuroaxis. At fast timescales, parallel engagement of pathways allows urgent behaviors while retaining flexibility. Overall, the proposed architecture enables context-dependent, adaptive behaviors spanning proximate to distant spatio-temporal scales. The framework promotes an integrative perspective and a distributed, heterarchical view of brain function.

The Entangled Brain.

The Entangled Brain (Pessoa, L., 2002. MIT Press) promotes the idea that we need to understand the brain as a complex, entangled system. Why does the complex systems perspective, one that entails emergent properties, matter for brain science? In fact, many neuroscientists consider these ideas a distraction. We discuss three principles of brain organization that inform the question of the interactional complexity of the brain: (1) massive combinatorial anatomical connectivity; (2) highly distributed functional coordination; and (3) networks/circuits as functional units. To motivate the challenges of mapping structure and function, we discuss neural circuits illustrating the high anatomical and functional interactional complexity typical in the brain. We discuss potential avenues for testing for network-level properties, including those relying on distributed computations across multiple regions. We discuss implications for brain science, including the need to characterize decentralized and heterarchical anatomical-functional organization. The view advocated has important implications for causation, too, because traditional accounts of causality provide poor candidates for explanation in interactionally complex systems like the brain given the distributed, mutual, and reciprocal nature of the interactions. Ultimately, to make progress understanding how the brain supports complex mental functions, we need to dissolve boundaries within the brain-those suggested to be associated with perception, cognition, action, emotion, motivation-as well as outside the brain, as we bring down the walls between biology, psychology, mathematics, computer science, philosophy, and so on.

Distributed and Multifaceted Effects of Threat and Safety.

In the present fMRI study, we examined how anxious apprehension is processed in the human brain. A central goal of the study was to test the prediction that a subset of brain regions would exhibit sustained response profiles during threat periods, including the anterior insula, a region implicated in anxiety disorders. A second important goal was to evaluate the responses in the amygdala and the bed nucleus of the stria terminals, regions that have been suggested to be involved in more transient and sustained threat, respectively. A total of 109 participants performed an experiment in which they encountered "threat" or "safe" trials lasting approximately 16 sec. During the former, they experienced zero to three highly unpleasant electrical stimulations, whereas in the latter, they experienced zero to three benign electrical stimulations (not perceived as unpleasant). The timing of the stimulation during trials was randomized, and as some trials contained no stimulation, stimulation delivery was uncertain. We contrasted responses during threat and safe trials that did not contain electrical stimulation, but only the potential that unpleasant (threat) or benign (safe) stimulation could occur. We employed Bayesian multilevel analysis to contrast responses to threat and safe trials in 85 brain regions implicated in threat processing. Our results revealed that the effect of anxious apprehension is distributed across the brain and that the temporal evolution of the responses is quite varied, including more transient and more sustained profiles, as well as signal increases and decreases with threat.

Learning brain dynamics for decoding and predicting individual differences.

Insights from functional Magnetic Resonance Imaging (fMRI), as well as recordings of large numbers of neurons, reveal that many cognitive, emotional, and motor functions depend on the multivariate interactions of brain signals. To decode brain dynamics, we propose an architecture based on recurrent neural networks to uncover distributed spatiotemporal signatures. We demonstrate the potential of the approach using human fMRI data during movie-watching data and a continuous experimental paradigm. The model was able to learn spatiotemporal patterns that supported 15-way movie-clip classification (∼90%) at the level of brain regions, and binary classification of experimental conditions (∼60%) at the level of voxels. The model was also able to learn individual differences in measures of fluid intelligence and verbal IQ at levels comparable to that of existing techniques. We propose a dimensionality reduction approach that uncovers low-dimensional trajectories and captures essential informational (i.e., classification related) properties of brain dynamics. Finally, saliency maps and lesion analysis were employed to characterize brain-region/voxel importance, and uncovered how dynamic but consistent changes in fMRI activation influenced decoding performance. When applied at the level of voxels, our framework implements a dynamic version of multivariate pattern analysis. Our approach provides a framework for visualizing, analyzing, and discovering dynamic spatially distributed brain representations during naturalistic conditions.

Effectiveness of soil conditioners to enhance salt extraction ability of Salicornia ramosissima in saline-sodic soil for different soil moisture contents.

Soil salinity is considered one of the main types of soil degradation in semiarid environments around the globe. This work aims to evaluate the effectiveness of soil conditioners to enhance the growth and salt extraction ability of Salicornia ramosíssima for different soil moisture contents. Salicornia plants were cultivated in pots in which the soils were treated with the following conditioners: control; gypsum + organic matter; elemental sulfur + organic matter; and gypsum + elemental sulfur + organic matter. Salicornia plants were subjected to two soil moisture rates - at 35 and 85% field capacity. Soil conditioners associated with higher contents of soil moisture promoted significant increases, compared to control, in fresh (6.20 - 11.13 g) and dry matter (1.20 - 2.07 g), relative biomass (100 - 179%) as well as significantly increased the concentrations of Na+ (56.09 - 65.64 mg kg-1) and Cl- (110.83 - 150.0 mg kg-1) in plant tissues. Soil conditioners significantly increased salt extraction ability under the two moisture levels, mainly by promoting higher values for both transfer factor and phytoremediation potential. The best performance of Salicornia in terms of plant yield and salt extraction, regardless of the moisture level, was the gypsum + organic matter.Novelty statementThere are no studies in the literature relating the use of conditioners as a strategy to enhance Salicornia's ability to extract salts.This work contributes to the management of salinized areas around the globe in two main aspects. The first is that many of these salt-degraded areas are desertified and through this study, it is possible to revegetate and recover them. The second one is that, since Salicornia is a plant with economic value, this can serve as an incentive for farmers to grow Salicornia in saline areas.

To pool or not to pool: Can we ignore cross-trial variability in FMRI?

In this work, we investigate the importance of explicitly accounting for cross-trial variability in neuroimaging data analysis. To attempt to obtain reliable estimates in a task-based experiment, each condition is usually repeated across many trials. The investigator may be interested in (a) condition-level effects, (b) trial-level effects, or (c) the association of trial-level effects with the corresponding behavior data. The typical strategy for condition-level modeling is to create one regressor per condition at the subject level with the underlying assumption that responses do not change across trials. In this methodology of complete pooling, all cross-trial variability is ignored and dismissed as random noise that is swept under the rug of model residuals. Unfortunately, this framework invalidates the generalizability from the confine of specific trials (e.g., particular faces) to the associated stimulus category ("face"), and may inflate the statistical evidence when the trial sample size is not large enough. Here we propose an adaptive and computationally tractable framework that meshes well with the current two-level pipeline and explicitly accounts for trial-by-trial variability. The trial-level effects are first estimated per subject through no pooling. To allow generalizing beyond the particular stimulus set employed, the cross-trial variability is modeled at the population level through partial pooling in a multilevel model, which permits accurate effect estimation and characterization. Alternatively, trial-level estimates can be used to investigate, for example, brain-behavior associations or correlations between brain regions. Furthermore, our approach allows appropriate accounting for serial correlation, handling outliers, adapting to data skew, and capturing nonlinear brain-behavior relationships. By applying a Bayesian multilevel model framework at the level of regions of interest to an experimental dataset, we show how multiple testing can be addressed and full results reported without arbitrary dichotomization. Our approach revealed important differences compared to the conventional method at the condition level, including how the latter can distort effect magnitude and precision. Notably, in some cases our approach led to increased statistical sensitivity. In summary, our proposed framework provides an effective strategy to capture trial-by-trial responses that should be of interest to a wide community of experimentalists.

Representational Organization of Novel Task Sets during Proactive Encoding.

Recent multivariate analyses of brain data have boosted our understanding of the organizational principles that shape neural coding. However, most of this progress has focused on perceptual visual regions (Connolly et al., 2012), whereas far less is known about the organization of more abstract, action-oriented representations. In this study, we focused on humans' remarkable ability to turn novel instructions into actions. While previous research shows that instruction encoding is tightly linked to proactive activations in frontoparietal brain regions, little is known about the structure that orchestrates such anticipatory representation. We collected fMRI data while participants (both males and females) followed novel complex verbal rules that varied across control-related variables (integrating within/across stimuli dimensions, response complexity, target category) and reward expectations. Using representational similarity analysis (Kriegeskorte et al., 2008), we explored where in the brain these variables explained the organization of novel task encoding, and whether motivation modulated these representational spaces. Instruction representations in the lateral PFC were structured by the three control-related variables, whereas intraparietal sulcus encoded response complexity and the fusiform gyrus and precuneus organized its activity according to the relevant stimulus category. Reward exerted a general effect, increasing the representational similarity among different instructions, which was robustly correlated with behavioral improvements. Overall, our results highlight the flexibility of proactive task encoding, governed by distinct representational organizations in specific brain regions. They also stress the variability of motivation-control interactions, which appear to be highly dependent on task attributes, such as complexity or novelty.SIGNIFICANCE STATEMENT In comparison with other primates, humans display a remarkable success in novel task contexts thanks to our ability to transform instructions into effective actions. This skill is associated with proactive task-set reconfigurations in frontoparietal cortices. It remains yet unknown, however, how the brain encodes in anticipation the flexible, rich repertoire of novel tasks that we can achieve. Here we explored cognitive control and motivation-related variables that might orchestrate the representational space for novel instructions. Our results showed that different dimensions become relevant for task prospective encoding, depending on the brain region, and that the lateral PFC simultaneously organized task representations following different control-related variables. Motivation exerted a general modulation upon this process, diminishing rather than increasing distances among instruction representations.

Comparing functional connectivity matrices: A geometry-aware approach applied to participant identification.

Understanding the correlation structure associated with multiple brain measurements informs about potential "functional groupings" and network organization. The correlation structure can be conveniently captured in a matrix format that summarizes the relationships among a set of brain measurements involving two regions, for example. Such functional connectivity matrix is an important component of many types of investigation focusing on network-level properties of the brain, including clustering brain states, characterizing dynamic functional states, performing participant identification (so-called "fingerprinting") understanding how tasks reconfigure brain networks, and inter-subject correlation analysis. In these investigations, some notion of proximity or similarity of functional connectivity matrices is employed, such as their Euclidean distance or Pearson correlation (by correlating the matrix entries). Here, we propose the use of a geodesic distance metric that reflects the underlying non-Euclidean geometry of functional correlation matrices. The approach is evaluated in the context of participant identification (fingerprinting): given a participant's functional connectivity matrix based on resting-state or task data, how effectively can the participant be identified? Using geodesic distance, identification accuracy was over 95% on resting-state data, and exceeded the Pearson correlation approach by 20%. For whole-cortex regions, accuracy improved on a range of tasks by between 2% and as much as 20%. We also investigated identification using pairs of subnetworks (say, dorsal attention plus default mode), and particular combinations improved accuracy over whole-cortex participant identification by over 10%. The geodesic distance also outperformed Pearson correlation when the former employed a fourth of the data as the latter. Finally, we suggest that low-dimensional distance visualizations based on the geodesic approach help uncover the geometry of task functional connectivity in relation to that during resting-state. We propose that the use of the geodesic distance is an effective way to compare the correlation structure of the brain across a broad range of studies.

Fighting or embracing multiplicity in neuroimaging? neighborhood leverage versus global calibration.

Neuroimaging faces the daunting challenge of multiple testing - an instance of multiplicity - that is associated with two other issues to some extent: low inference efficiency and poor reproducibility. Typically, the same statistical model is applied to each spatial unit independently in the approach of massively univariate modeling. In dealing with multiplicity, the general strategy employed in the field is the same regardless of the specifics: trust the local "unbiased" effect estimates while adjusting the extent of statistical evidence at the global level. However, in this approach, modeling efficiency is compromised because each spatial unit (e.g., voxel, region, matrix element) is treated as an isolated and independent entity during massively univariate modeling. In addition, the required step of multiple testing "correction" by taking into consideration spatial relatedness, or neighborhood leverage, can only partly recoup statistical efficiency, resulting in potentially excessive penalization as well as arbitrariness due to thresholding procedures. Moreover, the assigned statistical evidence at the global level heavily relies on the data space (whole brain or a small volume). The present paper reviews how Stein's paradox (1956) motivates a Bayesian multilevel (BML) approach that, rather than fighting multiplicity, embraces it to our advantage through a global calibration process among spatial units. Global calibration is accomplished via a Gaussian distribution for the cross-region effects whose properties are not a priori specified, but a posteriori determined by the data at hand through the BML model. Our framework therefore incorporates multiplicity as integral to the modeling structure, not a separate correction step. By turning multiplicity into a strength, we aim to achieve five goals: 1) improve the model efficiency with a higher predictive accuracy, 2) control the errors of incorrect magnitude and incorrect sign, 3) validate each model relative to competing candidates, 4) reduce the reliance and sensitivity on the choice of data space, and 5) encourage full results reporting. Our modeling proposal reverberates with recent proposals to eliminate the dichotomization of statistical evidence ("significant" vs. "non-significant"), to improve the interpretability of study findings, as well as to promote reporting the full gamut of results (not only "significant" ones), thereby enhancing research transparency and reproducibility.

Interactions between emotion and action in the brain.

A growing literature supports the existence of interactions between emotion and action in the brain, and the central participation of the anterior midcingulate cortex (aMCC) in this regard. In the present functional magnetic resonance imaging study, we sought to investigate the role of self-relevance during such interactions by varying the context in which threating pictures were presented (with guns pointed towards or away from the observer). Participants performed a simple visual detection task following exposure to such stimuli. Except for voxelwise tests, we adopted a Bayesian analysis framework which evaluated evidence for the hypotheses of interest, given the data, in a continuous fashion. Behaviorally, our results demonstrated a valence by context interaction such that there was a tendency of speeding up responses to targets after viewing threat pictures directed towards the participant. In the brain, interaction patterns that paralleled those observed behaviorally were observed most notably in the middle temporal gyrus, supplementary motor area, precentral gyrus, and anterior insula. In these regions, activity was overall greater during threat conditions relative to neutral ones, and this effect was enhanced in the directed towards context. A valence by context interaction was observed in the aMCC too, where we also observed a correlation (across participants) of evoked responses and reaction time data. Taken together, our study revealed the context-sensitive engagement of motor-related areas during emotional perception, thus supporting the idea that emotion and action interact in important ways in the brain.

Do applications of soil conditioner mixtures improve the salt extraction ability of Atriplex nummularia at early growth stage?

The reclamation of salt-affected soils is considered a slow process that urges the development of faster recovery strategies as a priority. The present article aimed at investigating the effect of mixing different chemical and organic conditioners on the growth of Atriplex and its salt extraction efficiency during its early growth stage. A pot experiment was conducted on saline-sodic Cambisol cultivated with Atriplex for 60 days and subjected to the following conditioner mixtures: Atriplex; Atriplex + gypsum + organic matter; Atriplex + elemental sulfur + organic matter; and Atriplex + gypsum + elemental sulfur + organic matter. The mixtures of conditioners did not affect the Atriplex growth, but caused significant increase in Na+ and Cl- contents in the dry matter of Atriplex plants. Additionally, there was an increase in the Atriplex's ability of extracting salt due to the application of the mixtures. Results suggest that the "gypsum + organic matter" mixture is preferable for a faster recovery of salts/sodium affected soils through phytoremediation by Atriplex plants, mainly due to a more significant increase in the efficiency of salt absorption.

Dynamic Threat Processing.

During real-life situations, multiple factors interact dynamically to determine threat level. In the current fMRI study involving healthy adult human volunteers, we investigated interactions between proximity, direction (approach vs. retreat), and speed during a dynamic threat-of-shock paradigm. As a measure of threat-evoked physiological arousal, skin conductance responses were recorded during fMRI scanning. Some brain regions tracked individual threat-related factors, and others were also sensitive to combinations of these variables. In particular, signals in the anterior insula tracked the interaction between proximity and direction where approach versus retreat responses were stronger when threat was closer compared with farther. A parallel proximity-by-direction interaction was also observed in physiological skin conductance responses. In the right amygdala, we observed a proximity by direction interaction, but intriguingly in the opposite direction as the anterior insula; retreat versus approach responses were stronger when threat was closer compared with farther. In the right bed nucleus of the stria terminalis, we observed an effect of threat proximity, whereas in the right periaqueductal gray/midbrain we observed an effect of threat direction and a proximity by direction by speed interaction (the latter was detected in exploratory analyses but not in a voxelwise fashion). Together, our study refines our understanding of the brain mechanisms involved during aversive anticipation in the human brain. Importantly, it emphasizes that threat processing should be understood in a manner that is both context-sensitive and dynamic.

Brain dynamics and temporal trajectories during task and naturalistic processing.

Human functional Magnetic Resonance Imaging (fMRI) data are acquired while participants engage in diverse perceptual, motor, cognitive, and emotional tasks. Although data are acquired temporally, they are most often treated in a quasi-static manner. Yet, a fuller understanding of the mechanisms that support mental functions necessitates the characterization of dynamic properties. Here, we describe an approach employing a class of recurrent neural networks called reservoir computing, and show the feasibility and potential of using it for the analysis of temporal properties of brain data. We show that reservoirs can be used effectively both for condition classification and for characterizing lower-dimensional "trajectories" of temporal data. Classification accuracy was approximately 90% for short clips of "social interactions" and around 70% for clips extracted from movie segments. Data representations with 12 or fewer dimensions (from an original space with over 300) attained classification accuracy within 5% of the full data. We hypothesize that such low-dimensional trajectories may provide "signatures" that can be associated with tasks and/or mental states. The approach was applied across participants (that is, training in one set of participants, and testing in a separate group), showing that representations generalized well to unseen participants. Taken together, we believe the present approach provides a promising framework to characterize dynamic fMRI information during both tasks and naturalistic conditions.

An integrative Bayesian approach to matrix-based analysis in neuroimaging.

Understanding the correlation structure associated with brain regions is a central goal in neuroscience, as it informs about interregional relationships and network organization. Correlation structure can be conveniently captured in a matrix that indicates the relationships among brain regions, which could involve electroencephalogram sensors, electrophysiology recordings, calcium imaging data, or functional magnetic resonance imaging (FMRI) data-We call this type of analysis matrix-based analysis, or MBA. Although different methods have been developed to summarize such matrices across subjects, including univariate general linear models (GLMs), the available modeling strategies tend to disregard the interrelationships among the regions, leading to "inefficient" statistical inference. Here, we develop a Bayesian multilevel (BML) modeling framework that simultaneously integrates the analyses of all regions, region pairs (RPs), and subjects. In this approach, the intricate relationships across regions as well as across RPs are quantitatively characterized. The adoption of the Bayesian framework allows us to achieve three goals: (a) dissolve the multiple testing issue typically associated with seeking evidence for the effect of each RP under the conventional univariate GLM; (b) make inferences on effects that would be treated as "random" under the conventional linear mixed-effects framework; and (c) estimate the effect of each brain region in a manner that indexes their relative "importance". We demonstrate the BML methodology with an FMRI dataset involving a cognitive-emotional task and compare it to the conventional GLM approach in terms of model efficiency, performance, and inferences. The associated program MBA is available as part of the AFNI suite for general use.

Embracing integration and complexity: placing emotion within a science of brain and behaviour.

The present paper addresses conceptual issues that are central to emotion research. What is emotion? What are its defining characteristics? The field struggles with questions like these almost constantly. I argue that definitions, and deciding what is the proper status of emotion, are not a requirement for scientific progress - in fact, they can hinder it. Therefore, "emotion" researchers should strive to develop a science of complex behaviours, and worry less about their exact nature. But for interesting behaviours, is most of the explaining that is needed present at the level of isolated systems (perception, cognition, etc.) or at the level of interactions between them? I suggest that the level of interactions is where most of the work is needed. Accordingly, I advocate that it is important to embrace integration, and not to strive to necessarily disentangle the multiple contributions underlying behaviours. More generally, it is argued that we need to revise models of causation adopted when reasoning about the mind and brain. Instead, a "complex systems" approach is required where the interactions between multiple components lead to system-level - emergent - properties that cannot be isolated or attributed to more elementary parts.

Screening Aedes aegypti (Diptera: Culicidae) Populations From Pernambuco, Brazil for Resistance to Temephos, Diflubenzuron, and Cypermethrin and Characterization of Potential Resistance Mechanisms.

Resistance to chemical insecticides detected in Aedes aegypti (L.) mosquitoes has been a problem for the National Dengue Control Program (PNCD) over the last years. In order to provide deeper knowledge of resistance to xenobiotics, our study evaluated the susceptibility profile of temephos, diflubenzuron, and cypermethrin insecticides in natural mosquito populations from the Pernambuco State, associating these results with the local historical use of such compounds. Furthermore, mechanisms that may be associated with this particular type of resistance were characterized. Bioassays with multiple temephos and diflubenzuron concentrations were performed to detect and quantify resistance. For cypermethrin, diagnostic dose assays were performed. Biochemical tests were carried out to quantify the activity of detoxification enzymes. In addition, a screening of mutations present in the voltage-gated sodium channel gene (NaV) was performed in samples previously submitted to bioassays with cypermethrin. The populations under study were resistant to temephos and showed a positive correlation between insecticide consumption and the resistance ratio (RR) to the compound. For diflubenzuron, the biological activity ratio (BAR) ranged from 1.3 to 4.7 times, when compared to the susceptible strain. All populations showed resistance to cypermethrin. Altered enzymatic profiles of alpha, p-nitrophenyl acetate (PNPA) esterases and glutathione-S-transferases were recorded in most of these samples. Molecular analysis demonstrated that Arcoverde was the only population that presented the mutated form 1016Ile/Ile. These findings show that the situation is critical vis-à-vis the effectiveness of mosquito control using chemical insecticides, since resistance to temephos and cypermethrin is widespread in Ae. aegypti from Pernambuco.

Brain networks for emotion and cognition: Implications and tools for understanding mental disorders and pathophysiology.

Understanding how structure maps to function in the brain in terms of large-scale networks is critical to elucidating the brain basis of mental phenomena and mental disorders. Given that this mapping is many-to-many, I argue that researchers need to shift to a multivariate brain and behavior characterization to fully unravel the contributions of brain processes to typical and atypical function.

LEICA: Laplacian eigenmaps for group ICA decomposition of fMRI data.

Independent component analysis (ICA) is a data-driven method that has been increasingly used for analyzing functional Magnetic Resonance Imaging (fMRI) data. However, generalizing ICA to multi-subject studies is non-trivial due to the high-dimensionality of the data, the complexity of the underlying neuronal processes, the presence of various noise sources, and inter-subject variability. Current group ICA based approaches typically use several forms of the Principal Component Analysis (PCA) method to extend ICA for generating group inferences. However, linear dimensionality reduction techniques have serious limitations including the fact that the underlying BOLD signal is a complex function of several nonlinear processes. In this paper, we propose an effective non-linear ICA-based model for extracting group-level spatial maps from multi-subject fMRI datasets. We use a non-linear dimensionality reduction algorithm based on Laplacian eigenmaps to identify a manifold subspace common to the group, such that this mapping preserves the correlation among voxels' time series as much as possible. These eigenmaps are modeled as linear mixtures of a set of group-level spatial features, which are then extracted using ICA. The resulting algorithm is called LEICA (Laplacian Eigenmaps for group ICA decomposition). We introduce a number of methods to evaluate LEICA using 100-subject resting state and 100-subject working memory task fMRI datasets from the Human Connectome Project (HCP). The test results show that the extracted spatial maps from LEICA are meaningful functional networks similar to those produced by some of the best known methods. Importantly, relative to state-of-the-art methods, our algorithm compares favorably in terms of the functional cohesiveness of the spatial maps generated, as well as in terms of the reproducibility of the results.