A Scalable Approach to Independent Vector Analysis by Shared Subspace Separation for Multi-Subject fMRI Analysis

[Abstract]: Joint blind source separation (JBSS) has wide applications in modeling latent structures across multiple related datasets. However, JBSS is computationally prohibitive with high-dimensional data, limiting the number of datasets that can be included in a tractable analysis. Furthermore, JBSS may not be effective if the data’s true latent dimensionality is not adequately modeled, where severe overparameterization may lead to poor separation and time performance. In this paper, we propose a scalable JBSS method by modeling and separating the “shared” subspace from the data. The shared subspace is defined as the subset of latent sources that exists across all datasets, represented by groups of sources that collectively form a low-rank structure. Our method first provides the efficient initialization of the independent vector analysis (IVA) with a multivariate Gaussian source prior (IVA-G) specifically designed to estimate the shared sources. Estimated sources are then evaluated regarding whether they are shared, upon which further JBSS is applied separately to the shared and non-shared sources. This provides an effective means to reduce the dimensionality of the problem, improving analyses with larger numbers of datasets. We apply our method to resting-state fMRI datasets, demonstrating that our method can achieve an excellent estimation performance with significantly reduced computational costs.The computational hardware used is part of the UMBC High Performance Computing Facility (HPCF), supported by the US NSF through the MRI and SCREMS programs (grants CNS-0821258, CNS-1228778, OAC-1726023, CNS-1920079, DMS-0821311), with additional substantial support from the University of Maryland, Baltimore County (UMBC). This work was supported by the grants NIH R01 MH118695, NIH R01 MH123610, and NIH R01 AG073949. Xunta de Galicia was supported by a postdoctoral grant No. ED481B 2022/012 and the Fulbright Program, sponsored by the US Department of State.Xunta de Galicia; ED481B 2022/01

Intrinsic Inter-Subject Variability in Functional Neuroimaging: Verification Using Blind Source Separation Features

The holy grail of brain imaging is the identification of a biomarker, which can identify an abnormality that can be used to diagnose disease and track the effectiveness of treatment and disease progression. Typically approaches that search for biomarkers start by identifying mean differences between groups of patients and healthy controls. However, combining data from different subjects and groups to be able to make meaningful inferences is not trivial. The structure of the brain in each individual is unique in size and shape as well as in the relative location of anatomical landmarks (e.g. sulci and gyri). When looking for mean differences in functional images, this issue is exacerbated by the presence of variability in functional localization, i.e. variability in the location of functional regions in the brain. This is notably an important reason to focus on looking for inter-individual differences or variability. Inter-subject variability in neuroimaging experiments is often viewed as noise. The analyses are setup in a manner to ignore this variability assuming that a global spatial normalization brings the data into the same space. Nonetheless, functional activation patterns can be impacted by variability in multiple ways for e.g., there could be spatial variability of the maps or variability in the spectral composition of the timecourses or variability in the connectivity between the activation patterns identified. The overarching problem this thesis seeks to contribute to, is seeking improved measures to quantify biologically significant spatial, spectral and connectivity based variability and to identify associated cognitive or behavioral differences in the distribution of brain networks. We have successfully shown that different (spatial and spectral) measures of variability in blind source separated functional activation patterns underline previously unexplained characteristics that help in discerning schizophrenia patients from healthy controls. Additionally, we show that variance measures in dynamic connectivity between networks in healthy controls can justify relationship between connectivity patterns and executive functioning abilities

Any-way and Sparse Analyses for Multimodal Fusion and Imaging Genomics

This dissertation aims to develop new algorithms that leverage sparsity and mutual information across data modalities built upon the independent component analysis (ICA) framework to improve the performance of current ICA-based multimodal fusion approaches. These algorithms are further applied to both simulated data and real neuroimaging and genomic data to examine their performance. The identified neuroimaging and genomic patterns can help better delineate the pathology of mental disorders or brain development. To alleviate the signal-background separation difficulties in infomax-decomposed sources for genomic data, we propose a sparse infomax by enhancing a robust sparsity measure, the Hoyer index. Hoyer index is scale-invariant and well suited for ICA frameworks since the scale of decomposed sources is arbitrary. Simulation results demonstrate that sparse infomax increases the component detection accuracy for situations where the source signal-to-background (SBR) ratio is low, particularly for single nucleotide polymorphism (SNP) data. The proposed sparse infomax is further extended into two data modalities as a sparse parallel ICA for applications to imaging genomics in order to investigate the associations between brain imaging and genomics. Simulation results show that sparse parallel ICA outperforms parallel ICA with improved accuracy for structural magnetic resonance imaging (sMRI)-SNP association detection and component spatial map recovery, as well as with enhanced sparsity for sMRI and SNP components under noisy cases. Applying the proposed sparse parallel ICA to fuse the whole-brain sMRI and whole-genome SNP data of 24985 participants in the UK biobank, we identify three stable and replicable sMRI-SNP pairs. The identified sMRI components highlight frontal, parietal, and temporal regions and associate with multiple cognitive measures (with different association strengths in different age groups for the temporal component). Top SNPs in the identified SNP factor are enriched in inflammatory disease and inflammatory response pathways, which also regulate gene expression, isoform percentage, transcription expression, or methylation level in the frontal region, and the regulation effects are significantly enriched. Applying the proposed sparse parallel ICA to imaging genomics in attention-deficit/hyperactivity disorder (ADHD), we identify and replicate one SNP component related to gray matter volume (GMV) alterations in superior and middle frontal gyri underlying working memory deficit in adults and adolescents with ADHD. The association is more significant in ADHD families than controls and stronger in adults and older adolescents than younger ones. The identified SNP component highlights SNPs in long non-coding RNAs (lncRNAs) in chromosome 5 and in several protein-coding genes that are involved in ADHD, such as MEF2C, CADM2, and CADPS2. Top SNPs are enriched in human brain neuron cells and regulate gene expression, isoform percentage, transcription expression, or methylation level in the frontal region. Moreover, to increase the flexibility and robustness in mining multimodal data, we propose aNy-way ICA, which optimizes the entire correlation structure of linked components across any number of modalities via the Gaussian independent vector analysis and simultaneously optimizes independence via separate (parallel) ICAs. Simulation results demonstrate that aNy-way ICA recover sources and loadings, as well as the true covariance patterns with improved accuracy compared to existing multimodal fusion approaches, especially under noisy conditions. Applying the proposed aNy-way ICA to integrate structural MRI, fractal n-back, and emotion identification task functional MRIs collected in the Philadelphia Neurodevelopmental Cohort (PNC), we identify and replicate one linked GMV-threat-2-back component, and the threat and 2-back components are related to intelligence quotient (IQ) score in both discovery and replication samples. Lastly, we extend the proposed aNy-way ICA with a reference constraint to enable prior-guided multimodal fusion. Simulation results show that aNy-way ICA with reference recovers the designed linkages between reference and modalities, cross-modality correlations, as well as loading and component matrices with improved accuracy compared to multi-site canonical correlation analysis with reference (MCCAR)+joint ICA under noisy conditions. Applying aNy-way ICA with reference to supervise structural MRI, fractal n-back, and emotion identification task functional MRIs fusion in PNC with IQ as the reference, we identify and replicate one IQ-related GMV-threat-2-back component, and this component is significantly correlated across modalities in both discovery and replication samples.Ph.D

Assisted Dictionary Learning for fMRI Data Analysis

Extracting information from functional magnetic resonance (fMRI) images has been a major area of research for more than two decades. The goal of this work is to present a new method for the analysis of fMRI data sets, that is capable to incorporate a priori available information, via an efficient optimization framework. Tests on synthetic data sets demonstrate significant performance gains over existing methods of this kind.Comment: 5 pages, 2 figure

In vivo widefield calcium imaging of the mouse cortex for analysis of network connectivity in health and brain disease

The organization of brain areas in functionally connected networks, their dynamic changes, and perturbations in disease states are subject of extensive investigations. Research on functional networks in humans predominantly uses functional magnetic resonance imaging (fMRI). However, adopting fMRI and other functional imaging methods to mice, the most widely used model to study brain physiology and disease, poses major technical challenges and faces important limitations. Hence, there is great demand for alternative imaging modalities for network characterization. Here, we present a refined protocol for in vivo widefield calcium imaging of both cerebral hemispheres in mice expressing a calcium sensor in excitatory neurons. We implemented a stringent protocol for minimizing anesthesia and excluding movement artifacts which both imposed problems in previous approaches. We further adopted a method for unbiased identification of functional cortical areas using independent component analysis (ICA) on resting-state imaging data. Biological relevance of identified components was confirmed using stimulus-dependent cortical activation. To explore this novel approach in a model of focal brain injury, we induced photothrombotic lesions of the motor cortex, determined changes in inter- and intrahemispheric connectivity at multiple time points up to 56 days post-stroke and correlated them with behavioral deficits. We observed a severe loss in interhemispheric connectivity after stroke, which was partially restored in the chronic phase and associated with corresponding behavioral motor deficits. Taken together, we present an improved widefield calcium imaging tool accounting for anesthesia and movement artifacts, adopting an advanced analysis pipeline based on human fMRI algorithms and with superior sensitivity to recovery mechanisms in mouse models compared to behavioral tests. This tool will enable new studies on interhemispheric connectivity in murine models with comparability to human imaging studies for a wide spectrum of neuroscience applications in health and disease