Evaluation of spatio-temporal decomposition techniques for group analysis of fMRI resting state data sets

Denoising brain networks using a fixed mathematical phase change in independent component analysis of magnitude-only fMRI data

Brain networks extracted by independent component analysis (ICA) from magnitude-only fMRI data are usually denoised using various amplitude-based thresholds. By contrast, spatial source phase (SSP) or the phase information of ICA brain networks extracted from complex-valued fMRI data, has provided a simple yet effective way to perform the denoising using a fixed phase change. In this work, we extend the approach to magnitude-only fMRI data to avoid testing various amplitude thresholds for denoising magnitude maps extracted by ICA, as most studies do not save the complex-valued data. The main idea is to generate a mathematical SSP map for a magnitude map using a mapping framework, and the mapping framework is built using complex-valued fMRI data with a known SSP map. Here we leverage the fact that the phase map derived from phase fMRI data has similar phase information to the SSP map. After verifying the use of the magnitude data of complex-valued fMRI, this framework is generalized to work with magnitude-only data, allowing use of our approach even without the availability of the corresponding phase fMRI datasets. We test the proposed method using both simulated and experimental fMRI data including complex-valued data from University of New Mexico and magnitude-only data from Human Connectome Project. The results provide evidence that the mathematical SSP denoising with a fixed phase change is effective for denoising spatial maps from magnitude-only fMRI data in terms of retaining more BOLD-related activity and fewer unwanted voxels, compared with amplitude-based thresholding. The proposed method provides a unified and efficient SSP approach to denoise ICA brain networks in fMRI data.

A technical review of canonical correlation analysis for neuroscience applications

Collecting comprehensive data sets of the same subject has become a standard in neuroscience research and uncovering multivariate relationships among collected data sets have gained significant attentions in recent years. Canonical correlation analysis (CCA) is one of the powerful multivariate tools to jointly investigate relationships among multiple data sets, which can uncover disease or environmental effects in various modalities simultaneously and characterize changes during development, aging, and disease progressions comprehensively. In the past 10 years, despite an increasing number of studies have utilized CCA in multivariate analysis, simple conventional CCA dominates these applications. Multiple CCA-variant techniques have been proposed to improve the model performance; however, the complicated multivariate formulations and not well-known capabilities have delayed their wide applications. Therefore, in this study, a comprehensive review of CCA and its variant techniques is provided. Detailed technical formulation with analytical and numerical solutions, current applications in neuroscience research, and advantages and limitations of each CCA-related technique are discussed. Finally, a general guideline in how to select the most appropriate CCA-related technique based on the properties of available data sets and particularly targeted neuroscience questions is provided.

Pretherapeutic resting-state fMRI profiles are associated with MR signature volumes after stereotactic radiosurgical thalamotomy for essential tremor

OBJECTIVEEssential tremor (ET) is the most common movement disorder. Drug-resistant ET can benefit from standard stereotactic deep brain stimulation or radiofrequency thalamotomy or, alternatively, minimally invasive techniques, including stereotactic radiosurgery (SRS) and high-intensity focused ultrasound, at the level of the ventral intermediate nucleus (Vim). The aim of the present study was to evaluate potential correlations between pretherapeutic interconnectivity (IC), as depicted on resting-state functional MRI (rs-fMRI), and MR signature volume at 1 year after Vim SRS for tremor, to be able to potentially identify hypo- and hyperresponders based only on pretherapeutic neuroimaging data.METHODSSeventeen consecutive patients with ET were included, who benefitted from left unilateral SRS thalamotomy (SRS-T) between September 2014 and August 2015. Standard tremor assessment and rs-fMRI were acquired pretherapeutically and 1 year after SRS-T. A healthy control group was also included (n = 12). Group-level independent component analysis (ICA; only n = 17 for pretherapeutic rs-fMRI) was applied. The mean MR signature volume was 0.125 ml (median 0.063 ml, range 0.002-0.600 ml). The authors correlated baseline IC with 1-year MR signatures within all networks. A 2-sample t-test at the level of each component was first performed in two groups: group 1 (n = 8, volume < 0.063 ml) and group 2 (n = 9, volume ≥ 0.063 ml). These groups did not statistically differ by age, duration of symptoms, baseline ADL score, ADL point decrease at 1 year, time to tremor arrest, or baseline tremor score on the treated hand (TSTH; p > 0.05). An ANOVA was then performed on each component, using individual subject-level maps and continuous values of 1-year MR signatures, correlated with pretherapeutic IC.RESULTSUsing 2-sample t-tests, two networks were found to be statistically significant: network 3, including the brainstem, motor cerebellum, bilateral thalamus, and left supplementary motor area (SMA) (pFWE = 0.004, cluster size = 94), interconnected with the red nucleus (MNI -2, -22, -32); and network 9, including the brainstem, posterior insula, bilateral thalamus, and left SMA (pFWE = 0.002, cluster size = 106), interconnected with the left SMA (MNI 24, -28, 44). Higher pretherapeutic IC was associated with higher MR volumes, in a network including the anterior default-mode network and bilateral thalamus (ANOVA, pFWE = 0.004, cluster size = 73), interconnected with cerebellar lobule V (MNI -12, -70, -22). Moreover, in the same network, radiological hyporesponders presented with negative IC values.CONCLUSIONSThese findings have clinical implications for predicting MR signature volumes after SRS-T. Here, using pretherapeutic MRI and data processing without prior hypothesis, the authors showed that pretherapeutic network interconnectivity strength predicts 1-year MR signature volumes following SRS-T.

Multiway canonical correlation analysis of brain data

Brain data recorded with electroencephalography (EEG), magnetoencephalography (MEG) and related techniques often have poor signal-to-noise ratios due to the presence of multiple competing sources and artifacts. A common remedy is to average responses over repeats of the same stimulus, but this is not applicable for temporally extended stimuli that are presented only once (speech, music, movies, natural sound). An alternative is to average responses over multiple subjects that were presented with identical stimuli, but differences in geometry of brain sources and sensors reduce the effectiveness of this solution. Multiway canonical correlation analysis (MCCA) brings a solution to this problem by allowing data from multiple subjects to be fused in such a way as to extract components common to all. This paper reviews the method, offers application examples that illustrate its effectiveness, and outlines the caveats and risks entailed by the method.

Automated iterative reclustering framework for determining hierarchical functional networks in resting state fMRI

To spatially cluster resting state-functional magnetic resonance imaging (rs-fMRI) data into potential networks, there are only a few general approaches that determine the number of networks/clusters, despite a wide variety of techniques proposed for clustering. For individual subjects, extraction of a large number of spatially disjoint clusters results in multiple small networks that are spatio-temporally homogeneous but irreproducible across subjects. Alternatively, extraction of a small number of clusters creates spatially large networks that are temporally heterogeneous but spatially reproducible across subjects. We propose a fully automatic, iterative reclustering framework in which a small number of spatially large, heterogeneous networks are initially extracted to maximize spatial reproducibility. Subsequently, the large networks are iteratively subdivided to create spatially reproducible subnetworks until the overall within-network homogeneity does not increase substantially. The proposed approach discovers a rich network hierarchy in the brain while simultaneously optimizing spatial reproducibility of networks across subjects and individual network homogeneity. We also propose a novel metric to measure the connectivity of brain regions, and in a simulation study show that our connectivity metric and framework perform well in the face of low signal to noise and initial segmentation errors. Experimental results generated using real fMRI data show that the proposed metric improves stability of network clusters across subjects, and generates a meaningful pattern for spatially hierarchical structure of the brain.

A Hybrid LDA+gCCA Model for fMRI Data Classification and Visualization

Linear predictive models are applied to functional MRI (fMRI) data to estimate boundaries that predict experimental task states for scans. These boundaries are visualized as statistical parametric maps (SPMs) and range from low to high spatial reproducibility across subjects (e.g., Strother , 2004; LaConte , 2003). Such inter-subject pattern reproducibility is an essential characteristic of interpretable SPMs that generalize across subjects. Therefore, we introduce a flexible hybrid model that optimizes reproducibility by simultaneously enhancing the prediction power and reproducibility. This hybrid model is formed by a weighted summation of the optimization functions of a linear discriminate analysis (LDA) model and a generalized canonical correlation (gCCA) model (Afshin-Pour , 2012). LDA preserves the model's ability to discriminate the fMRI scans of multiple brain states while gCCA finds a linear combination for each subject's scans such that the estimated boundary map is reproducible. The hybrid model is implemented in a split-half resampling framework (Strother , 2010) which provides reproducibility (r) and prediction (p) quality metrics. Then the model was compared with LDA, and Gaussian Naive Bayes (GNB). For simulated fMRI data, the hybrid model outperforms the other two techniques in terms of receiver operating characteristic (ROC) curves, particularly for detecting less predictable but spatially reproducible networks. These techniques were applied to real fMRI data to estimate the maps for two task contrasts. Our results indicate that compared to LDA and GNB, the hybrid model can provide maps with large increases in reproducibility for small reductions in prediction, which are jointly closer to the ideal performance point of (p=1, r=1).

Whole-brain connectivity analysis and classification of spinocerebellar ataxia type 7 by functional MRI

Background

Spinocerebellar ataxia type 7 (SCA7) is a genetic disorder characterized by degeneration of the motor and visual systems. Besides neural deterioration, these patients also show functional connectivity changes linked to the degenerated brain areas. However, it is not known if there are functional connectivity changes in regions not necessarily linked to the areas undergoing structural deterioration. Therefore, in this study we have explored the whole-brain functional connectivity of SCA7 patients in order to find the overall abnormal functional pattern of this disease. Twenty-six patients and age-and-gender-matched healthy controls were recruited. Whole-brain functional connectivity analysis was performed in both groups. A classification algorithm was used to find the discriminative power of the abnormal connections by classifying patients and healthy subjects.

Results

Nineteen abnormal functional connections involving cerebellar and cerebral regions were selected for the classification stage. Support vector machine classification reached 92.3% accuracy with 95% sensitivity and 89.6% specificity using a 10-fold cross-validation. Most of the selected regions were well known degenerated brain regions including cerebellar and visual cortices, but at the same time, our whole-brain connectivity analysis revealed new regions not previously reported involving temporal and prefrontal cortices.

Conclusion

Our whole-brain connectivity approach provided information that seed-based analysis missed due to its region-specific searching method. The high classification accuracy suggests that using resting state functional connectivity may be a useful biomarker in SCA 7.

Electronic supplementary material

The online version of this article (doi:10.1186/2053-8871-1-2) contains supplementary material, which is available to authorized users.