Improved reconstruction of crossing fibers in the mouse optic pathways with orientation distribution function fingerprinting

PURPOSE

The accuracy of diffusion MRI tractography reconstruction decreases in the white matter regions with crossing fibers. The optic pathways in rodents provide a challenging structure to test new diffusion tractography approaches because of the small crossing volume within the optic chiasm and the unbalanced 9:1 proportion between the contra- and ipsilateral neural projections from the retina to the lateral geniculate nucleus, respectively.

METHODS

Common approaches based on Orientation Distribution Function (ODF) peak finding or statistical inference were compared qualitatively and quantitatively to ODF Fingerprinting (ODF-FP) for reconstruction of crossing fibers within the optic chiasm using in vivo diffusion MRI (n = 18 $$ n=18 $$ healthy C57BL/6 mice). Manganese-Enhanced MRI (MEMRI) was obtained after intravitreal injection of manganese chloride and used as a reference standard for the optic pathway anatomy.

RESULTS

ODF-FP outperformed by over 100% all the tested methods in terms of the ratios between the contra- and ipsilateral segments of the reconstructed optic pathways as well as the spatial overlap between tractography and MEMRI.

CONCLUSION

In this challenging model system, ODF-Fingerprinting reduced uncertainty of diffusion tractography for complex structural formations of fiber bundles.

Tractography passes the test: Results from the diffusion-simulated connectivity (disco) challenge

Estimating structural connectivity from diffusion-weighted magnetic resonance imaging is a challenging task, partly due to the presence of false-positive connections and the misestimation of connection weights. Building on previous efforts, the MICCAI-CDMRI Diffusion-Simulated Connectivity (DiSCo) challenge was carried out to evaluate state-of-the-art connectivity methods using novel large-scale numerical phantoms. The diffusion signal for the phantoms was obtained from Monte Carlo simulations. The results of the challenge suggest that methods selected by the 14 teams participating in the challenge can provide high correlations between estimated and ground-truth connectivity weights, in complex numerical environments. Additionally, the methods used by the participating teams were able to accurately identify the binary connectivity of the numerical dataset. However, specific false positive and false negative connections were consistently estimated across all methods. Although the challenge dataset doesn't capture the complexity of a real brain, it provided unique data with known macrostructure and microstructure ground-truth properties to facilitate the development of connectivity estimation methods.

Stepwise Stochastic Dictionary Adaptation Improves Microstructure Reconstruction with Orientation Distribution Function Fingerprinting

Fitting of the multicompartment biophysical model of white matter is an ill-posed optimization problem. One approach to make it computationally tractable is through Orientation Distribution Function (ODF) Fingerprinting. However, the accuracy of this method relies solely on ODF dictionary generation mechanisms which either sample the microstructure parameters on a multidimensional grid or draw them randomly with a uniform distribution. In this paper, we propose a stepwise stochastic adaptation mechanism to generate ODF dictionaries tailored specifically to the diffusion-weighted images in hand. The results we obtained on a diffusion phantom and in vivo human brain images show that our reconstructed diffusivities are less noisy and the separation of a free water fraction is more pronounced than for the prior (uniform) distribution of ODF dictionaries.

Performance of orientation distribution function-fingerprinting with a biophysical multicompartment diffusion model

PURPOSE

Orientation Distribution Function (ODF) peak finding methods typically fail to reconstruct fibers crossing at shallow angles below 40°, leading to errors in tractography. ODF-Fingerprinting (ODF-FP) with the biophysical multicompartment diffusion model allows for breaking this barrier.

METHODS

A randomized mechanism to generate a multidimensional ODF-dictionary that covers biologically plausible ranges of intra- and extra-axonal diffusivities and fraction volumes is introduced. This enables ODF-FP to address the high variability of brain tissue. The performance of the proposed approach is evaluated on both numerical simulations and a reconstruction of major fascicles from high- and low-resolution in vivo diffusion images.

RESULTS

ODF-FP with the suggested modifications correctly identifies fibers crossing at angles as shallow as 10 degrees in the simulated data. In vivo, our approach reaches 56% of true positives in determining fiber directions, resulting in visibly more accurate reconstruction of pyramidal tracts, arcuate fasciculus, and optic radiations than the state-of-the-art techniques. Moreover, the estimated diffusivity values and fraction volumes in corpus callosum conform with the values reported in the literature.

CONCLUSION

The modified ODF-FP outperforms commonly used fiber reconstruction methods at shallow angles, which improves deterministic tractography outcomes of major fascicles. In addition, the proposed approach allows for linearization of the microstructure parameters fitting problem.

Free-breathing radial imaging using a pilot-tone radiofrequency transmitter for detection of respiratory motion

PURPOSE

To describe an approach for detection of respiratory signals using a transmitted radiofrequency (RF) reference signal called Pilot-Tone (PT) and to use the PT signal for creation of motion-resolved images based on 3D stack-of-stars imaging under free-breathing conditions.

METHODS

This work explores the use of a reference RF signal generated by a small RF transmitter, placed outside the MR bore. The reference signal is received in parallel to the MR signal during each readout. Because the received PT amplitude is modulated by the subject's breathing pattern, a respiratory signal can be obtained by detecting the strength of the received PT signal over time. The breathing-induced PT signal modulation can then be used for reconstructing motion-resolved images from free-breathing scans. The PT approach was tested in volunteers using a radial stack-of-stars 3D gradient echo (GRE) sequence with golden-angle acquisition.

RESULTS

Respiratory signals derived from the proposed PT method were compared to signals from a respiratory cushion sensor and k-space-center-based self-navigation under different breathing conditions. Moreover, the accuracy was assessed using a modified acquisition scheme replacing the golden-angle scheme by a zero-angle acquisition. Incorporating the PT signal into eXtra-Dimensional (XD) motion-resolved reconstruction led to improved image quality and clearer anatomical depiction of the lung and liver compared to k-space-center signal and motion-averaged reconstruction, when binned into 6, 8, and 10 motion states.

CONCLUSION

PT is a novel concept for tracking respiratory motion. Its small dimension (8 cm), high sampling rate, and minimal interaction with the imaging scan offers great potential for resolving respiratory motion.

Mapping brain-behavior networks using functional and structural connectome fingerprinting in the HCP dataset

INTRODUCTION

Connectome analysis of the human brain's structural and functional architecture provides a unique opportunity to understand the organization of the brain's functional architecture. In previous studies, connectome fingerprinting using brain functional connectivity profiles as an individualized trait was able to predict an individual's neurocognitive performance from the Human Connectome Project (HCP) neurocognitive datasets.

MATERIALS AND METHODS

In the present study, we extend connectome fingerprinting from functional connectivity (FC) to structural connectivity (SC), identifying multiple relationships between behavioral traits and brain connectivity. Higher-order neurocognitive tasks were found to have a weaker association with structural connectivity than its functional connectivity counterparts.

RESULTS

Neurocognitive tasks with a higher sensory footprint were, however, found to have a stronger association with structural connectivity than their functional connectivity counterparts. Language behavioral measurements had a particularly stronger correlation, especially between performance on the picture language test (Pic Vocab) and both FC (r = .28, p < .003) and SC (r = 0.27, p < .00077).

CONCLUSIONS

At the neural level, we found that the pattern of structural brain connectivity related to high-level language performance is consistent with the language white matter regions identified in presurgical mapping. We illustrate how this approach can be used to generalize the connectome fingerprinting framework to structural connectivity and how this can help understand the connections between cognitive behavior and the white matter connectome of the brain.

Using fMRI connectivity to define a treatment-resistant form of post-traumatic stress disorder

A mechanistic understanding of the pathology of psychiatric disorders has been hampered by extensive heterogeneity in biology, symptoms, and behavior within diagnostic categories that are defined subjectively. We investigated whether leveraging individual differences in information-processing impairments in patients with post-traumatic stress disorder (PTSD) could reveal phenotypes within the disorder. We found that a subgroup of patients with PTSD from two independent cohorts displayed both aberrant functional connectivity within the ventral attention network (VAN) as revealed by functional magnetic resonance imaging (fMRI) neuroimaging and impaired verbal memory on a word list learning task. This combined phenotype was not associated with differences in symptoms or comorbidities, but nonetheless could be used to predict a poor response to psychotherapy, the best-validated treatment for PTSD. Using concurrent focal noninvasive transcranial magnetic stimulation and electroencephalography, we then identified alterations in neural signal flow in the VAN that were evoked by direct stimulation of that network. These alterations were associated with individual differences in functional fMRI connectivity within the VAN. Our findings define specific neurobiological mechanisms in a subgroup of patients with PTSD that could contribute to the poor response to psychotherapy.

Noninvasive PET Imaging of CDK4/6 Activation in Breast Cancer

The cell cycle is a progression of 4 distinct phases (G1, S, G2, and M), with various cycle proteins being essential in regulating this process. We aimed to develop a radiolabeled cyclin-dependent kinase 4/6 (CDK4/6) inhibitor for breast cancer imaging. Our transfluorinated analog (¹⁸F-CDKi) was evaluated and validated as a novel PET imaging agent to quantify CDK4/6 expression in estrogen receptor (ER)-positive human epidermal growth factor receptor 2 (HER₂)-negative breast cancer. Methods: ¹⁸F-CDKi was synthesized and assayed against CDK4/6 kinases. ¹⁸F-CDKi was prepared with a 2-step automated synthetic strategy that yielded the final product with remarkable purity and molar activity. In vitro and in vivo biologic specificity was assessed in a MCF-7 cell line and in mice bearing MCF-7 breast tumors. Nonradioactive palbociclib was used as a blocking agent to investigate the binding specificity and selectivity of ¹⁸F-CDKi. Results: ¹⁸F-CDKi was obtained with an overall radiochemical uncorrected yield of 15% and radiochemical purity higher than 98%. The total time from the start of synthesis to the final injectable formulated tracer is 70 min. The retention time reported for ¹⁸F-CDKi and ¹⁹F-CDKi is 27.4 min as demonstrated by coinjection with ¹⁹F-CDKi in a high-pressure liquid chromatograph. In vivo blood half-life (weighted, 7.03 min) and octanol/water phase partition coefficient (1.91 ± 0.24) showed a mainly lipophilic behavior. ¹⁸F-CDKi is stable in vitro and in vivo (>98% at 4 h after injection) and maintained its potent targeting affinity to CDK4/6. Cellular uptake experiments performed on the MCF-7 breast cancer cell line (ER-positive and HER₂-negative) demonstrated specific uptake with a maximum intracellular concentration of about 65% as early as 10 min after incubation. The tracer uptake was reduced to less than 5% when cells were coincubated with a molar excess of palbociclib. In vivo imaging and ex vivo biodistribution of ER-positive, HER₂-negative MCF-7 breast cancer models showed a specific uptake of approximately 4% injected dose/g of tumor (reduced to ∼0.3% with a 50-fold excess of cold palbociclib). A comprehensive biodistribution analysis also revealed a significantly lower activation of CDK4/6 in nontargeting organs. Conclusion: ¹⁸F-CDKi represents the first ¹⁸F PET CDK4/6 imaging agent and a promising imaging agent for ER-positive, HER₂-negative breast cancer.

Fingerprinting Orientation Distribution Functions in diffusion MRI detects smaller crossing angles

Diffusion tractography is routinely used to study white matter architecture and brain connectivity in vivo. A key step for successful tractography of neuronal tracts is the correct identification of tract directions in each voxel. Here we propose a fingerprinting-based methodology to identify these fiber directions in Orientation Distribution Functions, dubbed ODF-Fingerprinting (ODF-FP). In ODF-FP, fiber configurations are selected based on the similarity between measured ODFs and elements in a pre-computed library. In noisy ODFs, the library matching algorithm penalizes the more complex fiber configurations. ODF simulations and analysis of bootstrapped partial and whole-brain in vivo datasets show that the ODF-FP approach improves the detection of fiber pairs with small crossing angles while maintaining fiber direction precision, which leads to better tractography results. Rather than focusing on the ODF maxima, the ODF-FP approach uses the whole ODF shape to infer fiber directions to improve the detection of fiber bundles with small crossing angle. The resulting fiber directions aid tractography algorithms in accurately displaying neuronal tracts and calculating brain connectivity.

Cardiorespiratory motion-tracking via self-refocused rosette navigators

PURPOSE

To develop a flexible method for tracking respiratory and cardiac motions throughout MR and PET-MR body examinations that requires no additional hardware and minimal sequence modification.

METHODS

The incorporation of a contrast-neutral rosette navigator module following the RF excitation allows for robust cardiorespiratory motion tracking with minimal impact on the host sequence. Spatial encoding gradients are applied to the FID signal and the desired motion signals are extracted with a blind source separation technique. This approach is validated with an anthropomorphic, PET-MR-compatible motion phantom as well as in 13 human subjects.

RESULTS

Both respiratory and cardiac motions were reliably extracted from the proposed rosette navigator in phantom and patient studies. In the phantom study, the MR-derived motion signals were additionally validated against the ground truth measurement of diaphragm displacement and left ventricle model triggering pulse.

CONCLUSION

The proposed method yields accurate respiratory and cardiac motion-state tracking, requiring only a short (1.76 ms) additional navigator module, which is self-refocusing and imposes minimal constraints on sequence design.

MRI Evidence of Altered Callosal Sodium in Mild Traumatic Brain Injury

BACKGROUND AND PURPOSE

Mild traumatic brain injury is a leading cause of death and disability worldwide with 42 million cases reported annually, increasing the need to understand the underlying pathophysiology because this could help guide the development of targeted therapy. White matter, particularly the corpus callosum, is susceptible to injury. Animal models suggest stretch-induced mechanoporation of the axonal membrane resulting in ionic shifts and altered sodium ion distribution. The purpose of this study was to compare the distribution of total sodium concentration in the corpus callosum between patients with mild traumatic brain injury and controls using sodium (²³Na) MR imaging.

MATERIALS AND METHODS

Eleven patients with a history of mild traumatic brain injury and 10 age- and sex-matched controls underwent sodium (²³Na) MR imaging using a 3T scanner. Total sodium concentration was measured in the genu, body, and splenium of the corpus callosum with 5-mm ROIs; total sodium concentration of the genu-to-splenium ratio was calculated and compared between patients and controls.

RESULTS

Higher total sodium concentration in the genu (49.28 versus 43.29 mmol/L, P = .01) and lower total sodium concentration in the splenium (which was not statistically significant; 38.35 versus 44.06 mmol/L, P = .08) was seen in patients with mild traumatic brain injury compared with controls. The ratio of genu total sodium concentration to splenium total sodium concentration was also higher in patients with mild traumatic brain injury (1.3 versus 1.01, P = .001).

CONCLUSIONS

Complex differences are seen in callosal total sodium concentration in symptomatic patients with mild traumatic brain injury, supporting the notion of ionic dysfunction in the pathogenesis of mild traumatic brain injury. The total sodium concentration appears to be altered beyond the immediate postinjury phase, and further work is needed to understand the relationship to persistent symptoms and outcome.

Improved Detection of Small Pulmonary Nodules Through Simultaneous MR/PET Imaging

Magnetic resonance (MR)/PET scanners provide an imaging platform that enables simultaneous acquisition of MR and PET data in perfect spatial and temporal registration. This feature allows improving image quality for the MR and PET images obtained during the course of an examination. In this work the authors demonstrate the use of prospective MR-based motion tracking information for removing motion blur in MR/PET images of small pulmonary nodules. The theoretical basis for the algorithms is presented alongside clinical examples of its use.

Low Rank plus Sparse decomposition of ODFs for improved detection of group-level differences and variable correlations in white matter

A novel approach is presented for group statistical analysis of diffusion weighted MRI datasets through voxelwise Orientation Distribution Functions (ODF). Recent advances in MRI acquisition make it possible to use high quality diffusion weighted protocols (multi-shell, large number of gradient directions) for routine in vivo study of white matter architecture. The dimensionality of these data sets is however often reduced to simplify statistical analysis. While these approaches may detect large group differences, they do not fully capitalize on all acquired image volumes. Incorporation of all available diffusion information in the analysis however risks biasing the outcome by outliers. Here we propose a statistical analysis method operating on the ODF, either the diffusion ODF or fiber ODF. To avoid outlier bias and reliably detect voxelwise group differences and correlations with demographic or behavioral variables, we apply the Low-Rank plus Sparse (L+S) matrix decomposition on the voxelwise ODFs which separates the sparse individual variability in the sparse matrix S whilst recovering the essential ODF features in the low-rank matrix L. We demonstrate the performance of this ODF L+S approach by replicating the established negative association between global white matter integrity and physical obesity in the Human Connectome dataset. The volume of positive findings p<0.01,227cm³, agrees with and expands on the volume found by TBSS (17 cm³), Connectivity based fixel enhancement (15 cm³) and Connectometry (212 cm³). In the same dataset we further localize the correlations of brain structure with neurocognitive measures such as fluid intelligence and episodic memory. The presented ODF L+S approach will aid in the full utilization of all acquired diffusion weightings leading to the detection of smaller group differences in clinically relevant settings as well as in neuroscience applications.

Accelerated radial diffusion spectrum imaging using a multi-echo stimulated echo diffusion sequence

PURPOSE

Diffusion spectrum imaging (DSI) provides us non-invasively and robustly with anatomical details of brain microstructure. To achieve sufficient angular resolution, DSI requires a large number of q-space samples, leading to long acquisition times. This need is mitigated here by combining the beneficial properties of Radial q-space sampling for DSI with a Multi-Echo Stimulated Echo Sequence (MESTIM).

METHODS

Full 2D k-spaces for each of several q-space samples, along the same radially outward line in q-space, are acquired in one readout train with one spin and three stimulated echoes. RF flip angles are carefully chosen to distribute spin magnetization over the echoes and the DSI reconstruction is adapted to account for differences in diffusion time among echoes.

RESULTS

Individual datasets and bootstrapped reproducibility analysis demonstrate image quality and SNR of the more-than-twofold-accelerated RDSI MESTIM sequence. Orientation distribution functions (ODF) and tractography results benefit from the longer diffusion times of the latter echoes in the echo train.

CONCLUSION

A MESTIM sequence can be used to shorten RDSI acquisition times significantly without loss of image or ODF quality. Further acceleration is possible by combination with simultaneous multi-slice techniques. Magn Reson Med 79:306-316, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Improved Detection of Small Pulmonary Nodules Through Simultaneous MR/PET Imaging

Magnetic resonance (MR)/PET scanners provide an imaging platform that enables simultaneous acquisition of MR and PET data in perfect spatial and temporal registration. This feature allows improving image quality for the MR and PET images obtained during the course of an examination. In this work the authors demonstrate the use of prospective MR-based motion tracking information for removing motion blur in MR/PET images of small pulmonary nodules. The theoretical basis for the algorithms is presented alongside clinical examples of its use.

Voxelwise analysis of simultaneously acquired and spatially correlated 18 F-fluorodeoxyglucose (FDG)-PET and intravoxel incoherent motion metrics in breast cancer

PURPOSE

Diffusion-weighted imaging (DWI) and ¹⁸ F-fluorodeoxyglucose-positron emission tomography (¹⁸ F-FDG-PET) independently correlate with malignancy in breast cancer, but the relationship between their structural and metabolic metrics is not completely understood. This study spatially correlates diffusion, perfusion, and glucose avidity in breast cancer with simultaneous PET/MR imaging and compares correlations with clinical prognostics.

METHODS

In this Health Insurance Portability and Accountability Act-compliant prospective study, with written informed consent and approval of the institutional review board and using simultaneously acquired FDG-PET and DWI, tissue diffusion (D_t ), and perfusion fraction (f_p ) from intravoxel incoherent motion (IVIM) analysis were registered to FDG-PET within 14 locally advanced breast cancers. Lesions were analyzed using 2D histograms and correlation coefficients between D_t , f_p , and standardized uptake value (SUV). Correlations were compared with prognostics from biopsy, metastatic burden from whole-body PET, and treatment history.

RESULTS

SUV||D_t correlation coefficient significantly distinguished treated (0.11 ± 0.24) from nontreated (-0.33 ± 0.26) patients (P = 0.005). SUV||f_p correlations were on average negative for the whole cohort (-0.17 ± 0.13).

CONCLUSION

Simultaneously acquired and registered FDG-PET/DWI allowed quantifiable descriptions of breast cancer microenvironments that may provide a framework for monitoring and predicting response to treatment. Magn Reson Med 78:1147-1156, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Dixon Sequence with Superimposed Model-Based Bone Compartment Provides Highly Accurate PET/MR Attenuation Correction of the Brain

Simultaneous PET/MR of the brain is a promising technology for characterizing patients with suspected cognitive impairment or epilepsy. Unlike CT, however, MR signal intensities do not correlate directly with PET photon attenuation correction (AC), and inaccurate radiotracer SUV estimation can limit future PET/MR clinical applications. We tested a novel AC method that supplements standard Dixon-based tissue segmentation with a superimposed model-based bone compartment.

METHODS

We directly compared SUV estimation between MR-based AC and reference CT AC in 16 patients undergoing same-day PET/CT and PET/MR with a single (18)F-FDG dose for suspected neurodegeneration. Three Dixon-based MR AC methods were compared with CT: standard Dixon 4-compartment segmentation alone, Dixon with a superimposed model-based bone compartment, and Dixon with a superimposed bone compartment and linear AC optimized specifically for brain tissue. The brain was segmented using a 3-dimensional T1-weighted volumetric MR sequence, and SUV estimations were compared with CT AC for whole-image, whole-brain, and 91 FreeSurfer-based regions of interest.

RESULTS

Modifying the linear AC value specifically for brain and superimposing a model-based bone compartment reduced the whole-brain SUV estimation bias of Dixon-based PET/MR AC by 95% compared with reference CT AC (P < 0.05), resulting in a residual -0.3% whole-brain SUVmean bias. Further, brain regional analysis demonstrated only 3 frontal lobe regions with an SUV estimation bias of 5% or greater (P < 0.05). These biases appeared to correlate with high individual variability in frontal bone thickness and pneumatization.

CONCLUSION

Bone compartment and linear AC modifications result in a highly accurate MR AC method in subjects with suspected neurodegeneration. This prototype MR AC solution appears equivalent to other recently proposed solutions and does not require additional MR sequences and scanning time. These data also suggest that exclusively model-based MR AC approaches may be adversely affected by common individual variations in skull anatomy.

Radial q-space sampling for DSI

PURPOSE

Diffusion spectrum imaging (DSI) has been shown to be an effective tool for noninvasively depicting the anatomical details of brain microstructure. Existing implementations of DSI sample the diffusion encoding space using a rectangular grid. Here we present a different implementation of DSI whereby a radially symmetric q-space sampling scheme for DSI is used to improve the angular resolution and accuracy of the reconstructed orientation distribution functions.

METHODS

Q-space is sampled by acquiring several q-space samples along a number of radial lines. Each of these radial lines in q-space is analytically connected to a value of the orientation distribution functions at the same angular location by the Fourier slice theorem.

RESULTS

Computer simulations and in vivo brain results demonstrate that radial diffusion spectrum imaging correctly estimates the orientation distribution functions when moderately high b-values (4000 s/mm2) and number of q-space samples (236) are used.

CONCLUSION

The nominal angular resolution of radial diffusion spectrum imaging depends on the number of radial lines used in the sampling scheme, and only weakly on the maximum b-value. In addition, the radial analytical reconstruction reduces truncation artifacts which affect Cartesian reconstructions. Hence, a radial acquisition of q-space can be favorable for DSI. Magn Reson Med 76:769-780, 2016. © 2015 Wiley Periodicals, Inc.

In vivo brain rosette spectroscopic imaging (RSI) with LASER excitation, constant gradient strength readout, and automated LCModel quantification for all voxels

PURPOSE

To optimize the Rosette trajectories for high-sensitivity in vivo brain spectroscopic imaging and reduced gradient demands.

METHODS

Using LASER localization, a rosette based sampling scheme for in vivo brain spectroscopic imaging data on a 3 Tesla (T) system is described. The two-dimensional (2D) and 3D rosette spectroscopic imaging (RSI) data were acquired using 20 × 20 in-plane resolution (8 × 8 mm(2) ), and 1 (2D) -18 mm (1.1 cc) or 12 (3D) -8 mm partitions (0.5 cc voxels). The performance of the RSI acquisition was compared with a conventional spectroscopic imaging (SI) sequence using LASER localization and 2D or 3D elliptical phase encoding (ePE). Quantification of the entire RSI data set was performed using an LCModel based pipeline.

RESULTS

The RSI acquisitions took 32 s for the 2D scan, and as short as 5 min for the 3D 20 × 20 × 12 scan, using a maximum gradient strength Gmax=5.8 mT/m and slew-rate Smax=45 mT/m/ms. The Bland-Altman agreement between RSI and ePE CSI, characterized by the 95% confidence interval for their difference (RSI-ePE), is within 13% of the mean (RSI+ePE)/2. Compared with the 3D ePE at the same nominal resolution, the effective RSI voxel size was three times smaller while the measured signal-to-noise ratio sensitivity, after normalization for differences in effective size, was 43% greater.

CONCLUSION

3D LASER-RSI is a fast, high-sensitivity spectroscopic imaging sequence, which can acquire medium-to-high resolution SI data in clinically acceptable scan times (5-10 min), with reduced stress on the gradient system. Magn Reson Med 76:380-390, 2016. © 2015 Wiley Periodicals, Inc.

Short-T2 imaging for quantifying concentration of sodium (23 Na) of bi-exponential T2 relaxation

PURPOSE

This work intends to demonstrate a new method for quantifying concentration of sodium (²³ Na) of bi-exponential T₂ relaxation in patients on MRI scanners at 3.0 Tesla.

THEORY AND METHODS

Two single-quantum (SQ) sodium images acquired at very-short and short echo times (TE = 0.5 and 5.0 ms) are subtracted to produce an image of the short-T₂ component of the bi-exponential (or bound) sodium. An integrated calibration on the SQ and short-T₂ images quantifies both total and bound sodium concentrations. Numerical models were used to evaluate signal response of the proposed method to the short-T₂ components. MRI scans on agar phantoms and brain tumor patients were performed to assess accuracy and performance of the proposed method, in comparison with a conventional method of triple-quantum filtering.

RESULTS

A good linear relation (R²  = 0.98) was attained between the short-T₂ image intensity and concentration of bound sodium. A reduced total scan time of 22 min was achieved under the SAR restriction for human studies in quantifying both total and bound sodium concentrations.

CONCLUSION

The proposed method is feasible for quantifying bound sodium concentration in routine clinical settings at 3.0 Tesla. Magn Reson Med 74:162-174, 2015. © 2014 Wiley Periodicals, Inc.

Molecular and metabolic pattern classification for detection of brain glioma progression

OBJECTIVES

The ability to differentiate between brain tumor progression and radiation therapy induced necrosis is critical for appropriate patient management. In order to improve the differential diagnosis, we combined fluorine-18 2-fluoro-deoxyglucose positron emission tomography ((18)F-FDG PET), proton magnetic resonance spectroscopy ((1)H MRS) and histological data to develop a multi-parametric machine-learning model.

METHODS

We enrolled twelve post-therapy patients with grade 2 and 3 gliomas that were suspicious of tumor progression. All patients underwent (18)F-FDG PET and (1)H MRS. Maximal standardized uptake value (SUVmax) of the tumors and reference regions were obtained. Multiple 2D maps of choline (Cho), creatine (Cr), and N-acetylaspartate (NAA) of the tumors were generated. A support vector machine (SVM) learning model was established to take imaging biomarkers and histological data as input vectors. A combination of clinical follow-up and multiple sequential MRI studies served as the basis for assessing the clinical outcome. All vector combinations were evaluated for diagnostic accuracy and cross validation. The optimal cutoff value of individual parameters was calculated using Receiver operating characteristic (ROC) plots.

RESULTS

The SVM and ROC analyses both demonstrated that SUVmax of the lesion was the most significant single diagnostic parameter (75% accuracy) followed by Cho concentration (67% accuracy). SVM analysis of all paired parameters showed SUVmax and Cho concentration in combination could achieve 83% accuracy. SUVmax of the lesion paired with SUVmax of the white matter as well as the tumor Cho paired with the tumor Cr both showed 83% accuracy. These were the most significant paired diagnostic parameters of either modality. Combining all four parameters did not improve the results. However, addition of two more parameters, Cho and Cr of brain parenchyma contralateral to the tumor, increased the accuracy to 92%.

CONCLUSION

This study suggests that SVM models may improve detection of glioma progression more accurately than single parametric imaging methods.

RESEARCH SUPPORT

National Cancer Institute, Cancer Center Support Grant Supplement Award, Imaging Response Assessment Teams.

Multi-slice parallel transmission three-dimensional tailored RF (PTX 3DTRF) pulse design for signal recovery in ultra high field functional MRI

T(2)(∗) weighted fMRI at high and ultra high field (UHF) is often hampered by susceptibility-induced, through-plane, signal loss. Three-dimensional tailored RF (3DTRF) pulses have been shown to be an effective approach for mitigating through-plane signal loss at UHF. However, the required RF pulse lengths are too long for practical applications. Recently, parallel transmission (PTX) has emerged as a very effective means for shortening the RF pulse duration for 3DTRF without sacrificing the excitation performance. In this article, we demonstrate a RF pulse design strategy for 3DTRF based on the use of multi-slice PTX 3DTRF to simultaneously and precisely recover signal with whole-brain coverage. Phantom and human experiments are used to demonstrate the effectiveness and robustness of the proposed method on three subjects using an eight-channel whole body parallel transmission system.

Repeatability of ultrashort echo time-based two-component T2* measurements on cartilages in human knee at 3 T

Repeatability of in vivo measurement of multicomponent T2* relaxation in articular cartialges in human knee is important to clinical use. This study evaluated the repeatability of two-component T2* relaxation on seven healthy human subjects. The left knee was scanned once a day in three consecutive days, on a clinical 3T MRI scanner with eight-channel knee coil and ultrashort echo time pulse sequence at 11 echo times=0.6-40 ms. The intrasubject and intersubject repeatability was evaluated via coefficient of variation (CV=standard deviation/mean) in four typical cartilage regions: patellar, anterior articular, femoral, and tibial regions. It was found that the intrasubject repeatability was good, with CV<10% for the short- and long-T2* relaxation time in the layered regions in the four cartilages (with one exception) and CV<13% for the component intensity fraction (with two exceptions). The intersubject repeatability was also good, with CV∼8% (range 1-15%) for the short- and long-T2* relaxation time and CV∼10% (range 2-20%) for the component intensity fraction. The long-T2* component showed significantly better repeatability (CV∼8%) than the short-T2* component (CV∼12%) (P<0.005). These CV values suggest that in vivo measurement of two-component T2* relaxation in the knee cartilages is repeatable on clinical scanner at 3 T, with a signal-to-noise ratio of 90.

Combined imaging biomarkers for therapy evaluation in glioblastoma multiforme: correlating sodium MRI and F-18 FLT PET on a voxel-wise basis

We evaluate novel magnetic resonance imaging (MRI) and positron emission tomography (PET) quantitative imaging biomarkers and associated multimodality, serial-time-point analysis methodologies, with the ultimate aim of providing clinically feasible, predictive measures for early assessment of response to cancer therapy. A focus of this work is method development and an investigation of the relationship between the information content of the two modalities. Imaging studies were conducted on subjects who were enrolled in glioblastoma multiforme (GBM) therapeutic clinical trials. Data were acquired, analyzed and displayed using methods that could be adapted for clinical use. Subjects underwent dynamic [(18)F]fluorothymidine (F-18 FLT) PET, sodium ((23)Na) MRI and 3-T structural MRI scans at baseline (before initiation of therapy), at an early time point after beginning therapy and at a late follow-up time point after therapy. Sodium MRI and F-18 FLT PET images were registered to the structural MRI. F-18 FLT PET tracer distribution volumes and sodium MRI concentrations were calculated on a voxel-wise basis to address the heterogeneity of tumor physiology. Changes in, and differences between, these quantities as a function of scan timing were tracked. While both modalities independently show a change in tissue status as a function of scan time point, results illustrate that the two modalities may provide complementary information regarding tumor progression and response. Additionally, tumor status changes were found to vary in different regions of tumor. The degree to which these methods are useful for GBM therapy response assessment and particularly for differentiating true progression from pseudoprogression requires additional patient data and correlation of these imaging biomarker changes with clinical outcome.

Sodium MRI and the assessment of irreversible tissue damage during hyper-acute stroke

Sodium MRI (sMRI) has undergone a tremendous amount of technical development during the last two decades that makes it a suitable tool for the study of human pathology in the acute setting within the constraints of a clinical environment. The salient role of the sodium ion during impaired ATP production during the course of brain ischemia makes sMRI an ideal tool for the study of ischemic tissue viability during stroke. In this paper, the current limitations of conventional MRI for the determination of tissue viability during evolving brain ischemia are discussed. This discussion is followed by a summary of the known findings about the dynamics of tissue sodium changes during brain ischemia. A mechanistic model for the explanation of these findings is presented together with the technical requirements for its investigation using clinical MRI scanners. An illustration of the salient features of the technique is also presented using a nonhuman primate model of reversible middle cerebral artery occlusion.

Comparison of proton magnetic resonance spectroscopy with fluorine-18 2-fluoro-deoxyglucose positron emission tomography for assessment of brain tumor progression

OBJECTIVES

We investigated the accuracy of high-field proton magnetic resonance spectroscopy ((1) H MRS) and fluorine-18 2-fluoro-deoxyglucose positron emission tomography ((18) F-FDG-PET) for diagnosis of glioma progression following tumor resection, stereotactic radiation, and chemotherapy.

METHODS

Twelve post-therapy patients with histology proven gliomas (six grade II and six grade III) presented with magnetic resonance imaging (MRI) and clinical symptoms suggestive but not conclusive of progression were entered into the study. (1) H MRS data were acquired and 3-dimensional volumetric maps of choline (Cho) over creatine (Cr) were generated. Intensity of (18) F-FDG uptake was evaluated on a semiquantitative scale.

RESULTS

The accuracy of (1) H MRS and (18) F-FDG-PET imaging for diagnosis of glioma progression was 75% and 83%, respectively. Classifying the tumors by grade improved accuracy of (18) F-FDG-PET to 100% in high-grade gliomas and accuracy of (1) H MRS to 80% in low-grade tumors. Spearman's analysis demonstrated a trend between (18) F-FDG uptake and tumor grading (ρ= .612, P-value = .272). The results of (18) F-FDG-PET and (1) H MRS were concordant in 75% (9/12) of cases.

CONCLUSION

The combination of (1) H MRS data and (18) F-FDG-PET imaging can enhance detection of glioma progression. (1) H MRS imaging was more accurate in low-grade gliomas and (18) F-FDG-PET provided better accuracy in high-grade gliomas.

High-definition fiber tracking for assessment of neurological deficit in a case of traumatic brain injury: finding, visualizing, and interpreting small sites of damage

For patients with traumatic brain injury (TBI), current clinical imaging methods generally do not provide highly detailed information about the location of axonal injury, severity of injury, or expected recovery. In a case of severe TBI, the authors applied a novel high-definition fiber tracking (HDFT) to directly visualize and quantify the degree of axonal fiber damage and predict functional deficits due to traumatic axonal injury and loss of cortical projections. This 32-year-old man sustained a severe TBI. Computed tomography and MRI revealed an area of hemorrhage in the basal ganglia with mass effect, but no specific information on the location of axonal injury could be obtained from these studies. Examinations of the patient at Week 3 and Week 8 after TBI revealed motor weaknesses of the left extremities. Four months postinjury, 257-direction diffusion spectrum imaging and HDFT analysis was performed to evaluate the degree of axonal damage in the motor pathway and quantify asymmetries in the left and right axonal pathways. High-definition fiber tracking was used to follow corticospinal and corona radiata pathways from the cortical surface to the midbrain and quantify projections from motor areas. Axonal damage was then localized by assessing the number of descending fibers at the level of the cortex, internal capsule, and midbrain. The motor deficit apparent in the clinical examinations correlated with the axonal losses visualized using HDFT. Fiber loss estimates at 4 months postinjury accurately predicted the nature of the motor deficits (severe, focal left-hand weakness) when other standard clinical imaging modalities did not. A repeat scan at 10 months postinjury, when edema and hemorrhage had receded, replicated the fiber loss. Using HDFT, the authors accurately identified the presence and location of damage to the underlying white matter in this patient with TBI. Detailed information of injury provided by this novel technique holds future potential for precise neuroimaging assessment of TBI.

High-resolution sodium imaging of human brain at 7 T

The feasibility of high-resolution sodium magnetic resonance imaging on human brain at 7 T was demonstrated in this study. A three-dimensional anisotropic resolution data acquisition was used to address the challenge of low signal-to-noise ratio associated with high resolution. Ultrashort echo-time sequence was used for the anisotropic data acquisition. Phantoms and healthy human brains were studied on a whole-body 7-T magnetic resonance imaging scanner. Sodium images were obtained at two high nominal in-plane resolutions (1.72 and 0.86 mm) at a slice thickness of 4 mm. Signal-to-noise ratio in the brain image (cerebrospinal fluid) was measured as 14.4 and 6.8 at the two high resolutions, respectively. The actual in-plane resolution was measured as 2.9 and 1.6 mm, 69-86% larger than their nominal values. The quantification of sodium concentration on the phantom and brain images enabled better accuracy at the high nominal resolutions than at the low nominal resolution of 3.44 mm (measured resolution 5.5 mm) due to the improvement of in-plane resolution.

High-resolution ultrashort echo time (UTE) imaging on human knee with AWSOS sequence at 3.0 T

PURPOSE

To demonstrate the technical feasibility of high-resolution (0.28-0.14 mm) ultrashort echo time (UTE) imaging on human knee at 3T with the acquisition-weighted stack of spirals (AWSOS) sequence.

MATERIALS AND METHODS

Nine human subjects were scanned on a 3T MRI scanner with an 8-channel knee coil using the AWSOS sequence and isocenter positioning plus manual shimming.

RESULTS

High-resolution UTE images were obtained on the subject knees at TE = 0.6 msec with total acquisition time of 5.12 minutes for 60 slices at an in-plane resolution of 0.28 mm and 10.24 minutes for 40 slices at an in-plane resolution of 0.14 mm. Isocenter positioning, manual shimming, and the 8-channel array coil helped minimize image distortion and achieve high signal-to-noise ratio (SNR).

CONCLUSION

It is technically feasible on a clinical 3T MRI scanner to perform UTE imaging on human knee at very high spatial resolutions (0.28-0.14 mm) within reasonable scan time (5-10 min) using the AWSOS sequence.

Improved large tip angle parallel transmission pulse design through a perturbation analysis of the Bloch equation

Parallel transmission has emerged as an efficient means for implementing multidimensional spatially selective radiofrequency excitation pulses. To date, most theoretical and experimental work on parallel transmission radiofrequency (RF) pulse design is based on the small-tip-angle approximation to the Bloch equation. The small-tip-angle, while mathematically compact, is not an exact solution and leads to significant errors when large-tip-angle pulses are designed. Methods have been proposed to overcome the limitations of the small-tip-angle using regularized least-square optimization or optimal control algorithms. These methods, however, are based on further approximations to the Bloch equation or require the use of general purpose algorithms that do not capitalize fully on the dynamics of the physical model at hand. In this article, a novel algorithm for large-tip-angle parallel transmission pulse design is proposed. The algorithm relies on a perturbation analysis of the Bloch equation and it depicts the relationship between the excited magnetization, its deviation from the target pattern and the desired pulses. Simulations and experiments are used to validate the proposed method on a 7 T 8-channel transmit array. The results demonstrate that the perturbation analysis algorithm provides a fast and accurate approach for multidimensional large-tip-angle pulse design, especially when large acceleration factors and/or echo-planar trajectories are used.

K⁺ dynamics in ischemic rat brain in vivo by ⁸⁷Rb MRI at 7 T

The aims of the present study were as follows: (i) to perform the first (87)Rb MRI in live rats with focal ischemic stroke; and (ii) to test the hypothesis that K(+) egress from the brain in this model is quantifiable in individual animals by high-field (7-T) K/Rb substitution MRI. Rats preloaded with dietary Rb(+) (resulting in Rb/(K + Rb) replacement ratios of 0.1-0.2 in the brain) were subjected to permanent occlusion of the middle cerebral artery, and (87)Rb MRI was implemented with 13-min temporal resolution using a dedicated RF coil and a spiral ultrashort-TE sequence (TR/TE = 3/0.07 ms). The ischemic core was localized by apparent diffusion coefficient mapping, by microtubule-associated protein-2 immunohistochemistry and by changes in surface reflectivity. [K], [Na] and [Rb] were determined independently in the micropunched samples by post-mortem flame photometry. Both techniques were generally in agreement in the nonischemic cortex; however, the MRI-assessed [K(+) + Rb(+)] drop in ischemic brain was less pronounced (average efflux rate of 4.8 ± 0.2 nEq/mm(3) /h versus 10 ± 1 nEq/mm(3)/h by flame photometry; p < 0.0001). The use of higher field gradients for better spatial resolution, and hence more accurate quantification, is suggested.

Multicomponent T2* mapping of knee cartilage: technical feasibility ex vivo

Disorganization of collagen fibers is a sign of early-stage cartilage degeneration in osteoarthritic knees. Water molecules trapped within well-organized collagen fibrils would be sensitive to collagen alterations. Multicomponent effective transverse relaxation (T2*) mapping with ultrashort echo time acquisitions is here proposed to probe short T(2) relaxations in those trapped water molecules. Six human tibial plateau explants were scanned on a 3T MRI scanner using a home-developed ultrashort echo time sequence with echo times optimized via Monte Carlo simulations. Time constants and component intensities of T2* decays were calculated at individual pixels, using the nonnegative least squares algorithm. Four T2*-decay types were found: 99% of cartilage pixels having mono-, bi-, or nonexponential decay, and 1% showing triexponential decay. Short T2* was mainly in 1-6 ms, while long T2* was ∼ 22 ms. A map of decay types presented spatial distribution of these T2* decays. These results showed the technical feasibility of multicomponent T2* mapping on human knee cartilage explants.

Whole-body MRI and ethnic differences in adipose tissue and skeletal muscle distribution in overweight black and white adolescent boys

It is unclear whether ethnic differences exist in adipose tissue (AT) and skeletal muscle (SM) distribution in black and white youth. Investigation into the pattern of AT and SM distribution in black versus white youth may provide insight into the previously reported health disparities between these ethnicities. Therefore, we examined total and regional AT and SM in overweight black and white boys. The study sample included overweight black (n = 19) and white (n = 21) boys (11-18 yr, BMI ≥ 85th) whose body composition was evaluated using whole-body MRI. Despite similar age, Tanner stage, and BMI, black boys had significantly (P < .05) less visceral AT than white boys and more (P < .05) total and lower-body subcutaneous AT (SAT) in both absolute (kg) and relative (%) terms. There was a main effect (P < .05) of ethnicity on the relationship between total and regional AT, such that for a given amount of total body AT (kg), black boys had a greater (P < .05) lower-body SAT and less visceral AT than their white peers. For a given amount of total SM, black boys had more (P < .05) SM in the thigh. Compared with overweight white boys, overweight black boys have less visceral fat, more subcutaneous fat, and more thigh skeletal muscle.

High-definition fiber tracking guidance for intraparenchymal endoscopic port surgery

The authors have applied high-definition fiber tracking (HDFT) to the resection of an intraparenchymal dermoid cyst by using a minimally invasive endoscopic port. The lesion was located within the mesial frontal lobe, septal area, hypothalamus, and suprasellar recess. Using high-dimensional (256 directions) diffusion imaging, more than 250,000 fiber tracts were imaged before and after surgery. Trajectory planning using HDFT in a computer model was used to facilitate cannulation of the cyst with the endoscopic port. Analysis of the proposed initial surgical route was overlaid onto the fiber tracts and was predicted to produce substantial disruption to prefrontal projection fibers (anterior limb of the internal capsule) and the cingulum. Adjustment of the cannulation entry point 1 cm medially was predicted to cross the corpus callosum instead of the anterior limb of the internal capsule or the cingulum. Following cyst resection performed using endoscopic port surgery, postoperative imaging demonstrated accurate cannulation of the lesion, with improved quantitative signal from both the anterior limb of the internal capsule and the cingulum. The observed fiber preservation from the cingulum and the anterior limb of the internal capsule, with minor injury to the corpus callosum, was in close agreement with preoperative trajectory modeling. Comparison of pre- and postoperative HDFT data facilitated quantification of the benefits and costs of the surgical trajectory. Future studies will help to determine whether HDFT combined with endoscopic port surgery facilitates anatomical and functional preservation in such challenging cases.

Rosette spectroscopic imaging: optimal parameters for alias-free, high sensitivity spectroscopic imaging

PURPOSE

To optimize the Rosette trajectories for fast, high sensitivity spectroscopic imaging experiments and to compare this acquisition technique with other chemical shift imaging (CSI) methods.

MATERIALS AND METHODS

A framework for comparing the sensitivity of the Rosette Spectroscopic Imaging (RSI) acquisition to other spectroscopic imaging experiments is outlined. Accounting for hardware constraints, trajectory parameters that provide for optimal sampling and minimal artifact production are found. Along with an analytical expression for the number of excitations to be used in an RSI experiment that is provided, the theoretical precompensation weights used for optimal image reconstruction are derived.

RESULTS

The spectral response function for RSI is shown to be approximately the same as the point spread function of standard Fourier reconstructions. While the signal-to-noise ratio (SNR) for an RSI experiment is reduced by the inherent nonuniform sampling of these trajectories, their circular k-space support and speed of spatial encoding leads to greater SNR efficiency and improvements in the total data acquisition time relative to the gold standard CSI approach with square k-space support and to similar efficiency to spiral CSI acquisitions. Numerical simulations and in vivo experimental data are presented to demonstrate the properties of this data acquisition technique.

CONCLUSION

This work demonstrates the use of Rosette trajectories and how to achieve improved efficiency for these trajectories in a two-dimensional spectroscopic imaging experiment.

Inhomogeneous sodium accumulation in the ischemic core in rat focal cerebral ischemia by 23Na MRI

PURPOSE

To test the hypotheses that (i) the regional heterogeneity of brain sodium concentration ([Na(+)](br)) provides a parameter for ischemic progression not available from apparent diffusion coefficient (ADC) data, and (ii) [Na(+)](br) increases more in ischemic cortex than in the caudate putamen (CP) with its lesser collateral circulation after middle cerebral artery occlusion in the rat.

MATERIALS AND METHODS

(23)Na twisted projection MRI was performed at 3 Tesla. [Na(+)](br) was independently determined by flame photometry. The ischemic core was localized by ADC, by microtubule-associated protein-2 immunohistochemistry, and by changes in surface reflectivity.

RESULTS

Within the ischemic core, the ADC ratio relative to the contralateral tissue was homogeneous (0.63 +/- 0.07), whereas the rate of [Na(+)](br) increase (slope) was heterogeneous (P < 0.005): 22 +/- 4%/h in the sites of maximum slope versus 14 +/- 1%/h elsewhere (here 100% is [Na(+)](br) in the contralateral brain). Maximum slopes in the cortex were higher than in CP (P < 0.05). In the ischemic regions, there was no slope/ADC correlation between animals and within the same brain (P > 0.1). Maximum slope was located at the periphery of ischemic core in 8/10 animals.

CONCLUSION

Unlike ADC, (23)Na MRI detected within-core ischemic lesion heterogeneity.

Sodium mapping in focal cerebral ischemia in the rat by quantitative (23)Na MRI

PURPOSE

To validate (23)Na twisted projection magnetic resonance imaging (MRI) as a quantitative technique to assess local brain sodium concentration ([Na(+)](br)) during rat focal ischemia every 5.3 minutes.

MATERIALS AND METHODS

The MRI protocol included an ultrashort echo-time (0.4 msec), a correction of radiofrequency (RF) inhomogeneities by B(1) mapping, and the use of 0-154 mM NaCl calibration standards. To compare MRI [Na(+)](br) values with those obtained by emission flame photometry in precision-punched brain samples of about 0.5 mm(3) size, MR images were aligned with a histological three-dimensional reconstruction of the punched brain and regions of interest (ROIs) were placed precisely over the punch voids.

RESULTS

The Bland-Altman analysis of [Na(+)](br) in normal and ischemic cortex and caudate putamen of seven rats quantitated by (23)Na MRI and flame photometry yielded a mean bias and limits of agreement (at +/-1.96 SD) of 2% and 43% of average, respectively. A linear increase in [Na(+)](br) was observed between 1 and 6 hours after middle cerebral artery occlusion.

CONCLUSION

(23)Na MRI provides accurate and reliable results within the whole range of [Na(+)](br) in ischemia with a temporal resolution of 5.3 minutes and precisely targeted submicroliter ROIs in selected brain structures.

Simultaneous z-shim method for reducing susceptibility artifacts with multiple transmitters

The signal loss susceptibility artifact is a major limitation in gradient-echo MRI applications. Various methods, including z-shim techniques and multidimensional tailored radio frequency (RF) pulses, have been proposed to mitigate the through-plane signal loss artifact, which is dominant in axial slices above the sinus region. Unfortunately, z-shim techniques require multiple steps and multidimensional RF methods are complex, with long pulse lengths. Parallel transmission methods were recently shown to be promising for improving B1 inhomogeneity and reducing the specific absorption rate. In this work, a novel method using time-shifted slice-select RF pulses is presented for reducing the through-plane signal loss artifact in parallel transmission applications. A simultaneous z-shim is obtained by concurrently applying unique time-shifted pulses on each transmitter. The method is shown to reduce the signal loss susceptibility artifact in gradient-echo images using a four-channel parallel transmission system at 3T.

Acquisition-weighted stack of spirals for fast high-resolution three-dimensional ultra-short echo time MR imaging

Ultra-short echo time (UTE) MRI requires both short excitation ( approximately 0.5 ms) and short acquisition delay (<0.2 ms) to minimize T(2)-induced signal decay. These requirements currently lead to low acquisition efficiency when high resolution (<1 mm) is pursued. A novel pulse sequence, acquisition-weighted stack of spirals (AWSOS), is proposed here to acquire high-resolution three-dimensional (3D) UTE images with short scan time ( approximately 72 s). The AWSOS sequence uses variable-duration slice encoding to minimize T(2) decay, separates slice thickness from in-plane resolution to reduce the number of slice encodings, and uses spiral trajectories to accelerate in-plane data collections. T(2)- and off-resonance induced slice widening and image blurring were calculated from 1.5 to 7 Tesla (T) through point spread function. Computer simulations were performed to optimize spiral interleaves and readout times. Phantom scans and in vivo experiments on human heads were implemented on a clinical 1.5T scanner (G(max) = 40 mT/m, S(max) = 150 T/m/s). Accounting for the limits on B(1) maximum, specific absorption rate (SAR), and the lowered amplitude of slab-select gradient, a sinc radiofrequency (RF) pulse of 0.8ms duration and 1.5 cycles was found to produce a flat slab profile. High in-plane resolution (0.86 mm) images were obtained for the human head using echo time (TE) = 0.608 ms and total shots = 720 (30 slice-encodings x 24 spirals). Compared with long-TE (10 ms) images, the ultrashort-TE AWSOS images provided clear visualization of short-T(2) tissues such as the nose cartilage, the eye optic nerve, and the brain meninges and parenchyma.

Reduction of transmitter B1 inhomogeneity with transmit SENSE slice-select pulses

Parallel transmitter techniques are a promising approach for reducing transmitter B1 inhomogeneity due to the potential for adjusting the spatial excitation profile with independent RF pulses. These techniques may be further improved with transmit sensitivity encoding (SENSE) methods because the sensitivity information in pulse design provides an excitation that is inherently compensated for transmitter B1 inhomogeneity. This paper presents a proof of this concept using transmit SENSE 3D tailored RF pulses designed for small flip angles. An eight-channel receiver coil was used to mimic parallel transmission for brain imaging at 3T. The transmit SENSE pulses were based on the fast-k(z) design and produced 5-mm-thick slices at a flip angle of 30 degrees with only a 4.3-ms pulse length. It was found that the transmit SENSE pulses produced more homogeneous images than those obtained from the complex sum of images from all receivers excited with a standard RF pulse.

Non-invasive assessment of tumor proliferation using triple quantum filtered 23/Na MRI: technical challenges and solutions

We address the development of triple-quantum-filtered sodium MRI as a non-invasive surrogate measure for cell proliferation in brain tumors. We demonstrate that through careful consideration of the theoretical description of the signal, triple-quantum-filtered sodium images of adequate signal-to-noise ratio (SNR) can be acquired in clinically acceptable imaging times.

Monitoring of brain potassium with rubidium flame photometry and MRI

An animal model was developed to monitor [K(+)] in the brain using partial K(+) replacement with Rb(+) and (87)Rb MRI. Fifty-one rats were given 0-80 mM of RbCl in the drinking water for up to 90 days. Focal cerebral ischemia was produced in 15 of the animals. Na, K, and Rb content in precision-guided submilligram samples of cortical brain were determined by emission flame photometry. Multinuclear (87)Rb/(23)Na/(1)H MRI was performed on phantoms and rats at 3T using a twisted projection imaging (TPI) scheme for (87)Rb/(23)Na, and custom-built surface or parallel cosine transmit/receive coils. Brain [Rb(+)] was safely brought up to 17-25 mEq/kg within 2-3 weeks of feeding. The characteristic patterns of [K(+)] decrease (with a sharp drop at 3-4 hr of ischemia) and [Na(+)] increase (at a rate of 31%/hr) observed previously in animals without Rb/K substitution were reproduced in ischemic cortex. The Rb/(Rb+K) ratio increased over time in ischemic areas (R = 0.91, P < 0.001), suggesting an additional index of ischemia progression. Preliminary (87)Rb MRI gave an estimate of 20-25 mEq Rb/kg brain weight (N = 2). In conclusion, brain Rb(+) is detectable by (87)Rb MRI and does not significantly interfere with ion dynamics in ischemic brain, which enables (87)Rb MRI studies of K(+) in ischemia.

The relation of white matter hyperintensities to cognitive performance in the normal old: education matters

This study examined whether the severity of cerebral white matter abnormalities (evident on MR images as white matter hyperintensities (WMH)) was related to the cognitive performance of 141 high-functioning older adults. The elderly showed the typical age decrement on measures of processing speed, working memory, and inhibition; however WMH severity was significantly related only to processing speed. The strength of this relationship was, however, influenced by the educational level of the participants, such that processing speed was more associated with WMH severity in less-educated than in well-educated participants. This is consistent with recent concepts of cognitive reserve, but does raise a question as to the underlying source of the cognitive decrement found in the sort of well-educated elders typically used in cognitive-aging studies.

Excitation UNFOLD (XUNFOLD) to improve the temporal resolution of multishot tailored RF pulses

An extension of the "UNaliasing by Fourier encoding the Overlaps using the temporaL Dimension" (UNFOLD) method to the excitation domain (XUNFOLD) is presented to improve the temporal resolution of multishot tailored RF (TRF) pulses. Multishot three-dimensional TRF pulses were designed to produce a time series of images with periodically aliased excitation profiles. The XUNFOLD method is shown to remove the excitation profile aliasing from the dynamic imaging data by filtering in the temporal frequency dimension. The technique is demonstrated to improve the temporal resolution of simulated functional MRI (fMRI) activation in a time series of brain images.

Fast-kz three-dimensional tailored radiofrequency pulse for reduced B1 inhomogeneity

This article presents a small-flip-angle, three-dimensional tailored RF pulse that excites thin slices with an adjustable quadratic in-plane spatial variation. The quadratic spatial variation helps to compensate for the loss in image uniformity using a volume coil at 3 T due to the wavelike properties of the RF field. The pulse is based on a novel "fast-kz" design that uses a series of slice-select subpulses along kz and phase encoding "blips" along kx-ky. The method is demonstrated by acquiring a series of 5-mm-thick T2-weighted images of the human brain at 3 T using pulses 4.8 ms in length with a 45 degrees flip angle.

Stroke onset time using sodium MRI in rat focal cerebral ischemia

BACKGROUND AND PURPOSE

Thrombolytic therapy with intravenous tPA must be administered within 3 hours after stroke onset. However, stroke onset time cannot be established in 20% to 45% of potential patients. We propose that the rate of increase of the brain concentration of sodium ([Na+]br) after stroke, monitored using sodium MRI in a rat model of cortical ischemia, is linear in each individual animal, can locate the ischemic region, and can be used to estimate onset time.

METHODS

After induction of focal cortical ischemia in rats under isoflurane anesthesia, [Na+]br time course maps were acquired continuously on a 3 T whole body scanner from 2 to 7 hours after occlusion followed by T2-weighted proton images. Microtubule-associated protein-2 immunostained brain sections were used to verify the location of the infarct.

RESULTS

The ischemic region identified with microtubule-associated protein-2 corresponded to the region of maximum [Na+]br increase (P<0.001; n=5), and all of the animals demonstrated high linearity. [Na+]br increased at a mean rate of 25+/-4.7%/h in ischemic tissue (P=0.013) but not in normal cortex (1.0+/-1.1%/h; P=0.42). The mean onset time error was 1+/-4 minutes (n=4).

CONCLUSIONS

These results of sodium MRI show that the region of maximum [Na+]br increase corresponds to the ischemic region. Although [Na+]br increases at a different rate in each animal, the increase is linear, and, therefore, onset time can be estimated. These findings suggest that this method can be used as a ticking clock to estimate time elapsed after vascular occlusion.