The impact of normalization and segmentation on resting state brain networks

Graph theory has recently received a lot of attention from the neuroscience community as a method to represent and characterize brain networks. Still, there is a lack of a gold standard for the methods that should be employed for the preprocessing of the data and the construction of the networks, as well as a lack of knowledge on how different methodologies can affect the metrics reported. We used graph theory analysis applied to resting-state functional Magnetic Resonance Imaging (rs-fMRI) to investigate the influence of different node-defining strategies and the effect of normalizing the functional acquisition on several commonly reported metrics used to characterize brain networks. The nodes of the network were defined using either the individual FreeSurfer segmentation of each subject or the FreeSurfer segmented MNI (Montreal National Institute) 152 template, using the Destrieux and sub-cortical atlas. The functional acquisition was either kept on the functional native space or normalized into MNI standard space. The comparisons were done at three levels: on the connections, on the edge properties and on the network properties levels. Our results reveal that different registration and brain parcellation strategies have a strong impact on all the levels of analysis, possibly favoring the use of individual segmentation strategies and conservative registration approaches. In conclusion, several technical aspects must be considered so that graph theoretical analysis of connectivity MRI data can provide a framework to understand brain pathologies.(undefined

Graph analysis of functional brain networks: practical issues in translational neuroscience

The brain can be regarded as a network: a connected system where nodes, or units, represent different specialized regions and links, or connections, represent communication pathways. From a functional perspective communication is coded by temporal dependence between the activities of different brain areas. In the last decade, the abstract representation of the brain as a graph has allowed to visualize functional brain networks and describe their non-trivial topological properties in a compact and objective way. Nowadays, the use of graph analysis in translational neuroscience has become essential to quantify brain dysfunctions in terms of aberrant reconfiguration of functional brain networks. Despite its evident impact, graph analysis of functional brain networks is not a simple toolbox that can be blindly applied to brain signals. On the one hand, it requires a know-how of all the methodological steps of the processing pipeline that manipulates the input brain signals and extract the functional network properties. On the other hand, a knowledge of the neural phenomenon under study is required to perform physiological-relevant analysis. The aim of this review is to provide practical indications to make sense of brain network analysis and contrast counterproductive attitudes

Learning and comparing functional connectomes across subjects

Functional connectomes capture brain interactions via synchronized fluctuations in the functional magnetic resonance imaging signal. If measured during rest, they map the intrinsic functional architecture of the brain. With task-driven experiments they represent integration mechanisms between specialized brain areas. Analyzing their variability across subjects and conditions can reveal markers of brain pathologies and mechanisms underlying cognition. Methods of estimating functional connectomes from the imaging signal have undergone rapid developments and the literature is full of diverse strategies for comparing them. This review aims to clarify links across functional-connectivity methods as well as to expose different steps to perform a group study of functional connectomes

Modeling brain dynamics in brain tumor patients using the virtual brain

Presurgical planning for brain tumor resection aims at delineating eloquent tissue in the vicinity of the lesion to spare during surgery. To this end, noninvasive neuroimaging techniques such as functional MRI and diffusion-weighted imaging fiber tracking are currently employed. However, taking into account this information is often still insufficient, as the complex nonlinear dynamics of the brain impede straightforward prediction of functional outcome after surgical intervention. Large-scale brain network modeling carries the potential to bridge this gap by integrating neuroimaging data with biophysically based models to predict collective brain dynamics. As a first step in this direction, an appropriate computational model has to be selected, after which suitable model parameter values have to be determined. To this end, we simulated large-scale brain dynamics in 25 human brain tumor patients and 11 human control participants using The Virtual Brain, an open-source neuroinformatics platform. Local and global model parameters of the Reduced Wong-Wang model were individually optimized and compared between brain tumor patients and control subjects. In addition, the relationship between model parameters and structural network topology and cognitive performance was assessed. Results showed (1) significantly improved prediction accuracy of individual functional connectivity when using individually optimized model parameters; (2) local model parameters that can differentiate between regions directly affected by a tumor, regions distant from a tumor, and regions in a healthy brain; and (3) interesting associations between individually optimized model parameters and structural network topology and cognitive performance

Estimating Time-Varying Effective Connectivity in High-Dimensional fMRI Data Using Regime-Switching Factor Models

Recent studies on analyzing dynamic brain connectivity rely on sliding-window analysis or time-varying coefficient models which are unable to capture both smooth and abrupt changes simultaneously. Emerging evidence suggests state-related changes in brain connectivity where dependence structure alternates between a finite number of latent states or regimes. Another challenge is inference of full-brain networks with large number of nodes. We employ a Markov-switching dynamic factor model in which the state-driven time-varying connectivity regimes of high-dimensional fMRI data are characterized by lower-dimensional common latent factors, following a regime-switching process. It enables a reliable, data-adaptive estimation of change-points of connectivity regimes and the massive dependencies associated with each regime. We consider the switching VAR to quantity the dynamic effective connectivity. We propose a three-step estimation procedure: (1) extracting the factors using principal component analysis (PCA) and (2) identifying dynamic connectivity states using the factor-based switching vector autoregressive (VAR) models in a state-space formulation using Kalman filter and expectation-maximization (EM) algorithm, and (3) constructing the high-dimensional connectivity metrics for each state based on subspace estimates. Simulation results show that our proposed estimator outperforms the K-means clustering of time-windowed coefficients, providing more accurate estimation of regime dynamics and connectivity metrics in high-dimensional settings. Applications to analyzing resting-state fMRI data identify dynamic changes in brain states during rest, and reveal distinct directed connectivity patterns and modular organization in resting-state networks across different states.Comment: 21 page

Comparison of longitudinal changes in resting state functional magnetic resonance imaging between alzheimer’s and healthy controls

Resting State Functional Magnetic Resonance Imaging (rs-fMRI) is a technique that is widely used for analyzing brain function using different approaches and methods. This study involves rs-fMRI analysis of Blood Oxygenation Level Dependent (BOLD) signals acquired from Alzheimer’s disease (AD) Patients and Healthy Controls (HC). Each subject in the study had both functional and anatomical images with at least one rs-fMRI scan with their Anatomical (T1) scans. Previous rs-fMRI studies have demonstrated that AD shows differences in Amplitude of Low Frequency (\u3c0.1 Hz) Fluctuations (ALFF), and Regional Homogeneity (ReHo) measures according to HCs. The aim of the study is to investigate individual and group level differences using ReHo and mALFF related measures in a longitudinal analysis. The hypothesis is that with the age and group (AD or HC) of the subject, it is possible to separate AD and HC subjects from each other using 3 different ROIs (DMN – MT – MV), These regions are known to show abnormalities in AD patients but clinical wise never been identified as neuroimaging biomarkers. This study tries to check these ROIs to see if there are significant differences between the AD patients and HCs using 3 different features

Retrieving the hemodynamic response function in resting state fMRI: methodology and application

In this paper we present a procedure to retrieve the hemodynamic response function (HRF) from resting state functional magnetic resonance imaging (fMRI) data. The fundamentals of the procedures are further validated by considering simultaneous electroencephalographic (EEG) recordings. The typical HRF shape at rest for a group of healthy subject is presented. Then we present the modifications to the shape of the HRF at rest following two physiological modulations: eyes open versus eyes closed and propofol-induced modulations of consciousness

Evaluation of atlas-based segmentation of hippocampi in healthy humans

Introduction and aim: Region of interest (ROI)-based functional magnetic resonance imaging (fMRI) data analysis relies on extracting signals from a specific area which is presumed to be involved in the brain activity being studied. The hippocampus is of interest in many functional connectivity studies for example in epilepsy as it plays an important role in epileptogenesis. In this context, ROI may be defined using different techniques. Our study aims at evaluating the spatial correspondence of hippocampal ROIs obtained using three brain atlases with hippocampal ROI obtained using an automatic segmentation algorithm dedicated to the hippocampus. Material and methods: High-resolution volumetric T1-weighted MR images of 18 healthy volunteers (five females) were acquired on a 3T scanner. Individual ROIs for both hippocampi of each subject were segmented from the MR images using an automatic hippocampus and amygdala segmentation software called SACHA providing the gold standard ROI for comparison with the atlas-derived results. For each subject, hippocampal ROIs were also obtained using three brain atlases: PickAtlas available as a commonly used software toolbox; automated anatomical labeling (AAL) atlas included as a subset of ROI into PickAtlas toolbox and a frequency-based brain atlas by Hammers et al. The levels of agreement between the SACHA results and those obtained using the atlases were assessed based on quantitative indices measuring volume differences and spatial overlap. The comparison was performed in standard Montreal Neurological Institute space, the registration being obtained with SPM5 (http://www.fil.ion.ucl.ac.uk/spm/). Results: The mean volumetric error across all subjects was 73% for hippocampal ROIs derived from AAL atlas; 20% in case of ROIs derived from the Hammers atlas and 107% for ROIs derived from PickAtlas. The mean false-positive and false-negative classification rates were 60% and 10% respectively for the AAL atlas; 16% and 32% for the Hammers atlas and 6% and 72% for the PickAtlas. Conclusion: Though atlas-based ROI definition may be convenient, the resulting ROIs may be poor representations of the hippocampus in some studies critical to under- or oversampling. Performance of the AAL atlas was inferior to that of the Hammers atlas. Hippocampal ROIs derived from PickAtlas are highly significantly smaller, and this results in the worst performance out of three atlases. It is advisable that the defined ROIs should be verified with knowledge of neuroanatomy before using it for further data analysis