Building connectomes using diffusion MRI: why, how and but

Why has diffusion MRI become a principal modality for mapping connectomes in vivo? How do different image acquisition parameters, fiber tracking algorithms and other methodological choices affect connectome estimation? What are the main factors that dictate the success and failure of connectome reconstruction? These are some of the key questions that we aim to address in this review. We provide an overview of the key methods that can be used to estimate the nodes and edges of macroscale connectomes, and we discuss open problems and inherent limitations. We argue that diffusion MRI-based connectome mapping methods are still in their infancy and caution against blind application of deep white matter tractography due to the challenges inherent to connectome reconstruction. We review a number of studies that provide evidence of useful microstructural and network properties that can be extracted in various independent and biologically-relevant contexts. Finally, we highlight some of the key deficiencies of current macroscale connectome mapping methodologies and motivate future developments

Using diffusion MRI to discriminate areas of cortical grey matter

Cortical area parcellation is a challenging problem that is often approached by combining structural imaging (e.g., quantitative T1, diffusion-based connectivity) with functional imaging (e.g., task activations, topological mapping, resting state correlations). Diffusion MRI (dMRI) has been widely adopted to analyse white matter microstructure, but scarcely used to distinguish grey matter regions because of the reduced anisotropy there. Nevertheless, differences in the texture of the cortical 'fabric' have long been mapped by histologists to distinguish cortical areas. Reliable area-specific contrast in the dMRI signal has previously been demonstrated in selected occipital and sensorimotor areas. We expand upon these findings by testing several diffusion-based feature sets in a series of classification tasks. Using Human Connectome Project (HCP) 3T datasets and a supervised learning approach, we demonstrate that diffusion MRI is sensitive to architectonic differences between a large number of different cortical areas defined in the HCP parcellation. By employing a surface-based cortical imaging pipeline, which defines diffusion features relative to local cortical surface orientation, we show that we can differentiate areas from their neighbours with higher accuracy than when using only fractional anisotropy or mean diffusivity. The results suggest that grey matter diffusion may provide a new, independent source of information for dividing up the cortex

Fractional order magnetic resonance fingerprinting in the human cerebral cortex

Mathematical models are becoming increasingly important in magnetic resonance imaging (MRI), as they provide a mechanistic approach for making a link between tissue microstructure and signals acquired using the medical imaging instrument. The Bloch equations, which describes spin and relaxation in a magnetic field, is a set of integer order differential equations with a solution exhibiting mono-exponential behaviour in time. Parameters of the model may be estimated using a non-linear solver, or by creating a dictionary of model parameters from which MRI signals are simulated and then matched with experiment. We have previously shown the potential efficacy of a magnetic resonance fingerprinting (MRF) approach, i.e. dictionary matching based on the classical Bloch equations, for parcellating the human cerebral cortex. However, this classical model is unable to describe in full the mm-scale MRI signal generated based on an heterogenous and complex tissue micro-environment. The time-fractional order Bloch equations has been shown to provide, as a function of time, a good fit of brain MRI signals. We replaced the integer order Bloch equations with the previously reported time-fractional counterpart within the MRF framework and performed experiments to parcellate human gray matter, which is cortical brain tissue with different cyto-architecture at different spatial locations. Our findings suggest that the time-fractional order parameters, {\alpha} and {\beta}, potentially associate with the effect of interareal architectonic variability, hypothetically leading to more accurate cortical parcellation

A Symmetry-Based Method to Infer Structural Brain Networks from Probabilistic Tractography Data

Recent progress in diffusion MRI and tractography algorithms as well as the launch of the Human Connectome Project (HCP) have provided brain research with an abundance of structural connectivity data. In this work, we describe and evaluate a method that can infer the structural brain network that interconnects a given set of Regions of Interest (ROIs) from probabilistic tractography data. The proposed method, referred to as Minimum Asymmetry Network Inference Algorithm (MANIA), does not determine the connectivity between two ROIs based on an arbitrary connectivity threshold. Instead, we exploit a basic limitation of the tractography process: the observed streamlines from a source to a target do not provide any information about the polarity of the underlying white matter, and so if there are some fibers connecting two voxels (or two ROIs) X and Y, tractography should be able in principle to follow this connection in both directions, from X to Y and from Y to X. We leverage this limitation to formulate the network inference process as an optimization problem that minimizes the (appropriately normalized) asymmetry of the observed network. We evaluate the proposed method using both the FiberCup dataset and based on a noise model that randomly corrupts the observed connectivity of synthetic networks. As a case-study, we apply MANIA on diffusion MRI data from 28 healthy subjects to infer the structural network between 18 corticolimbic ROIs that are associated with various neuropsychiatric conditions including depression, anxiety and addiction