Diffusion-weighted MRI of the lung at 3T evaluated using echo-planar-based and single-shot turbo spin-echo-based acquisition techniques for radiotherapy applications

Purpose

To compare single‐shot echo‐planar (SS‐EPI)‐based and turbo spin‐echo (SS‐TSE)‐based diffusion‐weighted imaging (DWI) in Non‐Small Cell Lung Cancer (NSCLC) patients and to characterize the distributions of apparent diffusion coefficient (ADC) values generated by the two techniques.

Methods

Ten NSCLC patients were enrolled in a prospective IRB‐approved study to compare and optimize DWI using EPI and TSE‐based techniques for radiotherapy planning. The imaging protocol included axial T2w, EPI‐based DWI and TSE‐based DWI on a 3 T Philips scanner. Both EPI‐based and TSE‐based DWI sequences used three b values (0, 400, and 800 s/mm²). The acquisition times for EPI‐based and TSE‐based DWI were 5 and 8 min, respectively. DW‐MR images were manually coregistered with axial T2w images, and tumor volume contoured on T2w images were mapped onto the DWI scans. A pixel‐by‐pixel fit of tumor ADC was calculated based on monoexponential signal behavior. Tumor ADC mean, standard deviation, kurtosis, and skewness were calculated and compared between EPI and TSE‐based DWI. Image distortion and ADC values between the two techniques were also quantified using fieldmap analysis and a NIST traceable ice‐water diffusion phantom, respectively.

Results

The mean ADC for EPI and TSE‐based DWI were 1.282 ± 0.42 × 10⁻³ and 1.211 ± 0.31 × 10⁻³ mm²/s. The average skewness and kurtosis were 0.14 ± 0.4 and 2.43 ± 0.40 for DWI‐EPI and −0.06 ± 0.69 and 2.89 ± 0.62 for DWI‐TSE. Fieldmap analysis showed a mean distortion of 13.72 ± 8.12 mm for GTV for DWI‐EPI and 0.61 ± 0.4 mm for DWI‐TSE. ADC values obtained using the diffusion phantom for the two techniques were within 0.03 × 10⁻³ mm²/s with respect to each other as well as the established values.

Conclusions

Diffusion‐weighted turbo spin‐echo shows better geometrical accuracy compared to DWI‐EPI. Mean ADC values were similar with both acquisitions but the shape of the histograms was different based on the skewness and kurtosis values. The impact of differences in respiratory technique on ADC values requires further investigation.

Navigator-Free EPI Ghost Correction With Structured Low-Rank Matrix Models: New Theory and Methods

Structured low-rank matrix models have previously been introduced to enable calibrationless MR image reconstruction from sub-Nyquist data, and such ideas have recently been extended to enable navigator-free echo-planar imaging (EPI) ghost correction. This paper presents a novel theoretical analysis which shows that, because of uniform subsampling, the structured low-rank matrix optimization problems for EPI data will always have either undesirable or non-unique solutions in the absence of additional constraints. This theory leads us to recommend and investigate problem formulations for navigator-free EPI that incorporate side information from either image-domain or k-space domain parallel imaging methods. The importance of using nonconvex low-rank matrix regularization is also identified. We demonstrate using phantom and in vivo data that the proposed methods are able to eliminate ghost artifacts for several navigator-free EPI acquisition schemes, obtaining better performance in comparison with the state-of-the-art methods across a range of different scenarios. Results are shown for both single-channel acquisition and highly accelerated multi-channel acquisition.

Tilted-CAIPI for highly accelerated distortion-free EPI with point spread function (PSF) encoding

PURPOSE

To develop a method for fast distortion- and blurring-free imaging.

THEORY

EPI with point-spread-function (PSF) mapping can achieve distortion- and blurring-free imaging at a cost of long acquisition time. In this study, an acquisition/reconstruction technique, termed "tilted-CAIPI," is proposed to achieve >20× acceleration for PSF-EPI. The proposed method systematically optimized the k-space sampling trajectory with B₀ -inhomogeneity-informed reconstruction, to exploit the inherent signal correlation in PSF-EPI and take full advantage of coil sensitivity. Susceptibility-induced phase accumulation is regarded as an additional encoding that is estimated by calibration data and integrated into reconstruction. Self-navigated phase correction was developed to correct shot-to-shot phase variation in diffusion imaging.

METHODS

Tilted-CAIPI was implemented at 3T, with incorporation of partial Fourier and simultaneous multislice to achieve further accelerations. T₂ -weighted, T₂ ^* -weighted, and diffusion-weighted imaging experiments were conducted to evaluate the proposed method.

RESULTS

The ability of tilted-CAIPI to provide highly accelerated imaging without distortion and blurring was demonstrated through in vivo brain experiments, where only 8 shots per simultaneous slice group were required to provide high-quality, high-SNR imaging at 0.8-1 mm resolution.

CONCLUSION

Tilted-CAIPI achieved fast distortion- and blurring-free imaging with high SNR. Whole-brain T₂ -weighted, T₂ ^* -weighted, and diffusion imaging can be obtained in just 15-60 s.

ROBUST AUTOCALIBRATED LORAKS FOR EPI GHOST CORRECTION

Nyquist ghosts are a longstanding problem in a variety of fast MRI experiments that use echo-planar imaging (EPI). Recently, several structured low-rank matrix modeling approaches have been proposed that achieve state-of-the-art ghost-elimination, although the performance of these approaches is still inadequate in several important scenarios. We present a new structured low-rank matrix recovery ghost correction method that we call Robust Autocalibrated LORAKS (RAC-LORAKS). RAC-LORAKS incorporates constraints from autocalibration data to avoid ill-posedness, but allows adaptation of these constraints to gain robustness against possible autocalibration imperfections. RAC-LORAKS is tested in two challenging scenarios: highly-undersampled multi-channel EPI of the brain, and cardiac EPI with a double-oblique slice orientation. Results show that RAC-LORAKS can provide substantial improvements over existing ghost correction methods, and potentially enables new imaging applications that were previously confounded by ghost artifacts.

Volumetric MRI thermometry using a three-dimensional stack-of-stars echo-planar imaging pulse sequence

PURPOSE

To measure temperature over a large brain volume with fine spatiotemporal resolution.

METHODS

A three-dimensional stack-of-stars echo-planar imaging sequence combining echo-planar imaging and radial sampling with golden angle spacing was implemented at 3T for proton resonance frequency-shift temperature imaging. The sequence acquires a 188x188x43 image matrix with 1.5x1.5x2.75 mm3 spatial resolution. Temperature maps were reconstructed using sensitivity encoding (SENSE) image reconstruction followed by the image domain hybrid method, and using the k-space hybrid method. In vivo temperature maps were acquired without heating to measure temperature precision in the brain, and in a phantom during high-intensity focused ultrasound sonication.

RESULTS

In vivo temperature standard deviation was less than 1°C at dynamic scan times down to 0.75 s. For a given frame rate, scanning at a minimum repetition time (TR) with minimum acceleration yielded the lowest standard deviation. With frame rates around 3 s, the scan was tolerant to a small number of receive coils, and temperature standard deviation was 48% higher than a standard two-dimensional Fourier transform temperature mapping scan, but provided whole-brain coverage. Phantom temperature maps with no visible aliasing were produced for dynamic scan times as short as 0.38 s. k-Space hybrid reconstructions were more tolerant to acceleration.

CONCLUSION

Three-dimensional stack-of-stars echo-planar imaging temperature mapping provides volumetric brain coverage and fine spatiotemporal resolution. Magn Reson Med 79:2003-2013, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Comparison of DWI Methods in the Pediatric Brain: PROPELLER Turbo Spin-Echo Imaging Versus Readout-Segmented Echo-Planar Imaging Versus Single-Shot Echo-Planar Imaging

OBJECTIVE

The purpose of this study is to compare DWI for pediatric brain evaluation using single-shot echo-planar imaging (EPI), periodically rotated overlapping parallel lines with enhanced reconstruction (Blade), and readout-segmented EPI (Resolve).

MATERIALS AND METHODS

Blade, Resolve, and single-shot EPI were performed for 27 pediatric patients (median age, 9 years), and three datasets were independently reviewed by two radiologists. Qualitative analyses were performed for perceptive coarseness, image distortion, susceptibility-related changes, motion artifacts, and lesion conspicuity using a 5-point Likert scale. Quantitative analyses were conducted for spatial distortion and signal uniformity of each sequence.

RESULTS

Mean scores were 2.13, 3.17, and 3.76 for perceptive coarseness; 4.85, 3.96, and 2.19 for image distortion; 4.76, 3.96, and 2.30 for susceptibility-related change; 4.96, 3.83, and 4.69 for motion artifacts; and 2.71, 3.75, and 1.92 for lesion conspicuity, for Blade, Resolve, and single-shot EPI, respectively. Blade and Resolve showed better quality than did single-shot EPI for image distortion, susceptibility-related changes, and lesion conspicuity. Blade showed less image distortion, fewer susceptibility-related changes, and fewer motion artifacts than did Resolve, whereas lesion conspicuity was better with Resolve. Blade showed increased signal variation compared with Resolve and single-shot EPI (coefficients of variation were 0.10, 0.08, and 0.05 for lateral ventricle; 0.13, 0.09, and 0.05 for centrum semiovale; and 0.16, 0.09, and 0.06 for pons in Blade, Resolve, and single-shot EPI, respectively).

CONCLUSION

DWI with Resolve or Blade yields better quality regarding distortion, susceptibility-related changes, and lesion conspicuity, compared with single-shot EPI. Blade is less susceptible to motion artifacts than is Resolve, whereas Resolve yields less noise and better lesion conspicuity than does Blade.

Proton change of parotid glands after gustatory stimulation examined by magnetic resonance imaging

The aim of this study was to investigate proton changes of the parotid gland after gustatory stimulation by semi-quantitative parameters and an empirical mathematical model (EMM) using high-temporal-resolution, double-echo, echo-planar imaging (EPI). Approved by a local institutional review board, this study examined 20 parotid glands from 10 healthy volunteers (male:female = 6: 4; age ± standard deviation =35.1 ± 14.1 years) with written informed consent obtained. All participants underwent 1.5-T, double-echo EPI with gustatory stimulation. Semi-quantitative parameters, including maximal drop ratio (MDR), time to peak (TTP), drop slope (DS), recovery slope (RS) and recovery ratio (RR), were calculated. The effect of temporal resolution on parotid functional parameters was evaluated. An EMM comprising an output function ( Sot=Aoe-kot+B) and an input function ( Sint=Ain1-e-kint) was also applied to fit all dynamic curves. Kruskal-Wallis test, Wilcoxon test, linear regression analysis and goodness of fit were used for statistical analysis. p < 0.05 was considered to be statistically significant. The signal intensity dropped significantly after gustatory stimulation on the proton density (PD) image (p < 0.01). MDR was 8.26% in the PD image. MDR and RR were negatively associated with time interval, whereas DS and TTP were significantly positively associated with time interval (all p < 0.05). EMM parametric values derived from PD-time curves of parotid glands were 12.04 ± 6.81%, 6.43 ± 4.23 min^-1 , 88.73 ± 6.18%, 8.41 ± 4.86 min^-1 and 1.09 ± 1.35 for A_o , k_o , B, A_in and k_in , respectively. Semi-quantitative functional parameters and EMM parameters using high-temporal-resolution, double-echo EPI allow the quantification of parotid proton changes after gustatory stimulation.

Connectome analysis for pre-operative brain mapping in neurosurgery

Object: Brain mapping has entered a new era focusing on complex network connectivity. Central to this is the search for the connectome or the brains ‘wiring diagram’. Graph theory analysis of the connectome allows understanding of the importance of regions to network function, and the consequences of their impairment or excision. Our goal was to apply connectome analysis in patients with brain tumours to characterise overall network topology and individual patterns of connectivity alterations.

Methods: Resting-state functional MRI data were acquired using multi-echo, echo planar imaging pre-operatively from five participants each with a right temporal–parietal–occipital glioblastoma. Complex networks analysis was initiated by parcellating the brain into anatomically regions amongst which connections were identified by retaining the most significant correlations between the respective wavelet decomposed time-series.

Results: Key characteristics of complex networks described in healthy controls were preserved in these patients, including ubiquitous small world organization. An exponentially truncated power law fit to the degree distribution predicted findings of general network robustness to injury but with a core of hubs exhibiting disproportionate vulnerability. Tumours produced a consistent reduction in local and long-range connectivity with distinct patterns of connection loss depending on lesion location.

Conclusions: Connectome analysis is a feasible and novel approach to brain mapping in individual patients with brain tumours. Applications to pre-surgical planning include identifying regions critical to network function that should be preserved and visualising connections at risk from tumour resection. In the future one could use such data to model functional plasticity and recovery of cognitive deficits.

The pulsatility volume index: an indicator of cerebrovascular compliance based on fast magnetic resonance imaging of cardiac and respiratory pulsatility

The influence of cardiac activity on the viscoelastic properties of intracranial tissue is one of the mechanisms through which brain-heart interactions take place, and is implicated in cerebrovascular disease. Cerebrovascular disease risk is not fully explained by current risk factors, including arterial compliance. Cerebrovascular compliance is currently estimated indirectly through Doppler sonography and magnetic resonance imaging (MRI) measures of blood velocity changes. In order to meet the need for novel cerebrovascular disease risk factors, we aimed to design and validate an MRI indicator of cerebrovascular compliance based on direct endogenous measures of blood volume changes. We implemented a fast non-gated two-dimensional MRI pulse sequence based on echo-planar imaging (EPI) with ultra-short repetition time (approx. 30-50 ms), which stepped through slices every approximately 20 s. We constrained the solution of the Bloch equations for spins moving faster than a critical speed to produce an endogenous contrast primarily dependent on spin volume changes, and an approximately sixfold signal gain compared with Ernst angle acquisitions achieved by the use of a 90° flip angle. Using cardiac and respiratory peaks detected on physiological recordings, average cardiac and respiratory MRI pulse waveforms in several brain compartments were obtained at 7 Tesla, and used to derive a compliance indicator, the pulsatility volume index (pVI). The pVI, evaluated in larger cerebral arteries, displayed significant variation within and across vessels. Multi-echo EPI showed the presence of significant pulsatility effects in both S0 and [Formula: see text] signals, compatible with blood volume changes. Lastly, the pVI dynamically varied during breath-holding compared with normal breathing, as expected for a compliance indicator. In summary, we characterized and performed an initial validation of a novel MRI indicator of cerebrovascular compliance, which might prove useful to investigate brain-heart interactions in cerebrovascular disease and other disorders.

Single-scan MRI with exceptional resilience to field heterogeneities

PURPOSE

Single-scan two-dimensional MRI has been generally constrained to acquisitions in high quality magnets. This study introduces a methodology, cross-term spatiotemporal encoding (xSPEN), that delivers such images under much poorer external field conditions.

METHODS

xSPEN departs from conventional k-space scanning, by relying on spatiotemporally encoding the image being sought. Unlike hitherto proposed SPEN methods, however, xSPEN's image readout does not take place using a field gradient along the direction being probed, but rather with the aid of an ancillary source of inhomogeneous frequency broadening. This ancillary dimension was here imposed by an orthogonal field gradient; for example, images along the "y" axis were read out by application of a "z" gradient. The principles and characteristics of this new approach, compatible with existing scanners and free from the need to collect auxiliary information such as field maps, are presented and discussed.

RESULTS

Single- and multi-slice in vitro, ex vivo, and in vivo MRI experiments, confirmed the unusual resilience of this new single-shot MRI method to multiple chemical sites on phantoms, animals and humans.

CONCLUSION

xSPEN can deliver single-scan MRI with good sensitivity and exceptional resilience to field inhomogeneities. This could enable investigations that have hitherto escaped from MRI's scope. Magn Reson Med 77:623-634, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Reference-free single-pass EPI Nyquist ghost correction using annihilating filter-based low rank Hankel matrix (ALOHA)

PURPOSE

MR measurements from an echo-planar imaging (EPI) sequence produce Nyquist ghost artifacts that originate from inconsistencies between odd and even echoes. Several reconstruction algorithms have been proposed to reduce such artifacts, but most of these methods require either additional reference scans or multipass EPI acquisition. This article proposes a novel and accurate single-pass EPI ghost artifact correction method that does not require any additional reference data.

THEORY AND METHODS

After converting a ghost correction problem into separate k-space data interpolation problems for even and odd phase encoding, our algorithm exploits an observation that the differential k-space data between the even and odd echoes is a Fourier transform of an underlying sparse image. Accordingly, we can construct a rank-deficient Hankel structured matrix, whose missing data can be recovered using an annihilating filter-based low rank Hankel structured matrix completion approach.

RESULTS

The proposed method was applied to EPI data for both single and multicoil acquisitions. Experimental results using in vivo data confirmed that the proposed method can completely remove ghost artifacts successfully without prescan echoes.

CONCLUSION

Owing to the discovery of the annihilating filter relationship from the intrinsic EPI image property, the proposed method successfully suppresses ghost artifacts without a prescan step. Magn Reson Med 76:1775-1789, 2016. © 2016 International Society for Magnetic Resonance in Medicine.

Accelerated 3D echo-planar imaging with compressed sensing for time-resolved hyperpolarized 13 C studies

PURPOSE

To enable large field-of-view, time-resolved volumetric coverage in hyperpolarized ¹³ C metabolic imaging by implementing a novel data acquisition and image reconstruction method based on the compressed sensing framework.

METHODS

A spectral-spatial pulse for single-resonance excitation followed by a symmetric echo-planar imaging (EPI) readout was implemented for encoding a 72 × 18 cm² field of view at 5 × 5 mm² resolution. Random undersampling was achieved with blipped z-gradients during the ramp portion of the echo-planar imaging readout. The sequence and reconstruction were tested with phantom studies and consecutive in vivo hyperpolarized ¹³ C scans in rats. Retrospectively and prospectively undersampled data were compared on the basis of structural similarity in the reconstructed images and the quantification of the lactate-to-pyruvate ratio in rat kidneys.

RESULTS

No artifacts or loss of resolution are evident in the compressed sensing reconstructed images acquired with the proposed sequence. Structural similarity analysis indicate that compressed sensing reconstructions can accurately recover spatial features in the metabolic images evaluated.

CONCLUSION

A novel z-blip acquisition sequence for compressed sensing accelerated hyperpolarized ¹³ C 3D echo-planar imaging was developed and demonstrated. The close agreement in lactate-to-pyruvate ratios from both retrospectively and prospectively undersampled data from rats shows that metabolic information is preserved with acceleration factors up to 3-fold with the developed method. Magn Reson Med 77:538-546, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Motion immune diffusion imaging using augmented MUSE for high-resolution multi-shot EPI

PURPOSE

To develop new techniques for reducing the effects of microscopic and macroscopic patient motion in diffusion imaging acquired with high-resolution multishot echo-planar imaging.

THEORY

The previously reported multiplexed sensitivity encoding (MUSE) algorithm is extended to account for macroscopic pixel misregistrations, as well as motion-induced phase errors in a technique called augmented MUSE (AMUSE). Furthermore, to obtain more accurate quantitative diffusion-tensor imaging measures in the presence of subject motion, we also account for the altered diffusion encoding among shots arising from macroscopic motion.

METHODS

MUSE and AMUSE were evaluated on simulated and in vivo motion-corrupted multishot diffusion data. Evaluations were made both on the resulting imaging quality and estimated diffusion tensor metrics.

RESULTS

AMUSE was found to reduce image blurring resulting from macroscopic subject motion compared to MUSE but yielded inaccurate tensor estimations when neglecting the altered diffusion encoding. Including the altered diffusion encoding in AMUSE produced better estimations of diffusion tensors.

CONCLUSION

The use of AMUSE allows for improved image quality and diffusion tensor accuracy in the presence of macroscopic subject motion during multishot diffusion imaging. These techniques should facilitate future high-resolution diffusion imaging.

Dual optimization method of radiofrequency and quasistatic field simulations for reduction of eddy currents generated on 7T radiofrequency coil shielding

PURPOSE

To optimize the design of radiofrequency (RF) shielding of transmit coils at 7T and reduce eddy currents generated on the RF shielding when imaging with rapid gradient waveforms.

METHODS

One set of a four-element, 2 × 2 Tic-Tac-Toe head coil structure was selected and constructed to study eddy currents on the RF coil shielding. The generated eddy currents were quantitatively studied in the time and frequency domains. The RF characteristics were studied using the finite difference time domain method. Five different kinds of RF shielding were tested on a 7T MRI scanner with phantoms and in vivo human subjects.

RESULTS

The eddy current simulation method was verified by the measurement results. Eddy currents induced by solid/intact and simple-structured slotted RF shielding significantly distorted the gradient fields. Echo-planar images, B1+ maps, and S matrix measurements verified that the proposed slot pattern suppressed the eddy currents while maintaining the RF characteristics of the transmit coil.

CONCLUSION

The presented dual-optimization method could be used to design RF shielding and reduce the gradient field-induced eddy currents while maintaining the RF characteristics of the transmit coil.

Interleaved diffusion-weighted improved by adaptive partial-Fourier and multiband multiplexed sensitivity-encoding reconstruction

PURPOSE

We report a series of techniques to reliably eliminate artifacts in interleaved echo-planar imaging (EPI) based diffusion-weighted imaging (DWI).

METHODS

First, we integrate the previously reported multiplexed sensitivity encoding (MUSE) algorithm with a new adaptive Homodyne partial-Fourier reconstruction algorithm, so that images reconstructed from interleaved partial-Fourier DWI data are free from artifacts even in the presence of either (a) motion-induced k-space energy peak displacement, or (b) susceptibility field gradient induced fast phase changes. Second, we generalize the previously reported single-band MUSE framework to multiband MUSE, so that both through-plane and in-plane aliasing artifacts in multiband multishot interleaved DWI data can be effectively eliminated.

RESULTS

The new adaptive Homodyne-MUSE reconstruction algorithm reliably produces high-quality and high-resolution DWI, eliminating residual artifacts in images reconstructed with previously reported methods. Furthermore, the generalized MUSE algorithm is compatible with multiband and high-throughput DWI.

CONCLUSION

The integration of the multiband and adaptive Homodyne-MUSE algorithms significantly improves the spatial-resolution, image quality, and scan throughput of interleaved DWI. We expect that the reported reconstruction framework will play an important role in enabling high-resolution DWI for both neuroscience research and clinical uses.

Readout-segmented echo-planar imaging for diffusion-weighted imaging in the pelvis at 3T-A feasibility study

RATIONALE AND OBJECTIVES

Diffusion-weighted imaging (DWI) of the pelvis at 3T is prone to artifacts that diminish the image quality. Readout-segmented echo-planar imaging (RS-EPI) is a new DWI technique that can reduce the artifacts associated with standard single-shot echo-planar imaging (SS-EPI) DWI. The purpose of this study was to evaluate the feasibility and image quality of RS-EPI in pelvic DWI compared to SS-EPI on a 3T imaging system.

MATERIALS AND METHODS

Thirty patients underwent pelvic DWI on a 3T scanner with SS-EPI and RS-EPI techniques. Two blinded readers independently assessed each set of images for geometric distortion, image blurring, ghosting artifacts, lesion conspicuity, and overall image quality on a 7-point scale. Qualitative image scores were compared using paired Wilcoxon signed rank test. Interreader correlation was assessed by Spearman rank correlation.

RESULTS

Geometric distortion, imaging blurring, ghosting artifacts, lesion conspicuity, and overall image quality were rated significantly better by both readers for RS-EPI technique (P < .01 for all parameters). There was moderate-high correlation between the readers (r = 0.649-0.752) for all parameters apart from lesion conspicuity (r = 0.351). Both readers preferred the RS-EPI set of DWI images in most of the cases (reader 1: 0.87, 95% CI 0.74-0.99; reader 2: 0.77, 95% CI 0.61-0.93). Mean difference and limits of agreement between apparent diffusion coefficient (ADC) values obtained from the two methods were 0.01 (-0.08, 0.10) × 10(-3) mm(2)/s.

CONCLUSIONS

RS-EPI DWI images showed improved image quality compared to SS-EPI technique at 3T. RS-EPI is a feasible technique in the pelvis for producing high-resolution DWI.

Frequency correction method for improved spatial correlation of hyperpolarized 13C metabolites and anatomy

Blip-reversed echo-planar imaging (EPI) is investigated as a method for measuring and correcting the spatial shifts that occur due to bulk frequency offsets in (13)C metabolic imaging in vivo. By reversing the k-space trajectory for every other time point, the direction of the spatial shift for a given frequency is reversed. Here, mutual information is used to find the 'best' alignment between images and thereby measure the frequency offset. Time-resolved 3D images of pyruvate/lactate/urea were acquired with 5 s temporal resolution over a 1 min duration in rats (N = 6). For each rat, a second injection was performed with the demodulation frequency purposely mis-set by +35 Hz, to test the correction for erroneous shifts in the images. Overall, the shift induced by the 35 Hz frequency offset was 5.9 ± 0.6 mm (mean ± standard deviation). This agrees well with the expected 5.7 mm shift based on the 2.02 ms delay between k-space lines (giving 30.9 Hz per pixel). The 0.6 mm standard deviation in the correction corresponds to a frequency-detection accuracy of 4 Hz. A method was presented for ensuring the spatial registration between (13)C metabolic images and conventional anatomical images when long echo-planar readouts are used. The frequency correction method was shown to have an accuracy of 4 Hz. Summing the spatially corrected frames gave a signal-to-noise ratio (SNR) improvement factor of 2 or greater, compared with the highest single frame.

Strategies for assessing diffusion anisotropy on the basis of magnetic resonance images: comparison of systematic errors

Diffusion weighted imaging uses the signal loss associated with the random thermal motion of water molecules in the presence of magnetic field gradients to derive a number of parameters that reflect the translational mobility of the water molecules in tissues. With a suitable experimental set-up, it is possible to calculate all the elements of the local diffusion tensor (DT) and derived parameters describing the behavior of the water molecules in each voxel. One of the emerging applications of the information obtained is an interpretation of the diffusion anisotropy in terms of the architecture of the underlying tissue. These interpretations can only be made provided the experimental data which are sufficiently accurate. However, the DT results are susceptible to two systematic error sources: On one hand, the presence of signal noise can lead to artificial divergence of the diffusivities. In contrast, the use of a simplified model for the interaction of the protons with the diffusion weighting and imaging field gradients (b matrix calculation), common in the clinical setting, also leads to deviation in the derived diffusion characteristics. In this paper, we study the importance of these two sources of error on the basis of experimental data obtained on a clinical magnetic resonance imaging system for an isotropic phantom using a state of the art single-shot echo planar imaging sequence. Our results show that optimal diffusion imaging require combining a correct calculation of the b-matrix and a sufficiently large signal to noise ratio.

Improved B0 -distortion correction in diffusion MRI using interlaced q-space sampling and constrained reconstruction

PURPOSE

To enable high-quality correction of susceptibility-induced geometric distortion artifacts in diffusion magnetic resonance imaging (MRI) images without increasing scan time.

THEORY AND METHODS

A new method for distortion correction is proposed based on subsampling a generalized version of the state-of-the-art reversed-gradient distortion correction method. Rather than acquire each q-space sample multiple times with different distortions (as in the conventional reversed-gradient method), we sample each q-space point once with an interlaced sampling scheme that measures different distortions at different q-space locations. Distortion correction is achieved using a novel constrained reconstruction formulation that leverages the smoothness of diffusion data in q-space.

RESULTS

The effectiveness of the proposed method is demonstrated with simulated and in vivo diffusion MRI data. The proposed method is substantially faster than the reversed-gradient method, and can also provide smaller intensity errors in the corrected images and smaller errors in derived quantitative diffusion parameters.

CONCLUSION

The proposed method enables state-of-the-art distortion correction performance without increasing data acquisition time.

Is there a cerebellar compensatory effort in first-episode, treatment-naive major depressive disorder at rest?

BACKGROUND

This study was undertaken to explore whether there is a cerebellar compensatory response in patients with first-episode, treatment-naive major depressive disorder (MDD). The cerebellar compensatory response is defined as a cerebellar hyperactivity which would be inversely correlated with both the activation of the functionally connected cerebral regions and the depression severity.

METHODS

Resting-state functional magnetic resonance imaging (fMRI) data of 24 patients with MDD and 24 healthy subjects were analyzed with the fractional amplitude of low-frequency fluctuations (fALFF) and functional connectivity (FC) methods. The structural images were processed with the voxel-based morphometry (VBM) method.

RESULTS

Compared to healthy controls, depressed patients had significantly increased fALFF in the left Crus I and the left cerebellar lobule VI. FC analysis of these two seeded regions found that depressed patients had increased FC between the left Crus I and the right hippocampus, but had decreased FC between the left Crus I and the left inferior parietal lobule (IPL), and between the left cerebellar lobule VI and bilateral inferior temporal gyrus. No correlation was observed between the abnormal fALFF of the seeds and their connected regions and the depression severity or the executive function. The VBM results did not show significant reduction in gray or white matter volume in any above-mentioned region.

CONCLUSIONS

Our findings suggest that increased cerebellar activity at resting state may be a disease state phenomenon but not a compensatory response to the dysfunction of the default mode network (DMN) in MDD.

Direct design of 2D RF pulses using matrix inversion

Multi-dimensional pulses are frequently used in MRI for applications such as targeted excitation, fat-water separation or metabolic imaging with hyperpolarised (13)C compounds. For the design, the problem is typically separated into the different dimensions. In this work, a method to directly design two-dimensional pulses using the small-tip angle approximation is introduced based on a direct matrix representation. The numerical problem is solved in a single step directly in two dimensions by matrix inversion. Exemplary spectral-spatial excitation and spatio-temporal encoding (SPEN) pulses are designed and validated. The main benefits of the direct design approach include a reduction of artefacts in case of spectral-spatial pulses, a simple and straightforward computer implementation and high flexibility in the pulse design.

Quasi-periodic patterns (QPP): large-scale dynamics in resting state fMRI that correlate with local infraslow electrical activity

Functional connectivity measurements from resting state blood-oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) are proving a powerful tool to probe both normal brain function and neuropsychiatric disorders. However, the neural mechanisms that coordinate these large networks are poorly understood, particularly in the context of the growing interest in network dynamics. Recent work in anesthetized rats has shown that the spontaneous BOLD fluctuations are tightly linked to infraslow local field potentials (LFPs) that are seldom recorded but comparable in frequency to the slow BOLD fluctuations. These findings support the hypothesis that long-range coordination involves low frequency neural oscillations and establishes infraslow LFPs as an excellent candidate for probing the neural underpinnings of the BOLD spatiotemporal patterns observed in both rats and humans. To further examine the link between large-scale network dynamics and infraslow LFPs, simultaneous fMRI and microelectrode recording were performed in anesthetized rats. Using an optimized filter to isolate shared components of the signals, we found that time-lagged correlation between infraslow LFPs and BOLD is comparable in spatial extent and timing to a quasi-periodic pattern (QPP) found from BOLD alone, suggesting that fMRI-measured QPPs and the infraslow LFPs share a common mechanism. As fMRI allows spatial resolution and whole brain coverage not available with electroencephalography, QPPs can be used to better understand the role of infraslow oscillations in normal brain function and neurological or psychiatric disorders.

Investigating univariate temporal patterns for intrinsic connectivity networks based on complexity and low-frequency oscillation: a test-retest reliability study

Intrinsic connectivity networks (ICNs) are composed of spatial components and time courses. The spatial components of ICNs were discovered with moderate-to-high reliability. So far as we know, few studies focused on the reliability of the temporal patterns for ICNs based their individual time courses. The goals of this study were twofold: to investigate the test-retest reliability of temporal patterns for ICNs, and to analyze these informative univariate metrics. Additionally, a correlation analysis was performed to enhance interpretability. Our study included three datasets: (a) short- and long-term scans, (b) multi-band echo-planar imaging (mEPI), and (c) eyes open or closed. Using dual regression, we obtained the time courses of ICNs for each subject. To produce temporal patterns for ICNs, we applied two categories of univariate metrics: network-wise complexity and network-wise low-frequency oscillation. Furthermore, we validated the test-retest reliability for each metric. The network-wise temporal patterns for most ICNs (especially for default mode network, DMN) exhibited moderate-to-high reliability and reproducibility under different scan conditions. Network-wise complexity for DMN exhibited fair reliability (ICC<0.5) based on eyes-closed sessions. Specially, our results supported that mEPI could be a useful method with high reliability and reproducibility. In addition, these temporal patterns were with physiological meanings, and certain temporal patterns were correlated to the node strength of the corresponding ICN. Overall, network-wise temporal patterns of ICNs were reliable and informative and could be complementary to spatial patterns of ICNs for further study.

The Importance of k-Space Trajectory on Off-Resonance Artifact in Segmented Echo-Planar Imaging

Segmented interleaved echo planar imaging (EPI) is a highly efficient data acquisition technique; however, EPI is sensitive to artifacts from off-resonance spins. The choice of k-space trajectories is important in determining how off-resonance spins contribute to image artifacts. Top-down and center-out trajectories are theoretically analyzed, simulated, implemented, and tested in phantom and volunteer experiments. Theoretical results show off-resonance artifact manifests as a simple positional shift for the top-down trajectory, while for the center-out trajectory off-resonance artifact manifests as a splitting of the object, which entails both shift and blurring. These results were validated using simulation and phantom scan data where a frequency-offset was introduced ranging from -300 Hz to +300 Hz. As predicted by the theoretical results, inferior image quality was observed for the center-out trajectory in a single volunteer. Off-resonance produces more severe and complex artifacts with the center-out trajectory than the top-down trajectory.

Robust GRAPPA-accelerated diffusion-weighted readout-segmented (RS)-EPI

Readout segmentation (RS-EPI) has been suggested as a promising variant to echo-planar imaging (EPI) for high-resolution imaging, particularly when combined with parallel imaging. This work details some of the technical aspects of diffusion-weighted (DW)-RS-EPI, outlining a set of reconstruction methods and imaging parameters that can both minimize the scan time and afford high-resolution diffusion imaging with reduced distortions. These methods include an efficient generalized autocalibrating partially parallel acquisition (GRAPPA) calibration for DW-RS-EPI data without scan time penalty, together with a variant for the phase correction of partial Fourier RS-EPI data. In addition, the role of pulsatile and rigid-body brain motion in DW-RS-EPI was assessed. Corrupt DW-RS-EPI data arising from pulsatile nonlinear brain motion had a prevalence of approximately 7% and were robustly identified via k-space entropy metrics. For DW-RS-EPI data corrupted by rigid-body motion, we showed that no blind overlap was required. The robustness of RS-EPI toward phase errors and motion, together with its minimized distortions compared with EPI, enables the acquisition of exquisite 3 T DW images with matrix sizes close to 512(2).