Scholarly article on topic 'Treatment envelope evaluation in transcranial magnetic resonance-guided focused ultrasound utilizing 3D MR thermometry'

Treatment envelope evaluation in transcranial magnetic resonance-guided focused ultrasound utilizing 3D MR thermometry Academic research paper on "Medical engineering"

Share paper
Academic journal
Journal of Therapeutic Ultrasound
OECD Field of science

Academic research paper on topic "Treatment envelope evaluation in transcranial magnetic resonance-guided focused ultrasound utilizing 3D MR thermometry"


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

Henrik Odeen1,2*, Joshua de Bever3, Scott Almquist3, Alexis Farrer4, Nick Todd1, Allison Payne1, John W Snell5,6, Douglas A Christensen4,7 and Dennis L Parker1

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.

Keywords: Treatment envelope, tcMRgFUS, MR thermometry, PRF, Brain

* Correspondence:

1Utah Center for Advanced Imaging Research, Department of Radiology, University of Utah, Salt Lake City, Utah 84108, USA

2Department of Physics and Astronomy, University of Utah, Salt Lake City, Utah 84112, USA

Fulllist of author information is available at the end of the article


Bio Med Central

© 2014 Odeen et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.


Focused ultrasound (FUS) is a promising treatment modality for a wide range of disorders due to its precise, non-invasive, and ionizing radiation-free nature. It is currently being performed under ultrasound [1-4] (US) and magnetic resonance (MR) guidance [5-7]. Both thermal and mechanical [8-11] effects of the US interaction with tissues are being utilized, and current targets of clinical interest include the prostate [12-14], liver [15-17], breast [18-20], uterine fibroid [21,22], and brain [5,23-25].

Neurological applications of FUS have been of particular research interest due to the unique possibility of non-invasive delivery of US energy through the intact skull bone to a precise location within the brain [26-28]. Irreversible surgical procedures utilizing high US intensities for tissue ablation [24,29,30] as well as reversible neuro-modulation and blood-brain barrier disruption procedures utilizing lower US intensities [31-35] have been investigated. Challenges to transcranial FUS have included the defocusing of the US beam due to phase aberrations by the skull bone [36], and the difficulty of non-invasively monitoring the procedures. Recent approaches to phase aberration correction [37-43] have greatly improved the ability to maintain a localized US focal spot even after propagation through the skull. Non-invasive treatment monitoring has been made possible with the use of magnetic resonance imaging (MRI). MRI offers the ability to acquire high-resolution anatomical images with various contrasts in arbitrary imaging planes, real-time feedback from MR thermometry, as well as methods for post-treatment evaluation [44-46]. For MR temperature imaging (MRTI), various temperature-sensitive MR parameters such as the T1 and T2 relaxation times can be utilized [47,48]; however, the current gold standard for temperature measurements in aqueous tissues is the proton resonance frequency (PRF) shift method [49,50]. The primary challenge in MR thermometry for transcranial applications is achieving accurate and precise temperature measurements with adequate spatial and temporal resolution over a sufficiently large field of view (FOV). Typical MRTI methods cover only one or a few two-dimensional (2D) slices in order to get high enough spatio-temporal resolution for accurate temperature measurements [23,24].

Despite recent improvements in phase aberration correction and non-invasive treatment monitoring, current clinical targets in transcranial MR-guided FUS (tcMRgFUS)

are all close to the geometric center of the skull convexity [23,24]. For current clinical tcMRgFUS systems, which utilize large hemispherical phased-array US transducers, this minimizes the challenges of focusing through the skull, as most of the ultrasound elements will be nearly orthogonal to the skull surface. Central focusing also results in less unintended off-focal heating (i.e., heating in near- and far-field tissue/bone interfaces, and heating due to grating lobes) than may occur when focusing close to the skull. However, to be able to treat a wider variety of neurological disorders and solid tumors, therapeutic levels of FUS will have to be delivered to non-central targets within the skull as well. This will lead to more challenging beam aberrations and more unintended off-focal heating. To investigate in what regions of the brain therapeutic levels of FUS can be delivered, Eames et al. [51] and Dallapiazza et al. [52] performed treatment envelope studies where the US focus was mechanically and electronically moved to various locations within human cadaver heads and phantom skulls while monitoring the heating with 2D MRTI.

In this work, we present methods to perform treatment envelope studies utilizing three-dimensional (3D) MRTI by combining k-space subsampling and a previously described constrained reconstruction method. The utilization of 3D MRTI enables monitoring of the fully insonified FOV and thereby simultaneous monitoring of the focal spot and any unintended near-field heating. The methodology of this proof-of-concept study is evaluated in both a phantom and a lamb model, and it is shown that although therapeutic levels of US can be delivered to an intracranial region, the amount of unintended heating occurring in the near-field may make it unsafe to do so.


All MR imaging was performed with a 3D gradient recalled echo pulse sequence with segmented echo planar imaging readout on a 3 T MRI scanner (TIM Trio, Siemens Medical Solutions, Erlangen, Germany). All imaging parameters used for the phantom and the lamb studies are given in Table 1. The repetition time (TR) was longer for the lamb study since a fat-saturation pre-pulse was used in this study, but not in the phantom study. k-space was evenly subsampled in the ky-phase encode direction, while the kx-read out and the kz-slice encode directions were both fully sampled in each

Table 1 MR scan parameters for the phantom and lamb studies, respectively

TR/TE (ms) FOV (mm) Voxel size (mm) BW (Hz/px) ETL ES (ms) FA (°) tacq (s) R

Phantom 25/11 192x192x110 1.2x1.2x2.5 744 9 1.59 15 3.30 6

Lamb 36/11 192x144x60 1.0x1.0x2.0 752 9 1.61 20 4.32 4

TR repetition time, TE echo time, FOV field of view, BW bandwidth (in readout direction), ETL echo train length, ES echo spacing, FA flip angle, taeq acquisition time, R k-space reduction/subsampling factor.

dynamic time frame. All subsampled data were reconstructed with a previously described temporally constrained reconstruction (TCR) algorithm [53,54]. In the TCR algorithm, the series of temperature images are reconstructed by iteratively minimizing a cost function consisting of a temporal constraint term and a data fidelity term using a steepest descent method, according to

m = argminml ||WFm'-d!

phase difference between two images is scaled according to

AT = -


where m are the artifact-free images, W is a binary matrix representing sampled phase encode lines, F is the fast Fourier transform, m' is the image estimate, d' is the acquired subsampled k-space data, a is a spatially varying weighting factor, the sum is over all N pixels in the data set, Vt is the temporal gradient operator, and II * ll2 is the L2 norm. Following TCR reconstruction, all data were zero-filled interpolated to 0.5-mm in-plane resolution and 1.0-mm slice thickness to minimize partial volume effects. MR temperature measurements were performed with the PRF shift method [49], where the

where AT is the change in temperature (°C), Ap is the phase difference between the two images (radians), y is the gyromagnetic ratio (Hz/T), a is a thermal coefficient (here assumed to be -0.01 ppm/°C), B0 is the main magnetic field strength (T), and TE is the echo time (s). All image reconstruction and post-processing was performed using MATLAB (R2013a, The MathWorks Inc., Natick, MA, USA).

For the phantom study, two in-house-built radio-frequency (RF) receive-only single-channel loop coils were used for signal detection, and for the lamb study, two in-house-built dual-channel receive-only RF coils [55] were used (Figure 1a,e).

All FUS heating was performed with an MR-compatible phased-array ultrasound transducer (256 elements, 1-MHz frequency, 13-cm radius of curvature, focal spot full-width- at-half maximum 2x2x8 mm, Imasonic, Voray-sur-l'Ognon, France) and accompanying hardware and software for mechanical positioning and electronic beam steering (Image Guided Therapy, Pessac, France). The

Figure 1 Experimental setup. (a) Coronalview of the experimentalsetup for the phantom study, where the yellow box indicates the 3D imaging volume and the white lines indicate the positions of the two 1-channelRF receive coils. Green, red, and blue lines indicate the positions of the three planes where the FUS was focused, and transverse views of these planes are shown in (b-d). In the transverse views of the three focal planes, (b-d), the solid red dots indicate the approximate mechanical focal positions, and solid green dots indicate the approximate electronically steered focal positions. Red and green rings indicate heatings "mirrored" from the left hemisphere to the right hemisphere. (e) Coronalview of the experimental setup for the lamb study, where the yellow box indicates the 3D imaging volume and the white lines indicate the positions of the two 2-channelRF receive coils. Green and red lines indicate the positions of the two planes where the FUS was focused, and transverse views of these planes are shown in (f) and (g). In the transverse views in (f) and (g), black crosses correspond to the transducer's physical position, and the green dots correspond to the position of the electronically steered focal spot.

transducer was coupled to the phantom and the lamb skull with a bath of de-ionized and de-gassed water.

The phantom used in this study was composed of a single-layer polyvinyl chloride (PVC) plastic (speed of sound 2,330 m/s, specific heat capacity 1,047 J/kg*K, thermal conductivity 0.19 W/m*K [56,57]) skullcap (model A20, 3B Scientific, Tucker, GA, USA) embedded in tissue-mimicking gelatin (speed of sound 1,553 m/s, specific heat capacity 3,635 J/kg*K, thermal conductivity 0.55 W/m*K) [58]. A total of 29 sonications were performed in three planes, spaced 1.5 mm apart, utilizing both mechanical and electronic steering (Figure 1a-d). To minimize the number of sonications in this proof-of-concept study, it was assumed that the right and the left hemispheres were symmetric, and sonications were only performed in one hemisphere and then "mirrored" to the other hemisphere for treatment envelope evaluations. Each sonication was performed at 52 acoustic watts for 19.8 s, and after each heating, the phantom was allowed to cool for approximately 10 min.

For the lamb study, a recently euthanized lamb was placed in a supine position with the US focused through the skullcap (Figure 1e-g). The Institutional Animal Care and Use Committee approved all experiments. The lamb model was used since the thin skull enabled getting enough energy through the skull to cause substantial focal spot heating with our available FUS system. For this proof-of-concept study, no additional information would have been obtained if the animal had been alive, which is why the lamb was euthanized before the study. To ensure good acoustic coupling to the head and avoid trapping any air, the head was treated with hair removal cream before being placed on the FUS system. A total of 32 sonications were performed in two planes, spaced 8 mm apart. The transducer was mechanically moved to five positions in the plane more distal to the transducer and to four positions in the plane more proximal to the transducer, and from each mechanical position, the US was electronically steered ±5-mm in x and y (i.e., to the four "corners"). Each sonication was performed at 32 acoustic watts for 17.3 s, with 10 min of cooling in between sonications. No phase aberration correction of the US phases was applied in either the phantom or the lamb study.

Two different types of treatment envelopes were derived from the phantom and the lamb data. The first type utilized only the focal spot temperatures and demonstrates the type of envelope that can be derived from 2D MRTI with proper slice placement. In these envelopes, the hottest focal spot temperatures from the 29/ 32 heatings were extrapolated to the full intracranial volume. For the second type of envelope, the full 3D MRTI information was utilized, and the relative near-field heating was calculated as the ratio between the hottest focal

spot voxel and the mean of the ten hottest near-field voxels inside the skull. Ten voxels were chosen as a compromise between reducing the standard deviation without introducing an underestimation of the temperature rise. For this second type of envelope, it is desirable that the ratio be as large as possible, i.e., for every degree of temperature rise occurring in the focal spot, only a fraction of that is occurring in the near field. These ratios were also extrapolated to the full intracranial volume.

Hydrophone scans (Model HNR-0500, ONDA, Sunnyvale, CA, USA) through the skullcaps of both the plastic skull (before it was embedded in the tissue-mimicking gelatin) and the lamb skull (post-treatment and with all soft tissue removed) were performed in a bath of de-ionized and degassed water to investigate the amount of beam aberration and intensity loss that were to be expected when focusing through the skulls. The hydrophone was attached to 2D stepper motors and a 20 x 20 mm FOV with a resolution of 0.25 x 0.25 mm was scanned both with and without the skullcaps in the beam path with a power of 2.7 acoustic watts.


Figure 2 shows results from the hydrophone scans with and without the plastic and lamb skulls in the US beam path. In (a-b) and (e-f), 2D plots are shown for the plastic skull and lamb skull studies when electronically steering the US focus 5 mm in x and y. In (c-d) and (g-h), line plots through the focal spot in the x- and y-directions for the two studies are shown. All plots are normalized to the peak intensity measured in corresponding scans in the water bath only (i.e., without the skull in the US beam path). Intensity losses of about 85% (for the phantom study) and 70% (for the lamb study) were observed when focusing through the skulls. Only slight beam aberrations were observed in both cases, but in the phantom, a slight positional shift of the focus (~1 mm) in the anterior-posterior direction was observed.

In Figure 3a,b, three orthogonal views of superpositions of the hottest time frame for all 29 and 32 sonica-tions in the phantom and the lamb studies are shown. The focal spots can be seen to experience only slight beam aberration. The near field can further be seen to have more substantial temperature rise in the lamb study than in the phantom study. This is also highlighted in Figure 4, which shows transverse 2D thin-slab maximum intensity projections (MIP) of the tissue in the near field inside the skull (i.e., not including any focal spot heating, but just near-field voxels). The phantom, which can be placed closer to the transducer due to the larger acoustic window, experiences temperature rises in the near field of approximately 2°C-4°C, spread over a relatively large area (approximately 6 x 5 cm). For the lamb, the energy entering the skull is spread over a considerably smaller

Figure 2 Hydrophone scans result. 2D plot of hydrophone scans without (a) and with (b) the plastic skullin the US beam path. Line plots in the x-direction (indicated by yellow arrows) and in the y-direction (green arrows) for the two cases are shown in (c) and (d), respectively. All intensities are normalized to the maximum value in the water-only scan. The US focus is steered 5 mm in the x- and y-directions, relative to the geometric focus. Beam aberration is not significant; however, an intensity drop of about 85% and slight shift in the y-direction are seen to occur when focusing through the skull. (e-h) Corresponding data from the lamb skullcap after all soft tissue was removed. In this case, an intensity drop of about 70% and only slight beam aberration are observed.

Figure 3 Superposition of temperatures. Three orthogonalviews of superpositions of temperatures for the 29 sonications in the phantom study (a) and the 32 sonications in the lamb study (b). In both (a) and (b), the temperature maps are overlayed on magnitude images. Heatings in the right hemisphere in (a) are mirrored from the left hemisphere (see Figure 1). In both the phantom and the lamb studies, the focal spots are only slightly aberrated, in agreement with the hydrophone results in Figure 2. The near-field heating is more severe in the lamb study than in the phantom study. In (a) and (b), the three views are color coded, with the colored frame and arrows indicating each view's position in the other two views.

Figure 4 Near-field MIPs. Transverse 2D thin-slab MIPs of near-field temperatures (i.e., not including the focal spot) inside the skull for the phantom study (a-b) and the lamb study (c-d). (a) and (c) show MIPs when focusing in the plane most distal to the transducer, and (b) and (d) show MIPs when focusing in planes more proximalto the transducer. In both cases, more severe near-field heating can be seen in the lamb study than in the phantom study, in agreement with Figure 3. All MPIs are overlayed on magnitude images.

area (approximately 1 x 1.5 cm) because the transducer was located further from the skull, resulting in temperature rises of approximately 20°C-30°C.

The temporal characteristics of the sonications are demonstrated in Figure 5 as the temperature rise as a function of time for the focal spot and the near field (mean of ten hottest voxels). Data for both the phantom and the lamb studies, and for both the most distal and the most proximal planes relative to the transducer, are shown.

Figures 6 and 7 show the two types of treatment envelopes derived in this study. Figure 6 is based only on the focal spot temperatures, and Figure 7 is based on the relative near-field heating as the ratio between the focal spot and corresponding near-field heating, made possible by 3D MRTI. In Figure 6, temperature rises greater than 10°C can be seen in large parts of the brains for both the phantom and the lamb, suggesting that therapeutic levels of FUS can be delivered to these areas. Figure 7 shows that relative near-field heating ratios of up to 3 can be achieved in the phantom study, whereas for the lamb study, the majority of the brain has a ratio less than 1, indicating that the near field will experience more heating than the focal spot when attempting to treat these areas.

The reconstruction time for 250 iterations of the TCR algorithm for the phantom study was 803 s on a laptop

with 2.5 GHz Intel Core i7 processor and 8 GB RAM (MacBook Pro, Apple Inc., Cupertino, CA, USA).


The methods presented in this study are transducer design independent, and the results highlight the need for treatment envelope studies to evaluate not only feasibility and efficacy, but safety as well. With 3D MRTI coupled with constrained reconstruction for improved spatial and temporal resolution, the fully insonified FOV can simultaneously be monitored and characterized. Substantially more information about possible unintended near-field heating, for improved safety, is available compared to envelope studies utilizing 2D MRTI.

The hydrophone scans show that substantial loss in focal spot intensity occurs when focusing through both the plastic skull and the lamb skull. This is explained by the high US attenuation in the plastic and in the skull, and reflections occurring at the water/plastic or water/bone interfaces. It can be further noted that since the hydrophone scans are performed in the absence of phantom material and soft tissue, the intensity loss in actual experiments will be greater than the reported 85% and 70%, which is due only to the skulls. The focal spot shape does not experience any substantial amount of distortion in the hydrophone experiments presented, which supports not

Figure 5 Temperature rise vs. time. Temperature rise as a function of time for (a) the phantom study and (b) the lamb study. Curves for both the hottest voxelin the focalspot (solid lines) and the mean of the ten hottest near-field voxels (dashed lines) are shown, both when focusing in the plane most distalto the transducer (red lines) and when focusing in a plane more proximalto the transducer (blue lines). For the phantom, the focalspot can be seen to heat more than the near field, whereas in the lamb study, the near field experiences greater temperature rises than at the focalspot. In both cases, the near field reaches the peak temperature after the US is turned off and cools slower than the focalspot. Green arrows indicate when the US was applied.

utilizing aberration correction for the FUS heating studies. The amount of focal distortion observed in both the phantom and the lamb is limited by the relative homogeneity of the skulls and the FUS setup used. CT scans (data not shown) show that the lamb skull is considerably more homogenous than a human skull (there is very little diploe in the lamb skull), so the variation in speed of sound for the individual beam paths through the skull is very limited. This is also the case for the plastic skull, which is composed of a single homogenous layer of PVC. For the lamb study, the transducer, which has a relatively small aperture, is placed far away from the skull. Hence, the beam paths from the individual elements all go through

similar parts of the skull. The plastic skull was placed closer to the transducer than the lamb skull, but due to the limited size of the transducer, the difference in beam paths is still fairly small. From Figure 3a, it can be seen that all the focal spots in the phantom study are well defined. For the lamb study in Figure 3b, however, focal spots extending over the lateral ventricles can be seen to experience a corrupted focal spot shape due to the movement that the cerebrospinal fluid (CSF) in the ventricles experiences from the acoustic radiation force.

The sagittal and coronal views of Figure 3 also indicate that more severe near-field heating occurs in the lamb study than in the phantom study. This is further highlighted

Figure 6 Treatment envelopes based on 2D MRTI. Three orthogonalviews of the treatment envelopes based only on the focalspot heating for the phantom (a) and lamb (b) studies. Near-field temperatures are not considered for this type of envelope. Temperature rises of at least 10°C can be seen in large parts of the intracranial volume in both cases. This is the type of envelope that can be derived from 2D MRTI. In both (a) and (b), the temperature maps are overlayed on magnitude images. In (a) and (b), the three views are color coded, with the colored frame and arrows indicating each view's position in the other two views.

by the 2D MIPs of the near-field heating in Figure 4, showing that near-field heating of 2°C-4°C over a relatively large area in the phantom occurs, while in the lamb study, near-field heating of 20°C-30°C over a smaller area occurs. The main reason for the discrepancy is likely the fact that the larger size of the phantom allowed it to be positioned closer to the transducer (Figure 1a,e), and therefore, the incident US intensity was spread over a larger area at the plastic/gelatin interface than in the lamb study. The thermal properties of the plastic skull compared to the lamb skull can also have a limited effect. The specific heat of the PVC is higher than that of skull bone (1.05 kJ/kg*K vs. 0.44 kJ/kg*K), meaning that more energy is needed to heat the plastic skull; however, its thermal conductivity is lower (0.19 W/m*K vs. 0.32 W/m*K),

meaning that it retains heat longer [56,57,59]. In the lamb study, the near-field heating is more severe when focusing more proximal to the transducer (Figure 4d), while this is not the case for the phantom study (Figure 4b). It can be hypothesized that the larger electronic steering of 1.5 cm between the different planes for the phantom study (compared to 0.8 cm for the lamb) results in greater intensity losses that counteract the fact that the power is spread over a smaller area.

The plots of temperature as a function of time in Figure 5 again show that more near-field heating occurs in the lamb study than in the phantom study. It can further be seen that the peak temperature of the near field occurs after the US is turned off and that the post-heating cooling is much slower than at the focal spot. These observations

Figure 7 Treatment envelopes based on 3D MRTI. Three orthogonalviews of treatment envelopes based on the relative near-field heating, calculated as the ratio between the hottest focalspot voxeland the mean of the ten hottest near-field voxels for the 29 and 32 heatings in the phantom (a) and lamb (b) studies. This type of envelope can only be derived from 3D MRTIdata. In the phantom study, ratios of up to 3 can be observed for the most centralregions. For the lamb study, the highest ratio observed is about 1, meaning that for the majority of the intracranial volume, the near-field region willbe heated more than the focalspot. For both (a) and (b), the temperature maps are overlayed on corresponding magnitude images. In (a) and (b), the three views are color coded, with the colored frame and arrows indicating each view's position in the other two views.

are in agreement with previously published results on tissue heating close to the bone [60,61] and are explained by the high US absorption and low specific heat and thermal conductivity of the bone. High US absorption and low specific heat result in the bone heating more easily than surrounding tissue, and the low thermal conductivity results in a delayed and prolonged heat delivery to the surrounding tissue.

Treatment envelopes that are based only on the focal spot temperatures (Figure 6), which is what could be derived from 2D MRTI, suggest that with the US power applied, at least a 10°C temperature rise can be achieved in a large part of the intracranial volume. Of course, by increasing the power, even higher temperatures could have been achieved, and therapeutic levels of FUS could have been delivered to an even larger part of the intra-cranial volume. If these plots had been derived from actual 2D MRTI data, information about near-field heating would only be available at points where the 2D slices intercepted the skull. Taking into account the full 3D MRTI data by looking at the relative near-field heating in Figure 7, it can be seen that areas where Figure 6 suggests therapeutic levels of FUS can be delivered can in fact not be safely treated. This is especially true for the lamb study, where the ratio is less than 1 for the majority of the intracranial volume. This indicates that the temperature rise at the focal spot will be less than that at the skull. It should be noted that neither blood perfusion nor CSF circulation was present in the phantom or in the lamb study, and both would act as heat sinks both at the focal spot and in the near field in clinical studies. The presented treatment envelopes emphasize the need for envelope studies to incorporate not only feasibility and efficacy, but safety as well. It is shown that safety can be incorporated by monitoring the fully insonified 3D FOV and that this leads to a very different envelope than if only feasibility and efficacy are considered.

The amount of near-field heating observed in this proof-of-concept study is almost certainly more severe than what will be observed with the clinical tcMRgFUS systems used in previous treatment envelope studies [51,52] and current clinical trials [24]. These clinical systems operate at lower frequencies of 0.2 and 0.7 MHz that are more optimal for transcranial transmission and brain absorption [62,63] than the 1-MHz transducer used in this study. Compared to the transducer used in this study, the elements of the clinical transducers are also spread over a larger hemispherical surface that surrounds the head. This spreads the US energy over a larger surface area of the skull and thereby minimizes the near-field heating.

As in all iterative reconstruction methods, the reconstruction time of the TCR algorithm is fairly long. A real-time version of the TCR algorithm [64] has been

presented and significantly decreases the reconstruction time to allow real-time availability of the reconstructed temperature maps. The reconstruction time is decreased by only applying the TCR algorithm over the focal spot region where the temperatures change most rapidly, and performs a significantly less computationally heavy sliding window reconstruction over areas where the temperature changes more gradually. The algorithm was implemented on a GPU to further speed up the reconstruction. However, for treatment envelope studies, realtime availability of the temperature maps is generally of secondary interest.


This study demonstrates the concept and importance of assessing treatment envelopes in terms of both safety and efficacy. The treatment envelopes derived in this study incorporate safety assessments in terms of near-field heating evaluations and show that due to unintended near-field heating, it is not necessarily safe to treat in all areas where therapeutic levels of FUS can be delivered. The methods described in this study can be utilized to maximize the treatment envelope when designing new transducers for tcMRgFUS, as well as to evaluate existing systems. Since this proof-of-concept study was performed utilizing a non-optimal transducer design, and only in a phantom and a lamb model, more studies with clinical tcMRgFUS systems and human cadaver skulls will be necessary to derive meaningful treatment envelopes for these systems. The results in this study further indicate that from a safety point of view, tcMRgFUS treatments are ideally monitored with 3D MRTI covering the fully insonified field of view.

Competing interests

The authors declare that they have no competing interests. Authors' contributions

HO and JdB performed the MRgFUS experiments. HO conducted the data reconstruction and analysis, and drafted the manuscript. SA and DAC were responsible for the hydrophone equipment, and SA and HO performed the hydrophone experiments. HO, NT, AP, JWS, and DLP designed the study. HO, NT, AP, and DLP developed and implemented the software algorithms used. AF manufactured the skullphantom. Allauthors read and approved the final manuscript.


The authors thank MarJanna Dahland Dr. Kurt H Albertine for help with the lamb study. This work was supported by The Focused Ultrasound Surgery Foundation, Siemens Healthcare, The Ben B. and Iris M. Margolis Foundation, and NIH grants F32 EB012917-02 and R01s EB013433, CA134599, and HL110002.

Author details

1Utah Center for Advanced Imaging Research, Department of Radiology, University of Utah, Salt Lake City, Utah 84108, USA. 2Department of Physics and Astronomy, University of Utah, Salt Lake City, Utah 84112, USA. 3School of Computing, University of Utah, Salt Lake City, Utah 84112, USA. 4Department of Bioengineering, University of Utah, Salt Lake City, Utah 84112, USA. 5Focused Ultrasound Foundation, Charlottesville, Virginia 22903, USA. 6Department of NeurologicalSurgery, University of Virginia,

Charlottesville, Virginia 22908, USA. 'Department of Electrical and Computer Engineering, University of Utah, Salt Lake City, Utah 84112, USA.

Received: 24 June 2014 Accepted: 17 September 2014 Published: 16 October 2014


1. Wu F, Wang Z-B, Chen W-Z, Zou J-Z, Bai J, Zhu H, Li K-Q, Xie F-L, Jin C-B, Su H-B, Gao G-W. Extracorporeal focused ultrasound surgery for treatment of human solid carcinomas: early Chinese clinical experience. Ultrasound Med Biol. 2004; 30:245-60.

2. Li S, Wu P-H. Magnetic resonance image-guided versus ultrasound-guided high-intensity focused ultrasound in the treatment of breast cancer. Chin J Cancer. 2013; 32:441 -52.

3. Gelet A, Chapelon J, Margonari J, Theillere Y, Gorry F, Souchon R, Bouvier R. High-intensity focused ultrasound experimentation on human benign prostatic hypertrophy. Eur Urol. 1993; 23:44-7.

4. Yang R, Sanghvi N, Rescorla F, Kopecky K, Grosfeld J. Liver cancer ablation with extracorporeal high-intensity focused ultrasound. Eur Urol. 1993; 23:17-22.

5. Medel R, Monteith SJ, Elias WJ, Eames M, Snell J, Sheehan JP, Wintermark M, Jolesz FA, Kassell NF. Magnetic resonance-guided focused ultrasound surgery: part 2: a review of current and future applications.

Neurosurgery. 2012; 71:755-63.

6. Schlesinger D, Benedict S, Diederich C, Gedroyc W, Klibanov A, Larner J. MR-guided focused ultrasound surgery, present and future. Med Phys. 2013; 40:080901.

7. Ellis S, Rieke V, Kohi M, Westphalen AC. Clinical applications for magnetic resonance guided high intensity focused ultrasound (MRgHIFU): present and future. J Med Imaging Radiat Oncol. 2013; 57:391 -9.

8. Wahab RA, Choi M, Liu Y, Krauthamer V, Zderic V, Myers MR. Mechanical bioeffects of pulsed high intensity focused ultrasound on a simple neural model. Med Phys. 2012; 39:4274-83.

9. Hempel CR, Hall TL, Cain CA, Fowlkes JB, Xu Z, Roberts WW. Histotripsy fractionation of prostate tissue: local effects and systemic response in a canine model. J Urol. 2011; 185:1484-9.

10. Wright C, Hynynen K, Goertz D. In vitro and in vivo high intensity focused ultrasound thrombolysis. Invest Radiol. 2012; 47:217-25.

11. McDannold N, Zhang Y-Z, Power C, Jolesz F, Vykhodtseva N. Nonthermal ablation with microbubble-enhanced focused ultrasound close to the optic tract without affecting nerve function. J Neurosurg. 2013; 119:1208-20.

12. Chopra R, Burtnyk M, N'djin WA, Bronskill M. MRI-controlled transurethral ultrasound therapy for localised prostate cancer. Int J Hyperth. 2010; 26:804-21.

13. Crouzet S, Murat FJ, Pasticier G, Cassier P, Chapelon JY, Gelet A. High intensity focused ultrasound (HIFU) for prostate cancer: current clinical status, outcomes and future perspectives. Int J Hyperth. 2010; 26:796-803.

14. Salgaonkar VA, Prakash P, Rieke V, Ozhinsky E, Plata J, Kurhanewicz J, Hsu I-CJ, Diederich CJ. Model-based feasibility assessment and evaluation of prostate hyperthermia with a commercial MR-guided endorectal HIFU ablation array. Med Phys. 2014; 41:033301.

15. Wijlemans J, Bartels L, Deckers R, Ries M, Mali W, Moonen C, van den Bosch M. Magnetic resonance-guided high-intensity focused ultrasound (MR-HIFU) ablation of liver tumours. Cancer Imaging. 2013; 12:387-94.

16. Illing RO, Kennedy JE, Wu F, ter Haar GR, Protheroe AS, Friend PJ, Gleeson FV, Cranston DW, Phillips RR, Middleton MR. The safety and feasibility of extracorporeal high-intensity focused ultrasound (HIFU) for the treatment of liver and kidney tumours in a Western population. Br J Cancer. 2005; 93:890-5.

17. Anzidei M, Napoli A, Sandolo F, Marincola B, Di Martino M, Berloco P, Bosco S, Bezzi M, Catalano C. Magnetic resonance-guided focused ultrasound ablation in abdominal moving organs: a feasibility study in selected cases of pancreatic and liver cancer Cardiovasc Intervent Radiol. 2014. doi:10.1007/s00270-014-0861 -x.

18. Schmitz AC, Gianfelice D, Daniel BL, Mali WPTM, van den Bosch MAAJ. Image-guided focused ultrasound ablation of breast cancer: current status, challenges, and future directions. Eur Radiol. 2008; 18:1431-41.

19. Payne A, Merrill R, Minalga E, Vyas U, de Bever J, Todd N, Hadley R, Dumont E, Neumayer L, Christensen D, Roemer R, Parker D. Design and characterization of a laterally mounted phased-array transducer breast-specific MRgHIFU device with integrated 11-channel receiver array. Med Phys. 2012; 39:1552-60.

20. Wu F, Wang Z, Zhu H, Chen W, Zou J, Bai J, Li K, Jin C, Xie F, Su H. Extracorporeal high intensity focused ultrasound treatment for patients with breast cancer. Breast Cancer Res Treat. 2005; 92:51-60.

21. Shen S-H, Fennessy F, McDannold N, Jolesz F, Tempany C. Image-guided thermal therapy of uterine fibroids. Semin Ultyrasound CT MR. 2009; 30:91-104.

22. Hesley GK, Gorny KR, Woodrum DA. MR-guided focused ultrasound for the treatment of uterine fibroids. Cardiovasc Intervent Radiol. 2013; 36:5-13.

23. Martin E, Jeanmonod D, Morel A, Zadicario E, Werner B. High-intensity focused ultrasound for noninvasive functional neurosurgery. Ann Neurol. 2009; 66:858-61.

24. Elias WJ, Huss D, Voss T, Loomba J, Khaled M, Zadicario E, Frysinger RC, Sperling SA, Wylie S, Monteith SJ, Druzgal J, Shah BB, Harrison M, Wintermark M. A pilot study of focused ultrasound thalamotomy for essential tremor. N Engl J Med. 2013; 369:640-8.

25. McDannold N, Clement G, Black P, Jolesz F, Hynynen K. Transcranial MRI-guided focused ultrasound surgery of brain tumors: initial findings in three patients. Neurosurgery. 2010; 66:323-32.

26. Lynn JG, Putnam TJ. Histology of cerebral lesions produced by focused ultrasound. Am J Pathol. 1944; 20:637-49.

27. Fry F. Transkull transmission of an intense focused ultrasonic beam. Ultrasound Med Biol. 1977; 3:179-84.

28. Lindstrom PA. Prefrontal ultrasonic irradiation - a substitute for lobotomy. AMA Arch Neurol Psychiatry. 1954; 72:399-425.

29. Fry WJ, Mosberg WH, Barnard JW, Fry FJ. Production of focal destructive lesions in the central nervous system with ultrasound. J Neurosurg. 1954; 11:471-8.

30. Monteith S, Sheehan J, Medel R, Wintermark M, Eames M, Snell J, Kassell NF, Elias WJ. Potential intracranial applications of magnetic resonance-guided focused ultrasound surgery. J Neurosurg. 2013; 118:215-21.

31. Fry JF, Ades HW, Fry WJ. Production of reversible changes in the central nervous system by ultrasound. Science (80-). 1958; 127:83-4.

32. Legon W, Sato TF, Opitz A, Mueller J, Barbour A, Williams A, Tyler WJ. Transcranial focused ultrasound modulates the activity of primary somatosensory cortex in humans. Nat Neurosci. 2014; 17:322-9.

33. Legon W, Rowlands A, Opitz A, Sato TF, Tyler WJ. Pulsed ultrasound differentially stimulates somatosensory circuits in humans as indicated by EEG and FMRI. PLoS One. 2012; 7:e51177.

34. Bystritsky A, Korb AS, Douglas PK, Cohen MS, Melega WP, Mulgaonkar AP, DeSalles A, Min B-K, Yoo S-S. A review of low-intensity focused ultrasound pulsation. Brain Stimul. 2011; 4:125-36.

35. Hynynen K, McDannold N, Vykhodtseva N, Jolesz FA. Non-invasive opening of BBB by focused ultrasound. Acta Neurochir Suppl. 2003; 86:555-8.

36. Fry FJ, Barger JE. Acoustical properties of the human skull. J Acoust Soc Am. 1978; 63:1576-90.

37. Clement GT, Hynynen K. A non-invasive method for focusing ultrasound through the human skull. Phys Med Biol. 2002; 47:1219-36.

38. Aubry J-F, Tanter M, Pernot M, Thomas J-L, Fink M. Experimental demonstration of noninvasive transskull adaptive focusing based on prior computed tomography scans. J Acoust Soc Am. 2003; 113:84.

39. Almquist S, Todd N, Parker DL, Christensen D. Correcting phase aberrations at multiple treatment sites in transcranial brain HIFU. In: Proc. 3rd Int. Symp. Focus. Bethesda, Md: Ultrasound Surg. Found; 2012: p. 148. P-91-BR.

40. Marsac L, Chauvet D, Larrat B, Pernot M, Robert B, Fink M, Boch AL, Aubry JF, Tanter M. MR-guided adaptive focusing of therapeutic ultrasound beams in the human head. Med Phys. 2012; 39:1141-9.

41. Kaye EA, Hertzberg Y, Marx M, Werner B, Navon G, Levoy M, Pauly KB. Application of Zernike polynomials towards accelerated adaptive focusing of transcranial high intensity focused ultrasound. Med Phys. 2012; 39:6254-63.

42. Larrat B, Pernot M, Montaldo G, Fink M, Tanter M. MR-guided adaptive focusing of ultrasound. IEEE Trans Ultrason Ferroelectr Freq Control. 2010; 57:1734-7.

43. Kyriakou A, Neufeld E, Werner B, Paulides M, Szekely G, Kuster N. A review of numerical and experimental compensation techniques for skull-induced phase aberrations in transcranial focused ultrasound. Int J Hyperth. 2014; 30:36-46.

44. Diakite M, Payne A, Todd N, Parker DL. Irreversible change in the T1 temperature dependence with thermal dose using the proton resonance frequency-T1 technique. Magn Reson Med. 2013; 69:1122-30.

45. Hectors S, Jacobs I, Strijkers G, Nicolay K. Identification of high intensity focused ultrasound treated tumor tissue using a multiparametric MRI

protocol and ISODATA analysis. Proc 20th ISMRM Annu Meet Melbourne, Aust. 2012; 1537.

46. Hectors SJCG, Moonen RPM, Strijkers GJ, Nicolay K. T1 p mapping for the evaluation of high intensity focused ultrasound tumor treatment. Magn Reson Med. 2014; 00:1-9.

47. Parker DL. Applications of NMR imaging in hyperthermia: an evaluation of the potential for localized tissue heating and noninvasive temperature monitoring. IEEE Trans Biomed Eng. 1984; 31:161-7.

48. Graham SJ, Bronskill MJ, Henkelman RM. Time and temperature dependence of MR parameters during thermal coagulation of ex vivo rabbit muscle. Magn Reson Med. 1998; 39:198-203.

49. Ishihara Y, Calderon A, Watanabe H, Okamoto K, Suzuki Y, Kuroda K. A precise and fast temperature mapping using water proton chemical shift. Magn Reson Med. 1995; 34:814-23.

50. De Poorter J, De Wagter C, De Deene Y, Thomsen C, Stahlberg F, Achten E. Noninvasive MRI thermometry with the proton resonance frequency (PRF) method: in vivo results in human muscle. Magn Reson Med. 1995; 33:74-81.

51. Eames M, Hananel A, Snell J. Transcranial MRg-FUS treatment envelope.

In: Int. Symp. Focus. Washington DC: Ultrasound; 2012: p. 124.

52. Dallapiazza RF, Eames M, Snell J, Khaled M, Elias WJ. A cadaveric study to define the treatment envelope for transcranial MR-guided focused ultrasound ablation. In: Int Symp Ther Ultrasound Las Vegas. 2014.

53. Adluru G, Awate SP, Tasdizen T, Whitaker RT, Dibella EVR. Temporally constrained reconstruction of dynamic cardiac perfusion MRI. Magn Reson Med. 2007; 57:1027-36.

54. Todd N, Adluru G, Payne A, DiBella EVR, Parker D. Temporally constrained reconstruction applied to MRI temperature data. Magn Reson Med. 2009; 62:406-19.

55. Hayes CE, Mathis CM, Yuan C. Surface coil phased arrays for highresolution imaging of the carotid arteries. J Magn Reson Imaging. 1996; 6:109-12.

56. Mark JE. Physical Properties of Polymers Handbook. New York: Springer Science & Business Media; 2007: p. 1096.

57. Engineering ToolBox. []. Accessed on 22 Aug 2014.

58. Farrer AI, de Bever J, Coats B, Christenson DA, Payne A. Fabrication and evaluation of tissue-mimicking phantoms for use with MR-ARFI and MRgFUS. In: Int. Symp. Ther. Las Vegas: Ultrasound; 2014: p. 232.

59. Tissue properties database - IT'IS Foundation. []. Accessed on 22 Aug 2014.

60. McDannold N, King RL, Hynynen K. MRI monitoring of heating produced by ultrasound absorption in the skull: in vivo study in pigs. Magn Reson Med. 2004; 51:1061 -5.

61. Webb TD, Bitton R, Ghanouni P. Butts Pauly K. Spatial and temporal characteristics of soft tissue heating in MR-HIFU treatment of bone metastasis In Int Soc Magn Reson Med. 2014; 22:2344.

62. White PJ, Clement GT, Hynynen K. Local frequency dependence in transcranial ultrasound transmission. Phys Med Biol. 2006; 51:2293-305.

63. Hayner M, Hynynen K. Numerical analysis of ultrasonic transmission and absorption of oblique plane waves through the human skull. J Acoust Soc Am. 2001; 110:3319.

64. Todd N, Prakash J, Odeen H, de Bever J, Payne A, Yalavarthy P, Parker DL. Toward real-time availability of 3D temperature maps created with temporally constrained reconstruction. Magn Reson Med. 2014; 71:1394-404.


Cite this article as: Odéen et al.: Treatment envelope evaluation in transcranial magnetic resonance-guided focused ultrasound utilizing 3D MR thermometry. Journal of Therapeutic Ultrasound 2014 2:19.

Submit your next manuscript to BioMed Central and take full advantage of:

• Convenient online submission

• Thorough peer review

• No space constraints or color figure charges

• Immediate publication on acceptance

• Inclusion in PubMed, CAS, Scopus and Google Scholar

• Research which is freely available for redistribution

Submit your manuscript at Tpntral >_ciiuai