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.

T2-weighted turbo spin-echo magnetic resonance imaging of canine brain anatomy at 1.5T, 3T, and 7T field strengths.

Post-mortem T2 weighted images of canine heads were acquired at 1.5T, 3T, and 7T. This study aimed to (1) identify anatomical structures of the canine brain using an ultra-high-field magnetic resonance imaging (MRI) (7T) to help to facilitate their localization on high field MRI images (3T and 1.5T), where these structures may appear less well defined and delineated and (2) evaluate the visibility of canine brain anatomical structures on 1.5T, 3T, and 7T MRI images for optimizing clinical utility. Our hypothesis was that the provided subjective image quality comparison at different field strengths may offer a general baseline for canine brain anatomy and may help clinicians evaluate MRI options better. Six canine heads were examined with 1.5T, 3T, and 7T MRI scanners. T2-weighted images were acquired in three orthogonal planes at each field strength using a turbo spin-echo sequence. Fifty neuroanatomic structures were identified and evaluated on the 7T MR images; subsequently, those were found on the 3T and 45 out of the 50 structures were detected on the 1.5T imaging. The structures that were not able to be identified on the 1.5T imaging included the septum pellucidum, oculomotor nucleus, substantia nigra, claustrum, and thalamic nucleus griseus. Images acquired at 7T were subjective of higher spatial and contrast resolution. However, the ultra-high-field images were prone to artifacts at the interface between tissues of different magnetic properties. In conclusion, 3T MR imaging appears to be the best comprise for evaluating canine brain anatomy on MRI with fewer imaging artifacts.

Blood Flow Velocity Pulsatility and Arterial Diameter Pulsatility Measurements of the Intracranial Arteries Using 4D PC-MRI.

4D phase contrast magnetic resonance imaging (PC-MRI) allows for the visualization and quantification of the cerebral blood flow. A drawback of software that is used to quantify the cerebral blood flow is that it oftentimes assumes a static arterial luminal area over the cardiac cycle. Quantifying the lumen area pulsatility index (aPI), i.e. the change in lumen area due to an increase in distending pressure over the cardiac cycle, can provide insight in the stiffness of the arteries. Arterial stiffness has received increased attention as a predictor in the development of cerebrovascular disease. In this study, we introduce software that allows for measurement of the aPI as well as the blood flow velocity pulsatility index (vPI) from 4D PC-MRI. The internal carotid arteries of seven volunteers were imaged using 7 T MRI. The aPI and vPI measurements from 4D PC-MRI were validated against measurements from 2D PC-MRI at two levels of the internal carotid arteries (C3 and C7). The aPI and vPI computed from 4D PC-MRI were comparable to those measured from 2D PC-MRI (aPI: mean difference: 0.03 (limits of agreement: -0.14 - 0.23); vPI: 0.03 (-0.17-0.23)). The measured blood flow rate for the C3 and C7 segments was similar, indicating that our proposed software correctly captures the variation in arterial lumen area and blood flow velocity that exists along the distal end of the carotid artery. Our software may potentially aid in identifying changes in arterial stiffness of the intracranial arteries caused by pathological changes to the vessel wall.

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.

Explaining RF induced current patterns on implantable medical devices during MRI using the transfer matrix.

PURPOSE: In this work a simulation study is performed to gain insights in the patterns of induced radiofrequency (RF) currents for various implant-like structures at 1.5 T. The previously introduced transfer matrix (TM) is used to determine why certain current patterns have a tendency to naturally occur. This can benefit current safety assessment techniques and may enable the identification of critical exposure conditions. THEORY AND METHODS: The induced current on an elongated implant can be determined by multiplication of the incident electric field along the implant with its TM. The eigenmode spectrum of the TMs for various lengths and various types of implants are determined. The eigenvector with the highest eigenvalue describes the incident electric field pattern that induces the highest current which in turn will lead to highest heating. Subsequently, a statistical probability analysis is performed using a wide range of potential incident electric field distributions in a representative human subject model during a 1.5 T MR exam which are determined by means of electromagnetic FDTD simulations. These incident electric field distributions and the resulting induced current patterns are projected onto eigenvectors of the TM to determine which eigenmodes of the implant dominate the current patterns. RESULTS: The eigenvectors of the TM of bare and insulated wires resemble sinusoidal harmonics of a string fixed at both ends similar to the natural-current distribution on thin antennas(1). The currents on implants shorter than 20 cm are generally dominated by the first harmonic (similar to half a sine wave). This is firstly because for these implant lengths (relative to the RF wavelength), the first eigenvalue is more than three times bigger than the second showing the ability of an implant to accommodate one eigenmode better than another. Secondly, the incident electric fields have a high likelihood (≳95,7%) to project predominantly on this first eigenmode. CONCLUSION: The eigenmode spectrum of the TM of an implant provides insight into the expected shape of induced current distributions and worst-case exposure conditions. For short implants, the first eigenvector is dominant. In addition, realistic incident electric field distributions project more heavily on this eigenvector. Both effects together cause significant currents to always resemble the dominant eigenmode of the TM for short implants at 1.5 T.

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.

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%).

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.

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.

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.

Microbleeds colocalize with enlarged juxtacortical perivascular spaces in amnestic mild cognitive impairment and early Alzheimer's disease: A 7 Tesla MRI study.

MRI-visible perivascular spaces (PVS) in the semioval centre are associated with cerebral amyloid angiopathy (CAA), but it is unknown if PVS co-localize with MRI markers of CAA. To examine this, we assessed the topographical association between cortical cerebral microbleeds (CMBs) - as an indirect marker of CAA - and dilatation of juxtacortical perivascular spaces (jPVS) in 46 patients with amnestic mild cognitive impairment (aMCI) or early Alzheimer's disease (eAD). The degree of dilatation of jPVS <1 cm around each cortical CMBs was compared with a similar reference site (no CMB) in the contralateral hemisphere, using a 4-point scale. Also, jPVS dilatation was compared between patients with and without cortical CMBs. Eleven patients (24%) had cortical CMBs [total=35, median=1, range=1-14] of whom five had >1 cortical CMBs. The degree of jPVS dilatation was higher around CMBs than at the reference sites [Wilcoxon signed rank test, Z = 2.2, p = 0.03]. Patients with >1 cortical CMBs had a higher degree of jPVS dilation [median=2.2, IQR = 1.8-2.3] than patients without cortical CMBs [median=1.4, IQR = 1.0-1.8], p = 0.02. We found a topographical association between a high degree of jPVS dilatation and cortical CMBs, supporting a common underlying pathophysiology - most likely CAA.

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.

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.

MR-MOTUS: model-based non-rigid motion estimation for MR-guided radiotherapy using a reference image and minimal k-space data.

Time-resolved motion estimation from MRI data has received an increasing amount of interest due to the advent of the MR-Linac. The combination of an MRI scanner and a linear accelerator enables radiation plan adaptation based on internal organ motion estimated from MRI data. However, time-resolved estimation of this motion from MRI data still remains a challenge. In light of this application, we propose MR-MOTUS, a framework to estimate non-rigid 3D motion from minimal k-space data. MR-MOTUS consists of two main components: (1) a signal model that explicitly relates the k-space signal of a deforming object to non-rigid motion-fields and a reference image, and (2) model-based reconstructions of the non-rigid motion-fields directly from k-space data. Using an a priori available reference image and the fact that internal body motion exhibits a high level of spatial correlation, we represent the motion-fields in a low-dimensional space and reconstruct them from minimal k-space data that can be acquired very rapidly. The signal model is validated through numerical experiments with a digital 3D phantom and motion-fields are reconstructed from retrospectively undersampled in vivo head and abdomen data using various undersampling strategies. A comparison is made with state-of-the-art image registration performed on images reconstructed from the same undersampled data. Results show that MR-MOTUS reconstructs in vivo 3D rigid head motion from 474-fold retrospectively downsampled k-space data, and in vivo non-rigid 3D respiratory motion from 63-fold retrospectively undersampled k-space data. Preliminary results on prospectively undersampled data acquired with a 2D golden angle acquisition during free-breathing demonstrate the practical feasibility of the method.

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.

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.

Opening a new window on MR-based Electrical Properties Tomography with deep learning.

In the radiofrequency (RF) range, the electrical properties of tissues (EPs: conductivity and permittivity) are modulated by the ionic and water content, which change for pathological conditions. Information on tissues EPs can be used e.g. in oncology as a biomarker. The inability of MR-Electrical Properties Tomography techniques (MR-EPT) to accurately reconstruct tissue EPs by relating MR measurements of the transmit RF field to the EPs limits their clinical applicability. Instead of employing electromagnetic models posing strict requirements on the measured MRI quantities, we propose a data driven approach where the electrical properties reconstruction problem can be casted as a supervised deep learning task (DL-EPT). DL-EPT reconstructions for simulations and MR measurements at 3 Tesla on phantoms and human brains using a conditional generative adversarial network demonstrate high quality EPs reconstructions and greatly improved precision compared to conventional MR-EPT. The supervised learning approach leverages the strength of electromagnetic simulations, allowing circumvention of inaccessible MR electromagnetic quantities. Since DL-EPT is more noise-robust than MR-EPT, the requirements for MR acquisitions can be relaxed. This could be a major step forward to turn electrical properties tomography into a reliable biomarker where pathological conditions can be revealed and characterized by abnormalities in tissue electrical properties.

Contradiction between amide-CEST signal and pH in breast cancer explained with metabolic MRI.

PURPOSE: Metabolic MRI is a noninvasive technique that can give new insights into understanding cancer metabolism and finding biomarkers to evaluate or monitor treatment plans. Using this technique, a previous study has shown an increase in pH during neoadjuvant chemotherapy (NAC) treatment, while recent observation in a different study showed a reduced amide proton transfer (APT) signal during NAC treatment (negative relation). These findings are counterintuitive, given the known intrinsic positive relation of APT signal to pH. METHODS: In this study we combined APT MRI and 31 P-MRSI measurements to unravel the relation between the APT signal and pH in breast cancer. Twenty-two breast cancer patients were scanned with a 7 T MRI before and after the first cycle of NAC treatment. pH was determined by the chemical shift of inorganic phosphate (Pi). RESULTS: While APT signals have a positive relation to pH and amide content, we observed a direct negative linear correlation between APT signals and pH in breast tumors in vivo. CONCLUSIONS: As differentiation of cancer stages was confirmed by observation of a linear correlation between cell proliferation marker PE/Pi (phosphoethanolamine over inorganic phosphate) and pH in the tumor, our data demonstrates that the concentration of mobile proteins likely supersedes the contribution of the exchange rate to the APT signal.

Early detection of changes in phospholipid metabolism during neoadjuvant chemotherapy in breast cancer patients using phosphorus magnetic resonance spectroscopy at 7T.

The purpose of this work was to investigate whether noninvasive early detection (after the first cycle) of response to neoadjuvant chemotherapy (NAC) in breast cancer patients was possible. 31 P-MRSI at 7 T was used to determine different phosphor metabolites ratios and correlate this to pathological response. 31 P-MRSI was performed in 12 breast cancer patients treated with NAC. 31 P spectra were fitted and aligned to the frequency of phosphoethanolamine (PE). Metabolic signal ratios for phosphomonoesters/phosphodiesters (PME/PDE), phosphocholine/glycerophosphatidylcholine (PC/GPtC), phosphoethanolamine/glycerophosphoethanolamine (PE/GPE) and phosphomonoesters/in-organic phosphate (PME/Pi) were determined from spectral fitting of the individual spectra and the summed spectra before and after the first cycle of NAC. Metabolic ratios were subsequently related to pathological response. Additionally, the correlation between the measured metabolic ratios and Ki-67 levels was determined using linear regression. Four patients had a pathological complete response after treatment, five patients a partial pathological response, and three patients did not respond to NAC. In the summed spectrum after the first cycle of NAC, PME/Pi and PME/PDE decreased by 18 and 13%, respectively. A subtle difference among the different response groups was observed in PME/PDE, where the nonresponders showed an increase and the partial and complete responders a decrease (P = 0.32). No significant changes in metabolic ratios were found. However, a significant association between PE/Pi and the Ki-67 index was found (P = 0.03). We demonstrated that it is possible to detect subtle changes in 31 P metabolites with a 7 T MR system after the first cycle of NAC treatment in breast cancer patients. Nonresponders showed different changes in metabolic ratios compared with partial and complete responders, in particular for PME/PDE; however, more patients need to be included to investigate its clinical value.