Clusterwise Independent Component Analysis (C-ICA):An R package for clustering subjects based on ICA patterns underlying three-way (brain) data

In many areas of science, like neuroscience, genomics and text mining, several important and challenging research questions imply the study of (subject) heterogeneity present in three-way data. In clinical neuroscience, for example, disclosing differences or heterogeneity between subjects in resting state networks (RSNs) underlying multi-subject fMRI data (i.e., time by voxel by subject three-way data) may advance the subtyping of psychiatric and mental diseases. Recently, the Clusterwise Independent Component Analysis (C-ICA) method was proposed that enables the disclosure of heterogeneity between subjects in RSNs that is present in multi-subject rs-fMRI data [1]. Up to now, however, no publicly available software exists that allows to fit C-ICA to empirical data at hand. The goal of this paper, therefore, is to present the CICA R package, which contains the necessary functions to estimate the C-ICA parameters and to interpret and visualize the analysis output. Further, the package also includes functions to select suitable initial values for the C-ICA model parameters and to determine the optimal number of clusters and components for a given empirical data set (i.e., model selection). The use of the main functions of the package is discussed and demonstrated with simulated data. Herewith, the necessary analytical choices that have to be made by the user (e.g., starting values) are explained and showed step by step. The rich functionality of the package is further illustrated by applying C-ICA to empirical rs-fMRI data from a group of Alzheimer patients and elderly control subjects and to multi-country stock market data. Finally, extensions regarding the C-ICA algorithm and procedures for model selection that could be implemented in future releases of the package are discussed

Analyzing hierarchical multi-view MRI data with StaPLR: An application to Alzheimer's disease classification

Multi-view data refers to a setting where features are divided into feature sets, for example because they correspond to different sources. Stacked penalized logistic regression (StaPLR) is a recently introduced method that can be used for classification and automatically selecting the views that are most important for prediction. We introduce an extension of this method to a setting where the data has a hierarchical multi-view structure. We also introduce a new view importance measure for StaPLR, which allows us to compare the importance of views at any level of the hierarchy. We apply our extended StaPLR algorithm to Alzheimer's disease classification where different MRI measures have been calculated from three scan types: structural MRI, diffusion-weighted MRI, and resting-state fMRI. StaPLR can identify which scan types and which derived MRI measures are most important for classification, and it outperforms elastic net regression in classification performance.Comment: 36 pages, 9 figures. Accepted manuscrip

Gray Matter Covariance Networks as Classifiers and Predictors of Cognitive Function in Alzheimer's Disease

The study of shared variation in gray matter morphology may define neurodegenerative diseases beyond what can be detected from the isolated assessment of regional brain volumes. We, therefore, aimed to (1) identify SCNs (structural covariance networks) that discriminate between Alzheimer's disease (AD) patients and healthy controls (HC), (2) investigate their diagnostic accuracy in comparison and above established markers, and (3) determine if they are associated with cognitive abilities. We applied a random forest algorithm to identify discriminating networks from a set of 20 SCNs. The algorithm was trained on a main sample of 104 AD patients and 104 age-matched HC and was then validated in an independent sample of 28 AD patients and 28 controls from another center. Only two of the 20 SCNs contributed significantly to the discrimination between AD and controls. These were a temporal and a secondary somatosensory SCN. Their diagnostic accuracy was 74% in the original cohort and 80% in the independent samples. The diagnostic accuracy of SCNs was comparable with that of conventional volumetric MRI markers including whole brain volume and hippocampal volume. SCN did not significantly increase diagnostic accuracy beyond that of conventional MRI markers. We found the temporal SCN to be associated with verbal memory at baseline. No other associations with cognitive functions were seen. SCNs failed to predict the course of cognitive decline over an average of 18 months. We conclude that SCNs have diagnostic potential, but the diagnostic information gain beyond conventional MRI markers is limited

Combining anatomical, diffusion, and resting state functional magnetic resonance imaging for individual classification of mild and moderate Alzheimer's disease

AbstractMagnetic resonance imaging (MRI) is sensitive to structural and functional changes in the brain caused by Alzheimer's disease (AD), and can therefore be used to help in diagnosing the disease. Improving classification of AD patients based on MRI scans might help to identify AD earlier in the disease's progress, which may be key in developing treatments for AD. In this study we used an elastic net classifier based on several measures derived from the MRI scans of mild to moderate AD patients (N=77) from the prospective registry on dementia study and controls (N=173) from the Austrian Stroke Prevention Family Study. We based our classification on measures from anatomical MRI, diffusion weighted MRI and resting state functional MRI. Our unimodal classification performance ranged from an area under the curve (AUC) of 0.760 (full correlations between functional networks) to 0.909 (grey matter density). When combining measures from multiple modalities in a stepwise manner, the classification performance improved to an AUC of 0.952. This optimal combination consisted of grey matter density, white matter density, fractional anisotropy, mean diffusivity, and sparse partial correlations between functional networks. Classification performance for mild AD as well as moderate AD also improved when using this multimodal combination. We conclude that different MRI modalities provide complementary information for classifying AD. Moreover, combining multiple modalities can substantially improve classification performance over unimodal classification

Alterations and test-retest reliability of functional connectivity network measures in cerebral small vessel disease

While structural network analysis consolidated the hypothesis of cerebral small vessel disease (SVD) being a disconnection syndrome, little is known about functional changes on the level of brain networks. In patients with genetically defined SVD (CADASIL,n= 41) and sporadic SVD (n= 46), we independently tested the hypothesis that functional networks change with SVD burden and mediate the effect of disease burden on cognitive performance, in particular slowing of processing speed. We further determined test-retest reliability of functional network measures in sporadic SVD patients participating in a high-frequency (monthly) serial imaging study (RUN DMC-InTENse, median: 8 MRIs per participant). Functional networks for the whole brain and major subsystems (i.e., default mode network, DMN;fronto-parietal task control network, FPCN;visual network, VN;hand somatosensory-motor network, HSMN) were constructed based on resting-state multi-band functional MRI. In CADASIL, global efficiency (a graph metric capturing network integration) of the DMN was lower in patients with high disease burden (standardized beta = -.44;p[corrected] = .035) and mediated the negative effect of disease burden on processing speed (indirect path: std. beta = -.20,p= .047;direct path: std. beta = -.19,p= .25;total effect: std. beta = -.39,p= .02). The corresponding analyses in sporadic SVD showed no effect. Intraclass correlations in the high-frequency serial MRI dataset of the sporadic SVD patients revealed poor test-retest reliability and analysis of individual variability suggested an influence of age, but not disease burden, on global efficiency. In conclusion, our results suggest that changes in functional connectivity networks mediate the effect of SVD-related brain damage on cognitive deficits. However, limited reliability of functional network measures, possibly due to age-related comorbidities, impedes the analysis in elderly SVD patients

Genetic variants associated with subjective well-being, depressive symptoms, and neuroticism identified through genome-wide analyses

Very few genetic variants have been associated with depression and neuroticism, likely because of limitations on sample size in previous studies. Subjective well-being, a phenotype that is genetically correlated with both of these traits, has not yet been studied with genome-wide data. We conducted genome-wide association studies of three phenotypes: subjective well-being (n = 298,420), depressive symptoms (n = 161,460), and neuroticism (n = 170,911). We identify 3 variants associated with subjective well-being, 2 variants associated with depressive symptoms, and 11 variants associated with neuroticism, including 2 inversion polymorphisms. The two loci associated with depressive symptoms replicate in an independent depression sample. Joint analyses that exploit the high genetic correlations between the phenotypes (|ρ^| ≈ 0.8) strengthen the overall credibility of the findings and allow us to identify additional variants. Across our phenotypes, loci regulating expression in central nervous system and adrenal or pancreas tissues are strongly enriched for association.</p

Author Correction:Study of 300,486 individuals identifies 148 independent genetic loci influencing general cognitive function

Christina M. Lill, who contributed to analysis of data, was inadvertently omitted from the author list in the originally published version of this article. This has now been corrected in both the PDF and HTML versions of the article

Clusterwise Independent Component Analysis (C-ICA):Using fMRI resting state networks to cluster subjects and find neurofunctional subtypes

BACKGROUND: FMRI resting state networks (RSNs) are used to characterize brain disorders. They also show extensive heterogeneity across patients. Identifying systematic differences between RSNs in patients, i.e. discovering neurofunctional subtypes, may further increase our understanding of disease heterogeneity. Currently, no methodology is available to estimate neurofunctional subtypes and their associated RSNs simultaneously.NEW METHOD: We present an unsupervised learning method for fMRI data, called Clusterwise Independent Component Analysis (C-ICA). This enables the clustering of patients into neurofunctional subtypes based on differences in shared ICA-derived RSNs. The parameters are estimated simultaneously, which leads to an improved estimation of subtypes and their associated RSNs.RESULTS: In five simulation studies, the C-ICA model is successfully validated using both artificially and realistically simulated data (N = 30-40). The successful performance of the C-ICA model is also illustrated on an empirical data set consisting of Alzheimer's disease patients and elderly control subjects (N = 250). C-ICA is able to uncover a meaningful clustering that partially matches (balanced accuracy = .72) the diagnostic labels and identifies differences in RSNs between the Alzheimer and control cluster.COMPARISON WITH OTHER METHODS: Both in the simulation study and the empirical application, C-ICA yields better results compared to competing clustering methods (i.e., a two step clustering procedure based on single subject ICA's and a Group ICA plus dual regression variant thereof) that do not simultaneously estimate a clustering and associated RSNs. Indeed, the overall mean adjusted Rand Index, a measure for cluster recovery, equals 0.65 for C-ICA and ranges from 0.27 to 0.46 for competing methods.CONCLUSIONS: The successful performance of C-ICA indicates that it is a promising method to extract neurofunctional subtypes from multi-subject resting state-fMRI data. This method can be applied on fMRI scans of patient groups to study (neurofunctional) subtypes, which may eventually further increase understanding of disease heterogeneity.</p

FMRI to probe sex-related differences in brain function with multitasking.

Although established as a general notion in society, there is no solid scientific foundation for the existence of sex-differences in multitasking. Reaction time and accuracy in dual task conditions have an inverse relationship relative to single task, independently from sex. While a more disseminated network, parallel to decreasing accuracy and reaction time has been demonstrated in dual task fMRI studies, little is known so far whether there exist respective sex-related differences in activation.We subjected 20 women (mean age = 25.45; SD = 5.23) and 20 men (mean age = 27.55; SD = 4.00) to a combined verbal and spatial fMRI paradigm at 3.0T to assess sex-related skills, based on the assumption that generally women better perform in verbal tasks while men do better in spatial tasks. We also obtained behavioral tests for verbal and spatial intelligence, attention, executive functions, and working memory.No differences between women and men were observed in behavioral measures of dual-tasking or cognitive performance. Generally, brain activation increased with higher task load, mainly in the bilateral inferior and prefrontal gyri, the anterior cingulum, thalamus, putamen and occipital areas. Comparing sexes, women showed increased activation in the inferior frontal gyrus in the verbal dual-task while men demonstrated increased activation in the precuneus and adjacent visual areas in the spatial task.Against the background of equal cognitive and behavioral dual-task performance in both sexes, we provide first evidence for sex-related activation differences in functional networks for verbal and spatial dual-tasking