Characterizing brain anatomical connections using diffusion weighted MRI and graph theory
Introduction
Random motion of water molecules inside the brain is influenced by the architectural properties of tissues. Water diffusion is known to be highly anisotropic in certain white matter regions, with preferential movement along the nervous fibers. A recent development of a non-invasive technique which quantifies water diffusion process, known as Diffusion Weighted Magnetic Resonance Imaging (DW-MRI), has allowed to obtain structural information about the intravoxel axon arrangement (Basser et al., 1994, LeBihan, 2003). Based on this information, fiber tractography arises as a crucial technique to attain a better in vivo anatomical characterization of the brain (Mori et al., 1999, Conturo et al., 1999, Tuch, 2002, Parker et al., 2002, Koch et al., 2002, Behrens et al., 2003). Also, quantification of the anatomical connectivity between different gray matter structures would be a significant contribution to the understanding of functional integration of the human brain (LeBihan et al., 2001, Koch et al., 2002, Ramnani et al., 2004, Sporns et al., 2005, Sotero et al., 2007).
Reconstruction of nervous fiber trajectories is an extensively treated topic. In the traditional Streamline Tractography (SLT) approach (Mori et al., 1999, Conturo et al., 1999, Basser et al., 2000), a continuous trajectory is traced tangential to the direction of the principal eigenvector of the diffusion tensor measured at each voxel using a discretization step smaller than the size of the voxel. This approach usually fails in voxels where fibers cross each other, merge, kiss or diverge, and it is very sensitive to the influence of MR signal noise (Basser and Pajevic, 2000, Lori et al., 2002). In those situations, traced path strays from the real trajectory of nervous fibers. To overcome these limitations, modified Streamline Tractography (mSLT) methods based on Diffusion Tensor Deflection (Weinstein et al., 1999, Lazar et al., 2003) and Probabilistic Monte-Carlo Method (Parker and Alexander, 2003) have been proposed. The former uses the entire diffusion tensor to deflect the propagation direction computed in the previous step. The latter uses the uncertainty of the estimated nervous fiber orientation to compute a large number of possible paths from the seed point; a quantity can be assigned to each path reflecting some connectivity relationship between seed and target points.
In recent years, many other mSLT methods have been proposed (Tuch, 2002, Tench et al., 2002, Behrens et al., 2003, Hagmann et al., 2003). Usually, they define the anatomical connection probability between seed and target voxels as the ratio between the number of shared paths and the number of generated paths.
In contrast to SLT and mSLT methods, Level Set-Based Fast Marching (FM) techniques (Parker et al., 2002, Staempfli et al., 2006) express the tractography in terms of a wave front that emanates from a source point and whose evolution is controlled by the diffusion data. FM methods have two advantages over the SLT and mSLT methods: 1, better performance in situations of branching and fiber crossing, and 2, direct estimation of the probability of white matter connectivity between two points (Jun Zhang et al., 2005).
In FM methods, front evolution speed and direction in a voxel depend on the measured diffusion tensor. Generally, all proposed FM algorithms have used only the principal eigenvector of the diffusion tensor, therefore these methods fail to reconstruct fiber pathways in those places where fibers cross, merge, kiss or diverge. For dealing with this limitation, recently Staempfli et al. (2006) proposed an advanced implementation of FM (aFM), combining the advantages of classical FM and the tensor deflection approach. The objective is to take into account the entire information contained in the diffusion tensor. As an intrinsic limitation, aFM needs an empirical threshold value to classify geometrically the diffusion tensor ellipsoid (i.e. prolate, oblate or spherical tensor) and therefore to set the corresponding speed function. Also, only four possible situations of voxel transitions are considered, which are those involving prolate and oblate cases. Thus, some combinations of more than two fibers crossing may be ignored.
On the other hand, although probability of connection between seed and target voxels has been previously used (Tuch, 2002, Parker et al., 2002, Koch et al., 2002, Parker et al., 2003, Behrens et al., 2003, Staempfli et al., 2006), the generalization of this concept to characterize anatomical connections between different brain gray matter structures is not straightforward. An initial approach (Iturria-Medina et al., 2005) was proposed to quantify the anatomical connection strength (ACS) between two gray matter structures using geometrical information from probabilistic fiber paths. ACS was considered proportional to the total area comprised by the fiber connector volume over the surfaces of the two connected structures. This was evaluated by counting the number of superficial voxels involved in the connection, where each voxel is weighted according to the validity of the paths that connect it with the second structure. A connectivity matrix estimated using the aforementioned approach was employed to couple several brain areas in a realistic neural mass model for the EEG generation, obtaining physiologically plausible results (Sotero et al., 2007).
In addition, recently Hagmann et al. (2006) proposed a technique based on graph theory to study the connectivity between small cortical areas. Nodes of a graph correspond to small cubic regions of interest (ROI) covering the brain gray matter. Fiber tractography is performed by initiating fibers over the whole brain and arc weight between any two ROIs is assigned according to the connection density between them. An unweighted version of this graph was constructed in order to analyze its small world and hierarchical properties.
In this work, our interest lies in the development of a DW-MRI-based methodology, capable of characterizing directly anatomical connections between brain gray matter structures, which can be defined according to cytoarchitectonic, histological or other sort of anatomical and functional information. In order to accomplish this, the graph framework is employed to introduce a new anatomical connectivity model. Firstly, each voxel of the cerebral volume is assumed to be a node of a non-directed weighted graph. In this case, the weight of an arc is considered to be proportional to the probability of the existence of a nervous fiber connecting its corresponding nodes. Probabilistic tissue segmentation and intravoxel white matter orientational distribution function (ODF) are combined to compute the arc weight. Secondly, an iterative algorithm is used to solve the most probable path problem between any two nodes in the graph, which we will indistinctly refer to as the most reliable connection route between these nodes. This approach allows to asses probabilistic anatomical connectivity maps between brain voxels. Finally, in order to assessing anatomical connectivity between K gray matter structures, the graph is partitioned in the corresponding K non-overlapped subsets and one subset containing the remaining nodes. This allowed for the definition of three different anatomical connectivity measures between any pair of gray matter structures: Anatomical Connection Strength (ACS), Anatomical Connection Density (ACD) and Anatomical Connection Probability (ACP).
Section snippets
Methods
This section will be devoted to present some basic elements of graph theory, as well as the principal steps of the proposed methodology: 1, definition of a Brain Graph, 2, introduction of an iterative fiber tracking algorithm and quantification of node–node connectivity and 3, definition of anatomical connectivity measures between gray matter areas. Details on experimental data to be used and its preprocessing will also appear.
Artificial data
The performance of the proposed methodology was explored using artificial DW-MRI data. Firstly, anatomical connections between different regions of interest (ROIs), defined for the straight crossing, curve crossing and Feuille trajectory were estimated (see Fig. 7). Each ROI consists of seven contiguous voxels representing a line perpendicular to the fiber tract being analyzed. Left column (Figs. 7a, c and e) shows the reconstructed connection routes (most probable paths) between ROI1 and the
Summary and discussions
In this work, we developed a diffusion imaging methodology capable to characterize anatomical connections between different brain gray matter structures. It consists of three basic steps: 1, definition of a Brain Graph model in which each voxel is considered as a node of a non-directed weighted graph; 2, the use of an iterative algorithm based on analysis of the voxels neighborhood to find the route of maximum probability between two nodes and the subsequent definition of the anatomical
Acknowledgments
Authors would like to thank Nelson Rubal Lorenzo, Agustín Lage-Castellanos and Erick C. Jones for their useful comments on the anatomical, statistical and algorithmic aspects of this work, respectively. Also, we thank to Roberto C. Sotero and Thomas Koenig, for the final revision of the manuscript.
References (51)
- et al.
Voxel-based morphometry-the methods
NeuroImage
(2000) - et al.
Estimation of the effective self-diffusion tensor from the NMR spin echo
J. Magn. Reson., Ser. B.
(1994) - et al.
DTI mapping of human brain connectivity: statistical fibre tracking and virtual dissection
NeuroImage
(2003) - et al.
Image distortion correction in fMRI: a quantitative evaluation
NeuroImage
(2002) - et al.
An investigation of functional and anatomical connectivity using magnetic resonance imaging
NeuroImage
(2002) - et al.
A probabilistic atlas of the human brain: theory and rationale for its development. The International Consortium for Brain Mapping (ICBM)
NeuroImage
(1995) - et al.
New approaches for exploring anatomical and functional connectivity in the human brain
Biol. Psychiatry
(2004) - et al.
Resolving fiber crossing using advanced fast marching tractography based on diffusion tensor imaging
NeuroImage
(2006) - Alemán-Gómez, Y., Melie-García, L., Valdes-Hernández, P., 2006. IBASPM: Toolbox for automatic parcellation of brain...
Multiple-fiber reconstruction algorithms for diffusion MRI
Ann. N. Y. Acad. Sci.
(2005)
Spatial transformations of diffusion tensor magnetic resonance images
IEEE Trans. Med. Imag.
Nonlinear spatial normalization using basis functions
Hum. Brain Mapp.
Statistical artifacts in diffusion tensor MRI (DT-MRI) caused by background noise
Magn. Reson. Med.
In vivo fiber tractography using DT-MRI data
Magn. Reson. Med.
Non-invasive mapping of connections between human thalamus and cortex using diffusion imaging
Nat. Neurosci.
Automatic 3D intersubject registration of MR volumetric data in standardized Talairach space
J. Comput. Assist. Tomogr.
Tracking neuronal fiber pathways in the living human brain
Proc. Natl. Acad. Sci. U. S. A.
A note on two problems in connexion with graphs
Numer. Math.
Large-scale simultaneous hypothesis testing: the choice of a null hypothesis
J. Am. Stat. Assoc.
Three-dimensional correlative imaging: applications in human brain mapping
Anatomía II. Folleto complementario
Graphs and Algorithms
Imaging the brain neuronal network with diffusion MRI: a way to understand its global architecture
Proc. ISMRM
Cited by (306)
Diagnostic of autism spectrum disorder based on structural brain MRI images using, grid search optimization, and convolutional neural networks
2023, Biomedical Signal Processing and ControlMultiple Classification of Brain MRI Autism Spectrum Disorder by Age and Gender Using Deep Learning
2024, Journal of Medical SystemsTopological alterations in white matter anatomical networks in cervical dystonia
2024, Research Square