Analysis of Dynamic Brain Imaging Data

Modern imaging techniques for probing brain function, including functional Magnetic Resonance Imaging, intrinsic and extrinsic contrast optical imaging, and magnetoencephalography, generate large data sets with complex content. In this paper we develop appropriate techniques of analysis and visualization of such imaging data, in order to separate the signal from the noise, as well as to characterize the signal. The techniques developed fall into the general category of multivariate time series analysis, and in particular we extensively use the multitaper framework of spectral analysis. We develop specific protocols for the analysis of fMRI, optical imaging and MEG data, and illustrate the techniques by applications to real data sets generated by these imaging modalities. In general, the analysis protocols involve two distinct stages: `noise' characterization and suppression, and `signal' characterization and visualization. An important general conclusion of our study is the utility of a frequency-based representation, with short, moving analysis windows to account for non-stationarity in the data. Of particular note are (a) the development of a decomposition technique (`space-frequency singular value decomposition') that is shown to be a useful means of characterizing the image data, and (b) the development of an algorithm, based on multitaper methods, for the removal of approximately periodic physiological artifacts arising from cardiac and respiratory sources.Comment: 40 pages; 26 figures with subparts including 3 figures as .gif files. Originally submitted to the neuro-sys archive which was never publicly announced (was 9804003

Brain connectivity analysis: a short survey

This short survey the reviews recent literature on brain connectivity studies. It encompasses all forms of static and dynamic connectivity whether anatomical, functional, or effective. The last decade has seen an ever increasing number of studies devoted to deduce functional or effective connectivity, mostly from functional neuroimaging experiments. Resting state conditions have become a dominant experimental paradigm, and a number of resting state networks, among them the prominent default mode network, have been identified. Graphical models represent a convenient vehicle to formalize experimental findings and to closely and quantitatively characterize the various networks identified. Underlying these abstract concepts are anatomical networks, the so-called connectome, which can be investigated by functional imaging techniques as well. Future studies have to bridge the gap between anatomical neuronal connections and related functional or effective connectivities

Spatiotemporal dynamics of low frequency fluctuations in bold fMRI

Traditional fMRI utilizes blood oxygenation level dependent (BOLD) contrast to map brain activity. BOLD signal is sensitive to the hemodynamic changes associated with brain activity, and gives an indirect measure of brain activity. Low frequency fluctuations (LFFs) have been observed in the BOLD signal even in the absence of any anesthetic agent, and the correlations between the fluctuations from different brain regions has been used to map functional connectivity in the brain. Most studies involving spontaneous fluctuations in the BOLD signal extract connectivity patterns that show relationships between brain areas that are maintained over the length of the scanning session. The research presented in this document investigates the spatiotemporal dynamics of the BOLD fluctuations to identify common spatiotemporal patterns within a scan. First, the presence of a visually detectable spatiotemporal propagation pattern is demonstrated by utilizing single-slice data with high spatial and temporal resolution. The pattern consists of lateral-medial propagation of BOLD signal, demonstrating the presence of time-varying features in spontaneous BOLD fluctuations. Further, a novel pattern finding algorithm is developed for detecting repeated spatiotemporal patterns in BOLD fMRI data. The algorithm is applied to high temporal resolution T2*-weighted multislice images obtained from rats and humans in the absence of any task or stimulation. In rats, the primary pattern consists of waves of high signal intensity, propagating in a lateral-medial direction across the cortex, replicating the results obtained using visual observation. In humans, the most common spatiotemporal pattern consisted of an alteration between activation of areas comprising the "default-mode" (e.g., posterior cingulate and anterior medial prefrontal cortices) and the "task-positive" (e.g., superior parietal and premotor cortices) networks. Signal propagation from focal starting points is also observed. The pattern finding algorithm is shown to be reasonably insensitive to the variation in user-defined parameters, and the results are consistent within and between subjects. This novel approach for probing the spontaneous network activity of the brain has implications for the interpretation of conventional functional connectivity studies, and may increase the amount of information that can be obtained from neuroimaging data.Ph.D.Committee Chair: Keilholz, Shella; Committee Member: Hu, Xiaoping; Committee Member: Jaeger, Dieter; Committee Member: Sathian, Krish; Committee Member: Schumacher, Eri

Tensor Analysis and Fusion of Multimodal Brain Images

Current high-throughput data acquisition technologies probe dynamical systems with different imaging modalities, generating massive data sets at different spatial and temporal resolutions posing challenging problems in multimodal data fusion. A case in point is the attempt to parse out the brain structures and networks that underpin human cognitive processes by analysis of different neuroimaging modalities (functional MRI, EEG, NIRS etc.). We emphasize that the multimodal, multi-scale nature of neuroimaging data is well reflected by a multi-way (tensor) structure where the underlying processes can be summarized by a relatively small number of components or "atoms". We introduce Markov-Penrose diagrams - an integration of Bayesian DAG and tensor network notation in order to analyze these models. These diagrams not only clarify matrix and tensor EEG and fMRI time/frequency analysis and inverse problems, but also help understand multimodal fusion via Multiway Partial Least Squares and Coupled Matrix-Tensor Factorization. We show here, for the first time, that Granger causal analysis of brain networks is a tensor regression problem, thus allowing the atomic decomposition of brain networks. Analysis of EEG and fMRI recordings shows the potential of the methods and suggests their use in other scientific domains.Comment: 23 pages, 15 figures, submitted to Proceedings of the IEE

Analysis of resting-state neurovascular coupling and locomotion-associated neural dynamics using wide-field optical mapping

Understanding the relationship between neural activity and cortical hemodynamics, or neurovascular coupling is the foundation to interpret neuroimaging signals such as functional magnetic resonance imaging (fMRI) which measure local changes in hemodynamics as a proxy for underlying neural activity. Even though the stereotypical stimulus-evoked hemodynamic response pattern with increased concentration of oxy- and total-hemoglobin and decrease in concentration of deoxy-hemoglobin has been well-recognized, the linearity of neurovascular coupling and its variances depending on brain state and tasks haven’t been thoroughly evaluated. To directly assess the cortical neurovascular coupling, simultaneous recordings of neural and hemodynamic activity were imaged by wide-field optical mapping (WFOM) over the bilateral dorsal surface of the mouse brain through a bilateral thinned-skull cranial window. Neural imaging is achieved through wide-field fluorescence imaging in animals expressing genetically encoded calcium sensor (Thy1-GCaMP). Hemodynamics are recorded via simultaneous imaging of multi-spectral reflectance. Significant hemodynamic crosstalk was found in the detected fluorescence signal and the physical model of the contamination, methods of correction as well as electrophysiological verification are presented. A linear model between neural and hemodynamic signals was used to fit spatiotemporal hemodynamics can be predicted by convolving local fluorescence changes with hemodynamic response functions derived through both deconvolution and gamma-variate fitting. Beyond confirming that the resting-state hemodynamics in the awake and anesthetized brain are coupled to underlying neural activity, the patterns of bilaterally symmetric spontaneous neural activity observed by WFOM emulate the functionally connected networks detected by fMRI. This result provides reassurance that resting-state functional connectivity has neural origins. With the access to cortical neural activity at mesoscopic level, we further explore the cortical neural representations preceding and during spontaneous locomotion