|Home | About | Journals | Submit | Contact Us | Français|
In this study, we propose a framework to map functional MRI (fMRI) activation signals using diffusion tensor imaging (DTI)-tractography. This framework, which we termed functional by structural hierarchical (FSH)mapping, takes into account the origin of fMRI brain activation to construct the “N-step reachable structural maps” - a series of vectors that denote signal propagation after the N-th step tracing of the structural connectivity matrix from the origin. Linear combinations of these N-step reachable maps are then used to predict the observed fMRI signals. Additionally, a utilization matrix, which numerically estimates whether the inclusion of a specific structural connection better predicts fMRI, can be constructed and estimated using simulated annealing. We tested the proposed framework on a visual fMRI task in a sample of body dysmorphic disorder (BDD) subjects and comparable healthy controls. Group differences were inferred by comparing the observed utilization differences against 10,000 permutations under the null hypothesis. Results revealed that BDD subjects under-utilized several key local connections in the visual system, which may help explain previously reported fMRI findings and further elucidate the underlying pathophysiology of BDD.
How the brain is organized into functional networks remains an elusive question. Despite widespread use of functional neuroimaging, inferences on network properties have largely been based on associative activation patterns; its limited temporal resolution hinders an understanding of dynamic interactions between spatially defined nodes. Advanced modeling techniques such as Dynamic Causal Modeling  provide estimates of complex neuronal dynamic interactions estimated from BOLD signal patterns, but may be practical for understanding only a small number of nodes.
Inferences about network dynamical interactions on a larger scale could potentially be derived using strategies involving structural to functional mapping. This is based on the assumption that a system's structural connectivity pattern provides a framework for determining information available as inputs from certain regions and its outputs and potential influence on other regions . Such an approach could utilize a combination of DTI-tractography to estimate brain white matter connectivity and functional magnetic resonance imaging (fMRI) to estimate the neuronal activity coupled to blood flow changes in anatomical regions that comprise nodes of the network.
Several structural to functional mapping approaches have been proposed [2-9]. Some have focused on specific, limited regional activation patterns [2, 3] while others have developed models for understanding functional connections in larger set of regions comprising networks or systems [4-9]. However, to our knowledge no study has used structural to functional mapping to estimate information transfer on a systems-wide level. We therefore developed a model to parsimoniously explain overall observed activation patterns with the assumption that the system is comprised of a set of interacting nodes.
To test this approach, we applied it to understanding activation patterns within the visual system. The visual system has been well characterized in monkeys and humans and can be reliably activated with functional neuroimaging paradigms. In addition, previous studies have provided models for a temporal sequence of activation, particularly for the “first wave” of information transfer (see, for example [10, 11]), from V1 to V2 and other visual cortical areas, followed by temporal and parietal regions. This allows us to make several simplifying assumptions about direction of information transfer. We tested our model in a dataset consisting of healthy human controls and individuals with body dysmorphic disorder (BDD). Individuals with BDD have perceptual distortions in which they misperceived defects in their appearances . Previous fMRI studies in BDD found abnormal visual processing in primary and secondary visual processing systems for face and inanimate object (house) stimuli [13-15]. To better characterize these findings, we thus developed a strategy for functional by structural hierarchical mapping to estimate step-wise activation patterns in the visual system. Here, the word “hierarchical” refers to the n-step structural maps used for model fitting. The fMRI portion of this dataset, which utilized well-characterized house visual stimuli, was previously analyzed for group differences using the general linear model . Here we demonstrate a model of hierarchical structural to functional mapping that, under several simplifying assumptions, elaborates on previous results by providing estimates of structural connection utilization within the visual system.
Eleven unmedicated participants with BDD and 13 healthy controls, aged 20 to 48 years, were scanned in a 3-T Siemens Allegra MRI scanner. Diffusion-weighted MR imaging (DWI) data were acquired using single-shot spin-echo echo-planar imaging (EPI) (field of view=240mm; voxel size=2.5×2.5×3.0mm, with 0.75 mm gap; TR/TE=7400/96ms; flip angle 9°). We collected 44 contiguous axial slices aligned to the AC-PC line along 34 gradient-sensitizing directions with b=1000s/mm2 and one minimally diffusion-weighted scan. DWI data were corrected for eddy current artifacts using FSL toolbox (http://www.fmrib.ox.ac.uk/fsl/fdt/fdt_eddy.html). fMRI data came from a previous object-processing fMRI study; subject selection, acquisition and preprocessing are as previously described . fMRI stimuli consisted of photographs of houses that contained only low spatial frequency (LSF) information , in order to more sensitively probe configural and holistic visual processing.
We computed whole-brain deterministic DTI tractography using Diffusion Toolkit with the Fiber Assignment by Continuous Tracking (FACT) algorithm  and a maximum turn angle of 35 degrees. Cortical regions of interest (ROIs) were defined using the Harvard Oxford Cortical probabilistic atlas. The masks were set to a liberal threshold of 10% to allow for the inclusion of tissue along the gray-white matter interface, where DTI tractography estimates are most reliable . To register these ROIs to each subject's DTI space, we used FSL's FLIRT program  to determine the optimal affine transformation between the MNI152 T1 average brain (in which the Harvard Oxford probabilistic atlases are based) and each subject's unique FA image, using a mutual information-based cost function. We applied the resulting transformation to register the 26 visual pathway related ROIs to each subject's DTI space. For each pair of ROIs, the number of fibers connecting them was counted to construct a structural connectivity matrix (a fiber was considered to connect two ROIs if it originated in the first ROI and terminated in the second, or vice versa). Thus, this matrix was symmetric, with the diagonal entries set to zero (no self connections).
To generate functional MRI signals for each region of interest, we computed the percentage signal change using the Featquery tool in FSL to estimate BOLD signal differences between the house stimuli and baseline.
Several assumptions and simplifications are needed in order to perform FSH mapping, outlined step-by-step as follows:
Now, we illustrate how to set up the FSH mapping using the BDD house processing example. In step 1, a total of 26 visual system ROIs were carefully selected from the group mean fMRI activation maps overlaid on the Harvard-Oxford cortical probabilistic atlas for house stimuli viewing vs. baseline, for combined BDD and healthy control groups. These 26 regions are the (bilateral) lateral occipital cortex superior division (nodes 1,2), lateral occipital cortex inferior division (3,4), intracalcarine cortex (5,6), precuneus cortex (7,8), cuneal cortex (9,10), parahippocampal gyrus posterior division (11,12), lingual gyrus (13,14), temporal fusiform cortex anterior division (15,16), temporal fusiform cortex posterior division (17,18), temporal occipital fusiform cortex (19,20), occipital fusiform cortex (21,22), supracalcarine cortex (23,24), and occipital pole (25,26).
As the fMRI task was a visual task, we thus assumed that cortical fMRI signals originated in the bilateral primary visual cortex (V1 or intracalcarine cortex; nodes 5 and 6). Mathematically, this is represented by a 26-element “activation” column vector A, whose entries are all zeros except for bilateral V1 (where A takes values of 1). In steps 2-4, we first observe that 0-step reachable map is A itself, and define the 1-step reachable map as (U C) (here the circle operator denotes the entry-wise or Hadamard product of two matrices of the same dimension). The n-step reachable structural maps (An) were then determined by the iterative application of the following procedure:
A0 = A, Co = C; for n = 0, 1, …,
Lastly, we fit the following linear equation so that the entire set of observed fMRI activations in the visual system, S, can be explained by linear combinations of An
Here, S is the fMRI signal column vector (26×1) and e is the error term. The optimization of the FSH mapping is achieved in two parts. First, U was initialized as a 26×26 matrix whose every entry took a value of 1 (i.e., all structural connections were utilized). We then fixed U and the linear equation above was fitted to determine the weights (w0, w1, w2, etc in the above equation) and to calculate the fitting errors. For the second part of the optimization, the weights were fixed, and data from all subjects in each group were pooled together to estimate a group-level optimal U using the simulated annealing algorithm (thus yielding two U matrices, one for the BDD and one for the control group). To this end, we perturbed U (from state i to state j) by randomly picking one element in U and changing its value (between 0 and 1). The acceptance criterion determined whether the new state j was accepted from the current state i by applying the following decision rule with respect to an artificial cooling temperature (c).
Where e(i) and e(j) are the linear fitting residuals. This perturbation was repeated numerous times until the solution space was adequately sampled and the global minimum reached.
As structural connection strengths vary across edges, utilization differences for connections that have on average higher fiber counts should be considered potentially more significant. As a result, we devised a statistic dU, by summing up all rows and columns of the difference matrix (subtracting the U of the control group from that of the BDD) weighted by the mean connectivity matrix C̄ (structural connectivity matrices averaged across all study participants).
Here, dU measures the overall group utilization differences, and permutation testing can be conducted to determine whether the observed dU reaches statistical significance. To this end, we performed 10,000 permutations by randomly re-assigning each subject's diagnosis (into two groups of size 11 and 13). At each permutation, the test statistic dU was re-computed and recorded, and the observed dU was then ranked against the re-sampled dUs. If the observed dU was relatively extreme and ranked among the top 5%, then we determined that the observed utilization group difference has reached statistical significance at a p value of 0.05.
Linear fitting errors vs. the level of n-step structural maps used for model fitting for one subject is shown in Figure 1. Figure 2 summarizes fitting results for all subjects in both groups. Linear combinations of 1-, 2-, and 3-step structural maps were tested. Three-level structural maps provided sufficient linear fitting, which was further improved with simulated annealing.
Group differences in the binary utilization matrices (BDD-control) weighted by the mean connectivity strength (note that this matrix is symmetrical, and thus its (i, j) and (j,i) elements indicate the utilization difference for the same structural connection) are shown in Figure 3. According to this model, the control group utilized several more structural connections than the BDD group (the most significant differences, marked A to C in Figure 3, are structural connections between the left intracalcarine cortex and left lingual gyrus, the right intracalcarine cortex and right lingual gyrus, and the right temporal occipital fusiform cortex and right temporal fusiform cortex, posterior division. These connections are visualized in Figure 4. Permutation tests confirmed that the observed utilization group difference reached statistical significance (p = 0.0329).
These results are consistent with previous fMRI results for house processing, which found relative hypoactivity in the BDD group vs. controls in the lingual gyrus, precuneus, and parahippocampal gyrus . Moreover, the use of FSH mapping provides additional information about utilization of structural connections that may explain observed fMRI differences; an inference from this is that individuals with BDD may have abnormally low information transfer between both primary and secondary visual cortical regions, and within higher-order temporal lobe visual processing systems. Results must be considered in light of multiple simplifying assumptions made to aid in this initial model development, and in light of small sample size. Future models that account for feedback, inhibitory interactions, and deactivation are necessary for more accurate estimations. Nevertheless, FSH may serve as a useful tool for integrating structural connectivity with fMRI data for improved understanding of functional activation patterns.