Hybrid referenceless and multibaseline subtraction MR thermometry for monitoring thermal therapies in moving organs

Technical advances in motion-robust MR thermometry

Proton resonance frequency shift (PRFS) MR thermometry is the most common method used in clinical thermal treatments because of its fast acquisition and high sensitivity to temperature. However, motion is the biggest obstacle in PRFS MR thermometry for monitoring thermal treatment in moving organs. This challenge arises because of the introduction of phase errors into the PRFS calculation through multiple methods, such as image misregistration, susceptibility changes in the magnetic field, and intraframe motion during MRI acquisition. Various approaches for motion correction have been developed for real-time, motion-robust, and volumetric MR thermometry. However, current technologies have inherent trade-offs among volume coverage, processing time, and temperature accuracy. These tradeoffs should be considered and chosen according to the thermal treatment application. In hyperthermia treatment, precise temperature measurements are of increased importance rather than the requirement for exceedingly high temporal resolution. In contrast, ablation procedures require robust temporal resolution to accurately capture a rapid temperature rise. This paper presents a comprehensive review of current cutting-edge MRI techniques for motion-robust MR thermometry, and recommends which techniques are better suited for each thermal treatment. We expect that this study will help discern the selection of motion-robust MR thermometry strategies and inspire the development of motion-robust volumetric MR thermometry for practical use in clinics.

A Non-Contrast Multi-Parametric MRI Biomarker for Assessment of MR-Guided Focused Ultrasound Thermal Therapies

OBJECTIVE

We present the development of a non-contrast multi-parametric magnetic resonance (MPMR) imaging biomarker to assess treatment outcomes for magnetic resonance-guided focused ultrasound (MRgFUS) ablations of localized tumors. Images obtained immediately following MRgFUS ablation were inputs for voxel-wise supervised learning classifiers, trained using registered histology as a label for thermal necrosis.

METHODS

VX2 tumors in New Zealand white rabbits quadriceps were thermally ablated using an MRgFUS system under 3 T MRI guidance. Animals were re-imaged three days post-ablation and euthanized. Histological necrosis labels were created by 3D registration between MR images and digitized H&E segmentations of thermal necrosis to enable voxel-wise classification of necrosis. Supervised MPMR classifier inputs included maximum temperature rise, cumulative thermal dose (CTD), post-FUS differences in T2-weighted images, and apparent diffusion coefficient, or ADC, maps. A logistic regression, support vector machine, and random forest classifier were trained in red a leave-one-out strategy in test data from four subjects.

RESULTS

In the validation dataset, the MPMR classifiers achieved higher recall and Dice than a clinically adopted 240 cumulative equivalent minutes at 43 °C (CEM 43) threshold (0.43) in all subjects. The average Dice scores of overlap with the registered histological label for the logistic regression (0.63) and support vector machine (0.63) MPMR classifiers were within 6% of the acute contrast-enhanced non-perfused volume (0.67).

CONCLUSIONS

Voxel-wise registration of MPMR data to histological outcomes facilitated supervised learning of an accurate non-contrast MR biomarker for MRgFUS ablations in a rabbit VX2 tumor model.

MR Thermometry during Transcranial MR Imaging-Guided Focused Ultrasound Procedures: A Review

Interest in transcranial MR imaging-guided focused ultrasound procedures has recently grown. These incisionless procedures enable precise focal ablation of brain tissue using real-time monitoring by MR thermometry. This article will provide an updated review on clinically applicable technical underpinnings and considerations of proton resonance frequency MR thermometry, the most common clinically used MR thermometry sequence.

Magnetic resonance fingerprinting based thermometry (MRFT): application toex vivoimaging near DBS leads

The purpose of this study is to demonstrate the first work ofT1-based magnetic resonance thermometry using magnetic resonance fingerprinting (dubbed MRFT). We compared temperature estimation of MRFT with proton resonance frequency shift (PRFS) thermometry onex vivobovine muscle. We demonstrated MRFT's feasibility in predicting temperature onex vivobovine muscles with deep brain stimulation (DBS) lead.B0maps generated from MRFT were compared with gold standardB0maps near the DBS lead. MRFT and PRFS estimated temperatures were compared in the presence of motion. All experiments were performed on a 3 Tesla whole-body GE Premier system with a 21-channel receive head coil (GE Healthcare, Milwaukee, WI). Four fluoroptic probes were used to measure the temperature at the center of a cold muscle (probe 1), the room temperature water bottle (probe 2), and the center and periphery of the heated muscle (probes 3 and 4). We selected regions of interest (ROIs) around the location of the probes and used simple linear regression to generate the temperature sensitivity calibration equations that convertT1maps and Δsmaps to temperature maps. We then repeated the same setup and compared MRFT and PRFS thermometry temperature estimation with gold standard probe measurements. For the MRFT experiment on DBS lead, we taped the probe to the tip of the DBS lead and used a turbo spin echo sequence to induce heating near the lead. We selected ROIs around the tip of the lead to compare MRFT temperature estimation with probe measurements and compared with PRFS temperature estimation. Vendor-suppliedB0mapping sequence was acquired to compare with MRFT-generatedB0maps. We found strong linear relationships (R2> 0.958) betweenT1and temperature and Δsand temperatures in our temperature sensitivity calibration experiment. MRFT and PRFS thermometry both accurately predict temperature (RMSE < 1.55 °C) compared to probe measurements. MRFT estimated temperature near DBS lead has a similar trend as the probe temperature. BothB0maps show inhomogeneities around the lead. MRFT estimated temperature is less sensitive to motion.

Appearance and modeling of bubble artifacts in intracranial magnetic resonance-guided laser interstitial thermal therapy (MRg-LITT) temperature images

PURPOSE

To understand if unexplained signal artifacts in MRg-LITT proton resonance frequency- (PRF-) shift thermometry images are caused by air bubbles or hemorrhages, and to characterize their effects on temperature measurements.

METHODS

Retrospective image data from an IRB-approved clinical trial of intracranial MRg-LITT were inspected for asymmetric distortions observed in phase data during ablations, which have been previously reported as likely hemorrhages. A total of eight patient cases were selected: seven with artifact occurrence and one without. Mathematical image models for air bubbles or hemorrhages were implemented to estimate the size of the air bubble or hemorrhage needed to explain the clinically observed phase artifacts. Correlations and Bland-Altman analyses were used to determine if an air bubble model or a hemorrhage model was better correlated to the clinical data. The model was used to inject bubbles into clean PRF phase data without artifacts to examine how temperature profile distortions change with slice orientation. The simulated air-bubble injected data were compared to clinical data containing artifacts to examine the bubbles' effects on temperature and thermal damage estimates.

RESULTS

The model demonstrated that air bubbles up to approximately 1 cm in diameter could explain the clinically observed phase artifacts. The bubble model predicts that a hemorrhage would have to be 2.2 times as large as an air bubble in order to explain the same extent of phase distortion observed in clinical data. Air bubbles had 16% percent higher correlations to the clinical PRF phase data than hemorrhages, even after rescaling the hemorrhage phases to better match the data. The air bubble model also explains how the phase artifacts lead to both large positive and large negative temperature errors, up to ±100 °C, which could cascade to damage estimate errors of several millimeters.

CONCLUSION

Results showed that the artifacts are likely caused by air bubbles rather than hemorrhages, which may be introduced before heating or appear during heating. Manufacturers and users of devices that rely upon PRF-shift thermometry should be aware these phase distortions from bubble artifacts can result in large temperature errors.

Detection of laser-associated heating in the brain during simultaneous fMRI and optogenetic stimulation

PURPOSE

To calculate temperatures from T2 *-weighted images collected during optogenetic fMRI based on proton resonance frequency (PRF) shift thermometry, to monitor confounding heating effects and determine appropriate light parameters for optogenetic stimulation.

METHODS

fMRI is mainly based on long-TE gradient-recalled echo acquisitions that are also suitable for measuring small temperature changes via the PRF shift. A motion- and respiration-robust processing pipeline was developed to calculate temperature changes based on the PRF shift directly from the T2 *-weighted images collected for fMRI with a two-shot 2D gradient-recalled echo-EPI sequence at 9.4T. Optogenetic fMRI protocols which differed in stimulation durations (3, 6 and 9 s) within a total block duration of 30 s were applied in a squirrel monkey to validate the methods with blue and green light (20 Hz, 30 mW) delivery interleaved between periods. General linear modeling was performed on the resulting time series temperature maps to verify if light delivery with each protocol resulted in significant heating in the brain around the optical fiber.

RESULTS

The temperature SD was 0.05°C with the proposed imaging protocol and processing. Statistical analysis showed that the optogenetic stimulation protocol with a 3 s stimulation duration did not result in significant temperature rises. Significant temperature rises up to 0.13°C (p < 0. 05) were observed with 6 and 9 s stimulation durations for blue and green light.

CONCLUSION

The proposed processing pipeline can be useful for the design of optogenetic stimulation protocols and for monitoring confounding heating effects.

Motion-robust, multi-slice, real-time MR thermometry for MR-guided thermal therapy in abdominal organs

PURPOSE

To develop an effective and practical reconstruction pipeline to achieve motion-robust, multi-slice, real-time MR thermometry for monitoring thermal therapy in abdominal organs.

METHODS

The application includes a fast spiral magnetic resonance imaging (MRI) pulse sequence and a real-time reconstruction pipeline based on multi-baseline proton resonance frequency shift (PRFS) method with visualization of temperature imaging. The pipeline supports multi-slice acquisition with minimal reconstruction lag. Simulations with a virtual motion phantom were performed to investigate the influence of the number of baselines and respiratory rate on the accuracy of temperature measurement. Phantom experiments with ultrasound heating were performed using a custom-made motion phantom to evaluate the performance of the pipeline. Lastly, experiments in healthy volunteers (N = 2) without heating were performed to evaluate the accuracy and stability of MR thermometry in abdominal organs (liver and kidney).

RESULTS

The multi-baseline approach with greater than 25 baselines resulted in minimal temperature errors in the simulation. Phantom experiments demonstrated a 713 ms update time for 3-slice acquisitions. Temperature maps with 30 baselines showed clear temperature distributions caused by ultrasound heating in the respiratory phantom. Finally, the pipeline was evaluated with physiologic motions in healthy volunteers without heating, which demonstrated the accuracy (root mean square error [RMSE]) of 1.23  ±   0.18 °C (liver) and 1.21  ±   0.17 °C (kidney) and precision of 1.13  ±   0.11 °C (liver) and 1.16  ±   0.15 °C (kidney) using 32 baselines.

CONCLUSIONS

The proposed real-time acquisition and reconstruction pipeline allows motion-robust, multi-slice, real-time temperature monitoring within the abdomen during free breathing.

A k-space-based method to measure and correct for temporal B0 field variations in MR temperature imaging

PURPOSE

Present a method to use change in phase in repeated Cartesian k-space measurements to monitor the change in magnetic field for dynamic MR temperature imaging.

METHODS

The method is applied to focused ultrasound heating experiments in a gelatin phantom and an ex vivo salt pork sample, without and with simulated respiratory motion.

RESULTS

In each experiment, phase variations due to B0 field drift and respiration were readily apparent in the measured phase difference. With correction, the SD of the temperature over time was reduced from 0.18°C to 0.14°C (no breathing) and from 0.81°C to 0.22°C (with breathing) for the gelatin phantom, and from 0.68°C to 0.13°C (no breathing) and from 1.06°C to 0.17°C (with breathing) for the pork sample. The accuracy in nonheated regions, assessed as the RMS error deviation from 0°C, improved from 1.70°C to 1.11°C (no breathing) and from 4.73°C to 1.47°C (with breathing) for the gelatin phantom, and from 5.95°C to 0.88°C (no breathing) and from 13.40°C to 1.73°C (with breathing) for the pork sample. The correction did not affect the temperature measurement accuracy in the heated regions.

CONCLUSION

This work demonstrates that phase changes resulting from variations in B0 due to drift and respiration, commonly seen in MR thermometry applications, can be measured directly from 3D Cartesian acquisition methods. The correction of temporal field variations using the presented technique improved temperature accuracy, reduced variability in nonheated regions, and did not reduce accuracy in heated regions.

Susceptibility artifact correction in MR thermometry for monitoring of mild radiofrequency hyperthermia using total field inversion

PURPOSE

MR temperature monitoring of mild radiofrequency hyperthermia (RF-HT) of cancer exploits the linear resonance frequency shift of water with temperature. Motion-induced susceptibility distribution changes cause artifacts that we correct here using the total field inversion (TFI) approach.

METHODS

The performance of TFI was compared to two background field removal (BFR) methods: Laplacian boundary value (LBV) and projection onto dipole fields (PDF). Data sets with spatial susceptibility change and B 0 -drift were simulated, phantom heating experiments were performed, four volunteer data sets at thermoneutral conditions as well as data from one cervical cancer, two sarcoma, and one seroma patients undergoing mild RF-HT were corrected using the proposed methods.

RESULTS

Simulations and phantom heating experiments revealed that using BFR or TFI preserves temperature-induced phase change, while removing susceptibility artifacts and B 0 -drift. TFI resulted in the least cumulative error for all four volunteers. Temperature probe information from four patient data sets were best depicted by TFI-corrected data in terms of accuracy and precision. TFI also performed best in case of the sarcoma treatment without temperature probe.

CONCLUSION

TFI outperforms previously suggested BFR methods in terms of accuracy and robustness. While PDF consistently overestimates susceptibility contribution, and LBV removes valuable pixel information, TFI is more robust and leads to more accurate temperature estimations.

Continuous cardiac thermometry via simultaneous catheter tracking and undersampled radial golden angle acquisition for radiofrequency ablation monitoring

The complexity of the MRI protocol is one of the factors limiting the clinical adoption of MR temperature mapping for real-time monitoring of cardiac ablation procedures and a push-button solution would ease its use. Continuous gradient echo golden angle radial acquisition combined with intra-scan motion correction and undersampled temperature determination could be a robust and more user-friendly alternative than the ultrafast GRE-EPI sequence which suffers from sensitivity to magnetic field susceptibility artifacts and requires ECG-gating. The goal of this proof-of-concept work is to establish the temperature uncertainty as well as the spatial and temporal resolutions achievable in an Agar-gel phantom and in vivo using this method. GRE radial golden angle acquisitions were used to monitor RF ablations in a phantom and in vivo in two sheep hearts with different slice orientations. In each case, 2D rigid motion correction based on catheter micro-coil signal, tracking its motion, was performed and its impact on the temperature imaging was assessed. The temperature uncertainty was determined for three spatial resolutions (1 × 1 × 3 mm³, 2 × 2 × 3 mm³, and 3 × 3 × 3 mm³) and three temporal resolutions (0.48, 0.72, and 0.97 s) with undersampling acceleration factors ranging from 2 to 17. The combination of radial golden angle GRE acquisition, simultaneous catheter tracking, intra-scan 2D motion correction, and undersampled thermometry enabled temperature monitoring in the myocardium in vivo during RF ablations with high temporal (< 1 s) and high spatial resolution. The temperature uncertainty ranged from 0.2 ± 0.1 to 1.8 ± 0.2 °C for the various temporal and spatial resolutions and, on average, remained superior to the uncertainty of an EPI acquisition while still allowing clinical monitoring of the RF ablation process. The proposed method is a robust and promising alternative to EPI acquisition to monitor in vivo RF cardiac ablations. Further studies remain required to improve the temperature uncertainty and establish its clinical applicability.

Deployable ultrasound applicators for endoluminal delivery of volumetric hyperthermia

PURPOSE

To investigate the design of an endoluminal deployable ultrasound applicator for delivering volumetric hyperthermia to deep tissue sites as a possible adjunct to radiation and chemotherapy.

METHOD

This study considers an ultrasound applicator consisting of two tubular transducers situated at the end of a catheter assembly, encased within a distensible conical shaped balloon-based reflector that redirects acoustic energy distally into the tissue. The applicator assembly can be inserted endoluminally or laparoscopically in a compact form and expanded after delivery to the target site. Comprehensive acoustic and biothermal simulations and parametric studies were employed in generalized 3D and patient-specific pancreatic head and body tumor models to characterize the acoustic performance and evaluate heating capabilities of the applicator by investigating the device at a range of operating frequencies, tissue acoustic and thermal properties, transducer configurations, power modulation, applicator positioning, and by analyzing the resultant 40, 41, and 43 °C isothermal volumes and penetration depth of the heating volume. Intensity distributions and volumetric temperature contours were calculated to define moderate hyperthermia boundaries.

RESULTS

Parametric studies demonstrated the frequency selection to control volume and depth of therapeutic heating from 62 to 22 cm³ and 4 to 2.6 cm as frequency ranges from 1 MHz to 4.7 MHz, respectively. Width of the heating profile tracks closely with the aperture. Water cooling within the reflector balloon was effective in controlling temperature to 37 °C maximum within the luminal wall. Patient-specific studies indicated that applicators with extended OD in the range of 3.6-6.2 cm with 0.5-1 cm long and 1 cm OD transducers can heat volumes of 1.1-7 cm³, 3-26 cm³, and 3.3-37.4 cm³ of pancreatic body and head tumors above 43, 41, and 40 °C, respectively.

CONCLUSION

In silico studies demonstrated the feasibility of combining endoluminal ultrasound with an integrated expandable balloon reflector for delivering volumetric hyperthermia in regions adjacent to body lumens and cavities.

AAPM Task Group 241: A medical physicist's guide to MRI-guided focused ultrasound body systems

Magnetic resonance-guided focused ultrasound (MRgFUS) is a completely non-invasive technology that has been approved by FDA to treat several diseases. This report, prepared by the American Association of Physicist in Medicine (AAPM) Task Group 241, provides background on MRgFUS technology with a focus on clinical body MRgFUS systems. The report addresses the issues of interest to the medical physics community, specific to the body MRgFUS system configuration, and provides recommendations on how to successfully implement and maintain a clinical MRgFUS program. The following sections describe the key features of typical MRgFUS systems and clinical workflow and provide key points and best practices for the medical physicist. Commonly used terms, metrics and physics are defined and sources of uncertainty that affect MRgFUS procedures are described. Finally, safety and quality assurance procedures are explained, the recommended role of the medical physicist in MRgFUS procedures is described, and regulatory requirements for planning clinical trials are detailed. Although this report is limited in scope to clinical body MRgFUS systems that are approved or currently undergoing clinical trials in the United States, much of the material presented is also applicable to systems designed for other applications.

Techniques for Temperature Monitoring of Myocardial Tissue Undergoing Radiofrequency Ablation Treatments: An Overview

Cardiac radiofrequency ablation (RFA) has received substantial attention for the treatment of multiple arrhythmias. In this scenario, there is an ever-growing demand for monitoring the temperature trend inside the tissue as it may allow an accurate control of the treatment effects, with a consequent improvement of the clinical outcomes. There are many methods for monitoring temperature in tissues undergoing RFA, which can be divided into invasive and non-invasive. This paper aims to provide an overview of the currently available techniques for temperature detection in this clinical scenario. Firstly, we describe the heat generation during RFA, then we report the principle of work of the most popular thermometric techniques and their features. Finally, we introduce their main applications in the field of cardiac RFA to explore the applicability in clinical settings of each method.

Deep correction of breathing-related artifacts in real-time MR-thermometry

Real-time MR-imaging has been clinically adapted for monitoring thermal therapies since it can provide on-the-fly temperature maps simultaneously with anatomical information. However, proton resonance frequency based thermometry of moving targets remains challenging since temperature artifacts are induced by the respiratory as well as physiological motion. If left uncorrected, these artifacts lead to severe errors in temperature estimates and impair therapy guidance. In this study, we evaluated deep learning for on-line correction of motion related errors in abdominal MR-thermometry. For this, a convolutional neural network (CNN) was designed to learn the apparent temperature perturbation from images acquired during a preparative learning stage prior to hyperthermia. The input of the designed CNN is the most recent magnitude image and no surrogate of motion is needed. During the subsequent hyperthermia procedure, the recent magnitude image is used as an input for the CNN-model in order to generate an on-line correction for the current temperature map. The method's artifact suppression performance was evaluated on 12 free breathing volunteers and was found robust and artifact-free in all examined cases. Furthermore, thermometric precision and accuracy was assessed for in vivo ablation using high intensity focused ultrasound. All calculations involved at the different stages of the proposed workflow were designed to be compatible with the clinical time constraints of a therapeutic procedure.

Systematic review of pre-clinical and clinical devices for magnetic resonance-guided radiofrequency hyperthermia

Clinical trials have demonstrated the therapeutic benefits of adding radiofrequency (RF) hyperthermia (HT) as an adjuvant to radio- and chemotherapy. However, maximum utilization of these benefits is hampered by the current inability to maintain the temperature within the desired range. RF HT treatment quality is usually monitored by invasive temperature sensors, which provide limited data sampling and are prone to infection risks. Magnetic resonance (MR) temperature imaging has been developed to overcome these hurdles by allowing noninvasive 3D temperature monitoring in the target and normal tissues. To exploit this feature, several approaches for inserting the RF heating devices into the MR scanner have been proposed over the years. In this review, we summarize the status quo in MR-guided RF HT devices and analyze trends in these hybrid hardware configurations. In addition, we discuss the various approaches, extract best practices and identify gaps regarding the experimental validation procedures for MR - RF HT, aimed at converging to a common standard in this process.

Novel fractionated ultrashort thermal exposures with MRI-guided focused ultrasound for treating tumors with thermosensitive drugs

Thermosensitive liposomes represent an important paradigm in oncology, where hyperthermia-mediated release coupled with thermal bioeffects enhance the effectiveness of chemotherapy. Their widespread clinical adoption hinges upon performing controlled targeted hyperthermia, and a leading candidate to achieve this is temperature-based magnetic resonance imaging (MRI)-guided focused ultrasound (MRgFUS). However, the current approach to hyperthermia involves exposures lasting tens of minutes to hours, which is not possible to achieve in many circumstances because of blood vessel cooling and respiratory motion. Here, we investigate a novel approach to overcome these limitations: to use fractionated ultrashort (~30 s) thermal exposures (~41° to 45°C) to release doxorubicin from a thermosensitive liposome. This is first demonstrated in a dorsal chamber tumor model using two-photon microscopy. Thermal exposures were then conducted with a rabbit tumor model using a custom MRgFUS system incorporating temperature feedback control. Drug release was confirmed, and longitudinal experiments demonstrated profoundly enhanced tumor growth inhibition and survival.

Field drift correction of proton resonance frequency shift temperature mapping with multichannel fast alternating nonselective free induction decay readouts

PURPOSE

To demonstrate that proton resonance frequency shift MR thermometry (PRFS-MRT) acquisition with nonselective free induction decay (FID), combined with coil sensitivity profiles, allows spatially resolved B0 drift-corrected thermometry.

METHODS

Phantom experiments were performed at 1.5T and 3T. Acquisition of PRFS-MRT and FID were performed during MR-guided high-intensity focused ultrasound heating. The phase of the FIDs was used to estimate the change in angular frequency δωdrift per coil element. Two correction methods were investigated: (1) using the average δωdrift over all coil elements (0th-order) and (2) using coil sensitivity profiles for spatially resolved correction. Optical probes were used for independent temperature verification. In-vivo feasibility of the methods was evaluated in the leg of 1 healthy volunteer at 1.5T.

RESULTS

In 30 minutes, B0 drift led to an apparent temperature change of up to -18°C and -98°C at 1.5T and 3T, respectively. In the sonicated area, both corrections had a median error of 0.19°C at 1.5T and -0.54°C at 3T. At 1.5T, the measured median error with respect to the optical probe was -1.28°C with the 0th-order correction and improved to 0.43°C with the spatially resolved correction. In vivo, without correction the spatiotemporal median of the apparent temperature was at -4.3°C and interquartile range (IQR) of 9.31°C. The 0th-order correction had a median of 0.75°C and IQR of 0.96°C. The spatially resolved method had the lowest median at 0.33°C and IQR of 0.80°C.

CONCLUSION

FID phase information from individual receive coil elements allows spatially resolved B0 drift correction in PRFS-based MRT.

Simulation-based design and characterization of a microwave applicator for MR-guided hyperthermia experimental studies in small animals

Purpose

The objective of this study was to design and characterize a 2.45 GHz microwave hyperthermia applicator for delivering hyperthermia in experimental small animals to 2 - 4 mm diameter targets located 1 - 3 mm from the skin surface, with minimal heating of the surrounding tissue, under 14.1 T MRI real-time monitoring and feedback control.

Materials and methods

An experimentally validated 3D computational model was employed to design and characterize a non-invasive directional water-cooled microwave hyperthermia applicator. We assessed the effects of: reflector geometry, monopole shape, cooling water temperature, and flow rate on spatial-temperature profiles. The system was integrated with real-time MR thermometry and feedback control to monitor and maintain temperature elevations in the range of 4 - 5 °C at 1 - 3 mm from the applicator surface. The quality of heating was quantified by determining the fraction of the target volume heated to the desired temperature, and the extent of heating in non-targeted regions.

Results

Model-predicted hyperthermic profiles were in good agreement with experimental measurements (Dice Similarity Coefficient of 0.95 - 0.99). Among the four considered criteria, a reflector aperture angle of 120 °, S-shaped monopole antenna with 0.6 mm displacement, and coolant flow rate of 150 ml/min were selected as the end result of the applicator design. The temperature of circulating water and input power were identified as free variables, allowing considerable flexibility in heating target sizes within varying distances from the applicator surface. 2 - 4 mm diameter targets positioned 1 - 3 mm from the applicator surface were heated to hyperthermic temperatures, with target coverage ratio ranging between 76 - 93 % and 11 - 26 % of non-targeted tissue heated.

Conclusion

We have designed an experimental platform for MR-guided hyperthermia, incorporating a microwave applicator integrated with temperature-based feedback control to heat deep-seated targets for experimental studies in small animals.

Recent technological advancements in thermometry

Thermometry is the key factor for achieving successful thermal therapy. Although invasive thermometry with a probe has been used for more than four decades, this method can only detect the local temperature within the probing volume. Noninvasive temperature imaging using a tomographic technique is ideal for monitoring hot-spot formation in the human body. Among various techniques, such as X-ray computed tomography, microwave tomography, echo sonography, and magnetic resonance (MR) imaging, the proton resonance frequency shift method of MR thermometry is the only method currently available for clinical practice because its temperature sensitivity is consistent in most aqueous tissues and can be easily observed using common clinical scanners. New techniques are being proposed to improve the robustness of this method against tissue motion. MR techniques for fat thermometry were also developed based on relaxation times. One of the latest non-MR techniques to attract attention is photoacoustic imaging.

Considerations for ultrasound exposure during transcranial MR acoustic radiation force imaging

The aim of this study was to improve the sensitivity of magnetic resonance-acoustic radiation force imaging (MR-ARFI) to minimize pressures required to localize focused ultrasound (FUS) beams, and to establish safe FUS localization parameters for ongoing ultrasound neuromodulation experiments in living non-human primates. We developed an optical tracking method to ensure that the MR-ARFI motion-encoding gradients (MEGs) were aligned with a single-element FUS transducer and that the imaged slice was prescribed at the optically tracked location of the acoustic focus. This method was validated in phantoms, which showed that MR-ARFI-derived displacement sensitivity is maximized when the MR-ARFI MEGs were maximally aligned with the FUS propagation direction. The method was then applied in vivo to acquire displacement images in two healthy macaque monkeys (M fascicularis) which showed the FUS beam within the brain. Temperature images were acquired using MR thermometry to provide an estimate of in vivo brain temperature changes during MR-ARFI, and pressure and thermal simulations of the acoustic pulses were performed using the k-Wave package which showed no significant heating at the focus of the FUS beam. The methods presented here will benefit the multitude of transcranial FUS applications as well as future human applications.

Assessment of left ventricle magnetic resonance temperature stability in patients in the presence of arrhythmias

BACKGROUND

Magnetic resonance (MR) thermometry allows visualization of lesion formation in real-time during cardiac radiofrequency (RF) ablation. The present study was performed to evaluate the precision of MR thermometry without RF heating in patients exhibiting cardiac arrhythmia in a clinical setting. The evaluation relied on quantification of changes in temperature measurements caused by noise and physiological motion.

METHODS

Fourteen patients referred for cardiovascular magnetic resonance imaging underwent an extra sequence to test the temperature mapping stability during free-breathing acquisition. Phase images were acquired using a multi-slice, cardiac-triggered, single-shot echo planar imaging sequence. Temperature maps were calculated and displayed in real-time while the electrocardiogram (ECG) was recorded. The precision of temperature measurement was assessed by measuring the temporal standard deviation and temporal mean of consecutive temperature maps over a period of three minutes. The cardiac cycle was analyzed from ECG recordings to quantify the impact of arrhythmia events on the precision of temperature measurement. Finally, two retrospective strategies were tested to remove acquisition dynamics related either to arrhythmia events or sudden breathing motion.

RESULTS

ECG synchronization allowed categorization of inter-beat intervals (RR) into distinct beat morphologies. Five patients were in stable sinus rhythm, while nine patients showed irregular RR intervals due to ectopic beats. An average temporal standard deviation of temperature of 1.6°C was observed in patients under sinus rhythm with a frame rate corresponding to the heart rate of the patient. The temporal standard deviation rose to 2.5°C in patients with arrhythmia. The retrospective rejection strategies increased the temperature precision measurement while maintaining a sufficient frame rate.

CONCLUSIONS

Our results indicated that real-time cardiac MR thermometry shows good precision in patients under clinical conditions, even in the presence of arrhythmia. By providing real-time visualization of temperature distribution within the myocardium during RF delivery, MR thermometry could prevent insufficient or excessive heating and thus improve safety and efficacy.

Reducing temperature errors in transcranial MR-guided focused ultrasound using a reduced-field-of-view sequence

PURPOSE

To reduce temperature errors due to water motion in transcranial MR-guided focused ultrasound (tcMRgFUS) ablation.

THEORY AND METHODS

In tcMRgFUS, water is circulated in the transducer bowl around the patient's head for acoustic coupling and heat removal. The water moves during sonications that are monitored by MR thermometry, which causes it to alias into the brain and create temperature errors. To reduce these errors, a two-dimensional excitation pulse was implemented in a gradient-recalled echo thermometry sequence. The pulse suppresses water signal by selectively exciting the brain only, which reduces the imaging FOV. Improvements in temperature precision compared to the conventional full-FOV scan were evaluated in healthy subject scans outside the tcMRgFUS system, gel phantom scans in the system with heating, and in 2×-accelerated head phantom scans in the system without heating.

RESULTS

In vivo temperature precision (standard deviation of temperature errors) outside the tcMRgFUS system was improved 43% on average, due to the longer TR and TE of the reduced-FOV sequence. In the phantom heating experiments, the hot spot was less distorted in the reduced-FOV scans, and background temperature precision was improved 59% on average. In the accelerated head phantom temperature reconstructions, temperature precision was improved 89% using the reduced-FOV sequence.

CONCLUSIONS

Reduced-FOV temperature imaging alleviates temperature errors due to water bath motion in tcMRgFUS, and enables accelerated temperature mapping with greater precision.

Simultaneous multislice MRI thermometry with a single coil using incoherent blipped-controlled aliasing

PURPOSE

To increase volume coverage in real-time MR thermometry for transcranial MR-guided focused ultrasound (tcMRgFUS) ablation, without multiple receive coils.

THEORY AND METHODS

Multiband excitation and incoherent blipped-controlled aliasing were implemented in a 2DFT pulse sequence used clinically for tcMRgFUS, and an extended k-space hybrid reconstruction was developed that recovers slice-separated temperature maps assuming that heating is focal, given slice-separated pretreatment images. Simulations were performed to characterize slice leakage, the number of slices that can be simultaneously imaged with low-temperature error, and robustness across random slice-phase k-space permutations. In vivo experiments were performed using a single receive coil without heating to measure temperature precision, and gel phantom FUS experiments were performed to test the method with heating and with a water bath.

RESULTS

Simulations showed that with large hot spots and identical magnitude images on each slice, up to three slices can be simultaneously imaged with less than 1 ∘ C temperature root-mean-square error. They also showed that hot spots do not alias coherently between slices, and that an average 86% of random slice-phase k-space permutations yielded less than 1 ∘ C temperature error. Temperature precision was not degraded compared to single-slice imaging in the in vivo SMS scans, and the gel phantom SMS temperature maps closely tracked single-slice temperature in the hot spot, with no coherent aliasing to other slices.

CONCLUSIONS

Incoherent controlled aliasing SMS enables accurate reconstruction of focal heating maps from two or three slices simultaneously, using a single receive coil and a sparsity-promoting temperature reconstruction.

A planning strategy for combined motion-assisted/gated MR guided focused ultrasound treatment of the pancreas

Objective: To develop and evaluate a combined motion-assisted/gated MRHIFU heating strategy designed to accelerate the treatment procedure by reducing the required number of sonications to ablate a target volume in the pancreas. Methods: A planning method for combined motion-assisted/gated MRHIFU using 4D-MRI and motion characterization is introduced. Six healthy volunteers underwent 4D-MRI for target motion characterization on a 3.0-T clinical scanner. Using displacement patterns, simulations were performed for all volunteers for three sonication approaches: gated, combined motion-assisted/gated, and static. The number of sonications needed to ablate the pancreas head was compared. The influence of displacement amplitude and target volume size was investigated. Spherical target volumes (8, 15, 20 and 34 mL) and displacement amplitudes ranging from 5 to 25 mm were evaluated. For this case, the number of sonications required to ablate the whole target was determined. Results: The number of required sonications was lowest for a static target, 62 on average (range 49-78). The gated approach required most sonications, 126 (range 97-159). The combined approach was almost as efficient as the hypothetical static case, with an average of 78 (range 53-123). Simulations showed that with a 5-mm displacement amplitude, the target could be treated by making use of motion-assisted MRHIFU sonications only. In that case, this approach allowed the lowest number of sonication, while for 10 mm and above, the number of required sonications increased. Conclusion: The use of a combined motion-assisted/gated MRHIFU strategy may accelerate tumor ablation in the pancreas when respiratory-induced displacement amplitudes are between 5 and 10 mm.

A variable flip angle golden-angle-ordered 3D stack-of-radial MRI technique for simultaneous proton resonant frequency shift and T1 -based thermometry

PURPOSE

To develop and evaluate a variable-flip-angle golden-angle-ordered 3D stack-of-radial MRI technique for simultaneous proton resonance frequency shift (PRF) and T₁ -based thermometry in aqueous and adipose tissues, respectively.

METHODS

The proposed technique acquires multiecho radial k-space data in segments with alternating flip angles to measure 3D temperature maps dynamically on the basis of PRF and T₁ . A sliding-window k-space weighted image contrast filter is used to increase temporal resolution. PRF is measured in aqueous tissues and T₁ in adipose tissues using fat/water masks. The accuracy for T₁ quantification was evaluated in a reference T₁ /T₂ phantom. In vivo nonheating experiments were conducted in healthy subjects to evaluate the stability of PRF and T₁ in the brain, prostate, and breast. The proposed technique was used to monitor high-intensity focused ultrasound (HIFU) ablation in ex vivo porcine fat/muscle tissues and compared to temperature probe readings.

RESULTS

The proposed technique achieved 3D coverage with 1.1-mm to 1.3-mm in-plane resolution and 2-s to 5-s temporal resolution. During 20 to 30 min of nonheating in vivo scans, the temporal coefficient of variation for T₁ was <5% in the brain, prostate, and breast fatty tissues, while the standard deviation of relative PRF temperature change was within 3°C in aqueous tissues. During ex vivo HIFU ablation, the temperatures measured by PRF and T₁ were consistent with temperature probe readings, with an absolute mean difference within 2°C.

CONCLUSION

The proposed technique achieves simultaneous PRF and T₁ -based dynamic 3D MR temperature mapping in aqueous and adipose tissues. It may be used to improve MRI-guided thermal procedures.

A rapid beam simulation framework for transcranial focused ultrasound

Transcranial focused ultrasound is a non-invasive therapeutic modality that can be used to treat essential tremor. Beams of energy are focused into a small spot in the thalamus, resulting in tissue heating and ablation. Here, we report on a rapid 3D numeric simulation framework that can be used to predict focal spot characteristics prior to the application of ultrasound. By comparing with magnetic resonance proton resonance frequency shift thermometry (MR thermometry) data acquired during treatments of essential tremor, we verified that our simulation framework can be used to predict focal spot position, and with patient-specific calibration, predict focal spot temperature rise. Preliminary data suggests that lateral smearing of the focal spot can be simulated. The framework may also be relevant for other therapeutic ultrasound applications such as blood brain barrier opening and neuromodulation.

Magnetic resonance thermometry and its biological applications - Physical principles and practical considerations

Most parameters that influence the magnetic resonance imaging (MRI) signal experience a temperature dependence. The fact that MRI can be used for non-invasive measurements of temperature and temperature change deep inside the human body has been known for over 30 years. Today, MR temperature imaging is widely used to monitor and evaluate thermal therapies such as radio frequency, microwave, laser, and focused ultrasound therapy. In this paper we cover the physical principles underlying the biological applications of MR temperature imaging and discuss practical considerations and remaining challenges. For biological tissue, the MR signal of interest comes mostly from hydrogen protons of water molecules but also from protons in, e.g., adipose tissue and various metabolites. Most of the discussed methods, such as those using the proton resonance frequency (PRF) shift, T₁, T₂, and diffusion only measure temperature change, but measurements of absolute temperatures are also possible using spectroscopic imaging methods (taking advantage of various metabolite signals as internal references) or various types of contrast agents. Currently, the PRF method is the most used clinically due to good sensitivity, excellent linearity with temperature, and because it is largely independent of tissue type. Because the PRF method does not work in adipose tissues, T₁- and T₂-based methods have recently gained interest for monitoring temperature change in areas with high fat content such as the breast and abdomen. Absolute temperature measurement methods using spectroscopic imaging and contrast agents often offer too low spatial and temporal resolution for accurate monitoring of ablative thermal procedures, but have shown great promise in monitoring the slower and usually less spatially localized temperature change observed during hyperthermia procedures. Much of the current research effort for ablative procedures is aimed at providing faster measurements, larger field-of-view coverage, simultaneous monitoring in aqueous and adipose tissues, and more motion-insensitive acquisitions for better precision measurements in organs such as the heart, liver, and kidneys. For hyperthermia applications, larger coverage, motion insensitivity, and simultaneous aqueous and adipose monitoring are also important, but great effort is also aimed at solving the problem of long-term field drift which gets interpreted as temperature change when using the PRF method.

Breath-hold MR-HIFU hyperthermia: phantom and in vivo feasibility

Background: The use of magnetic resonance imaging-guided high-intensity focused ultrasound (MR-HIFU) to deliver mild hyperthermia requires stable temperature mapping for long durations. This study evaluates the effects of respiratory motion on MR thermometry precision in pediatric subjects and determines the in vivo feasibility of circumventing breathing-related motion artifacts by delivering MR thermometry-controlled HIFU mild hyperthermia during repeated forced breath holds.Materials and methods: Clinical and preclinical studies were conducted. Clinical studies were conducted without breath-holds. In phantoms, breathing motion was simulated by moving an aluminum block towards the phantom along a sinusoidal trajectory using an MR-compatible motion platform. In vivo experiments were performed in ventilated pigs. MR thermometry accuracy and stability were evaluated.Results: Clinical data confirmed acceptable MR thermometry accuracy (0.12-0.44 °C) in extremity tumors, but not in the tumors in the chest/spine and pelvis. In phantom studies, MR thermometry accuracy and stability improved to 0.37 ± 0.08 and 0.55 ± 0.18 °C during simulated breath-holds. In vivo MR thermometry accuracy and stability in porcine back muscle improved to 0.64 ± 0.22 and 0.71 ± 0.25 °C during breath-holds. MR-HIFU hyperthermia delivered during intermittent forced breath holds over 10 min duration heated an 18-mm diameter target region above 41 °C for 10.0 ± 1.0 min, without significant overheating. For a 10-min mild hyperthermia treatment, an optimal treatment effect (TIR > 9 min) could be achieved when combining 36-60 s periods of forced apnea with 60-155.5 s free-breathing.Conclusion: MR-HIFU delivery during forced breath holds enables stable control of mild hyperthermia in targets adjacent to moving anatomical structures.

Transurethral high-intensity ultrasound for treatment of stress urinary incontinence (SUI): simulation studies with patient-specific models

BACKGROUND

Stress urinary incontinence (SUI) is prevalent in adult women, attributed to weakened endopelvic supporting tissues, and typically treated using drugs and invasive surgical procedures. The objective of this in silico study is to explore transurethral high-intensity ultrasound for delivery of precise thermal therapy to the endopelvic tissues adjacent to the mid-urethra, to induce thermal remodeling as a potential minimally invasive treatment alternative.

METHODS

3D acoustic (Rayleigh-Sommerfeld) and biothermal (Pennes bioheat) models of the ultrasound applicator and surrounding tissues were devised. Parametric studies over transducer configuration [frequency, radius-of-curvature (ROC)] and treatment settings (power, duration) were performed, and select cases on patient-specific models were used for further evaluation. Transient temperature and thermal dose distributions were calculated, and temperature and dose metrics reported.

RESULTS

Configurations using a 5-MHz curvilinear transducer (3.5 × 10 mm, 28 mm ROC) with single 90 s sonication can create heated zones with 11 mm penetration (>50 °C) while sparing the inner 1.8 mm (<45 °C) radial depth of the urethral mucosa. Sequential and discrete applicator rotations can sweep out bilateral coagulation volumes (1.4 W power, 15° rotations, 600 s total time), produce large volumetric (1124 mm³ above 60 EM43 °C) and wide angular (∼50.5° per lateral sweep) coverage, with up to 15.6 mm thermal penetration and at least 1.6 mm radial urethral protection (<5 EM43 °C).

CONCLUSION

Transurethral applicators with curvilinear ultrasound transducers can deliver spatially selective temperature elevations to lateral mid-urethral targets as a possible means to tighten the endopelvic fascia and adjacent tissues.

Multi-echo MR thermometry using iterative separation of baseline water and fat images

PURPOSE

To perform multi-echo water/fat separated proton resonance frequency (PRF)-shift temperature mapping.

METHODS

State-of-the-art, iterative multi-echo water/fat separation algorithms produce high-quality water and fat images in the absence of heating but are not suitable for real-time imaging due to their long compute times and potential errors in heated regions. Existing fat-referenced PRF-shift temperature reconstruction methods partially address these limitations but do not address motion or large time-varying and spatially inhomogeneous B₀ shifts. We describe a model-based temperature reconstruction method that overcomes these limitations by fitting a library of separated water and fat images measured before heating directly to multi-echo data measured during heating, while accounting for the PRF shift with temperature.

RESULTS

Simulations in a mixed water/fat phantom with focal heating showed that the proposed algorithm reconstructed more accurate temperature maps in mixed tissues compared to a fat-referenced thermometry method. In a porcine phantom experiment with focused ultrasound heating at 1.5 Tesla, temperature maps were accurate to within 1^∘ C of fiber optic probe temperature measurements and were calculated in 0.47 s per time point. Free-breathing breast and liver imaging experiments demonstrated motion and off-resonance compensation. The algorithm can also accurately reconstruct water/fat separated temperature maps from a single echo during heating.

CONCLUSIONS

The proposed model-based water/fat separated algorithm produces accurate PRF-shift temperature maps in mixed water and fat tissues in the presence of spatiotemporally varying off-resonance and motion.

Motion compensation using principal component analysis and projection onto dipole fields for abdominal magnetic resonance thermometry

PURPOSE

High intensity focused ultrasound (HIFU) has the potential to locally and non-invasively treat cancer with fewer side effects than alternative therapies. However, motion and tissue heterogeneity in the abdomen can compromise the HIFU focus and confound current thermometry methods.

METHODS

The proposed thermometry method combines principal component analysis (PCA), as a multi-baseline technique, and projection onto dipole fields (PDF), as a near-referenceless method. PCA forgoes tracking tools by projecting incoming images onto a subspace spanning the motion history. PDF is subsequently used to synthesize the naturally feasible components of the residual phase using a magnetic dipole model. This leaves only the phase shifts that are induced by HIFU.

RESULTS

With in vivo measurements, in porcine and human kidneys, the mean pixel-wise temperature SD was 0.86 ± 0.41°C in selected regions of interest (ROIs) across all data sets, without any user-interaction or supplementary tracking tools. This is an improvement over a benchmark hybrid method, which scored 1.36 ± 1.20°C on the same data. Uncorrected subtraction of the data yielded a score of 3.02 ± 2.87°C.

CONCLUSION

The PCA-PDF hybrid method achieves superior artifact correction by exploiting the motion history and intrinsic magnetic susceptibility of the underlying tissue.

Motion Correction in Proton Resonance Frequency-based Thermometry in the Liver

The unique ability of magnetic resonance imaging to measure temperature noninvasively, in vivo, makes it an attractive tool for monitoring interventional procedures, such as radiofrequency or microwave ablation in real-time. The most frequently used approach for magnetic resonance-based temperature measurement is proton resonance frequency (PRF) thermometry. Although it has many advantages, including tissue-independence and real-time capability, the main drawback is its motion sensitivity. This is likely the reason PRF thermometry in moving organs, such as the liver, is not commonly used in the clinical arena. In recent years, however, several developments suggest that motion-corrected thermometry in the liver is achievable. The present article summarizes the diverse attempts to correct thermometry in the liver. Therefore, the physical principle of PRF is introduced, with additional references for necrosis zone estimation and how to deal with fat phase modulation, and main magnetic field drifts. The primary categories of motion correction are presented, including general methods for motion compensation and library-based approaches, and referenceless thermometry and hybrid methods. Practical validation of the described methods in larger patient groups will be necessary to establish accurate motion-corrected thermometry in the clinical arena, with the goal of complete liver tumor ablation.

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.

Referenced MR thermometry using three-echo phase-based fat water separation method

A three-point image reconstruction method for internally referenced MR thermometry was developed. The technique exploits the fact that temperature-induced changes in the water resonance frequency are small relative to the chemical shift difference between water and fat signals. This property enabled the use of small angle approximations to derive an analytic phase-based fat-water separation method for MR thermometry. Ethylene glycol and cream cool-down experiments were performed to validate measurement technique. Over a cool-down temperature range of 20 °C, maximum deviation between probe and MR measurement (averaged over 1.3 cm³ region surrounding probe) was 0.6 °C and 1.1 °C for ethylene glycol and cream samples, respectively.

MR techniques for guiding high-intensity focused ultrasound (HIFU) treatments

To make full use of the ability of magnetic resonance (MR) to guide high-intensity focused ultrasound (HIFU) treatment, effort has been made to improve techniques for thermometry, motion tracking, and sound beam visualization. For monitoring rapid temperature elevation with proton resonance frequency (PRF) shift, data acquisition and processing can be accelerated with parallel imaging and/or sparse sampling in conjunction with appropriate signal processing methods. Thermometry should be robust against tissue motion, motion-induced magnetic field variation, and susceptibility change. Thus, multibaseline, referenceless, or hybrid techniques have become important. In cases with adipose or bony tissues, for which PRF shift cannot be used, thermometry with relaxation times or signal intensity may be utilized. Motion tracking is crucial not only for thermometry but also for targeting the focus of an ultrasound in moving organs such as the liver, kidney, or heart. Various techniques for motion tracking, such as those based on an anatomical image atlas with optical-flow displacement detection, a navigator echo to seize the diaphragm position, and/or rapid imaging to track vessel positions, have been proposed. Techniques for avoiding the ribcage and near-field heating have also been examined. MR acoustic radiation force imaging (MR-ARFI) is an alternative to thermometry that can identify the location and shape of the focal spot and sound beam path. This technique could be useful for treating heterogeneous tissue regions or performing transcranial therapy. All of these developments, which will be discussed further in this review, expand the applicability of HIFU treatments to a variety of clinical targets while maintaining safety and precision.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 4 J. Magn. Reson. Imaging 2018;47:316-331.

MRI use for atrial tissue characterization in arrhythmias and for EP procedure guidance

We review the utilization of magnetic resonance imaging methods for classifying atrial tissue properties that act as a substrate for common cardiac arrhythmias, such as atrial fibrillation. We then review state-of-the-art methods for mapping this substrate as a predicate for treatment, as well as methods used to ablate the electrical pathways that cause arrhythmia and restore patients to sinus rhythm.

Fluid filling of the digestive tract for improved proton resonance frequency shift-based MR thermometry in the pancreas

PURPOSE

To demonstrate that fluid filling of the digestive tract improves the performance of respiratory motion-compensated proton resonance frequency shift (PRFS)-based magnetic resonance (MR) thermometry in the pancreas.

MATERIALS AND METHODS

In seven volunteers (without heating), we evaluated PRFS thermometry in the pancreas with and without filling of the surrounding digestive tract. All data acquisition was performed at 1.5T, then all datasets were analyzed and compared with three different PRFS respiratory motion-compensated thermometry methods: gating, multibaseline, and referenceless. The temperature precision of the different methods was evaluated by assessing temperature standard deviation over time, while a simulation experiment was used to study the accuracy of the methods.

RESULTS

Without fluid intake, errors in temperature precision in the pancreas up to 10°C were observed for all evaluated methods. After liquid intake, temperature precision improved to median values between 1.8 and 2.9°C. The simulations showed that gating had the lowest accuracy, with errors up to 7°C. Multibaseline and referenceless thermometry performed better, with a median error in the pancreas between -3 and +3°C after fluid intake, for all volunteers.

CONCLUSION

Preparation of the digestive tract near the pancreas by filling it with fluid improved MR thermometry precision and accuracy for all common respiratory motion-compensated methods evaluated. These improvements are attributed to reducing field inhomogeneity in the pancreas.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2018;47:692-701.

Spatially-segmented undersampled MRI temperature reconstruction for transcranial MR-guided focused ultrasound

Background

Volumetric thermometry with fine spatiotemporal resolution is desirable to monitor MR-guided focused ultrasound (MRgFUS) procedures in the brain, but requires some form of accelerated imaging. Accelerated MR temperature imaging methods have been developed that undersample k-space and leverage signal correlations over time to suppress the resulting undersampling artifacts. However, in transcranial MRgFUS treatments, the water bath surrounding the skull creates signal variations that do not follow those correlations, leading to temperature errors in the brain due to signal aliasing.

Methods

To eliminate temperature errors due to the water bath, a spatially-segmented iterative reconstruction method was developed. The method fits a k-space hybrid signal model to reconstruct temperature changes in the brain, and a conventional MR signal model in the water bath. It was evaluated using single-channel 2DFT Cartesian, golden angle radial, and spiral data from gel phantom heating, and in vivo 8-channel 2DFT data from a FUS thalamotomy. Water bath signal intensity in phantom heating images was scaled between 0-100% to investigate its effect on temperature error. Temperature reconstructions of retrospectively undersampled data were performed using the spatially-segmented method, and compared to conventional whole-image k-space hybrid (phantom) and SENSE (in vivo) reconstructions.

Results

At 100% water bath signal intensity, 3 ×-undersampled spatially-segmented temperature reconstruction error was nearly 5-fold lower than the whole-image k-space hybrid method. Temperature root-mean square error in the hot spot was reduced on average by 27 × (2DFT), 5 × (radial), and 12 × (spiral) using the proposed method. It reduced in vivo error 2 × in the brain for all acceleration factors, and between 2 × and 3 × in the temperature hot spot for 2-4 × undersampling compared to SENSE.

Conclusions

Separate reconstruction of brain and water bath signals enables accelerated MR temperature imaging during MRgFUS procedures with low errors due to undersampling using Cartesian and non-Cartesian trajectories. The spatially-segmented method benefits from multiple coils, and reconstructs temperature with lower error compared to measurements from SENSE-reconstructed images. The acceleration can be applied to increase volumetric coverage and spatiotemporal resolution.

Dual-echo Z-shimmed proton resonance frequency-shift magnetic resonance thermometry near metallic ablation probes: Technique and temperature precision

PURPOSE

To improve the precision of proton resonance frequency-shift magnetic resonance thermometry near ablation probes by recovering near-probe image signals that are typically lost due to magnetic susceptibility-induced field distortions.

METHODS

A dual-echo gradient-recalled echo sequence was implemented, in which the first echo was under- or over-refocused in the slice dimension to recover image signal and temperature precision near a probe, and the second echo was fully refocused to obtain image signal everywhere else in the slice. A penalized maximum likelihood algorithm was implemented to estimate a single temperature map from both echoes. Agar phantom and ex vivo experiments with and without microwave heating at 3 T evaluated how much temperature precision was improved near a microwave ablator compared to a conventional single-echo scan as a function of slice and needle orientation in the magnet.

RESULTS

The number of near-probe voxels with temperature standard deviation σ>1°C was decreased by 51% in the phantom experiment, averaged across orientations, and by 31% in the pork. Temperature maps near the probe were more smoother and more complete in all orientations.

CONCLUSION

Dual-echo z-shimmed temperature imaging can recover image signal for more precise temperature mapping near metallic ablation probes. Magn Reson Med 78:2299-2306, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Feasibility of real-time MR thermal dose mapping for predicting radiofrequency ablation outcome in the myocardium in vivo

BACKGROUND

Clinical treatment of cardiac arrhythmia by radiofrequency ablation (RFA) currently lacks quantitative and precise visualization of lesion formation in the myocardium during the procedure. This study aims at evaluating thermal dose (TD) imaging obtained from real-time magnetic resonance (MR) thermometry on the heart as a relevant indicator of the thermal lesion extent.

METHODS

MR temperature mapping based on the Proton Resonance Frequency Shift (PRFS) method was performed at 1.5 T on the heart, with 4 to 5 slices acquired per heartbeat. Respiratory motion was compensated using navigator-based slice tracking. Residual in-plane motion and related magnetic susceptibility artifacts were corrected online. The standard deviation of temperature was measured on healthy volunteers (N = 5) in both ventricles. On animals, the MR-compatible catheter was positioned and visualized in the left ventricle (LV) using a bSSFP pulse sequence with active catheter tracking. Twelve MR-guided RFA were performed on three sheep in vivo at various locations in left ventricle (LV). The dimensions of the thermal lesions measured on thermal dose images, on 3D T1-weighted (T1-w) images acquired immediately after the ablation and at gross pathology were correlated.

RESULTS

MR thermometry uncertainty was 1.5 °C on average over more than 96% of the pixels covering the left and right ventricles, on each volunteer. On animals, catheter repositioning in the LV with active slice tracking was successfully performed and each ablation could be monitored in real-time by MR thermometry and thermal dosimetry. Thermal lesion dimensions on TD maps were found to be highly correlated with those observed on post-ablation T1-w images (R = 0.87) that also correlated (R = 0.89) with measurements at gross pathology.

CONCLUSIONS

Quantitative TD mapping from real-time rapid CMR thermometry during catheter-based RFA is feasible. It provides a direct assessment of the lesion extent in the myocardium with precision in the range of one millimeter. Real-time MR thermometry and thermal dosimetry may improve safety and efficacy of the RFA procedure by offering a reliable indicator of therapy outcome during the procedure.

Kalman Filtered Bio Heat Transfer Model Based Self-adaptive Hybrid Magnetic Resonance Thermometry

To develop a self-adaptive and fast thermometry method by combining the original hybrid magnetic resonance thermometry method and the bio heat transfer equation (BHTE) model. The proposed Kalman filtered Bio Heat Transfer Model Based Self-adaptive Hybrid Magnetic Resonance Thermometry, abbreviated as KalBHT hybrid method, introduced the BHTE model to synthesize a window on the regularization term of the hybrid algorithm, which leads to a self-adaptive regularization both spatially and temporally with change of temperature. Further, to decrease the sensitivity to accuracy of the BHTE model, Kalman filter is utilized to update the window at each iteration time. To investigate the effect of the proposed model, computer heating simulation, phantom microwave heating experiment and dynamic in-vivo model validation of liver and thoracic tumor were conducted in this study. The heating simulation indicates that the KalBHT hybrid algorithm achieves more accurate results without adjusting λ to a proper value in comparison to the hybrid algorithm. The results of the phantom heating experiment illustrate that the proposed model is able to follow temperature changes in the presence of motion and the temperature estimated also shows less noise in the background and surrounding the hot spot. The dynamic in-vivo model validation with heating simulation demonstrates that the proposed model has a higher convergence rate, more robustness to susceptibility problem surrounding the hot spot and more accuracy of temperature estimation. In the healthy liver experiment with heating simulation, the RMSE of the hot spot of the proposed model is reduced to about 50% compared to the RMSE of the original hybrid model and the convergence time becomes only about one fifth of the hybrid model. The proposed model is able to improve the accuracy of the original hybrid algorithm and accelerate the convergence rate of MR temperature estimation.

Accurate MR thermometry by hyperpolarized 129 Xe

PURPOSE

To investigate the temperature dependence of the resonance frequency of lipid-dissolved xenon (LDX) and to assess the accuracy of LDX-based MR thermometry.

METHODS

The chemical shift temperature dependence of water protons, methylene protons, and LDX was measured from samples containing tissues with varying fat contents using a high-resolution NMR spectrometer. LDX results were then used to acquire relative and absolute temperature maps in vivo and the results were compared with PRF-based MR thermometry.

RESULTS

The temperature dependence of proton resonance frequency (PRF) is strongly affected by the specific distribution of water and fat. A redistribution of water and fat compartments can reduce the apparent temperature dependence of the water chemical shift from -0.01 ppm/°C to -0.006 ppm, whereas the LDX chemical shift shows a consistent temperature dependence of -0.21 ppm/°C. The use of the methylene protons resonance frequency as internal reference improves the accuracy of LDX-based MR thermometry, but degrades that of PRF-based MR thermometry, as microscopic susceptibility gradients affected lipid and water spins differently.

CONCLUSION

The LDX resonance frequency, with its higher temperature dependence, provides more accurate and precise temperature measurements, both in vitro and in vivo. More importantly, the resonance frequency of nearby methylene protons can be used to extract absolute temperature information. Magn Reson Med 78:1070-1079, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Open-source, small-animal magnetic resonance-guided focused ultrasound system

BACKGROUND

MR-guided focused ultrasound or high-intensity focused ultrasound (MRgFUS/MRgHIFU) is a non-invasive therapeutic modality with many potential applications in areas such as cancer therapy, drug delivery, and blood-brain barrier opening. However, the large financial costs involved in developing preclinical MRgFUS systems represent a barrier to research groups interested in developing new techniques and applications. We aim to mitigate these challenges by detailing a validated, open-source preclinical MRgFUS system capable of delivering thermal and mechanical FUS in a quantifiable and repeatable manner under real-time MRI guidance.

METHODS

A hardware and software package was developed that includes closed-loop feedback controlled thermometry code and CAD drawings for a therapy table designed for a preclinical MRI scanner. For thermal treatments, the modular software uses a proportional integral derivative controller to maintain a precise focal temperature rise in the target given input from MR phase images obtained concurrently. The software computes the required voltage output and transmits it to a FUS transducer that is embedded in the delivery table within the magnet bore. The delivery table holds the FUS transducer, a small animal and its monitoring equipment, and a transmit/receive RF coil. The transducer is coupled to the animal via a water bath and is translatable in two dimensions from outside the magnet. The transducer is driven by a waveform generator and amplifier controlled by real-time software in Matlab. MR acoustic radiation force imaging is also implemented to confirm the position of the focus for mechanical and thermal treatments.

RESULTS

The system was validated in tissue-mimicking phantoms and in vivo during murine tumor hyperthermia treatments. Sonications were successfully controlled over a range of temperatures and thermal doses for up to 20 min with minimal temperature overshoot. MR thermometry was validated with an optical temperature probe, and focus visualization was achieved with acoustic radiation force imaging.

CONCLUSIONS

We developed an MRgFUS platform for small-animal treatments that robustly delivers accurate, precise, and controllable sonications over extended time periods. This system is an open source and could increase the availability of low-cost small-animal systems to interdisciplinary researchers seeking to develop new MRgFUS applications and technology.

Procedural sedation and analgesia for respiratory-gated MR-HIFU in the liver: a feasibility study

BACKGROUND

Previous studies demonstrated both pre-clinically and clinically the feasibility of magnetic resonance-guided high-intensity focused ultrasound (MR-HIFU) ablations in the liver. To overcome the associated problem of respiratory motion of the ablation area, general anesthesia (GA) and mechanical ventilation was used in conjunction with either respiratory-gated energy delivery or energy delivery during induced apnea. However, clinical procedures requiring GA are generally associated with increased mortality, morbidity, and complication rate compared to procedural sedation and analgesia (PSA). Furthermore, PSA is associated with faster recovery and an increased eligibility for non- and mini-invasive interventions.

METHODS

In this study, we investigate both in an animal model and on a small patient group the kinetics of the diaphragm during free-breathing, when a tailored remifentanil/propofol-based PSA protocol inducing partial respiratory depression is used. Subsequently, we demonstrate in an animal study the compatibility of the resulting respiratory pattern of the PSA protocol with a gated HIFU ablation in the liver by direct comparison with gated ablations conducted under GA. Wilcoxon signed-rank tests were performed for statistical analysis of non-perfused and necrosed tissue volumes. Duty cycles (ratio or percentage of the breathing cycle with the diaphragm in its resting position, such that acoustic energy delivery with MR-HIFU was allowed) were statistically compared for both GA and PSA using student's t tests.

RESULTS

In both animal and human experiments, the breathing frequency was decreased below 9/min, while maintaining stable vital functions. Furthermore an end-exhalation resting phase was induced by this PSA protocol during which the diaphragm is virtually immobile. Median non-perfused volumes, non-viable volumes based on NADH staining, and duty cycles were larger under PSA than under GA or equal.

CONCLUSIONS

We conclude that MR-HIFU ablations of the liver under PSA are feasible and potentially increase the non-invasive nature of this type of intervention.

Drift correction for accurate PRF-shift MR thermometry during mild hyperthermia treatments with MR-HIFU

There is growing interest in performing hyperthermia treatments with clinical magnetic resonance imaging-guided high-intensity focused ultrasound (MR-HIFU) therapy systems designed for tissue ablation. During hyperthermia treatment, however, due to the narrow therapeutic window (41-45 °C), careful evaluation of the accuracy of proton resonant frequency (PRF) shift MR thermometry for these types of exposures is required.

PURPOSE

The purpose of this study was to evaluate the accuracy of MR thermometry using a clinical MR-HIFU system equipped with a hyperthermia treatment algorithm.

METHODS

Mild heating was performed in a tissue-mimicking phantom with implanted temperature sensors using the clinical MR-HIFU system. The influence of image-acquisition settings and post-acquisition correction algorithms on the accuracy of temperature measurements was investigated. The ability to achieve uniform heating for up to 40 min was evaluated in rabbit experiments.

RESULTS

Automatic centre-frequency adjustments prior to image-acquisition corrected the image-shifts in the order of 0.1 mm/min. Zero- and first-order phase variations were observed over time, supporting the use of a combined drift correction algorithm. The temperature accuracy achieved using both centre-frequency adjustment and the combined drift correction algorithm was 0.57° ± 0.58 °C in the heated region and 0.54° ± 0.42 °C in the unheated region.

CONCLUSION

Accurate temperature monitoring of hyperthermia exposures using PRF shift MR thermometry is possible through careful implementation of image-acquisition settings and drift correction algorithms. For the evaluated clinical MR-HIFU system, centre-frequency adjustment eliminated image shifts, and a combined drift correction algorithm achieved temperature measurements with an acceptable accuracy for monitoring and controlling hyperthermia exposures.

Fast temperature estimation from undersampled k-space with fully-sampled center for MR guided microwave ablation

PURPOSE

This study aims to accelerate MR temperature imaging using the proton resonance frequency (PRF) shift method for real time temperature monitoring during thermal ablation.

MATERIALS AND METHODS

The proposed method estimates temperature changes from undersampled k-space with a fully sampled center. This proposed algorithm is based on the hybrid multi-baseline and referenceless treatment image model and can be seen as an extension of the conventional k-space-based hybrid thermometry. The parameters of hybrid model are acquired by utilizing information from low resolution images which are obtained from fully-sampled centers of k-space. Registration is used to correct temperature errors due to the displacement of the subject. Phantom heating simulations, motion simulations, phantom heating and in-vivo experiments were performed to investigate the efficiency of the proposed method. SPIRiT and the conventional k-space estimation reconstruction thermometry were implemented for comparison using the same sampling pattern.

RESULTS

The phantom heating simulations showed that the proposed method results in lower RMSEs than the conventional k-space hybrid thermometry and SPIRiT at various reduction factors tested. The motion simulations indicated the robustness of the proposed method to displacement of the subject. Phantom heating experiment further demonstrated the ability of the method to reconstruct temperature maps with less computation time and higher accuracy (RMSEs lower than 0.4°C) at a net reduction factor of 3.5 in the presence of large noise caused by a microwave needle. In-vivo experiments validated the feasibility of the proposed method to estimate temperature changes from undersampled k-space (net reduction factor 4.3) in presence of respiratory motion and complicated anatomical structure, while reducing computation time as much as 10-fold compared with the conventional k-space method.

CONCLUSION

The proposed method accelerates the PRF-shift MR thermometry and provides more accurate temperature maps in presence of motion with relatively short computation time, which may make real time imaging for MR-guided microwave ablation possible.

One-lung flooding reduces the ipsilateral diaphragm motion during mechanical ventilation

BACKGROUND

Diaphragm motion during spontaneous or mechanical respiration hinders image-guided percutaneous interventions of tumours in lung and upper abdomen. Motion-tracking methods can be applied but increase procedure complexity and procedure time. One-lung flooding (OLF) generates a suitable acoustic pathway to lung tumours and likely suppress diaphragm motion. The aim of this study was to quantify the effect of OLF on ipsilateral diaphragm motion during contralateral one-lung ventilation.

METHODS

To measure the diaphragm motion, M-mode ultrasonography of the right hemidiaphragm was performed during spontaneous breathing and mechanical ventilation, as well as after right-side lung flooding, in three pigs. Diaphragm motion was analysed using magnetic resonance images during left-side lung flooding and mechanical ventilation, in four pigs.

RESULTS

Double-lung ventilation increased the diaphragm movement in comparison with spontaneous breathing (17.8 ± 4.4 vs. 12.2 ± 3.4 mm, p = 0.014). Diaphragm movement on the flooded side during contralateral one-lung ventilation was significantly reduced compared to that during double-lung ventilation (3.9 ± 1.0 vs. 17.8 ± 4.4 mm, p = 0.041). By analysing the magnetic resonance images, the hemidiaphragm on the flooded side showed an average displacement of 4.2 mm, a maximum displacement of 15 mm close to the ventilated lung and no displacement at the lateral side.

CONCLUSION

OLF leads to a drastic reduction of diaphragm motion on the ipsilateral side which implies that targeting and motion compensation algorithms for interventions like high-intensity focused ultrasound ablation of intrapulmonary and hepatic lesions might not be required.

Improved cardiac magnetic resonance thermometry and dosimetry for monitoring lesion formation during catheter ablation

PURPOSE

A new real-time MR-thermometry pipeline was developed to measure multiple temperature images per heartbeat with 1.6×1.6×3 mm³ spatial resolution. The method was evaluated on 10 healthy volunteers and during radiofrequency ablation (RFA) in sheep.

METHODS

Multislice, electrocardiogram-triggered, echo-planar imaging was combined with parallel imaging, under free breathing conditions. In-plane respiratory motion was corrected on magnitude images by an optical flow algorithm. Motion-related susceptibility artifacts were compensated on phase images by an algorithm based on Principal Component Analysis. Correction of phase drift and temporal filter were included in the pipeline implemented in the Gadgetron framework. Contact electrograms were recorded simultaneously with MR thermometry by an MR-compatible ablation catheter.

RESULTS

The temporal standard deviation of temperature in the left ventricle remained below 2 °C on each volunteer. In sheep, focal heated regions near the catheter tip were observed on temperature images (maximal temperature increase of 38 °C) during RFA, with contact electrograms of acceptable quality. Thermal lesion dimensions at gross pathology were in agreement with those observed on thermal dose images.

CONCLUSION

This fully automated MR thermometry pipeline (five images/heartbeat) provides direct assessment of lesion formation in the heart during catheter-based RFA, which may improve treatment of cardiac arrhythmia by ablation. Magn Reson Med 77:673-683, 2017. © 2016 International Society for Magnetic Resonance in Medicine.