Lipid suppression in CSI with spatial priors and highly undersampled peripheral k-space

Mapping 1H brain metabolites using chemical shift imaging is hampered by the presence of subcutaneous lipid signals, which contaminate the metabolites by ringing due to limited spatial resolution. Even though chemical shift imaging at spatial resolution high enough to mitigate the lipid artifacts is infeasible due to signal-to-noise constraints on the metabolites, the lipid signals have orders of magnitude of higher concentration, which enables the collection of high-resolution lipid maps with adequate signal-to-noise. The previously proposed dual-density approach exploits this high signal-to-noise property of the lipid layer to suppress truncation artifacts using high-resolution lipid maps. Another recent approach for lipid suppression makes use of the fact that metabolite and lipid spectra are approximately orthogonal, and seeks sparse metabolite spectra when projected onto lipid-basis functions. This work combines and extends the dual-density approach and the lipid-basis penalty, while estimating the high-resolution lipid image from 2-average k-space data to incur minimal increase on the scan time. Further, we exploit the spectral-spatial sparsity of the lipid ring and propose to estimate it from substantially undersampled (acceleration R=10 in the peripheral k-space) 2-average in vivo data using compressed sensing and still obtain improved lipid suppression relative to using dual-density or lipid-basis penalty alone.

Local SAR in parallel transmission pulse design

The management of local and global power deposition in human subjects (specific absorption rate, SAR) is a fundamental constraint to the application of parallel transmission (pTx) systems. Even though the pTx and single channel have to meet the same SAR requirements, the complex behavior of the spatial distribution of local SAR for transmission arrays poses problems that are not encountered in conventional single-channel systems and places additional requirements on pTx radio frequency pulse design. We propose a pTx pulse design method which builds on recent work to capture the spatial distribution of local SAR in numerical tissue models in a compressed parameterization in order to incorporate local SAR constraints within computation times that accommodate pTx pulse design during an in vivo magnetic resonance imaging scan. Additionally, the algorithm yields a protocol-specific ultimate peak in local SAR, which is shown to bound the achievable peak local SAR for a given excitation profile fidelity. The performance of the approach was demonstrated using a numerical human head model and a 7 Tesla eight-channel transmit array. The method reduced peak local 10 g SAR by 14-66% for slice-selective pTx excitations and 2D selective pTx excitations compared to a pTx pulse design constrained only by global SAR. The primary tradeoff incurred for reducing peak local SAR was an increase in global SAR, up to 34% for the evaluated examples, which is favorable in cases where local SAR constraints dominate the pulse applications.

Correlation chemical shift imaging with low-power adiabatic pulses and constant-density spiral trajectories

In this work we introduce the concept of correlation chemical shift imaging (CCSI). Novel CCSI pulse sequences are demonstrated on clinical scanners for two-dimensional Correlation Spectroscopy (COSY) and Total Correlation Spectroscopy (TOCSY) imaging experiments. To date there has been limited progress reported towards a feasible and robust multivoxel 2D COSY. Localized 2D TOCSY imaging is shown for the first time in this work. Excitation with adiabatic GOIA-W(16,4) pulses (Gradient Offset Independent Adiabaticity Wurst modulation) provides minimal chemical shift displacement error, reduced lipid contamination from subcutaneous fat, uniform optimal flip angles, and efficient mixing for coupled spins, while enabling short repetition times due to low power requirements. Constant-density spiral readout trajectories are used to acquire simultaneously two spatial dimensions and f(2) frequency dimension in (k(x),k(y),t(2)) space in order to speed up data collection, while f(1) frequency dimension is encoded by consecutive time increments of t(1) in (k(x),k(y),t(1),t(2)) space. The efficient spiral sampling of the k-space enables the acquisition of a single-slice 2D COSY dataset with an 8 × 8 matrix in 8:32  min on 3 T clinical scanners, which makes it feasible for in vivo studies on human subjects. Here we present the first results obtained on phantoms, human volunteers and patients with brain tumors. The patient data obtained by us represent the first clinical demonstration of a feasible and robust multivoxel 2D COSY. Compared to the 2D J-resolved method, 2D COSY and TOCSY provide increased spectral dispersion which scales up with increasing main magnetic field strength and may have improved ability to unambiguously identify overlapping metabolites. It is expected that the new developments presented in this work will facilitate in vivo application of 2D chemical shift correlation MRS in basic science and clinical studies.

Denoising sparse images from GRAPPA using the nullspace method

To accelerate magnetic resonance imaging using uniformly undersampled (nonrandom) parallel imaging beyond what is achievable with generalized autocalibrating partially parallel acquisitions (GRAPPA) alone, the DEnoising of Sparse Images from GRAPPA using the Nullspace method is developed. The trade-off between denoising and smoothing the GRAPPA solution is studied for different levels of acceleration. Several brain images reconstructed from uniformly undersampled k-space data using DEnoising of Sparse Images from GRAPPA using the Nullspace method are compared against reconstructions using existing methods in terms of difference images (a qualitative measure), peak-signal-to-noise ratio, and noise amplification (g-factors) as measured using the pseudo-multiple replica method. Effects of smoothing, including contrast loss, are studied in synthetic phantom data. In the experiments presented, the contrast loss and spatial resolution are competitive with existing methods. Results for several brain images demonstrate significant improvements over GRAPPA at high acceleration factors in denoising performance with limited blurring or smoothing artifacts. In addition, the measured g-factors suggest that DEnoising of Sparse Images from GRAPPA using the Nullspace method mitigates noise amplification better than both GRAPPA and L1 iterative self-consistent parallel imaging reconstruction (the latter limited here by uniform undersampling).

MRI estimates of brain iron concentration in normal aging using quantitative susceptibility mapping

Quantifying tissue iron concentration in vivo is instrumental for understanding the role of iron in physiology and in neurological diseases associated with abnormal iron distribution. Herein, we use recently-developed Quantitative Susceptibility Mapping (QSM) methodology to estimate the tissue magnetic susceptibility based on MRI signal phase. To investigate the effect of different regularization choices, we implement and compare ℓ1 and ℓ2 norm regularized QSM algorithms. These regularized approaches solve for the underlying magnetic susceptibility distribution, a sensitive measure of the tissue iron concentration, that gives rise to the observed signal phase. Regularized QSM methodology also involves a pre-processing step that removes, by dipole fitting, unwanted background phase effects due to bulk susceptibility variations between air and tissue and requires data acquisition only at a single field strength. For validation, performances of the two QSM methods were measured against published estimates of regional brain iron from postmortem and in vivo data. The in vivo comparison was based on data previously acquired using Field-Dependent Relaxation Rate Increase (FDRI), an estimate of MRI relaxivity enhancement due to increased main magnetic field strength, requiring data acquired at two different field strengths. The QSM analysis was based on susceptibility-weighted images acquired at 1.5 T, whereas FDRI analysis used Multi-Shot Echo-Planar Spin Echo images collected at 1.5 T and 3.0 T. Both datasets were collected in the same healthy young and elderly adults. The in vivo estimates of regional iron concentration comported well with published postmortem measurements; both QSM approaches yielded the same rank ordering of iron concentration by brain structure, with the lowest in white matter and the highest in globus pallidus. Further validation was provided by comparison of the in vivo measurements, ℓ1-regularized QSM versus FDRI and ℓ2-regularized QSM versus FDRI, which again yielded perfect rank ordering of iron by brain structure. The final means of validation was to assess how well each in vivo method detected known age-related differences in regional iron concentrations measured in the same young and elderly healthy adults. Both QSM methods and FDRI were consistent in identifying higher iron concentrations in striatal and brain stem ROIs (i.e., caudate nucleus, putamen, globus pallidus, red nucleus, and substantia nigra) in the older than in the young group. The two QSM methods appeared more sensitive in detecting age differences in brain stem structures as they revealed differences of much higher statistical significance between the young and elderly groups than did FDRI. However, QSM values are influenced by factors such as the myelin content, whereas FDRI is a more specific indicator of iron content. Hence, FDRI demonstrated higher specificity to iron yet yielded noisier data despite longer scan times and lower spatial resolution than QSM. The robustness, practicality, and demonstrated ability of predicting the change in iron deposition in adult aging suggest that regularized QSM algorithms using single-field-strength data are possible alternatives to tissue iron estimation requiring two field strengths.

Phase-based regional oxygen metabolism (PROM) using MRI

Venous oxygen saturation (Y(v) ) in cerebral veins and the cerebral metabolic rate of oxygen (CMRO(2)) are important indicators for brain function and disease. Although MRI has been used for global measurements of these parameters, currently there is no recognized technique to quantify regional Y(v) and CMRO(2) using noninvasive imaging. This article proposes a technique to quantify CMRO(2) from independent MRI estimates of Y(v) and cerebral blood flow. The approach uses standard gradient-echo and arterial spin labeling acquisitions to make these measurements. Using MR susceptometry on gradient-echo phase images, Y(v) was quantified for candidate vein segments in gray matter that approximate a long cylinder parallel to the main magnetic field. Local cerebral blood flow for the identified vessel was determined from a corresponding region in the arterial spin labeling perfusion map. Fick's principle of arteriovenous difference was then used to quantify CMRO(2) locally around each vessel. Application of this method in young, healthy subjects provided gray matter averages of 59.6% ± 2.3% for Y(v), 51.7 ± 6.4 mL/100 g/min for cerebral blood flow, and 158 ± 18 μmol/100 g/min for CMRO(2) (mean ± SD, n = 12), which is consistent with values previously reported by positron emission tomography and MRI.

QUantitative Imaging of eXtraction of oxygen and TIssue consumption (QUIXOTIC) using venular-targeted velocity-selective spin labeling

While oxygen extraction fraction (OEF) and cerebral metabolic rate of oxygen (CMRO(2)) are fundamental parameters of brain health and function, a robust MRI-based mapping of OEF and CMRO(2) amenable to functional MRI (fMRI) has not been established. To address this issue, a novel method called QUantitative Imaging of eXtraction of Oxygen and TIssue Consumption, or QUIXOTIC, is introduced. The key innovation in QUIXOTIC is the use of velocity-selective spin labeling to isolate MR signal exclusively from postcapillary venular blood on a voxel-by-voxel basis. Measuring the T(2) of this venular-targeted blood allows calibration to venular oxygen saturation (Y(v)) via theoretical and experimental T(2) versus blood oxygen saturation relationships. Y(v) is converted to OEF, and baseline CMRO(2) is subsequently estimated from OEF and additional cerebral blood flow and hematocrit measurements. Theory behind the QUIXOTIC technique is presented, and implications of cutoff velocity (V(CUTOFF)) and outflow time parameters are discussed. Cortical gray matter values obtained with QUIXOTIC in 10 healthy volunteers are Y(v) = 0.73 ± 0.02, OEF = 0.26 ± 0.02, and CMRO(2) = 125 ± 15 μmol/100 g min. Results are compared to global measures obtained with the T(2) relaxation under spin tagging (TRUST) technique. The preliminary data presented suggest that QUIXOTIC will be useful for mapping Y(v), OEF, and CMRO(2), in both clinical and functional MRI settings.

Multi-contrast reconstruction with Bayesian compressed sensing

Clinical imaging with structural MRI routinely relies on multiple acquisitions of the same region of interest under several different contrast preparations. This work presents a reconstruction algorithm based on Bayesian compressed sensing to jointly reconstruct a set of images from undersampled k-space data with higher fidelity than when the images are reconstructed either individually or jointly by a previously proposed algorithm, M-FOCUSS. The joint inference problem is formulated in a hierarchical Bayesian setting, wherein solving each of the inverse problems corresponds to finding the parameters (here, image gradient coefficients) associated with each of the images. The variance of image gradients across contrasts for a single volumetric spatial position is a single hyperparameter. All of the images from the same anatomical region, but with different contrast properties, contribute to the estimation of the hyperparameters, and once they are found, the k-space data belonging to each image are used independently to infer the image gradients. Thus, commonality of image spatial structure across contrasts is exploited without the problematic assumption of correlation across contrasts. Examples demonstrate improved reconstruction quality (up to a factor of 4 in root-mean-square error) compared with previous compressed sensing algorithms and show the benefit of joint inversion under a hierarchical Bayesian model.

Simulating Magnetic Nanoparticle Behavior in Low-field MRI under Transverse Rotating Fields and Imposed Fluid Flow

In the presence of alternating-sinusoidal or rotating magnetic fields, magnetic nanoparticles will act to realign their magnetic moment with the applied magnetic field. The realignment is characterized by the nanoparticle's time constant, τ. As the magnetic field frequency is increased, the nanoparticle's magnetic moment lags the applied magnetic field at a constant angle for a given frequency, Ω, in rad/s. Associated with this misalignment is a power dissipation that increases the bulk magnetic fluid's temperature which has been utilized as a method of magnetic nanoparticle hyperthermia, particularly suited for cancer in low-perfusion tissue (e.g., breast) where temperature increases of between 4°C and 7°C above the ambient in vivo temperature cause tumor hyperthermia. This work examines the rise in the magnetic fluid's temperature in the MRI environment which is characterized by a large DC field, B(0). Theoretical analysis and simulation is used to predict the effect of both alternating-sinusoidal and rotating magnetic fields transverse to B(0). Results are presented for the expected temperature increase in small tumors (~1 cm radius) over an appropriate range of magnetic fluid concentrations (0.002 to 0.01 solid volume fraction) and nanoparticle radii (1 to 10 nm). The results indicate that significant heating can take place, even in low-field MRI systems where magnetic fluid saturation is not significant, with careful The goal of this work is to examine, by means of analysis and simulation, the concept of interactive fluid magnetization using the dynamic behavior of superparamagnetic iron oxide nanoparticle suspensions in the MRI environment. In addition to the usual magnetic fields associated with MRI, a rotating magnetic field is applied transverse to the main B(0) field of the MRI. Additional or modified magnetic fields have been previously proposed for hyperthermia and targeted drug delivery within MRI. Analytical predictions and numerical simulations of the transverse rotating magnetic field in the presence of B(0) are investigated to demonstrate the effect of Ω, the rotating field frequency, and the magnetic field amplitude on the fluid suspension magnetization. The transverse magnetization due to the rotating transverse field shows strong dependence on the characteristic time constant of the fluid suspension, τ. The analysis shows that as the rotating field frequency increases so that Ωτ approaches unity, the transverse fluid magnetization vector is significantly non-aligned with the applied rotating field and the magnetization's magnitude is a strong function of the field frequency. In this frequency range, the fluid's transverse magnetization is controlled by the applied field which is determined by the operator. The phenomenon, which is due to the physical rotation of the magnetic nanoparticles in the suspension, is demonstrated analytically when the nanoparticles are present in high concentrations (1 to 3% solid volume fractions) more typical of hyperthermia rather than in clinical imaging applications, and in low MRI field strengths (such as open MRI systems), where the magnetic nanoparticles are not magnetically saturated. The effect of imposed Poiseuille flow in a planar channel geometry and changing nanoparticle concentration is examined. The work represents the first known attempt to analyze the dynamic behavior of magnetic nanoparticles in the MRI environment including the effects of the magnetic nanoparticle spin-velocity. It is shown that the magnitude of the transverse magnetization is a strong function of the rotating transverse field frequency. Interactive fluid magnetization effects are predicted due to non-uniform fluid magnetization in planar Poiseuille flow with high nanoparticle concentrations.

Heating in the MRI environment due to superparamagnetic fluid suspensions in a rotating magnetic field

In the presence of alternating-sinusoidal or rotating magnetic fields, magnetic nanoparticles will act to realign their magnetic moment with the applied magnetic field. The realignment is characterized by the nanoparticle's time constant, τ. As the magnetic field frequency is increased, the nanoparticle's magnetic moment lags the applied magnetic field at a constant angle for a given frequency, Ω, in rad/s. Associated with this misalignment is a power dissipation that increases the bulk magnetic fluid's temperature which has been utilized as a method of magnetic nanoparticle hyperthermia, particularly suited for cancer in low-perfusion tissue (e.g., breast) where temperature increases of between 4°C and 7°C above the ambient in vivo temperature cause tumor hyperthermia. This work examines the rise in the magnetic fluid's temperature in the MRI environment which is characterized by a large DC field, B(0). Theoretical analysis and simulation is used to predict the effect of both alternating-sinusoidal and rotating magnetic fields transverse to B(0). Results are presented for the expected temperature increase in small tumors (~1 cm radius) over an appropriate range of magnetic fluid concentrations (0.002 to 0.01 solid volume fraction) and nanoparticle radii (1 to 10 nm). The results indicate that significant heating can take place, even in low-field MRI systems where magnetic fluid saturation is not significant, with careful selection of the rotating or sinusoidal field parameters (field frequency and amplitude). The work indicates that it may be feasible to combine low-field MRI with a magnetic hyperthermia system using superparamagnetic iron oxide nanoparticles.

Measuring SPIO and Gd contrast agent magnetization using 3 T MRI

Traditional methods of measuring magnetization in magnetic fluid samples, such as vibrating sample magnetometry (VSM), are typically limited to maximum field strengths of about 1 T. This work demonstrates the ability of MRI to measure the magnetization associated with two commercial MRI contrast agents at 3 T by comparing analytical solutions to experimental imaging results for the field pattern associated with agents in cylindrical vials. The results of the VSM and fitted MRI data match closely. The method represents an improvement over VSM measurements since results are attainable at imaging field strengths. The agents investigated are Feridex, a superparamagnetic iron oxide suspension used primarily for liver imaging, and Magnevist, a paramagnetic, gadolinium-based compound used for tumors, inflammation and vascular lesions. MR imaging of the agents took place in sealed cylindrical vials in the presence of a surrounding volume of deionized water where the effects of the contrast agents had a measurable effect on the water's magnetization in the vicinity of the compartment of contrast agent. A pair of phase images were used to reconstruct a B(0) fieldmap. The resultant B(0) maps in the water region, corrected for shimming and container edge effects, were used to predict the agent's magnetization at 3 T. The results were compared with the results from VSM measurements up to 1.2 T and close correlation was observed. The technique should be of interest to those seeking quantification of the magnetization associated with magnetic suspensions beyond the traditional scope of VSM. The magnetization needs to be sufficiently strong (M(s) >or= 50 Am(2)/kg Fe for Feridex and X(m) >or=5 x 10(-5) m(3)/kg Gd for Magnevist) for a measurable dipole field in the surrounding water. For this reason, the technique is mostly suitable for undiluted agents.

Broadband slab selection with B1+ mitigation at 7T via parallel spectral-spatial excitation

Chemical shift imaging benefits from signal-to-noise ratio (SNR) and chemical shift dispersion increases at stronger main field such as 7 Tesla, but the associated shorter radiofrequency (RF) wavelengths encountered require B1+ mitigation over both the spatial field of view (FOV) and a specified spectral bandwidth. The bandwidth constraint presents a challenge for previously proposed spatially tailored B1+ mitigation methods, which are based on a type of echovolumnar trajectory referred to as "spokes" or "fast-kz". Although such pulses, in conjunction with parallel excitation methodology, can efficiently mitigate large B1+ inhomogeneities and achieve relatively short pulse durations with slice-selective excitations, they exhibit a narrow-band off-resonance response and may not be suitable for applications that require B1+ mitigation over a large spectral bandwidth. This work outlines a design method for a general parallel spectral-spatial excitation that achieves a target-error minimization simultaneously over a bandwidth of frequencies and a specified spatial-domain. The technique is demonstrated for slab-selective excitation with in-plane B1+ mitigation over a 600-Hz bandwidth. The pulse design method is validated in a water phantom at 7T using an eight-channel transmit array system. The results show significant increases in the pulse's spectral bandwidth, with no additional pulse duration penalty and only a minor tradeoff in spatial B1+ mitigation compared to the standard spoke-based parallel RF design.

Seven-Tesla proton magnetic resonance spectroscopic imaging in adult X-linked adrenoleukodystrophy

BACKGROUND

Adults with X-linked adrenoleukodystrophy (X-ALD) remain at risk for progressive neurological deterioration. Phenotypes vary in their pathology, ranging from axonal degeneration to inflammatory demyelination. The severity of symptoms is poorly explained by conventional imaging.

OBJECTIVE

To test the hypothesis that neurochemistry in normal-appearing brains differs in adult phenotypes of X-ALD and that neurochemical changes correlate with the severity of symptoms.

PATIENTS AND METHODS

Using a 7-Tesla scanner, we performed structural and proton magnetic resonance spectroscopic imaging in 13 adult patients with X-ALD: 4 patients with adult cerebral ALD, 5 patients with adrenomyeloneuropathy (AMN), and 4 female heterozygotes. Nine healthy controls were included.

RESULTS

Among adult X-ALD phenotypes, the myo-inositol to creatine ratio was 46% higher and the choline to creatine ratio was 21% higher in normal-appearing white matter of those with adult cerebral ALD compared with those with AMN (P < .05). Both N-acetylaspartate to creatine (P = .03) and glutamate to creatine (P = .04) ratios were lower in AMN patients than in controls. There were no significant differences between patients with AMN and female heterozygotes. In the cortex, patients with adult cerebral ALD had lower N-acetylaspartate to creatine ratios compared with female heterozygotes and controls (P = .02). The global myo-inositol to creatine ratio demonstrated a significant association with Expanded Disability Status Scale score (Spearman rho = 0.66, P = .04).

CONCLUSIONS

Seven-Tesla proton magnetic resonance spectroscopic imaging reveals differences in the neurochemistry of adult cerebral ALD but cannot distinguish AMN patients from female heterozygotes. Myo-inositol to creatine ratio correlates with the severity of the symptoms and may be a meaningful biomarker in adult X-ALD.

Slice-selective RF pulses for in vivo B1+ inhomogeneity mitigation at 7 tesla using parallel RF excitation with a 16-element coil

Slice-selective RF waveforms that mitigate severe B1+ inhomogeneity at 7 Tesla using parallel excitation were designed and validated in a water phantom and human studies on six subjects using a 16-element degenerate stripline array coil driven with a butler matrix to utilize the eight most favorable birdcage modes. The parallel RF waveform design applied magnitude least-squares (MLS) criteria with an optimized k-space excitation trajectory to significantly improve profile uniformity compared to conventional least-squares (LS) designs. Parallel excitation RF pulses designed to excite a uniform in-plane flip angle (FA) with slice selection in the z-direction were demonstrated and compared with conventional sinc-pulse excitation and RF shimming. In all cases, the parallel RF excitation significantly mitigated the effects of inhomogeneous B1+ on the excitation FA. The optimized parallel RF pulses for human B1+ mitigation were only 67% longer than a conventional sinc-based excitation, but significantly outperformed RF shimming. For example the standard deviations (SDs) of the in-plane FA (averaged over six human studies) were 16.7% for conventional sinc excitation, 13.3% for RF shimming, and 7.6% for parallel excitation. This work demonstrates that excitations with parallel RF systems can provide slice selection with spatially uniform FAs at high field strengths with only a small pulse-duration penalty.

High-flip-angle slice-selective parallel RF transmission with 8 channels at 7 T

At high magnetic field, B(1)(+) non-uniformity causes undesired inhomogeneity in SNR and image contrast. Parallel RF transmission using tailored 3D k-space trajectory design has been shown to correct for this problem and produce highly uniform in-plane magnetization with good slice selection profile within a relatively short excitation duration. However, at large flip angles the excitation k-space based design method fails. Consequently, several large-flip-angle parallel transmission designs have recently been suggested. In this work, we propose and demonstrate a large-flip-angle parallel excitation design for 90 degrees and 180 degrees spin-echo slice-selective excitations that mitigate severe B(1)(+) inhomogeneity. The method was validated on an 8-channel transmit array at 7T using a water phantom with B(1)(+) inhomogeneity similar to that seen in human brain in vivo. Slice-selective excitations with parallel RF systems offer means to implement conventional high-flip excitation sequences without a severe pulse-duration penalty, even at very high B(0) field strengths where large B(1)(+) inhomogeneity is present.

Specific absorption rate studies of the parallel transmission of inner-volume excitations at 7T

PURPOSE

To investigate the behavior of whole-head and local specific absorption rate (SAR) as a function of trajectory acceleration factor and target excitation pattern due to the parallel transmission (pTX) of spatially tailored excitations at 7T.

MATERIALS AND METHODS

Finite-difference time domain (FDTD) simulations in a multitissue head model were used to obtain B(1) (+) and electric field maps of an eight-channel transmit head array. Local and average SAR produced by 2D-spiral-trajectory excitations were examined as a function of trajectory acceleration factor, R, and a variety of target excitation parameters when pTX pulses are designed for constant root-mean-square excitation pattern error.

RESULTS

Mean and local SAR grow quadratically with flip angle and more than quadratically with R, but the ratio of local to mean SAR is not monotonic with R. SAR varies greatly with target position, exhibiting different behaviors as a function of target shape and size for small and large R. For example, exciting large regions produces less SAR than exciting small ones for R >or=4, but the opposite trend occurs when R <4. Furthermore, smoother and symmetric patterns produce lower SAR.

CONCLUSION

Mean and local SAR vary by orders of magnitude depending on acceleration factor and excitation pattern, often exhibiting complex, nonintuitive behavior. To ensure safety compliance, it seems that model-based validation of individual target patterns and corresponding pTX pulses is necessary.

Sparsity-enforced slice-selective MRI RF excitation pulse design

We introduce a novel algorithm for the design of fast slice-selective spatially-tailored magnetic resonance imaging (MRI) excitation pulses. This method, based on sparse approximation theory, uses a second-order cone optimization to place and modulate a small number of slice-selective sinc-like radio-frequency (RF) pulse segments ("spokes") in excitation k-space, enforcing sparsity on the number of spokes allowed while simultaneously encouraging those that remain to be placed and modulated in a way that best forms a user-defined in-plane target magnetization. Pulses are designed to mitigate B(1) inhomogeneity in a water phantom at 7 T and to produce highly-structured excitations in an oil phantom on an eight-channel parallel excitation system at 3 T. In each experiment, pulses generated by the sparsity-enforced method outperform those created via conventional Fourier-based techniques, e.g., when attempting to produce a uniform magnetization in the presence of severe B(1) inhomogeneity, a 5.7-ms 15-spoke pulse generated by the sparsity-enforced method produces an excitation with 1.28 times lower root mean square error than conventionally-designed 15-spoke pulses. To achieve this same level of uniformity, the conventional methods need to use 29-spoke pulses that are 7.8 ms long.

Fast slice-selective radio-frequency excitation pulses for mitigating B+1 inhomogeneity in the human brain at 7 Tesla

A novel radio-frequency (RF) pulse design algorithm is presented that generates fast slice-selective excitation pulses that mitigate B+1 inhomogeneity present in the human brain at high field. The method is provided an estimate of the B+1 field in an axial slice of the brain and then optimizes the placement of sinc-like "spokes" in kz via an L1-norm penalty on candidate (kx, ky) locations; an RF pulse and gradients are then designed based on these weighted points. Mitigation pulses are designed and demonstrated at 7T in a head-shaped water phantom and the brain; in each case, the pulses mitigate a significantly nonuniform transmit profile and produce nearly uniform flip angles across the field of excitation (FOX). The main contribution of this work, the sparsity-enforced spoke placement and pulse design algorithm, is derived for conventional single-channel excitation systems and applied in the brain at 7T, but readily extends to lower field systems, nonbrain applications, and multichannel parallel excitation arrays.

Magnitude least squares optimization for parallel radio frequency excitation design demonstrated at 7 Tesla with eight channels

Spatially tailored radio frequency (RF) excitations accelerated with parallel transmit systems provide the opportunity to create shaped volume excitations or mitigate inhomogeneous B(1) excitation profiles with clinically relevant pulse lengths. While such excitations are often designed as a least-squares optimized approximation to a target magnitude and phase profile, adherence to the target phase profile is usually not important as long as the excitation phase is slowly varying compared with the voxel dimension. In this work, we demonstrate a method for a magnitude least squares optimization of the target magnetization profile for multichannel parallel excitation to improve the magnitude profile and reduce the RF power at the cost of a less uniform phase profile. The method enables the designer to trade off the allowed spatial phase variation for the improvement in magnitude profile and reduction in RF power. We validate the method with simulation studies and demonstrate its performance in fourfold accelerated two-dimensional spiral excitations, as well as for uniform in-plane slice selective parallel excitations using an eight-channel transmit array on a 7T human MRI scanner. The experimental results are in good agreement with the simulations, which show significant improvement in the magnitude profile and reductions in the required RF power while still maintaining negligible intravoxel phase variation.

Development and resolution of brain lesions caused by pyrithiamine- and dietary-induced thiamine deficiency and alcohol exposure in the alcohol-preferring rat: a longitudinal magnetic resonance imaging and spectroscopy study

Wernicke's encephalopathy (WE) is characterized by lesions in thalamus, hypothalamus (including mammillary nuclei), and inferior colliculi, results in serious disabilities, has an etiology of thiamine deficiency, is treatable with thiamine, and occurs most commonly with alcoholism. Despite decades of study, whether alcohol exposure exacerbates the neuropathology or retards its resolution remains controversial. To examine patterns of brain damage and recovery resulting from thiamine deprivation with and without alcohol exposure, we conducted in vivo magnetic resonance imaging (MRI) and magnetic resonance spectroscopy (MRS) at 3 T in alcohol-preferring (P) rats, which had voluntarily consumed large amounts of alcohol before thiamine manipulation. A total of 18 adult male P rats (nine alcohol-exposed) received a thiamine-deficient diet for 2 weeks: 10 (five alcohol-exposed) received intraperitoneal (i.p.) pyrithiamine (PT) and eight (four alcohol-exposed) received i.p. thiamine supplementation. Neurological signs developed by day 14. Rats were scanned before thiamine depletion and 18 and 35 days after thiamine repletion. Two-dimensional J-resolved MRS single-voxel spectra with water reference were collected in a voxel subtending the thalamus; metabolite quantification was corrected for voxel tissue content. MRI identified significant enlargement of dorsal ventricles and increase in signal intensities in thalamus, inferior colliculi, and mammillary nuclei of PT compared with thiamine-treated (TT) groups from MRI 1-2, followed by significant normalization from MRI 2-3 in thalamus and colliculi, but not mammillary nuclei and lateral ventricles. Voxel-by-voxel analysis revealed additional hyperintense signal clusters in the dorsal and ventral hippocampus and enlargement of the fourth ventricle. MRS showed a significant decline and then partial recovery in thalamic N-acetylaspartate, a marker of neuronal integrity, in PT compared with TT rats, with no change detected in creatine, choline, or glutamate. PT rats with prior alcohol exposure exhibited attenuated recovery in the thalamus and arrested growth of the corpus callosum; further, two of the five alcohol-exposed PT rats died prematurely. Parenchymal and ventricular changes with thiamine manipulation concur with human radiological signs of WE. The enduring macrostructural and neurochemical abnormalities involving critical nodes of Papez circuit carry liabilities for development of amnesia and incomplete recovery from other cognitive and motor functions subserved by the affected neural systems.

In vivo fiber tracking in the rat brain on a clinical 3T MRI system using a high strength insert gradient coil

In vivo neuroimaging methods permit longitudinal quantitative examination of the dynamic course of neurodegenerative conditions in humans and animal models and enable assessment of therapeutic efforts in mitigating disease effects on brain systems. The study of conditions affecting white matter, such as multiple sclerosis, demyelinating conditions, and drug and alcohol dependence, can be accomplished with diffusion tensor imaging (DTI), a technique uniquely capable of probing the microstructural integrity of white matter fibers in the living brain. We used a 3T clinical MR scanner equipped with an insert gradient coil that yields an order of magnitude increase in performance over the whole-body hardware to acquire in vivo DTI images of rat brain. The resolution allowed for fiber tracking evaluation of fractional anisotropy (FA) and apparent diffusion coefficients in the genu and splenium of the corpus callosum. A comparison of short (46 min) and long (92 min) acquisition time DTI protocols indicated low but adequate signal-to-noise ratio (SNR=6.2) of the shorter protocol to conduct quantitative fiber tracking enhanced by multiple acquisitions. As observed in human studies, FA in the rat splenium was higher than in the genu. Advantages of this technology include the use of similar user interface, pulse sequences, and field strength for preclinical animal and clinical human research, enhancing translational capabilities. An additional benefit of scanning at lower field strength, such as 3 T, is the reduction of artifacts due to main field inhomogeneity relative to higher field animal systems.

Degenerate mode band-pass birdcage coil for accelerated parallel excitation

An eight-rung, 3T degenerate birdcage coil (DBC) was constructed and evaluated for accelerated parallel excitation of the head with eight independent excitation channels. Two mode configurations were tested. In the first, each of the eight loops formed by the birdcage was individually excited, producing an excitation pattern similar to a loop coil array. In the second configuration a Butler matrix transformed this "loop coil" basis set into a basis set representing the orthogonal modes of the birdcage coil. In this case the rung currents vary sinusoidally around the coil and only four of the eight modes have significant excitation capability (the other four produce anticircularly polarized (ACP) fields). The lowest useful mode produces the familiar uniform B(1) field pattern, and the higher-order modes produce center magnitude nulls and azimuthal phase variations. The measured magnitude and phase excitation profiles of the individual modes were used to generate one-, four-, six-, and eightfold-accelerated spatially tailored RF excitations with 2D and 3D k-space excitation trajectories. Transmit accelerations of up to six-fold were possible with acceptable levels of spatial artifact. The orthogonal basis set provided by the Butler matrix was found to be advantageous when an orthogonal subset of these modes was used to mitigate B(1) transmit inhomogeneities using parallel excitation.

In vivo quantification of ethanol kinetics in rat brain

Proton magnetic resonance spectroscopy was used at 3T to measure the uptake and clearance of brain ethanol in rats after bolus intraperitoneal (i.p.) or intragastric (i.g.) alcohol injection, and to estimate the effects of acute alcohol on brain metabolites. The observation duration was 1-1.5 h with temporal resolution of alcohol sampling ranging from 4 s-4 min. The observed time course of alcohol brain concentration followed a consistent pattern characterized by a rapid absorption, an intermediate distribution, and a slower clearance that approached a linear decay. In a sample of eight healthy Wistar rats, the intercept of the linear clearance term, extrapolated back to the time of injection, correlated well with the administered dose per unit of lean body mass. Alcohol concentration estimation based on spectroscopically measured clearance was compared with blood alcohol levels from blood samples at the end of observation, and were in good agreement with the administered dose. Serial proton spectroscopy measurements provide a valid in vivo method for quantifying brain alcohol uptake and elimination kinetics in real time.

Diffusion tensor imaging with quantitative fibre tracking in HIV infection and alcoholism comorbidity: synergistic white matter damage

A substantial proportion of individuals infected with human immunodeficiency virus (HIV) also abuse alcohol. Given that each condition can disrupt brain structural integrity, with a predilection for white matter, we used MR diffusion tensor imaging (DTI) and quantitative fibre tracking to examine the separate and combined effects on the microstructure of the corpus callosum. Subjects were men and women with alcoholism alone (n = 87), HIV infection alone (n = 42), alcoholism and HIV infection comorbidity (n = 52) and non-affected controls (n = 88). The two alcoholism groups had similar lifetime alcohol consumption histories; the two HIV-infected groups had similar CD4+ counts and viral loads; all groups were matched in body mass index, and no participant was demented. Compared with controls, all patient groups had lower fractional anisotropy (FA) and higher mean diffusivity (MD) in callosal regions and fibre bundles coursing through the genu and splenium, but these effects were only significant in the two groups with alcoholism, which exhibited 0.65-1.2 SD abnormalities in FA and MD. The callosal regions were differentially affected by alcoholism, with the genu more affected than the splenium, a pattern even more pronounced in the fibre tracks. When the HIV-infected groups were divided by disease severity defined as an acquired immunodeficiency syndrome (AIDS)-defining event or low CD4+ counts (<200) and alcoholism comorbidity, the HIV-infected subgroup with AIDS and alcoholism exhibited approximately 2 SD FA and MD abnormalities in the callosal sectors and fibres, abnormalities that were more than twice the effect sizes observed in the other three HIV-infected subgroups. Degradation of the callosal microstructure was consistently associated with alcoholism, with evidence for compounded alcoholism-HIV effects. Functional relevance of the microstructural abnormalities was supported by associations between motor deficits and low FA or high MD within the diagnostic groups. The high prevalence of alcoholism in HIV-infected individuals and the interfering effect of alcohol on HIV pharmacological response and therapy compliance underscore the need to recognize the independent and synergistic contributions of each condition to brain structure and function.

Parallel RF transmission with eight channels at 3 Tesla

Spatially selective RF waveforms were designed and demonstrated for parallel excitation with a dedicated eight-coil transmit array on a modified 3T human MRI scanner. Measured excitation profiles of individual coils in the array were used in a low-flip-angle pulse design to achieve desired spatial target profiles with two- (2D) and three-dimensional (3D) k-space excitation with simultaneous transmission of RF on eight channels. The 2D pulse excited a high-resolution spatial pattern in-plane, while the 3D trajectory produced high-quality slice selection with a uniform in-plane excitation despite the highly nonuniform individual spatial profiles of the coil array. The multichannel parallel RF excitation was used to accelerate the 2D excitation by factors of 2-8, and experimental results were in excellent agreement with simulations based on the measured coil maps. Parallel RF transmission may become critical for robust and routine human studies at very high field strengths where B(1) inhomogeneity is commonly severe.