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.

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.

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.

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.

The YOUth cohort study: MRI protocol and test-retest reliability in adults.

The YOUth cohort study is a unique longitudinal study on brain development in the general population. As part of the YOUth study, 2000 children will be included at 8, 9 or 10 years of age and planned to return every three years during adolescence. Magnetic resonance imaging (MRI) brain scans are collected, including structural T1-weighted imaging, diffusion-weighted imaging (DWI), resting-state functional MRI and task-based functional MRI. Here, we provide a comprehensive report of the MR acquisition in YOUth Child & Adolescent including the test-retest reliability of brain measures derived from each type of scan. To measure test-retest reliability, 17 adults were scanned twice with a week between sessions using the full YOUth MRI protocol. Intraclass correlation coefficients were calculated to quantify reliability. Global brain measures derived from structural T1-weighted and DWI scans were reliable. Resting-state functional connectivity was moderately reliable, as well as functional brain measures for both the inhibition task (stop versus go) and the emotion task (face versus house). Our results complement previous studies by presenting reliability results of regional brain measures collected with different MRI modalities. YOUth facilitates data sharing and aims for reliable and high-quality data. Here we show that using the state-of-the art YOUth MRI protocol brain measures can be estimated reliably.

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

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.

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.

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.

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.

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.

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.