Evaluating metabolites in patients with major depressive disorder who received mindfulness-based cognitive therapy and healthy controls using short echo MRSI at 7 Tesla

Objectives

Our aim was to evaluate differences in metabolite levels between unmedicated patients with major depressive disorder (MDD) and healthy controls, to assess changes in metabolites in patients after they completed an 8-week course of mindfulness-based cognitive therapy (MBCT), and to exam the correlation between metabolites and depression severity.

Materials and methods

Sixteen patients with MDD and ten age- and gender-matched healthy controls were studied using 3D short echo-time (20 ms) magnetic resonance spectroscopic imaging (MRSI) at 7 Tesla. Relative metabolite ratios were estimated in five regions of interest corresponding to insula, anterior cingulate cortex (ACC), caudate, putamen, and thalamus.

Results

In all cases, MBCT reduced severity of depression. The ratio of total choline-containing compounds/total creatine (tCr) in the right caudate was significantly increased compared to that in healthy controls, while ratios of N-acetyl aspartate (NAA)/tCr in the left ACC, myo-inositol/tCr in the right insula, and glutathione/tCr in the left putamen were significantly decreased. At baseline, the severity of depression was negatively correlated with my-inositol/tCr in the left insula and putamen. The improvement in depression severity was significantly associated with changes in NAA/tCr in the left ACC.

Conclusions

This study has successfully evaluated regional differences in metabolites for patients with MDD who received MBCT treatment and in controls using 7 Tesla MRSI.

     

Non-gadolinium dynamic angiography of the neurovasculature using arterial spin labeling MRI: preliminary experience in children

Objective

We demonstrate the potential clinical utility of a 4D non-gadolinium dynamic angiography technique based on arterial spin-labeling called contrast inherent inflow enhanced multi-phase angiography (CINEMA) in pediatric patients.

Materials and Methods

CINEMA was qualitatively compared to conventional time-of-flight (TOF) angiography in a cohort of 31 pediatric patients at 3 Tesla.

Results

CINEMA data were successfully acquired and reconstructed in all patients with no image artifacts. There were no cases where CINEMA was rated inferior to TOF in depicting intracranial vessel conspicuity. In 19 cases, CINEMA was rated equivalent to TOF and in the 12 remaining cases CINEMA was rated superior to TOF.

Conclusion

There is a steadily rising concern in adults and children over the potential effects of intracranial deposition of gadolinium. CINEMA is therefore a viable alternative in dynamic neurovascular imaging.

     

Phosphorus magnetic resonance spectroscopic imaging using flyback echo planar readout trajectories

Object

To present and evaluate a fast phosphorus magnetic resonance spectroscopic imaging (MRSI) sequence using echo planar spectroscopic imaging with flyback readout gradient trajectories.

Materials and Methods

Waveforms were designed and implemented using a 3 Tesla MRI system. ³¹P spectra were acquired with 2 × 2 cm² and 3 × 3 cm² resolution over a 20- and 21-cm field of view and spectral bandwidths up to 1923 Hz. The sequence was first tested using a 20-cm-diameter phosphate phantom, and subsequent in vivo tests were performed on healthy human calf muscles and brains from five volunteers.

Results

Flyback EPSI achieved 10× and 7× reductions in acquisition time, with 68.0 ± 1.2 and 69.8 ± 2.2% signal-to-noise ratio (SNR) per unit of time efficiency (theoretical SNR efficiency was 74.5 and 76.4%) for the in vivo experiments, compared to conventional phase-encoded MRSI for the 2 × 2 cm² and 3 × 3 cm² resolution waveforms, respectively. Statistical analysis showed no difference in the quantification of most metabolites. Time savings and SNR comparisons were consistent across phantom, leg and brain experiments.

Conclusion

EPSI using flyback readout trajectories was found to be a reliable alternative for acquiring ³¹P-MRSI data in a shorter acquisition time.

     

A comparison of fitting algorithms for diffusion-weighted MRI data analysis using an intravoxel incoherent motion model

Object

The objective of this study is to propose a modified VARiable PROjection (VARPRO) algorithm specifically tailored for fitting the intravoxel incoherent motion (IVIM) model to diffusion-weighted magnetic resonance imaging (DW-MRI) data from locally advanced rectal cancer (LARC).

Materials and methods

The proposed algorithm is compared with classical non-linear least squares (NLLS) analysis using the Levenberg-Marquardt (LM) algorithm and with two recently proposed algorithms for ‘segmented’ analysis. These latter two comprise two consecutive steps: first, a subset of parameters is estimated using a portion of data; second, the remaining parameters are estimated using the whole data and the previous estimates. The comparison between the algorithms was based on the \(R^2\) goodness-of-fit measure: performance analysis was carried out on real data obtained by DW-MRI on 40 LARC patients.

Results

The performance of the proposed algorithm was higher than that of LM in 64 % of cases; ‘segmented’ methods were poorer than our algorithm in 100 % of cases.

Conclusion

The proposed modified VARPRO algorithm can lead to better fit of the IVIM model to LARC DW-MRI data compared to other techniques.

     

Theory of MRI contrast in the annulus fibrosus of the intervertebral disc

Objective

Here we develop a three-dimensional analytic model for MR image contrast of collagen lamellae in the annulus fibrosus of the intervertebral disc of the spine, based on the dependence of the MRI signal on collagen fiber orientation.

Materials and methods

High-resolution MRI scans were performed at 1.5 and 7 T on intact whole disc specimens from ovine, bovine, and human spines. An analytic model that approximates the three-dimensional curvature of the disc lamellae was developed to explain inter-lamellar contrast and intensity variations in the annulus. The model is based on the known anisotropic dipolar relaxation of water in tissues with ordered collagen.

Results

Simulated MRI data were generated that reproduced many features of the actual MRI data. The calculated inter-lamellar image contrast demonstrated a strong dependence on the collagen fiber angle and on the circumferential location within the annulus.

Conclusion

This analytic model may be useful for interpreting MR images of the disc and for predicting experimental conditions that will optimize MR image contrast in the annulus fibrosus.

     

ESMRMB 2017, 34th Annual Scientific Meeting,Barcelona, ES,October 19–October 21: Abstracts,Thursday

Objective

Current magnetic resonance imaging (MRI) axon diameter measurements rely on the pulsed gradient spin-echo sequence, which is unable to provide diffusion times short enough to measure small axon diameters. This study combines the AxCaliber axon diameter fitting method with data generated from Monte Carlo simulations of oscillating gradient spin-echo sequences (OGSE) to infer micron-sized axon diameters, in order to determine the feasibility of using MRI to infer smaller axon diameters in brain tissue.

Materials and methods

Monte Carlo computer simulation data were synthesized from tissue geometries of cylinders of different diameters using a range of gradient frequencies in the cosine OGSE sequence . Data were fitted to the AxCaliber method modified to allow the new pulse sequence. Intra- and extra-axonal water were studied separately and together.

Results

The simulations revealed the extra-axonal model to be problematic. Rather than change the model, we found that restricting the range of gradient frequencies such that the measured apparent diffusion coefficient was constant over that range resulted in more accurate fitted diameters. Thus a careful selection of frequency ranges is needed for the AxCaliber method to correctly model extra-axonal water, or adaptations to the method are needed. This restriction helped reduce the necessary gradient strengths for measurements that could be performed with parameters feasible for a Bruker BG6 gradient set. For these experiments, the simulations inferred diameters as small as 0.5 μm on square-packed and randomly packed cylinders. The accuracy of the inferred diameters was found to be dependent on the signal-to-noise ratio (SNR), with smaller diameters more affected by noise, although all diameter distributions were distinguishable from one another for all SNRs tested.

Conclusion

The results of this study indicate the feasibility of using MRI with OGSE on preclinical scanners to infer small axon diameters.

     

Distortion correction of echo planar images applying the concept of finite rate of innovation to point spread function mapping (FRIP)

Objective

Point spread function (PSF) mapping enables estimating the displacement fields required for distortion correction of echo planar images. Recently, a highly accelerated approach was introduced for estimating displacements from the phase slope of under-sampled PSF mapping data. Sampling schemes with varying spacing were proposed requiring stepwise phase unwrapping. To avoid unwrapping errors, an alternative approach applying the concept of finite rate of innovation to PSF mapping (FRIP) is introduced, using a pattern search strategy to locate the PSF peak, and the two methods are compared.

Materials and methods

Fully sampled PSF data was acquired in six subjects at 3.0 T, and distortion maps were estimated after retrospective under-sampling. The two methods were compared for both previously published and newly optimized sampling patterns. Prospectively under-sampled data were also acquired. Shift maps were estimated and deviations relative to the fully sampled reference map were calculated.

Results

The best performance was achieved when using FRIP with a previously proposed sampling scheme. The two methods were comparable for the remaining schemes. The displacement field errors tended to be lower as the number of samples or their spacing increased.

Conclusion

A robust method for estimating the position of the PSF peak has been introduced.

     

Spectroscopic sampling of the left side of long-TE spin echoes: a free lunch?

Objective

Use of spectroscopically-acquired spin echoes typically involves Fourier transformation of the right side of the echo while largely neglecting the left side. For sufficiently long echo times, the left side may have enough spectral resolution to offer some utility. Since the acquisition of this side is “free”, we deemed it worthy of attention and investigated the spectral properties and information content of this data.

Materials and methods

Theoretical expressions for left- and right-side spectra were derived assuming Lorentzian frequency distributions. For left-side spectra, three regimes were identified based upon the relative magnitudes of reversible and irreversible transverse relaxation rates, R ₂′ and R ₂, respectively. Point-resolved spectroscopy (PRESS) data from muscle, fat deposit and bone marrow were acquired at 1.5 T to test aspects of the theoretical expressions.

Results

For muscle water or methylene marrow resonances, left-side signals were substantially or moderately larger than right-side signals but were similar in magnitude for muscle choline and creatine resonances. Left- versus right-side spectral-peak amplitude ratios depend sensitively on the relative values of R ₂ and R ₂′, which can be estimated given this ratio and a right-side linewidth measurement.

Conclusion

Left-side spectra can be used to augment signal-to-noise and to estimate spectral R ₂ and R ₂′ values under some circumstances.

     

<Emphasis Type="Italic">K</Emphasis>-space data processing for magnetic resonance elastography (MRE)

Objective

Magnetic resonance elastography (MRE) requires substantial data processing based on phase image reconstruction, wave enhancement, and inverse problem solving. The objective of this study is to propose a new, fast MRE method based on MR raw data processing, particularly adapted to applications requiring fast MRE measurement or high elastogram update rate.

Materials and methods

The proposed method allows measuring tissue elasticity directly from raw data without prior phase image reconstruction and without phase unwrapping. Experimental feasibility is assessed both in a gelatin phantom and in the liver of a porcine model in vivo. Elastograms are reconstructed with the raw MRE method and compared to those obtained using conventional MRE. In a third experiment, changes in elasticity are monitored in real-time in a gelatin phantom during its solidification by using both conventional MRE and raw MRE.

Results

The raw MRE method shows promising results by providing similar elasticity values to the ones obtained with conventional MRE methods while decreasing the number of processing steps and circumventing the delicate step of phase unwrapping. Limitations of the proposed method are the influence of the magnitude on the elastogram and the requirement for a minimum number of phase offsets.

Conclusion

This study demonstrates the feasibility of directly reconstructing elastograms from raw data.

     

Automatic gallbladder segmentation using combined 2D and 3D shape features to perform volumetric analysis in native and secretin-enhanced MRCP sequences

Objectives

We aimed to develop the first fully automated 3D gallbladder segmentation approach to perform volumetric analysis in volume data of magnetic resonance (MR) cholangiopancreatography (MRCP) sequences. Volumetric gallbladder analysis is performed for non-contrast-enhanced and secretin-enhanced MRCP sequences.

Materials and methods

Native and secretin-enhanced MRCP volume data were produced with a 1.5-T MR system. Images of coronal maximum intensity projections (MIP) are used to automatically compute 2D characteristic shape features of the gallbladder in the MIP images. A gallbladder shape space is generated to derive 3D gallbladder shape features, which are then combined with 2D gallbladder shape features in a support vector machine approach to detect gallbladder regions in MRCP volume data. A region-based level set approach is used for fine segmentation. Volumetric analysis is performed for both sequences to calculate gallbladder volume differences between both sequences.

Results

The approach presented achieves segmentation results with mean Dice coefficients of 0.917 in non-contrast-enhanced sequences and 0.904 in secretin-enhanced sequences.

Conclusion

This is the first approach developed to detect and segment gallbladders in MR-based volume data automatically in both sequences. It can be used to perform gallbladder volume determination in epidemiological studies and to detect abnormal gallbladder volumes or shapes. The positive volume differences between both sequences may indicate the quantity of the pancreatobiliary reflux.

     

Reproducibility of locus coeruleus and substantia nigra imaging with neuromelanin sensitive MRI

Objectives

The purpose of this study was to assess the reproducibility of substantia nigra pars compacta (SNpc) and locus coeruleus (LC) delineation and measurement with neuromelanin-sensitive MRI.

Materials and methods

Eleven subjects underwent two neuromelanin-sensitive MRI scans. SNpc and LC volumes were extracted for each scan. Reproducibility of volume and magnetization transfer contrast measurements in SNpc and LC was assessed using intraclass correlation coefficients (ICC) and dice similarity coefficients (DSC).

Results

SNpc and LC volume measurements showed excellent reproducibility (SNpc-ICC: 0.94, p < 0.001; LC-ICC: 0.96, p < 0.001). SNpc and LC were accurately delineated between scans (SNpc-DSC: 0.80 ± 0.03; LC-DSC: 0.63 ± 0.07).

Conclusion

Neuromelanin-sensitive MRI can consistently delineate SNpc and LC.

     

In vivo functional connectome of human brainstem nuclei of the ascending arousal,autonomic, and motor systems by high spatial resolution 7-Tesla fMRI

Objective

Our aim was to map the in vivo human functional connectivity of several brainstem nuclei with the rest of the brain by using seed-based correlation of ultra-high magnetic field functional magnetic resonance imaging (fMRI) data.

Materials and methods

We used the recently developed template of 11 brainstem nuclei derived from multi-contrast structural MRI at 7 Tesla as seed regions to determine their connectivity to the rest of the brain. To achieve this, we used the increased contrast-to-noise ratio of 7-Tesla fMRI compared with 3 Tesla and time-efficient simultaneous multi-slice imaging to cover the brain with high spatial resolution (1.1-mm isotropic nominal resolution) while maintaining a short repetition time (2.5 s).

Results

The delineated Pearson’s correlation-based functional connectivity diagrams (connectomes) of 11 brainstem nuclei of the ascending arousal, motor, and autonomic systems from 12 controls are presented and discussed in the context of existing histology and animal work.

Conclusion

Considering that the investigated brainstem nuclei play a crucial role in several vital functions, the delineated preliminary connectomes might prove useful for future in vivo research and clinical studies of human brainstem function and pathology, including disorders of consciousness, sleep disorders, autonomic disorders, Parkinson’s disease, and other motor disorders.

     

Improving the clinical potential of ultra-high field fMRI using a model-free analysis method based on response consistency

Objective

To develop an analysis method that is sensitive to non-model-conform responses often encountered in ultra-high field presurgical planning fMRI. Using the consistency of time courses over a number of experiment repetitions, it should exclude low quality runs and generate activation maps that reflect the reliability of responses.

Materials and methods

7 T fMRI data were acquired from six healthy volunteers: three performing purely motor tasks and three a visuomotor task. These were analysed with the proposed approach (UNBIASED) and the GLM.

Results

UNBIASED results were generally less affected by false positive results than the GLM. Runs that were identified as being of low quality were confirmed to contain little or no activation. In two cases, regions were identified as activated in UNBIASED but not GLM results. Signal changes in these areas were time-locked to the task, but were delayed or transient.

Conclusion

UNBIASED is shown to be a reliable means of identifying consistent task-related signal changes regardless of response timing. In presurgical planning, UNBIASED could be used to rapidly generate reliable maps of the consistency with which eloquent brain regions are activated without recourse to task timing and despite modified hemodynamics.

     

The traveling heads: multicenter brain imaging at 7 Tesla

Objective

This study evaluates the inter-site and intra-site reproducibility of 7 Tesla brain imaging and compares it to literature values for other field strengths.

Materials and methods

The same two subjects were imaged at eight different 7 T sites. MP2RAGE, TSE, TOF, SWI, EPI as well as B1 and B0 field maps were analyzed quantitatively to assess inter-site reproducibility. Intra-site reproducibility was measured with rescans at three sites.

Results

Quantitative measures of MP2RAGE scans showed high agreement. Inter-site and intra-site reproducibility errors were comparable to 1.5 and 3 T. Other sequences also showed high reproducibility between the sites, but differences were also revealed. The different RF coils used were the main source for systematic differences between the sites.

Conclusion

Our results show for the first time that multi-center brain imaging studies of the supratentorial brain can be performed at 7 T with high reproducibility and similar reliability as at 3T. This study develops the basis for future large-scale 7 T multi-site studies.

     

3D true-phase polarity recovery with independent phase estimation using three-tier stacks based region growing (3D-TRIPS)

Objectives

A postprocessing technique termed 3D true-phase polarity recovery with independent phase estimation using three-tier stacks based region growing (3D-TRIPS) was developed, which directly reconstructs phase-sensitive inversion-recovery images without acquisition of phase-reference images. The utility of this technique is demonstrated in myocardial late gadolinium enhancement (LGE) imaging.

Materials and methods

A data structure with three tiers of stacks was used for 3D-TRIPS to directly achieve reliable region growing for successful background-phase estimation. Fifteen patients undergoing postgadolinium 3D phase-sensitive inversion recovery (PSIR) cardiac LGE magnetic resonance imaging (MRI) were recruited, and 3D-TRIPS LGE reconstructions were compared with standard PSIR. Objective voxel-by-voxel comparison was performed. Additionally, blinded review by two radiologists compared scar visibility, clinical acceptability, voxel polarity error, or groups and blurring.

Results

3D-TRIPS efficiently reconstructed postcontrast phase-sensitive myocardial LGE images. Objective analysis showed an average 95% voxel-by-voxel agreement between 3D-TRIPS and PSIR images. Blinded radiologist review demonstrated similar image quality between 3D-TRIPS and PSIR reconstruction.

Conclusion

3D-TRIPS provided similar image quality to PSIR for phase-sensitive myocardial LGE MRI reconstruction. 3D-TRIPS does not require acquisition of a reference image and can therefore be used to accelerate phase-sensitive LGE imaging.

     

The effect of diffusion gradient direction number on corticospinal tractography in the human brain: an along-tract analysis

Objectives

We evaluated diffusion imaging measures of the corticospinal tract obtained with a probabilistic tractography algorithm applied to data of two acquisition protocols based on different numbers of diffusion gradient directions (NDGDs).

Materials and methods

The corticospinal tracts (CST) of 18 healthy subjects were delineated using 22 and 66-NDGD data. An along-tract analysis of diffusion metrics was performed to detect possible local differences due to NDGD.

Results

FA values at 22-NDGD showed an increase along the central portion of the CST. The mean of partial volume fraction of the orientation of the second fiber (f2) was higher at 66-NDGD bilaterally, because for 66-NDGD data the algorithm more readily detects dominant fiber directions beyond the first, thus the increase in FA at 22-NDGD is due to a substantially reduced detection of crossing fiber volume. However, the good spatial correlation between the tracts drawn at 22 and 66 NDGD shows that the extent of the tract can be successfully defined even at lower NDGD.

Conclusions

Given the spatial tract localization obtained even at 22-NDGD, local analysis of CST can be performed using a NDGD compatible with clinical protocols. The probabilistic approach was particularly powerful in evaluating crossing fibers when present.

     

Assessment of pharmacokinetics for microvessel proliferation by DCE-MRI for early detection of physeal bone bridge formation in an animal model

Objectives

Bone bridge formation occurs after physeal lesions and can lead to growth arrest if not reversed. Previous investigations on the underlying mechanisms of this formation used histological methods. Therefore, this study aimed to apply a minimally invasive method using dynamic contrast-enhanced MRI (DCE-MRI).

Materials and methods

Changes in functional parameters related to the microvessel system were assessed in a longitudinal study of a cohort of an animal model applying a reference region model. The development of morphology of the injured physis was investigated with 3D high-resolution MRI. To acquire complementary information for MRI-related findings qRT-PCR and immunohistochemical data were acquired for a second cohort of the animal model.

Results

The evaluation of the pharmacokinetic parameters showed a first rise of the transfer coefficient 7 days post-lesion and a maximum 42 days after operation. The analysis of the complementary data showed a connection of the first rise to microvessel proliferation while the maximum value was linked to bone remodeling.

Conclusion

The pharmacokinetic analysis of DCE-MRI provides information on a proliferation of microvessels during the healing process as a sign for bone bridge formation. Thereby, DCE-MRI could identify details, which up to now required analyses of highly invasive methods.

     

Optimization of diffusion-weighted single-refocused spin-echo EPI by reducing eddy-current artifacts and shortening the echo time

Objective

The purpose of this work was to optimize the acquisition of diffusion-weighted (DW) single-refocused spin-echo (srSE) data without intrinsic eddy-current compensation (ECC) for an improved performance of ECC postprocessing. The rationale is that srSE sequences without ECC may yield shorter echo times (TE) and thus higher signal-to-noise ratios (SNR) than srSE or twice-refocused spin-echo (trSE) schemes with intrinsic ECC.

Materials and methods

The proposed method employs dummy scans with DW gradients to drive eddy currents into a steady state before data acquisition. Parameters of the ECC postprocessing algorithm were also optimized. Simulations were performed to obtain minimum TE values for the proposed sequence and sequences with intrinsic ECC. Experimentally, the proposed method was compared with standard DW-trSE imaging, both in vitro and in vivo.

Results

Simulations showed substantially shorter TE for the proposed method than for methods with intrinsic ECC when using shortened echo readouts. Data of the proposed method showed a marked increase in SNR. A dummy scan duration of at least 1.5 s improved performance of the ECC postprocessing algorithm.

Conclusion

Changes proposed for the DW-srSE sequence and for the parameter setting of the postprocessing ECC algorithm considerably reduced eddy-current artifacts and provided a higher SNR.

     

Single-breath-hold 3-D CINE imaging of the left ventricle using Cartesian sampling

Objectives

Our objectives were to evaluate a single-breath-hold approach for Cartesian 3-D CINE imaging of the left ventricle with a nearly isotropic resolution of \(1.9 \times 1.9 \times 2.5\,{\text {mm}^3}\) and a breath-hold duration of \(\sim \)19 s against a standard stack of 2-D CINE slices acquired in multiple breath-holds. Validation is performed with data sets from ten healthy volunteers.

Materials and methods

A Cartesian sampling pattern based on the spiral phyllotaxis and a compressed sensing reconstruction method are proposed to allow 3-D CINE imaging with high acceleration factors. The fully integrated reconstruction uses multiple graphics processing units to speed up the reconstruction. The 2-D CINE and 3-D CINE are compared based on ventricular function parameters, contrast-to-noise ratio and edge sharpness measurements.

Results

Visual comparisons of corresponding short-axis slices of 2-D and 3-D CINE show an excellent match, while 3-D CINE also allows reformatting to other orientations. Ventricular function parameters do not significantly differ from values based on 2-D CINE imaging. Reconstruction times are below 4 min.

Conclusion

We demonstrate single-breath-hold 3-D CINE imaging in volunteers and three example patient cases, which features fast reconstruction and allows reformatting to arbitrary orientations.

     

Fitting interrelated datasets: metabolite diffusion and general lineshapes

Objective

Simultaneous modeling of true 2-D spectroscopy data, or more generally, interrelated spectral datasets has been described previously and is useful for quantitative magnetic resonance spectroscopy applications. In this study, a combined method of reference-lineshape enhanced model fitting and two-dimensional prior-knowledge fitting for the case of diffusion weighted MR spectroscopy is presented.

Materials and methods

Time-dependent field distortions determined from a water reference are applied to the spectral bases used in linear-combination modeling of interrelated spectra. This was implemented together with a simultaneous spectral and diffusion model fitting in the previously described Fitting Tool for Arrays of Interrelated Datasets (FiTAID), where prior knowledge conditions and restraints can be enforced in two dimensions.

Results

The benefit in terms of increased accuracy and precision of parameters is illustrated with examples from Monte Carlo simulations, in vitro and in vivo human brain scans for one- and two-dimensional datasets from 2-D separation, inversion recovery and diffusion-weighted spectroscopy (DWS). For DWS, it was found that acquisitions could be substantially shortened.

Conclusion

It is shown that inclusion of a measured lineshape into modeling of interrelated MR spectra is beneficial and can be combined also with simultaneous spectral and diffusion modeling.