PURPOSE: Diffusion-weighted MRI is sensitive to incoherent tissue motion, which may confound the measured signal and subsequent analysis. We propose a "motion-compensated" gradient waveform design for tensor-valued diffusion encoding that negates the effects bulk motion and incoherent motion in the ballistic regime. METHODS: Motion compensation was achieved by constraining the magnitude of gradient waveform moment vectors. The constraint was incorporated into a numerical optimization framework, along with existing constraints that account for b-tensor shape, hardware restrictions, and concomitant field gradients. We evaluated the efficacy of encoding and motion compensation in simulations, and we demonstrated the approach by linear and planar b-tensor encoding in a healthy heart in vivo. RESULTS: The optimization framework produced asymmetric motion-compensated waveforms that yielded b-tensors of arbitrary shape with improved efficiency compared with previous designs for tensor-valued encoding, and equivalent efficiency to previous designs for linear (conventional) encoding. Technical feasibility was demonstrated in the heart in vivo, showing vastly improved data quality when using motion compensation. The optimization framework is available online in open source. CONCLUSION: Our gradient waveform design is both more flexible and efficient than previous methods, facilitating tensor-valued diffusion encoding in tissues in which motion would otherwise confound the signal. The proposed design exploits asymmetric encoding times, a single refocusing pulse or multiple refocusing pulses, and integrates compensation for concomitant gradient effects throughout the imaging volume.
Multidimensional MRI Core Publications
PURPOSE: Reperfusion therapy enables effective treatment of ischemic stroke presenting within 4-6 hours. However, tissue progression from ischemia to infarction is variable, and some patients benefit from treatment up until 24 hours. Improved imaging techniques are needed to identify these patients. Here, it was hypothesized that time dependence in diffusion MRI may predict tissue outcome in ischemic stroke. METHODS: Diffusion MRI data were acquired with multiple diffusion times in five non-reperfused patients at 2, 9, and 100 days after stroke onset. Maps of "rate of kurtosis change" (k), mean kurtosis, ADC, and fractional anisotropy were derived. The ADC maps defined lesions, normal-appearing tissue, and the lesion tissue that would either be infarcted or remain viable by day 100. Diffusion parameters were compared (1) between lesions and normal-appearing tissue, and (2) between lesion tissue that would be infarcted or remain viable. RESULTS: Positive values of k were observed within stroke lesions on day 2 (P = .001) and on day 9 (P = .023), indicating diffusional exchange. On day 100, high ADC values indicated infarction of 50 ± 20% of the lesion volumes. Tissue infarction was predicted by high k values both on day 2 (P = .026) and on day 9 (P = .046), by low mean kurtosis values on day 2 (P = .043), and by low fractional anisotropy values on day 9 (P = .029), but not by low ADC values. CONCLUSIONS: Diffusion time dependence predicted tissue outcome in ischemic stroke more accurately than the ADC, and may be useful for predicting reperfusion benefit.
Probing the cellular structure of in vivo biological tissue is a fundamental problem in biomedical imaging and medical science. This work introduces an approach for analyzing diffusion magnetic resonance imaging data acquired by the novel tensor-valued encoding technique for characterizing tissue microstructure. Our approach first uses a signal model to estimate the variance and skewness of the distribution of apparent diffusion tensors modeling the underlying tissue. Then several novel imaging indices, such as weighted microscopic anisotropy and microscopic skewness, are derived to characterize different ensembles of diffusion processes that are indistinguishable by existing techniques. The contributions of this work also include a theoretical proof that shows that, to estimate the skewness of a diffusion tensor distribution, the encoding protocol needs to include full-rank tensor diffusion encoding. This proof provides a guideline for the application of this technique. The properties of the proposed indices are illustrated using both synthetic data and in vivo data acquired from a human brain.
Language and theory of mind (ToM) are the cognitive capacities that allow for the successful interpretation and expression of meaning. While functional MRI investigations are able to consistently localize language and ToM to specific cortical regions, diffusion MRI investigations point to an inconsistent and sometimes overlapping set of white matter tracts associated with these two cognitive domains. To further examine the white matter tracts that may underlie these domains, we use a two-tensor tractography method to investigate the white matter microstructure of 809 participants from the Human Connectome Project. 20 association white matter tracts (10 in each hemisphere) are uniquely identified by leveraging a neuroanatomist-curated automated white matter tract atlas. The mean fractional anisotropy (FA), mean diffusivity (MD), and number of streamlines (NoS) are measured for each white matter tract. Performance on neuropsychological assessments of semantic memory (NIH Toolbox Picture Vocabulary Test, TPVT) and emotion perception (Penn Emotion Recognition Test, PERT) are used to measure critical subcomponents of the language and ToM networks, respectively. Regression models are constructed to examine how structural measurements of left and right white matter tracts influence performance across these two assessments. We find that semantic memory performance is influenced by the number of streamlines of the left superior longitudinal fasciculus III (SLF-III), and emotion perception performance is influenced by the number of streamlines of the right SLF-III. Additionally, we find that performance on both semantic memory & emotion perception is influenced by the FA of the left arcuate fasciculus (AF). The results point to multiple, overlapping white matter tracts that underlie the cognitive domains of language and ToM. Results are discussed in terms of hemispheric dominance and concordance with prior investigations.
BACKGROUND: Preoperative radiological assessment of meningioma characteristics is of value for pre- and post-operative patient management, counselling, and surgical approach. PURPOSE: To investigate whether tensor-valued diffusion MRI can add to the preoperative prediction of meningioma consistency, grade and type. MATERIALS AND METHODS: 30 patients with intracranial meningiomas (22 WHO grade I, 8 WHO grade II) underwent MRI prior to surgery. Diffusion MRI was performed with linear and spherical b-tensors with b-values up to 2000 s/mm2. The data were used to estimate mean diffusivity (MD), fractional anisotropy (FA), mean kurtosis (MK) and its components-the anisotropic and isotropic kurtoses (MKA and MKI). Meningioma consistency was estimated for 16 patients during resection based on ultrasonic aspiration intensity, ease of resection with instrumentation or suction. Grade and type were determined by histopathological analysis. The relation between consistency, grade and type and dMRI parameters was analyzed inside the tumor ("whole-tumor") and within brain tissue in the immediate periphery outside the tumor ("rim") by histogram analysis. RESULTS: Lower 10th percentiles of MK and MKA in the whole-tumor were associated with firm consistency compared with pooled soft and variable consistency (n = 7 vs 9; U test, p = 0.02 for MKA 10 and p = 0.04 for MK10) and lower 10th percentile of MD with variable against soft and firm (n = 5 vs 11; U test, p = 0.02). Higher standard deviation of MKI in the rim was associated with lower grade (n = 22 vs 8; U test, p = 0.04) and in the MKI maps we observed elevated rim-like structure that could be associated with grade. Higher median MKA and lower median MKI distinguished psammomatous type from other pooled meningioma types (n = 5 vs 25; U test; p = 0.03 for MKA 50 and p = 0.03 and p = 0.04 for MKI 50). CONCLUSION: Parameters from tensor-valued dMRI can facilitate prediction of consistency, grade and type.
PURPOSE: To accelerate the acquisition of relaxation-diffusion imaging by integrating time-division multiplexing (TDM) with simultaneous multi-slice (SMS) for EPI and evaluate imaging quality and diffusion measures. METHODS: The time-division multiplexing (TDM) technique and SMS method were integrated to achieve a high slice-acceleration (e.g., 6×) factor for acquiring relaxation-diffusion MRI. Two variants of the sequence, referred to as TDM3e-SMS and TDM2s-SMS, were developed to simultaneously acquire slice groups with three distinct TEs and two slice groups with the same TE, respectively. Both sequences were evaluated on a 3T scanner with in vivo human brains and compared with standard single-band (SB) -EPI and SMS-EPI using diffusion measures and tractography results. RESULTS: Experimental results showed that the TDM3e-SMS sequence with total slice acceleration of 6 (multiplexing factor (MP) = 3 × multi-band factor (MB) = 2) provided similar image intensity and microstructure measures compared to standard SMS-EPI with MB = 2, and yielded less bias in intensity compared to standard SMS-EPI with MB = 4. The three sequences showed a similar positive correlation between TE and mean kurtosis (MK) and a negative correlation between TE and mean diffusivity (MD) in white matter. Multi-fiber tractography also shows consistency of results in TE-dependent measures between different sequences. The TDM2s-SMS sequence (MP = 2, MB = 2) also provided imaging measures similar to standard SMS-EPI sequences (MB = 2) for single-TE diffusion imaging. CONCLUSIONS: The TDM-SMS sequence can provide additional 2x to 3x acceleration to SMS without degrading imaging quality. With the significant reduction in scan time, TDM-SMS makes joint relaxation-diffusion MRI a feasible technique in neuroimaging research to investigate new markers of brain disorders.
Diffusion magnetic resonance imaging (dMRI) tractography is an advanced imaging technique that enables in vivo reconstruction of the brain’s white matter connections at macro scale. It provides an important tool for quantitative mapping of the brain’s structural connectivity using measures of connectivity or tissue microstructure. Over the last two decades, the study of brain connectivity using dMRI tractography has played a prominent role in the neuroimaging research landscape. In this paper, we provide a high-level overview of how tractography is used to enable quantitative analysis of the brain’s structural connectivity in health and disease. We focus on two types of quantitative analyses of tractography, including: 1) tract-specific analysis that refers to research that is typically hypothesis-driven and studies particular anatomical fiber tracts, and 2) connectome-based analysis that refers to research that is more data-driven and generally studies the structural connectivity of the entire brain. We first provide a review of methodology involved in three main processing steps that are common across most approaches for quantitative analysis of tractography, including methods for tractography correction, segmentation and quantification. For each step, we aim to describe methodological choices, their popularity, and potential pros and cons. We then review studies that have used quantitative tractography approaches to study the brain’s white matter, focusing on applications in neurodevelopment, aging, neurological disorders, mental disorders, and neurosurgery. We conclude that, while there have been considerable advancements in methodological technologies and breadth of applications, there nevertheless remains no consensus about the "best" methodology in quantitative analysis of tractography, and researchers should remain cautious when interpreting results in research and clinical applications.
Matrix metalloproteinases 9 (MMP9) are enzymes involved in regulating neuroplasticity in the hippocampus. This, combined with evidence for disrupted hippocampal structure and function in schizophrenia, has prompted our current investigation into the relationship between MMP9 and hippocampal volumes in schizophrenia. 34 healthy individuals (mean age = 32.50, male = 21, female = 13) and 30 subjects with schizophrenia (mean age = 33.07, male = 19, female = 11) underwent a blood draw and T1-weighted magnetic resonance imaging. The hippocampus was automatically segmented utilizing FreeSurfer. MMP9 plasma levels were measured with ELISA. ANCOVAs were conducted to compare MMP9 plasma levels (corrected for age and sex) and hippocampal volumes between groups (corrected for age, sex, total intracranial volume). Spearman correlations were utilized to investigate the relationship between symptoms, medication, duration of illness, number of episodes, and MMP9 plasma levels in patients. Last, we explored the correlation between MMP9 levels and hippocampal volumes in patients and healthy individuals separately. Patients displayed higher MMP9 plasma levels than healthy individuals (F(1, 60) = 21.19, p < 0.0001). MMP9 levels correlated with negative symptoms in patients (R = 0.39, p = 0.035), but not with medication, duration of illness, or the number of episodes. Further, patients had smaller left (F(1,59) = 9.12, p = 0.0040) and right (F(1,59) = 6.49, p = 0.013) hippocampal volumes. Finally, left (R = -0.39, p = 0.034) and right (R = -0.37, p = 0.046) hippocampal volumes correlated negatively with MMP9 plasma levels in patients. We observe higher MMP9 plasma levels in SCZ, associated with lower hippocampal volumes, suggesting involvement of MMP9 in the pathology of SCZ. Future studies are needed to investigate how MMP9 influences the pathology of SCZ over the lifespan, whether the observed associations are specific for schizophrenia, and if a therapeutic modulation of MMP9 promotes neuroprotective effects in SCZ.
Monitoring time dependence with diffusion MRI yields observables sensitive to compartment sizes (restricted diffusion) and membrane permeability (water exchange). However, restricted diffusion and exchange have opposite effects on the diffusion-weighted signal, which can lead to errors in parameter estimates. In this work, we propose a signal representation that incorporates the effects of both restricted diffusion and exchange up to second order in b-value and is compatible with gradient waveforms of arbitrary shape. The representation features mappings from a gradient waveform to two scalars that separately control the sensitivity to restriction and exchange. We demonstrate that these scalars span a two-dimensional space that can be used to choose waveforms that selectively probe restricted diffusion or exchange, eliminating the correlation between the two phenomena. We found that waveforms with specific but unconventional shapes provide an advantage over conventional pulsed and oscillating gradient acquisitions. We also show that parametrization of waveforms into a two-dimensional space can be used to understand protocols from other approaches that probe restricted diffusion and exchange. For example, we found that the variation of mixing time in filter-exchange imaging corresponds to variation of our exchange-weighting scalar at a fixed value of the restriction-weighting scalar. The proposed signal representation was evaluated using Monte Carlo simulations in identical parallel cylinders with hexagonal and random packing as well as parallel cylinders with gamma-distributed radii. Results showed that the approach is sensitive to sizes in the interval 4-12 μm $$ \upmu \mathrmm $$ and exchange rates in the simulated range of 0 to 20 s - 1 $$ \mathrms^-1 $$ , but also that there is a sensitivity to the extracellular geometry. The presented theory constitutes a simple and intuitive description of how restricted diffusion and exchange influence the signal as well as a guide to protocol design capable of separating the two effects.
OBJECTIVES: Approximately 50% of comatose patients after cardiac arrest never regain consciousness. Cerebral ischaemia may lead to cytotoxic and/or vasogenic oedema, which can be detected by diffusion tensor imaging (DTI). Here, we evaluate the potential value of free water corrected mean diffusivity (MD) and fractional anisotropy (FA) based on DTI, for the prediction of neurological recovery of comatose patients after cardiac arrest. METHODS: A total of 50 patients after cardiac arrest were included in this prospective cohort study in two Dutch hospitals. DTI was obtained 2-4 days after cardiac arrest. Outcome was assessed at 6 months, dichotomised as poor (cerebral performance category 3-5; n = 20) or good (n = 30) neurological outcome. We calculated the whole brain mean MD and FA and compared between patients with good and poor outcomes. In addition, we compared a preliminary prediction model based on clinical parameters with or without the addition of MD and FA. RESULTS: We found significant differences between patients with good and poor outcome of mean MD (good: 726 [702-740] × 10-6 mm2/s vs. poor: 663 [575-736] × 10-6 mm2/s; p = 0.01) and mean FA (0.30 ± 0.03 vs. 0.28 ± 0.03; p = 0.03). An exploratory prediction model combining clinical parameters, MD and FA increased the sensitivity for reliable prediction of poor outcome from 60 to 85%, compared to the model containing clinical parameters only, but confidence intervals are overlapping. CONCLUSIONS: Free water-corrected MD and FA discriminate between patients with good and poor outcomes after cardiac arrest and hold the potential to add to multimodal outcome prediction. KEY POINTS: • Whole brain mean MD and FA differ between patients with good and poor outcome after cardiac arrest. • Free water-corrected MD can better discriminate between patients with good and poor outcome than uncorrected MD. • A combination of free water-corrected MD (sensitive to grey matter abnormalities) and FA (sensitive to white matter abnormalities) holds potential to add to the prediction of outcome.