Diffusion tensor fiber tracking of human brain connectivity: aquisition methods, reliability analysis and biological results.

We present a description, biological results and a reliability analysis for the method of diffusion tensor tracking (DTT) of white matter fiber pathways. In DTT, diffusion-tensor MRI (DT-MRI) data are collected and processed to visualize the line trajectories of fiber bundles within white matter (WM) pathways of living humans. A detailed description of the data acquisition is given. Technical aspects and experimental results are illustrated for the geniculo-calcarine tract with broad projections to visual cortex, occipital and parietal U-fibers, and the temporo-calcarine ventral pathway. To better understand sources of error and to optimize the method, accuracy and precision were analyzed by computer simulations. In the simulations, noisy DT-MRI data were computed that would be obtained for a WM pathway having a helical trajectory passing through gray matter. The error vector between the real and ideal track was computed, and random errors accumulated with the square root of track length consistent with a random-walk process. Random error was most dependent on signal-to-noise ratio, followed by number of averages, pathway anisotropy and voxel size, in decreasing order. Systematic error only occurred for a few conditions, and was most dependent on the stepping algorithm, anisotropy of the surrounding tissue, and non-equal voxel dimensions. Both random and systematic errors were typically below the voxel dimension. Other effects such as track rebound and track recovery also depended on experimental conditions. The methods, biological results and error analysis herein may improve the understanding and optimization of DTT for use in various applications in neuroscience and medicine.

Adaptive changes in early and late blind: a fMRI study of Braille reading.

Braille reading depends on remarkable adaptations that connect the somatosensory system to language. We hypothesized that the pattern of cortical activations in blind individuals reading Braille would reflect these adaptations. Activations in visual (occipital-temporal), frontal-language, and somatosensory cortex in blind individuals reading Braille were examined for evidence of differences relative to previously reported studies of sighted subjects reading print or receiving tactile stimulation. Nine congenitally blind and seven late-onset blind subjects were studied with fMRI as they covertly performed verb generation in response to reading Braille embossed nouns. The control task was reading the nonlexical Braille string "######". This study emphasized image analysis in individual subjects rather than pooled data. Group differences were examined by comparing magnitudes and spatial extent of activated regions first determined to be significant using the general linear model. The major adaptive change was robust activation of visual cortex despite the complete absence of vision in all subjects. This included foci in peri-calcarine, lingual, cuneus and fusiform cortex, and in the lateral and superior occipital gyri encompassing primary (V1), secondary (V2), and higher tier (VP, V4v, LO and possibly V3A) visual areas previously identified in sighted subjects. Subjects who never had vision differed from late blind subjects in showing even greater activity in occipital-temporal cortex, provisionally corresponding to V5/MT and V8. In addition, the early blind had stronger activation of occipital cortex located contralateral to the hand used for reading Braille. Responses in frontal and parietal cortex were nearly identical in both subject groups. There was no evidence of modifications in frontal cortex language areas (inferior frontal gyrus and dorsolateral prefrontal cortex). Surprisingly, there was also no evidence of an adaptive expansion of the somatosensory or primary motor cortex dedicated to the Braille reading finger(s). Lack of evidence for an expected enlargement of the somatosensory representation may have resulted from balanced tactile stimulation and gross motor demands during Braille reading of nouns and the control fields. Extensive engagement of visual cortex without vision is discussed in reference to the special demands of Braille reading. It is argued that these responses may represent critical language processing mechanisms normally present in visual cortex.

Normal brain maturation during childhood: developmental trends characterized with diffusion-tensor MR imaging.

PURPOSE: To characterize the maturational changes in water diffusion within central gray matter nuclei and central white matter pathways of the human brain by using diffusion-tensor magnetic resonance (MR) imaging. MATERIALS AND METHODS: Retrospective analysis of normal MR examination findings in 153 subjects (age range, 1 day to 11 years) referred for clinical neuroimaging was performed. All studies included diffusion tensor-encoded echo-planar MR imaging. Isotropic diffusion coefficient (D) and diffusion anisotropy (A(sigma)) were measured in the corpus callosum, internal capsule, caudate nucleus, lentiform nucleus, and thalamus. RESULTS: exhibited biexponential decay with age in gray and white matter regions, except for monoexponential decay in the genu of the corpus callosum. There was a steep nonlinear increase of A(sigma) in white matter tracts that paralleled the time course of the decline in D. In basal ganglia, only a small linear increase in A(sigma) was observed in patients. A(sigma) changes in the thalamus were intermediate between basal ganglia and white matter structures. CONCLUSION: Changes in magnitude and anisotropy of water diffusion follow stereotypical time courses during brain development that can be empirically described with multiexponential regression models, which suggests that quantitative scalar parameters derived from diffusion-tensor MR imaging may provide clinically useful developmental milestones for brain maturity.

Direct comparison of prefrontal cortex regions engaged by working and long-term memory tasks.

Neuroimaging studies have suggested the involvement of ventrolateral, dorsolateral, and frontopolar prefrontal cortex (PFC) regions in both working (WM) and long-term memory (LTM). The current study used functional magnetic resonance imaging (fMRI) to directly compare whether these PFC regions show selective activation associated with one memory domain. In a within-subjects design, subjects performed the n-back WM task (two-back condition) as well as LTM encoding (intentional memorization) and retrieval (yes-no recognition) tasks. Additionally, each task was performed with two different types of stimulus materials (familiar words, unfamiliar faces) in order to determine the influence of material-type vs task-type. A bilateral region of dorsolateral PFC (DL-PFC; BA 46/9) was found to be selectively activated during the two-back condition, consistent with a hypothesized role for this region in active maintenance and/or manipulation of information in WM. Left frontopolar PFC (FP-PFC) was also found to be selectively engaged during the two-back. Although FP-PFC activity has been previously associated with retrieval from LTM, no frontopolar regions were found to be selectively engaged by retrieval. Finally, lateralized ventrolateral PFC (VL-PFC) regions were found to be selectively engaged by material-type, but uninfluenced by task-type. These results highlight the importance of examining PFC activity across multiple memory domains, both for functionally differentiating PFC regions (e.g., task-selectivity vs material-selectivity in DL-PFC and VL-PFC) and for testing the applicability of memory domain-specific theories (e.g., FP-PFC in LTM retrieval).

The emotional modulation of cognitive processing: an fMRI study.

The functional neuroanatomy of visual processing of surface features of emotionally valenced pictorial stimuli was examined in normal human subjects using functional magnetic resonance imaging (fMRI). Pictorial stimuli were of two types: emotionally negative and neutral pictures. Task performance was slower for the negatively valenced than for the neutral pictures. Significant blood oxygen level dependent (BOLD) increases occurred in the medial and dorsolateral prefrontal cortex, midbrain, substantia innominata, and/or amygdala, and in the posterior cortical visual areas for both stimulus types. Increases were greater for the negatively valenced stimuli. While there was a small but significant BOLD decrease in the subgenual prefrontal cortex, which was larger in response to the negatively valenced pictures, there was an almost complete absence of other decreases prominently seen during the performance of demanding cognitive tasks [Shulman, G. L., Fiez, J. A., Corbetta, M., Buckner, R. L., Miezin, F. M., Raichle, M. E., & Petersen, S. E. (1997). Common blood flow changes across visual tasks: II. Decreases in cerebral cortex. Journal of Cognitive Neuroscience, 9, 648--663]. These results provide evidence that the emotional valence and arousing nature of stimuli used during the performance of an attention-demanding cognitive task are reflected in discernable, quantitative changes in the functional anatomy associated with task performance.

Differences between gray matter and white matter water diffusion in stroke: diffusion-tensor MR imaging in 12 patients.

PURPOSE: To investigate differences in water diffusion between white matter and gray matter in acute to early subacute stroke with diffusion-tensor magnetic resonance (MR) imaging. MATERIALS AND METHODS: Twelve patients with unilateral middle cerebral arterial infarcts were examined with diffusion tensor-encoded echo-planar MR imaging 17 hours to 5 days after stroke onset. Isotropic diffusion coefficient (D) and diffusion anisotropy (A(sigma)) images were computed. (D) values were measured in ischemic and contralateral gray matter and white matter by using A(sigma) images to differentiate white matter from gray matter. (D) images were compared with unidirectional and directionally averaged diffusion-weighted images. RESULTS: In all patients, (D) images showed two distinct levels of diffusion reduction in the infarct; more severe reduction occurred exclusively in white matter. (D) values were significantly less in infarcted white matter than in infarcted gray matter, whereas (D) values in the contralateral white matter and gray matter were not significantly different. Relative to the contralateral side, (D) values in the infarct were reduced by 46% in white matter and by 31% in gray matter (P <.001). Diffusion-weighted imaging caused underestimation of the magnitude and, in some cases, the spatial extent of the white matter diffusion abnormality. CONCLUSION: Isotropic diffusion is more reduced in white matter than in gray matter in acute to early subacute middle cerebral arterial stroke. Diffusion-tensor imaging may be more sensitive than diffusion-weighted imaging to white matter ischemia.

Areas involved in encoding and applying directional expectations to moving objects.

Two experiments used functional magnetic resonance imaging (fMRI) to examine the cortical areas involved in establishing an expectation about the direction of motion of an upcoming object and applying that expectation to the analysis of the object. In Experiment 1, subjects saw a stationary cue that either indicated the direction of motion of a subsequent test stimulus (directional cue) or provided no directional information (neutral cue). Their task was to detect the presence of coherent motion in the test stimulus. The stationary directional cue produced larger modulations than the neutral cue, with respect to a passive viewing baseline, both in motion-sensitive areas such as left MT+ and the anterior intraparietal sulcus, as well as motion-insensitive areas such as the posterior intraparietal sulcus and the junction of the left medial precentral sulcus and superior frontal sulcus. Experiment 2 used an event-related fMRI technique to separate signals during the cue period, in which the expectation was encoded and maintained, from signals during the subsequent test period, in which the expectation was applied to the test object. Cue period activations from a stationary, directional cue included many of the same motion-sensitive and -insensitive areas from Experiment 1 that produced directionally specific modulations. Prefrontal activations were not observed during the cue period, even though the stationary cue information had to be translated into a format appropriate for influencing motion detection, and this format was maintained for the duration of the cue period (approximately 5 sec).

Tracking neuronal fiber pathways in the living human brain.

Functional imaging with positron emission tomography and functional MRI has revolutionized studies of the human brain. Understanding the organization of brain systems, especially those used for cognition, remains limited, however, because no methods currently exist for noninvasive tracking of neuronal connections between functional regions [Crick, F. & Jones, E. (1993) Nature (London) 361, 109-110]. Detailed connectivities have been studied in animals through invasive tracer techniques, but these invasive studies cannot be done in humans, and animal results cannot always be extrapolated to human systems. We have developed noninvasive neuronal fiber tracking for use in living humans, utilizing the unique ability of MRI to characterize water diffusion. We reconstructed fiber trajectories throughout the brain by tracking the direction of fastest diffusion (the fiber direction) from a grid of seed points, and then selected tracks that join anatomically or functionally (functional MRI) defined regions. We demonstrate diffusion tracking of fiber bundles in a variety of white matter classes with examples in the corpus callosum, geniculo-calcarine, and subcortical association pathways. Tracks covered long distances, navigated through divergences and tight curves, and manifested topological separations in the geniculo-calcarine tract consistent with tracer studies in animals and retinotopy studies in humans. Additionally, previously undescribed topologies were revealed in the other pathways. This approach enhances the power of modern imaging by enabling study of fiber connections among anatomically and functionally defined brain regions in individual human subjects.

Quantitative diffusion-tensor anisotropy brain MR imaging: normative human data and anatomic analysis.

PURPOSE: To obtain normative human cerebral data and evaluate the anatomic information in quantitative diffusion anisotropy magnetic resonance (MR) imaging. MATERIALS AND METHODS: Quantitative diffusion anisotropy MR images were obtained in 13 healthy adults by using single-shot echo-planar MR imaging and a combination of tetrahedral and orthogonal gradient encoding (whole-brain coverage in about 1 minute). White matter (WM) anatomy was assessed at visual inspection, and values were measured in various brain regions. Different anisotropy measures, including total anisotropy (A sigma), were compared on the basis of information content, rotational invariance, and susceptibility to noise. Partial volume and noise effects were simulated. RESULTS: Anisotropy MR images depicted WM features not typically seen on conventional MR images (e.g., external capsule, thalamic substructures, basal ganglia, occipital WM, thickness of the internal capsule). Statistically significant anisotropy differences occurred across brain regions, which were reproducible within and across subjects. A sigma was highest in commissural WM and progressively lower in projection and association WM. This order paralleled that of known resistance to spread of vasogenic edema, which suggested that anisotropy may be sensitive to WM histologic structure. Gray matter (GM) A sigma data were consistent with zero anisotropy, and partial volume WM-GM effects were approximately linear. A sigma image quality could be effectively improved by means of averaging. CONCLUSION: Quantitative diffusion anisotropy images can be obtained rapidly and demonstrate subtle WM anatomy. Different histologic types of WM have significant and reproducible anisotropy differences.

Direct comparison of episodic encoding and retrieval of words: an event-related fMRI study.

Functional magnetic resonance imaging (fMRI) was used to compare directly episodic encoding and retrieval. During encoding, subjects studied visually presented words and reported via keypress whether each word represented a pleasant or unpleasant concept (intentional, deep encoding). During the retrieval phase, subjects indicated (via keypress) whether visually presented words had previously been studied. No reliable differences were found during the recognition phase for words that had been previously studied and those that had not been studied. Areas preferentially active during encoding (relative to retrieval) included left superior frontal cortex, medial frontal cortex, left superior temporal cortex, posterior cingulate, left parahippocampal gyrus, and left inferior frontal gyrus. Regions more active in retrieval than encoding included bilateral inferior parietal cortex, bilateral precuneus, right frontal polar cortex, right dorsolateral prefrontal cortex, and right inferior frontal/insular cortex.

A common network of functional areas for attention and eye movements.

Functional magnetic resonance imaging (fMRI) and surface-based representations of brain activity were used to compare the functional anatomy of two tasks, one involving covert shifts of attention to peripheral visual stimuli, the other involving both attentional and saccadic shifts to the same stimuli. Overlapping regional networks in parietal, frontal, and temporal lobes were active in both tasks. This anatomical overlap is consistent with the hypothesis that attentional and oculomotor processes are tightly integrated at the neural level.

Normal brain in human newborns: apparent diffusion coefficient and diffusion anisotropy measured by using diffusion tensor MR imaging.

PURPOSE: To establish quantitative standards for the directionally averaged water apparent diffusion coefficient (D) and quantitative diffusion anisotropy (A sigma) of normal brains in newborns by using diffusion tensor magnetic resonance (MR) imaging. MATERIALS AND METHODS: Diffusion tensor MR imaging was performed during the first 36 hours of life in 22 newborns (gestational age range, 31-41 weeks). Values of D and A sigma were measured in regions of interest chosen in the cortical gray matter, centrum semiovale, caudate nuclei, lentiform nuclei, thalami, internal capsules, and cerebellar hemispheres. RESULTS: The D values in the gray and white matter in newborns are considerably higher than those in adults. There is a striking correlation between gestational age and D, with D decreasing as gestational age increases. The A sigma values in the white matter in newborns are lower than those in adults. Values of A sigma show statistically significant correlations with gestational age only in the white matter of the centrum semiovale, in which A sigma values increase sharply near term. CONCLUSION: The D values primarily reflect overall brain water content. The A sigma values are more sensitive to tissue microstructure (e.g., white matter packing and myelination). The D and A sigma images reveal information and not apparent on T1- and T2-weighted images.

Functional MRI studies of word-stem completion: reliability across laboratories and comparison to blood flow imaging with PET.

Functional magnetic resonance imaging (fMRI) based on blood oxygen level-dependent (BOLD) contrast has become an increasingly popular technique for mapping the brain. The relationship between BOLD-fMRI imaging and imaging of blood flow activation with positron emission tomography (PET) remains unclear. Moreover, BOLD imaging strategies and analysis procedures vary widely across laboratories. To examine the relationship between these different methods, we compared brain activation maps of a word-stem completion task obtained both using PET and using fMRI across two separate institutions (Washington University and Massachusetts General Hospital) with different acquisitions (gradient-refocused echo and asymmetric spin echo) and different analysis techniques. Overall, activation maps were highly similar across both fMRI methods and PET. A set of activated brain areas, in consistent locations in Talairach atlas space, were identified across all three studies, including visual striate and extrastriate, left prefrontal, supplementary motor area (SMA), and right cerebellar areas. Decreases in activation were also consistently observed in medial parietal, posterior insular, and medial inferior frontal areas. Some differences were noted that may be related to the silent performance of the task with fMRI. The largely consistent results suggest that comparisons can be made appropriately across imaging modalities and laboratory methods. A further implication of the consistencies, which extended to both increases and decreases in signal, is that the underlying brain physiology leading to BOLD contrast may be more similar to blood flow than originally appreciated.

Hemispheric specialization in human dorsal frontal cortex and medial temporal lobe for verbal and nonverbal memory encoding.

The involvement of dorsal frontal and medial temporal regions during the encoding of words, namable line-drawn objects, and unfamiliar faces was examined using functional magnetic resonance imaging (fMRI). Robust dorsal frontal activations were observed in each instance, but lateralization was strongly dependent on the materials being encoded. Encoding of words produced left-lateralized dorsal frontal activation, whereas encoding of unfamiliar faces produced homologous right-lateralized activation. Encoding of namable objects, which are amenable to both verbal and nonverbal encoding, yielded bilateral dorsal frontal activation. A similar pattern of results was observed in the medial temporal lobe. These results indicate that regions in both hemispheres underlie human long-term memory encoding, and these regions can be engaged differentially according to the nature of the material being encoded.

Contrast-agent phase effects: an experimental system for analysis of susceptibility, concentration, and bolus input function kinetics.

A system is presented for experimental arterial input function (AIF) simulation and for accurate measurement of the concentration, susceptibility effects, and magnetic moment of paramagnetic MR contrast agents. Signal effects of contrast agents are evaluated with a stable, well-characterized, and precise experimental setup. A cylindrical phantom and a closed-loop circulating flow system were designed for AIF simulation, assessment of the physical determinants of contrast-agent phase effects, and measurement of contrast-agent properties under controlled conditions. A mathematical model of the AIF dynamics is proposed. From the experimental phase shift (delta phi), either the concentration or molar susceptibility, chiM, is determined. The linear dependence of delta phi on concentration and echo time (TE), the orientation dependence, and the lack of dependence on T1, T2, and diffusion time are proven precisely for water solutions under a wide variety of conditions. The measured effective magnetic moment of Gd+3, mu(eff), was 7.924 +/- 0.015 Bohr magnetons in agreement with the theoretical value of 7.937.

Anatomic localization and quantitative analysis of gradient refocused echo-planar fMRI susceptibility artifacts.

Functional magnetic resonance imaging (fMRI) techniques, such as echo-planar imaging, can permit rapid, sensitive, whole-brain measurements of local blood flow-induced MR signal changes seen during cognitive paradigms. Changes in blood oxygenation due to mismatch of flow and oxygen metabolism cause dynamic variations in microscopic susceptibility effects, leading to the blood oxygenation level-dependent (BOLD) signal measured by fMRI techniques. A related static macroscopic susceptibility effect is known to cause artifacts that attenuate the MR signal, leading to "blind spots" in some regions of brain adjacent to bone and air sinuses. The anatomical location, spatial extent, and magnitude of signal loss artifact are quantitated for a common whole-brain fMRI technique. Resting gradient-echo EPI studies were obtained in four healthy volunteers. Signal loss was primarily localized to inferior frontal regions (medial orbital gyri and gyrus rectus) and to inferior lateral temporal lobe (including part of fusiform gyrus) bilaterally. Increased echo time (TE) uniformly produced larger artifacts. The orientation of acquired slices and choice of phase-encoding direction influenced the location, shape, and extent of the artifacts. Regions of the brain with severe artifact may have attenuated activation signal, with potential implications for the design and interpretation of fMRI studies targeting activations in these areas.

Arterial input functions from MR phase imaging.

A method is presented for obtaining high-sensitivity arterial input functions following bolus intravenous contrast agent administration. Arterial contrast agent is monitored by phase reconstruction of single-shot echo-planar images. During bolus injections of a gadolinium (Gd) agent in a baboon, data were acquired at the mid-abdominal aorta, and magnitude and phase-shift images were reconstructed. Pairwise image subtraction was used to minimize phase aliasing. The phase-based method is shown to have a significant potential improvement in sensitivity compared to the magnitude approach. The phase method also has a general linear response to concentration. This method may have potential utility in quantitative imaging of blood flow and contrast agent kinetics.

Encoding of anisotropic diffusion with tetrahedral gradients: a general mathematical diffusion formalism and experimental results.

A diffusion imaging method with a tetrahedral sampling pattern has been developed for high-sensitivity diffusion analysis. The tetrahedral gradient pattern consists of four different combinations of x, y, and z gradients applied simultaneously at full strength to uniformly measure diffusion in four different directions. Signal-to-noise can be increased by up to a factor of about three using this approach, compared with diffusion measurements made using separately applied x, y, and z gradients. A mathematical formalism is presented describing six fundamental parameters: the directionally averaged diffusion coefficient D and diffusion element anisotropies eta and epsilon which are rotationally invariant, and diffusion ellipsoid orientation angles theta, phi, and omega which are rotationally variant. These six parameters contain all the information in the symmetric diffusion tensor D. Principal diffusion coefficients, reduced anisotropies, and other rotational invariants are further defined. It is shown that measurement of off-diagonal tensor elements is essential to assess anisotropy and orientation, and that the only parameter which can be measured with the orthogonal method is D. In cases of axial diffusion symmetry (e.g., fibers), the four tetrahedral diffusion measurements efficiently enable determination of D, eta, theta, and phi which contain all the diffusion information. From these four parameters, the diffusion parallel and perpendicular to the symmetry axis (D and D) and the axial anisotropy A can be determined. In more general cases, the six fundamental parameters can be determined with two additional diffusion measurements. Tetrahedral diffusion sequences were implemented on a clinical MR system. A muscle phantom demonstrates orientation independence of D, D, D, and A for large changes in orientation angles. Sample background gradients and diffusion gradient imbalances were directly measured and found to be insignificant in most cases.

Diffusion MRI: precision, accuracy and flow effects.

After a decade of evolution and application of diffusion imaging, a large body of literature has been accumulated. It is in this context that the accuracy and precision of diffusion-weighted and quantitative diffusion MRI are reviewed. The emphasis of the review is on practical methods for clinical human imaging, particularly in the brain. The requirements for accuracy and precision are reviewed for various clinical and basic science applications. The methods of measuring and calculating diffusion effects with MRI are reviewed. The pulse gradient spin echo (PGSE) methods are emphasized as these methods are used most commonly in the clinical setting. Processing of PGSE data is reviewed. Various PGSE encoding schemes are also reviewed in terms of the accuracy and precision of isotropic and anisotropic diffusion measurements. The broad range of factors impacting the accuracy of the PGSE methods and other encoding schemes is then considered. Firstly, system inaccuracies such as background imaging gradients, gradient linearity, refocusing RF pulses, eddy currents, image misregistration, noise and dynamic range are considered. A second class of inaccuracies is contributed by the bulk effects of the imaged object, and include sample background gradients, subject motion of cerebrospinal fluid and organs, and aperiodic organ motion. A final category of potential inaccuracies is classified as being contributed by microscopic, biophysical tissue properties and include partial volume effects, anisotropy, restriction, diffusion distance, compartmentation, exchange, multiexponential diffusion decay, T2 weighting and microvascular perfusion. Finally, the application of diffusion methods to studies of blood flow in the microvasculature (i.e. the arterioles, capillaries and venules) are reviewed in detail, particularly in terms of feasibility and the stringent accuracy and precision requirements. Recent provocative studies examining the use of PGSE approaches to suppress microvascular signals in brain functional MRI (fMRI) are also reviewed.