In vivo diffusion MRI of the human heart using a 300 mT/m gradient system.

PURPOSE: This work reports for the first time on the implementation and application of cardiac diffusion-weighted MRI on a Connectom MR scanner with a maximum gradient strength of 300 mT/m. It evaluates the benefits of the increased gradient performance for the investigation of the myocardial microstructure. METHODS: Cardiac diffusion-weighted imaging (DWI) experiments were performed on 10 healthy volunteers using a spin-echo sequence with up to second- and third-order motion compensation ( M 2 $$ {M}_2 $$ and M 3 $$ {M}_3 $$ ) and b = 100 , 450 $$ b=100,450 $$ , and 1000 s / m m 2 $$ \mathrm{s}/\mathrm{m}{\mathrm{m}}^2 $$ (twice the b max $$ {b}_{\mathrm{max}} $$ commonly used on clinical scanners). Mean diffusivity (MD), fractional anisotropy (FA), helix angle (HA), and secondary eigenvector angle (E2A) were calculated for b = [100, 450] s / m m 2 $$ \mathrm{s}/\mathrm{m}{\mathrm{m}}^2 $$ and b = [100, 1000] s / m m 2 $$ \mathrm{s}/\mathrm{m}{\mathrm{m}}^2 $$ for both M 2 $$ {M}_2 $$ and M 3 $$ {M}_3 $$ . RESULTS: The MD values with M 3 $$ {M}_3 $$ are slightly higher than with M 2 $$ {M}_2 $$ with Δ MD = 0 . 05 ± 0 . 05 [ × 1 0 - 3 mm 2 / s ] ( p = 4 e - 5 ) $$ \Delta \mathrm{MD}=0.05\pm 0.05\kern0.3em \left[\times 1{0}^{-3}\kern0.3em {\mathrm{mm}}^2/\mathrm{s}\right]\kern0.3em \left(p=4e-5\right) $$ for b max = 450 s / mm 2 $$ {b}_{\mathrm{max}}=450\kern0.3em \mathrm{s}/{\mathrm{mm}}^2 $$ and Δ MD = 0 . 03 ± 0 . 03 [ × 1 0 - 3 mm 2 / s ] ( p = 4 e - 4 ) $$ \Delta \mathrm{MD}=0.03\pm 0.03\kern0.3em \left[\times \kern0.3em 1{0}^{-3}\kern0.3em {\mathrm{mm}}^2/\mathrm{s}\right]\kern0.3em \left(p=4e-4\right) $$ for b max = 1000 s / mm 2 $$ {b}_{\mathrm{max}}=1000\kern0.3em \mathrm{s}/{\mathrm{mm}}^2 $$ . A reduction in MD is observed by increasing the b max $$ {b}_{\mathrm{max}} $$ from 450 to 1000 s / mm 2 $$ \mathrm{s}/{\mathrm{mm}}^2 $$ ( Δ MD = 0 . 06 ± 0 . 04 [ × 1 0 - 3 mm 2 / s ] ( p = 1 . 6 e - 9 ) $$ \Delta \mathrm{MD}=0.06\pm 0.04\kern0.3em \left[\times \kern0.3em 1{0}^{-3}\kern0.3em {\mathrm{mm}}^2/\mathrm{s}\right]\kern0.3em \left(p=1.6e-9\right) $$ for M 2 $$ {M}_2 $$ and Δ MD = 0 . 08 ± 0 . 05 [ × 1 0 - 3 mm 2 / s ] ( p = 1 e - 9 ) $$ \Delta \mathrm{MD}=0.08\pm 0.05\kern0.3em \left[\times \kern0.3em 1{0}^{-3}\kern0.3em {\mathrm{mm}}^2/\mathrm{s}\right]\kern0.3em \left(p=1e-9\right) $$ for M 3 $$ {M}_3 $$ ). The difference between FA, E2A, and HA was not significant in different schemes ( p > 0 . 05 $$ p>0.05 $$ ). CONCLUSION: This work demonstrates cardiac DWI in vivo with higher b-value and higher order of motion compensated diffusion gradient waveforms than is commonly used. Increasing the motion compensation order from M 2 $$ {M}_2 $$ to M 3 $$ {M}_3 $$ and the maximum b-value from 450 to 1000 s / mm 2 $$ \mathrm{s}/{\mathrm{mm}}^2 $$ affected the MD values but FA and the angular metrics (HA and E2A) remained unchanged. Our work paves the way for cardiac DWI on the next-generation MR scanners with high-performance gradient systems.

Compartmental anisotropy of filtered exchange imaging (FEXI) in human white matter: What is happening in FEXI?

PURPOSE: To investigate the effects of compartmental anisotropy on filtered exchange imaging (FEXI) in white matter (WM). THEORY AND METHODS: FEXI signals were measured using multiple combinations of diffusion filter and detection directions in five healthy volunteers. Additional filters, including a trace-weighted diffusion filter with trapezoidal gradients, a spherical b-tensor encoded diffusion filter, and a T2 filter, were tested with trace-weighted diffusion detection. RESULTS: A large range of apparent exchange rates (AXR) and both positive and negative filter efficiencies (σ) were found depending on the mutual orientation of the filter and detection gradients relative to WM fiber orientation. The data demonstrated that the fast-diffusion compartment suppressed by diffusional filtering is not exclusively extra-cellular, but also intra-cellular. While not comprehensive, a simple two-compartment diffusion tensor model with water exchange was able to account qualitatively for the trends in positive and negative filtering efficiencies, while standard model imaging (SMI) without exchange could not. This two-compartment diffusion tensor model also demonstrated smaller AXR variances across subjects. When employing trace-weighted diffusion detection, AXR values were on the order of the R1 (=1/T1) of water at 3T for crossing fibers, while being less than R1 for parallel fibers. CONCLUSION: Orientation-dependent AXR and σ values were observed when using multi-orientation filter and detection gradients in FEXI, indicating that WM FEXI models need to account for compartmental anisotropy. When using trace-weighted detection, AXR values were on the order of or less than R1, complicating the interpretation of FEXI results in WM in terms of biological exchange properties. These findings may contribute toward better understanding of FEXI results in WM.

Spherical convolutional neural networks can improve brain microstructure estimation from diffusion MRI data.

Diffusion magnetic resonance imaging is sensitive to the microstructural properties of brain tissue. However, estimating clinically and scientifically relevant microstructural properties from the measured signals remains a highly challenging inverse problem that machine learning may help solve. This study investigated if recently developed rotationally invariant spherical convolutional neural networks can improve microstructural parameter estimation. We trained a spherical convolutional neural network to predict the ground-truth parameter values from efficiently simulated noisy data and applied the trained network to imaging data acquired in a clinical setting to generate microstructural parameter maps. Our network performed better than the spherical mean technique and multi-layer perceptron, achieving higher prediction accuracy than the spherical mean technique with less rotational variance than the multi-layer perceptron. Although we focused on a constrained two-compartment model of neuronal tissue, the network and training pipeline are generalizable and can be used to estimate the parameters of any Gaussian compartment model. To highlight this, we also trained the network to predict the parameters of a three-compartment model that enables the estimation of apparent neural soma density using tensor-valued diffusion encoding.

Assessment of Precision and Accuracy of Brain White Matter Microstructure using Combined Diffusion MRI and Relaxometry.

Joint modeling of diffusion and relaxation has seen growing interest due to its potential to provide complementary information about tissue microstructure. For brain white matter, we designed an optimal diffusion-relaxometry MRI protocol that samples multiple b-values, B-tensor shapes, and echo times (TE). This variable-TE protocol (27 min) has as subsets a fixed-TE protocol (15 min) and a 2-shell dMRI protocol (7 min), both characterizing diffusion only. We assessed the sensitivity, specificity and reproducibility of these protocols with synthetic experiments and in six healthy volunteers. Compared with the fixed-TE protocol, the variable-TE protocol enables estimation of free water fractions while also capturing compartmental T2 relaxation times. Jointly measuring diffusion and relaxation offers increased sensitivity and specificity to microstructure parameters in brain white matter with voxelwise coefficients of variation below 10%.

Diffusion-weighted MR spectroscopy: Consensus, recommendations, and resources from acquisition to modeling.

Brain cell structure and function reflect neurodevelopment, plasticity, and aging; and changes can help flag pathological processes such as neurodegeneration and neuroinflammation. Accurate and quantitative methods to noninvasively disentangle cellular structural features are needed and are a substantial focus of brain research. Diffusion-weighted MRS (dMRS) gives access to diffusion properties of endogenous intracellular brain metabolites that are preferentially located inside specific brain cell populations. Despite its great potential, dMRS remains a challenging technique on all levels: from the data acquisition to the analysis, quantification, modeling, and interpretation of results. These challenges were the motivation behind the organization of the Lorentz Center workshop on "Best Practices & Tools for Diffusion MR Spectroscopy" held in Leiden, the Netherlands, in September 2021. During the workshop, the dMRS community established a set of recommendations to execute robust dMRS studies. This paper provides a description of the steps needed for acquiring, processing, fitting, and modeling dMRS data, and provides links to useful resources.

Tensor-valued diffusion MRI of human acute stroke.

PURPOSE: Tensor-valued diffusion encoding can disentangle orientation dispersion and subvoxel anisotropy, potentially offering insight into microstructural changes after cerebral ischemia. The purpose was to evaluate tensor-valued diffusion MRI in human acute ischemic stroke, assess potential confounders from diffusion time dependencies, and compare to Monte Carlo diffusion simulations of axon beading. METHODS: Linear (LTE) and spherical (STE) b-tensor encoding with inherently different effective diffusion times were acquired in 21 acute ischemic stroke patients between 3 and 57 h post-onset at 3 T in 2.5 min. In an additional 10 patients, STE with 2 LTE yielding different effective diffusion times were acquired for comparison. Diffusional variance decomposition (DIVIDE) was used to estimate microscopic anisotropy (µFA), as well as anisotropic, isotropic, and total diffusional variance (MKA , MKI , MKT ). DIVIDE parameters, and diffusion tensor imaging (DTI)-derived mean diffusivity and fractional anisotropy (FA) were compared in lesion versus contralateral white matter. Monte Carlo diffusion simulations of various cylindrical geometries for all b-tensor protocols were used to interpret parameter measurements. RESULTS: MD was ˜40% lower in lesions for all LTE/STE protocols. The DIVIDE parameters varied with effective diffusion time: higher µFA and MKA in lesion versus contralateral white matter for STE with longer effective diffusion time LTE, whereas the shorter effective diffusion time LTE protocol yielded lower µFA and MKA in lesions. Both protocols, regardless of diffusion time, were consistent with simulations of greater beading amplitude and intracellular volume fraction. CONCLUSION: DIVIDE parameters depend on diffusion time in acute stroke but consistently indicate neurite beading and larger intracellular volume fraction.

Human prostate MRI at ultrahigh-performance gradient: A feasibility study.

PURPOSE: To demonstrate the technical feasibility and the value of ultrahigh-performance gradient in imaging the prostate in a 3T MRI system. METHODS: In this local institutional review board-approved study, prostate MRI was performed on 4 healthy men. Each subject was scanned in a prototype 3T MRI system with a 42-cm inner-diameter gradient coil that achieves a maximum gradient amplitude of 200 mT/m and slew rate of 500 T/m/s. PI-RADS V2.1-compliant axial T2 -weighted anatomical imaging and single-shot echo planar DWI at standard gradient of 70 mT/m and 150 T/m/s were obtained, followed by DWI at maximum performance (i.e., 200 mT/m and 500 T/m/s). In comparison to state-of-the-art clinical whole-body MRI systems, the high slew rate improved echo spacing from 1020 to 596 µs and, together with a high gradient amplitude for diffusion encoding, TE was reduced from 55 to 36 ms. RESULTS: In all 4 subjects (waist circumference = 81-91 cm, age = 45-65 years), no peripheral nerve stimulation sensation was reported during DWI. Reduced image distortion in the posterior peripheral zone prostate gland and higher signal intensity, such as in the surrounding muscle of high-gradient DWI, were noted. CONCLUSION: Human prostate MRI at simultaneously high gradient amplitude of 200 mT/m and slew rate of 500 T/m/s is feasible, demonstrating that improved gradient performance can address image distortion and T2 decay-induced SNR issues for in vivo prostate imaging.

Diffusion MRI with free gradient waveforms on a high-performance gradient system: Probing restriction and exchange in the human brain.

The dependence of the diffusion MRI signal on the diffusion time carries signatures of restricted diffusion and exchange. Here we seek to highlight these signatures in the human brain by performing experiments using free gradient waveforms designed to be selectively sensitive to the two effects. We examine six healthy volunteers using both strong and ultra-strong gradients (80, 200 and 300 mT/m). In an experiment featuring a large set of 150 gradient waveforms with different sensitivities to restricted diffusion and exchange, our results reveal unique and different time-dependence signatures in grey and white matter. Grey matter was characterised by both restricted diffusion and exchange and white matter predominantly by restricted diffusion. Exchange in grey matter was at least twice as fast as in white matter, across all subjects and all gradient strengths. The cerebellar cortex featured relatively short exchange times (115 ms). Furthermore, we show that gradient waveforms with tailored designs can be used to map exchange in the human brain. We also assessed the feasibility of clinical applications of the method used in this work and found that the exchange-related contrast obtained with a 25-minute protocol at 300 mT/m was preserved in a 4-minute protocol at 300 mT/m and a 10-minute protocol at 80 mT/m. Our work underlines the utility of free waveforms for detecting time dependence signatures due to restricted diffusion and exchange in vivo, which may potentially serve as a tool for studying diseased tissue.

Impact of velocity- and acceleration-compensated encodings on signal dropout and black-blood state in diffusion-weighted magnetic resonance liver imaging at clinical TEs.

PURPOSE: The study aims to develop easy-to-implement concomitant field-compensated gradient waveforms with varying velocity-weighting (M1) and acceleration-weighting (M2) levels and to evaluate their efficacy in correcting signal dropouts and preserving the black-blood state in liver diffusion-weighted imaging. Additionally, we seek to determine an optimal degree of compensation that minimizes signal dropouts while maintaining blood signal suppression. METHODS: Numerically optimized gradient waveforms were adapted using a novel method that allows for the simultaneous tuning of M1- and M2-weighting by changing only one timing variable. Seven healthy volunteers underwent diffusion-weighted magnetic resonance imaging (DWI) with five diffusion encoding schemes (monopolar, velocity-compensated (M1 = 0), acceleration-compensated (M1 = M2 = 0), 84%-M1-M2-compensated, 67%-M1-M2-compensated) at b-values of 50 and 800 s/mm2 at a constant echo time of 70 ms. Signal dropout correction and apparent diffusion coefficients (ADCs) were quantified using regions of interest in the left and right liver lobe. The blood appearance was evaluated using two five-point Likert scales. RESULTS: Signal dropout was more pronounced in the left lobe (19%-42% less signal than in the right lobe with monopolar scheme) and best corrected by acceleration-compensation (8%-10% less signal than in the right lobe). The black-blood state was best with monopolar encodings and decreased significantly (p < 0.001) with velocity- and/or acceleration-compensation. The partially M1-M2-compensated encoding schemes could restore the black-blood state again. Strongest ADC bias occurred for monopolar encodings (difference between left/right lobe of 0.41 µm2/ms for monopolar vs. < 0.12 µm2/ms for the other encodings). CONCLUSION: All of the diffusion encodings used in this study demonstrated suitability for routine DWI application. The results indicate that a perfect value for the level of M1-M2-compensation does not exist. However, among the examined encodings, the 84%-M1-M2-compensated encodings provided a suitable tradeoff.

Probing brain tissue microstructure with MRI: principles, challenges, and the role of multidimensional diffusion-relaxation encoding.

Diffusion MRI uses the random displacement of water molecules to sensitize the signal to brain microstructure and to properties such as the density and shape of cells. Microstructure modeling techniques aim to estimate these properties from acquired data by separating the signal between virtual tissue 'compartments' such as the intra-neurite and the extra-cellular space. A key challenge is that the diffusion MRI signal is relatively featureless compared with the complexity of brain tissue. Another challenge is that the tissue microstructure is wildly different within the gray and white matter of the brain. In this review, we use results from multidimensional diffusion encoding techniques to discuss these challenges and their tentative solutions. Multidimensional encoding increases the information content of the data by varying not only the b-value and the encoding direction but also additional experimental parameters such as the shape of the b-tensor and the echo time. Three main insights have emerged from such encoding. First, multidimensional data contradict common model assumptions on diffusion and T2 relaxation, and illustrates how the use of these assumptions cause erroneous interpretations in both healthy brain and pathology. Second, many model assumptions can be dispensed with if data are acquired with multidimensional encoding. The necessary data can be easily acquired in vivo using protocols optimized to minimize Cramér-Rao lower bounds. Third, microscopic diffusion anisotropy reflects the presence of axons but not dendrites. This insight stands in contrast to current 'neurite models' of brain tissue, which assume that axons in white matter and dendrites in gray matter feature highly similar diffusion. Nevertheless, as an axon-based contrast, microscopic anisotropy can differentiate gray and white matter when myelin alterations confound conventional MRI contrasts.

Coregistered histology sections with diffusion tensor imaging data at 200 µm resolution in meningioma tumors.

A significant problem in diffusion MRI (dMRI) is the lack of understanding regarding which microstructural features account for the variability in the diffusion tensor imaging (DTI) parameters observed in meningioma tumors. A common assumption is that mean diffusivity (MD) and fractional anisotropy (FA) from DTI are inversely proportional to cell density and proportional to tissue anisotropy, respectively. Although these associations have been established across a wide range of tumors, they have been challenged for interpreting within-tumor variations where several additional microstructural features have been suggested as contributing to MD and FA. To facilitate the investigation of the biological underpinnings of DTI parameters, we performed ex-vivo DTI at 200 µm isotropic resolution on sixteen excised meningioma tumor samples. The samples exhibit a variety of microstructural features because the dataset includes meningiomas of six different meningioma types and two different grades. Diffusion-weighted signal (DWI) maps, DWI maps averaged over all directions for given b-value, signal intensities without diffusion encoding (S0) as well as DTI parameters: MD, FA, in-plane FA (FAIP), axial diffusivity (AD) and radial diffusivity (RD), were coregistered to Hematoxylin & Eosin- (H&E) and Elastica van Gieson-stained (EVG) histological sections by a non-linear landmark-based approach. Here, we provide DWI signal and DTI maps coregistered to histology sections and describe the pipeline for processing the raw DTI data and the coregistration. The raw, processed, and coregistered data are hosted by Analytic Imaging Diagnostics Arena (AIDA) data hub registry, and software tools for processing are provided via GitHub. We hope that data can be used in research and education concerning the link between the meningioma microstructure and parameters obtained by DTI.

Ultra-strong diffusion-weighted MRI reveals cerebellar grey matter abnormalities in movement disorders.

Structural brain MRI has proven invaluable in understanding movement disorder pathophysiology. However, most work has focused on grey/white matter volumetric (macrostructural) and white matter microstructural effects, limiting understanding of frequently implicated grey matter microstructural differences. Using ultra-strong spherical tensor encoding diffusion-weighted MRI, a persistent MRI signal was seen in healthy cerebellar grey matter even at high diffusion-weightings (b ​≥ 10,000 s/mm2). Quantifying the proportion of this signal (denoted fs), previously ascertained to originate from inside small spherical spaces, provides a potential proxy for cell body density. In this work, this approach was applied for the first time to a clinical cohort, including patients with diagnosed movement disorders in which the cerebellum has been implicated in symptom pathophysiology. Five control participants (control group 1, median age 24.5 years (20-39 years), imaged at two timepoints, demonstrated consistency in measurement of all three measures - MD (Mean Diffusivity) fs, and Ds (dot diffusivity)- with intraclass correlation coefficients (ICC) of 0.98, 0.86 and 0.76, respectively. Comparison with an older control group (control group 2 (n = 5), median age 51 years (43-58 years)) found no significant differences, neither with morphometric nor microstructural (MD (p = 0.36), fs (p = 0.17) and Ds (p = 0.22)) measures. The movement disorder cohort (Parkinson's Disease, n = 5, dystonia, n = 5. Spinocerebellar Ataxia 6, n = 5) when compared to the age-matched control cohort (Control Group 2) identified significantly lower MD (p < 0.0001 and p < 0.0001) and higher fs values (p < 0.0001 and p < 0.0001) in SCA6 and dystonia cohorts respectively. Lobar division of the cerebellum found these same differences in the superior and inferior posterior lobes, while no differences were seen in either the anterior lobes or with Ds measurements. In contrast to more conventional measures from diffusion tensor imaging, this framework provides enhanced specificity to differences in restricted spherical spaces in grey matter (including small cells) by eliminating signals from cerebrospinal fluid and axons. In the context of human and animal histopathology studies, these findings potentially implicate the cerebellar Purkinje and granule cells as contributors to the observed signal differences, with both cell types having been implicated in several neurological disorders through both postmortem and animal model studies. This novel microstructural imaging approach shows promise for improving movement disorder diagnosis, prognosis, and treatment.

Diffusion MRI with free gradient waveforms on a high-performance gradient system: Probing restriction and exchange in the human brain.

The dependence of the diffusion MRI signal on the diffusion time carries signatures of restricted diffusion and exchange. Here we seek to highlight these signatures in the human brain by performing experiments using free gradient waveforms that are selectively sensitive to the two effects. We examine six healthy volunteers using both strong and ultra-strong gradients (80, 200 and 300 mT/m). In an experiment featuring a large set of gradient waveforms with different sensitivities to restricted diffusion and exchange (150 samples), our results reveal unique time-dependence signatures in grey and white matter, where the former is characterised by both restricted diffusion and exchange and the latter predominantly exhibits restricted diffusion. Furthermore, we show that gradient waveforms with independently varying sensitivities to restricted diffusion and exchange can be used to map exchange in the human brain. We consistently find that exchange in grey matter is at least twice as fast as in white matter, across all subjects and all gradient strengths. The shortest exchange times observed in this study were in the cerebellar cortex (115 ms). We also assess the feasibility of future clinical applications of the method used in this work, where we find that the grey-white matter exchange contrast obtained with a 25-minute 300 mT/m protocol is preserved by a 4-minute 300 mT/m and a 10-minute 80 mT/m protocol. Our work underlines the utility of free waveforms for detecting time-dependence signatures due to restricted diffusion and exchange in vivo, which may potentially serve as a tool for studying diseased tissue.

Meningioma microstructure assessed by diffusion MRI: An investigation of the source of mean diffusivity and fractional anisotropy by quantitative histology.

BACKGROUND: Mean diffusivity (MD) and fractional anisotropy (FA) from diffusion MRI (dMRI) have been associated with cell density and tissue anisotropy across tumors, but it is unknown whether these associations persist at the microscopic level. PURPOSE: To quantify the degree to which cell density and anisotropy, as determined from histology, account for the intra-tumor variability of MD and FA in meningioma tumors. Furthermore, to clarify whether other histological features account for additional intra-tumor variability of dMRI parameters. MATERIALS AND METHODS: We performed ex-vivo dMRI at 200 µm isotropic resolution and histological imaging of 16 excised meningioma tumor samples. Diffusion tensor imaging (DTI) was used to map MD and FA, as well as the in-plane FA (FAIP). Histology images were analyzed in terms of cell nuclei density (CD) and structure anisotropy (SA; obtained from structure tensor analysis) and were used separately in a regression analysis to predict MD and FAIP, respectively. A convolutional neural network (CNN) was also trained to predict the dMRI parameters from histology patches. The association between MRI and histology was analyzed in terms of out-of-sample (R2OS) on the intra-tumor level and within-sample R2 across tumors. Regions where the dMRI parameters were poorly predicted from histology were analyzed to identify features apart from CD and SA that could influence MD and FAIP, respectively. RESULTS: Cell density assessed by histology poorly explained intra-tumor variability of MD at the mesoscopic level (200 µm), as median R2OS = 0.04 (interquartile range 0.01-0.26). Structure anisotropy explained more of the variation in FAIP (median R2OS = 0.31, 0.20-0.42). Samples with low R2OS for FAIP exhibited low variations throughout the samples and thus low explainable variability, however, this was not the case for MD. Across tumors, CD and SA were clearly associated with MD (R2 = 0.60) and FAIP (R2 = 0.81), respectively. In 37% of the samples (6 out of 16), cell density did not explain intra-tumor variability of MD when compared to the degree explained by the CNN. Tumor vascularization, psammoma bodies, microcysts, and tissue cohesivity were associated with bias in MD prediction based solely on CD. Our results support that FAIP is high in the presence of elongated and aligned cell structures, but low otherwise. CONCLUSION: Cell density and structure anisotropy account for variability in MD and FAIP across tumors but cell density does not explain MD variations within the tumor, which means that low or high values of MD locally may not always reflect high or low tumor cell density. Features beyond cell density need to be considered when interpreting MD.

Cardiac q-space trajectory imaging by motion-compensated tensor-valued diffusion encoding in human heart in vivo.

PURPOSE: Tensor-valued diffusion encoding can probe more specific features of tissue microstructure than what is available by conventional diffusion weighting. In this work, we investigate the technical feasibility of tensor-valued diffusion encoding at high b-values with q-space trajectory imaging (QTI) analysis, in the human heart in vivo. METHODS: Ten healthy volunteers were scanned on a 3T scanner. We designed time-optimal gradient waveforms for tensor-valued diffusion encoding (linear and planar) with second-order motion compensation. Data were analyzed with QTI. Normal values and repeatability were investigated for the mean diffusivity (MD), fractional anisotropy (FA), microscopic FA (µFA), isotropic, anisotropic and total mean kurtosis (MKi, MKa, and MKt), and orientation coherence (Cc ). A phantom, consisting of two fiber blocks at adjustable angles, was used to evaluate sensitivity of parameters to orientation dispersion and diffusion time. RESULTS: QTI data in the left ventricular myocardium were MD = 1.62 ± 0.07 µm2 /ms, FA = 0.31 ± 0.03, µFA = 0.43 ± 0.07, MKa = 0.20 ± 0.07, MKi = 0.13 ± 0.03, MKt = 0.33 ± 0.09, and Cc  = 0.56 ± 0.22 (mean ± SD across subjects). Phantom experiments showed that FA depends on orientation dispersion, whereas µFA was insensitive to this effect. CONCLUSION: We demonstrated the first tensor-valued diffusion encoding and QTI analysis in the heart in vivo, along with first measurements of myocardial µFA, MKi, MKa, and Cc . The methodology is technically feasible and provides promising novel biomarkers for myocardial tissue characterization.

A preliminary choroid plexus volumetric study in individuals with psychosis.

The choroid plexus (ChP) is part of the blood-cerebrospinal fluid barrier, regulating brain homeostasis and the brain's response to peripheral events. Its upregulation and enlargement are considered essential in psychosis. However, the timing of the ChP enlargement has not been established. This study introduces a novel magnetic resonance imaging-based segmentation method to examine ChP volumes in two cohorts of individuals with psychosis. The first sample consists of 41 individuals with early course psychosis (mean duration of illness = 1.78 years) and 30 healthy individuals. The second sample consists of 30 individuals with chronic psychosis (mean duration of illness = 7.96 years) and 34 healthy individuals. We utilized manual segmentation to measure ChP volumes. We applied ANCOVAs to compare normalized ChP volumes between groups and partial correlations to investigate the relationship between ChP, LV volumes, and clinical characteristics. Our segmentation demonstrated good reliability (.87). We further showed a significant ChP volume increase in early psychosis (left: p < .00010, right: p < .00010) and a significant positive correlation between higher ChP and higher LV volumes in chronic psychosis (left: r = .54, p = .0030, right: r = .68; p < .0010). Our study suggests that ChP enlargement may be a marker of acute response around disease onset. It might also play a modulatory role in the chronic enlargement of lateral ventricles, often reported in psychosis. Future longitudinal studies should investigate the dynamics of ChP enlargement as a promising marker for novel therapeutic strategies.

Optimal experimental design and estimation for q-space trajectory imaging.

Tensor-valued diffusion encoding facilitates data analysis by q-space trajectory imaging. By modeling the diffusion signal of heterogeneous tissues with a diffusion tensor distribution (DTD) and modulating the encoding tensor shape, this novel approach allows disentangling variations in diffusivity from microscopic anisotropy, orientation dispersion, and mixtures of multiple isotropic diffusivities. To facilitate the estimation of the DTD parameters, a parsimonious acquisition scheme coupled with an accurate and precise estimation of the DTD is needed. In this work, we create two precision-optimized acquisition schemes: one that maximizes the precision of the raw DTD parameters, and another that maximizes the precision of the scalar measures derived from the DTD. The improved precision of these schemes compared to a naïve sampling scheme is demonstrated in both simulations and real data. Furthermore, we show that the weighted linear least squares (WLLS) estimator that uses the squared reciprocal of the noisy signal as weights can be biased, whereas the iteratively WLLS estimator with the squared reciprocal of the predicted signal as weights outperforms the conventional unweighted linear LS and nonlinear LS estimators in terms of accuracy and precision. Finally, we show that the use of appropriate constraints can considerably increase the precision of the estimator with only a limited decrease in accuracy.

Diffusion MRI with pulsed and free gradient waveforms: Effects of restricted diffusion and exchange.

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 and exchange rates in the simulated range of 0 to 20 s - 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.

Tensor-valued diffusion magnetic resonance imaging in a radiotherapy setting.

Background and purpose: Diagnostic information about cell density variations and microscopic tissue anisotropy can be gained from tensor-valued diffusion magnetic resonance imaging (MRI). These properties of tissue microstructure have the potential to become novel imaging biomarkers for radiotherapy response. However, tensor-valued diffusion encoding is more demanding than conventional encoding, and its compatibility with MR scanners that are dedicated to radiotherapy has not been established. Thus, our aim was to investigate the feasibility of tensor-valued diffusion MRI with radiotherapy dedicated MR equipment. Material and methods: A tensor-valued diffusion protocol was implemented, and five healthy volunteers were scanned with different resolutions using conventional head coil and radiotherapy coil setup with fixation masks. Signal-to-noise-ratio (SNR) was evaluated to assess the risk of signal bias due to rectified noise floor. We also evaluated the repeatability and reproducibility of the microstructure parameters. One patient with brain metastasis was scanned to investigate the image quality and the transferability of the setup to diseased tissue. Results: A resolution of 3 × 3 × 3 mm3 provided images with SNR > 3 for 93 % of the voxels using radiotherapy coil setup. The parameter maps and repeatability characteristics were comparable to those observed with a conventional head coil. The patient evaluation demonstrated successful parameter analysis also in tumor tissue, with SNR > 3 for 93 % of the voxels. Conclusion: We demonstrate that tensor-valued diffusion MRI is compatible with radiotherapy fixation masks and coil setup for investigations of microstructure parameters. The reported reproducibility may be used to plan future investigations of imaging biomarkers in brain cancer radiotherapy.

MR Fingerprinting with b-Tensor Encoding for Simultaneous Quantification of Relaxation and Diffusion in a Single Scan.

PURPOSE: Although both relaxation and diffusion imaging are sensitive to tissue microstructure, studies have reported limited sensitivity and robustness of using relaxation or conventional diffusion alone to characterize tissue microstructure. Recently, it has been shown that tensor-valued diffusion encoding and joint relaxation-diffusion quantification enable more reliable quantification of compartment-specific microstructural properties. However, scan times to acquire such data can be prohibitive. Here, we aim to simultaneously quantify relaxation and diffusion using MR fingerprinting (MRF) and b-tensor encoding in a clinically feasible time. METHODS: We developed multidimensional MRF scans (mdMRF) with linear and spherical b-tensor encoding (LTE and STE) to simultaneously quantify T1, T2, and ADC maps from a single scan. The image quality, accuracy, and scan efficiency were compared between the mdMRF using LTE and STE. Moreover, we investigated the robustness of different sequence designs to signal errors and their impact on the maps. RESULTS: T1 and T2 maps derived from the mdMRF scans have consistently high image quality, while ADC maps are sensitive to different sequence designs. Notably, the fast imaging steady state precession (FISP)-based mdMRF scan with peripheral pulse gating provides the best ADC maps that are free of image distortion and shading artifacts. CONCLUSION: We demonstrated the feasibility of quantifying T1, T2, and ADC maps simultaneously from a single mdMRF scan in around 24 s/slice. The map quality and quantitative values are consistent with the reference scans.