Image Registration to Compensate for EPI Distortion in Patients with Brain Tumors: An Evaluation of Tract-Specific Effects.

BACKGROUND AND PURPOSE: Diffusion magnetic resonance imaging (dMRI) provides preoperative maps of neurosurgical patients' white matter tracts, but these maps suffer from echo-planar imaging (EPI) distortions caused by magnetic field inhomogeneities. In clinical neurosurgical planning, these distortions are generally not corrected and thus contribute to the uncertainty of fiber tracking. Multiple image processing pipelines have been proposed for image-registration-based EPI distortion correction in healthy subjects. In this article, we perform the first comparison of such pipelines in neurosurgical patient data. METHODS: Five pipelines were tested in a retrospective clinical dMRI dataset of 9 patients with brain tumors. Pipelines differed in the choice of fixed and moving images and the similarity metric for image registration. Distortions were measured in two important tracts for neurosurgery, the arcuate fasciculus and corticospinal tracts. RESULTS: Significant differences in distortion estimates were found across processing pipelines. The most successful pipeline used dMRI baseline and T2-weighted images as inputs for distortion correction. This pipeline gave the most consistent distortion estimates across image resolutions and brain hemispheres. CONCLUSIONS: Quantitative results of mean tract distortions on the order of 1-2 mm are in line with other recent studies, supporting the potential need for distortion correction in neurosurgical planning. Novel results include significantly higher distortion estimates in the tumor hemisphere and greater effect of image resolution choice on results in the tumor hemisphere. Overall, this study demonstrates possible pitfalls and indicates that care should be taken when implementing EPI distortion correction in clinical settings.

Dual-Pathway sequences for MR thermometry: When and where to use them.

PURPOSE: Dual-pathway sequences have been proposed to help improve the temperature-to-noise ratio (TNR) in MR thermometry. The present work establishes how much of an improvement these so-called "PSIF-FISP" sequences may bring in various organs and tissues. METHODS: Simulations and TNR calculations were validated against analytical equations, phantom, abdomen, and brain scans. Relative TNRs for PSIF-FISP, as compared to a dual-FISP reference standard, were calculated for flip angle (FA) = 1 to 85 º and repetition time (TR) = 6 to 60 ms, for gray matter, white matter, cervix, endometrium, myometrium, prostate, kidney medulla and cortex, bone marrow, pancreas, spleen, muscle, and liver tissues. RESULTS: PSIF-FISP was TNR superior in the kidney, pelvis, spleen, or gray matter at most tested TR and FA settings, and benefits increased at shorter TRs. PSIF-FISP was TNR superior in other tissues, e.g., liver, muscle, pancreas, for only short TR settings (20 ms or less). The TNR benefits of PSIF-FISP increased slightly with FA, and strongly with decreasing TR. Up to two- to three-fold reductions in TR with 20% TNR gains were achievable. In any given tissue, TNR performance is expected to further improve with heating, due to changes in relaxation rates. CONCLUSION: Dual-pathway PSIF-FISP can improve TNR and acquisition speed over standard gradient-recalled echo sequences, but optimal acquisition parameters are tissue dependent. Magn Reson Med 77:1193-1200, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Automated white matter fiber tract identification in patients with brain tumors.

We propose a method for the automated identification of key white matter fiber tracts for neurosurgical planning, and we apply the method in a retrospective study of 18 consecutive neurosurgical patients with brain tumors. Our method is designed to be relatively robust to challenges in neurosurgical tractography, which include peritumoral edema, displacement, and mass effect caused by mass lesions. The proposed method has two parts. First, we learn a data-driven white matter parcellation or fiber cluster atlas using groupwise registration and spectral clustering of multi-fiber tractography from healthy controls. Key fiber tract clusters are identified in the atlas. Next, patient-specific fiber tracts are automatically identified using tractography-based registration to the atlas and spectral embedding of patient tractography. Results indicate good generalization of the data-driven atlas to patients: 80% of the 800 fiber clusters were identified in all 18 patients, and 94% of the 800 fiber clusters were found in 16 or more of the 18 patients. Automated subject-specific tract identification was evaluated by quantitative comparison to subject-specific motor and language functional MRI, focusing on the arcuate fasciculus (language) and corticospinal tracts (motor), which were identified in all patients. Results indicate good colocalization: 89 of 95, or 94%, of patient-specific language and motor activations were intersected by the corresponding identified tract. All patient-specific activations were within 3mm of the corresponding language or motor tract. Overall, our results indicate the potential of an automated method for identifying fiber tracts of interest for neurosurgical planning, even in patients with mass lesions.

Characterization of gradient echo signal decays in healthy and cancerous prostate at 3T improves with a Gaussian augmentation of the mono-exponential (GAME) model.

A biomarker of cancer aggressiveness, such as hypoxia, could substantially impact treatment decisions in the prostate, especially radiation therapy, by balancing treatment morbidity (urinary incontinence, erectile dysfunction, etc.) against mortality. R2 (*) mapping with Mono-Exponential (ME) decay modeling has shown potential for identifying areas of prostate cancer hypoxia at 1.5T. However, Gaussian deviations from ME decay have been observed in other tissues at 3T. The purpose of this study is to assess whether gradient-echo signal decays are better characterized by a standard ME decay model, or a Gaussian Augmentation of the Mono-Exponential (GAME) decay model, in the prostate at 3T. Multi-gradient-echo signals were acquired on 20 consecutive patients with a clinical suspicion of prostate cancer undergoing MR-guided prostate biopsies. Data were fitted with both ME and GAME models. The information contents of these models were compared using Akaike's information criterion (second order, AICC ), in skeletal muscle, the prostate central gland (CG), and peripheral zone (PZ) regions of interest (ROIs). The GAME model had higher information content in 30% of the prostate on average (across all patients and ROIs), covering up to 67% of cancerous PZ ROIs, and up to 100% of cancerous CG ROIs (in individual patients). The higher information content of GAME became more prominent in regions that would be assumed hypoxic using ME alone, reaching 50% of the PZ and 70% of the CG as ME R2 (*) approached 40 s(-1) . R2 (*) mapping may have important applications in MRI; however, information lost due to modeling could mask differences in parameters due to underlying tissue anatomy or physiology. The GAME model improves characterization of signal behavior in the prostate at 3T, and may increase the potential for determining correlates of fit parameters with biomarkers, for example of oxygenation status.

Navigated DENSE strain imaging for post-radiofrequency ablation lesion assessment in the swine left atria.

AIMS: Prior work has demonstrated that magnetic resonance imaging (MRI) strain can separate necrotic/stunned myocardium from healthy myocardium in the left ventricle (LV). We surmised that high-resolution MRI strain, using navigator-echo-triggered DENSE, could differentiate radiofrequency ablated tissue around the pulmonary vein (PV) from tissue that had not been damaged by radiofrequency energy, similarly to navigated 3D myocardial delayed enhancement (3D-MDE). METHODS AND RESULTS: A respiratory-navigated 2D-DENSE sequence was developed, providing strain encoding in two spatial directions with 1.2 × 1.0 × 4 mm(3) resolution. It was tested in the LV of infarcted sheep. In four swine, incomplete circumferential lesions were created around the right superior pulmonary vein (RSPV) using ablation catheters, recorded with electro-anatomic mapping, and imaged 1 h later using atrial-diastolic DENSE and 3D-MDE at the left atrium/RSPV junction. DENSE detected ablation gaps (regions with >12% strain) in similar positions to 3D-MDE (2D cross-correlation 0.89 ± 0.05). Low-strain (<8%) areas were, on average, 33% larger than equivalent MDE regions, so they include both injured and necrotic regions. Optimal DENSE orientation was perpendicular to the PV trunk, with high shear strain in adjacent viable tissue appearing as a sensitive marker of ablation lesions. CONCLUSIONS: Magnetic resonance imaging strain may be a non-contrast alternative to 3D-MDE in intra-procedural monitoring of atrial ablation lesions.

Non-invasive quantification of absolute cerebral blood volume during functional activation applicable to the whole human brain.

PURPOSE: Cerebral blood volume (CBV) changes in many diverse pathologic conditions, and in response to functional challenges along with changes in blood flow, blood oxygenation, and the cerebral metabolic rate of oxygen. The feasibility of a new method for non-invasive quantification of absolute cerebral blood volume that can be applicable to the whole human brain was investigated. METHODS: Multi-slice data were acquired at 3 T using a novel inversion recovery echo planar imaging (IR-EPI) pulse sequence with varying contrast weightings and an efficient rotating slice acquisition order, at rest and during visual activation. A biophysical model was used to estimate absolute cerebral blood volume at rest and during activation, and oxygenation during activation, on data from 13 normal human subjects. RESULTS: Cerebral blood volume increased by 21.7% from 6.6 ± 0.8 mL/100 mL of brain parenchyma at rest to 8.0 ± 1.3 mL/100 mL of brain parenchyma in the occipital cortex during visual activation, with average blood oxygenation of 84 ± 2.1% during activation, comparing well with literature. CONCLUSION: The method is feasible, and could foster improved understanding of the fundamental physiological relationship between neuronal activity, hemodynamic changes, and metabolism underlying brain activation; complement existing methods for estimating compartmental changes; and potentially find utility in evaluating vascular health.

Noninvasive MRI measurement of the absolute cerebral blood volume-cerebral blood flow relationship during visual stimulation in healthy humans.

PURPOSE: The relationship between cerebral blood volume (CBV) and cerebral blood flow (CBF) underlies blood oxygenation level-dependent functional MRI signal. This study investigates the potential for improved characterization of the CBV-CBF relationship in humans, and examines sex effects as well as spatial variations in the CBV-CBF relationship. METHODS: Healthy subjects were imaged noninvasively at rest and during visual stimulation, constituting the first MRI measurement of the absolute CBV-CBF relationship in humans with complete coverage of the functional areas of interest. RESULTS: CBV and CBF estimates were consistent with the literature, and their relationship varied both spatially and with sex. In a region of interest with stimulus-induced activation in CBV and CBF at a significance level of the P < 0.05, a power function fit resulted in CBV = 2.1 CBF(0.32) across all subjects, CBV = 0.8 CBF(0.51) in females and CBV = 4.4 CBF(0.15) in males. Exponents decreased in both sexes as ROIs were expanded to include less significantly activated regions. CONCLUSION: Consideration for potential sex-related differences, as well as regional variations under a range of physiological states, may reconcile some of the variation across literature and advance our understanding of the underlying cerebrovascular physiology.

O-space imaging: Highly efficient parallel imaging using second-order nonlinear fields as encoding gradients with no phase encoding.

Recent improvements in parallel imaging have been driven by the use of greater numbers of independent surface coils placed so as to minimize aliasing along the phase-encode direction(s). However, gains from increasing the number of coils diminish as coil coupling problems begin to dominate and the ratio of acceleration gain to expense for multiple receiver chains becomes prohibitive. In this work, we redesign the spatial-encoding strategy in order to gain efficiency, achieving a gradient encoding scheme that is complementary to the spatial encoding provided by the receiver coils. This approach leads to "O-space" imaging, wherein the gradient shapes are tailored to an existing surface coil array, making more efficient use of the spatial information contained in the coil profiles. In its simplest form, for each acquired echo the Z2 spherical harmonic is used to project the object onto sets of concentric rings, while the X and Y gradients are used to offset this projection within the imaging plane. The theory is presented, an algorithm is introduced for image reconstruction, and simulations reveal that O-space encoding achieves high encoding efficiency compared to sensitivity encoding (SENSE) radial k-space trajectories, and parallel imaging technique with localized gradients (PatLoc), suggesting that O-space imaging holds great potential for accelerated scanning.