Irreversible change in the T1 temperature dependence with thermal dose using the proton resonance frequency-T1 technique

Simultaneous proton resonance frequency T1 - MR shear wave elastography for MR-guided focused ultrasound multiparametric treatment monitoring

PURPOSE

To develop an efficient MRI pulse sequence to simultaneously measure multiple parameters that have been shown to correlate with tissue nonviability following thermal therapies.

METHODS

A 3D segmented EPI pulse sequence was used to simultaneously measure proton resonance frequency shift (PRFS) MR thermometry (MRT), T1 relaxation time, and shear wave velocity induced by focused ultrasound (FUS) push pulses. Experiments were performed in tissue mimicking gelatin phantoms and ex vivo bovine liver. Using a carefully designed FUS triggering scheme, a heating duty cycle of approximately 65% was achieved by interleaving FUS ablation pulses with FUS push pulses to induce shear waves in the tissue.

RESULTS

In phantom studies, temperature increases measured with PRFS MRT and increases in T1 correlated with decreased shear wave velocity, consistent with material softening with increasing temperature. During ablation in ex vivo liver, temperature increase measured with PRFS MRT initially correlated with increasing T1 and decreasing shear wave velocity, and after tissue coagulation with decreasing T1 and increasing shear wave velocity. This is consistent with a previously described hysteresis in T1 versus PRFS curves and increased tissue stiffness with tissue coagulation.

CONCLUSION

An efficient approach for simultaneous and dynamic measurements of PRSF, T1 , and shear wave velocity during treatment is presented. This approach holds promise for providing co-registered dynamic measures of multiple parameters, which correlates to tissue nonviability during and following thermal therapies, such as FUS.

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.

A T1-based correction method for proton resonance frequency shift thermometry in breast tissue

PURPOSE

Develop and evaluate the effectiveness of a T1-based correction method for errors in proton resonant frequency shift thermometry due to non-local field effects caused by heating in fatty breast tissues.

METHODS

Computational models of human breast tissue were created by segmenting MRI data from a healthy human volunteer. MR-guided focused ultrasound (MRgFUS) heating and MR thermometry measurements were simulated in several locations in the heterogeneous segmented breast models. A T1-based correction method for PRF thermometry errors was applied and the maximum positive and negative errors and the root mean squared error (RMSE) in a region around each heating location was evaluated with and without correction. The method uses T1 measurements to estimate the temperature change in fatty tissues and correct for their influence. Experimental data from a heating study in cadaver breast tissue were analyzed, and the expected PRFS error computed.

RESULTS

The simulated MR thermometry had maximum single voxel errors ranging between 10% and 18% when no correction was applied. Applying the correction led to a considerable improvement, lowering the maximum error range to 2%-5%. The 5th to 95th percentile interval of the temperature error distribution was also lowered with correction, from approximately 3.5 to 1°C. This correction worked even when T1 times were uniformly raised or lowered by 5%-10%. The experimental data showed predicted errors of 15%.

CONCLUSIONS

This simulation study demonstrates that the T1-based correction method reduces MR thermometry errors due to non-local effects from heating in fatty tissues, potentially improving the accuracy of thermometry measurements during MRgFUS treatments. The presented correction method is reliant on having a patient-specific 3D model of the breast, and may be limited by the accuracy of the fat temperatures which in turn may be limited by noise or bias present in the T1 measurements.

Iron oxide magnetic nanoparticles based low-field MR thermometry

This paper reports on a highly accurate approach of magnetic resonance (MR) thermometry using iron oxide magnetic nanoparticles (MNPs) as temperature sensors. An empirical model for the description of the temperature dependent R ₂ relaxation rate is proposed by taking into account the temperature sensitivity of the MNP magnetization. The temperature sensitivity of the MNP magnetization (η) and the temperature sensitivity of the R ₂ relaxation rate (κ) are simulated with the proposed empirical models to investigate their dependence on the magnetic field and the particle size. Simulation results show the existence of optimal magnetic fields H_oη and H_oκ that maximize the temperature sensitivities η and κ. Furthermore, simulations and experiments demonstrate that the optimal magnetic field H_oη (H_oκ ) decreases with increasing the particle size. Experiments on temperature dependent R ₂ relaxation rate are performed at different magnetic fields for MNP samples with different iron concentrations. Experimental results show that the proposed MR thermometry using MNPs as temperature sensors allows a temperature estimation accuracy of about 0.05 °C. We believe that the achieved approach of highly accurate MR thermometry is of great interest and significance to biomedicine and biology.

2D multislice MP2RAGE sequence for fast T1 mapping at 7 T: Application to mouse imaging and MR thermometry

PURPOSE

To develop a 2D radial multislice MP2RAGE sequence for fast and reliable T₁ mapping at 7 T in mice and for MR thermometry.

METHODS

The 2D-MP2RAGE sequence was performed with the following parameters: TI₁ -TI₂ -MP2RAGE_TR = 1000-3000-9000 ms. The multiple dead times within the sequence were used for interleaved multislice acquisition, enabling one to acquire six slices in 9 seconds. The excitation pulse shape, inversion selectivity, and interslice gap were optimized. In vitro comparison with the inversion-recovery sequence was performed. The T₁ variations with temperature were measured on tubes with T₁ ranging from 800 ms to 2000 ms. The sequence was used to acquire T₁ maps continuously during 30 minutes on the brain and abdomen of healthy mice.

RESULTS

A three-lobe cardinal sine excitation pulse, combined with an inversion slice thickness and an interslice gap of respectively 150% and 50% of the imaging slice thickness, led to a SD and bias of the T₁ measurements below 1% and 2%, respectively. A linear dependence of T₁ with temperature was measured between 10°C and 60°C. In vivo, less than 1% variation was measured between successive T₁ maps in the mouse brain. In the abdomen, no obvious in-plane motion artifacts were observed but respiratory motion in the slice dimension led to 6% T₁ underestimation.

CONCLUSION

The multislice MP2RAGE sequence could be used for fast whole-body T₁ mapping and MR thermometry. Its reconstruction method would enable on-the-fly reconstruction.

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.

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.

Simultaneous proton resonance frequency shift thermometry and T1 measurements using a single reference variable flip angle T1 method

PURPOSE

Implement simultaneous proton resonance frequency (PRF) shift and T₁ measurements with equivalent temporal resolution using a single reference variable flip angle method. This novel method allows for simultaneous thermometry in both aqueous and fatty tissue.

METHODS

This method acquires a single reference image at the lower flip angle and all dynamic images at the higher angle. T₁ is calculated using a single reference variable flip angle method, which accounts for the reference image temperature remaining constant. Monte Carlo simulations determined the optimal dynamic flip angle for combined PRF and T₁ measurements. This method was evaluated in MR-guided focused ultrasound heating experiments using a gelatin phantom and human cadaver breasts. In vivo measurement precision was demonstrated in healthy female volunteers under nonheating conditions.

RESULTS

Temperature rise during MR-guided focused ultrasound heating was measured in aqueous tissue with both PRF and T₁ . Both measures show good qualitative agreement in both space and time in aqueous tissue. The T₁ change due to temperature increase was measured in fat, demonstrating the expected temporal response. The dynamic flip angle that produces optimal SNR for PRF measurements is lower than the optimal angle for T₁ measurements, necessitating the selection of a compromise angle.

CONCLUSION

The single reference variable flip angle method provides a reliable way to simultaneously measure PRF temperature and T₁ change and overcomes PRF's inability to simultaneously monitor temperature in aqueous and adipose tissues. Future work will calibrate T₁ change to temperature, enabling real-time temperature in fat and increasing patient safety and treatment efficacy during thermal interventional treatments.

New experimental model for single liver lobe hyperthermia in small animals using non-directional microwaves

PURPOSE

Our aim was to develop a new experimental model for in vivo hyperthermia using non-directional microwaves, applicable to small experimental animals. We present an affordable approach for targeted microwave heat delivery to an isolated liver lobe in rat, which allows rapid, precise and stable tissue temperature control.

MATERIALS AND METHODS

A new experimental model is proposed. We used a commercial available magnetron generating 2450 MHz, with 4.4V and 14A in the filament and 4500V anodic voltage. Modifications were required in order to adjust tissue heating such as to prevent overheating and to allow for fine adjustments according to real-time target temperature. The heating is controlled using a virtual instrument application implemented in LabView® and responds to 0.1° C variations in the target. Ten healthy adult male Wistar rats, weighing 250-270 g were used in this study. The middle liver lobe was the target for controlled heating, while the rest of the living animal was protected.

RESULTS

In vivo microwave delivery using our experimental setting is safe for the animals. Target tissue temperature rises from 30°C to 40°C with 3.375°C / second (R2 = 0.9551), while the increment is lower it the next two intervals (40-42°C and 42-44°C) with 0.291°C/ s (R2 = 0.9337) and 0.136°C/ s (R2 = 0.7894) respectively, when testing in sequences. After reaching the desired temperature, controlled microwave delivery insures a very stable temperature during the experiments.

CONCLUSIONS

We have developed an inexpensive and easy to manufacture system for targeted hyperthermia using non-directional microwave radiation. This system allows for fine and stable temperature adjustments within the target tissue and is ideal for experimental models testing below or above threshold hyperthermia.

Technical Note: The effect of 2D excitation profile on T1 measurement accuracy using the variable flip angle method with an average flip angle assumption

PURPOSE

To study the accuracy and precision of T₁ estimates using the Variable Flip Angle (VFA) method in 2D and 3D acquisitions.

METHODS

Excitation profiles were simulated using numerical implementation of the Bloch equations for Hamming-windowed sinc excitation pulses with different time-bandwidth products (TBP) of 2, 6, and 10 and for T₁ values of 295 ms and 1045 ms. Experimental data were collected in 5° increments from 5° to 90° for the same T₁ and TBP values. T₁ was calculated for every combination of flip angle with and without a correction for B₁ and slice profile variation. Calculations were also made for flat slice profile such as obtained in 3D acquisition. Monte Carlo simulations were performed to obtain T₁ measurement uncertainty.

RESULTS

VFA T₁ measurements in 2D without correction can result in a 40-80% underestimation of true T₁ . Flip angle correction can reduce the underestimation, but results in accurate measurements of T₁ only within a narrow band of flip angle combinations. The narrow band of accuracy increases with TBP, but remains too narrow for any practical range of T₁ values or B₁ variation. Simulated noisy VFA T₁ measurements in 3D were accurate as long as the two angles chosen are on either side of the Ernst angle.

CONCLUSIONS

Accurate T1 estimates from VFA 2D acquisitions are possible, but only a narrow range of T1 values within a narrow range of flip angle combinations can be accurately calculated using a 2D slice. Unless a better flip angle correction method is used, these results demonstrate that accurate measurements of T1 in 2D cannot be obtained robustly enough for practical use and are more likely obtained by a thin slab 3D VFA acquisition than from multiple-slice 2D acquisitions. VFA T₁ measurements in 3D are accurate for wide ranges of flip angle combinations and T₁ values.

Optothermal profile of an ablation catheter with integrated microcoil for MR-thermometry during Nd:YAG laser interstitial thermal therapies of the liver—an in-vitro experimental and theoretical study

PURPOSE

Flexible microcoils integrated with ablation catheters can improve the temperature accuracy during local MR-thermometry in Nd:YAG laser interstitial thermal therapies. Here, the authors are concerned with obtaining a preliminary confirmation of the clinical utility of the modified catheter. They investigate whether the thin-film substrate and copper tracks of the printed coil inductor affect the symmetry of the thermal profile, and hence of the lesion produced.

METHODS

Transmission spectroscopy in the near infrared was performed to test for the attenuation at 1064 nm through the 25 μm thick Kapton substrate of the microcoil. The radial transmission profile of an infrared high-power, light emitting diode with >80% normalized power at 1064 nm was measured through a cross section of the modified applicator to assess the impact of the copper inductor on the optical profile. The measurements were performed in air, as well as with the applicator surrounded by two types of scattering media; crystals of NaCl and a layer of liver-mimicking gel phantom. A numerical model based on Huygens-Fresnel principle and finite element simulations, using a commercially available package (COMSOL Multiphysics), were employed to compare with the optical measurements. The impact of the modified optical profile on the thermal symmetry was assessed by examining the high resolution microcoil derived thermal maps from a Nd:YAG laser ablation performed on a liver-mimicking gel phantom.

RESULTS

Less than 30% attenuation through the Kapton film was verified. Shadowing behind the copper tracks was observed in air and the measured radial irradiation correlated well with the diffraction pattern calculated numerically using the Huygens-Fresnel principle. Both optical experiments and simulations, demonstrate that shadowing is mitigated by the scattering properties of a turbid medium. The microcoil derived thermal maps at the end of a Nd:YAG laser ablation performed on a gel phantom in a 3 T scanner confirm that the modified irradiation pattern does not disrupt the thermal symmetry, even though, unlike tissue, the gel is minimally scattering.

CONCLUSIONS

The results from this initial assessment indicate that microcoils can be safely integrated with ablation catheters and ensure that the complete necrosis of the liver tumor can still be achieved.

Treatment envelope evaluation in transcranial magnetic resonance-guided focused ultrasound utilizing 3D MR thermometry

BACKGROUND

Current clinical targets for transcranial magnetic resonance-guided focused ultrasound (tcMRgFUS) are all located close to the geometric center of the skull convexity, which minimizes challenges related to focusing the ultrasound through the skull bone. Non-central targets will have to be reached to treat a wider variety of neurological disorders and solid tumors. Treatment envelope studies utilizing two-dimensional (2D) magnetic resonance (MR) thermometry have previously been performed to determine the regions in which therapeutic levels of FUS can currently be delivered. Since 2D MR thermometry was used, very limited information about unintended heating in near-field tissue/bone interfaces could be deduced.

METHODS

In this paper, we present a proof-of-concept treatment envelope study with three-dimensional (3D) MR thermometry monitoring of FUS heatings performed in a phantom and a lamb model. While the moderate-sized transducer used was not designed for transcranial geometries, the 3D temperature maps enable monitoring of the entire sonication field of view, including both the focal spot and near-field tissue/bone interfaces, for full characterization of all heating that may occur. 3D MR thermometry is achieved by a combination of k-space subsampling and a previously described temporally constrained reconstruction method.

RESULTS

We present two different types of treatment envelopes. The first is based only on the focal spot heating-the type that can be derived from 2D MR thermometry. The second type is based on the relative near-field heating and is calculated as the ratio between the focal spot heating and the near-field heating. This utilizes the full 3D MR thermometry data achieved in this study.

CONCLUSIONS

It is shown that 3D MR thermometry can be used to improve the safety assessment in treatment envelope evaluations. Using a non-optimal transducer, it is shown that some regions where therapeutic levels of FUS can be delivered, as suggested by the first type of envelope, are not necessarily safely treated due to the amount of unintended near-field heating occurring. The results presented in this study highlight the need for 3D MR thermometry in tcMRgFUS.

In vivo evaluation of multi-echo hybrid PRF/T1 approach for temperature monitoring during breast MR-guided focused ultrasound surgery treatments

PURPOSE

To evaluate the precision of in vivo temperature measurements in adipose and glandular breast tissue using a multi-echo hybrid PRF/T1 pulse sequence.

METHODS

A high-bandwidth, multi-echo hybrid PRF/T1 sequence was developed for monitoring temperature changes simultaneously in fat- and water-based tissues. The multiple echoes were combined with the optimal weightings for magnitude and phase images, allowing for precise measurement of both T1 and the proton resonance frequency (PRF) shift. The sequence was tested during in vivo imaging of 10 healthy volunteers in a breast-specific MR-guided focused ultrasound system and also during focused ultrasound heating of excised breast adipose tissue.

RESULTS

The in vivo results indicated that the sequence can measure PRF temperatures with 1.25 × 1.25 × 3.5 mm resolution, 1.9 s temporal resolution, and 1.0°C temperature precision, and can measure T1 values with 3.75 × 3.75 × 3.5 mm resolution, 3.8 s temporal resolution, and 2.5%-4.8% precision. The excised tissue heating experiments demonstrate the sequence's ability to monitor temperature changes simultaneously in water- and fat-based tissues.

CONCLUSION

The addition of a high-bandwidth, multi-echo readout to the hybrid PRF/T1 sequence improves the precision of each measurement, providing a sequence that will be beneficial to several MR-guided thermal therapies.

Phase reconstruction from multiple coil data using a virtual reference coil

PURPOSE

This study develops a method to obtain optimal estimates of absolute magnetization phase from multiple-coil MRI data.

THEORY AND METHODS

The element-specific phases of a multi-element receiver coil array are accounted for by using the phase of a real or virtual reference coil that is sensitive over the entire imaged volume. The virtual-reference coil is generated as a weighted combination of measurements from all receiver coils. The phase-corrected multiple coil complex images are combined using the inverse covariance matrix. These methods are tested on images of an agar phantom, an in vivo breast, and an anesthetized rabbit obtained using combinations of four, nine, and three receiver channels, respectively.

RESULTS

The four- and three-channel acquisitions require formation of a virtual-reference receiver coil while one channel of the nine-channel receive array has a sensitivity profile covering the entire imaged volume. Referencing to a real or virtual coil gives receiver phases that are essentially identical except for the individual receiver channel noise. The resulting combined images, which account for receiver channel noise covariance, show the expected reduction in phase variance.

CONCLUSION

The proposed virtual reference coil method determines a phase distribution for each coil from which an optimal phase map can be obtained.