Kestutis Barkauskas

Non-Cartesian Data Reconstruction Using GRAPPA Operator Gridding (GROG)

A novel approach that uses the concepts of parallel imaging to grid data sampled along a non‐Cartesian trajectory using GRAPPA operator gridding (GROG) is described. GROG shifts any acquired data point to its nearest Cartesian location, thereby converting non‐Cartesian to Cartesian data. Unlike other parallel imaging methods, GROG synthesizes the net weight for a shift in any direction from a single basis set of weights along the logical k‐space directions. Given the vastly reduced size of the basis set, GROG calibration and reconstruction requires fewer operations and less calibration data than other parallel imaging methods for gridding. Instead of calculating and applying a density compensation function (DCF), GROG requires only local averaging, as the reconstructed points fall upon the Cartesian grid. Simulations are performed to demonstrate that the root mean square error (RMSE) values of images gridded with GROG are similar to those for images gridded using the gold‐standard convolution gridding. Finally, GROG is compared to the convolution gridding technique using data sampled along radial, spiral, rosette, and BLADE (a.k.a. periodically rotated overlapping parallel lines with enhanced reconstruction [PROPELLER]) trajectories. Magn Reson Med, 2007.

Variation correction algorithm: Analysis of phase suppression and thermal profile fidelity for proton resonance frequency magnetic resonance thermometry at 0.2 T

To develop and analyze the performance of the variation correction algorithm (VCA), a phase correction technique that mitigates the contribution of background phase variations by combining accurate alignment of echoes, K‐space‐based phase correction (as opposed to spatial polynomials), and extraction of alias‐free phase difference images.

Focal transient CNS vessel leak provides a tissue niche for sequential immune cell accumulation during the asymptomatic phase of EAE induction

Peripheral immune cells are critical to the pathogenesis of neurodegenerative diseases including multiple sclerosis (MS) (Hendriks et al., 2005; Kasper and Shoemaker, 2010). However, the precise sequence of tissue events during the early asymptomatic induction phase of experimental autoimmune encephalomyelitis (EAE) pathogenesis remains poorly defined. Due to the spatial-temporal constrains of traditional methods used to study this disease, most studies had been performed in the spine during peak clinical disease; thus the debate continues as to whether tissue changes such as vessel disruption represent a cause or a byproduct of EAE pathophysiology in the cortex. Here, we provide dynamic, high-resolution information on the evolving structural and cellular processes within the gray matter of the mouse cortex during the first 12 asymptomatic days of EAE induction. We observed that transient focal vessel disruptions precede microglia activation, followed by infiltration of and directed interaction between circulating dendritic cells and T cells. Histamine antagonist minimizes but not completely ameliorates blood vessel leaks. Histamine H1 receptor blockade prevents early microglia function, resulting in subsequent reduction in immune cell accumulation, disease incidence and clinical severity.

Magnetic resonance imaging and model prediction for thermal ablation of tissue

To monitor and predict tissue temperature distributions and lesion boundaries during thermal ablation by combining MRI and thermal modeling methods.

Quantification of left ventricular functional parameter values using 3D spiral bSSFP and through-time Non-Cartesian GRAPPA

BackgroundThe standard clinical acquisition for left ventricular functional parameter analysis with cardiovascular magnetic resonance (CMR) uses a multi-breathhold multi-slice segmented balanced SSFP sequence. Performing multiple long breathholds in quick succession for ventricular coverage in the short-axis orientation can lead to fatigue and is challenging in patients with severe cardiac or respiratory disorders. This study combines the encoding efficiency of a six-fold undersampled 3D stack of spirals balanced SSFP sequence with 3D through-time spiral GRAPPA parallel imaging reconstruction. This 3D spiral method requires only one breathhold to collect the dynamic data.MethodsTen healthy volunteers were recruited for imaging at 3 T. The 3D spiral technique was compared against 2D imaging in terms of systolic left ventricular functional parameter values (Bland-Altman plots), total scan time (Welch’s t-test) and qualitative image rating scores (Wilcoxon signed-rank test).ResultsSystolic left ventricular functional values were not significantly different (i.e. 3D-2D) between the methods. The 95% confidence interval for ejection fraction was −0.1 ± 1.6% (mean ± 1.96*SD). The total scan time for the 3D spiral technique was 48 s, which included one breathhold with an average duration of 14’s for the dynamic scan, plus 34’s to collect the calibration data under free-breathing conditions. The 2D method required an average of 5min40s for the same coverage of the left ventricle. The difference between 3D and 2D image rating scores was significantly different from zero (Wilcoxon signed-rank test, p < 0.05); however, the scores were at least 3 (i.e. average) or higher for 3D spiral imaging.ConclusionThe 3D through-time spiral GRAPPA method demonstrated equivalent systolic left ventricular functional parameter values, required significantly less total scan time and yielded acceptable image quality with respect to the 2D segmented multi-breathhold standard in this study. Moreover, the 3D spiral technique used just one breathhold for dynamic imaging, which is anticipated to reduce patient fatigue as part of the complete cardiac examination in future studies that include patients.

Dynamics of MRI-Guided Thermal Ablation of VX2 Tumor in Paraspinal Muscle of Rabbits

This study combines fast magnetic resonance imaging (MRI) and model simulation of tissue thermal ablation for monitoring and predicting the dynamics of lesion size for tumor destruction. In vivo experiments were conducted using radiofrequency (RF) thermal ablation in paraspinal muscle of rabbit with a VX2 tumor. Before ablation, turbo-spin echo (TSE) images visualized the 3-D tumor (necrotic core and tumor periphery) and surrounding normal tissue. MR gradient-recalled echo (GRE) phase and magnitude images were acquired repeatedly in 3.3 s at 30-s intervals during and after thermal ablation to follow tissue temperature distribution dynamics and lesion development in tumor and surrounding normal tissue. Final lesion sizes estimated from GRE magnitude, post-ablation TSE, and stained histologic images were compared. Model simulations of temperature distribution and lesion development dynamics closely corresponded to the experimental data from MR images in tumor and normal tissue. The combined use of MR image monitoring and model simulation has the potential for improving pretreatment planning and real-time prediction of lesion-size dynamics for guidance of thermal ablation of tumors.

Accelerated 2D multi-slice first-pass contrast-enhanced myocardial perfusion using through-time radial GRAPPA

Background Non-Cartesian parallel imaging is a promising approach for reducing the scan time in multi-slice first-pass myocardial perfusion imaging, allowing increasing volumetric coverage in comparison with standard techniques. Through-time radial GRAPPA has previously been demonstrated for real-time functional cardiac imaging [Seiberlich, et al. MRM 2011 Feb;65 (2):492-505]. In this work, the through-time radial GRAPPA technique is applied for the acquisition of fifteen slices per heartbeat during a contrast-enhanced myocardial perfusion examination.

Single breath-hold 3D cardiac T 1 mapping using through-time spiral GRAPPA

The quantification of cardiac T1 relaxation time holds great potential for the detection of various cardiac diseases. However, as a result of both cardiac and respiratory motion, only one two‐dimensional T1 map can be acquired in one breath‐hold with most current techniques, which limits its application for whole heart evaluation in routine clinical practice. In this study, an electrocardiogram (ECG)‐triggered three‐dimensional Look–Locker method was developed for cardiac T1 measurement. Fast three‐dimensional data acquisition was achieved with a spoiled gradient‐echo sequence in combination with a stack‐of‐spirals trajectory and through‐time non‐Cartesian generalized autocalibrating partially parallel acquisition (GRAPPA) acceleration. The effects of different magnetic resonance parameters on T1 quantification with the proposed technique were first examined by simulating data acquisition and T1 map reconstruction using Bloch equation simulations. Accuracy was evaluated in studies with both phantoms and healthy subjects. These results showed that there was close agreement between the proposed technique and the reference method for a large range of T1 values in phantom experiments. In vivo studies further demonstrated that rapid cardiac T1 mapping for 12 three‐dimensional partitions (spatial resolution, 2 × 2 × 8 mm3) could be achieved in a single breath‐hold of ~12 s. The mean T1 values of myocardial tissue and blood obtained from normal volunteers at 3 T were 1311 ± 66 and 1890 ± 159 ms, respectively. In conclusion, a three‐dimensional T1 mapping technique was developed using a non‐Cartesian parallel imaging method, which enables fast and accurate T1 mapping of cardiac tissues in a single short breath‐hold.