Publications

We propose a technique to simultaneously estimate the local fiber orientations and perform multi-fiber tractography. Existing techniques estimate the local fiber orientation at each voxel independently so there is no running knowledge of confidence in the measured signal or estimated fiber orientation. Further, to overcome noise, many algorithms use a filter as a post-processing step to obtain a smooth trajectory. We formulate fiber tracking as causal estimation: at each step of tracing the fiber, the current estimate of the signal is guided by the previous. To do this, we model the signal as a discrete mixture of Watson directional functions and perform tractography within a filtering framework. Starting from a seed point, each fiber is traced to its termination using an unscented Kalman filter to simultaneously fit the signal and propagate in the most consistent direction. Despite the presence of noise and uncertainty, this provides an accurate estimate of the local structure at each point along the fiber. We choose the Watson function since it provides a compact representation of the signal parameterized by the principal diffusion direction and a scaling parameter describing anisotropy, and also allows analytic reconstruction of the oriented diffusion function from those parameters. Using a mixture of two and three components (corresponding to two-fiber and three-fiber models) we demonstrate in synthetic experiments that this approach reduces signal reconstruction error and significantly improves the angular resolution at crossings and branchings. In vivo experiments examine the corpus callosum and internal capsule and confirm the ability to trace through regions known to contain such crossing and branching while providing inherent path regularization.

Water diffusion in nerve fibers is strongly influenced by axon architecture. In this study, fractional diffusion anisotropy and transverse and longitudinal diffusion coefficients were measured in excised human cervical spinal cord with MR line-scan diffusion imaging, at 625 microm in-plane resolution and 3 mm slice thickness. A pixel-based comparison of fractional diffusion anisotropy, transverse diffusion coefficient, and longitudinal diffusion coefficient data with axon packing parameters derived from corresponding stained histological sections was performed for four slices. The axon packing parameters, axon density, axon area-fraction, and average axon size for entire specimen cross-sections were calculated by computerized segmentation of optical microscopy data obtained at 0.53 microm resolution. Salient features could be recognized on fractional diffusion anisotropy, transverse diffusion coefficient, axon density, axon area fraction, and average axon size maps. For white matter regions only, the average correlation coefficients for fractional diffusion anisotropy compared to histology-based parameters axon density and axon area fraction were 0.37 and 0.21, respectively. For transverse diffusion coefficient compared to axon density and axon area fraction, they were -0.40 and -0.36, and for longitudinal diffusion coefficient compared to axon density and axon area fraction, -0.14 and -0.30. All average correlation coefficients for average axon size were low. Correlation coefficients for collectively analyzed white and gray matter regions were significantly higher than correlation coefficients derived from analysis of white matter regions only.

Nyquist ghosts are an inherent artifact in echo planar imaging acquisitions. An approach to robustly eliminate Nyquist ghosts is presented that integrates two previous Nyquist ghost correction techniques: temporal domain encoding (phase labeling for additional coordinate encoding: PLACE and spatial domain encoding (phased array ghost elimination: PAGE). Temporal encoding modulates the echo planar imaging acquisition trajectory from frame to frame, enabling one to interleave data to remove inconsistencies that occur between sampling on positive and negative gradient readouts. With PLACE, one can coherently combine the interleaved data to cancel residual Nyquist ghosts. If the level of ghosting varies significantly from image to image, however, the signal cancellation that occurs with PLACE can adversely affect SNR-sensitive applications such as perfusion imaging with arterial spin labeling. This work proposes integrating PLACE into a PAGE-based reconstruction process to yield significantly better Nyquist ghost correction that is more robust than PLACE or PAGE alone. The robustness of this method is demonstrated in the presence of magnetic field drift with an in-vivo arterial spin labeling perfusion experiment.

The quantification of brain asymmetries may provide biomarkers for presurgical localization of language function and can improve our understanding of neural structure-function relationships in health and disease. We propose a new method for studying the asymmetry of the white matter tracts in the entire brain, and we apply it to a preliminary study of normal subjects across the handedness spectrum. Methods for quantifying white matter asymmetry using diffusion MRI tractography have thus far been based on comparing numbers of fibers or volumes of a single fiber tract across hemispheres. We propose a generalization of such methods, where the "number of fibers" laterality measurement is extended to the entire brain using a soft fiber comparison metric. We summarize the distribution of fiber laterality indices over the whole brain in a histogram, and we measure properties of the distribution such as its skewness, median, and inter-quartile range. The whole-brain fiber laterality histogram can be measured in an exploratory fashion without hypothesizing asymmetries only in particular structures. We demonstrate an overall difference in white matter asymmetry in consistent- and inconsistent-handers: the skewness of the fiber laterality histogram is significantly different across handedness groups.

As the number and complexity of partially sampled dynamic imaging methods continue to increase, reliable strategies to evaluate performance may prove most useful. In the present work, an analytical framework to evaluate given reconstruction methods is presented. A perturbation algorithm allows the proposed evaluation scheme to perform robustly without requiring knowledge about the inner workings of the method being evaluated. A main output of the evaluation process consists of a two-dimensional modulation transfer function, an easy-to-interpret visual rendering of a method’s ability to capture all combinations of spatial and temporal frequencies. Approaches to evaluate noise properties and artifact content at all spatial and temporal frequencies are also proposed. One fully sampled phantom and three fully sampled cardiac cine datasets were subsampled (R = 4 and 8) and reconstructed with the different methods tested here. A hybrid method, which combines the main advantageous features observed in our assessments, was proposed and tested in a cardiac cine application, with acceleration factors of 3.5 and 6.3 (skip factors of 4 and 8, respectively). This approach combines features from methods such as k-t sensitivity encoding, unaliasing by Fourier encoding the overlaps in the temporal dimension-sensitivity encoding, generalized autocalibrating partially parallel acquisition, sensitivity profiles from an array of coils for encoding and reconstruction in parallel, self, hybrid referencing with unaliasing by Fourier encoding the overlaps in the temporal dimension and generalized autocalibrating partially parallel acquisition, and generalized autocalibrating partially parallel acquisition-enhanced sensitivity maps for sensitivity encoding reconstructions.

The pulsed-field gradient (PFG) MR experiment enables one to measure particle displacements, velocities, and even higher moments of complex fluid motions. In diffusion-weighted MRI (DWI) in living tissue, where the PFG MRI experiment is used to measure diffusion, Brownian motion is assumed to dominate the displacements causing the observed signal loss. However, motions of water molecules caused by various active biological processes occurring at different length and time scales may also cause additional dephasing of magnetization and signal loss. To help understand their relative effects on the DWI signal attenuation, we used an integrated experimental and theoretical framework: a Rheo-NMR, which served as an experimental model system to precisely prescribe a microscopic velocity distribution; and a mathematical model that relates the DW signal intensity in the Rheo-NMR to experimental parameters that characterize the impressed velocity field. A technical innovation reported here is our use of ’natural’ (in this case, polar) coordinates both to simplify the description the fluid motion within the Couette cell of the Rheo-NMR, as well as to acquire and reconstruct magnitude and phase MR images obtained within it. We use this integrated model system to demonstrate how shear flows appears as pseudo-diffusion in magnitude DW MR signals obtained using PFG spin-echo (PGSE) NMR and MRI sequences. Our results lead us to reinterpret the possible causes of signal loss in DWI in vivo, in particular to revise and generalize the previous notion of intra-voxel incoherent motion (IVIM) in order to describe activity driven flows that appear as pseudo-diffusion over multiple length and time scales in living tissues.

The present study examined the relationship between hand preference degree and direction, functional language lateralization in Broca’s and Wernicke’s areas, and structural measures of the arcuate fasciculus. Results revealed an effect of degree of hand preference on arcuate fasciculus structure, such that consistently-handed individuals, regardless of the direction of hand preference, demonstrated the most asymmetric arcuate fasciculus, with larger left versus right arcuate, as measured by DTI. Functional language lateralization in Wernicke’s area, measured via fMRI, was related to arcuate fasciculus volume in consistent-left-handers only, and only in people who were not right hemisphere lateralized for language; given the small sample size for this finding, future investigation is warranted. Results suggest handedness degree may be an important variable to investigate in the context of neuroanatomical asymmetries.

BACKGROUND: While the neuroanatomical underpinnings of the functional brain disconnectivity observed in patients with schizophrenia (SZ) remain elusive, white matter fiber bundles of the brain are a likely candidate, given that they represent the infrastructure for long-distance neural communication. METHODS: This study investigated for diffusion abnormalities in 19 patients with chronic SZ, relative to 19 matched control subjects, across tractography-defined segments of the corpus callosum. Diffusion-weighted images were acquired with 51 noncollinear gradients on a 3T scanner (1.7 mm isotropic voxels). The corpus callosum was extracted by means of whole-brain tractography and automated fiber clustering and was parcelled into six segments on the basis of fiber trajectories. The diffusion indexes of fractional anisotropy (FA) and mode were calculated for each segment. RESULTS: Relative to the healthy control subjects, the SZ patients exhibited mode increases in the parietal fibers, suggesting a relative absence of crossing fibers. Schizophrenia patients also exhibited FA reductions in the frontal fibers, which were underpinned by increases in radial diffusivity, consistent with myelin abnormalities. Significant correlations were observed between patients’ degree of reality distortion and their FA and radial diffusivity, such that the most severely psychotic patients were the least abnormal in terms of their frontal fiber diffusivity. CONCLUSIONS: The SZ patients exhibited a variety of diffusion abnormalities in the corpus callosum, which were related to the severity of their psychotic symptoms. To the extent that diffusion abnormalities influence axonal transmission velocities, these results provide support for those theories that emphasize neural timing abnormalities in the etiology of schizophrenia.

Aging is associated with a variety of structural and pathological brain changes. Memory, or the ability to store and retrieve information, declines significantly during aging. In order to characterize the brain micro-structural correlates of memory performance, 52 healthy subjects, age 25-82 years, completed a computerized non-verbal memory test and were scanned using magnetic resonance diffusion tensor imaging. Partial correlation was conducted between DTI parameters and memory performance (accuracy and reaction time (RT) for different learning stages) controlling for age. A similar correlation pattern was found for apparent diffusion coefficient (ADC), FA, and radial and axial diffusivities, but correlations between ADC and memory performance were the most informative and are therefore reported here. While ADC was correlated with accuracy mainly in temporal and frontal cortical regions, it was correlated with RT in temporal and frontal white matter pathways, including the inferior longitudinal fasciculus and uncinate fasciculus. The task was repeated four times, performance in the first repetition was correlated with ADC in frontal white matter and in the fourth repetition with ADC in temporal gray matter structures mainly the parahippocapus and in the middle temporal gyrus. The localization of the correlations of ADC with the different task parameters is in line with previous studies. Thus, inter-subject variability in memory performance and tissue morphology, as expressed by diffusion tensor magnetic resonance imaging, can be used to relate tissue microstructures with cognitive performance, and to provide information to corroborate other functional localization techniques.

NMR observable nuclei undergoing restricted diffusion within confining pores are important reporters for microstructural features of porous media including, inter-alia, biological tissues, emulsions and rocks. Diffusion NMR, and especially the single-pulsed field gradient (s-PFG) methodology, is one of the most important noninvasive tools for studying such opaque samples, enabling extraction of important microstructural information from diffusion-diffraction phenomena. However, when the pores are not monodisperse and are characterized by a size distribution, the diffusion-diffraction patterns disappear from the signal decay, and the relevant microstructural information is mostly lost. A recent theoretical study predicted that the diffusion-diffraction patterns in double-PFG (d-PFG) experiments have unique characteristics, such as zero-crossings, that make them more robust with respect to size distributions. In this study, we theoretically compared the signal decay arising from diffusion in isolated cylindrical pores characterized by lognormal size distributions in both s-PFG and d-PFG methodologies using a recently presented general framework for treating diffusion in NMR experiments. We showed the gradual loss of diffusion-diffraction patterns in broadening size distributions in s-PFG and the robustness of the zero-crossings in d-PFG even for very large standard deviations of the size distribution. We then performed s-PFG and d-PFG experiments on well-controlled size distribution phantoms in which the ground-truth is well-known a priori. We showed that the microstructural information, as manifested in the diffusion-diffraction patterns, is lost in the s-PFG experiments, whereas in d-PFG experiments the zero-crossings of the signal persist from which relevant microstructural information can be extracted. This study provides a proof of concept that d-PFG may be useful in obtaining important microstructural features in samples characterized by size distributions.