Characterization of central manifestations in patients with Niemann-Pick disease type C

PURPOSE

Niemann-Pick disease type C (NPC) is a rare lysosomal storage disease characterized by progressive neurodegeneration and neuropsychiatric symptoms. This study investigated pathophysiological mechanisms underlying motor deficits, particularly speech production, and cognitive impairment.

METHODS

We prospectively phenotyped 8 adults with NPC and age-sex-matched healthy controls using a comprehensive assessment battery, encompassing clinical presentation, plasma biomarkers, hand-motor skills, speech production, cognitive tasks, and (micro-)structural and functional central nervous system properties through magnetic resonance imaging.

RESULTS

Patients with NPC demonstrated deficits in fine-motor skills, speech production timing and coordination, and cognitive performance. Magnetic resonance imaging revealed reduced cortical thickness and volume in cerebellar subdivisions (lobule VI and crus I), cortical (frontal, temporal, and cingulate gyri) and subcortical (thalamus and basal ganglia) regions, and increased choroid plexus volumes in NPC. White matter fractional anisotropy was reduced in specific pathways (intracerebellar input and Purkinje tracts), whereas diffusion tensor imaging graph theory analysis identified altered structural connectivity. Patients with NPC exhibited altered activity in sensorimotor and cognitive processing hubs during resting-state and speech production. Canonical component analysis highlighted the role of cerebellar-cerebral circuitry in NPC and its integration with behavioral performance and disease severity.

CONCLUSION

This deep phenotyping approach offers a comprehensive systems neuroscience understanding of NPC motor and cognitive impairments, identifying potential central nervous system biomarkers.

Focused ultrasound-mediated delivery of viral neuronal tracers in marmosets

In this issue of Cell Reports Methods, Parks et al. present an approach to non-invasively deliver adeno-associated viruses in a marmoset model using focused ultrasound (FUS) for neuronal tracing. The optimization of this technique in this non-human primate model is highly valuable for future FUS-mediated drug delivery studies.

A consensus protocol for functional connectivity analysis in the rat brain

Task-free functional connectivity in animal models provides an experimental framework to examine connectivity phenomena under controlled conditions and allows for comparisons with data modalities collected under invasive or terminal procedures. Currently, animal acquisitions are performed with varying protocols and analyses that hamper result comparison and integration. Here we introduce StandardRat, a consensus rat functional magnetic resonance imaging acquisition protocol tested across 20 centers. To develop this protocol with optimized acquisition and processing parameters, we initially aggregated 65 functional imaging datasets acquired from rats across 46 centers. We developed a reproducible pipeline for analyzing rat data acquired with diverse protocols and determined experimental and processing parameters associated with the robust detection of functional connectivity across centers. We show that the standardized protocol enhances biologically plausible functional connectivity patterns relative to previous acquisitions. The protocol and processing pipeline described here is openly shared with the neuroimaging community to promote interoperability and cooperation toward tackling the most important challenges in neuroscience.

Targeted manipulation of pain neural networks: The potential of focused ultrasound for treatment of chronic pain

Focused ultrasound (FUS) is a promising technology for facilitating treatment of brain diseases including chronic pain. Focused ultrasound is a unique modality for delivering therapeutic levels of energy into the body, including the central nervous system (CNS). It is non-invasive and can target spatially localized effects through the intact skull to cortical or subcortical regions of the brain. FUS can achieve three different mechanisms of action in the brain that are relevant for chronic pain treatment: (1) localized thermal ablation of neural tissue; (2) localized and transient disruption of the blood-brain barrier for targeted drug delivery to CNS structures; and (3) inhibition or stimulation of neuronal activity in targeted regions. This review provides an in-depth look at the technology of FUS with emphasis placed on applications to CNS-based treatments of chronic pain. While still in the early stages of clinical translation and with some technical challenges remaining, we suggest that FUS has great potential as a novel approach for manipulating CNS networks involved in pain treatment.

Secondary effects on brain physiology caused by focused ultrasound-mediated disruption of the blood-brain barrier

Focused ultrasound (FUS) combined with microbubbles is a non-invasive method for targeted, reversible disruption of the blood-brain barrier (FUS-BBB opening). This approach holds great promise for improving delivery of therapeutics to the brain. In order to achieve this clinically important goal, the approach necessarily breaks a protective barrier, temporarily, which plays a fundamental role in maintaining a homeostatic environment in the brain. Preclinical and clinical research has identified a set of treatment parameters under which this can be performed safely, whereby the BBB is disrupted to the point of being permeable to normally non-penetrant agents without causing significant acute damage to endothelial or neuronal cells. Much of the early work in this field focused on engineering questions around how to achieve optimal delivery of therapeutics via BBB disruption. However, there is increasing interest in addressing biological questions related to whether and how various aspects of neurophysiology might be affected when this fundamental protective barrier is compromised by the specific mechanisms of FUS-BBB opening. Improving our understanding of these secondary effects is becoming vital now that FUS-BBB opening treatments have entered clinical trials. Such information would help to safely expand FUS-BBB opening protocols into a wider range of drug delivery applications and may even lead to new types of treatments. In this paper, we will critically review our current knowledge of the secondary effects caused by FUS-BBB opening on brain physiology, identify areas that remain understudied, and discuss how a better understanding of these processes can be used to safely advance FUS-BBB opening into a wider range of clinical applications.

The Use of Portable Computers: A Survey of Visually Impaired College Students

This survey of 10 visually impaired college students found that the students used portable (laptop) computers to produce assignments and to perform other tasks, as well as to take notes in class. The use of laptops seemed to foster the students’ independence and gave them sufficient time to complete their work.

The neurovascular response is attenuated by focused ultrasound-mediated disruption of the blood-brain barrier

Focused ultrasound (FUS)-induced disruption of the blood-brain barrier (BBB) is a non-invasive method to target drug delivery to specific brain areas that is now entering into the clinic. Recent studies have shown that the method has several secondary effects on local physiology and brain function beyond making the vasculature permeable to normally non-BBB penetrant molecules. This study uses functional MRI methods to investigate how FUS BBB opening alters the neurovascular response in the rat brain. Nine rats underwent actual and sham FUS induced BBB opening targeted to the right somatosensory cortex (SI) followed by four runs of bilateral electrical hind paw stimulus-evoked fMRI. The neurovascular response was quantified using measurements of the blood oxygen level dependent (BOLD) signal and cerebral blood flow (CBF). An additional three rats underwent the same FUS-BBB opening followed by stimulus-evoked fMRI with high resolution BOLD imaging and BOLD imaging of a carbogen-breathing gas challenge. BOLD and CBF measurements at two different stimulus durations demonstrate that the neurovascular response to the stimulus is attenuated in both amplitude and duration in the region targeted for FUS-BBB opening. The carbogen results show that the attenuation in response amplitude, but not duration, is still present when the signaling mechanism originates from changes in blood oxygenation instead of stimulus-induced neuronal activity. There is some evidence of non-local effects, including a possible global decrease in baseline CBF. All effects are resolved by 24 h after FUS-BBB opening. Taken together, these results suggest that FUS-BBB opening alters that state of local brain neurovascular physiology in such a way that hinders its ability to respond to demands for increased blood flow to the region. The mechanisms for this effect need to be elucidated.

Modulation of brain function by targeted delivery of GABA through the disrupted blood-brain barrier

The technology of transcranial focused ultrasound (FUS) enables a novel approach to neuromodulation, a tool for selective manipulation of brain function to be used in neurobiology research and with potential applications in clinical treatment. The method uses transcranial focused ultrasound to non-invasively open the blood-brain barrier (BBB) in a localized region such that a systemically injected neurotransmitter chemical can be delivered to the targeted brain site. The approach modulates the chemical signaling that occurs in and between neurons, making it complimentary to most other neuromodulation techniques that affect the electrical properties of neuronal activity. Here, we report delivering the inhibitory neurotransmitter GABA to the right somatosensory cortex of the rat brain during bilateral hind paw electrical stimulation and measure the inhibition of activation using functional MRI (fMRI). In a 2 × 2 factorial design, we evaluated conditions of BBB Closed vs BBB Open and No GABA vs GABA. Results from fMRI measurements of the blood oxygen level-dependent (BOLD) signal show: 1) intravenous GABA injection without FUS-mediated BBB opening does not have an effect on the BOLD response; 2) FUS-mediated BBB opening alone significantly alters the BOLD signal response to the stimulus, both in amplitude and shape of the time course; 3) the combination of FUS-mediated BBB opening and GABA injection further reduces the peak amplitude and spatial extent of the BOLD signal response to the stimulus. The data support the thesis that FUS-mediated opening of the BBB can be used to achieve non-invasive delivery of neuroactive substances for targeted manipulation of brain function.

Optimizing Data for Modeling Neuronal Responses

In this technical note, we address an unresolved challenge in neuroimaging statistics: how to determine which of several datasets is the best for inferring neuronal responses. Comparisons of this kind are important for experimenters when choosing an imaging protocol-and for developers of new acquisition methods. However, the hypothesis that one dataset is better than another cannot be tested using conventional statistics (based on likelihood ratios), as these require the data to be the same under each hypothesis. Here we present Bayesian data comparison (BDC), a principled framework for evaluating the quality of functional imaging data, in terms of the precision with which neuronal connectivity parameters can be estimated and competing models can be disambiguated. For each of several candidate datasets, neuronal responses are modeled using Bayesian (probabilistic) forward models, such as General Linear Models (GLMs) or Dynamic Casual Models (DCMs). Next, the parameters from subject-specific models are summarized at the group level using a Bayesian GLM. A series of measures, which we introduce here, are then used to evaluate each dataset in terms of the precision of (group-level) parameter estimates and the ability of the data to distinguish similar models. To exemplify the approach, we compared four datasets that were acquired in a study evaluating multiband fMRI acquisition schemes, and we used simulations to establish the face validity of the comparison measures. To enable people to reproduce these analyses using their own data and experimental paradigms, we provide general-purpose Matlab code via the SPM software.

Accurate modeling of temporal correlations in rapidly sampled fMRI time series

Rapid imaging techniques are increasingly used in functional MRI studies because they allow a greater number of samples to be acquired per unit time, thereby increasing statistical power. However, temporal correlations limit the increase in functional sensitivity and must be accurately accounted for to control the false‐positive rate. A common approach to accounting for temporal correlations is to whiten the data prior to estimating fMRI model parameters. Models of white noise plus a first‐order autoregressive process have proven sufficient for conventional imaging studies, but more elaborate models are required for rapidly sampled data. Here we show that when the “FAST” model implemented in SPM is used with a well‐controlled number of parameters, it can successfully prewhiten 80% of grey matter voxels even with volume repetition times as short as 0.35 s. We further show that the temporal signal‐to‐noise ratio (tSNR), which has conventionally been used to assess the relative functional sensitivity of competing imaging approaches, can be augmented to account for the temporal correlations in the time series. This amounts to computing the t‐score testing for the mean signal. We show in a visual perception task that unlike the tSNR weighted by the number of samples, the t‐score measure is directly related to the t‐score testing for activation when the temporal correlations are correctly modeled. This score affords a more accurate means of evaluating the functional sensitivity of different data acquisition options.

Focused ultrasound induced opening of the blood-brain barrier disrupts inter-hemispheric resting state functional connectivity in the rat brain

Focused ultrasound (FUS) is a technology capable of delivering therapeutic levels of energy through the intact skull to a tightly localized brain region. Combining the FUS pressure wave with intravenously injected microbubbles creates forces on blood vessel walls that open the blood-brain barrier (BBB). This noninvasive and localized opening of the BBB allows for targeted delivery of pharmacological agents into the brain for use in therapeutic development. It is possible to use FUS power levels such that the BBB is opened without damaging local tissues. However, open questions remain related to the effects that FUS-induced BBB opening has on brain function including local physiology and vascular hemodynamics. We evaluated the effects that FUS-induced BBB opening has on resting state functional magnetic resonance imaging (rs-fMRI) metrics. Data from rs-fMRI was acquired in rats that underwent sham FUS BBB vs. FUS BBB opening targeted to the right primary somatosensory cortex hindlimb region (S1HL). FUS BBB opening reduced the functional connectivity between the right S1HL and other sensorimotor regions, including statistically significant reduction of connectivity to the homologous region in the left hemisphere (left S1HL). The effect was observed in all three metrics analyzed: functional connectivity between anatomically defined regions, whole brain voxel-wise correlation maps based on anatomical seeds, and spatial patterns from independent component analysis. Connectivity metrics for other regions where the BBB was not perturbed were not affected. While it is not clear whether the effect is vascular or neuronal in origin, these results suggest that even safe levels of FUS BBB opening have an effect on the physiological processes that drive the signals measured by BOLD fMRI. As such these effects must be accounted for when carrying out studies using fMRI to evaluate the effects of pharmacological agents delivered via FUS-induced BBB opening.

Functional Sensitivity of 2D Simultaneous Multi-Slice Echo-Planar Imaging: Effects of Acceleration on g-factor and Physiological Noise

Accelerated data acquisition with simultaneous multi-slice (SMS) imaging for functional MRI studies leads to interacting and opposing effects that influence the sensitivity to blood oxygen level-dependent (BOLD) signal changes. Image signal to noise ratio (SNR) is decreased with higher SMS acceleration factors and shorter repetition times (TR) due to g-factor noise penalties and saturation of longitudinal magnetization. However, the lower image SNR is counteracted by greater statistical power from more samples per unit time and a higher temporal Nyquist frequency that allows for better removal of spurious non-BOLD high frequency signal content. This study investigated the dependence of the BOLD sensitivity on these main driving factors and their interaction, and provides a framework for evaluating optimal acceleration of SMS-EPI sequences. functional magnetic resonance imaging (fMRI) data from a scenes/objects visualization task was acquired in 10 healthy volunteers at a standard neuroscience resolution of 3 mm on a 3T MRI scanner. SMS factors 1, 2, 4, and 8 were used, spanning TRs of 2800 ms to 350 ms. Two data processing methods were used to equalize the number of samples over the SMS factors. BOLD sensitivity was assessed using g-factors maps, temporal SNR (tSNR), and t-score metrics. tSNR results show a dependence on SMS factor that is highly non-uniform over the brain, with outcomes driven by g-factor noise amplification and the presence of high frequency noise. The t-score metrics also show a high degree of spatial dependence: the lower g-factor noise area of V1 shows significant improvements at higher SMS factors; the moderate-level g-factor noise area of the parahippocampal place area shows only a trend of improvement; and the high g-factor noise area of the ventral-medial pre-frontal cortex shows a trend of declining t-scores at higher SMS factors. This spatial variability suggests that the optimal SMS factor for fMRI studies is region dependent. For task fMRI studies done with similar parameters as were used here (3T scanner, 32-channel RF head coil, whole brain coverage at 3 mm isotropic resolution), we recommend SMS accelerations of 4x (conservative) to 8x (aggressive) for most studies and a more conservative acceleration of 2x for studies interested in anterior midline regions.

Preclinical evaluation of a low-frequency transcranial MRI-guided focused ultrasound system in a primate model

This study investigated thermal ablation and skull-induced heating with a 230 kHz transcranial MRI-guided focused ultrasound (TcMRgFUS) system in nonhuman primates. We evaluated real-time acoustic feedback and aimed to understand whether cavitation contributed to the heating and the lesion formation. In four macaques, we sonicated thalamic targets at acoustic powers of 34-560 W (896-7590 J). Tissue effects evaluated with MRI and histology were compared to MRI-based temperature and thermal dose measurements, acoustic emissions recorded during the experiments, and acoustic and thermal simulations. Peak temperatures ranged from 46 to 57 °C, and lesions were produced in 5/8 sonicated targets. A linear relationship was observed between the applied acoustic energy and both the focal and brain surface heating. Thermal dose thresholds were 15-50 cumulative equivalent minutes at 43 °C, similar to prior studies at higher frequencies. Histology was also consistent with earlier studies of thermal effects in the brain. The system successfully controlled the power level and maintained a low level of cavitation activity. Increased acoustic emissions observed in 3/4 animals occurred when the focal temperature rise exceeded approximately 16 °C. Thresholds for thermally-significant subharmonic and wideband emissions were 129 and 140 W, respectively, corresponding to estimated pressure amplitudes of 2.1 and 2.2 MPa. Simulated focal heating was consistent with the measurements for sonications without thermally-significant acoustic emissions; otherwise it was consistently lower than the measurements. Overall, these results suggest that the lesions were produced by thermal mechanisms. The detected acoustic emissions, however, and their association with heating suggest that cavitation might have contributed to the focal heating. Compared to earlier work with a 670 kHz TcMRgFUS system, the brain surface heating was substantially reduced and the focal heating was higher with this 230 kHz system, suggesting that a reduced frequency can increase the treatment envelope for TcMRgFUS and potentially reduce the risk of skull heating.

Prior image based temporally constrained reconstruction algorithm for magnetic resonance guided high intensity focused ultrasound

PURPOSE

A prior image based temporally constrained reconstruction (PITCR) algorithm was developed for obtaining accurate temperature maps having better volume coverage, and spatial, and temporal resolution than other algorithms for highly undersampled data in magnetic resonance (MR) thermometry.

METHODS

The proposed PITCR approach is an algorithm that gives weight to the prior image and performs accurate reconstruction in a dynamic imaging environment. The PITCR method is compared with the temporally constrained reconstruction (TCR) algorithm using pork muscle data.

RESULTS

The PITCR method provides superior performance compared to the TCR approach with highly undersampled data. The proposed approach is computationally expensive compared to the TCR approach, but this could be overcome by the advantage of reconstructing with fewer measurements. In the case of reconstruction of temperature maps from 16% of fully sampled data, the PITCR approach was 1.57× slower compared to the TCR approach, while the root mean square error using PITCR is 0.784 compared to 2.815 with the TCR scheme.

CONCLUSIONS

The PITCR approach is able to perform more accurate reconstructions of temperature maps compared to the TCR approach with highly undersampled data in MR guided high intensity focused ultrasound.

Evaluation of a three-dimensional MR acoustic radiation force imaging pulse sequence using a novel unbalanced bipolar motion encoding gradient

PURPOSE

MR guided focused ultrasound procedures require accurate focal spot localization in three dimensions. This study presents a three-dimensional (3D) pulse sequence for acoustic radiation force imaging (ARFI) that efficiently localizes the focal spot by means of ultrasound induced tissue displacement over a large field-of-view.

METHODS

A novel unbalanced bipolar motion encoding gradient was implemented to maximize time available for motion encoding, reduce echo times, and allow for longer echo train lengths. Two advanced features, kz reduction factor (KZRF) and kz -level interleaving, were implemented to reduce tissue heating. Studies in gelatin phantoms compared the location of peak displacement and temperature measured by 3D MR thermometry. MR-ARFI induced tissue heating was evaluated through a parametric study of sequence parameters and MR thermometry measurements during repeated application of ARFI sonication patterns. Sequence performance was characterized in the presence of respiration and tissue inhomogeneity.

RESULTS

The location of peak displacement and temperature rise agreed within 0.2 ± 0.1 mm and 0.5 ± 0.3 mm in the transverse and longitudinal direction, respectively. The 3D displacement maps were acquired safely, and the KZRF and kz -level interleaving features reduced tissue heating by 51%. High quality displacement maps were obtained despite respiration and tissue inhomogeneities.

CONCLUSION

This sequence provides a safe, accurate, and simple approach to localizing the focal spot in three dimensions with a single scan. Magn Reson Med 76:803-813, 2016. © 2015 Wiley Periodicals, Inc.

Sampling strategies for subsampled segmented EPI PRF thermometry in MR guided high intensity focused ultrasound

PURPOSE

To investigate k-space subsampling strategies to achieve fast, large field-of-view (FOV) temperature monitoring using segmented echo planar imaging (EPI) proton resonance frequency shift thermometry for MR guided high intensity focused ultrasound (MRgHIFU) applications.

METHODS

Five different k-space sampling approaches were investigated, varying sample spacing (equally vs nonequally spaced within the echo train), sampling density (variable sampling density in zero, one, and two dimensions), and utilizing sequential or centric sampling. Three of the schemes utilized sequential sampling with the sampling density varied in zero, one, and two dimensions, to investigate sampling the k-space center more frequently. Two of the schemes utilized centric sampling to acquire the k-space center with a longer echo time for improved phase measurements, and vary the sampling density in zero and two dimensions, respectively. Phantom experiments and a theoretical point spread function analysis were performed to investigate their performance. Variable density sampling in zero and two dimensions was also implemented in a non-EPI GRE pulse sequence for comparison. All subsampled data were reconstructed with a previously described temporally constrained reconstruction (TCR) algorithm.

RESULTS

The accuracy of each sampling strategy in measuring the temperature rise in the HIFU focal spot was measured in terms of the root-mean-square-error (RMSE) compared to fully sampled "truth." For the schemes utilizing sequential sampling, the accuracy was found to improve with the dimensionality of the variable density sampling, giving values of 0.65 °C, 0.49 °C, and 0.35 °C for density variation in zero, one, and two dimensions, respectively. The schemes utilizing centric sampling were found to underestimate the temperature rise, with RMSE values of 1.05 °C and 1.31 °C, for variable density sampling in zero and two dimensions, respectively. Similar subsampling schemes with variable density sampling implemented in zero and two dimensions in a non-EPI GRE pulse sequence both resulted in accurate temperature measurements (RMSE of 0.70 °C and 0.63 °C, respectively). With sequential sampling in the described EPI implementation, temperature monitoring over a 192×144×135 mm3 FOV with a temporal resolution of 3.6 s was achieved, while keeping the RMSE compared to fully sampled "truth" below 0.35 °C.

CONCLUSIONS

When segmented EPI readouts are used in conjunction with k-space subsampling for MR thermometry applications, sampling schemes with sequential sampling, with or without variable density sampling, obtain accurate phase and temperature measurements when using a TCR reconstruction algorithm. Improved temperature measurement accuracy can be achieved with variable density sampling. Centric sampling leads to phase bias, resulting in temperature underestimations.

An evaluation of prospective motion correction (PMC) for high resolution quantitative MRI

Quantitative imaging aims to provide in vivo neuroimaging biomarkers with high research and diagnostic value that are sensitive to underlying tissue microstructure. In order to use these data to examine intra-cortical differences or to define boundaries between different myelo-architectural areas, high resolution data are required. The quality of such measurements is degraded in the presence of motion hindering insight into brain microstructure. Correction schemes are therefore vital for high resolution, whole brain coverage approaches that have long acquisition times and greater sensitivity to motion. Here we evaluate the use of prospective motion correction (PMC) via an optical tracking system to counter intra-scan motion in a high resolution (800 μm isotropic) multi-parameter mapping (MPM) protocol. Data were acquired on six volunteers using a 2 × 2 factorial design permuting the following conditions: PMC on/off and motion/no motion. In the presence of head motion, PMC-based motion correction considerably improved the quality of the maps as reflected by fewer visible artifacts and improved consistency. The precision of the maps, parameterized through the coefficient of variation in cortical sub-regions, showed improvements of 11-25% in the presence of deliberate head motion. Importantly, in the absence of motion the PMC system did not introduce extraneous artifacts into the quantitative maps. The PMC system based on optical tracking offers a robust approach to minimizing motion artifacts in quantitative anatomical imaging without extending scan times. Such a robust motion correction scheme is crucial in order to achieve the ultra-high resolution required of quantitative imaging for cutting edge in vivo histology applications.

Prospective motion correction of 3D echo-planar imaging data for functional MRI using optical tracking

We evaluated the performance of an optical camera based prospective motion correction (PMC) system in improving the quality of 3D echo-planar imaging functional MRI data. An optical camera and external marker were used to dynamically track the head movement of subjects during fMRI scanning. PMC was performed by using the motion information to dynamically update the sequence's RF excitation and gradient waveforms such that the field-of-view was realigned to match the subject's head movement. Task-free fMRI experiments on five healthy volunteers followed a 2 × 2 × 3 factorial design with the following factors: PMC on or off; 3.0 mm or 1.5 mm isotropic resolution; and no, slow, or fast head movements. Visual and motor fMRI experiments were additionally performed on one of the volunteers at 1.5 mm resolution comparing PMC on vs PMC off for no and slow head movements. Metrics were developed to quantify the amount of motion as it occurred relative to k-space data acquisition. The motion quantification metric collapsed the very rich camera tracking data into one scalar value for each image volume that was strongly predictive of motion-induced artifacts. The PMC system did not introduce extraneous artifacts for the no motion conditions and improved the time series temporal signal-to-noise by 30% to 40% for all combinations of low/high resolution and slow/fast head movement relative to the standard acquisition with no prospective correction. The numbers of activated voxels (p < 0.001, uncorrected) in both task-based experiments were comparable for the no motion cases and increased by 78% and 330%, respectively, for PMC on versus PMC off in the slow motion cases. The PMC system is a robust solution to decrease the motion sensitivity of multi-shot 3D EPI sequences and thereby overcome one of the main roadblocks to their widespread use in fMRI studies.

•

Data quality was assessed in terms of tSNR, RMSE, and task fMRI t-statistics.

•

Correction improved tSNR by 30% to 40%.

•

Percent reduction in RMSE as a function of motion level was characterized.

•

Number of active voxels for motor and visual tasks increased significantly.

Model predictive filtering MR thermometry: Effects of model inaccuracies, k-space reduction factor, and temperature increase rate

PURPOSE

Evaluate effects of model parameter inaccuracies (thermal conductivity, k, and ultrasound power deposition density, Q), k-space reduction factor (R), and rate of temperature increase ( T˙) in a thermal model-based reconstruction for MR-thermometry during focused-ultrasound heating.

METHODS

Simulations and ex vivo experiments were performed to investigate the accuracy of the thermal model and the model predictive filtering (MPF) algorithm for varying R and T˙, and their sensitivity to errors in k and Q. Ex vivo data was acquired with a segmented EPI pulse sequence to achieve large field-of-view (192 × 162 × 96 mm) four-dimensional temperature maps with high spatiotemporal resolution (1.5 × 1.5 × 2.0 mm, 1.7 s).

RESULTS

In the simulations, 50% errors in k and Q resulted in maximum temperature root mean square errors (RMSE) of 6 °C for model only and 3 °C for MPF. Using recently developed methods, estimates of k and Q were accurate to within 3%. The RMSE between MPF and true temperature increased with R and T˙. In the ex vivo study the RMSE remained below 0.7 °C for R ranging from 4 to 12 and T˙ of 0.28-0.75 °C/s.

CONCLUSION

Errors in MPF temperatures occur due to errors in k and Q. These MPF temperature errors increase with increase in R and T˙, but are smaller than those obtained using the thermal model alone.

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.

Approaching the subject of violence: a response-based approach to working with men who have abused others

Traditional effects-based approaches to therapeutic work with men who have abused others often attempt to intervene by correcting personal deficits assumed to be causing the violence. This not only creates a hierarchical counseling relationship but also can inadvertently excuse aggressive actions. In this article, we outline a response-based alternative that emphasizes questions of choice, agency, and volition within a collaborative therapeutic relationship. Rather than impose external correction, we pay attention to details of how men describe their violent acts and position themselves as agents of those acts as we work toward supporting them in their own acts of self-correction.

In vivo evaluation of a breast-specific magnetic resonance guided focused ultrasound system in a goat udder model

PURPOSE

This work further evaluates the functionality, efficacy, and safety of a new breast-specific magnetic resonance guided high intensity focused ultrasound (MRgFUS) system in an in vivo goat udder model.

METHODS

Eight female goats underwent an MRgFUS ablation procedure using the breast-specific MRgFUS system. Tissue classification was achieved through the 3D magnetic resonance imaging (MRI) acquisition of several contrasts (T1w, T2w, PDw, 3-point Dixon). The MRgFUS treatment was performed with a grid trajectory executed in one or two planes within the glandular tissue of the goat udder. Temperature was monitored using a 3D proton resonance frequency (PRF) MRI technique. Delayed contrast enhanced-MR images were acquired immediately and 14 days post MRgFUS treatment. A localized tissue excision was performed in one animal and histological analysis was performed. Animals were available for adoption at the conclusion of the study.

RESULTS

The breast-specific MRgFUS system was able to ablate regions ranging in size from 0.4 to 3.6 cm(3) in the goat udder model. Tissue damage was confirmed through the correlation of thermal dose measurements obtained with realtime 3D MR thermometry to delayed contrast enhanced-MR images immediately after the treatment and 14 days postablation. In general, lesions were longer in the ultrasound propagation direction, which is consistent with the dimensions of the ultrasound focal spot. Thermal dose volumes had better agreement with nonenhancing areas of the DCE-MRI images obtained 14 days after the MRgFUS treatment.

CONCLUSIONS

The system was able to successfully ablate lesions up to 3.6 cm(3). The thermal dose volume was found to correlate better with the 14-day postablation nonenhancing delayed contrast enhanced-MR image volumes. While the goat udder is not an ideal model for the human breast, this study has proven the feasibility of using this system on a wide variety of udder shapes and sizes, demonstrating the flexibility that would be required in order to treat human subjects.

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.

Effects of MRTI sampling characteristics on estimation of HIFU SAR and tissue thermal diffusivity

While the non-invasive and three-dimensional nature of magnetic-resonance temperature imaging (MRTI) makes it a valuable tool for high-intensity focused ultrasound (HIFU) treatments, random and systematic errors in MRTI measurements may propagate into temperature-based parameter estimates used for pretreatment planning. This study assesses the MRTI effects of zero-mean Gaussian noise (SD = 0.0-2.0 °C), temporal sampling (tacq = 1.0-8.0 s), and spatial averaging (Res = 0.5-2.0 mm isotropic) on HIFU temperature measurements and temperature-based estimates of the amplitude and full width half maximum (FWHM) of the HIFU specific absorption rate and of tissue thermal diffusivity. The ultrasound beam used in simulations and ex vivo pork loin experiments has lateral and axial FWHM dimensions of 1.4 mm and 7.9 mm respectively. For spatial averaging simulations, beams with lateral FWHM varying from 1.2-2.2 mm are also assessed. Under noisy conditions, parameter estimates are improved by fitting to data from larger voxel regions. Varying the temporal sampling results in minimal changes in measured temperatures (<2% change) and parameter estimates (<5% change). For the HIFU beams studied, a spatial resolution of 1 × 1 × 3 mm(3) or smaller is required to keep errors in temperature and all estimated parameters less than 10%. By quantifying the errors associated with these sampling characteristics, this work provides researchers with appropriate MRTI conditions for obtaining estimates of parameters essential to pretreatment modeling of HIFU 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.

Toward real-time temperature monitoring in fat and aqueous tissue during magnetic resonance-guided high-intensity focused ultrasound using a three-dimensional proton resonance frequency T1 method

PURPOSE

To present a three-dimensional (3D) segmented echoplanar imaging (EPI) pulse sequence implementation that provides simultaneously the proton resonance frequency shift temperature of aqueous tissue and the longitudinal relaxation time (T1 ) of fat during thermal ablation.

METHODS

The hybrid sequence was implemented by combining a 3D segmented flyback EPI sequence, the extended two-point Dixon fat and water separation, and the double flip angle T1 mapping techniques. High-intensity focused ultrasound (HIFU) heating experiments were performed at three different acoustic powers on excised human breast fat embedded in ex vivo porcine muscle. Furthermore, T1 calibrations with temperature in four different excised breast fat samples were performed, yielding an estimate of the average and variation of dT1 /dT across subjects.

RESULTS

The water only images were used to mask the complex original data before computing the proton resonance frequency shift. T1 values were calculated from the fat-only images. The relative temperature coefficients were found in five fat tissue samples from different patients and ranged from 1.2% to 2.6%/°C.

CONCLUSION

The results demonstrate the capability of real-time simultaneous temperature mapping in aqueous tissue and T1 mapping in fat during HIFU ablation, providing a potential tool for treatment monitoring in organs with large fat content, such as the breast.

Design and characterization of a laterally mounted phased-array transducer breast-specific MRgHIFU device with integrated 11-channel receiver array

PURPOSE

This work presents the design and preliminary evaluation of a new laterally mounted phased-array MRI-guided high-intensity focused ultrasound (MRgHIFU) system with an integrated 11-channel phased-array radio frequency (RF) coil intended for breast cancer treatment. The design goals for the system included the ability to treat the majority of tumor locations, to increase the MR image's signal-to-noise ratio (SNR) throughout the treatment volume and to provide adequate comfort for the patient.

METHODS

In order to treat the majority of the breast volume, the device was designed such that the treated breast is suspended in a 17-cm diameter treatment cylinder. A laterally shooting 1-MHz, 256-element phased-array ultrasound transducer with flexible positioning is mounted outside the treatment cylinder. This configuration achieves a reduced water volume to minimize RF coil loading effects, to position the coils closer to the breast for increased signal sensitivity, and to reduce the MR image noise associated with using water as the coupling fluid. This design uses an 11-channel phased-array RF coil that is placed on the outer surface of the cylinder surrounding the breast. Mechanical positioning of the transducer and electronic steering of the focal spot enable placement of the ultrasound focus at arbitrary locations throughout the suspended breast. The treatment platform allows the patient to lie prone in a face-down position. The system was tested for comfort with 18 normal volunteers and SNR capabilities in one normal volunteer and for heating accuracy and stability in homogeneous phantom and inhomogeneous ex vivo porcine tissue.

RESULTS

There was a 61% increase in mean relative SNR achieved in a homogeneous phantom using the 11-channel RF coil when compared to using only a single-loop coil around the chest wall. The repeatability of the system's energy delivery in a single location was excellent, with less than 3% variability between repeated temperature measurements at the same location. The execution of a continuously sonicated, predefined 48-point, 8-min trajectory path resulted in an ablation volume of 8.17 cm(3), with one standard deviation of 0.35 cm(3) between inhomogeneous ex vivo tissue samples. Comfort testing resulted in negligible side effects for all volunteers.

CONCLUSIONS

The initial results suggest that this new device will potentially be suitable for MRgHIFU treatment in a wide range of breast sizes and tumor locations.

Toward real-time availability of 3D temperature maps created with temporally constrained reconstruction

PURPOSE

To extend the previously developed temporally constrained reconstruction (TCR) algorithm to allow for real-time availability of three-dimensional (3D) temperature maps capable of monitoring MR-guided high intensity focused ultrasound applications.

METHODS

A real-time TCR (RT-TCR) algorithm is developed that only uses current and previously acquired undersampled k-space data from a 3D segmented EPI pulse sequence, with the image reconstruction done in a graphics processing unit implementation to overcome computation burden. Simulated and experimental data sets of HIFU heating are used to evaluate the performance of the RT-TCR algorithm.

RESULTS

The simulation studies demonstrate that the RT-TCR algorithm has subsecond reconstruction time and can accurately measure HIFU-induced temperature rises of 20°C in 15 s for 3D volumes of 16 slices (RMSE = 0.1°C), 24 slices (RMSE = 0.2°C), and 32 slices (RMSE = 0.3°C). Experimental results in ex vivo porcine muscle demonstrate that the RT-TCR approach can reconstruct temperature maps with 192 × 162 × 66 mm 3D volume coverage, 1.5 × 1.5 × 3.0 mm resolution, and 1.2-s scan time with an accuracy of ±0.5°C.

CONCLUSION

The RT-TCR algorithm offers an approach to obtaining large coverage 3D temperature maps in real-time for monitoring MR-guided high intensity focused ultrasound treatments.

An 11-channel radio frequency phased array coil for magnetic resonance guided high-intensity focused ultrasound of the breast

In this study, a radio frequency phased array coil was built to image the breast in conjunction with a magnetic resonance guided high-intensity focused ultrasound (MRgHIFU) device designed specifically to treat the breast in a treatment cylinder with reduced water volume. The MRgHIFU breast coil was comprised of a 10-channel phased array coil placed around an MRgHIFU treatment cylinder where nearest-neighbor decoupling was achieved with capacitive decoupling in a shared leg. In addition a single loop coil was placed at the chest wall making a total of 11 channels. The radio frequency coil array design presented in this work was chosen based on ease of implementation, increased visualization into the treatment cylinder, image reconstruction speed, temporal resolution, and resulting signal-to-noise ratio profiles. This work presents a dedicated 11-channel coil for imaging of the breast tissue in the MRgHIFU setup without obstruction of the ultrasound beam and, specifically, compares its performance in signal-to-noise, overall imaging time, and temperature measurement accuracy to that of the standard single chest-loop coil typically used in breast MRgHIFU.

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

Denaturation of macromolecules within the tissues is believed to be the major factor contributing to the damage of tissues upon hyperthermia. As a result, the value of the spin-lattice relaxation time T1 of the tissue water, which is related to the translational and rotational rates of water, represents an intrinsic probe for investigating structural changes in tissues at high temperature. Therefore, the goal of this work is to investigate whether the simultaneous measurement of temperature and T1 using a hybrid proton resonance frequency (PRF)-T1 measurement technique can be used to detect irreversible changes in T1 that might be indicative of tissue damage. A new hybrid PRF-T1 sequence was implemented based on the variable flip angle driven-equilibrium single-pulse observation (DESPOT)1 method from a standard three dimensional segmented echo-planar imaging sequence by alternating two flip angles from measurement to measurement. The structural changes of the heated tissue volumes were analyzed based on the derived T1 values and the corresponding PRF temperatures. Using the hybrid PRF-T1 technique, we demonstrate that the change of spin lattice relaxation time T1 is reversible with temperature for low thermal dose (thermal dose ≤ 240 cumulative equivalent minutes [CEM] 43°C) and irreversible with temperature after significant accumulation of thermal dose in ex vivo chicken breast tissue. These results suggest that the hybrid PRF-T1 method may be a potentially powerful tool to investigate the extent and mechanism of heat damage of biological tissues.

Hybrid proton resonance frequency/T1 technique for simultaneous temperature monitoring in adipose and aqueous tissues

Thermal therapy procedures being carried out under MR guidance would be safer if temperature changes could be accurately monitored in both water-based and fat-based tissues. To this end, we present a hybrid proton resonance frequency (PRF)/T(1) approach for simultaneously measuring PRF shift temperatures in water-based tissues and T(1) changes in fat-based tissues. The hybrid PRF/T(1) sequence is a standard radiofrequency spoiled gradient echo sequence executed in a dynamic mode with two flip angles alternating every time frame. The PRF information is extracted every time frame using the image phase in the standard approach, and the T(1) information is extracted every two time frames using a variable flip angle approach. Simulation studies, ex vivo high intensity focused ultrasound heating experiments, and in vivo stability experiments were performed to test the feasibility of the approach. The results indicate that the hybrid PRF/T(1) approach provides PRF temperature maps of the same quality as those obtained by traditional PRF methods while simultaneously being able to track T(1) changes in fat-based tissues. Although several potential error sources exist for the T(1) measurements, the approach is a promising start toward realizing quantitative temperature measurements in both water-based and fat-based tissues.

The effect of electronically steering a phased array ultrasound transducer on near-field tissue heating

PURPOSE

This study presents the results obtained from both simulation and experimental techniques that show the effect of mechanically or electronically steering a phased array transducer on proximal tissue heating.

METHODS

The thermal response of a nine-position raster and a 16-mm diameter circle scanning trajectory executed through both electronic and mechanical scanning was evaluated in computer simulations and experimentally in a homogeneous tissue-mimicking phantom. Simulations were performed using power deposition maps obtained from the hybrid angular spectrum (HAS) method and applying a finite-difference approximation of the Pennes' bioheat transfer equation for the experimentally used transducer and also for a fully sampled transducer to demonstrate the effect of acoustic window, ultrasound beam overlap and grating lobe clutter on near-field heating.

RESULTS

Both simulation and experimental results show that electronically steering the ultrasound beam for the two trajectories using the 256-element phased array significantly increases the thermal dose deposited in the near-field tissues when compared with the same treatment executed through mechanical steering only. In addition, the individual contributions of both beam overlap and grating lobe clutter to the near-field thermal effects were determined through comparing the simulated ultrasound beam patterns and resulting temperature fields from mechanically and electronically steered trajectories using the 256-randomized element phased array transducer to an electronically steered trajectory using a fully sampled transducer with 40 401 phase-adjusted sample points.

CONCLUSIONS

Three distinctly different three distinctly different transducers were simulated to analyze the tradeoffs of selected transducer design parameters on near-field heating. Careful consideration of design tradeoffs and accurate patient treatment planning combined with thorough monitoring of the near-field tissue temperature will help to ensure patient safety during an MRgHIFU treatment.

Reconstruction of fully three-dimensional high spatial and temporal resolution MR temperature maps for retrospective applications

Many areas of MR-guided thermal therapy research would benefit from temperature maps with high spatial and temporal resolution that cover a large three-dimensional volume. This article describes an approach to achieve these goals, which is suitable for research applications where retrospective reconstruction of the temperature maps is acceptable. The method acquires undersampled data from a modified three-dimensional segmented echo-planar imaging sequence and creates images using a temporally constrained reconstruction algorithm. The three-dimensional images can be zero-filled to arbitrarily small voxel spacing in all directions and then converted into temperature maps using the standard proton resonance frequency shift technique. During high intensity focused ultrasound heating experiments, the proposed method was used to obtain temperature maps with 1.5 mm × 1.5 mm × 3.0 mm resolution, 288 mm × 162 mm × 78 mm field of view, and 1.7 s temporal resolution. The approach is validated to demonstrate that it can accurately capture the spatial characteristics and time dynamics of rapidly changing high intensity focused ultrasound-induced temperature distributions. Example applications from MR-guided high intensity focused ultrasound research are shown to demonstrate the benefits of the large coverage fully three-dimensional temperature maps, including characterization of volumetric heating trajectories and near- and far-field heating.

Ultrasound-mediated tumor imaging and nanotherapy using drug loaded, block copolymer stabilized perfluorocarbon nanoemulsions

Perfluorocarbon nanoemulsions can deliver lipophilic therapeutic agents to solid tumors and simultaneously provide for monitoring nanocarrier biodistribution via ultrasonography and/or (19)F MRI. In the first generation of block copolymer stabilized perfluorocarbon nanoemulsions, perfluoropentane (PFP) was used as the droplet forming compound. Although manifesting excellent therapeutic and ultrasound imaging properties, PFP nanoemulsions were unstable at storage, difficult to handle, and underwent hard to control phenomenon of irreversible droplet-to-bubble transition upon injection. To solve the above problems, perfluoro-15-crown-5-ether (PFCE) was used as a core forming compound in the second generation of block copolymer stabilized perfluorocarbon nanoemulsions. PFCE nanodroplets manifest both ultrasound and fluorine ((19)F) MR contrast properties, which allows using multimodal imaging and (19)F MR spectroscopy for monitoring nanodroplet pharmacokinetics and biodistribution. In the present paper, acoustic, imaging, and therapeutic properties of unloaded and paclitaxel (PTX) loaded PFCE nanoemulsions are reported. As manifested by the (19)F MR spectroscopy, PFCE nanodroplets are long circulating, with about 50% of the injected dose remaining in circulation 2h after the systemic injection. Sonication with 1-MHz therapeutic ultrasound triggered reversible droplet-to-bubble transition in PFCE nanoemulsions. Microbubbles formed by acoustic vaporization of nanodroplets underwent stable cavitation. The nanodroplet size (200nm to 350nm depending on a type of the shell and conditions of emulsification) as well as long residence in circulation favored their passive accumulation in tumor tissue that was confirmed by ultrasonography. In the breast and pancreatic cancer animal models, ultrasound-mediated therapy with paclitaxel-loaded PFCE nanoemulsions showed excellent therapeutic properties characterized by tumor regression and suppression of metastasis. Anticipated mechanisms of the observed effects are discussed.

The effects of spatial sampling choices on MR temperature measurements

The purpose of this article is to quantify the effects that spatial sampling parameters have on the accuracy of magnetic resonance temperature measurements during high intensity focused ultrasound treatments. Spatial resolution and position of the sampling grid were considered using experimental and simulated data for two different types of high intensity focused ultrasound heating trajectories (a single point and a 4-mm circle) with maximum measured temperature and thermal dose volume as the metrics. It is demonstrated that measurement accuracy is related to the curvature of the temperature distribution, where regions with larger spatial second derivatives require higher resolution. The location of the sampling grid relative temperature distribution has a significant effect on the measured values. When imaging at 1.0 × 1.0 × 3.0 mm(3) resolution, the measured values for maximum temperature and volume dosed to 240 cumulative equivalent minutes (CEM) or greater varied by 17% and 33%, respectively, for the single-point heating case, and by 5% and 18%, respectively, for the 4-mm circle heating case. Accurate measurement of the maximum temperature required imaging at 1.0 × 1.0 × 3.0 mm(3) resolution for the single-point heating case and 2.0 × 2.0 × 5.0 mm(3) resolution for the 4-mm circle heating case.

Model predictive filtering for improved temporal resolution in MRI temperature imaging

A novel method for reconstructing MRI temperature maps from undersampled data is presented. The method, model predictive filtering, combines temperature predictions from a preidentified thermal model with undersampled k-space data to create temperature maps in near real time. The model predictive filtering algorithm was implemented in three ways: using retrospectively undersampled k-space data from a fully sampled two-dimensional gradient echo (GRE) sequence (reduction factors R = 2.7 to R = 7.1), using actually undersampled data from a two-dimensional GRE sequence (R = 4.8), and using actually undersampled data from a three-dimensional GRE sequence (R = 12.1). Thirty-nine high-intensity focused ultrasound heating experiments were performed under MRI monitoring to test the model predictive filtering technique against the current gold standard for MR temperature mapping, the proton resonance frequency shift method. For both of the two-dimensional implementations, the average error over the five hottest voxels from the hottest time frame remained between +/-0.8 degrees C and the temperature root mean square error over a 24 x 7 x 3 x 25-voxel region of interest remained below 0.35 degrees C. The largest errors for the three-dimensional implementation were slightly worse: -1.4 degrees C for the mean error of the five hottest voxels and 0.61 degrees C for the temperature root mean square error.

Temporally constrained reconstruction applied to MRI temperature data

The monitoring of thermal ablation procedures would benefit from an acceleration in the rate at which MRI temperature maps are acquired. Constrained reconstruction techniques have been shown to be capable of generating high quality images using only a fraction of the k-space data. Here, we present a temporally constrained reconstruction (TCR) algorithm applied to proton resonance frequency shift (PRF) data. The algorithm generates images from undersampled data by iteratively minimizing a cost function. The unique challenges of using an iterative constrained reconstruction technique to generate real-time images were addressed. For a set of eight heating experiments on ex vivo porcine tissue, a maximum reduction factor of 4 was achieved while keeping the root mean square error (RMSE) of the temperature below 0.5 degrees C. For a set of three heating experiments on in vivo canine muscle tissue, the maximum reduction factor achieved was 3 while keeping the temperature RMSE below 1.0 degrees C. At these reduction factors, the TCR algorithm underpredicted the thermal dose by an average of 6% for the ex vivo data and 28% for the in vivo data. Compared with sliding window and low resolution reconstructions, the RMSE of the TCR algorithm was significantly lower (P < 0.05 in all cases).