Generalized relative entropy in functional magnetic resonance imaging

SummaryHumans show a remarkable ability to discriminate others' gaze direction, even though a given direction can be conveyed by many physically dissimilar configurations of different eye positions and head views. For example, eye contact can be signaled by a rightward glance in a left-turned head or by direct gaze in a front-facing head. Such acute gaze discrimination implies considerable perceptual invariance. Previous human research found that superior temporal sulcus (STS) responds preferentially to gaze shifts [1], but the underlying representation that supports such general responsiveness remains poorly understood. Using multivariate pattern analysis (MVPA) of human functional magnetic resonance imaging (fMRI) data, we tested whether STS contains a higher-order, head view-invariant code for gaze direction. The results revealed a finely graded gaze direction code in right anterior STS that was invariant to head view and physical image features. Further analyses revealed similar gaze effects in left anterior STS and precuneus. Our results suggest that anterior STS codes the direction of another's attention regardless of how this information is conveyed and demonstrate how high-level face areas carry out fine-grained, perceptually relevant discrimination through invariance to other face features.

A new adaptive signal-preserving technique for noise suppression in event-related functional magnetic resonance imaging (fMRI) data is proposed based on spectral subtraction. The proposed technique estimates a parametric model for the power spectrum of random noise from the acquired data based on the characteristics of the Rician statistical model. This model is subsequently used to estimate a noise-suppressed power spectrum for any given pixel time course by simple subtraction of power spectra. The new technique is tested using computer simulations and real data from event-related fMRI experiments. The results show the potential of the new technique in suppressing noise while preserving the other deterministic components in the signal. Moreover, we demonstrate that further analysis using principal component analysis and independent component analysis shows a significant improvement in both convergence and clarity of results when the new technique is used. Given its simple form, the new method does not change the statistical characteristics of the signal or cause correlated noise to be present in the processed signal. This suggests the value of the new technique as a useful preprocessing step for fMRI data analysis.

Whether signals from different sensory modalities converge and interact within primary cortices in humans is unresolved, despite emerging evidence in animals. This is partially because of debates concerning the appropriate analyses of functional magnetic resonance imaging (fMRI) data in response to multisensory phenomena. Using event-related fMRI, we observed that simple auditory stimuli (noise bursts) activated primary visual cortices and that simple visual stimuli (checkerboards) activated primary auditory cortices, indicative of multisensory convergence. Moreover, analyses of blood oxygen level-dependent response dynamics revealed facilitation of hemodynamic response peak latencies and slopes for multisensory auditory-visual stimuli versus either unisensory condition, indicative of multisensory interactions within primary sensory cortices. Neural processing at the lowest cortical levels can be modulated by interactions between the senses. Temporal information in fMRI data can reveal these modulations and overcome analytic and interpretational challenges of more traditional procedures. In addition to providing an essential translational link with animal models, these results suggest that longstanding notions of cortical organization need to be revised to include multisensory interactions as an inherent component of functional brain organization.

Previous studies have shown that hemodynamic response overlap severely limits the maximum presentation rate with event-related functional MRI (fMRI) using fixed intertrial experimental designs. Here we demonstrate that the use of randomized experimental designs can largely overcome this limitation, thereby allowing for event-related fMRI experiments with extremely rapid presentation rates. In the first experiment, fMRI time courses were simulated using a fixed intertrial interval design with intervals of 16, 3, and 1 s, and using a randomized design having the same mean intertrial intervals. We found that using fixed intertrial interval designs the transient information decreased with decreasing intertrial intervals, whereas using randomized designs the transient information increased with decreasing mean intertrial intervals. In a second experiment, fMRI data were collected from two subjects using a randomized paradigm with visual hemifield stimuli presented randomly every 500 ms. Robust event-related activation maps and hemodynamic response estimates were obtained. These results demonstrate the feasibility of performing event-related fMRI experiments with rapid, randomized paradigms identical to those used in electrophysiological and behavioral studies, thereby expanding the applicability of event-related fMRI to a whole new range of cognitive neurosciences questions and paradigms.

Background We are considering the statistical analysis of functional magnetic resonance imaging (fMRI) data. As demonstrated in previous work, grouping voxels into regions (of interest) and carrying out a multiple test for signal detection on the basis of these regions typically leads to a higher sensitivity when compared with voxel-wise multiple testing approaches. Methods In the case of a multi-subject study, we propose to define the regions for each subject separately based on their individual brain anatomy, represented, e.g., by regional labels. The aggregation of the subject-specific evidence for the presence of signals in the different regions is then performed by means of a combination function for p -values. We validate the proposed methodology with simulated data and apply it to real fMRI data of a hypothesis-driven approach towards identifying brain regions involved in understanding software code. Results The results of our simulated data indicate that the proposed approach improves power relative to voxel-wise inference and performs comparably to a cluster-based baseline. Testing our method on real fMRI data, we found that our approach yields overlapping results with a two-stage approach for which two independent experiments are needed, one for defining the regions and one for actual signal detection. Conclusions In this paper, we overall demonstrate that our method of utilizing anatomical information is a candidate to provide a more sensitive analysis of fMRI data.

P300 recordings during event-related fMRI: a feasibility study

Investigating the enhancement of template-free activation detection of event-related fMRI data using wavelet shrinkage and figures of merit

Independent component analysis (ICA) has proved to be a powerful method for exploratory analysis of functional magnetic resonance imaging (fMRI) data. It has been used to uncover unexpected activations in fMRI data derived from brain activation. ICA has been used to characterize other sources of variability in the fMRI signal besides task-related activity, as well as challenging some of the assumptions inherent in other fMRI analysis methods. As a data-driven fMRI analysis technique, the philosophy of ICA is often in disagreement with hypothesis-driven methods. By exploiting the fact that much of fMRI data has deterministic spatial-temporal structure, a scheme employing ICA denoising and least squares (LS) estimation of the evoked hemodynamic response (HDR) is proposed. Simulations suggest that the method is more robust to different noise models compared to naive application of LS. The result is a considerably increased level of significance of activation for a given voxel but still qualitatively similar spatial distribution of activations over all voxels. This suggests that the proposed method has the potential to substantially reduce total scanning time requirements to achieve the same level of statistically significant activation.

In recent years, the discovery of brain effective connectivity (EC) networks through computational analysis of functional magnetic resonance imaging (fMRI) data has gained prominence in neuroscience and neuroimaging. However, owing to the influence of diverse factors during data collection and processing, fMRI data typically exhibits high noise and limited sample characteristics, consequently leading to suboptimal performance of current methods. In this paper, we propose a novel brain effective connectivity discovery method based on meta-reinforcement learning, called MetaRLEC. The method mainly consists of three modules: actor, critic, and meta-critic. MetaRLEC first employs an encoder-decoder framework: the encoder utilizing a Transformer, converts noisy fMRI data into a state embedding; the decoder employing bidirectional LSTM, discovers brain region dependencies from the state and generates actions (EC networks). Then a critic network evaluates these actions, incentivizing the actor to learn higher-reward actions amidst the high-noise setting. Finally, a meta-critic framework facilitates online learning of historical state-action pairs, integrating an action-value neural network and supplementary training losses to enhance the model's adaptability to small-sample fMRI data. We conduct comprehensive experiments on both simulated and real-world data to demonstrate the efficacy of our proposed method.

Analysis of time-evolving data is crucial to understand the functioning of dynamic systems such as the brain. For instance, analysis of functional magnetic resonance imaging (fMRI) data collected during a task may reveal spatial regions of interest, and how they evolve during the task. However, capturing underlying spatial patterns as well as their change in time is challenging. The traditional approach in fMRI data analysis is to assume that underlying spatial regions of interest are static. In this article, using fractional amplitude of low-frequency fluctuations (fALFF) as an effective way to summarize the variability in fMRI data collected during a task, we arrange time-evolving fMRI data as a subjects by voxels by time windows tensor, and analyze the tensor using a tensor factorization-based approach called a PARAFAC2 model to reveal spatial dynamics. The PARAFAC2 model jointly analyzes data from multiple time windows revealing subject-mode patterns, evolving spatial regions (also referred to as networks) and temporal patterns. We compare the PARAFAC2 model with matrix factorization-based approaches relying on independent components, namely, joint independent component analysis (ICA) and independent vector analysis (IVA), commonly used in neuroimaging data analysis. We assess the performance of the methods in terms of capturing evolving networks through extensive numerical experiments demonstrating their modeling assumptions. In particular, we show that (i) PARAFAC2 provides a compact representation in all modes, i.e., subjects, time, and voxels, revealing temporal patterns as well as evolving spatial networks, (ii) joint ICA is as effective as PARAFAC2 in terms of revealing evolving networks but does not reveal temporal patterns, (iii) IVA's performance depends on sample size, data distribution and covariance structure of underlying networks. When these assumptions are satisfied, IVA is as accurate as the other methods, (iv) when subject-mode patterns differ from one time window to another, IVA is the most accurate. Furthermore, we analyze real fMRI data collected during a sensory motor task, and demonstrate that a component indicating statistically significant group difference between patients with schizophrenia and healthy controls is captured, which includes primary and secondary motor regions, cerebellum, and temporal lobe, revealing a meaningful spatial map and its temporal change.

A new adaptive signal-preserving technique for noise suppression in functional magnetic resonance imaging (fMRI) data is proposed based on spectrum subtraction. The proposed technique estimates a model for the power spectrum of random noise from the acquired data. This model is used to estimate a noise-suppressed power spectrum for any given pixel time course by simple subtraction of power spectra. The new technique is tested using computer simulations and real data for event-related fMRI experiments. The results show the potential of the new technique in suppressing noise while preserving the other deterministic components. Moreover, further analysis using principal component analysis (PCA) and independent component analysis (ICA) shows a significant improvement in both convergence and clarity of results when the new technique is used. This suggests the value of the new technique as a useful preprocessing step for this type of signal.

Function magnetic resonance imaging (fMRI) data are typically contaminated by noise introduced by head motion, physiological noise, and thermal noise. To mitigate noise artifact in fMRI data, a variety of denoising methods have been developed by removing noise factors derived from the whole time series of fMRI data and therefore are not applicable to real-time fMRI data analysis. In the present study, we develop a generally applicable, deep learning based fMRI denoising method to generate noise-free realistic individual fMRI volumes (time points). Particularly, we develop a fully data-driven 3D convolutional encapsulated Long Short-Term Memory (3DConv-LSTM) approach to generate noise-free fMRI volumes regularized by an adversarial network that makes the generated fMRI volumes more realistic by fooling a critic network. The 3DConv-LSTM model also integrates a gate-controlled self-attention model to memorize short-term dependency and historical information within a memory pool. We have evaluated our method based on both task and resting state fMRI data. Both qualitative and quantitative results have demonstrated that the proposed method outperformed state-of-the-art alternative deep learning methods.

Separating 4D multi-task fMRI data of multiple subjects by independent component analysis with projection

A group comparison in fMRI data using a semiparametric model under shape invariance

Analysis of functional magnetic resonance imaging (fMRI) data has been performed using both model-driven (parametric) methods and data-driven methods. An advantage of model-driven methods is incorporation of prior knowledge of spatial and temporal properties of the hemodynamic response (HDR). A novel analytical framework for fMRI data has been developed that identifies multi-voxel regions of activation through iterative segmentation-based optimization over HDR estimates for both individual voxels and regional groupings. Simulations using synthetic activation embedded in autoregressive integrated moving average (ARIMA) noise reveal the proposed procedure to be more sensitive and selective than conventional fMRI analysis methods (reference set: principle component analysis, PCA; independent component analysis, ICA; k-means clustering, k=100; univariate t-test) in identification of active regions over the range of average contrast-to-noise ratios of 0.5 to 4.0. Results of analysis of extant human data (for which the average contrast-to-noise ratio is unknown) are further suggestive of greater statistical detection power. Refinement of this new procedure is expected to reduce both false positive and negative rates, without resorting to filtering that can reduce the effective spatial resolution.