Vespa: Integrated applications for RF pulse design, spectral simulation and MRS data analysis.

PURPOSE: The Vespa package (Versatile Simulation, Pulses, and Analysis) is described and demonstrated. It provides workflows for developing and optimizing linear combination modeling (LCM) fitting for 1 H MRS data using intuitive graphical user interface interfaces for RF pulse design, spectral simulation, and MRS data analysis. Command line interfaces for embedding workflows in MR manufacturer platforms and utilities for synthetic dataset creation are included. Complete provenance is maintained for all steps in workflows. THEORY AND METHODS: Vespa is written in Python for compatibility across operating systems. It embeds the PyGAMMA spectral simulation library for spectral simulation. Multiprocessing methods accelerate processing and visualization. Applications use the Vespa database for results storage and cross-application access. Three projects demonstrate pulse, sequence, simulation, and data analysis workflows: (1) short TE semi-LASER single-voxel spectroscopy (SVS) LCM fitting, (2) optimizing MEGA-PRESS (MEscher-GArwood Point RESolved Spectroscopy) flip angle and LCM fitting, and (3) creating a synthetic short TE dataset. RESULTS: The LCM workflows for in vivo basis set creation and spectral analysis showed reasonable results for both the short TE semi-LASER and MEGA-PRESS. Examples of pulses, simulations, and data fitting are shown in Vespa application interfaces for various steps to demonstrate the interactive workflow. CONCLUSION: Vespa provides an efficient and extensible platform for characterizing RF pulses, pulse design, spectral simulation optimization, and automated LCM fitting via an interactive platform. Modular design and command line interface make it easy to embed in other platforms. As open source, it is free to the MRS community for use and extension. Vespa source code and documentation are available through GitHub.

Influence of editing pulse flip angle on J-difference MR spectroscopy.

PURPOSE: To investigate the editing-pulse flip angle (FA) dependence of editing efficiency and ultimately to maximize the edited signal of commonly edited MR spectroscopy (MRS) signals, such as gamma-aminobutyric acid (GABA) and lactate. METHODS: Density-matrix simulations were performed for a range of spin systems to find the editing-pulse FA for maximal editing efficiency. Simulations were confirmed by phantom experiments and in vivo measurements in 10 healthy participants using a 3T Philips scanner. Four MEGA-PRESS in vivo measurements targeting GABA+ and lactate were performed, comparing the conventional editing-pulse FA (FA = 180°) to the optimal one suggested by simulations (FA = 210°). RESULTS: Simulations and phantom experiments show that edited GABA and lactate signals are maximal at FA = 210°. Compared to conventional editing (FA = 180°), in vivo signals from GABA+ and lactate signals increase on average by 8.5% and 9.3%, respectively. CONCLUSION: Increasing the FA of editing-pulses in the MEGA-PRESS experiment from 180° to 210° increases the edited signals from GABA+ and lactate by about 9% in vivo.

NIfTI-MRS: A standard data format for magnetic resonance spectroscopy.

PURPOSE: Multiple data formats in the MRS community currently hinder data sharing and integration. NIfTI-MRS is proposed as a standard spectroscopy data format, implemented as an extension to the Neuroimaging informatics technology initiative (NIfTI) format. This standardized format can facilitate data sharing and algorithm development as well as ease integration of MRS analysis alongside other imaging modalities. METHODS: A file format using the NIfTI header extension framework incorporates essential spectroscopic metadata and additional encoding dimensions. A detailed description of the specification is provided. An open-source command-line conversion program is implemented to convert single-voxel and spectroscopic imaging data to NIfTI-MRS. Visualization of data in NIfTI-MRS is provided by development of a dedicated plugin for FSLeyes, the FMRIB Software Library (FSL) image viewer. RESULTS: Online documentation and 10 example datasets in the proposed format are provided. Code examples of NIfTI-MRS readers are implemented in common programming languages. Conversion software, spec2nii, currently converts 14 formats where data is stored in image-space to NIfTI-MRS, including Digital Imaging and Communications in Medicine (DICOM) and vendor proprietary formats. CONCLUSION: NIfTI-MRS aims to solve issues arising from multiple data formats being used in the MRS community. Through a single conversion point, processing and analysis of MRS data are simplified, thereby lowering the barrier to use of MRS. Furthermore, it can serve as the basis for open data sharing, collaboration, and interoperability of analysis programs. Greater standardization and harmonization become possible. By aligning with the dominant format in neuroimaging, NIfTI-MRS enables the use of mature tools present in the imaging community, demonstrated in this work by using a dedicated imaging tool, FSLeyes, for visualization.

Lumbar intervertebral disc diurnal deformations and T2 and T1rho relaxation times vary by spinal level and disc region.

PURPOSE: Magnetic resonance imaging (MRI) is routinely used to evaluate spine pathology; however, standard imaging findings weakly correlate to low back pain. Abnormal disc mechanical function is implicated as a cause of back pain but is not assessed using standard clinical MRI. Our objective was to utilize our established MRI protocol for measuring disc function to quantify disc mechanical function in a healthy cohort. METHODS: We recruited young, asymptomatic volunteers (6 male/6 female; age 18-30 years; BMI < 30) and used MRI to determine how diurnal deformations in disc height, volume, and perimeter were affected by spinal level, disc region, MRI biomarkers of disc health (T2, T1rho), and Pfirrmann grade. RESULTS: Lumbar discs deformed by a mean of -6.1% (95% CI: -7.6%, -4.7%) to -8.0% (CI: -10.6%, -5.4%) in height and -5.4% (CI: -7.6%, -3.3%) to -8.5% (CI: -11.0%, -6.0%) in volume from AM to PM across spinal levels. Regional deformations were more uniform in cranial lumbar levels and concentrated posteriorly in the caudal levels, reaching a maximum of 13.1% at L5-S1 (CI:-16.1%, -10.2%). T2 and T1rho relaxation times were greatest in the nucleus and varied circumferentially within the annulus. T2 relaxation times were greatest at the most cranial spinal levels and decreased caudally. In this young healthy cohort, we identified a weak association between nucleus T2 and the diurnal change in the perimeter. CONCLUSIONS: Spinal level is a key factor in determining regional disc deformations. Interestingly, deformations were concentrated in the posterior regions of caudal discs where disc herniation is most prevalent.

Advanced magnetic resonance spectroscopic neuroimaging: Experts' consensus recommendations.

Magnetic resonance spectroscopic imaging (MRSI) offers considerable promise for monitoring metabolic alterations associated with disease or injury; however, to date, these methods have not had a significant impact on clinical care, and their use remains largely confined to the research community and a limited number of clinical sites. The MRSI methods currently implemented on clinical MRI instruments have remained essentially unchanged for two decades, with only incremental improvements in sequence implementation. During this time, a number of technological developments have taken place that have already greatly benefited the quality of MRSI measurements within the research community and which promise to bring advanced MRSI studies to the point where the technique becomes a true imaging modality, while making the traditional review of individual spectra a secondary requirement. Furthermore, the increasing use of biomedical MR spectroscopy studies has indicated clinical areas where advanced MRSI methods can provide valuable information for clinical care. In light of this rapidly changing technological environment and growing understanding of the value of MRSI studies for biomedical studies, this article presents a consensus from a group of experts in the field that reviews the state-of-the-art for clinical proton MRSI studies of the human brain, recommends minimal standards for further development of vendor-provided MRSI implementations, and identifies areas which need further technical development.

Across-vendor standardization of semi-LASER for single-voxel MRS at 3T.

The semi-adiabatic localization by adiabatic selective refocusing (sLASER) sequence provides single-shot full intensity signal with clean localization and minimal chemical shift displacement error and was recommended by the international MRS Consensus Group as the preferred localization sequence at high- and ultra-high fields. Across-vendor standardization of the sLASER sequence at 3 tesla has been challenging due to the B1 requirements of the adiabatic inversion pulses and maximum B1 limitations on some platforms. The aims of this study were to design a short-echo sLASER sequence that can be executed within a B1 limit of 15 µT by taking advantage of gradient-modulated RF pulses, to implement it on three major platforms and to evaluate the between-vendor reproducibility of its perfomance with phantoms and in vivo. In addition, voxel-based first and second order B0 shimming and voxel-based B1 adjustments of RF pulses were implemented on all platforms. Amongst the gradient-modulated pulses considered (GOIA, FOCI and BASSI), GOIA-WURST was identified as the optimal refocusing pulse that provides good voxel selection within a maximum B1 of 15 µT based on localization efficiency, contamination error and ripple artifacts of the inversion profile. An sLASER sequence (30 ms echo time) that incorporates VAPOR water suppression and 3D outer volume suppression was implemented with identical parameters (RF pulse type and duration, spoiler gradients and inter-pulse delays) on GE, Philips and Siemens and generated identical spectra on the GE 'Braino' phantom between vendors. High-quality spectra were consistently obtained in multiple regions (cerebellar white matter, hippocampus, pons, posterior cingulate cortex and putamen) in the human brain across vendors (5 subjects scanned per vendor per region; mean signal-to-noise ratio > 33; mean water linewidth between 6.5 Hz to 11.4 Hz). The harmonized sLASER protocol is expected to produce high reproducibility of MRS across sites thereby allowing large multi-site studies with clinical cohorts.

Contribution of macromolecules to brain 1 H MR spectra: Experts' consensus recommendations.

Proton MR spectra of the brain, especially those measured at short and intermediate echo times, contain signals from mobile macromolecules (MM). A description of the main MM is provided in this consensus paper. These broad peaks of MM underlie the narrower peaks of metabolites and often complicate their quantification but they also may have potential importance as biomarkers in specific diseases. Thus, separation of broad MM signals from low molecular weight metabolites enables accurate determination of metabolite concentrations and is of primary interest in many studies. Other studies attempt to understand the origin of the MM spectrum, to decompose it into individual spectral regions or peaks and to use the components of the MM spectrum as markers of various physiological or pathological conditions in biomedical research or clinical practice. The aim of this consensus paper is to provide an overview and some recommendations on how to handle the MM signals in different types of studies together with a list of open issues in the field, which are all summarized at the end of the paper.

Methodological consensus on clinical proton MRS of the brain: Review and recommendations.

Proton MRS (1 H MRS) provides noninvasive, quantitative metabolite profiles of tissue and has been shown to aid the clinical management of several brain diseases. Although most modern clinical MR scanners support MRS capabilities, routine use is largely restricted to specialized centers with good access to MR research support. Widespread adoption has been slow for several reasons, and technical challenges toward obtaining reliable good-quality results have been identified as a contributing factor. Considerable progress has been made by the research community to address many of these challenges, and in this paper a consensus is presented on deficiencies in widely available MRS methodology and validated improvements that are currently in routine use at several clinical research institutions. In particular, the localization error for the PRESS localization sequence was found to be unacceptably high at 3 T, and use of the semi-adiabatic localization by adiabatic selective refocusing sequence is a recommended solution. Incorporation of simulated metabolite basis sets into analysis routines is recommended for reliably capturing the full spectral detail available from short TE acquisitions. In addition, the importance of achieving a highly homogenous static magnetic field (B0 ) in the acquisition region is emphasized, and the limitations of current methods and hardware are discussed. Most recommendations require only software improvements, greatly enhancing the capabilities of clinical MRS on existing hardware. Implementation of these recommendations should strengthen current clinical applications and advance progress toward developing and validating new MRS biomarkers for clinical use.

A convolutional neural network to filter artifacts in spectroscopic MRI.

PURPOSE: Proton MRSI is a noninvasive modality capable of generating volumetric maps of in vivo tissue metabolism without the need for ionizing radiation or injected contrast agent. Magnetic resonance spectroscopic imaging has been shown to be a viable imaging modality for studying several neuropathologies. However, a key hurdle in the routine clinical adoption of MRSI is the presence of spectral artifacts that can arise from a number of sources, possibly leading to false information. METHODS: A deep learning model was developed that was capable of identifying and filtering out poor quality spectra. The core of the model used a tiled convolutional neural network that analyzed frequency-domain spectra to detect artifacts. RESULTS: When compared with a panel of MRS experts, our convolutional neural network achieved high sensitivity and specificity with an area under the curve of 0.95. A visualization scheme was implemented to better understand how the convolutional neural network made its judgement on single-voxel or multivoxel MRSI, and the convolutional neural network was embedded into a pipeline capable of producing whole-brain spectroscopic MRI volumes in real time. CONCLUSION: The fully automated method for assessment of spectral quality provides a valuable tool to support clinical MRSI or spectroscopic MRI studies for use in fields such as adaptive radiation therapy planning.

AutoVOI: real-time automatic prescription of volume-of-interest for single voxel spectroscopy.

PURPOSE: To develop a fast and automated volume-of-interest (VOI) prescription pipeline (AutoVOI) for single-voxel MRS that removes the need for manual VOI placement, allows flexible VOI planning in any brain region, and enables high inter- and intra-subject consistency of VOI prescription. METHODS: AutoVOI was designed to transfer pre-defined VOIs from an atlas to the 3D anatomical data of the subject during the scan. The AutoVOI pipeline was optimized for consistency in VOI placement (precision), enhanced coverage of the targeted tissue (accuracy), and fast computation speed. The tool was evaluated against manual VOI placement using existing T1 -weighted data sets and corresponding VOI prescriptions. Finally, it was implemented on 2 scanner platforms to acquire MRS data from clinically relevant VOIs that span the cerebrum, cerebellum, and the brainstem. RESULTS: The AutoVOI pipeline includes skull stripping, non-linear registration of the atlas to the subject's brain, and computation of the VOI coordinates and angulations using a minimum oriented bounding box algorithm. When compared against manual prescription, AutoVOI showed higher intra- and inter-subject spatial consistency, as quantified by generalized Dice coefficients (GDC), lower intra- and inter-subject variability in tissue composition (gray matter, white matter, and cerebrospinal fluid) and higher or equal accuracy, as quantified by GDC of prescribed VOI with targeted tissues. High quality spectra were obtained on Siemens and Philips 3T systems from 6 automatically prescribed VOIs by the tool. CONCLUSION: Robust automatic VOI prescription is feasible and can help facilitate clinical adoption of MRS by avoiding operator dependence of manual selection.

Global brain metabolic quantification with whole-head proton MRS at 3 T.

Total N-acetyl-aspartate + N-acetyl-aspartate-glutamate (NAA), total creatine (Cr) and total choline (Cho) proton MRS (1 H-MRS) signals are often used as surrogate markers in diffuse neurological pathologies, but spatial coverage of this methodology is limited to 1%-65% of the brain. Here we wish to demonstrate that non-localized, whole-head (WH) 1 H-MRS captures just the brain's contribution to the Cho and Cr signals, ignoring all other compartments. Towards this end, 27 young healthy adults (18 men, 9 women), 29.9 ± 8.5 years old, were recruited and underwent T1 -weighted MRI for tissue segmentation, non-localizing, approximately 3 min WH 1 H-MRS (TE /TR /TI  = 5/10/940 ms) and 30 min 1 H-MR spectroscopic imaging (MRSI) (TE /TR  = 35/2100 ms) in a 360 cm3 volume of interest (VOI) at the brain's center. The VOI absolute NAA, Cr and Cho concentrations, 7.7 ± 0.5, 5.5 ± 0.4 and 1.3 ± 0.2 mM, were all within 10% of the WH: 8.6 ± 1.1, 6.0 ± 1.0 and 1.3 ± 0.2 mM. The mean NAA/Cr and NAA/Cho ratios in the WH were only slightly higher than the "brain-only" VOI: 1.5 versus 1.4 (7%) and 6.6 versus 5.9 (11%); Cho/Cr were not different. The brain/WH volume ratio was 0.31 ± 0.03 (brain ≈ 30% of WH volume). Air-tissue susceptibility-driven local magnetic field changes going from the brain outwards showed sharp gradients of more than 100 Hz/cm (1 ppm/cm), explaining the skull's Cr and Cho signal losses through resonance shifts, line broadening and destructive interference. The similarity of non-localized WH and localized VOI NAA, Cr and Cho concentrations and their ratios suggests that their signals originate predominantly from the brain. Therefore, the fast, comprehensive WH-1 H-MRS method may facilitate quantification of these metabolites, which are common surrogate markers in neurological disorders.

Correlation between quantitative whole-body muscle magnetic resonance imaging and clinical muscle weakness in Pompe disease.

INTRODUCTION: Previous examination of whole-body muscle involvement in Pompe disease has been limited to physical examination and/or qualitative magnetic resonance imaging (MRI). In this study we assess the feasibility of quantitative proton-density fat-fraction (PDFF) whole-body MRI in late-onset Pompe disease (LOPD) and compare the results with manual muscle testing. METHODS: Seven LOPD patients and 11 disease-free controls underwent whole-body PDFF MRI. Quantitative MR muscle group assessments were compared with physical testing of muscle groups. RESULTS: The 95% upper limits of confidence intervals for muscle groups were 4.9-12.6% in controls and 6.8-76.4% in LOPD patients. LOPD patients showed severe and consistent tongue and axial muscle group involvement, with less marked involvement of peripheral musculature. MRI was more sensitive than physical examination for detection of abnormality in multiple muscle groups. CONCLUSION: This integrated, quantitative approach to muscle assessment provides more detailed data than physical examination and may have clinical utility for monitoring disease progression and treatment response.

Clinical proton MR spectroscopy in central nervous system disorders.

A large body of published work shows that proton (hydrogen 1 [(1)H]) magnetic resonance (MR) spectroscopy has evolved from a research tool into a clinical neuroimaging modality. Herein, the authors present a summary of brain disorders in which MR spectroscopy has an impact on patient management, together with a critical consideration of common data acquisition and processing procedures. The article documents the impact of (1)H MR spectroscopy in the clinical evaluation of disorders of the central nervous system. The clinical usefulness of (1)H MR spectroscopy has been established for brain neoplasms, neonatal and pediatric disorders (hypoxia-ischemia, inherited metabolic diseases, and traumatic brain injury), demyelinating disorders, and infectious brain lesions. The growing list of disorders for which (1)H MR spectroscopy may contribute to patient management extends to neurodegenerative diseases, epilepsy, and stroke. To facilitate expanded clinical acceptance and standardization of MR spectroscopy methodology, guidelines are provided for data acquisition and analysis, quality assessment, and interpretation. Finally, the authors offer recommendations to expedite the use of robust MR spectroscopy methodology in the clinical setting, including incorporation of technical advances on clinical units.

Automated whole-brain N-acetylaspartate proton MRS quantification.

Concentration of the neuronal marker, N-acetylaspartate (NAA), a quantitative metric for the health and density of neurons, is currently obtained by integration of the manually defined peak in whole-head proton ((1) H)-MRS. Our goal was to develop a full spectral modeling approach for the automatic estimation of the whole-brain NAA concentration (WBNAA) and to compare the performance of this approach with a manual frequency-range peak integration approach previously employed. MRI and whole-head (1) H-MRS from 18 healthy young adults were examined. Non-localized, whole-head (1) H-MRS obtained at 3 T yielded the NAA peak area through both manually defined frequency-range integration and the new, full spectral simulation. The NAA peak area was converted into an absolute amount with phantom replacement and normalized for brain volume (segmented from T1 -weighted MRI) to yield WBNAA. A paired-sample t test was used to compare the means of the WBNAA paradigms and a likelihood ratio test used to compare their coefficients of variation. While the between-subject WBNAA means were nearly identical (12.8 ± 2.5 mm for integration, 12.8 ± 1.4 mm for spectral modeling), the latter's standard deviation was significantly smaller (by ~50%, p = 0.026). The within-subject variability was 11.7% (±1.3 mm) for integration versus 7.0% (±0.8 mm) for spectral modeling, i.e., a 40% improvement. The (quantifiable) quality of the modeling approach was high, as reflected by Cramer-Rao lower bounds below 0.1% and vanishingly small (experimental - fitted) residuals. Modeling of the whole-head (1) H-MRS increases WBNAA quantification reliability by reducing its variability, its susceptibility to operator bias and baseline roll, and by providing quality-control feedback. Together, these enhance the usefulness of the technique for monitoring the diffuse progression and treatment response of neurological disorders.

Influence of running on femoroacetabular joint bone-to-bone distances.

There is limited data quantifying the influence of running on hip cartilage mechanics. The goal of this investigation was to quantify changes in hip joint bone-to-bone distance in response to a 3-mile treadmill run. We acquired magnetic resonance (MR) images of the dominant hip of eight young, asymptomatic runners (five males, three females) before and immediately after they ran 3 miles at a self-selected pace on a level treadmill. The femoral heads and acetabula were semiautomatically segmented from the pre- and post-exercise MR images to generate three-dimensional models of each participant's hip that were used to compute changes in the bone-to-bone distances incurred by the running exercise. We observed a significant 3% decrease in bone-to-bone distance from 3.47 ± 0.20 to 3.36 ± 0.22 mm between the femoral head and acetabulum after a 3-mile treadmill run (mean ± 95% confidence interval; p = 0.03). These findings provide new baseline data describing how running impacts the hip joint in young, asymptomatic runners.

Spin-echo magnetic resonance spectroscopic imaging at 7 T with frequency-modulated refocusing pulses.

Two approaches to high-resolution SENSE-encoded magnetic resonance spectroscopic imaging (MRSI) of the human brain at 7 Tesla (T) with whole-slice coverage are described. Both sequences use high-bandwidth radiofrequency pulses to reduce chemical shift displacement artifacts, SENSE-encoding to reduce scan time, and dual-band water and lipid suppression optimized for 7 T. Simultaneous B0 and transmit B1 mapping was also used for both sequences to optimize field homogeneity using high-order shimming and determine optimum radiofrequency transmit level, respectively. One sequence ("Hahn-MRSI") used reduced flip angle (90°) refocusing pulses for lower radiofrequency power deposition, while the other sequence used adiabatic fast passage refocusing pulses for improved sensitivity and reduced signal dependence on the transmit-B1 level. In four normal subjects, adiabatic fast passage-MRSI showed a signal-to-noise ratio improvement of 3.2±0.5 compared to Hahn-MRSI at the same spatial resolution, pulse repetition time, echo time, and SENSE-acceleration factor. An interleaved two-slice Hahn-MRSI sequence is also demonstrated to be experimentally feasible.

Improved initial value estimation for short echo time magnetic resonance spectroscopy spectral analysis using short T2 signal attenuation.

Robust spectral analysis of magnetic resonance spectroscopy data frequently uses a spectral model with prior metabolite signal information within a nonlinear least squares optimization algorithm. Starting values for the spectral model greatly influence the final results. Short echo time magnetic resonance spectroscopy contains broad signals that overlap with metabolite signals, complicating the estimation of starting values. We describe a method for more robust initial value estimation using a filter to attenuate short T(2) signal contributions (e.g., macromolecules or residual lipids). The method attenuates signals by truncating early points in the data set. Metabolite peak estimation is simplified by the removal of broad, short T(2) signals, and corrections for metabolite signal truncation are described. Short echo time simulated Monte Carlo data and in vivo data were used to validate the method. Areas for metabolite signals in the Monte Carlo data with singlet (N-acetylaspartate, creatine, choline) and singlet-like (myo-inositol) resonances were estimated within 10% of actual value for various metabolite line widths, signal-to-noise ratios, and underlying broad signal contributions. Initial value estimates of in vivo magnetic resonance spectroscopy data were within 14% of metabolite area ratios relative to the creatine peak fitted using established magnetic resonance spectroscopy spectral analysis software.

Bulbar muscle weakness and fatty lingual infiltration in glycogen storage disorder type IIIa.

Glycogen storage disorder type III (GSD III) is a rare autosomal recessive disorder resulting from a deficiency of glycogen debranching enzyme, critical in cytosolic glycogen degradation. GSD IIIa, the most common form of GSD III, primarily affects the liver, cardiac muscle, and skeletal muscle. Although skeletal muscle weakness occurs commonly in GSD IIIa, bulbar muscle involvement has not been previously reported. Here we present three GSD IIIa patients with clinical evidence of bulbar weakness based on instrumental assessment of lingual strength. Dysarthria and/or dysphagia, generally mild in severity, were evident in all three individuals. One patient also underwent correlative magnetic resonance imaging (MRI) which was remarkable for fatty infiltration at the base of the intrinsic tongue musculature, as well as abnormal expansion of the fibro-fatty lingual septum. Additionally, we provide supportive evidence of diffuse glycogen infiltration of the tongue at necropsy in a naturally occurring canine model of GSD IIIa. While further investigation in a larger group of patients with GSD III is needed to determine the incidence of bulbar muscle involvement in this condition and whether it occurs in GSD IIIb, clinical surveillance of lingual strength is recommended.

Comprehensive analysis of the Cramer-Rao bounds for magnetic resonance temperature change measurement in fat-water voxels using multi-echo imaging.

OBJECT: The aim of this paper is to characterize the noise propagation for MRI temperature change measurement with emphasis on finding the best echo time combinations that yield the lowest temperature noise. MATERIALS AND METHODS: A Cramer-Rao lower-bound (CRLB) calculation was used to estimate the temperature noise for a model of the MR signal in fat-water voxels. The temperature noise CRLB was then used to find a set of echo times that gave the lowest temperature change noise for a range of fat-water frequency differences, temperature changes, fat/water signal ratios, and T2* values. CRLB estimates were verified by Monte Carlo simulation and in phantoms using images acquired in a 1.5 T magnet. RESULTS: Results show that regions exist where the CRLB predicts minimal temperature variation as a function of the other variables. The results also indicate that the CRLB values calculated in this paper provide excellent guidance for predicting the variation of temperature measurements due to changes in the signal parameters. For three echo scans, the best noise characteristics are seen for TE values of 20.71, 23.71, and 26.71 ms. Results for five and seven echo scans are also presented in the text. CONCLUSION: The results present a comprehensive analysis of the effects of different scan parameters on temperature noise, potentially benefiting the selection of scan parameters for clinical MRI thermometry.