Volumetric measurements of weak current-induced magnetic fields in the human brain at high resolution

PURPOSE

Clinical use of transcranial electrical stimulation (TES) requires accurate knowledge of the injected current distribution in the brain. MR current density imaging (MRCDI) uses measurements of the TES-induced magnetic fields to provide this information. However, sufficient sensitivity and image quality in humans in vivo has only been documented for single-slice imaging.

METHODS

A recently developed, optimally spoiled, acquisition-weighted, gradient echo-based 2D-MRCDI method has now been advanced for volume coverage with densely or sparsely distributed slices: The 3D rectilinear sampling (3D-DENSE) and simultaneous multislice acquisition (SMS-SPARSE) were optimized and verified by cable-loop experiments and tested with 1-mA TES experiments for two common electrode montages.

RESULTS

Comparisons between the volumetric methods against the 2D-MRCDI showed that relatively long acquisition times of 3D-DENSE using a single slab with six slices hindered the expected sensitivity improvement in the current-induced field measurements but improved sensitivity by 61% in the Laplacian of the field, on which some MRCDI reconstruction methods rely. Also, SMS-SPARSE acquisition of three slices, with a factor 2 CAIPIRINHA (controlled aliasing in parallel imaging results in higher acceleration) acceleration, performed best against the 2D-MRCDI with sensitivity improvements for the∆ B z , c $$ \Delta {B}_{z,c} $$ and Laplacian noise floors of 56% and 78% (baseline without current flow) as well as 43% and 55% (current injection into head). SMS-SPARSE reached a sensitivity of 67 pT for three distant slices at 2 × 2 × 3 mm3 resolution in 10 min of total scan time, and consistently improved image quality.

CONCLUSION

Volumetric MRCDI measurements with high sensitivity and image quality are well suited to characterize the TES field distribution in the human brain.

A head template for computational dose modelling for transcranial focused ultrasound stimulation

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.

Optimized flip angle schemes for the split acquisition of fast spin-echo signals (SPLICE) sequence and application to diffusion-weighted imaging

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 T 2 $$ {}_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.

Decomposition-based framework for tumor classification and prediction of treatment response from longitudinal MRI

Objective.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.Approach.Two previously published MRI datasets with multiple measurements during and after radiotherapy, were used for development and testing:T₂-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.Main Results.The framework was able to classify the two pancreatic tumor types with an area under curve (AUC) of 0.999,P< 0.001 and predict the tumor volume change with a correlation coefficient of 0.513,P= 0.034. A classification of the human brain metastases into responders and non-responders resulted in an AUC of 0.74,P= 0.065.Significance.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.

Tracking of rigid head motion during MRI using an EEG system

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.

Data-driven separation of MRI signal components for tissue characterization

PURPOSE

MRI can be utilized for quantitative characterization of tissue. To assess e.g. water fractions or diffusion coefficients for compartments in the brain, a decomposition of the signal is necessary. Imposing standard models carries the risk of estimating biased parameters if model assumptions are violated. This work introduces a data-driven multicomponent analysis, the monotonous slope non-negative matrix factorization (msNMF), tailored to extract data features expected in MR signals.

METHODS

The msNMF was implemented by extending the standard NMF with monotonicity constraints on the signal profiles and their first derivatives. The method was validated using simulated data, and subsequently applied to both ex vivo DWI data and in vivo relaxometry data. Reproducibility of the method was tested using the latter.

RESULTS

The msNMF recovered the multi-exponential signals in the simulated data and showed superiority to standard NMF (based on the explained variance, area under the ROC curve, and coefficient of variation). Diffusion components extracted from the DWI data reflected the cell density of the underlying tissue. The relaxometry analysis resulted in estimates of edema water fractions (EWF) highly correlated with published results, and demonstrated acceptable reproducibility.

CONCLUSION

The msNMF can robustly separate MR signals into components with relation to the underlying tissue composition, and may potentially be useful for e.g. tumor tissue characterization.

On the reconstruction of magnetic resonance current density images of the human brain: Pitfalls and perspectives

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.

Sensitivity and resolution improvement for in vivo magnetic resonance current-density imaging of the human brain

PURPOSE

Magnetic resonance current-density imaging (MRCDI) combines MRI with low-intensity transcranial electrical stimulation (TES; 1-2 mA) to map current flow in the brain. However, usage of MRCDI is still hampered by low measurement sensitivity and image quality.

METHODS

Recently, a multigradient-echo-based MRCDI approach has been introduced that presently has the best-documented efficiency. This MRCDI approach has now been advanced in three directions and has been validated by phantom and in vivo experiments. First, the importance of optimum spoiling for brain imaging was verified. Second, the sensitivity and spatial resolution were improved by using acquisition weighting. Third, navigators were added as a quality control measure for tracking physiological noise. Combining these advancements, the optimized MRCDI method was tested by using 1 mA TES for two different injection profiles.

RESULTS

For a session duration of 4:20 min, the new MRCDI method was able to detect TES-induced magnetic fields at a sensitivity level of 84 picotesla, representing a twofold efficiency increase against our original method. A comparison between measurements and simulations based on personalized head models showed a consistent increase in the coefficient of determination of ΔR² = 0.12 for the current-induced magnetic fields and ΔR² = 0.22 for the current flow reconstructions. Interestingly, some of the simulations still clearly deviated from the measurements despite the strongly improved measurement quality. This highlights the utility of MRCDI to improve head models for TES simulations.

CONCLUSION

The achieved sensitivity improvement is an important step from proof-of-concept studies toward a broader application of MRCDI in clinical and basic neuroscience research.

Elevated body weight modulates subcortical volume change and associated clinical response following electroconvulsive therapy

BACKGROUND

Obesity is a frequent somatic comorbidity of major depression, and it has been associated with worse clinical outcomes and brain structural abnormalities. Converging evidence suggests that electroconvulsive therapy (ECT) induces both clinical improvements and increased subcortical grey matter volume in patients with depression. However, it remains unknown whether increased body weight modulates the clinical response and structural neuroplasticity that occur with ECT.

METHODS

To address this question, we conducted a longitudinal investigation of structural MRI data from the Global ECT-MRI Research Collaboration (GEMRIC) in 223 patients who were experiencing a major depressive episode (10 scanning sites). Structural MRI data were acquired before and after ECT, and we assessed change in subcortical grey matter volume using FreeSurfer and Quarc.

RESULTS

Higher body mass index (BMI) was associated with a significantly lower increase in subcortical grey matter volume following ECT. We observed significant negative associations between BMI and change in subcortical grey matter volume, with pronounced effects in the thalamus and putamen, where obese participants showed increases in grey matter volume that were 43.3% and 49.6%, respectively, of the increases found in participants with normal weight. As well, BMI significantly moderated the association between subcortical grey matter volume change and clinical response to ECT. We observed no significant association between BMI and clinical response to ECT.

LIMITATIONS

Because only baseline BMI values were available, we were unable to study BMI changes during ECT and their potential association with clinical and grey matter volume change.

CONCLUSION

Future studies should take into account the relevance of body weight as a modulator of structural neuroplasticity during ECT treatment and aim to further explore the functional relevance of this novel finding.

Safety evaluation of a new setup for transcranial electric stimulation during magnetic resonance imaging

BACKGROUND

Transcranial electric stimulation during MR imaging can introduce safety issues due to coupling of the RF field with the stimulation electrodes and leads.

OBJECTIVE

To optimize the stimulation setup for MR current density imaging (MRCDI) and increase maximum stimulation current, a new low-conductivity (σ = 29.4 S/m) lead wire is designed and tested.

METHOD

The antenna effect was simulated to investigate the effect of lead conductivity. Subsequently, specific absorption rate (SAR) simulations for realistic lead configurations with low-conductivity leads and two electrode types were performed at 128 MHz and 298 MHz being the Larmor frequencies of protons at 3T and 7T. Temperature measurements were performed during MRI using high power deposition sequences to ensure that the electrodes comply with MRI temperature regulations.

RESULTS

The antenna effect was found for copper leads at ¼ RF wavelength and could be reliably eliminated using low-conductivity leads. Realistic lead configurations increased the head SAR and the local head SAR at the electrodes only minimally. The highest temperatures were measured on the rings of center-surround electrodes, while circular electrodes showed little heating. No temperature increase above the safety limit of 39 °C was observed.

CONCLUSION

Coupling to the RF field can be reliably prevented by low-conductivity leads, enabling cable paths optimal for MRCDI. Compared to commercial copper leads with safety resistors, the low-conductivity leads had lower total impedance, enabling the application of higher currents without changing stimulator design. Attention must be paid to electrode pads.

Multi-site benchmarking of clinical 13C RF coils at 3T

A quality assurance protocol for RF coils is proposed, which can be used for volume (Tx/Rx) and surface (Rx) coils. Following this protocol, a benchmarking of seven coils (from three different MR sites) dedicated to ¹³C MRI at 3T is reported. Coil performance is particularly important for 3T MRI at the ¹³C frequency, since the coil-to-sample noise ratio is typically high. The coils are evaluated experimentally using the proposed protocol based on MR spectroscopic imaging performed with two different phantoms: one head-shaped, and one with cylindrical shape and nearly twice the volume of the first one. To achieve an unbiased SNR comparison of volume and array coils, coil combination was done using sensitivity profiles extracted from the data. SNR, noise correlation matrices and example g-factor maps are reported. For globally calibrated, equal excitation angles, the measured SNR shows large differences for the volume coils of up to 115% at the phantom center for a head phantom. The arrays show lower differences in superficial SNR. The sample surface depth at which the volume coils outperform the arrays is estimated to 7 cm, and SNR furthest away from the coil surface is 28% lower for the best array compared to the best volume coil. A broad set of coils for ¹³C at 3T have been benchmarked. The results reported, and the method used to benchmark them, should guide the ¹³C community to choose the most suitable coil for a given experiment.

Brain Changes Induced by Electroconvulsive Therapy Are Broadly Distributed

BACKGROUND

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.

METHODS

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.

RESULTS

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.

CONCLUSIONS

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.

Safety of transcranial focused ultrasound stimulation: A systematic review of the state of knowledge from both human and animal studies

BACKGROUND

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.

OBJECTIVE

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.

METHODS

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.

RESULTS

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.

CONCLUSION

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.

[Magnetic resonance imaging for the diagnosis of non-alcoholic fatty liver disease]

Non-alcoholic fatty liver disease (NAFLD) is one of the leading causes of chronic liver disease with an estimated overall prevalence of 25% in the global adult population. Liver biopsy is the gold standard for the diagnosis of NAFLD. However, the risk of complications and collection of only 1/50,000 of the total liver volume, limits this diagnostic method in an unselected population. Non-invasive diag-nostic methods are warranted, and magnetic resonance imaging of the liver for NAFLD has shown promising results.

Development of a Symmetric Echo-Planar Spectroscopy Imaging Framework for Hyperpolarized 13C Imaging in a Clinical PET/MR Scanner

Here, we developed a symmetric echo-planar spectroscopic imaging (EPSI) sequence for hyperpolarized 13C imaging on a clinical hybrid positron emission tomography/magnetic resonance imaging system. The pulse sequence uses parallel reconstruction pipelines to separately reconstruct data from odd-and-even gradient echoes to reduce artifacts from gradient imbalances. The ramp-sampled data in the spatiotemporal frequency space are regridded to compensate for the chemical-shift displacements. Unaliasing of nonoverlapping peaks outside of the sampled spectral width was performed to double the effective spectral width. The sequence was compared with conventional phase-encoded chemical-shift imaging (CSI) in phantoms, and it was evaluated in a canine cancer patient with ameloblastoma after injection of hyperpolarized [1-13C]pyruvate. The relative signal-to-noise ratio of EPSI with respect to CSI was 0.88, which is consistent with the decrease in sampling efficiency due to ramp sampling. Data regridding in the spatiotemporal frequency space significantly reduced spatial blurring compared with direct fast Fourier transform. EPSI captured the spatial distributions of both metabolites and their temporal dynamics in vivo with an in-plane spatial resolution of 5 × 9 mm2 and a temporal resolution of 3 seconds. Significantly higher spatial and temporal resolution for delineating anatomical structures in vivo was achieved for EPSI metabolic maps than for CSI maps, which suffered spatiotemporal blurring. The EPSI sequence showed promising results in terms of short acquisition time and sufficient spectral bandwidth of 500 Hz, allowing to adjust the trade-off between signal-to-noise ratio and encoding speed.

Gamma-aminobutyric acid edited echo-planar spectroscopic imaging (EPSI) with MEGA-sLASER at 7T

PURPOSE

For rapid spatial mapping of gamma-aminobutyric acid (GABA) at the increased sensitivity and spectral separation for ultra-high magnetic field strength (7 tesla [T]), an accelerated edited magnetic resonance spectroscopic imaging technique was developed and optimized for the human brain at 7 T.

METHODS

A MEGA-sLASER sequence was used for GABA editing and volume selection to maximize editing efficiency and minimize chemical shift displacement errors. To accommodate the high bandwidth requirements at 7 T, a single-shot echo planar readout was used for rapid simultaneous encoding of the temporal dimension and 1 spatial. B₀ and B₁ field aspects specific for 7 T were studied together with correction procedures, and feasibility of the EPSI MEGA-sLASER technique was tested in vivo in 5 healthy subjects.

RESULTS

Localized edited spectra could be measured in all subjects giving spatial GABA signal distributions over a central brain region, having 45- to 50-Hz spatial intervoxel B₀ field variations and up to 30% B₁ field deviations. MEGA editing was found unaffected by the B₀ inhomogeneities for the optimized sequence. The correction procedures reduced effects of intervoxel B₀ inhomogeneities, corrected for spatial editing efficiency variations, and compensated for GABA resonance phase and frequency shifts from subtle motion and acquisition instabilities. The optimized oscillating echo-planar gradient scheme permitted full spectral acquisition at 7 T and exhibited minimal spectral-spatial ghosting effects for the selected brain region.

CONCLUSION

The EPSI MEGA-sLASER technique was shown to provide time-efficient mapping of regional variations in cerebral GABA in a central volume of interest with spatial B₁ and B₀ field variations typical for 7 T.

Human in-vivo brain magnetic resonance current density imaging (MRCDI)

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 ΔB_z,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 ΔB_z,c for the reconstruction of the current flow and ohmic conductivity. Here, we demonstrate reliable ΔB_z,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 ΔB_z,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 ΔB_z,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 ΔB_z,c fields as weak as 1 nT, caused by currents of 1 mA strength. Comparison of the ΔB_z,c measurements with simulated ΔB_z,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 ΔB_z,c data. Reconstructions from non-corrected ΔB_z,c measurements systematically overestimate the current densities. Comparing the current densities reconstructed from corrected ΔB_z,c measurements and from simulated ΔB_z,c images reveals an average coefficient of determination R² 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.

Effect of a treat-to-target strategy based on methotrexate and intra-articular betamethasone with or without additional cyclosporin on MRI-assessed synovitis, osteitis, tenosynovitis, bone erosion, and joint space narrowing in early rheumatoid arthritis: results from a 2-year randomized double-blind placebo-controlled trial (CIMESTRA)

OBJECTIVES

To investigate whether a treat-to-target strategy based on methotrexate (MTX) and intra-articular (IA) betamethasone suppresses magnetic resonance imaging (MRI)-determined measures of disease activity and reduces joint destruction in early rheumatoid arthritis (eRA) patients, and to investigate whether concomitant cyclosporin A (CyA) provides an additional effect.

METHOD

In the 2-year randomized, double-blind, treat-to-target trial CIMESTRA, 160 patients with eRA (< 6 months) were randomized to MTX, intra-articular betamethasone and CyA, or placebo CyA. A total of 129 patients participated in the MRI substudy, and had contrast-enhanced MR images of the non-dominant hand at months 0, 6, 12, and 24. MR images were evaluated for osteitis, synovitis, tenosynovitis, bone erosion, and joint space narrowing (JSN), using validated scoring methods.

RESULTS

Significant reductions were seen at 6 months in all inflammatory parameters [synovitis, mean change -1.6 (p < 0.001, Wilcoxon), tenosynovitis, -3.5 (p < 0.001), and osteitis, -1.3 (p < 0.05)] and at 12/24 months in synovitis and tenosynovitis [-1.6/-2.2 and -3.6/-3.8, respectively; all p < 0.001]. MRI signs of inflammation were not fully eliminated, and increases in erosion and JSN scores were observed at 6 months [0.4 (p < 0.01)/0.1 (p < 0.05)], 12 months [0.8 (p < 0.001)/0.3 (p < 0.01)], and 24 months [1.0 (p < 0.001)/0.4 (p < 0.001)]. Clinical measures decreased significantly (p < 0.001) at all time points. There were no consistent statistically significant differences between treatment groups.

CONCLUSIONS

In this eRA treat-to-target trial, MTX and intra-articular glucocorticoids markedly reduced, but did not eliminate, MRI osteitis, synovitis, and tenosynovitis. Accordingly, minimal but statistically significant increases in bone erosion and JSN were observed. No additional effect of CyA was demonstrated.

Regional brain volumes, diffusivity, and metabolite changes after electroconvulsive therapy for severe depression

OBJECTIVE

To investigate the role of hippocampal plasticity in the antidepressant effect of electroconvulsive therapy (ECT).

METHOD

We used magnetic resonance (MR) imaging including diffusion tensor imaging (DTI) and proton MR spectroscopy (¹ H-MRS) to investigate hippocampal volume, diffusivity, and metabolite changes in 19 patients receiving ECT for severe depression. Other regions of interest included the amygdala, dorsolateral prefrontal cortex (DLPFC), orbitofrontal cortex, and hypothalamus. Patients received a 3T MR scan before ECT (TP1), 1 week (TP2), and 4 weeks (TP3) after ECT.

RESULTS

Hippocampal and amygdala volume increased significantly at TP2 and continued to be increased at TP3. DLPFC exhibited a transient volume reduction at TP2. DTI revealed a reduced anisotropy and diffusivity of the hippocampus at TP2. We found no significant post-ECT changes in brain metabolite concentrations, and we were unable to identify a spectral signature at ≈1.30 ppm previously suggested to reflect neurogenesis induced by ECT. None of the brain imaging measures correlated to the clinical response.

CONCLUSION

Our findings show that ECT causes a remodeling of brain structures involved in affective regulation, but due to their lack of correlation with the antidepressant effect, this remodeling does not appear to be directly underlying the antidepressant action of ECT.

Measuring motion-induced B0 -fluctuations in the brain using field probes

PURPOSE

Fluctuations of the background magnetic field (B0 ) due to body and breathing motion can lead to significant artifacts in brain imaging at ultrahigh field. Corrections based on real-time sensing using external field probes show great potential. This study evaluates different aspects of field interpolation from these probes into the brain which is implicit in such methods. Measurements and simulations were performed to quantify how well B0 -fluctuations in the brain due to body and breathing motion are reflected in external field probe measurements.

METHODS

Field probe measurements were compared with scanner acquired B0 -maps from experiments with breathing and shoulder movements. A realistic simulation of B0 -fluctuations caused by breathing was performed, and used for testing different sets of field probe positions.

RESULTS

The B0 -fluctuations were well reflected in the field probe measurements in the shoulder experiments, while the breathing experiments showed only moderate correspondence. The simulations showed the importance of the probe positions, and that performing full 3(rd) order corrections based on 16 field probes is not recommended.

CONCLUSION

Methods for quantitative assessment of the field interpolation problem were developed and demonstrated. Field corrections based on external field measurements show great potential, although potential pitfalls were identified.

Blood-brain barrier permeability of normal appearing white matter in relapsing-remitting multiple sclerosis

BACKGROUND

Multiple sclerosis (MS) affects the integrity of the blood-brain barrier (BBB). Contrast-enhanced T1 weighted magnetic resonance imaging (MRI) is widely used to characterize location and extent of BBB disruptions in focal MS lesions. We employed quantitative T1 measurements before and after the intravenous injection of a paramagnetic contrast agent to assess BBB permeability in the normal appearing white matter (NAWM) in patients with relapsing-remitting MS (RR-MS).

METHODOLOGY/PRINCIPAL FINDINGS

Fifty-nine patients (38 females) with RR-MS undergoing immunomodulatory treatment and nine healthy controls (4 females) underwent quantitative T1 measurements at 3 tesla before and after injection of a paramagnetic contrast agent (0.2 mmol/kg Gd-DTPA). Mean T1 values were calculated for NAWM in patients and total cerebral white matter in healthy subjects for the T1 measurements before and after injection of Gd-DTPA. The pre-injection baseline T1 of NAWM (945±55 [SD] ms) was prolonged in RR-MS relative to healthy controls (903±23 ms, p = 0.028). Gd-DTPA injection shortened T1 to a similar extent in both groups. Mean T1 of NAWM was 866±47 ms in the NAWM of RR-MS patients and 824±13 ms in the white matter of healthy controls. The regional variability of T1 values expressed as the coefficient of variation (CV) was comparable between the two groups at baseline, but not after injection of the contrast agent. After intravenous Gd-DTPA injection, T1 values in NAWM were more variable in RR-MS patients (CV = 0.198±0.046) compared to cerebral white matter of healthy controls (CV = 0.166±0.018, p = 0.046).

CONCLUSIONS/SIGNIFICANCE

We found no evidence of a global BBB disruption within the NAWM of RR-MS patients undergoing immunomodulatory treatment. However, the increased variation of T1 values in NAWM after intravenous Gd-DTPA injection points to an increased regional inhomogeneity of BBB function in NAWM in relapsing-remitting MS.

Brain N-acetylaspartate levels correlate with motor function in metachromatic leukodystrophy

BACKGROUND

Late infantile metachromatic leukodystrophy (MLD) is an autosomal recessive lysosomal storage disorder that causes severe demyelination of the nervous system. The neuronal metabolite N-acetylaspartate (NAA) serves as a source of acetyl groups for myelin lipid synthesis in oligodendrocytes and is known as a marker for neuronal and axonal loss. NAA and other metabolite levels measured by proton magnetic resonance spectroscopy (MRS) correlate with performance of the brain in normal children. There is a need for sensitive measures of disease progression in patients with MLD to enable development of future treatments.

METHODS

A cross-section of 13 children with late infantile MLD were examined by proton MRS. Signals from NAA, total choline, and total creatine in the deep white matter were measured and correlated with the results of cognitive and motor function tests.

RESULTS

The NAA signal decreased as the disease process advanced. Motor function, measured by the Gross Motor Function Measure-88, varied from 13 (only head movement in the supine position) to 180 (able to walk) across the study cohort, demonstrating a wide range in functional status. Similarly, varied decreases were observed in cognitive function. We report strong positive correlations between standardized measures of motor and cognitive function and NAA levels in the deep white matter.

CONCLUSIONS

We suggest that NAA levels could serve as a sensitive biomarker in children with MLD. Proton MRS may provide a valuable tool for measuring the effects of treatment interventions in this disorder.

A Graphical Simulator for Teaching Basic and Advanced MR Imaging Techniques

Teaching of magnetic resonance (MR) imaging techniques typically involves considerable handwaving, literally, to explain concepts such as resonance, rotating frames, dephasing, refocusing, sequences, and imaging. A proper understanding of MR contrast and imaging techniques is crucial for radiologists, radiographers, and technical staff alike, but it is notoriously challenging to explain spin dynamics by using traditional teaching tools. The author developed a freely available graphical simulator based on the Bloch equations to aid in the teaching of topics ranging from precession and relaxation to advanced concepts such as stimulated echoes, spin tagging, and k-space-methods. A graphical user interface provides the user with a three-dimensional view of spin isochromates that can be manipulated by selecting radiofrequency pulses and gradient events. Even complicated sequences can be visualized in an intuitive way. The cross-platform software is primarily designed for use in lectures, but is also useful for self studies and student assignments. Movies available at http://radiographics.rsnajnls.org/cgi/content/full/e27/DC1.

Encoding of electrophysiology and other signals in MR images

PURPOSE

To develop a gradient insensitive, generic technique for recording of non-MR signals by use of surplus scanner bandwidth.

MATERIALS AND METHODS

Relatively simple battery driven hardware is used to transform one or more signals into radio waves detectable by the MR scanner. Similar to the "magstripe" technique used for encoding of soundtracks in motion pictures, the electrical signals are in this way encoded as artifacts appearing in the MR images or spectra outside the region of interest. The encoded signals are subsequently reconstructed from the signal recorded by the scanner.

RESULTS

Electrophysiological (EP) eye and heart muscular recording (electrooculography [EOG] and electrocardiography [ECG]) during fast echo planar imaging (EPI) is demonstrated with an expandable, modular 8-channel prototype implementation. The gradient artifacts that would normally be dominating EOG are largely eliminated.

CONCLUSION

The method provides relatively inexpensive sampling with inherent microsecond synchronization and it reduces gradient artifacts in physiological recordings significantly. When oversampling is employed, the method is compatible with all MR reconstruction and postprocessing techniques.

The relationship between MRI and PET changes and cognitive disturbances in MS

Cognitive dysfunction in multiple sclerosis (MS) is present in approximately 50% of the patients. Only moderate correlations have been found between cognitive dysfunction and T(2) lesion load, black holes or atrophy. Cognitive dysfunction in MS is probably related to the overall disease burden of the brain including abnormalities in normal appearing white matter (NAWM) and cortical grey matter, which is undetected with conventional magnetic resonance imaging (MRI). Hence, imaging techniques that embrace such abnormalities are needed to achieve better correlation with cognitive dysfunction. MR spectroscopy (MRS) performed with multi-slice echo planar spectroscopic imaging (EPSI) and PET measurements of brain metabolism as the cortical cerebral metabolic rate of glucose are imaging methods that are able to provide information on axonal loss or dysfunction in both MS lesions and in NAWM and cortical grey matter. Measurements of global NAA using multi-slice EPSI is a new promising method for measurement of the global neuron capacity and can be repeated with only little discomfort and without any risk for the patient.

Magnetic resonance imaging at 3.0 tesla detects more lesions in acute optic neuritis than at 1.5 tesla

OBJECTIVE

We sought to assess whether magnetic resonance imaging (MRI) at 3.0 T detects more brain lesions in acute optic neuritis (ON) than MRI at 1.5 T.

MATERIALS AND METHODS

Twenty-eight patients with acute ON were scanned at both field-strengths using fast-fluid-attenuated inversion recovery (FLAIR), proton density and T2-weighted turbo spin echo, and T1-weighted spin echo after contrast. In addition, magnetization-prepared rapid acquisition gradient echo (MPRAGE) was obtained after contrast at 3.0 T. Lesion number and volumes were assessed by an observer blind to patient identity and field strength.

RESULTS

Scans at 3.0 T showed a significantly increase in number of lesions detected on FLAIR images (P = 0.002) relative to scanning at 1.5 T. MPRAGE proved to be suitable for detecting enhancing lesions in ON.

CONCLUSION

The MRI protocol at 3.0 T was more sensitive to hyperintense brain lesions in ON than the standard MRI protocol at 1.5 T.

Correlation of Global N-Acetyl Aspartate With Cognitive Impairment in Multiple Sclerosis

BACKGROUND

Whole-brain N-acetyl aspartate (NAA), a measure of neuronal function, can be assessed by multislice echo-planar spectroscopic imaging.

OBJECTIVE

To test the hypothesis that the global brain NAA/creatine (Cr) ratio is a better predictor of cognitive dysfunction in multiple sclerosis than conventional magnetic resonance imaging measures.

DESIGN

Survey.

SETTING

Research-oriented hospitals.

PATIENTS

Twenty patients, 16 women and 4 men (mean age, 36 years), with early relapsing-remitting multiple sclerosis (mean Expanded Disability Status Scale score, 2.5).

MAIN OUTCOME MEASURES

Correlation between the global NAA/Cr ratio and a cognitive dysfunction factor comprising 16 measures from an extensive neuropsychological test battery that best distinguished patients with multiple sclerosis from healthy control subjects.

RESULTS

A significant partial correlation between the global NAA/Cr ratio and the cognitive dysfunction factor was found (partial r = 0.62, P = .01), and 9 cognitively impaired patients had significantly lower global NAA/Cr ratios than 11 unimpaired patients (P = .04). No significant correlations were found between the cognitive dysfunction factor and conventional magnetic resonance imaging measures (ie, brain parenchymal fraction and lesion volume).

CONCLUSIONS

Multislice echo-planar spectroscopic imaging provides global metabolic measures that distinguish between cognitively impaired and unimpaired patients with multiple sclerosis and correlate with a global cognitive measure. Standardization of the technique is needed, and larger-scale studies that include healthy controls are suggested.

Multi-slice echo-planar spectroscopic MR imaging provides both global and local metabolite measures in multiple sclerosis

MR spectroscopy (MRS) provides information about neuronal loss or dysfunction by measuring decreases in N-acetyl aspartate (NAA), a metabolite widely believed to be a marker of neuronal viability. In multiple sclerosis (MS), whole-brain NAA (WBNAA) has been suggested as a marker of disease progression and treatment efficacy in treatment trials, and the ability to measure NAA loss in specific brain regions early in the evolution of this disease may have prognostic value. Most spectroscopic studies to date have been limited to single voxels or nonlocalized measurements of WBNAA only, and longitudinal studies have often been hampered by standardization and reproducibility problems. Multi-slice echo-planar spectroscopic imaging (EPSI) is presented as a promising alternative to single-voxel or nonlocalized spectroscopy for obtaining global metabolite estimates in MS. In the same session, measurements of metabolites in specific brain areas chosen after image acquisition (e.g., normal-appearing white matter (NAWM), gray matter (GM), and lesions) can be obtained. The identification and exclusion of regions that are inadequate for spectroscopic evaluation in global assessments can significantly improve quality and reproducibility, as demonstrated by a low within-subject variance in healthy controls. The reproducibility of the technique makes it a promising tool for future longitudinal spectroscopic studies of MS.

Hyperpolarized 3He MRI and 81mKr SPECT in chronic obstructive pulmonary disease

PURPOSE

During recent years, magnetic resonance imaging (MRI) using hyperpolarised (HP) 3He gas has emerged as a promising new method for the imaging of lung ventilation. However, systematic comparisons with nuclear medicine techniques have not yet been performed. The aim of this study was to compare ventilation imaging methods in 26 patients with chronic obstructive pulmonary disease (COPD) and nine lung healthy volunteers.

METHODS

HP 3He MRI, 81mKr single-photon emission computed tomography (SPECT), high-resolution computed tomography (HRCT) and pulmonary function tests were performed. The three scans were scored visually as percentage of non-ventilated/diseased lung, and a computer-based objective measure of the ventilated volume in HP 3He MRI and 81mKr SPECT and an emphysema index in HRCT were calculated.

RESULTS

We found a good correlation between HP 3He MRI and 81mKr SPECT for both visual defect score (r=0.80, p<0.0001) and objective estimate of ventilation (r=0.45, p=0.0157). In addition, both scans were well correlated with reference methods for the diagnosis of emphysema (pulmonary function test and HRCT). The defect scores were largest on 81mKr SPECT (the score on HP 3He MRI was one-third less than that on 81mKr SPECT), but the difference was reduced after normalisation for different breathing depths (HP 3He MRI at total lung capacity; 81mKr SPECT at tidal breathing at functional residual capacity).

CONCLUSION

HP 3He MRI provides detailed ventilation distribution images and defect scores are comparable on HP 3He MRI and 81mKr SPECT. Additionally, new insights into the regional pulmonary microstructure via the apparent diffusion coefficient measurements are provided by HP 3He MRI. HP 3He MRI is a promising new diagnostic tool for the assessment of ventilation distribution.

Cortical deactivation induced by visual stimulation in human slow-wave sleep

It has previously been demonstrated that sleeping and sedated young children respond with a paradoxical decrease in the blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signal in the rostro-medial occipital visual cortex during visual stimulation. It is unresolved whether this negative BOLD response pattern is of developmental neurobiological origin particular to a given age or to a general effect of sleep or sedative drugs. To further elucidate this issue, we used fMRI and positron emission tomography (PET) to study the brain activation pattern during visual stimulation in spontaneously sleeping adult volunteers. In five sleeping volunteers fMRI studies confirmed a robust signal decrease during stimulation in the rostro-medial occipital cortex. A similar relative decrease at the same location was found during visual stimulation and polysomnographically verified slow-wave sleep in a separate group of six subjects using H(2)(15)O PET measures of the regional cerebral blood flow (rCBF). This decrease was more rostro-dorsal compared to the relative rCBF increase along the calcarine sulcus found during visual stimulation in the awake state. This study reconfirms the previously described paradoxical stimulation-correlated negative BOLD signal change in the rostro-medial occipital cortex, expanding this response mode to an age spectrum ranging from the newborn to the adult. Further, the use of complementary brain mapping techniques suggests that this decrease was secondary to a relative rCBF decrease. Possible mechanisms for the paradoxical response pattern during sleep include an active inhibition of the visual cortex or a disruption of an energy-consuming process.

Perfusion quantification using Gaussian process deconvolution

The quantification of perfusion using dynamic susceptibility contrast MRI (DSC-MRI) requires deconvolution to obtain the residual impulse response function (IRF). In this work, a method using the Gaussian process for deconvolution (GPD) is proposed. The fact that the IRF is smooth is incorporated as a constraint in the method. The GPD method, which automatically estimates the noise level in each voxel, has the advantage that model parameters are optimized automatically. The GPD is compared to singular value decomposition (SVD) using a common threshold for the singular values, and to SVD using a threshold optimized according to the noise level in each voxel. The comparison is carried out using artificial data as well as data from healthy volunteers. It is shown that GPD is comparable to SVD with a variable optimized threshold when determining the maximum of the IRF, which is directly related to the perfusion. GPD provides a better estimate of the entire IRF. As the signal-to-noise ratio (SNR) increases or the time resolution of the measurements increases, GPD is shown to be superior to SVD. This is also found for large distribution volumes.

Evaluation of CA1 damage using single-voxel 1H-MRS and un-biased stereology: Can non-invasive measures of N-acetyl-asparate following global ischemia be used as a reliable measure of neuronal damage?

Global brain ischemia provoked by transient occlusion of the carotid arteries (2VO) in gerbils results in a severe loss of neurons in the hippocampal CA1 region. We measured the concentration of the neuron specific N-acetyl-aspartate, [NAA], in the gerbil dorsal hippocampus by proton MR spectroscopy (1H-MRS) in situ, and HPLC, 4 days after global ischemia. The [NAA] was correlated with graded hippocampus damage scoring and stereologically determined neuronal density. A basal hippocampal [NAA] of 8.37+/-0.10 and 9.81+/-0.44 mmol/l were found from HPLC and 1H-MRS, respectively. HPLC measurements of [NAA] obtained from hippocampus 4 days after 2VO showed a 20% reduction in the [NAA] following 4 min of ischemia (P<0.001). 1H-MRS measurements on gerbils subjected to 4 or 8 min of ischemia showed a similar 24% decline in the [NAA] (P<0.05). Thus, there was correlation between the HPLC and 1H-MRS determined NAA decline. There was also a significant correlation between 1H-MRS [NAA] and the corresponding reduction in CA1 neuronal density (P<0.004). In summary our findings show that single voxel 1H-MRS can be used as a supplement to histological evaluation of neuronal injury in studies after global brain ischemia. Accordingly, volume selective spectroscopy has a potential for assessment of neuroprotective therapeutic compounds/strategies with respect to neuronal rescue for delayed ischemic brain damage.

Reconstruction strategy for echo planar spectroscopy and its application to partially undersampled imaging

The most commonly encountered form of echo planar spectroscopy involves oscillating gradients in one spatial dimension during readout. Data are consequently not sampled on a Cartesian grid. A fast gridding algorithm applicable to this particular situation is presented. The method is optimal, i.e., it performs as well as the full discrete Fourier transform for band limited signals while allowing for use of the fast Fourier transform. The method is demonstrated for reconstruction of data that are partially undersampled in the time domain. The advantages of undersampling are lower hardware requirements or fewer interleaves per acquisition. The method is of particular interest when large bandwidths are needed (e.g., for high field scanning) and for scanners with limited gradient performance. The unavoidable artifacts resulting from undersampling are demonstrated to be acceptable for spectroscopy with long echo times.

Optimal voxel size for measuring global gray and white matter proton metabolite concentrations using chemical shift imaging

Quantification of gray and white matter levels of spectroscopically visible metabolites can provide important insights into brain development and pathological conditions. Chemical shift imaging offers a gain in efficiency for estimation of global gray and white matter metabolite concentrations compared to single voxel methods. In the present study, the optimal voxel size is calculated from segmented human brain data and accompanying field maps. The optimal voxel size is found to be approximately 8 cc, but a wide range of values, 4-64 cc, can be chosen with little increase in estimated concentration error (<15%). Magn Reson Med 44:10-18, 2000.