Relating structural connectivity with functional activity is fundamentally important to understanding the brain's physiology. The thalamocortical system serves as a good model system for exploring structure/function relationships because of its well-documented anatomical connectivity. Here we performed functional and structural magnetic resonance mapping of the human thalamocortical system using intrinsic brain activity and diffusion-weighted imaging. The accuracy of these imaging techniques is tested by comparison with human histology registered to common anatomical space and connectional anatomy derived from nonhuman primates. In general, there is good overall concordance among structural, functional, and histological results which suggests that a simple model of direct anatomical connectivity between the cerebral cortex and the thalamus is capable of explaining much of the observed correlations in neuronal activity. However, important differences between structural and functional mapping results are also manifest which suggests a more complex interpretation and emphasizes the unique contributions from structural and functional mapping.
The functional repertoire of any system is ultimately determined by its structural composition. Equally important in the brain, the underlying structure is continually reshaped by function in relation to experience. This interdependence underscores the importance of exploring the relationship between structure and function in neuroscience. Cellular electrophysiology provides a window on function but only with sparse sampling. Likewise, invasive tract tracing, especially performed postmortem, is similarly limited in its areal coverage. Recent developments have enabled whole brain characterization of both neuronal activity using functional magnetic resonance imaging (fMRI) and large fiber tracts using diffusion tensor imaging (DTI). The combination of these two imaging modalities permits a systems level interrogation of functional interactions and their structural underpinnings. However, most of the existing dual imaging studies are limited in the scope of their functional measures by reliance on task performance to elicit measurable neuronal activity (Toosy et al. 2004; Baird et al. 2005) (for a review of this extensive literature; see Rykhlevskaia et al. 2008). First, each task activates a very restricted set of brain regions, and second, a “wired together” interpretation of coactivations must be dissociated from the confound of stimulus driven synchrony. If the goal is a global characterization of structure/function relationships in the brain, reliance on task performance may not be the most efficient approach.
An alternative to the task-based approach to studying the brain's functional organization is to examine the organization of its ongoing functional activity. A factor motivating this approach is the fact that the majority of the brain's enormous energy budget is devoted to its ongoing functional activity rather than evoked activity (Raichle and Mintun 2006). In fact, less than 5% of its energy budget is devoted to momentary task demands (Raichle and Mintun 2006). Ongoing fluctuations in the blood oxygen level–dependent (BOLD) fMRI signal provide a unique view of the organization of intrinsic activity. The ongoing activity revealed by the fMRI BOLD signal persists across different states of consciousness including wakefulness (Biswal et al. 1995; Lowe et al. 1998; Cordes et al. 2000; Greicius et al. 2003), sleep (Horovitz et al. 2008; Larson-Prior et al. 2009), anesthesia induced unresponsiveness (Johnston et al. 2008), and coma (Boly et al. 2008). Ongoing activity exhibits correlated signaling within neuroanatomically defined systems that replicate with remarkable precision responses to experimentally administered tasks (Vincent et al. 2007). In this way, multiple functional networks can be delineated from a single scan without the need for any task performance (Fox and Raichle 2007).
The focus of the current study is to map structural and functional connections between the human thalamus and cortex. The thalamus is a core structure of the brain that contains primary relay nuclei that are highly specific in their connections with distinct zones of the cerebral cortex. A seminal study has previously investigated the structural connectivity of the human thalamocortical system using DTI and probabilistic tractography (Behrens, Johansen-Berg, et al. 2003). A specific question we were interested in was how well structural information derived from DTI compares with functional information derived from fMRI of intrinsic brain activity (Zhang et al. 2008). Furthermore, how well does structural and functional MR mapping compare with “gold standard” histological parcellation of the human thalamus and connectional anatomy inferred from non-human primates? We hypothesized that, in general, correlations in intrinsic neuronal activity should correspond well with connectional anatomy derived from DTI, but important differences may be observed given that complex neuronal interactions occur via polysynaptic pathways and given that the strength of an anatomical connection is not always reflective of the importance of its functional role (Van Horn et al. 2000; Sherman and Guillery 2006). Therefore, in considering thalamocortical “connectivity”, defining a functional connectivity map may be just as important as defining a structural connectivity map of the thalamocortical system, each providing unique information relevant to the physiology of the system.
Materials and Methods
Intrinsic brain activity was assessed with BOLD fMRI in a group of 17 subjects (3T Siemens Allegra MR scanner [Erlangen, Germany], 4-mm isotropic voxels, time echo [TE] = 25 ms, time repetition [TR] = 2.16 s, 776 frames per subject). Subjects visually fixated on a projected cross-hair. No other task was imposed except to remain still and not fall asleep. Structural fiber tracks were assessed using diffusion-weighted MRI (DWI) acquired in a group of 12 subjects using a 1.5T Siemens Sonata MR scanner. Diffusion was measured in 48 directions plus 4 unsensitized volumes (2.0- to 2.5-mm isotropic voxels, TE = 113 ms, TR = 7 s, b value = 400–1200 s/mm2). Five DTI data sets were acquired in each individual. Total acquisition time was approximately 90 min per subject. All data were acquired in the course of previously published and unpublished studies (Shimony et al. 2006; Fox et al. 2007). Although structural and functional imaging were performed in two separate groups, it is unlikely that the current study would benefit from dual modality imaging on the same individuals due to methodological variability in imaging a small subcortical structure at the single subject level (see Supplemental Materials for an extended discussion). Dual modality imaging in a separate population of 3 subjects supports this assertion (data not shown). In the current results, all data presented were analyzed at the population level. Anatomical images acquired for definitive atlas transformation included a T1-weighted magnetization–prepared rapid gradient echo (MP–RAGE) and a T2-weighted (T2W) fast spin echo scan.
Cortical Regions of Interest Definition
Five cortical regions were manually drawn on the cortical surface using CARET (Van Essen 2005; Van Essen et al. 2001) (see Zhang et al. 2008) on the basis of major sulcal landmarks, largely following previous work taking into account the known anatomical connectivity of the thalamus and the cortex (Behrens, Johansen-Berg, et al. 2003; Zhang et al. 2008). These cortical regions of interest (ROIs) were defined as 1) fronto-polar and frontal cortex including the orbital surface and anterior cingulate; 2) motor and premotor cortex (Brodmann areas 6 and 4) excluding adjacent portions of cingulate cortex; 3) somatosensory cortex (Brodmann areas 3, 1, 2, 5, and parts of 40); 4) parietal and occipital cortex including posterior cingulate and lingual gyrus; 5) temporal cortex including the lateral surface, temporal pole and parahippocampal areas (Fig. 1a). For our functional connectivity MRI (fcMRI) study, BOLD data were first registered to our standard atlas space which ensured correct alignment with our cortical ROIs. For our DWI study, tractography was performed in native acquisition space. Therefore, our cortical ROIs were first deformed to each subject's anatomy using CARET's surface-based registration (Van Essen and Drury 1997; Van Essen et al. 2001), and then converted to volume space with an assigned cortical thickness of 3 mm (1.5 mm above and below fiducial surface corresponding to “layer IV”).
Resting-State fcMRI Preprocessing and Analysis
Preprocessing included compensation for systematic, slice-dependent time shifts, elimination of systematic odd-even slice intensity differences due to interleaved acquisition, rigid body correction for head motion within and across runs, and normalization of the signal intensity across each run (skipping the first four frames) to obtain a whole brain mode value of 1000 (Ojemann et al. 1997). Atlas transformation was achieved by composition of affine transforms connecting the first functional volume (averaged over all fMRI runs after cross-run realignment) with the T2W and T1-weighted structural images. Common mode image registration (e.g., to correct for head motion within and across fMRI runs and affine warping of T1-weighted structural images to the atlas-representative target image) was performed using in-house versions of standard algorithms (Woods et al. 1998). Cross-modal registration (e.g., T2W to T1-weighted images) was performed using the vector gradient measure maximization (Rowland et al. 2005). Our atlas representative template includes MP–RAGE data from 12 normal individuals and was made to conform to the 1988 Talairach atlas (Talairach and Tournoux 1988; Lancaster et al. 1995). Linear trends across runs were removed voxelwise and the data were low-passed to retain frequencies below 0.1 Hz. Several sources of spurious variance were removed by regression of the following nuisance variables along with their first derivatives: 1) the six parameters resulting from rigid body correction for head motion; 2) a signal from a ventricular region of interest; and 3) a signal from a white matter region of interest (Fox et al. 2005). To prepare our BOLD data for the present analyses, each fMRI run was resampled to 2-mm3 voxels in atlas space.
To map correlations in intrinsic brain activity between the cortical ROIs and the thalamus, the following analysis was performed. First, the average BOLD time course was extracted from each cortical ROI. Then, the partial correlation was computed between each of these cortical ROI time course and the time course extracted from each thalamic voxel, eliminating the shared variance among all the cortical ROI time courses (Zhang et al. 2008). For purposes of calculating statistical significance, partial correlation coefficients were converted to a normal distribution using Fisher's r-to-z transform. Because consecutive BOLD frame are not independent but autocorrelated, it was necessary to correct the degrees of freedom according to Bartlett's theory (Jenkins and Watts 1968). Accordingly, for purposes of significance testing under a fixed-effects model, the z transformed partial correlation values were converted to Z scores by multiplying by , where n is the degrees of freedom (number of frames divided by 3.183) and k is the order of the partial correlation (in the present analysis k = 4) (Weatherburn 1949). Z score maps were combined across subjects using a fixed-effects analysis. Thresholding for displayed Figures 1 and 4 was set to half of the maximum partial correlation value for the given thalamic map. Thus, figures are labeled as 50% and 100% corresponding to half maximum and maximum partial correlation values. To reduce bias due to possible outlier values, the maximum partial correlation value was defined as the 99th percentile value instead of the 100th percentile value. Voxels displayed using a 50% threshold all pass a multiple comparisons criteria (false positive error rate of 5%) empirically determined using a bootstrap procedure (see Zhang et al. 2008 for details).
DWI Preprocessing and Analysis
Each DWI data set was assembled by collating slices from two interleaved scans. Head motion correction was performed according to the methods of referenced work (Shimony et al. 2006). All DWI computations were performed in native acquisition space thereby minimizing resampling of the data. Subsequent analysis was performed with FMRIB's Diffusion Toolbox (FDT) as part of the FSL package (FSL v4.1.2, FDT v2.0, Oxford, UK) (Behrens, Woolrich, et al. 2003; Smith et al. 2004; Behrens et al. 2007; Woolrich et al. 2009). First, eddy current correction was applied. Then, the local probability density function of a diffusion tensor model was calculated using an algorithm that models intravoxel crossing fibers (Behrens et al. 2007). Probabilistic tractography was performed at a global level starting from each voxel in the thalamus and then tracked to its ipsilateral cortex (Behrens, Johansen-Berg, et al. 2003; Behrens, Woolrich, et al. 2003). The probability of connectivity from a thalamic voxel to a cortical ROI was calculated to be the number of fiber tracts successfully propagated to a given cortical ROI divided by the total number of tracts propagated to all cortical ROIs. This last preprocessing step is not normally performed in probabilistic tractography (see Behrens et al. 2007; Behrens, Johansen-Berg, et al. 2003; Johansen-Berg et al. 2005, for examples). However, it should be noted that inclusion of this step results in higher overlap for all cortical ROI measures both in fcMRI and in DTI results (i.e., all values in Table 1). Therefore, there is methodological value to this step as well as the possibility that the ratios of probabilistic tractography values are, to a degree, reflected in the cofluctuation of neuronal activity.
|Cortical partition||Percentage of overlapping voxels|
|fcMRI map||DTI map|
|Cortical partition||Percentage of overlapping voxels|
|fcMRI map||DTI map|
Post-tractography intersubject spatial alignment was performed for group level analysis. The spatial relationships among diffusion volumes, atlas representative target, T1-weighted and T2W images were computed using previously published methods (Shimony et al. 2006). Population level probability maps were calculated using a fixed effects model. The same thresholding methodology applied to functional mapping results was applied to the population level DWI results. Thus, DWI maps display values between approximately 50% and 100% of maximum calculated probability value. For display purposes, voxel boundaries for both fcMRI and DTI results were smoothed using fourfold interpolation and then displayed using in-house software written on the MATLAB platform (The MathWorks; Natick, MA). The T1-weighted anatomical image displayed behind our MR results is the Montréal Neurological Institute's (MNI) 152 subject average registered to common space using non-linear transform and distributed with FSL (Jenkinson and Smith 2001; Smith et al. 2004; Woolrich et al. 2009). We performed affine transform to register the MNI152 template to our standard atlas target.
3D Reconstruction of the Human Brain using 2D Histological Slices
The Atlas of the Human Brain in Stereotaxic Space by Mai, Paxinos, Voxx (hereafter referred to as the “Mai atlas”) provided serial coronal sections of a 24-year-old male's brain (Mai et al. 2008). In the generation of this atlas, a variety of morphometric measures and staining techniques were utilized to partition and label various subcortical as well as cortical structures. We reconstructed the 3-dimensional brain volume from these 2-dimensional histological slices using 1-dimensional interpolation with a lookup table of slice coordinates in mm units along the y-axis (implemented using MATLAB's interpolation algorithm, The Mathworks). To facilitate comparison our functional and DWI MR results with this atlas, we spatially registered all images to a common atlas representative target. Registration of functional and DWI images are described in the above respective sections. Cross modal registration of the Mai atlas with our T1 atlas was performed using vector gradient measure maximization (Rowland et al. 2005). Spatial relationships between the source and target images were defined in terms of affine transforms.
Monte Carlo Calculation of Significance of Overlap between fcMRI and DTI Results
The simulation consisted of recreating all thalamic voxels, assigning each voxel a binary value (Fig. 1D). The proportion of voxels assigned a value of 1 instead of 0 is determined by the number of voxels that exceeded the 50% threshold in the MR results (Fig. 1B,C). In this way, the thalamus was simulated for both the fcMRI results and the DTI results. Overlap between fcMRI and DTI maps was calculated as the number of simulated voxels that contained a value of 1 in both maps. The simulation was repeated 10 000 times, each time randomly permuting which voxels of the thalamus were assigned a value of 1. This null distribution was then compared with the empirical number of voxels in the overlap from Figure 1D to calculate statistical significance. The above steps were performed 5 separate times, once for each cortical ROI.
Five cortical regions shown in Figure 1A were defined based on major sulcal landmarks (Van Essen 2005; Van Essen et al. 2001), taking into account the anatomical connectivity between the thalamus and the cortex. For our functional study, the average BOLD signal was extracted from each cortical region and compared with the BOLD signal at each voxel in the thalamus using partial correlation to eliminate shared variance among all cortical regions (Zhang et al. 2008). Figure 1B demonstrates the specificity of these functional connections. Each cortical region was specifically correlated with distinct, spatially restricted zones within the thalamus. For our structural study using DTI, we adopted a parallel approach of relating the thalamus to its ipsilateral cortex using the implementation of Behrens, Johansen-Berg, et al. (2003; Behrens et al. 2007). Distinct zones of thalamic connectivity were observed for each cortical region in a distribution similar to that of Behrens, Johansen-Berg, et al. (2003) (Fig. 1C). The overlay between fcMRI and DTI results is shown in Figure 1D. Even though the physical basis of the two MR measurements is vastly different, there is significant overlap in their connectivity profiles (P < .00001). Table 1 shows the percentage of thresholded voxels that overlap between the two methods. Additionally, multiple thalamic maps were condensed into a single winner-take-all projection (Fig. 2) in which every thalamic voxel is labeled with the color of the cortical ROI that generated the highest partial correlation or highest tractography probability value using fcMRI and probabilistic DTI, respectively.
Evaluation of Connectivity Overlap
For an evaluation of the accuracy of our imaging results compared with histology, we constructed a 3-dimensional interpolation of a coronal atlas of thalamic histology created by Mai et al. (2008). This atlas was registered to our standard anatomical space by performing cross modal matching of gyral and sulcal contours which allowed for comparisons at a voxelwise level with our MR imaging results, notwithstanding subject-to-subject variability from an atlas characteristic of a single individual. Figure 3A shows a representative axial slice of labeled thalamic histology. In comparing our MR results with histology, we first evaluated the localization of the overlap area between the two methods (Fig. 1D, purple). The overlap between motor/premotor fcMRI and DTI results traced on top of the Mai atlas (Fig. 3B) highlights the ventral lateral (VL) and ventral lateral posterior (VLP) thalamic nuclei, which corresponds well with documented connectional neuroanatomy (Ray and Price 1993; Jones 2007). Thalamic connectivity with other cortical regions similarly demonstrated close correspondence with neuroanatomy. The area of overlap between the two methods with prefrontal cortex localized to mediodorsal (MD) nucleus as well as ventral anterior nuclei and parts of the anterior group of thalamic nuclei (Fig. 3C). Overlap with temporal cortex localized well with the medial geniculate nucleus (MGN, Fig. 3D) and medial pulvinar (Supplemental Fig. 1). Overlap with parietal and occipital cortex localized predominantly with lateral pulvinar (Fig. 3E). fcMRI results also localized to the lateral geniculate nucleus (LGN), whereas DTI localization to LGN was not as pronounced (see Supplemental Fig. 1 for localization as well as a discussion of this difference). Both methods demonstrated connectivity between the somatosensory cortex and relevant parts of the thalamus. The overlap between the two methods localized predominantly to anterior pulvinar, which projects mainly to Brodmann area 5 (BA5) (Fig. 3F). DTI results also highlighted ventral posterior lateral (VPL) nuclei, whereas fcMRI results were localized more medially to just anterior pulvinar (Supplemental Fig. 1).
Evaluation of Connectivity Differences
The greatest difference (by total voxel count) in the present thalamic localization results was observed in the lateral extent of the MD nucleus (Fig. 4, crosshairs). Although neuronal activity in motor/premotor cortex was correlated with nuclei VL/VLP, the strongest functional correlations were observed more medially in nucleus MD (Fig. 4A,C). In contrast, MD showed highest probability of DTI connectivity with prefrontal cortex (Fig. 4B–D), which is concordant with connectional anatomy derived from invasive track tracing studies in monkeys (Jones 2007). This difference is not the consequence of the level of thresholding and changes in boundary size with different levels of thresholding because the strongest functional correlations occur at a locus of the thalamus that is a transition zone in the tractography results from high connectivity with motor/premotor to high connectivity with prefrontal cortex.
We investigated this difference from the perspective of our cortical ROIs. Because each anatomically defined cortical ROI includes extended regions of gray matter that encompass multiple functional systems, it is likely that the average signal extracted from each ROI does not uniformly correlate with all voxels within that ROI. The same partial correlation analysis that generated maps in Figure 1B was now extended to the cortex. Figure 4E shows the results of this analysis. As one can see, most of the maps show cortical correlations that are well-restricted to within the cortical ROIs that generated the map. However, the prefrontal ROI correlation map not only highlights the prefrontal cortex but clearly includes a relatively intact default mode network (DMN), the only functional network to prominently retain its long range functional connections. This is to be contrasted with prefrontal connectivity with DTI tractography which predominantly tracks to lateral and not medial prefrontal cortex (i.e., outside of the DMN) because of crossing fiber issues from the corpus callosum. Motor/premotor correlations also show hints of functional system contribution in highlighting the parietal component of the dorsal attention system, though this effect was much weaker. It should be noted that the areas of largest difference between fcMRI and DTI tractography (Fig. 4A,B) concern the prefrontal and motor/premotor cortex.
Areas of highest probability in the DTI results tend to localize near the surface of the thalamus where anisotropy is highest, whereas functional localization does not preferentially localize to the outer thalamic surface (Fig. 1B,C). No consistent systematic shifts were observed between the fcMRI and DTI results, ruling out the possibility of systematic biases due to anatomical registration. Connectivity between LGN and parietal/occipital cortex was clearly seen with fcMRI but less clear with DTI, as noted above (Supplemental Fig. 1). On the other hand, DTI tractography defined connections between motor/premotor cortex and posterior portions of VLP, whereas functional connectivity results were more restricted to anterior and medial parts of the thalamus (Fig. 1B,C, third column). As noted above, somatosensory connectivity with nuclei VPL was observed with structural imaging but absent with functional imaging (Fig. 1B,C, fourth column). Functional correlations with somatosensory cortex were predominantly circumscribed to the anterior pulvinar (Fig. 1B, fourth column). A closer look at the fcMRI results in region VPL showed noticeably weak functional correlations with all cortical ROIs, although motor/premotor cortex still demonstrated statistically significant values in this region (Fig. 2).
The high overall correspondence between the structural and functional results suggests two important implications. First, the substantial overlap between observed structure and function along with their correspondence with histology supports the feasibility and validity of the two imaging approaches in mapping thalamocortical connectivity. In the current work, measurements of structural and functional connectivity within the thalamocortical system were performed independently. This is to be contrasted with a large body of existing work where task-evoked fMRI is used to define a priori ROIs for use in constraining tractography results, in essence asking the binary question of whether a structural connection exists given a set of functional activation regions (see Rykhlevskaia et al. 2008) for a review of this literature). Here, the pattern of connectivity generated with one method is not influenced by the other method. It is only at the postprocessing stage that the results from both approaches are combined for comparison purposes. It should be noted that our present results extend the findings of a previous study which demonstrated the correspondence between DTI tractography and task-evoked fMRI in select thalamic nuclei (Johansen-Berg et al. 2005).
The second implication of the substantial overlap in observed structure and function is that a straightforward model of neuronal activity reflecting more direct physical connectivity is capable of explaining much of the mapping results. In correlation mapping of fMRI activity, it is usually an open question whether the observed correlation is a result of more direct (fewer synaptic junctions) or less direct connections (more synaptic junctions). In the present analysis, we can directly modify the relative contribution of these two parameters by using partial correlation to eliminate shared variance among the five cortical ROIs, thereby mitigating the contributions of cortico-cortical interactions in the generation of thalamic correlation maps. At a global level, the average BOLD signal across all voxels in the brain is significantly correlated with much of gray matter, including the thalamus (Fox et al. 2009). Partial correlation also eliminates globally shared signals, whether of neuronal or non-neuronal origin (Birn et al. 2006). The result of this analysis step is a delineation of neuronal activity that is reflective of more direct physical connectivity and therefore corresponds better with DTI tractography results which predominantly track long distance fiber pathways with few intervening synapses.
Interpretation of Connectivity Differences
Areas of observed differences clearly exist in the present results, suggesting a more complex interpretation. The greatest difference by voxel count is in connectivity with motor/premotor cortex (Fig. 4). There is a shift in the approximate boundary between motor/premotor regions and prefrontal regions of the thalamus, as well as a large associated shift in peak foci of localization (Fig. 4A,B). Figure 4E showed that this effect may be related to spatial nonuniformity of correlations with each cortical ROI. Despite the use of an anatomical parcellation scheme, the DMN emerged as the principal functional system to retain its long range functional connections. The robustness of this functional system may be related to graph theoretic descriptions of hub-like properties of some of its canonical nodes (Hagmann et al. 2008).
The functional cohesiveness of the DMN shown here suggests the distinctiveness of this functional network and also suggests the presence of polysynaptic contributions to our functional results. From a cortical parcellation perspective, these polysynaptic influences would likely be better controlled for in our partial correlation strategy if the size of our ROIs were decreased (thus possibly increasing the correspondence between fcMRI and DTI tractography results), although we would lose whole brain coverage if the number of our ROIs were not dramatically increased (note that increasing the ROI number to double digits would not be ideal given a partial correlation strategy because of the reduction in signal to noise; see Zhang et al. 2008). Other polysynaptic interactions also likely contribute to our observations in ways that we do not control for in the present analysis. For example there are known subcortico-subcortical interactions such as those mediated by the thalamic reticular nucleus (Guillery and Harting 2003) and the brainstem reticular formation (Jones 2008). In contrast, DTI tractography results generated in the current study predominantly reflect more direct white matter pathways between the thalamus and cortex. Inference on less direct pathways such as those involving cortico-cortical connections can be performed with multi-ROI/graph theoretic measures, but the accuracy of these reconstructed polysynaptic pathways scales inversely with the number of synapses (i.e., graph nodes) involved in the pathway and the acquisition resolution. Thus, in a comparison of the two methods, one should be aware of inherent limitations in inferring functional connectivity strictly on the basis of tractography alone and vice versa. However, the limitation of each lone standing method highlights the complementary nature of the two imaging modalities. Empirically derived measures of functional connectivity can serve as a physiological constraint on the extensive permutation of possible polysynaptic anatomical connections derived from DTI tractography. Likewise, DTI tractography can often be used to generate the most plausible anatomical tract connecting two brain regions that have correlated neuronal activity (see Honey et al. 2009 for an excellent example of a graph theoretic approach to combining fcMRI and DTI tractography results).
Methodological idiosyncrasies also contribute to observed differences. In the present analysis, we utilized probabilistic tractography which is able to track through areas of the brain near gray matter (Behrens, Johansen-Berg, et al. 2003). This differs from the very limited number of studies that have combined fcMRI with diffusion weighted imaging (Koch et al. 2002; Skudlarski et al. 2008; van den Heuvel et al. 2008, 2009; Greicius et al. 2009; Honey et al. 2009), all of which have used conventional streamline tractography which produce reliable results only in brain areas where anisotropy is high (Behrens, Johansen-Berg, et al. 2003). This distinction in tractography methodology is an important point of consideration and represents a crucial step in the feasibility of thalamocortical structural mapping (Behrens, Johansen-Berg, et al. 2003). However, even with this technical advance in the science of tractography, there are still well known and well-described current limitations (Behrens et al. 2007). The probabilistic algorithm utilized, in part, accounts for crossing fibers (Behrens et al. 2007), but areas with major fiber bundles still tend to suppress tracking of smaller fiber bundles crossing perpendicular to them. For example, the MGN can be distinguished from LGN based on differences in their connectivity profiles but showing that the connections from MGN reliably and accurately traverse to auditory cortex has been a much harder task due to crossing fibers from the optic radiations (Devlin et al. 2006; Behrens et al. 2007). In the current work, functional mapping from temporal cortex appeared to generate more specific localization to MGN than its DTI counterpart (Fig. 1B–D, far right panel), but the best definition of MGN from an anatomical point of view comes from the area of overlap between the two methods (Fig. 3D).
Methodological variability can be seen to have a pronounced effect in the winner-take-all projections (Fig. 2). Because this data reduction technique performs segmentation based only on the highest values, the boundaries of the WTA maps are expected to be more variable between the two methods and more prone to methodological idiosyncrasies. Notably, borders in DTI WTA maps are qualitatively more stripe-like, whereas fcMRI WTA borders better resemble irregular ellipsoids. This pattern is somewhat expected given that probabilistic tracts from deep thalamic nuclei have to physically pass through the outer surface of the thalamus on their way to the cortex, thus contributing to strip-like patterns. In this sense, inner and outer layers of the thalamus are less dependent on each other using functional imaging. Instead, the BOLD signal is strongest in gray matter and significantly attenuated in white matter structures.
The current choice of thresholding was chosen based on a percentage of the maximum correlation/probability value to better accommodate the differences in unit measurements from the two disparate methods. One advantage of our percentage thresholding approach is that it does not artificially place a constraint on the number of voxels displayed from each method. Nevertheless, for our fcMRI results, all displayed voxels satisfied a multiple comparisons statistical significance test at P < 0.05 determined by a bootstrap estimation on the data (Zhang et al. 2008). We, however, performed an alternative thresholding scheme in which the functional data were thresholded at a multiple comparisons corrected significance of P = 0.05 and then the DTI data were thresholded at a level that produced an equal number of displayed voxels. This alternative strategy produced results qualitatively similar to those presented in Figure 1 (data not shown), only with less degrees of freedom in the display parameters.
In the current work, we have demonstrated that structural and functional mapping of the thalamocortical system, performed independently of each other, demonstrate good overall correspondence. Areas of overlap between the two methods match well with histological partitioning of the human thalamus and connectional anatomy derived from non-human primates. Within the scope of the thalamocortical system, the correspondence in functional activity with more direct anatomical connectivity reinforces the notion that the thalamus is a relay structure of the brain that faithfully performs the duty of signal transmission. More generally, fcMRI and DTI will be increasingly utilized in combined structure/function studies. In this respect, the current results provide a validation of combined imaging approaches with comparison to a histological “gold standard”.
Differences between the two techniques do exist, and both results may be equally valid given the two approaches measure different aspects of “connectivity”. In light of these differences, caution is warranted when inferring functional interactions strictly on the basis of DTI tractography, and vice versa. Each technique provides unique and complementary information that should be interpreted in their corresponding context. Here, we have demonstrated the feasibility of combining these two approaches within the thalamocortical system. In the larger scope, an obvious complementarity exists between fcMRI and probabilistic DTI tractography toward a global understanding of the brain's physiology, both in health and disease.
National Institutes of Health (NIH) (NS06833) to M.E.R.; National Institutes of Mental Health (NIMH) (F30MH083483) to D.Z.; and NIH (K23 HD053212) to J.S.S.
Conflict of Interest: None declared.