Low-intensity transcranial focused ultrasound stimulation (TFUS) holds great promise as a highly focal technique for transcranial stimulation even for deep brain areas. Yet, knowledge about the ...safety of this novel technique is still limited.
To systematically review safety related aspects of TFUS. The review covers the mechanisms-of-action by which TFUS may cause adverse effects and the available data on the possible occurrence of such effects in animal and human studies.
Initial screening used key term searches in PubMed and bioRxiv, and a review of the literature lists of relevant papers. We included only studies where safety assessment was performed, and this results in 33 studies, both in humans and animals.
Adverse effects of TFUS were very rare. At high stimulation intensity and/or rate, TFUS may cause haemorrhage, cell death or damage, and unintentional blood-brain barrier (BBB) opening. TFUS may also unintentionally affect long-term neural activity and behaviour. A variety of methods was used mainly in rodents to evaluate these adverse effects, including tissue staining, magnetic resonance imaging, temperature measurements and monitoring of neural activity and behaviour. In 30 studies, adverse effects were absent, even though at least one Food and Drug Administration (FDA) safety index was frequently exceeded. Two studies reported microhaemorrhages after long or relatively intense stimulation above safety limits. Another study reported BBB opening and neuronal damage in a control condition, which intentionally and substantially exceeded the safety limits.
Most studies point towards a favourable safety profile of TFUS. Further investigations are warranted to establish a solid safety framework for the therapeutic window of TFUS to reliably avoid adverse effects while ensuring neural effectiveness. The comparability across studies should be improved by a more standardized reporting of TFUS parameters.
•TFUS is an emerging non-invasive stimulation method with excellent focality•We summarize the safety-related data from the available literature.•Adverse effects were absent in 30 studies, and were reported in 3 studies.•Many studies used parameters outside safety limits for diagnostic ultrasound.•Further studies are warranted to establish the safety margin for TFUS.
•A head template for acoustic and thermal simulations of TUS is introduced.•The template was built from calibrated CT scans of normal adults (age < 50).•Acoustic and thermal simulations based on the ...template represent the group median well.•The template is useful for treatment planning and optimization in young adults.
Transcranial focused Ultrasound Stimulation (TUS) at low intensities is emerging as a novel non-invasive brain stimulation method with higher spatial resolution than established transcranial stimulation methods and the ability to selectively stimulate also deep brain areas. Accurate control of the focus position and strength of the TUS acoustic waves is important to enable a beneficial use of the high spatial resolution and to ensure safety. As the human skull causes strong attenuation and distortion of the waves, simulations of the transmitted waves are needed to accurately determine the TUS dose distribution inside the cranial cavity. The simulations require information of the skull morphology and its acoustic properties. Ideally, they are informed by computed tomography (CT) images of the individual head. However, suited individual imaging data is often not readily available. For this reason, we here introduce and validate a head template that can be used to estimate the average effects of the skull on the TUS acoustic wave in the population.
The template was created from CT images of the heads of 29 individuals of different ages (between 20–50 years), gender and ethnicity using an iterative non-linear co-registration procedure. For validation, we compared acoustic and thermal simulations based on the template to the average of the simulation results of all 29 individual datasets. Acoustic simulations were performed for a model of a focused transducer driven at 500 kHz, placed at 24 standardized positions by means of the EEG 10–10 system. Additional simulations at 250 kHz and 750 kHz at 16 of the positions were used for further confirmation. The amount of ultrasound-induced heating at 500 kHz was estimated for the same 16 transducer positions. Our results show that the template represents the median of the acoustic pressure and temperature maps from the individuals reasonably well in most cases. This underpins the usefulness of the template for the planning and optimization of TUS interventions in studies of healthy young adults. Our results further indicate that the amount of variability between the individual simulation results depends on the position. Specifically, the simulated ultrasound-induced heating inside the skull exhibited strong interindividual variability for three posterior positions close to the midline, caused by a high variability of the local skull shape and composition. This should be taken into account when interpreting simulation results based on the template.
Hippocampal enlargements are commonly reported after electroconvulsive therapy (ECT). To clarify mechanisms, we examined if ECT-induced hippocampal volume change relates to dose (number of ECT ...sessions and electrode placement) and acts as a biomarker of clinical outcome.
Longitudinal neuroimaging and clinical data from 10 independent sites participating in the Global ECT-Magnetic Resonance Imaging Research Collaboration (GEMRIC) were obtained for mega-analysis. Hippocampal volumes were extracted from structural magnetic resonance images, acquired before and after patients (n = 281) experiencing a major depressive episode completed an ECT treatment series using right unilateral and bilateral stimulation. Untreated nondepressed control subjects (n = 95) were scanned twice.
The linear component of hippocampal volume change was 0.28% (SE 0.08) per ECT session (p < .001). Volume change varied by electrode placement in the left hippocampus (bilateral, 3.3 ± 2.2%, d = 1.5; right unilateral, 1.6 ± 2.1%, d = 0.8; p < .0001) but not the right hippocampus (bilateral, 3.0 ± 1.7%, d = 1.8; right unilateral, 2.7 ± 2.0%, d = 1.4; p = .36). Volume change for electrode placement per ECT session varied similarly by hemisphere. Individuals with greater treatment-related volume increases had poorer outcomes (Montgomery–Åsberg Depression Rating Scale change –1.0 SE 0.35, per 1% volume increase, p = .005), although the effects were not significant after controlling for ECT number (slope –0.69 SE 0.38, p = .069).
The number of ECT sessions and electrode placement impacts the extent and laterality of hippocampal enlargement, but volume change is not positively associated with clinical outcome. The results suggest that the high efficacy of ECT is not explained by hippocampal enlargement, which alone might not serve as a viable biomarker for treatment outcome.
Electroconvulsive therapy (ECT) is associated with volumetric enlargements of corticolimbic brain regions. However, the pattern of whole-brain structural alterations following ECT remains unresolved. ...Here, we examined the longitudinal effects of ECT on global and local variations in gray matter, white matter, and ventricle volumes in patients with major depressive disorder as well as predictors of ECT-related clinical response.
Longitudinal magnetic resonance imaging and clinical data from the Global ECT-MRI Research Collaboration (GEMRIC) were used to investigate changes in white matter, gray matter, and ventricle volumes before and after ECT in 328 patients experiencing a major depressive episode. In addition, 95 nondepressed control subjects were scanned twice. We performed a mega-analysis of single subject data from 14 independent GEMRIC sites.
Volumetric increases occurred in 79 of 84 gray matter regions of interest. In total, the cortical volume increased by mean ± SD of 1.04 ± 1.03% (Cohen’s d = 1.01, p < .001) and the subcortical gray matter volume increased by 1.47 ± 1.05% (d = 1.40, p < .001) in patients. The subcortical gray matter increase was negatively associated with total ventricle volume (Spearman’s rank correlation ρ = −.44, p < .001), while total white matter volume remained unchanged (d = −0.05, p = .41). The changes were modulated by number of ECTs and mode of electrode placements. However, the gray matter volumetric enlargements were not associated with clinical outcome.
The findings suggest that ECT induces gray matter volumetric increases that are broadly distributed. However, gross volumetric increases of specific anatomically defined regions may not serve as feasible biomarkers of clinical response.
Purpose
The diffusion‐weighted SPLICE (split acquisition of fast spin‐echo signals) sequence employs split‐echo rapid acquisition with relaxation enhancement (RARE) readout to provide images almost ...free of geometric distortions. However, due to the varying T2$$ {}_2 $$‐weighting during k‐space traversal, SPLICE suffers from blurring. This work extends a method for controlling the spatial point spread function (PSF) while optimizing the signal‐to‐noise ratio (SNR) achieved by adjusting the flip angles in the refocusing pulse train of SPLICE.
Methods
An algorithm based on extended phase graph (EPG) simulations optimizes the flip angles by maximizing SNR for a flexibly chosen predefined target PSF that describes the desired k‐space density weighting and spatial resolution. An optimized flip angle scheme and a corresponding post‐processing correction filter which together achieve the target PSF was tested by healthy subject brain imaging using a clinical 1.5 T scanner.
Results
Brain images showed a clear and consistent improvement over those obtained with a standard constant flip angle scheme. SNR was increased and apparent diffusion coefficient estimates were more accurate. For a modified Hann k‐space weighting example, considerable benefits resulted from acquisition weighting by flip angle control.
Conclusion
The presented flexible method for optimizing SPLICE flip angle schemes offers improved MR image quality of geometrically accurate diffusion‐weighted images that makes the sequence a strong candidate for radiotherapy planning or stereotactic surgery.
Purpose
To demonstrate a novel method for tracking of head movements during MRI using electroencephalography (EEG) hardware for recording signals induced by native imaging gradients.
Theory and ...Methods
Gradient switching during simultaneous EEG–fMRI induces distortions in EEG signals, which depend on subject head position and orientation. When EEG electrodes are interconnected with high‐impedance carbon wire loops, the induced voltages are linear combinations of the temporal gradient waveform derivatives. We introduce head tracking based on these signals (CapTrack) involving 3 steps: (1) phantom scanning is used to characterize the target sequence and a fast calibration sequence; (2) a linear relation between changes of induced signals and head pose is established using the calibration sequence; and (3) induced signals recorded during target sequence scanning are used for tracking and retrospective correction of head movement without prolonging the scan time of the target sequence. Performance of CapTrack is compared directly to interleaved navigators.
Results
Head‐pose tracking at 27.5 Hz during echo planar imaging (EPI) was demonstrated with close resemblance to rigid body alignment (mean absolute difference: 0.14 0.38 0.15‐mm translation, 0.30 0.27 0.22‐degree rotation). Retrospective correction of 3D gradient‐echo imaging shows an increase of average edge strength of 12%/−0.39% for instructed/uninstructed motion with CapTrack pose estimates, with a tracking interval of 1561 ms and high similarity to interleaved navigator estimates (mean absolute difference: 0.13 0.33 0.12 mm, 0.28 0.15 0.22 degrees).
Conclusion
Motion can be estimated from recordings of gradient switching with little or no sequence modification, optionally in real time at low computational burden and synchronized to image acquisition, using EEG equipment already found at many research institutions.
Magnetic resonance current density imaging (MRCDI) of the human brain aims to reconstruct the current density distribution caused by transcranial electric stimulation from MR-based measurements of ...the current-induced magnetic fields. So far, the MRCDI data acquisition achieves only a low signal-to-noise ratio, does not provide a full volume coverage and lacks data from the scalp and skull regions. In addition, it is only sensitive to the component of the current-induced magnetic field parallel to the scanner field. The reconstruction problem thus involves coping with noisy and incomplete data, which makes it mathematically challenging. Most existing reconstruction methods have been validated using simulation studies and measurements in phantoms with simplified geometries. Only one reconstruction method, the projected current density algorithm, has been applied to human in-vivo data so far, however resulting in blurred current density estimates even when applied to noise-free simulated data.
We analyze the underlying causes for the limited performance of the projected current density algorithm when applied to human brain data. In addition, we compare it with an approach that relies on the optimization of the conductivities of a small number of tissue compartments of anatomically detailed head models reconstructed from structural MR data. Both for simulated ground truth data and human in-vivo MRCDI data, our results indicate that the estimation of current densities benefits more from using a personalized volume conductor model than from applying the projected current density algorithm. In particular, we introduce a hierarchical statistical testing approach as a principled way to test and compare the quality of reconstructed current density images that accounts for the limited signal-to-noise ratio of the human in-vivo MRCDI data and the fact that the ground truth of the current density is unknown for measured data. Our results indicate that the statistical testing approach constitutes a valuable framework for the further development of accurate volume conductor models of the head. Our findings also highlight the importance of tailoring the reconstruction approaches to the quality and specific properties of the available data.
Magnetic resonance current density imaging (MRCDI) and MR electrical impedance tomography (MREIT) are two emerging modalities, which combine weak time-varying currents injected via surface electrodes ...with magnetic resonance imaging (MRI) to acquire information about the current flow and ohmic conductivity distribution at high spatial resolution. The injected current flow creates a magnetic field in the head, and the component of the induced magnetic field ΔBz,c parallel to the main scanner field causes small shifts in the precession frequency of the magnetization. The measured MRI signal is modulated by these shifts, allowing to determine ΔBz,c for the reconstruction of the current flow and ohmic conductivity.
Here, we demonstrate reliable ΔBz,c measurements in-vivo in the human brain based on multi-echo spin echo (MESE) and steady-state free precession free induction decay (SSFP-FID) sequences. In a series of experiments, we optimize their robustness for in-vivo measurements while maintaining a good sensitivity to the current-induced fields. We validate both methods by assessing the linearity of the measured ΔBz,c with respect to the current strength. For the more efficient SSFP-FID measurements, we demonstrate a strong influence of magnetic stray fields on the ΔBz,c images, caused by non-ideal paths of the electrode cables, and validate a correction method. Finally, we perform measurements with two different current injection profiles in five subjects. We demonstrate reliable recordings of ΔBz,c fields as weak as 1 nT, caused by currents of 1 mA strength. Comparison of the ΔBz,c measurements with simulated ΔBz,c images based on FEM calculations and individualized head models reveals significant linear correlations in all subjects, but only for the stray field-corrected data. As final step, we reconstruct current density distributions from the measured and simulated ΔBz,c data. Reconstructions from non-corrected ΔBz,c measurements systematically overestimate the current densities. Comparing the current densities reconstructed from corrected ΔBz,c measurements and from simulated ΔBz,c images reveals an average coefficient of determination R2 of 71%. In addition, it shows that the simulations underestimated the current strength on average by 24%.
Our results open up the possibility of using MRI to systematically validate and optimize numerical field simulations that play an important role in several neuroscience applications, such as transcranial brain stimulation, and electro- and magnetoencephalography.
In the field of radiation oncology, the benefit of MRI goes beyond that of providing high soft-tissue contrast images for staging and treatment planning. With the recent clinical introduction of ...hybrid MRI linear accelerators it has become feasible to map physiological parameters describing diffusion, perfusion, and relaxation during the entire course of radiotherapy, for example. However, advanced data analysis tools are required for extracting qualified prognostic and predictive imaging biomarkers from longitudinal MRI data. In this study, we propose a new prediction framework tailored to exploit temporal dynamics of tissue features from repeated measurements. We demonstrate the framework using a newly developed decomposition method for tumor characterization.
Two previously published MRI datasets with multiple measurements during and after radiotherapy, were used for development and testing:
-weighted multi-echo images obtained for two mouse models of pancreatic cancer, and diffusion-weighted images for patients with brain metastases. Initially, the data was decomposed using the novel monotonous slope non-negative matrix factorization (msNMF) tailored for MR data. The following processing consisted of a tumor heterogeneity assessment using descriptive statistical measures, robust linear modelling to capture temporal changes of these, and finally logistic regression analysis for stratification of tumors and volumetric outcome.
The framework was able to classify the two pancreatic tumor types with an area under curve (AUC) of 0.999,
< 0.001 and predict the tumor volume change with a correlation coefficient of 0.513,
= 0.034. A classification of the human brain metastases into responders and non-responders resulted in an AUC of 0.74,
= 0.065.
A general data processing framework for analyses of longitudinal MRI data has been developed and applications were demonstrated by classification of tumor type and prediction of radiotherapy response. Further, as part of the assessment, the merits of msNMF for tumor tissue decomposition were demonstrated.