RESUMO
Real-time brain functional MRI (rt-fMRI) allows in vivo non-invasive monitoring of neural networks. The use of multivariate data-driven analysis methods such as independent component analysis (ICA) offers an attractive trade-off between data interpretability and information extraction, and can be used during both task-based and rest experiments. The purpose of this study was to assess the effectiveness of different ICA-based procedures to monitor in real-time a target IC defined from a functional localizer which also used ICA. Four novel methods were implemented to monitor ongoing brain activity in a sliding window approach. The methods differed in the ways in which a priori information, derived from ICA algorithms, was used to monitor a target independent component (IC). We implemented four different algorithms, all based on ICA. One Back-projection method used ICA to derive static spatial information from the functional localizer, off-line, which was then back-projected dynamically during the real-time acquisition. The other three methods used real-time ICA algorithms that dynamically exploited temporal, spatial, or spatial-temporal priors during the real-time acquisition. The methods were evaluated by simulating a rt-fMRI experiment that used real fMRI data. The performance of each method was characterized by the spatial and/or temporal correlation with the target IC component monitored, computation time, and intrinsic stochastic variability of the algorithms. In this study the Back-projection method, which could monitor more than one IC of interest, outperformed the other methods. These results are consistent with a functional task that gives stable target ICs over time. The dynamic adaptation possibilities offered by the other ICA methods proposed may offer better performance than the Back-projection in conditions where the functional activation shows higher spatial and/or temporal variability.
RESUMO
Independent component analysis (ICA) techniques offer a data-driven possibility to analyze brain functional MRI data in real-time. Typical ICA methods used in functional magnetic resonance imaging (fMRI), however, have been until now mostly developed and optimized for the off-line case in which all data is available. Real-time experiments are ill-posed for ICA in that several constraints are added: limited data, limited analysis time and dynamic changes in the data and computational speed. Previous studies have shown that particular choices of ICA parameters can be used to monitor real-time fMRI (rt-fMRI) brain activation, but it is unknown how other choices would perform. In this rt-fMRI simulation study we investigate and compare the performance of 14 different publicly available ICA algorithms systematically sampling different growing window lengths (WLs), model order (MO) as well as a priori conditions (none, spatial or temporal). Performance is evaluated by computing the spatial and temporal correlation to a target component as well as computation time. Four algorithms are identified as best performing (constrained ICA, fastICA, amuse, and evd), with their corresponding parameter choices. Both spatial and temporal priors are found to provide equal or improved performances in similarity to the target compared with their off-line counterpart, with greatly reduced computation costs. This study suggests parameter choices that can be further investigated in a sliding-window approach for a rt-fMRI experiment.
RESUMO
BACKGROUND: In the absence of overt stimuli, the brain shows correlated fluctuations in functionally related brain regions. Approximately ten largely independent resting state networks (RSNs) showing this behaviour have been documented to date. Recent studies have reported the existence of an RSN in the basal ganglia - albeit inconsistently and without the means to interpret its function. Using two large study groups with different resting state conditions and MR protocols, the reproducibility of the network across subjects, behavioural conditions and acquisition parameters is assessed. Independent Component Analysis (ICA), combined with novel analyses of temporal features, is applied to establish the basis of signal fluctuations in the network and its relation to other RSNs. Reference to prior probabilistic diffusion tractography work is used to identify the basal ganglia circuit to which these fluctuations correspond. RESULTS: An RSN is identified in the basal ganglia and thalamus, comprising the pallidum, putamen, subthalamic nucleus and substantia nigra, with a projection also to the supplementary motor area. Participating nuclei and thalamo-cortical connection probabilities allow this network to be identified as the motor control circuit of the basal ganglia. The network was reproducibly identified across subjects, behavioural conditions (fixation, eyes closed), field strength and echo-planar imaging parameters. It shows a frequency peak at 0.025 +/- 0.007 Hz and is most similar in spectral composition to the Default Mode (DM), a network of regions that is more active at rest than during task processing. Frequency features allow the network to be classified as an RSN rather than a physiological artefact. Fluctuations in this RSN are correlated with those in the task-positive fronto-parietal network and anticorrelated with those in the DM, whose hemodynamic response it anticipates. CONCLUSION: Although the basal ganglia RSN has not been reported in most ICA-based studies using a similar methodology, we demonstrate that it is reproducible across subjects, common resting state conditions and imaging parameters, and show that it corresponds with the motor control circuit. This characterisation of the basal ganglia network opens a potential means to investigate the motor-related neuropathologies in which the basal ganglia are involved.