Human brain mapping: a systematic comparison of parcellation methods for the human cerebral cortex

The macro-connectome elucidates the pathways through which brain regions are structurally connected or functionally coupled to perform a specific cognitive task. It embodies the notion of representing and understanding all connections within the brain as a network, while the subdivision of the brain into interacting functional units is inherent in its architecture. As a result, the definition of network nodes is one of the most critical steps in connectivity network analysis. Although brain atlases obtained from cytoarchitecture or anatomy have long been used for this task, connectivity-driven methods have arisen only recently, aiming to delineate more homogeneous and functionally coherent regions. This study provides a systematic comparison between anatomical, connectivity-driven and random parcellation methods proposed in the thriving field of brain parcellation. Using resting-state functional MRI data from the Human Connectome Project and a plethora of quantitative evaluation techniques investigated in the literature, we evaluate 10 subject-level and 24 groupwise parcellation methods at different resolutions. We assess the accuracy of parcellations from four different aspects: (1) reproducibility across different acquisitions and groups, (2) fidelity to the underlying connectivity data, (3) agreement with fMRI task activation, myelin maps, and cytoarchitectural areas, and (4) network analysis. This extensive evaluation of different parcellations generated at the subject and group level highlights the strengths and shortcomings of the various methods and aims to provide a guideline for the choice of parcellation technique and resolution according to the task at hand. The results obtained in this study suggest that there is no optimal method able to address all the challenges faced in this endeavour simultaneously

Disentangling causal webs in the brain using functional Magnetic Resonance Imaging: A review of current approaches

In the past two decades, functional Magnetic Resonance Imaging has been used to relate neuronal network activity to cognitive processing and behaviour. Recently this approach has been augmented by algorithms that allow us to infer causal links between component populations of neuronal networks. Multiple inference procedures have been proposed to approach this research question but so far, each method has limitations when it comes to establishing whole-brain connectivity patterns. In this work, we discuss eight ways to infer causality in fMRI research: Bayesian Nets, Dynamical Causal Modelling, Granger Causality, Likelihood Ratios, LiNGAM, Patel's Tau, Structural Equation Modelling, and Transfer Entropy. We finish with formulating some recommendations for the future directions in this area

Functional Brain Organization in Space and Time

The brain is a network functionally organized at many spatial and temporal scales. To understand how the brain processes information, controls behavior and dynamically adapts to an ever-changing environment, it is critical to have a comprehensive description of the constituent elements of this network and how relationships between these elements may change over time. Decades of lesion studies, anatomical tract-tracing, and electrophysiological recording have given insight into this functional organization. Recently, however, resting state functional magnetic resonance imaging (fMRI) has emerged as a powerful tool for whole-brain non-invasive measurement of spontaneous neural activity in humans, giving ready access to macroscopic scales of functional organization previously much more difficult to obtain. This thesis aims to harness the unique combination of spatial and temporal resolution provided by functional MRI to explore the spatial and temporal properties of the functional organization of the brain. First, we establish an approach for defining cortical areas using transitions in correlated patterns of spontaneous BOLD activity (Chapter 2). We then propose and apply measures of internal and external validity to evaluate the credibility of the areal parcellation generated by this technique (Chapter 3). In chapter 4, we extend the study of functional brain organization to a highly sampled individual. We describe the idiosyncratic areal and systems-level organization of the individual relative to a standard group-average description. Further, we develop a model describing the reliability of BOLD correlation estimates across days that accounts for relevant sources of variability. Finally, in Chapter 5, we examine whether BOLD correlations meaningfully vary over the course of single resting-state scans

Investigating microstructural variation in the human hippocampus using non-negative matrix factorization

In this work we use non-negative matrix factorization to identify patterns of microstructural variance in the human hippocampus. We utilize high-resolution structural and diffusion magnetic resonance imaging data from the Human Connectome Project to query hippocampus microstructure on a multivariate, voxelwise basis. Application of non-negative matrix factorization identifies spatial components (clusters of voxels sharing similar covariance patterns), as well as subject weightings (individual variance across hippocampus microstructure). By assessing the stability of spatial components as well as the accuracy of factorization, we identified 4 distinct microstructural components. Furthermore, we quantified the benefit of using multiple microstructural metrics by demonstrating that using three microstructural metrics (T1-weighted/T2-weighted signal, mean diffusivity and fractional anisotropy) produced more stable spatial components than when assessing metrics individually. Finally, we related individual subject weightings to demographic and behavioural measures using a partial least squares analysis. Through this approach we identified interpretable relationships between hippocampus microstructure and demographic and behavioural measures. Taken together, our work suggests non-negative matrix factorization as a spatially specific analytical approach for neuroimaging studies and advocates for the use of multiple metrics for data-driven component analyses

Support vector classification analysis of resting state functional connectivity fMRI

Since its discovery in 1995 resting state functional connectivity derived from functional MRI data has become a popular neuroimaging method for study psychiatric disorders. Current methods for analyzing resting state functional connectivity in disease involve thousands of univariate tests, and the specification of regions of interests to employ in the analysis. There are several drawbacks to these methods. First the mass univariate tests employed are insensitive to the information present in distributed networks of functional connectivity. Second, the null hypothesis testing employed to select functional connectivity dierences between groups does not evaluate the predictive power of identified functional connectivities. Third, the specification of regions of interests is confounded by experimentor bias in terms of which regions should be modeled and experimental error in terms of the size and location of these regions of interests. The objective of this dissertation is to improve the methods for functional connectivity analysis using multivariate predictive modeling, feature selection, and whole brain parcellation. A method of applying Support vector classification (SVC) to resting state functional connectivity data was developed in the context of a neuroimaging study of depression. The interpretability of the obtained classifier was optimized using feature selection techniques that incorporate reliability information. The problem of selecting regions of interests for whole brain functional connectivity analysis was addressed by clustering whole brain functional connectivity data to parcellate the brain into contiguous functionally homogenous regions. This newly developed famework was applied to derive a classifier capable of correctly seperating the functional connectivity patterns of patients with depression from those of healthy controls 90% of the time. The features most relevant to the obtain classifier match those previously identified in previous studies, but also include several regions not previously implicated in the functional networks underlying depression.Ph.D.Committee Chair: Hu, Xiaoping; Committee Co-Chair: Vachtsevanos, George; Committee Member: Butera, Robert; Committee Member: Gurbaxani, Brian; Committee Member: Mayberg, Helen; Committee Member: Yezzi, Anthon

Adaptive cortical parcellations for source reconstructed EEG/MEG connectomes.

There is growing interest in the rich temporal and spectral properties of the functional connectome of the brain that are provided by Electro- and Magnetoencephalography (EEG/MEG). However, the problem of leakage between brain sources that arises when reconstructing brain activity from EEG/MEG recordings outside the head makes it difficult to distinguish true connections from spurious connections, even when connections are based on measures that ignore zero-lag dependencies. In particular, standard anatomical parcellations for potential cortical sources tend to over- or under-sample the real spatial resolution of EEG/MEG. By using information from cross-talk functions (CTFs) that objectively describe leakage for a given sensor configuration and distributed source reconstruction method, we introduce methods for optimising the number of parcels while simultaneously minimising the leakage between them. More specifically, we compare two image segmentation algorithms: 1) a split-and-merge (SaM) algorithm based on standard anatomical parcellations and 2) a region growing (RG) algorithm based on all the brain vertices with no prior parcellation. Interestingly, when applied to minimum-norm reconstructions for EEG/MEG configurations from real data, both algorithms yielded approximately 70 parcels despite their different starting points, suggesting that this reflects the resolution limit of this particular sensor configuration and reconstruction method. Importantly, when compared against standard anatomical parcellations, resolution matrices of adaptive parcellations showed notably higher sensitivity and distinguishability of parcels. Furthermore, extensive simulations of realistic networks revealed significant improvements in network reconstruction accuracies, particularly in reducing false leakage-induced connections. Adaptive parcellations therefore allow a more accurate reconstruction of functional EEG/MEG connectomes

Performing group-level functional image analyses based on homologous functional regions mapped in individuals

Functional MRI (fMRI) studies have traditionally relied on intersubject normalization based on global brain morphology, which cannot establish proper functional correspondence between subjects due to substantial intersubject variability in functional organization. Here, we reliably identified a set of discrete, homologous functional regions in individuals to improve intersubject alignment of fMRI data. These functional regions demonstrated marked intersubject variability in size, position, and connectivity. We found that previously reported intersubject variability in functional connectivity maps could be partially explained by variability in size and position of the functional regions. Importantly, individual differences in network topography are associated with individual differences in task-evoked activations, suggesting that these individually specified regions may serve as the localizer to improve the alignment of task-fMRI data. We demonstrated that aligning task-fMRI data using the regions derived from resting state fMRI may lead to increased statistical power of task-fMRI analyses. In addition, resting state functional connectivity among these homologous regions is able to capture the idiosyncrasies of subjects and better predict fluid intelligence (gF) than connectivity measures derived from group-level brain atlases. Critically, we showed that not only the connectivity but also the size and position of functional regions are related to human behavior. Collectively, these findings suggest that identifying homologous functional regions across individuals can benefit a wide range of studies in the investigation of connectivity, task activation, and brain-behavior associations. Author summary No two individuals are alike. The size, shape, position, and connectivity patterns of brain functional regions can vary drastically between individuals. While interindividual differences in functional organization are well recognized, to date, standard procedures for functional neuroimaging research still rely on aligning different subjects' data to a nominal average brain based on global brain morphology. We developed an approach to reliably identify homologous functional regions in each individual and demonstrated that aligning data based on these homologous functional regions can significantly improve the study of resting state functional connectivity, task-fMRI activations, and brain-behavior associations. Moreover, we showed that individual differences in size, position, and connectivity of brain functional regions are dissociable, and each can provide nonredundant information in explaining human behavior