Andrew M. Michael

A Baseline for the Multivariate Comparison of Resting-State Networks

As the size of functional and structural MRI datasets expands, it becomes increasingly important to establish a baseline from which diagnostic relevance may be determined, a processing strategy that efficiently prepares data for analysis, and a statistical approach that identifies important effects in a manner that is both robust and reproducible. In this paper, we introduce a multivariate analytic approach that optimizes sensitivity and reduces unnecessary testing. We demonstrate the utility of this mega-analytic approach by identifying the effects of age and gender on the resting-state networks (RSNs) of 603 healthy adolescents and adults (mean age: 23.4 years, range: 12–71 years). Data were collected on the same scanner, preprocessed using an automated analysis pipeline based in SPM, and studied using group independent component analysis. RSNs were identified and evaluated in terms of three primary outcome measures: time course spectral power, spatial map intensity, and functional network connectivity. Results revealed robust effects of age on all three outcome measures, largely indicating decreases in network coherence and connectivity with increasing age. Gender effects were of smaller magnitude but suggested stronger intra-network connectivity in females and more inter-network connectivity in males, particularly with regard to sensorimotor networks. These findings, along with the analysis approach and statistical framework described here, provide a useful baseline for future investigations of brain networks in health and disease.

Investigation of relationships between fMRI brain networks in the spectral domain using ICA and Granger causality reveals distinct differences between schizophrenia patients and healthy controls

Functional network connectivity (FNC) is an approach that examines the relationships between brain networks (as opposed to functional connectivity (FC) that focuses upon the relationships between single voxels). FNC may help explain the complex relationships between distributed cerebral sites in the brain and possibly provide new understanding of neurological and psychiatric disorders such as schizophrenia. In this paper, we use independent component analysis (ICA) to extract the time courses of spatially independent components and then use these in Granger causality test (GCT) to investigate causal relationships between brain activation networks. We present results using both simulations and fMRI data of 155 subjects obtained during two different tasks. Unlike previous research, causal relationships are presented over different portions of the frequency spectrum in order to differentiate high and low-frequency effects and not merged in a scalar. The results obtained using Sternberg item recognition paradigm (SIRP) and auditory oddball (AOD) tasks showed FNC differentiations between schizophrenia and control groups, and explained how the two groups differed during these tasks. During the SIRP task, secondary visual and cerebellum activation networks served as hubs and included most complex relationships between the activated regions. Secondary visual and temporal lobe activations replaced these components during the AOD task.

Does function follow form?: Methods to fuse structural and functional brain images show decreased linkage in schizophrenia

When both structural magnetic resonance imaging (sMRI) and functional MRI (fMRI) data are collected they are typically analyzed separately and the joint information is not examined. Techniques that examine joint information can help to find hidden traits in complex disorders such as schizophrenia. The brain is vastly interconnected, and local brain morphology may influence functional activity at distant regions. In this paper we introduce three methods to identify inter-correlations among sMRI and fMRI voxels within the whole brain. We apply these methods to examine sMRI gray matter data and fMRI data derived from an auditory sensorimotor task from a large study of schizophrenia. In Method 1 the sMRI-fMRI cross-correlation matrix is reduced to a histogram and results show that healthy controls (HC) have stronger correlations than do patients with schizophrenia (SZ). In Method 2 the spatial information of sMRI-fMRI correlations is retained. Structural regions in the cerebellum and frontal regions show more positive and more negative correlations, respectively, with functional regions in HC than in SZ. In Method 3 significant sMRI-fMRI inter-regional links are detected, with regions in the cerebellum showing more significant positive correlations with functional regions in HC relative to SZ. Results from all three methods indicate that the linkage between gray matter and functional activation is stronger in HC than SZ. The methods introduced can be easily extended to comprehensively correlate large data sets.

Multifaceted genomic risk for brain function in schizophrenia.

Recently, deriving candidate endophenotypes from brain imaging data has become a valuable approach to study genetic influences on schizophrenia (SZ), whose pathophysiology remains unclear. In this work we utilized a multivariate approach, parallel independent component analysis, to identify genomic risk components associated with brain function abnormalities in SZ. 5157 candidate single nucleotide polymorphisms (SNPs) were derived from genome-wide array based on their possible connections with SZ and further investigated for their associations with brain activations captured with functional magnetic resonance imaging (fMRI) during a sensorimotor task. Using data from 92 SZ patients and 116 healthy controls, we detected a significant correlation (r=0.29; p=2.41 × 10(-5)) between one fMRI component and one SNP component, both of which significantly differentiated patients from controls. The fMRI component mainly consisted of precentral and postcentral gyri, the major activated regions in the motor task. On average, higher activation in these regions was observed in participants with higher loadings of the linked SNP component, predominantly contributed to by 253 SNPs. 138 identified SNPs were from known coding regions of 100 unique genes. 31 identified SNPs did not differ between groups, but moderately correlated with some other group-discriminating SNPs, indicating interactions among alleles contributing toward elevated SZ susceptibility. The genes associated with the identified SNPs participated in four neurotransmitter pathways: GABA receptor signaling, dopamine receptor signaling, neuregulin signaling and glutamate receptor signaling. In summary, our work provides further evidence for the complexity of genomic risk to the functional brain abnormality in SZ and suggests a pathological role of interactions between SNPs, genes and multiple neurotransmitter pathways.

Preserving subject variability in group fMRI analysis: performance evaluation of GICA vs. IVA

Independent component analysis (ICA) is a widely applied technique to derive functionally connected brain networks from fMRI data. Group ICA (GICA) and Independent Vector Analysis (IVA) are extensions of ICA that enable users to perform group fMRI analyses; however a full comparison of the performance limits of GICA and IVA has not been investigated. Recent interest in resting state fMRI data with potentially higher degree of subject variability makes the evaluation of the above techniques important. In this paper we compare component estimation accuracies of GICA and an improved version of IVA using simulated fMRI datasets. We systematically change the degree of inter-subject spatial variability of components and evaluate estimation accuracy over all spatial maps (SMs) and time courses (TCs) of the decomposition. Our results indicate the following: (1) at low levels of SM variability or when just one SM is varied, both GICA and IVA perform well, (2) at higher levels of SM variability or when more than one SMs are varied, IVA continues to perform well but GICA yields SM estimates that are composites of other SMs with errors in TCs, (3) both GICA and IVA remove spatial correlations of overlapping SMs and introduce artificial correlations in their TCs, (4) if number of SMs is over estimated, IVA continues to perform well but GICA introduces artifacts in the varying and extra SMs with artificial correlations in the TCs of extra components, and (5) in the absence or presence of SMs unique to one subject, GICA produces errors in TCs and IVA estimates are accurate. In summary, our simulation experiments (both simplistic and realistic) and our holistic analyses approach indicate that IVA produces results that are closer to ground truth and thereby better preserves subject variability. The improved version of IVA is now packaged into the GIFT toolbox (http://mialab.mrn.org/software/gift).

Heritability of Multivariate Gray Matter Measures in Schizophrenia

Structural brain measures are employed as endophenotypes in the search for schizophrenia susceptibility genes. We analyzed two independent structural imaging datasets with voxel-based morphometry and with source-based morphometry, a multivariate, independent components analysis, to determine the stability and heritability of regional gray matter concentration abnormalities in schizophrenia. The samples comprised 209 and 102 patients with schizophrenia and 208 and 96 healthy volunteers, respectively. The second sample additionally included non-ill siblings of participants with and without schizophrenia. A standard voxel-based analysis showed reproducible regional gray matter deficits in the affected participants compared with unrelated, unaffected controls in both datasets: patients showed significant gray matter concentration deficits in cortical frontal, temporal, and insular lobes. Source-based morphometry (SBM) was applied to the gray matter images of the entire sample to determine the effects of diagnosis on networks of covarying structures. The SBM analysis extracted 24 significant sets of covarying regions (components). Four of these components showed significantly lower gray matter concentrations in patients (p < .05). We determined the familiality of the observed SBM components based on 66 sibling pairs (25 discordant for schizophrenia). Two components, one including the medial frontal, insular, inferior frontal, and temporal lobes, and the other including the posterior occipital lobe, showed significant familiality (p < .05). We conclude that structural brain deficits in schizophrenia are replicable, and that SBM can extract unique familial and likely heritable components. SBM provides a useful data reduction technique that can provide measures that may serve as endophenotypes for schizophrenia.

Wavelet-based fMRI analysis: 3-D denoising, signal separation, and validation metrics

We present a novel integrated wavelet-domain based framework (w-ICA) for 3-D denoising functional magnetic resonance imaging (fMRI) data followed by source separation analysis using independent component analysis (ICA) in the wavelet domain. We propose the idea of a 3-D wavelet-based multi-directional denoising scheme where each volume in a 4-D fMRI data set is sub-sampled using the axial, sagittal and coronal geometries to obtain three different slice-by-slice representations of the same data. The filtered intensity value of an arbitrary voxel is computed as an expected value of the denoised wavelet coefficients corresponding to the three viewing geometries for each sub-band. This results in a robust set of denoised wavelet coefficients for each voxel. Given the de-correlated nature of these denoised wavelet coefficients, it is possible to obtain more accurate source estimates using ICA in the wavelet domain. The contributions of this work can be realized as two modules: First, in the analysis module we combine a new 3-D wavelet denoising approach with signal separation properties of ICA in the wavelet domain. This step helps obtain an activation component that corresponds closely to the true underlying signal, which is maximally independent with respect to other components. Second, we propose and describe two novel shape metrics for post-ICA comparisons between activation regions obtained through different frameworks. We verified our method using simulated as well as real fMRI data and compared our results against the conventional scheme (Gaussian smoothing+spatial ICA: s-ICA). The results show significant improvements based on two important features: (1) preservation of shape of the activation region (shape metrics) and (2) receiver operating characteristic curves. It was observed that the proposed framework was able to preserve the actual activation shape in a consistent manner even for very high noise levels in addition to significant reduction in false positive voxels.

ICA-fNORM: Spatial Normalization of fMRI Data Using Intrinsic Group-ICA Networks

A common pre-processing challenge associated with group level fMRI analysis is spatial registration of multiple subjects to a standard space. Spatial normalization, using a reference image such as the Montreal Neurological Institute brain template, is the most common technique currently in use to achieve spatial congruence across multiple subjects. This method corrects for global shape differences preserving regional asymmetries, but does not account for functional differences. We propose a novel approach to co-register task-based fMRI data using resting state group-ICA networks. We posit that these intrinsic networks (INs) can provide to the spatial normalization process with important information about how each individual’s brain is organized functionally. The algorithm is initiated by the extraction of single subject representations of INs using group level independent component analysis (ICA) on resting state fMRI data. In this proof-of-concept work two of the robust, commonly identified, networks are chosen as functional templates. As an estimation step, the relevant INs are utilized to derive a set of normalization parameters for each subject. Finally, the normalization parameters are applied individually to a different set of fMRI data acquired while the subjects performed an auditory oddball task. These normalization parameters, although derived using rest data, generalize successfully to data obtained with a cognitive paradigm for each subject. The improvement in results is verified using two widely applied fMRI analysis methods: the general linear model and ICA. Resulting activation patterns from each analysis method show significant improvements in terms of detection sensitivity and statistical significance at the group level. The results presented in this article provide initial evidence to show that common functional domains from the resting state brain may be used to improve the group statistics of task-fMRI data.

An ICA with reference approach in identification of genetic variation and associated brain networks

To address the statistical challenges associated with genome-wide association studies, we present an independent component analysis (ICA) with reference approach to target a specific genetic variation and associated brain networks. First, a small set of single nucleotide polymorphisms (SNPs) are empirically chosen to reflect a feature of interest and these SNPs are used as a reference when applying ICA to a full genomic SNP array. After extracting the genetic component maximally representing the characteristics of the reference, we test its association with brain networks in functional magnetic resonance imaging (fMRI) data. The method was evaluated on both real and simulated datasets. Simulation demonstrates that ICA with reference can extract a specific genetic factor, even when the variance accounted for by such a factor is so small that a regular ICA fails. Our real data application from 48 schizophrenia patients (SZs) and 40 healthy controls (HCs) include 300K SNPs and fMRI images in an auditory oddball task. Using SNPs with allelic frequency difference in two groups as a reference, we extracted a genetic component that maximally differentiates patients from controls (p < 4 × 10−17), and discovered a brain functional network that was significantly associated with this genetic component (p < 1 × 10−4). The regions in the functional network mainly locate in the thalamus, anterior and posterior cingulate gyri. The contributing SNPs in the genetic factor mainly fall into two clusters centered at chromosome 7q21 and chromosome 5q35. The findings from the schizophrenia application are in concordance with previous knowledge about brain regions and gene function. All together, the results suggest that the ICA with reference can be particularly useful to explore the whole genome to find a specific factor of interest and further study its effect on brain.

A Data-Driven Investigation of Gray Matter–Function Correlations in Schizophrenia during a Working Memory Task

The brain is a vastly interconnected organ and methods are needed to investigate its long range structure(S)–function(F) associations to better understand disorders such as schizophrenia that are hypothesized to be due to distributed disconnected brain regions. In previous work we introduced a methodology to reduce the whole brain S–F correlations to a histogram and here we reduce the correlations to brain clusters. The application of our approach to sMRI [gray matter (GM) concentration maps] and functional magnetic resonance imaging data (general linear model activation maps during Encode and Probe epochs of a working memory task) from patients with schizophrenia (SZ, n = 100) and healthy controls (HC, n = 100) presented the following results. In HC the whole brain correlation histograms for GM–Encode and GM–Probe overlap for Low and Medium loads and at High the histograms separate, but in SZ the histograms do not overlap for any of the load levels and Medium load shows the maximum difference. We computed GM–F differential correlation clusters using activation for Probe Medium, and they included regions in the left and right superior temporal gyri, anterior cingulate, cuneus, middle temporal gyrus, and the cerebellum. Inter-cluster GM–Probe correlations for Medium load were positive in HC but negative in SZ. Within group inter-cluster GM–Encode and GM–Probe correlation comparisons show no differences in HC but in SZ differences are evident in the same clusters where HC vs. SZ differences occurred for Probe Medium, indicating that the S–F integrity during Probe is aberrant in SZ. Through a data-driven whole brain analysis approach we find novel brain clusters and show how the S–F differential correlation changes during Probe and Encode at three memory load levels. Structural and functional anomalies have been extensively reported in schizophrenia and here we provide evidences to suggest that evaluating S–F associations can provide important additional information.