Tier-based formalism for safety assessment of custom-built radio-frequency transmit coils.

The purpose of this work is to propose a tier-based formalism for safety assessment of custom-built radio-frequency (RF) coils that balances validation effort with the effort put in determinating the safety factor. The formalism has three tier levels. Higher tiers require increased effort when validating electromagnetic simulation results but allow for less conservative safety factors. In addition, we propose a new method to calculate modeling uncertainty between simulations and measurements and a new method to propagate uncertainties in the simulation into a safety factor that minimizes the risk of underestimating the peak specific absorption rate (SAR). The new safety assessment procedure was completed for all tier levels for an eight-channel dipole array for prostate imaging at 7 T and an eight-channel dipole array for head imaging at 10.5 T, using data from two different research sites. For the 7 T body array, the validation procedure resulted in a modeling uncertainty of 77% between measured and simulated local SAR distributions. For a situation where RF shimming is performed on the prostate, average power limits of 2.4 and 4.5 W/channel were found for tiers 2 and 3, respectively. When the worst-case peak SAR among all phase settings was calculated, power limits of 1.4 and 2.7 W/channel were found for tiers 2 and 3, respectively. For the 10.5 T head array, a modeling uncertainty of 21% was found based on B1 + mapping. For the tier 2 validation, a power limit of 2.6 W/channel was calculated. The demonstrated tier system provides a strategy for evaluating modeling inaccuracy, allowing for the rapid translation of novel coil designs with conservative safety factors and the implementation of less conservative safety factors for frequently used coil arrays at the expense of increased validation effort.

Real-time assessment of potential peak local specific absorption rate value without phase monitoring: Trigonometric maximization method for worst-case local specific absorption rate determination.

PURPOSE: Multi-transmit MRI systems are typically equipped with dedicated hardware to sample the reflected/lost power in the transmit channels. After extensive calibration, the amplitude and phase of the signal at the feed of each array element can be accurately determined. However, determining the phase is more difficult and monitoring errors can lead to a hazardous peak local specific absorption rate (pSAR10g ) underestimation. For this purpose, methods were published for online maximum potential pSAR10g estimation without relying on phase monitoring, but these methods produce considerable overestimation. We present a trigonometric maximization method to determine the actual worst-case pSAR10g without any overestimation. THEORY AND METHOD: The proposed method takes advantage of the sinusoidal relation between the SAR10g in each voxel and the phases of input signals, to return the maximum achievable SAR10g in a few iterations. The method is applied to determine the worst-case pSAR10g for three multi-transmit array configurations at 7T: (1) body array with eight fractionated dipoles; (2) head array with eight fractionated dipoles; (3) head array with eight rectangular loops. The obtained worst-case pSAR10g values are compared with the pSAR10g values determined with a commonly used method and with a more efficient method based on reference-phases. RESULTS: For each voxel, the maximum achievable SAR10g is determined in less than 0.1 ms. Compared to the reference-phases-based method, the proposed method reduces the mean overestimation of the actual pSAR10g up to 52%, while never underestimating the true pSAR10g . CONCLUSION: The proposed method can widely improve the performance of parallel transmission MRI systems without phase monitoring.

SAR and temperature distributions in a database of realistic human models for 7 T cardiac imaging.

PURPOSE: To investigate inter-subject variability of B1+ , SAR and temperature rise in a database of human models using a local transmit array for 7 T cardiac imaging. METHODS: Dixon images were acquired of 14 subjects and segmented in dielectric models with an eight-channel local transmit array positioned around the torso for cardiac imaging. EM simulations were done to calculate SAR distributions. Based on the SAR distributions, temperature simulations were performed for exposure times of 6 min and 30 min. Peak local SAR and temperature rise levels were calculated for different RF shim settings. A statistical analysis of the resulting peak local SAR and temperature rise levels was performed to arrive at safe power limits. RESULTS: For RF shim vectors with random phase and uniformly distributed power, a safe average power limit of 35.7 W was determined (first level controlled mode). When RF amplitude and phase shimming was performed on the heart, a safe average power limit of 35.0 W was found. According to Pennes' model, our numerical study suggests a very low probability of exceeding the absolute local temperature limit of 40 °C for a total exposure time of 6 min and a peak local SAR of 20 W/kg. For a 30 min exposure time at 20 W/kg, it was shown that the absolute temperature limit can be exceeded in the case where perfusion does not change with temperature. CONCLUSION: Safe power constraints were found for 7 T cardiac imaging with an eight-channel local transmit array, while considering the inter-subject variability of B1+ , SAR and temperature rise.

Volume increase in the dentate gyrus after electroconvulsive therapy in depressed patients as measured with 7T.

Electroconvulsive therapy (ECT) is the most effective treatment for depression, yet its working mechanism remains unclear. In the animal analog of ECT, neurogenesis in the dentate gyrus (DG) of the hippocampus is observed. In humans, volume increase of the hippocampus has been reported, but accurately measuring the volume of subfields is limited with common MRI protocols. If the volume increase of the hippocampus in humans is attributable to neurogenesis, it is expected to be exclusively present in the DG, whereas other processes (angiogenesis, synaptogenesis) also affect other subfields. Therefore, we acquired an optimized MRI scan at 7-tesla field strength allowing sensitive investigation of hippocampal subfields. A further increase in sensitivity of the within-subjects measurements is gained by automatic placement of the field of view. Patients receive two MRI scans: at baseline and after ten bilateral ECT sessions (corresponding to a 5-week interval). Matched controls are also scanned twice, with a similar 5-week interval. A total of 31 participants (23 patients, 8 controls) completed the study. A large and significant increase in DG volume was observed after ECT (M = 75.44 mm3, std error = 9.65, p < 0.001), while other hippocampal subfields were unaffected. We note that possible type II errors may be present due to the small sample size. In controls no changes in volume were found. Furthermore, an increase in DG volume was related to a decrease in depression scores, and baseline DG volume predicted clinical response. These findings suggest that the volume change of the DG is related to the antidepressant properties of ECT, and may reflect neurogenesis.

Validating faster DENSE measurements of cardiac-induced brain tissue expansion as a potential tool for investigating cerebral microvascular pulsations.

Displacement Encoding with Stimulated Echoes (DENSE) has recently shown potential for measuring cardiac-induced cerebral volumetric strain in the human brain. As such, it may provide a powerful tool for investigating the cerebral small vessels. However, further development and validation are necessary. This study aims, first, to validate a retrospectively-gated implementation of the DENSE method for assessing brain tissue pulsations as a physiological marker, and second, to use the acquired measurements to explore intracranial volume dynamics. We acquired repeated measurements of cerebral volumetric strain in 8 healthy subjects, and internally validated these measurements by comparing them to spinal CSF stroke volumes obtained in the same scan session. Peak volumetric strain was found to be highly repeatable between scan sessions. First/second measured peak volumetric strains were: (6.4 â€‹± â€‹1.7)x10-4/(6.7 â€‹± â€‹1.6)x10-4 for whole brain, (9.5 â€‹± â€‹2.5)x10-4/(9.6 â€‹± â€‹2.4)x10-4 for grey matter, and (4.4 â€‹± â€‹1.7)x10-4/(4.1 â€‹± â€‹0.8)x10-4 for white matter. Grey matter showed significantly higher peak strain (p â€‹< â€‹0.001) and earlier time-to-peak strain (p â€‹< â€‹0.02) than white matter. An approximately linear relationship was found between CSF and brain tissue volume pulsations over the cardiac cycle (mean slope and R2 of 0.88 â€‹± â€‹0.23 and 0.89 â€‹± â€‹0.07, respectively). The close similarity between CSF and brain tissue volume pulsations implies limited contributions from large intracranial vessel pulsations, providing further evidence for venous compression as an additional mechanism for maintaining stable intracranial pressures over the cardiac cycle. Cerebral pulsatility showed consistent inter-subject peak values in healthy subjects, and was strongly correlated to CSF stroke volumes. These results strengthen the potential of brain tissue volumetric strain as a means for investigating the intracranial dynamics of the ageing brain in normal or diseased states.

MRI-based transfer function determination through the transfer matrix by jointly fitting the incident and scattered B1+ field.

PURPOSE: A purely experimental method for MRI-based transfer function (TF) determination is presented. A TF characterizes the potential for radiofrequency heating of a linear implant by relating the incident tangential electric field to a scattered electric field at its tip. We utilize the previously introduced transfer matrix (TM) to determine transfer functions solely from the MR measurable quantities, that is, the B1+ and transceive phase distributions. This technique can extend the current practice of phantom-based TF assessment with dedicated experimental setup toward MR-based methods that have the potential to assess the TF in more realistic situations. THEORY AND METHODS: An analytical description of the B1+ magnitude and transceive phase distribution around a wire-like implant was derived based on the TM. In this model, the background field is described using a superposition of spherical and cylindrical harmonics while the transfer matrix is parameterized using a previously introduced attenuated wave model. This analytical description can be used to estimate the transfer matrix and transfer function based on the measured B1+ distribution. RESULTS: The TF was successfully determined for 2 mock-up implants: a 20-cm bare copper wire and a 20-cm insulated copper wire with 10 mm of insulation stripped at both endings in respectively 4 and 3 different trajectories. The measured TFs show a strong correlation with a reference determined from simulations and between the separate experiments with correlation coefficients above 0.96 between all TFs. Compared to the simulated TF, the maximum deviation in the estimated tip field is 9.4% and 12.2% for the bare and insulated wire, respectively. CONCLUSIONS: A method has been developed to measure the TF of medical implants using MRI experiments. Jointly fitting the incident and scattered B1+ distributions with an analytical description based on the transfer matrix enables accurate determination of the TF of 2 test implants. The presented method no longer needs input from simulated data and can therefore, in principle, be used to measure TF's in test animals or corpses.

Conditional safety margins for less conservative peak local SAR assessment: A probabilistic approach.

PURPOSE: The introduction of a linear safety factor to address peak local specific absorption rate (pSAR10g ) uncertainties (eg, intersubject variation, modeling inaccuracies) bears one considerable drawback: It often results in over-conservative scanning constraints. We present a more efficient approach to define a variable safety margin based on the conditional probability density function of the effectively obtained pSAR10g value, given the estimated pSAR10g value. METHODS: The conditional probability density function can be estimated from previously simulated data. A representative set of true and estimated pSAR10g samples was generated by means of our database of 23 subject-specific models with an 8-fractionated dipole array for prostate imaging at 7 T. The conditional probability density function was calculated for each possible estimated pSAR10g value and used to determine the corresponding safety margin with an arbitrary low probability of underestimation. This approach was applied to five state-of-the-art local SAR estimation methods, namely: (1) using just the generic body model "Duke"; (2) using our model library to assess the maximum pSAR10g value over all models; (3) using the most representative "local SAR model"; (4) using the five most representative local SAR models; and (5) using a recently developed deep learning-based method. RESULTS: Compared with the more conventional safety factor, the conditional safety-margin approach results in lower (up to 30%) mean overestimation for all investigated local SAR estimation methods. CONCLUSION: The proposed probabilistic approach for pSAR10g correction allows more accurate local SAR assessment with much lower overestimation, while a predefined level of underestimation is accepted (eg, 0.1%).

Introduction of the snake antenna array: Geometry optimization of a sinusoidal dipole antenna for 10.5T body imaging with lower peak SAR.

PURPOSE: To improve imaging performance for body MRI with a local transmit array at 10.5T, the geometry of a dipole antenna was optimized to achieve lower peak specific absorption rate (SAR) levels and a more uniform transmit profile. METHODS: Electromagnetic simulations on a phantom were used to evaluate the SAR and B 1 + -performance of different dipole antenna geometries. The best performing antenna (the snake antenna) was simulated on human models in a 12-channel array configuration for safety assessment and for comparison to a previous antenna design. This 12-channel array was constructed after which electromagnetic simulations were validated by B 1 + -maps and temperature measurements. After obtaining approval by the Food and Drug Administration to scan with the snake antenna array, in vivo imaging was performed on 2 volunteers. RESULTS: Simulation results on a phantom indicate a lower SAR and a higher transmit efficiency for the snake antenna compared to the fractionated dipole array. Similar results are found on a human body model: when comparing the trade-off between uniformity and peak SAR, the snake antenna performs better for all imaging targets. Simulations and measurements are in good agreement. Preliminary imaging result were acquired in 2 volunteers with the 12-channel snake antenna array. CONCLUSION: By optimizing the geometry of a dipole antenna, peak SAR levels were lowered while achieving a more uniform transmit field as demonstrated in simulations on a phantom and a human body model. The array was constructed, validated, and successfully used to image 2 individuals at 10.5T.

Correcting time-intensity curves in dynamic contrast-enhanced breast MRI for inhomogeneous excitation fields at 7T.

PURPOSE: Inhomogeneous excitation at ultrahigh field strengths (7T and above) compromises the reliability of quantified dynamic contrast-enhanced breast MRI. This can hamper the introduction of ultrahigh field MRI into the clinic. Compensation for this non-uniformity effect can consist of both hardware improvements and post-acquisition corrections. This paper investigated the correctable radiofrequency transmit ( B1+ ) range post-acquisition in both simulations and patient data for 7T MRI. METHODS: Simulations were conducted to determine the minimum B1+ level at which corrections were still beneficial because of noise amplification. Two correction strategies leading to differences in noise amplification were tested. The effect of the corrections on a 7T patient data set (N = 38) with a wide range of B1+ levels was investigated in terms of time-intensity curve types as well as washin, washout and peak enhancement values. RESULTS: In simulations assuming a common amount of T1 saturation, the lowest B1+ level at which the SNR of the corrected images was at least that of the original precontrast image was 43% of the nominal angle. After correction, time-intensity curve types changed in 24% of included patients, and the distribution of curve types corresponded better to the distribution found in literature. Additionally, the overlap between the distributions of washin, washout, and peak enhancement values for grade 1 and grade 2 tumors was slightly reduced. CONCLUSION: Although the correctable range varies with the amount of T1 saturation, post-acquisition correction for inhomogeneous excitation was feasible down to B1+ levels of 43% of the nominal angle in vivo.

Accelerating implant RF safety assessment using a low-rank inverse update method.

PURPOSE: Patients who have medical metallic implants, e.g. orthopaedic implants and pacemakers, often cannot undergo an MRI exam. One of the largest risks is tissue heating due to the radio frequency (RF) fields. The RF safety assessment of implants is computationally demanding. This is due to the large dimensions of the transmit coil compared to the very detailed geometry of an implant. METHODS: In this work, we explore a faster computational method for the RF safety assessment of implants that exploits the small geometry. The method requires the RF field without an implant as a basis and calculates the perturbation that the implant induces. The inputs for this method are the incident fields and a library matrix that contains the RF field response of every edge an implant can occupy. Through a low-rank inverse update, using the Sherman-Woodbury-Morrison matrix identity, the EM response of arbitrary implants can be computed within seconds. We compare the solution from full-wave simulations with the results from the presented method, for two implant geometries. RESULTS: From the comparison, we found that the resulting electric and magnetic fields are numerically equivalent (maximum error of 1.35%). However, the computation was between 171 to 2478 times faster than the corresponding GPU accelerated full-wave simulation. CONCLUSIONS: The presented method enables for rapid and efficient evaluation of the RF fields near implants and might enable situation-specific scanning conditions.

High-resolution in vivo MR-STAT using a matrix-free and parallelized reconstruction algorithm.

MR-STAT is a recently proposed framework that allows the reconstruction of multiple quantitative parameter maps from a single short scan by performing spatial localisation and parameter estimation on the time-domain data simultaneously, without relying on the fast Fourier transform (FFT). To do this at high resolution, specialized algorithms are required to solve the underlying large-scale nonlinear optimisation problem. We propose a matrix-free and parallelized inexact Gauss-Newton based reconstruction algorithm for this purpose. The proposed algorithm is implemented on a high-performance computing cluster and is demonstrated to be able to generate high-resolution (1 mm × 1 mm in-plane resolution) quantitative parameter maps in simulation, phantom, and in vivo brain experiments. Reconstructed T1 and T2 values for the gel phantoms are in agreement with results from gold standard measurements and, for the in vivo experiments, the quantitative values show good agreement with literature values. In all experiments, short pulse sequences with robust Cartesian sampling are used, for which MR fingerprinting reconstructions are shown to fail.

Pushing functional MRI spatial and temporal resolution further: High-density receive arrays combined with shot-selective 2D CAIPIRINHA for 3D echo-planar imaging at 7 T.

To be able to examine dynamic and detailed brain functions, the spatial and temporal resolution of 7 T MRI needs to improve. In this study, it was investigated whether submillimeter multishot 3D EPI fMRI scans, acquired with high-density receive arrays, can benefit from a 2D CAIPIRINHA sampling pattern, in terms of noise amplification (g-factor), temporal SNR and fMRI sensitivity. High-density receive arrays were combined with a shot-selective 2D CAIPIRINHA implementation for multishot 3D EPI sequences at 7 T. In this implementation, in contrast to conventional inclusion of extra kz gradient blips, specific EPI shots are left out to create a CAIPIRINHA shift and reduction of scan time. First, the implementation of the CAIPIRINHA sequence was evaluated with a standard receive setup by acquiring submillimeter whole brain T2 *-weighted anatomy images. Second, the CAIPIRINHA sequence was combined with high-density receive arrays to push the temporal resolution of submillimeter 3D EPI fMRI scans of the visual cortex. Results show that the shot-selective 2D CAIPIRINHA sequence enables a reduction in scan time for 0.5 mm isotropic 3D EPI T2 *-weighted anatomy scans by a factor of 4 compared with earlier reports. The use of the 2D CAIPIRINHA implementation in combination with high-density receive arrays, enhances the image quality of submillimeter 3D EPI scans of the visual cortex at high acceleration as compared to conventional SENSE. Both the g-factor and temporal SNR improved, resulting in a method that is more sensitive to the fMRI signal. Using this method, it is possible to acquire submillimeter single volume 3D EPI scans of the visual cortex in a subsecond timeframe. Overall, high-density receive arrays in combination with shot-selective 2D CAIPIRINHA for 3D EPI scans prove to be valuable for reducing the scan time of submillimeter MRI acquisitions.

Metabolite cycled liver 1 H MRS on a 7 T parallel transmit system.

INTRODUCTION: Single-voxel 1 H MRS in body applications often suffers from respiratory and other motion induced phase and frequency shifts, which lead to incoherent averaging and hence to suboptimal results. METHODS: Here we show the application of metabolite cycling (MC) for liver STEAM-localized 1 H MRS on a 7 T parallel transmit system, using eight transmit-receive fractionated dipole antennas with 16 additional, integrated receive loops. MC-STEAM measurements were made in six healthy, lean subjects and compared with STEAM measurements using VAPOR water suppression. Measurements were performed during free breathing and during synchronized breathing, for which the subjects did breathe in between the MRS acquisitions. Both intra-session repeatability and inter-session reproducibility of liver lipid quantification with MC-STEAM and VAPOR-STEAM were determined. RESULTS: The preserved water signal in MC-STEAM allowed for robust phase and frequency correction of individual acquisitions before averaging, which resulted in in vivo liver spectra that were of equal quality when measurements were made with free breathing or synchronized breathing. Intra-session repeatability and inter-session reproducibility of liver lipid quantification were better for MC-STEAM than for VAPOR-STEAM. This may also be explained by the more robust phase and frequency correction of the individual MC-STEAM acquisitions as compared with the VAPOR-STEAM acquisitions, for which the low-signal-to-noise ratio lipid signals had to be used for the corrections. CONCLUSION: Non-water-suppressed MC-STEAM on a 7 T system with parallel transmit is a promising approach for 1 H MRS applications in the body that are affected by motion, such as in the liver, and yields better repeatability and reproducibility compared with water-suppressed measurements.

Higher Pulsatility in Cerebral Perforating Arteries in Patients With Small Vessel Disease Related Stroke, a 7T MRI Study.

Background and Purpose- Cerebral small vessel disease (SVD) is a major cause of stroke and dementia, but underlying disease mechanisms are still largely unknown, partly because of the difficulty in assessing small vessel function in vivo. We developed a method to measure blood flow velocity pulsatility in perforating arteries in the basal ganglia and semioval center. We aimed to determine whether this novel method could detect functional abnormalities at the level of the small vessels in patients with stroke attributable to SVD. Methods- We investigated 10 patients with lacunar infarction (mean age 61 years, 80% men), 11 patients with deep intracerebral hemorrhage (ICH) considered to be caused by SVD (ICH, mean age 58 years, 82% men) and 18 healthy controls that were age- and sex-matched. We performed 2-dimensional phase contrast magnetic resonance imaging at 7 T to measure time-resolved blood flow velocity in cerebral perforating arteries of the semioval center and the basal ganglia. We compared the number of detected arteries, pulsatility index and mean velocity between the patient groups and controls. Results- In the basal ganglia, the number of detected perforators was lower in lacunar infarction (26±9, P=0.01) and deep ICH patients (28±6, P=0.02) than in controls (35±7). The pulsatility index in the basal ganglia was higher in lacunar infarction (1.07±0.13, P=0.03), and deep ICH patients (1.02±0.11, P=0.11), than in controls (0.94±0.10). Observations in the semioval center were similar. Number of detected perforators was lower in lacunar infarction (32±18, P=0.06), and deep ICH patients (28±18, P=0.02), than in controls (45±16). The pulsatility index was higher in lacunar infarction (1.18±0.15, P=0.02), and deep ICH patients (1.17±0.14, P=0.045) than in controls (1.08±0.07). No velocity differences were detected. Conclusions- This exploratory study shows that SVD can be expressed in terms of functional measures, such as pulsatility index, which are derived directly from the small vessels themselves. Future studies may use this technique to further unravel the mechanisms underlying SVD.

Potential acceleration performance of a 256-channel whole-brain receive array at 7 T.

PURPOSE: Assess the potential gain in acceleration performance of a 256-channel versus 32-channel receive coil array at 7 T in combination with a 2D CAIPIRINHA sequence for 3D data sets. METHODS: A 256-channel receive setup was simulated by placing 2 small 16-channel high-density receive arrays at 2 × 8 different locations on the head of healthy participants. Multiple consecutive measurements were performed and coil sensitivity maps were combined to form a complete 256-channel data set. This setup was compared with a standard 32-channel head coil, in terms of SNR, noise correlation, and acceleration performance (g-factor). RESULTS: In the periphery of the brain, the receive SNR was on average a factor 1.5 higher (ranging up to a factor 2.7 higher) than the 32-channel coil; in the center of the brain the SNR was comparable or lower, depending on the size of the region of interest, with a factor 1.0 on average (ranging from 0.7 up to a factor of 1.6). The average noise correlation between coil elements was 3% for the 256-channel coil, and 5% for the 32-channel coil. At acceptable g-factors (< 2), the achievable acceleration factor using SENSE and 2D CAIPIRINHA was 24 and 28, respectively, versus 9 and 12 for the 32-channel coil. CONCLUSION: The receive performance of the simulated 256 channel array was better than the 32-channel reference. Combined with 2D CAIPIRINHA, a peak acceleration factor of 28 was assessed, showing great potential for high-density receive arrays.

Intersubject specific absorption rate variability analysis through construction of 23 realistic body models for prostate imaging at 7T.

PURPOSE: For ultrahigh field (UHF) MRI, the expected local specific absorption rate (SAR) distribution is usually calculated by numerical simulations using a limited number of generic body models and adding a safety margin to take into account intersubject variability. Assessment of this variability with a large model database would be desirable. In this study, a procedure to create such a database with accurate subject-specific models is presented. Using 23 models, intersubject variability is investigated for prostate imaging at 7T with an 8-channel fractionated dipole antenna array with 16 receive loops. METHOD: From Dixon images of a volunteer acquired at 1.5T with a mockup array in place, an accurate dielectric model is built. Following this procedure, 23 subject-specific models for local SAR assessment at 7T were created enabling an extensive analysis of the intersubject B1+ and peak local SAR variability. RESULTS: For the investigated setup, the maximum possible peak local SAR ranges from 2.6 to 4.6 W/kg for 8 × 1 W input power. The expected peak local SAR values represent a Gaussian distribution (µ/σ=2.29/0.29 W/kg) with realistic prostate-shimmed phase settings and a gamma distribution Γ(24,0.09) with multidimensional radiofrequency pulses. Prostate-shimmed phase settings are similar for all models. Using 1 generic phase setting, average B1+ reduction is 7%. Using only 1 model, the required safety margin for intersubject variability is 1.6 to 1.8. CONCLUSION: The presented procedure allows for the creation of a customized model database. The results provide valuable insights into B1+ and local SAR variability. Recommended power thresholds per channel are 3.1 W with phase shimming on prostate or 2.6 W for multidimensional pulses.

Reducing distortions in echo-planar breast imaging at ultrahigh field with high-resolution off-resonance maps.

PURPOSE: DWI is a promising modality in breast MRI, but its clinical acceptance is slow. Analysis of DWI is hampered by geometric distortion artifacts, which are caused by off-resonant spins in combination with the low phase-encoding bandwidth of the EPI sequence used. Existing correction methods assume smooth off-resonance fields, which we show to be invalid in the human breast, where high discontinuities arise at tissue interfaces. METHODS: We developed a distortion correction method that incorporates high-resolution off-resonance maps to better solve for severe distortions at tissue interfaces. The method was evaluated quantitatively both ex vivo in a porcine tissue phantom and in vivo in 5 healthy volunteers. The added value of high-resolution off-resonance maps was tested using a Wilcoxon signed rank test comparing the quantitative results obtained with a low-resolution off-resonance map with those obtained with a high-resolution map. RESULTS: Distortion correction using low-resolution off-resonance maps corrected most of the distortions, as expected. Still, all quantitative comparison metrics showed increased conformity between the corrected EPI images and a high-bandwidth reference scan for both the ex vivo and in vivo experiments. All metrics showed a significant improvement when a high-resolution off-resonance map was used (P < 0.05), in particular at tissue boundaries. CONCLUSION: The use of off-resonance maps of a resolution higher than EPI scans significantly improves upon existing distortion correction techniques, specifically by superior correction at glandular tissue boundaries.

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.

Quantifying cardiac-induced brain tissue expansion using DENSE.

Brain tissue undergoes viscoelastic deformation and volumetric strain as it expands over the cardiac cycle due to blood volume changes within the underlying microvasculature. Volumetric strain measurements may therefore provide insights into small vessel function and tissue viscoelastic properties. Displacement encoding via stimulated echoes (DENSE) is an MRI technique that can quantify the submillimetre displacements associated with brain tissue motion. Despite previous studies reporting brain tissue displacements using DENSE and other MRI techniques, a complete picture of brain tissue volumetric strain over the cardiac cycle has not yet been obtained. To address this need we implemented 3D cine-DENSE at 7 T and 3 T to investigate the feasibility of measuring cardiac-induced volumetric strain as a marker for small vessel blood volume changes. Volumetric strain over the entire cardiac cycle was computed for the whole brain and for grey and white matter tissue separately in six healthy human subjects. Signal-to-noise ratio (SNR) measurements were used to determine the voxel-wise volumetric strain noise. Mean peak whole brain volumetric strain at 7 T (mean ± SD) was (4.5 ± 1.0) × 10-4 (corresponding to a volume expansion of 0.48 ± 0.1 mL), which is in agreement with literature values for cerebrospinal fluid that is displaced into the spinal canal to maintain a stable intracranial pressure. The peak volumetric strain ratio of grey to white matter was 4.4 ± 2.8, reflecting blood volume and tissue stiffness differences between these tissue types. The mean peak volumetric strains of grey and white matter tissue were found to be significantly different (p < 0.001). The mean SNR at 7 T and 3 T of the DENSE measurements was 22.0 ± 7.3 and 7.0 ± 2.8 respectively, which currently limits a voxel-wise strain analysis at both field strengths. We demonstrate that tissue specific quantification of volumetric strain is feasible with DENSE. This metric holds potential for studying blood volume pulsations in the ageing brain in healthy and diseased states.

SNR optimized 31 P functional MRS to detect mitochondrial and extracellular pH change during visual stimulation.

Energy metabolism of the human visual cortex was investigated by performing 31 P functional MRS. INTRODUCTION: The human brain is known to be the main glucose demanding organ of the human body and neuronal activity can increase this energy demand. In this study we investigate whether alterations in pH during activation of the brain can be observed with MRS, focusing on the mitochondrial inorganic phosphate (Pi) pool as potential marker of energy demand. METHODS: Six participants were scanned with 16 consecutive 31 P-MRSI scans, which were divided in 4 blocks of 8:36 minutes of either rest or visual stimulation. Since the signals from the mitochondrial compartments of Pi are low, multiple approaches to achieve high SNR 31 P measurements were combined. This included: a close fitting 31 P RF coil, a 7 T-field strength, Ernst angle acquisitions and a stimulus with a large visual angle allowing large spectroscopy volumes containing activated tissue. RESULTS: The targeted resonance downfield of the main Pi peak could be distinguished, indicating the high SNR of the 31 P spectra. The peak downfield of the main Pi peak is believed to be connected to mitochondrial performance. In addition, a BOLD effect in the PCr signal was observed as a signal increase of 2-3% during visual stimulation as compared to rest. When averaging data over multiple volunteers, a small subtle shift of about 0.1 ppm of the downfield Pi peak towards the main Pi peak could be observed in the first 4 minutes of visual stimulation, but no longer in the 4 to 8 minute scan window. Indications of a subtle shift during visual stimulation were found, but this effect remains small and should be further validated. CONCLUSION: Overall, the downfield peak of Pi could be observed, revealing opportunities and considerations to measure specific acidity (pH) effects in the human visual cortex.