Fast image reconstruction with L2-regularization

Noise-reduction techniques for 1H-FID-MRSI at 14.1 T: Monte Carlo validation and in vivo application

Proton magnetic resonance spectroscopic imaging (1H-MRSI) is a powerful tool that enables the multidimensional non-invasive mapping of the neurochemical profile at high resolution over the entire brain. The constant demand for higher spatial resolution in 1H-MRSI has led to increased interest in post-processing-based denoising methods aimed at reducing noise variance. The aim of the present study was to implement two noise-reduction techniques, Marchenko-Pastur principal component analysis (MP-PCA) based denoising and low-rank total generalized variation (LR-TGV) reconstruction, and to test their potential with and impact on preclinical 14.1 T fast in vivo 1H-FID-MRSI datasets. Since there is no known ground truth for in vivo metabolite maps, additional evaluations of the performance of both noise-reduction strategies were conducted using Monte Carlo simulations. Results showed that both denoising techniques increased the apparent signal-to-noise ratio (SNR) while preserving noise properties in each spectrum for both in vivo and Monte Carlo datasets. Relative metabolite concentrations were not significantly altered by either method and brain regional differences were preserved in both synthetic and in vivo datasets. Increased precision of metabolite estimates was observed for the two methods, with inconsistencies noted for lower-concentration metabolites. Our study provided a framework for how to evaluate the performance of MP-PCA and LR-TGV methods for preclinical 1H-FID MRSI data at 14.1 T. While gains in apparent SNR and precision were observed, concentration estimations ought to be treated with care, especially for low-concentration metabolites.

Deep learning pipeline for quality filtering of MRSI spectra

With the rise of novel 3D magnetic resonance spectroscopy imaging (MRSI) acquisition protocols in clinical practice, which are capable of capturing a large number of spectra from a subject's brain, there is a need for an automated preprocessing pipeline that filters out bad-quality spectra and identifies contaminated but salvageable spectra prior to the metabolite quantification step. This work introduces such a pipeline based on an ensemble of deep-learning classifiers. The dataset consists of 36,338 spectra from one healthy subject and five brain tumor patients, acquired with an EPSI variant, which implemented a novel type of spectral editing named SLOtboom-Weng (SLOW) editing on a 7T MR scanner. The spectra were labeled manually by an expert into four classes of spectral quality as follows: (i) noise, (ii) spectra greatly influenced by lipid-related artifacts (deemed not to contain clinical information), (iii) spectra containing metabolic information slightly contaminated by lipid signals, and (iv) good-quality spectra. The AI model consists of three pairs of networks, each comprising a convolutional autoencoder and a multilayer perceptron network. In the classification step, the encoding half of the autoencoder is kept as a dimensionality reduction tool, while the fully connected layers are added to its output. Each of the three pairs of networks is trained on different representations of spectra (real, imaginary, or both), aiming at robust decision-making. The final class is assigned via a majority voting scheme. The F1 scores obtained on the test dataset for the four previously defined classes are 0.96, 0.93, 0.82, and 0.90, respectively. The arguably lower value of 0.82 was reached for the least represented class of spectra mildly influenced by lipids. Not only does the proposed model minimise the required user interaction, but it also greatly reduces the computation time at the metabolite quantification step (by selecting a subset of spectra worth quantifying) and enforces the display of only clinically relevant information.

Deep learning unlocks label-free viability assessment of cancer spheroids in microfluidics

Despite recent advances in cancer treatment, refining therapeutic agents remains a critical task for oncologists. Precise evaluation of drug effectiveness necessitates the use of 3D cell culture instead of traditional 2D monolayers. Microfluidic platforms have enabled high-throughput drug screening with 3D models, but current viability assays for 3D cancer spheroids have limitations in reliability and cytotoxicity. This study introduces a deep learning model for non-destructive, label-free viability estimation based on phase-contrast images, providing a cost-effective, high-throughput solution for continuous spheroid monitoring in microfluidics. Microfluidic technology facilitated the creation of a high-throughput cancer spheroid platform with approximately 12 000 spheroids per chip for drug screening. Validation involved tests with eight conventional chemotherapeutic drugs, revealing a strong correlation between viability assessed via LIVE/DEAD staining and phase-contrast morphology. Extending the model's application to novel compounds and cell lines not in the training dataset yielded promising results, implying the potential for a universal viability estimation model. Experiments with an alternative microscopy setup supported the model's transferability across different laboratories. Using this method, we also tracked the dynamic changes in spheroid viability during the course of drug administration. In summary, this research integrates a robust platform with high-throughput microfluidic cancer spheroid assays and deep learning-based viability estimation, with broad applicability to various cell lines, compounds, and research settings.

7 Tesla magnetic resonance spectroscopic imaging predicting IDH status and glioma grading

INTRODUCTION

With the application of high-resolution 3D 7 Tesla Magnetic Resonance Spectroscopy Imaging (MRSI) in high-grade gliomas, we previously identified intratumoral metabolic heterogeneities. In this study, we evaluated the potential of 3D 7 T-MRSI for the preoperative noninvasive classification of glioma grade and isocitrate dehydrogenase (IDH) status. We demonstrated that IDH mutation and glioma grade are detectable by ultra-high field (UHF) MRI. This technique might potentially optimize the perioperative management of glioma patients.

METHODS

We prospectively included 36 patients with WHO 2021 grade 2-4 gliomas (20 IDH mutated, 16 IDH wildtype). Our 7 T 3D MRSI sequence provided high-resolution metabolic maps (e.g., choline, creatine, glutamine, and glycine) of these patients' brains. We employed multivariate random forest and support vector machine models to voxels within a tumor segmentation, for classification of glioma grade and IDH mutation status.

RESULTS

Random forest analysis yielded an area under the curve (AUC) of 0.86 for multivariate IDH classification based on metabolic ratios. We distinguished high- and low-grade tumors by total choline (tCho) / total N-acetyl-aspartate (tNAA) ratio difference, yielding an AUC of 0.99. Tumor categorization based on other measured metabolic ratios provided comparable accuracy.

CONCLUSIONS

We successfully classified IDH mutation status and high- versus low-grade gliomas preoperatively based on 7 T MRSI and clinical tumor segmentation. With this approach, we demonstrated imaging based tumor marker predictions at least as accurate as comparable studies, highlighting the potential application of MRSI for pre-operative tumor classifications.

Recommended implementation of quantitative susceptibility mapping for clinical research in the brain: A consensus of the ISMRM electro-magnetic tissue properties study group

This article provides recommendations for implementing QSM for clinical brain research. It is a consensus of the International Society of Magnetic Resonance in Medicine, Electro-Magnetic Tissue Properties Study Group. While QSM technical development continues to advance rapidly, the current QSM methods have been demonstrated to be repeatable and reproducible for generating quantitative tissue magnetic susceptibility maps in the brain. However, the many QSM approaches available have generated a need in the neuroimaging community for guidelines on implementation. This article outlines considerations and implementation recommendations for QSM data acquisition, processing, analysis, and publication. We recommend that data be acquired using a monopolar 3D multi-echo gradient echo (GRE) sequence and that phase images be saved and exported in Digital Imaging and Communications in Medicine (DICOM) format and unwrapped using an exact unwrapping approach. Multi-echo images should be combined before background field removal, and a brain mask created using a brain extraction tool with the incorporation of phase-quality-based masking. Background fields within the brain mask should be removed using a technique based on SHARP or PDF, and the optimization approach to dipole inversion should be employed with a sparsity-based regularization. Susceptibility values should be measured relative to a specified reference, including the common reference region of the whole brain as a region of interest in the analysis. The minimum acquisition and processing details required when reporting QSM results are also provided. These recommendations should facilitate clinical QSM research and promote harmonized data acquisition, analysis, and reporting.

DF-QSM: Data Fidelity based Hybrid Approach for Improved Quantitative Susceptibility Mapping of the Brain

Quantitative Susceptibility Mapping (QSM) is an advanced magnetic resonance imaging (MRI) technique to quantify the magnetic susceptibility of the tissue under investigation. Deep learning methods have shown promising results in deconvolving the susceptibility distribution from the measured local field obtained from the MR phase. Although existing deep learning based QSM methods can produce high-quality reconstruction, they are highly biased toward training data distribution with less scope for generalizability. This work proposes a hybrid two-step reconstruction approach to improve deep learning based QSM reconstruction. The susceptibility map prediction obtained from the deep learning methods has been refined in the framework developed in this work to ensure consistency with the measured local field. The developed method was validated on existing deep learning and model-based deep learning methods for susceptibility mapping of the brain. The developed method resulted in improved reconstruction for MRI volumes obtained with different acquisition settings, including deep learning models trained on constrained (limited) data settings.

Proton Free Induction Decay MRSI at 7T in the Human Brain Using an Egg-Shaped Modified Rosette K-Space Trajectory

Purpose

1.1 Proton ( 1 H)-MRSI via spatial-spectral encoding poses high demands on gradient hardware at ultra-high fields and high-resolutions. Rosette trajectories help alleviate these problems, but at reduced SNR-efficiency due to their k-space densities not matching any desired k-space filter. We propose modified rosette trajectories, which more closely match a Hamming filter, and thereby improve SNR performance while still staying within gradient hardware limitations and without prolonging scan time.

Methods

1.2Analytical and synthetic simulations were validated with phantom and in vivo measurements at 7 T. The rosette and modified rosette trajectories were measured in five healthy volunteers in six minutes in a 2D slice in the brain. A 3D sequence was measured in one volunteer within 19 minutes. The SNR, linewidth, CRLBs, lipid contamination and data quality of the proposed modified rosette trajectory were compared to the rosette trajectory.

Results

1.3Using the modified rosette trajectories, an improved k-space weighting function was achieved resulting in an increase of up to 12% in SNR compared to rosette's dependent on the two additional trajectory parameters. Similar results were achieved for the theoretical SNR calculation based on k-space densities, as well as when using the pseudo-replica method for simulated, in-vivo and phantom data. The CRLBs improved slightly, but non-significantly for the modified rosette trajectories, while the linewidths and lipid contamination remained similar.

Conclusion

1.4By improving the rosette trajectory's shape, modified rosette trajectories achieved higher SNR at the same scan time and data quality.

Identifying virulence factors using graph transformer autoencoder with ESMFold-predicted structures

With the increasing resistance of bacterial pathogens to conventional antibiotics, antivirulence strategies targeting virulence factors (VFs) have become an effective new therapy for the treatment of pathogenic bacterial infections. Therefore, the identification and prediction of VFs can provide ideal candidate targets for the implementation of antivirulence strategies in treating infections caused by pathogenic bacteria. Currently, the existing computational models predominantly rely on the amino acid sequences of virulence proteins while overlooking structural information. Here, we propose a novel graph transformer autoencoder for VF identification (GTAE-VF), which utilizes ESMFold-predicted 3D structures and converts the VF identification problem into a graph-level prediction task. In an encoder-decoder framework, GTAE-VF adaptively learns both local and global information by integrating a graph convolutional network and a transformer to implement all-pair message passing, which can better capture long-range correlations and potential relationships. Extensive experiments on an independent test dataset demonstrate that GTAE-VF achieves reliable and robust prediction accuracy with an AUC of 0.963, which is consistently better than that of other structure-based and sequence-based approaches. We believe that GTAE-VF has the potential to emerge as a valuable tool for assessing VFs and devising antivirulence strategies.

Model resolution-based deconvolution for improved quantitative susceptibility mapping

Quantitative susceptibility mapping (QSM) utilizes the relationship between the measured local field and the unknown susceptibility map to perform dipole deconvolution. The aim of this work is to introduce and systematically evaluate the model resolution-based deconvolution for improved estimation of the susceptibility map obtained using the thresholded k-space division (TKD). A two-step approach has been proposed, wherein the first step involves the TKD susceptibility map computation and the second step involves the correction of this susceptibility map using the model-resolution matrix. The TKD-estimated susceptibility map can be expressed as the weighted average of the true susceptibility map, where the weights are determined by the rows of the model-resolution matrix, and hence a deconvolution of the TKD susceptibility map using the model-resolution matrix yields a better approximation to the true susceptibility map. The model resolution-based deconvolution is realized using closed-form, iterative, and sparsity-regularized implementations. The proposed approach was compared with L2 regularization, TKD, rescaled TKD in superfast dipole inversion, the modulated closed-form method, and iterative dipole inversion, as well as sparsity-regularized dipole inversion. It was observed that the proposed approach showed a substantial reduction in the streaking artifacts across 94 test volumes considered in this study. The proposed approach also showed better error reduction and edge preservation compared with other approaches. The proposed model resolution-based deconvolution compensates for the truncation of zero coefficients in the dipole kernel at the magic angle and hence provides a closer approximation to the true susceptibility map compared with other direct methods.

High-resolution proton metabolic mapping of the human brain at 7 T using free induction decay rosette spectroscopic imaging

Magnetic resonance spectroscopic imaging (MRSI) provides information about the spatial distribution of metabolites in the brain. These metabolite maps can be valuable in diagnosing central nervous system pathology. However, MRSI generally suffers from a long acquisition time, poor spatial resolution, and a low metabolite signal-to-noise ratio (SNR). Ultrahigh field strengths (≥ 7 T) can benefit MRSI with an improved SNR and allow high-resolution metabolic mapping. Non-Cartesian spatial-spectral encoding techniques, such as rosette spectroscopic imaging, can efficiently sample spatial and temporal domains, which significantly reduces the imaging time and enables high-resolution metabolic mapping in a clinically relevant scan time. In the current study, high-resolution (in-plane resolution of 2 × 2 mm2 ) mapping of proton (1 H) metabolites in the human brain at 7 T, is demonstrated. Five healthy subjects participated in the study. Using a time-efficient rosette trajectory and short TR/TE free induction decay MRSI, high-resolution maps of 1 H metabolites were obtained in a clinically relevant imaging time (6 min). Suppression of the water signal was achieved with an optimized water suppression enhanced through T1 effects approach and lipid removal was performed using L2 -regularization in the postprocessing. Spatial distributions of N-acetyl-aspartate, total choline, creatine, N-acetyl-aspartyl glutamate, myo-inositol, and glutamate were generated with Cramer-Rao lower bounds of less than 20%.

Precision Brain Morphometry Using Cluster Scanning

Measurement error limits the statistical power to detect group differences and longitudinal change in structural MRI morphometric measures (e.g., hippocampal volume, prefrontal thickness). Recent advances in scan acceleration enable extremely fast T1-weighted scans (~1 minute) to achieve morphometric errors that are close to the errors in longer traditional scans. As acceleration allows multiple scans to be acquired in rapid succession, it becomes possible to pool estimates to increase measurement precision, a strategy known as "cluster scanning." Here we explored brain morphometry using cluster scanning in a test-retest study of 40 individuals (12 younger adults, 18 cognitively unimpaired older adults, and 10 adults diagnosed with mild cognitive impairment or Alzheimer's Dementia). Morphometric errors from a single compressed sensing (CS) 1.0mm scan with 6x acceleration (CSx6) were, on average, 12% larger than a traditional scan using the Alzheimer's Disease Neuroimaging Initiative (ADNI) protocol. Pooled estimates from four clustered CSx6 acquisitions led to errors that were 34% smaller than ADNI despite having a shorter total acquisition time. Given a fixed amount of time, a gain in measurement precision can thus be achieved by acquiring multiple rapid scans instead of a single traditional scan. Errors were further reduced when estimates were pooled from eight CSx6 scans (51% smaller than ADNI). Neither pooling across a break nor pooling across multiple scan resolutions boosted this benefit. We discuss the potential of cluster scanning to improve morphometric precision, boost statistical power, and produce more sensitive disease progression biomarkers.

Metabolic Insights into Iron Deposition in Relapsing-Remitting Multiple Sclerosis via 7 T Magnetic Resonance Spectroscopic Imaging

OBJECTIVE

To investigate the metabolic pattern of different types of iron accumulation in multiple sclerosis (MS) lesions, and compare metabolic alterations within and at the periphery of lesions and newly emerging lesions in vivo according to iron deposition.

METHODS

7 T MR spectroscopic imaging and susceptibility-weighted imaging was performed in 31 patients with relapsing-remitting MS (16 female/15 male; mean age, 36.9 ± 10.3 years). Mean metabolic ratios of four neuro-metabolites were calculated for regions of interest (ROI) of normal appearing white matter (NAWM), "non-iron" (lesion without iron accumulation on SWI), and three distinct types of iron-laden lesions ("rim": distinct rim-shaped iron accumulation; "area": iron deposition across the entire lesions; "transition": transition between "area" and "rim" accumulation shape), and for lesion layers of "non-iron" and "rim" lesions. Furthermore, newly emerging "non-iron" and "iron" lesions were compared longitudinally, as measured before their appearance and one year later.

RESULTS

Thirty-nine of 75 iron-containing lesions showed no distinct paramagnetic rim. Of these, "area" lesions exhibited a 65% higher mIns/tNAA (p = 0.035) than "rim" lesions. Comparing lesion layers of both "non-iron" and "rim" lesions, a steeper metabolic gradient of mIns/tNAA ("non-iron" +15%, "rim" +40%) and tNAA/tCr ("non-iron" -15%, "rim" -35%) was found in "iron" lesions, with the lesion core showing +22% higher mIns/tNAA (p = 0.005) and -23% lower tNAA/tCr (p = 0.048) in "iron" compared to "non-iron" lesions. In newly emerging lesions, 18 of 39 showed iron accumulation, with the drop in tNAA/tCr after lesion formation remaining significantly lower compared to pre-lesional tissue over time in "iron" lesions (year 0: p = 0.013, year 1: p = 0.041) as opposed to "non-iron" lesions (year 0: p = 0.022, year 1: p = 0.231).

CONCLUSION

7 T MRSI allows in vivo characterization of different iron accumulation types each presenting with a distinct metabolic profile. Furthermore, the larger extent of neuronal damage in lesions with a distinct iron rim was reconfirmed via reduced tNAA/tCr concentrations, but with metabolic differences in lesion development between (non)-iron-containing lesions. This highlights the ability of MRSI to further investigate different types of iron accumulation and suggests possible implications for disease monitoring.

Reproducibility of 3D MRSI for imaging human brain glucose metabolism using direct (2H) and indirect (1H) detection of deuterium labeled compounds at 7T and clinical 3T

INTRODUCTION

Deuterium metabolic imaging (DMI) and quantitative exchange label turnover (QELT) are novel MR spectroscopy techniques for non-invasive imaging of human brain glucose and neurotransmitter metabolism with high clinical potential. Following oral or intravenous administration of non-ionizing [6,6'-2H2]-glucose, its uptake and synthesis of downstream metabolites can be mapped via direct or indirect detection of deuterium resonances using 2H MRSI (DMI) and 1H MRSI (QELT), respectively. The purpose of this study was to compare the dynamics of spatially resolved brain glucose metabolism, i.e., estimated concentration enrichment of deuterium labeled Glx (glutamate+glutamine) and Glc (glucose) acquired repeatedly in the same cohort of subjects using DMI at 7T and QELT at clinical 3T.

METHODS

Five volunteers (4 m/1f) were scanned in repeated sessions for 60 min after overnight fasting and 0.8 g/kg oral [6,6'-2H2]-glucose administration using time-resolved 3D 2H FID-MRSI with elliptical phase encoding at 7T and 3D 1H FID-MRSI with a non-Cartesian concentric ring trajectory readout at clinical 3T.

RESULTS

One hour after oral tracer administration regionally averaged deuterium labeled Glx4 concentrations and the dynamics were not significantly different over all participants between 7T 2H DMI and 3T 1H QELT data for GM (1.29±0.15 vs. 1.38±0.26 mM, p=0.65 & 21±3 vs. 26±3 µM/min, p=0.22) and WM (1.10±0.13 vs. 0.91±0.24 mM, p=0.34 & 19±2 vs. 17±3 µM/min, p=0.48). Also, the observed time constants of dynamic Glc6 data in GM (24±14 vs. 19±7 min, p=0.65) and WM (28±19 vs. 18±9 min, p=0.43) dominated regions showed no significant differences. Between individual 2H and 1H data points a weak to moderate negative correlation was observed for Glx4 concentrations in GM (r=-0.52, p<0.001), and WM (r=-0.3, p<0.001) dominated regions, while a strong negative correlation was observed for Glc6 data GM (r=-0.61, p<0.001) and WM (r=-0.70, p<0.001).

CONCLUSION

This study demonstrates that indirect detection of deuterium labeled compounds using 1H QELT MRSI at widely available clinical 3T without additional hardware is able to reproduce absolute concentration estimates of downstream glucose metabolites and the dynamics of glucose uptake compared to 2H DMI data acquired at 7T. This suggests significant potential for widespread application in clinical settings especially in environments with limited access to ultra-high field scanners and dedicated RF hardware.

Brain morphometry in older adults with and without dementia using extremely rapid structural scans

T1-weighted structural MRI is widely used to measure brain morphometry (e.g., cortical thickness and subcortical volumes). Accelerated scans as fast as one minute or less are now available but it is unclear if they are adequate for quantitative morphometry. Here we compared the measurement properties of a widely adopted 1.0 mm resolution scan from the Alzheimer's Disease Neuroimaging Initiative (ADNI = 5'12'') with two variants of highly accelerated 1.0 mm scans (compressed-sensing, CSx6 = 1'12''; and wave-controlled aliasing in parallel imaging, WAVEx9 = 1'09'') in a test-retest study of 37 older adults aged 54 to 86 (including 19 individuals diagnosed with a neurodegenerative dementia). Rapid scans produced highly reliable morphometric measures that largely matched the quality of morphometrics derived from the ADNI scan. Regions of lower reliability and relative divergence between ADNI and rapid scan alternatives tended to occur in midline regions and regions with susceptibility-induced artifacts. Critically, the rapid scans yielded morphometric measures similar to the ADNI scan in regions of high atrophy. The results converge to suggest that, for many current uses, extremely rapid scans can replace longer scans. As a final test, we explored the possibility of a 0'49'' 1.2 mm CSx6 structural scan, which also showed promise. Rapid structural scans may benefit MRI studies by shortening the scan session and reducing cost, minimizing opportunity for movement, creating room for additional scan sequences, and allowing for the repetition of structural scans to increase precision of the estimates.

Recommended Implementation of Quantitative Susceptibility Mapping for Clinical Research in The Brain: A Consensus of the ISMRM Electro-Magnetic Tissue Properties Study Group

This article provides recommendations for implementing quantitative susceptibility mapping (QSM) for clinical brain research. It is a consensus of the ISMRM Electro-Magnetic Tissue Properties Study Group. While QSM technical development continues to advance rapidly, the current QSM methods have been demonstrated to be repeatable and reproducible for generating quantitative tissue magnetic susceptibility maps in the brain. However, the many QSM approaches available give rise to the need in the neuroimaging community for guidelines on implementation. This article describes relevant considerations and provides specific implementation recommendations for all steps in QSM data acquisition, processing, analysis, and presentation in scientific publications. We recommend that data be acquired using a monopolar 3D multi-echo GRE sequence, that phase images be saved and exported in DICOM format and unwrapped using an exact unwrapping approach. Multi-echo images should be combined before background removal, and a brain mask created using a brain extraction tool with the incorporation of phase-quality-based masking. Background fields should be removed within the brain mask using a technique based on SHARP or PDF, and the optimization approach to dipole inversion should be employed with a sparsity-based regularization. Susceptibility values should be measured relative to a specified reference, including the common reference region of whole brain as a region of interest in the analysis, and QSM results should be reported with - as a minimum - the acquisition and processing specifications listed in the last section of the article. These recommendations should facilitate clinical QSM research and lead to increased harmonization in data acquisition, analysis, and reporting.

SLIPMAT: a pipeline for extracting tissue-specific spectral profiles from 1H MR spectroscopic imaging data

¹H Magnetic Resonance Spectroscopy (MRS) is an important non-invasive tool for measuring brain metabolism, with numerous applications in the neuroscientific and clinical domains. In this work we present a new analysis pipeline (SLIPMAT), designed to extract high-quality, tissue-specific, spectral profiles from MR spectroscopic imaging data (MRSI). Spectral decomposition is combined with spatially dependant frequency and phase correction to yield high SNR white and grey matter spectra without partial-volume contamination. A subsequent series of spectral processing steps are applied to reduce unwanted spectral variation, such as baseline correction and linewidth matching, before direct spectral analysis with machine learning and traditional statistical methods. The method is validated using a 2D semi-LASER MRSI sequence, with a 5-minute duration, from data acquired in triplicate across 8 healthy participants. Reliable spectral profiles are confirmed with principal component analysis, revealing the importance of total-choline and scyllo-inositol levels in distinguishing between individuals - in good agreement with our previous work. Furthermore, since the method allows the simultaneous measurement of metabolites in grey and white matter, we show the strong discriminative value of these metabolites in both tissue types for the first time. In conclusion, we present a novel and time efficient MRSI acquisition and processing pipeline, capable of detecting reliable neuro-metabolic differences between healthy individuals, and suitable for the sensitive neurometabolic profiling of in-vivo brain tissue.

Noninvasive 3-Dimensional 1 H-Magnetic Resonance Spectroscopic Imaging of Human Brain Glucose and Neurotransmitter Metabolism Using Deuterium Labeling at 3T : Feasibility and Interscanner Reproducibility

OBJECTIVES

Noninvasive, affordable, and reliable mapping of brain glucose metabolism is of critical interest for clinical research and routine application as metabolic impairment is linked to numerous pathologies, for example, cancer, dementia, and depression. A novel approach to map glucose metabolism noninvasively in the human brain has been presented recently on ultrahigh-field magnetic resonance (MR) scanners (≥7T) using indirect detection of deuterium-labeled glucose and downstream metabolites such as glutamate, glutamine, and lactate. The aim of this study was to demonstrate the feasibility to noninvasively detect deuterium-labeled downstream glucose metabolites indirectly in the human brain via 3-dimensional (3D) proton ( 1 H) MR spectroscopic imaging on a clinical 3T MR scanner without additional hardware.

MATERIALS AND METHODS

This prospective, institutional review board-approved study was performed in 7 healthy volunteers (mean age, 31 ± 4 years, 5 men/2 women) after obtaining written informed consent. After overnight fasting and oral deuterium-labeled glucose administration, 3D metabolic maps were acquired every ∼4 minutes with ∼0.24 mL isotropic spatial resolution using real-time motion-, shim-, and frequency-corrected echo-less 3D 1 H-MR spectroscopic Imaging on a clinical routine 3T MR system. To test the interscanner reproducibility of the method, subjects were remeasured on a similar 3T MR system. Time courses were analyzed using linear regression and nonparametric statistical tests. Deuterium-labeled glucose and downstream metabolites were detected indirectly via their respective signal decrease in dynamic 1 H MR spectra due to exchange of labeled and unlabeled molecules.

RESULTS

Sixty-five minutes after deuterium-labeled glucose administration, glutamate + glutamine (Glx) signal intensities decreased in gray/white matter (GM/WM) by -1.63 ± 0.3/-1.0 ± 0.3 mM (-13% ± 3%, P = 0.02/-11% ± 3%, P = 0.02), respectively. A moderate to strong negative correlation between Glx and time was observed in GM/WM ( r = -0.64, P < 0.001/ r = -0.54, P < 0.001), with 60% ± 18% ( P = 0.02) steeper slopes in GM versus WM, indicating faster metabolic activity. Other nonlabeled metabolites showed no significant changes. Excellent intrasubject repeatability was observed across scanners for static results at the beginning of the measurement (coefficient of variation 4% ± 4%), whereas differences were observed in individual Glx dynamics, presumably owing to physiological variation of glucose metabolism.

CONCLUSION

Our approach translates deuterium metabolic imaging to widely available clinical routine MR scanners without specialized hardware, offering a safe, affordable, and versatile (other substances than glucose can be labeled) approach for noninvasive imaging of glucose and neurotransmitter metabolism in the human brain.

Reproducibility of 3D MRSI for imaging human brain glucose metabolism using direct ( 2 H) and indirect ( 1 H) detection of deuterium labeled compounds at 7T and clinical 3T

Introduction

Deuterium metabolic imaging (DMI) and quantitative exchange label turnover (QELT) are novel MR spectroscopy techniques for non-invasive imaging of human brain glucose and neurotransmitter metabolism with high clinical potential. Following oral or intravenous administration of non-ionizing [6,6'- ² H ₂ ]-glucose, its uptake and synthesis of downstream metabolites can be mapped via direct or indirect detection of deuterium resonances using ² H MRSI (DMI) and ¹ H MRSI (QELT), respectively. The purpose of this study was to compare the dynamics of spatially resolved brain glucose metabolism, i.e., estimated concentration enrichment of deuterium labeled Glx (glutamate+glutamine) and Glc (glucose) acquired repeatedly in the same cohort of subjects using DMI at 7T and QELT at clinical 3T.

Methods

Five volunteers (4m/1f) were scanned in repeated sessions for 60 min after overnight fasting and 0.8g/kg oral [6,6'- ² H ₂ ]-glucose administration using time-resolved 3D ² H FID-MRSI with elliptical phase encoding at 7T and 3D ¹ H FID-MRSI with a non-Cartesian concentric ring trajectory readout at clinical 3T.

Results

One hour after oral tracer administration regionally averaged deuterium labeled Glx ₄ concentrations and the dynamics were not significantly different over all participants between 7T ² H DMI and 3T ¹ H QELT data for GM (1.29±0.15 vs. 1.38±0.26 mM, p=0.65 & 21±3 vs. 26±3 µM/min, p=0.22) and WM (1.10±0.13 vs. 0.91±0.24 mM, p=0.34 & 19±2 vs. 17±3 µM/min, p=0.48). Also, the observed time constants of dynamic Glc ₆ data in GM (24±14 vs. 19±7 min, p=0.65) and WM (28±19 vs. 18±9 min, p=0.43) dominated regions showed no significant differences. Between individual ² H and ¹ H data points a weak to moderate negative correlation was observed for Glx ₄ concentrations in GM (r=-0.52, p<0.001), and WM (r=-0.3, p<0.001) dominated regions, while a strong negative correlation was observed for Glc ₆ data GM (r=- 0.61, p<0.001) and WM (r=-0.70, p<0.001).

Conclusion

This study demonstrates that indirect detection of deuterium labeled compounds using ¹ H QELT MRSI at widely available clinical 3T without additional hardware is able to reproduce absolute concentration estimates of downstream glucose metabolites and the dynamics of glucose uptake compared to ² H DMI data acquired at 7T. This suggests significant potential for widespread application in clinical settings especially in environments with limited access to ultra-high field scanners and dedicated RF hardware.

AST-GIN: Attribute-Augmented Spatiotemporal Graph Informer Network for Electric Vehicle Charging Station Availability Forecasting

Electric Vehicle (EV) charging demand and charging station availability forecasting is one of the challenges in the intelligent transportation system. With accurate EV station availability prediction, suitable charging behaviors can be scheduled in advance to relieve range anxiety. Many existing deep learning methods have been proposed to address this issue; however, due to the complex road network structure and complex external factors, such as points of interest (POIs) and weather effects, many commonly used algorithms can only extract the historical usage information and do not consider the comprehensive influence of external factors. To enhance the prediction accuracy and interpretability, the Attribute-Augmented Spatiotemporal Graph Informer (AST-GIN) structure is proposed in this study by combining the Graph Convolutional Network (GCN) layer and the Informer layer to extract both the external and internal spatiotemporal dependence of relevant transportation data. The external factors are modeled as dynamic attributes by the attributeaugmented encoder for training. The AST-GIN model was tested on the data collected in Dundee City, and the experimental results showed the effectiveness of our model considering external factors' influence on various horizon settings compared with other baselines.

Vertical Nystagmus Recognition Based on Deep Learning

Vertical nystagmus is a common neuro-ophthalmic sign in vestibular medicine. Vertical nystagmus not only reflects the functional state of vertical semicircular canal but also reflects the effect of otoliths. Medical experts can take nystagmus symptoms as the key factor to determine the cause of dizziness. Traditional observation (visual observation conducted by medical experts) may be biased subjectively. Visual examination also requires medical experts to have enough experience to make an accurate diagnosis. With the development of science and technology, the detection system for nystagmus can be realized by using artificial intelligence technology. In this paper, a vertical nystagmus recognition method is proposed based on deep learning. This method is mainly composed of a dilated convolution layer module, a depthwise separable convolution module, a convolution attention module, a Bilstm-GRU module, etc. The average recognition accuracy of the proposed method is 91%. Using the same training dataset and test set, the recognition accuracy of this method for vertical nystagmus was 2% higher than other methods.

Lesion-Specific Metabolic Alterations in Relapsing-Remitting Multiple Sclerosis Via 7 T Magnetic Resonance Spectroscopic Imaging

BACKGROUND

Magnetic resonance spectroscopic imaging (MRSI) of the brain enables in vivo assessment of metabolic alterations in multiple sclerosis (MS). This provides complementary insights into lesion pathology that cannot be obtained via T1- and T2-weighted conventional magnetic resonance imaging (cMRI).

PURPOSE

The aims of this study were to assess focal metabolic alterations inside and at the periphery of lesions that are visible or invisible on cMRI, and to correlate their metabolic changes with T1 hypointensity and the distance of lesions to cortical gray matter (GM).

METHODS

A 7 T MRSI was performed on 51 patients with relapsing-remitting MS (30 female/21 male; mean age, 35.4 ± 9.9 years). Mean metabolic ratios were calculated for segmented regions of interest (ROIs) of normal-appearing white matter, white matter lesions, and focal regions of increased mIns/tNAA invisible on cMRI. A subgroup analysis was performed after subdividing based on T1 relaxation and distance to cortical GM. Metabolite ratios were correlated with T1 and compared between different layers around cMRI-visible lesions.

RESULTS

Focal regions of, on average, 2.8-fold higher mIns/tNAA than surrounding normal-appearing white matter and with an appearance similar to that of MS lesions were found, which were not visible on cMRI (ie, ~4% of metabolic hotspots). T1 relaxation was positively correlated with mIns/tNAA ( P ≤ 0.01), and negatively with tNAA/tCr ( P ≤ 0.01) and tCho/tCr ( P ≤ 0.01). mIns/tCr was increased outside lesions, whereas tNAA/tCr distributions resembled macroscopic tissue damage inside the lesions. mIns/tCr was -21% lower for lesions closer to cortical GM ( P ≤ 0.05).

CONCLUSIONS

7 T MRSI allows in vivo visualization of focal MS pathology not visible on cMRI and the assessment of metabolite levels in the lesion center, in the active lesion periphery and in cortical lesions. This demonstrated the potential of MRSI to image mIns as an early biomarker in lesion development.

Multi-echo dipole inversion for magnetic susceptibility mapping

PURPOSE

Reconstructing tissue magnetic susceptibility (QSM) from MRI phase data involves solving multiple consecutive ill-posed inverse problems such as phase unwrapping, background field removal, and field-to-source inversion. Multi-echo acquisitions present an additional challenge, as the magnetization field is typically computed from the multiple phase data prior to reconstructing the susceptibility map. Processing the multiple phase data introduces errors during the field estimation, violating assumptions of the subsequent inverse problems, manifesting as streaking artifacts in the susceptibility map. To address this challenge, we propose a multi-echo field-to-source forward model that forgoes the field estimation step. Moreover, we propose a fully general underestimation correction step to recover susceptibility sources that were regularized away during the field-to-source inversion.

METHODS

The multi-echo forward model and correction step were validated on the QSM Challenge 2.0 datasets and compared to the standard single field-to-source model in in vivo human brains using different types of deconvolution algorithms.

RESULTS

On the QSM Challenge 2.0 datasets the multi-echo forward model and correction step attain state-of-the-art results on all metrics by a wide margin. Experiments in in vivo brains show that the multi-echo model is in agreement with the single field-to-source model and that the proposed forward model and correction step can be used with any available dipole inversion method.

CONCLUSION

A multi-echo field-to-source forward model forgoes the need to fit multi-echo phase data and achieves state-of-the-art results on the QSM Challenge 2.0 data. Underestimated low-frequency susceptibility distributions can be partially recovered using a correction step.

Proton metabolic mapping of the brain at 7 T using a two-dimensional free induction decay-echo-planar spectroscopic imaging readout with lipid suppression

The increased signal-to-noise ratio (SNR) and chemical shift dispersion at high magnetic fields (≥7 T) have enabled neuro-metabolic imaging at high spatial resolutions. To avoid very long acquisition times with conventional magnetic resonance spectroscopic imaging (MRSI) phase-encoding schemes, solutions such as pulse-acquire or free induction decay (FID) sequences with short repetition time and inner volume selection methods with acceleration (echo-planar spectroscopic imaging [EPSI]), have been proposed. With the inner volume selection methods, limited spatial coverage of the brain and long echo times may still impede clinical implementation. FID-MRSI sequences benefit from a short echo time and have a high SNR per time unit; however, contamination from strong extra-cranial lipid signals remains a problem that can hinder correct metabolite quantification. L2-regularization can be applied to remove lipid signals in cases with high spatial resolution and accurate prior knowledge. In this work, we developed an accelerated two-dimensional (2D) FID-MRSI sequence using an echo-planar readout and investigated the performance of lipid suppression by L2-regularization, an external crusher coil, and the combination of these two methods to compare the resulting spectral quality in three subjects. The reduction factor of lipid suppression using the crusher coil alone varies from 2 to 7 in the lipid region of the brain boundary. For the combination of the two methods, the average lipid area inside the brain was reduced by 2% to 38% compared with that of unsuppressed lipids, depending on the subject's region of interest. 2D FID-EPSI with external lipid crushing and L2-regularization provides high in-plane coverage and is suitable for investigating brain metabolite distributions at high fields.

Test-retest reproducibility of human brain multi-slice 1 H FID-MRSI data at 9.4T after optimization of lipid regularization, macromolecular model, and spline baseline stiffness

PURPOSE

This study analyzes the effects of retrospective lipid suppression, a simulated macromolecular prior knowledge and different spline baseline stiffness values on 9.4T multi-slice proton FID-MRSI data spanning the whole cerebrum of human brain and the reproducibility of respective metabolite ratio to total creatine (/tCr) maps for 10 brain metabolites.

METHODS

Measurements were performed twice on 5 volunteers using a short TR and TE FID MRSI 2D sequence at 9.4T. The effects of retrospective lipid L2-regularization, macromolecular spectrum and different LCModel baseline flexibilities on SNR, FWHM, fitting residual, Cramér-Rao lower bound, and metabolite ratio maps were investigated. Intra-subject, inter-session coefficient of variation and the test-retest reproducibility of the mean metabolite ratios (/tCr) of each slice was calculated.

RESULTS

Transversal, sagittal, and coronal slices of many metabolite ratio maps correspond to the anatomically expected concentration relations in gray and white matter for the majority of the cerebrum when using a flexible baseline in LCModel fit. Results from the second measurements of the same subjects show that slice positioning and data quality correlate significantly to the first measurement. L2-regularization provided effective suppression of lipid-artifacts, but should be avoided if no artifacts are detected.

CONCLUSION

Reproducible concentration ratio maps (/tCr) for 4 metabolites (total choline, N-acetylaspartate, glutamate, and myoinositol) spanning the majority of the cerebrum and 6 metabolites (N-acetylaspartylglutamate, γ-aminobutyric acid, glutathione, taurine, glutamine, and aspartate) covering 32 mm in the upper part of the brain were acquired at 9.4T using multi-slice FID MRSI with retrospective lipid suppression, a macromolecular spectrum and a flexible LCModel baseline.

Short-term repeatability and long-term reproducibility of quantitative MR imaging biomarkers in a single centre longitudinal study

Quantitative imaging biomarkers (QIBs) can be defined as objective measures that are sensitive and specific to changes in tissue physiology. Provided the acquired QIBs are not affected by scanner changes, they could play an important role in disease diagnosis, prognosis, management, and treatment monitoring. The precision of selected QIBs was assessed from data collected on a 3-T scanner in four healthy participants over a 5-year period. Inevitable scanner changes and acquisition protocol revisions occurred during this time. Standard and custom processing pipelines were used to calculate regional brain volume, cortical thickness, T2, T2*, quantitative susceptibility, cerebral blood flow, axial, radial and mean diffusivity, peak width of skeletonized mean diffusivity, and fractional anisotropy from the acquired images. Coefficient of variation (CoV) and intra-class correlation (ICC) indices were determined in the short-term (i.e., repeatable over three acquisitions within 4 weeks) and in the long-term (i.e., reproducible over four acquisition sessions in 5 years). Precision indices varied based on acquisition technique, processing pipeline, and anatomical region. Good repeatability (average CoV=2.40% and ICC=0.78) and reproducibility (average CoV=8.86 % and ICC=0.72) were found over all QIBs. The best performance indices were obtained for diffusion derived biomarkers (CoV∼0.96% and ICCs=0.87); conversely, the poorest indices were found for the cerebral blood flow biomarker (CoV>10% and ICC<0.5). These results demonstrate that changes in protocol, along with hardware and software upgrades, did not affect the estimates of the selected biomarkers and their precision. Further characterization of the QIB is necessary to understand meaningful changes in the biomarkers in longitudinal studies of normal brain aging and translation to clinical research.

Sampling Possible Reconstructions of Undersampled Acquisitions in MR Imaging With a Deep Learned Prior

Undersampling the k-space during MR acquisitions saves time, however results in an ill-posed inversion problem, leading to an infinite set of images as possible solutions. Traditionally, this is tackled as a reconstruction problem by searching for a single "best" image out of this solution set according to some chosen regularization or prior. This approach, however, misses the possibility of other solutions and hence ignores the uncertainty in the inversion process. In this paper, we propose a method that instead returns multiple images which are possible under the acquisition model and the chosen prior to capture the uncertainty in the inversion process. To this end, we introduce a low dimensional latent space and model the posterior distribution of the latent vectors given the acquisition data in k-space, from which we can sample in the latent space and obtain the corresponding images. We use a variational autoencoder for the latent model and the Metropolis adjusted Langevin algorithm for the sampling. We evaluate our method on two datasets; with images from the Human Connectome Project and in-house measured multi-coil images. We compare to five alternative methods. Results indicate that the proposed method produces images that match the measured k-space data better than the alternatives, while showing realistic structural variability. Furthermore, in contrast to the compared methods, the proposed method yields higher uncertainty in the undersampled phase encoding direction, as expected.

7T HR FID-MRSI Compared to Amino Acid PET: Glutamine and Glycine as Promising Biomarkers in Brain Tumors

Simple Summary

Magnetic resonance spectroscopic imaging is an imaging method that can map the distribution of multiple biochemicals in the human brain in one scan. Using stronger magnetic fields, such as 7 Tesla, allows for higher resolution images and more biochemical maps. To test these results, we compared it to positron emission tomography, the established clinical standard for metabolic imaging. This comparison mainly looked at the overlap between regions with increased signal between both methods. We found that the molecules glutamine and glycine, only mappable at 7 Tesla, corresponded better to positron emission tomography than the commonly used choline.

Abstract

(1) Background: Recent developments in 7T magnetic resonance spectroscopic imaging (MRSI) made the acquisition of high-resolution metabolic images in clinically feasible measurement times possible. The amino acids glutamine (Gln) and glycine (Gly) were identified as potential neuro-oncological markers of importance. For the first time, we compared 7T MRSI to amino acid PET in a cohort of glioma patients. (2) Methods: In 24 patients, we co-registered 7T MRSI and routine PET and compared hotspot volumes of interest (VOI). We evaluated dice similarity coefficients (DSC), volume, center of intensity distance (CoI), median and threshold values for VOIs of PET and ratios of total choline (tCho), Gln, Gly, myo-inositol (Ins) to total N-acetylaspartate (tNAA) or total creatine (tCr). (3) Results: We found that Gln and Gly ratios generally resulted in a higher correspondence to PET than tCho. Using cutoffs of 1.6-times median values of a control region, DSCs to PET were 0.53 ± 0.36 for tCho/tNAA, 0.66 ± 0.40 for Gln/tNAA, 0.57 ± 0.36 for Gly/tNAA, and 0.38 ± 0.31 for Ins/tNAA. (4) Conclusions: Our 7T MRSI data corresponded better to PET than previous studies at lower fields. Our results for Gln and Gly highlight the importance of future research (e.g., using Gln PET tracers) into the role of both amino acids.

Exploring Hierarchical Auditory Representation via a Neural Encoding Model

By integrating hierarchical feature modeling of auditory information using deep neural networks (DNNs), recent functional magnetic resonance imaging (fMRI) encoding studies have revealed the hierarchical neural auditory representation in the superior temporal gyrus (STG). Most of these studies adopted supervised DNNs (e.g., for audio classification) to derive the hierarchical feature representation of external auditory stimuli. One possible limitation is that the extracted features could be biased toward discriminative features while ignoring general attributes shared by auditory information in multiple categories. Consequently, the hierarchy of neural acoustic processing revealed by the encoding model might be biased toward classification. In this study, we explored the hierarchical neural auditory representation via an fMRI encoding framework in which an unsupervised deep convolutional auto-encoder (DCAE) model was adopted to derive the hierarchical feature representations of the stimuli (naturalistic auditory excerpts in different categories) in fMRI acquisition. The experimental results showed that the neural representation of hierarchical auditory features is not limited to previously reported STG, but also involves the bilateral insula, ventral visual cortex, and thalamus. The current study may provide complementary evidence to understand the hierarchical auditory processing in the human brain.

Emerging methods and applications of ultra-high field MR spectroscopic imaging in the human brain

Magnetic Resonance Spectroscopic Imaging (MRSI) of the brain enables insights into the metabolic changes and fluxes in diseases such as tumors, multiple sclerosis, epilepsy, or hepatic encephalopathy, as well as insights into general brain functionality. However, the routine application of MRSI is mostly hampered by very low signal-to-noise ratios (SNR) due to the low concentrations of metabolites, about 10000 times lower than water. Furthermore, MRSI spectra have a dense information content with many overlapping metabolite resonances, especially for proton MRSI. MRI scanners at ultra-high field strengths, like 7 T or above, offer the opportunity to increase SNR, as well as the separation between resonances, thus promising to solve both challenges. Yet, MRSI at ultra-high field strengths is challenged by decreased B0- and B1-homogeneity, shorter T2 relaxation times, stronger chemical shift displacement errors, and aggravated lipid contamination. Therefore, to capitalize on the advantages of ultra-high field strengths, these challenges must be overcome. This review focuses on the challenges MRSI of the human brain faces at ultra-high field strength, as well as the possible applications to this date.

Extensive Brain Pathologic Alterations Detected with 7.0-T MR Spectroscopic Imaging Associated with Disability in Multiple Sclerosis

Background MR spectroscopic imaging (MRSI) allows in vivo assessment of brain metabolism and is of special interest in multiple sclerosis (MS), where morphologic MRI cannot depict major parts of disease activity. Purpose To evaluate the ability of 7.0-T MRSI to depict and visualize pathologic alterations in the normal-appearing white matter (NAWM) and cortical gray matter (CGM) in participants with MS and to investigate their relation to disability. Materials and Methods Free-induction decay MRSI was performed at 7.0 T. Participants with MS and age- and sex-matched healthy controls were recruited prospectively between January 2016 and December 2017. Metabolic ratios were obtained in white matter lesions, NAWM, and CGM regions. Subgroup analysis for MS-related disability based on Expanded Disability Status Scale (EDSS) scores was performed using analysis of covariance. Partial correlations were applied to explore associations between metabolic ratios and disability. Results Sixty-five participants with MS (mean age ± standard deviation, 34 years ± 9; 34 women) and 20 age- and sex-matched healthy controls (mean age, 32 years ± 7; 11 women) were evaluated. Higher signal intensity of myo-inositol (mI) with and without reduced signal intensity of N-acetylaspartate (NAA) was visible on metabolic images in the NAWM of participants with MS. A higher ratio of mI to total creatine (tCr) was observed in the NAWM of the centrum semiovale of all MS subgroups, including participants without disability (marginal mean ± standard error, healthy controls: 0.78 ± 0.04; EDSS 0-1: 0.86 ± 0.03 [P = .02]; EDSS 1.5-3: 0.95 ± 0.04 [P < .001]; EDSS ≥3.5: 0.94 ± 0.04 [P = .001]). A lower ratio of NAA to tCr was found in MS subgroups with disabilities, both in their NAWM (marginal mean ± standard error, healthy controls: 1.46 ± 0.04; EDSS 1.5-3: 1.33 ± 0.03 [P = .03]; EDSS ≥3.5: 1.30 ± 0.04 [P = .01]) and CGM (marginal mean ± standard error, healthy controls: 1.42 ± 0.05; EDSS ≥3.5: 1.23 ± 0.05 [P = .006]). mI/NAA correlated with EDSS (NAWM of centrum semiovale: r = 0.47, P < .001; parietal NAWM: r = 0.43, P = .002; frontal NAWM: r = 0.34, P = .01; frontal CGM: r = 0.37, P = .004). Conclusion MR spectroscopic imaging at 7.0 T allowed in vivo visualization of multiple sclerosis pathologic findings not visible at T1- or T2-weighted MRI. Metabolic abnormalities in the normal-appearing white matter and cortical gray matter were associated with disability. © RSNA, 2022 Online supplemental material is available for this article. See also the editorial by Barker in this issue.

BioSeq-BLM: a platform for analyzing DNA, RNA and protein sequences based on biological language models

In order to uncover the meanings of ‘book of life’, 155 different biological language models (BLMs) for DNA, RNA and protein sequence analysis are discussed in this study, which are able to extract the linguistic properties of ‘book of life’. We also extend the BLMs into a system called BioSeq-BLM for automatically representing and analyzing the sequence data. Experimental results show that the predictors generated by BioSeq-BLM achieve comparable or even obviously better performance than the exiting state-of-the-art predictors published in literatures, indicating that BioSeq-BLM will provide new approaches for biological sequence analysis based on natural language processing technologies, and contribute to the development of this very important field. In order to help the readers to use BioSeq-BLM for their own experiments, the corresponding web server and stand-alone package are established and released, which can be freely accessed at http://bliulab.net/BioSeq-BLM/.

QSM reconstruction challenge 2.0: Design and report of results

PURPOSE

The aim of the second quantitative susceptibility mapping (QSM) reconstruction challenge (Oct 2019, Seoul, Korea) was to test the accuracy of QSM dipole inversion algorithms in simulated brain data.

METHODS

A two-stage design was chosen for this challenge. The participants were provided with datasets of multi-echo gradient echo images synthesized from two realistic in silico head phantoms using an MR simulator. At the first stage, participants optimized QSM reconstructions without ground truth data available to mimic the clinical setting. At the second stage, ground truth data were provided for parameter optimization. Submissions were evaluated using eight numerical metrics and visual ratings.

RESULTS

A total of 98 reconstructions were submitted for stage 1 and 47 submissions for stage 2. Iterative methods had the best quantitative metric scores, followed by deep learning and direct inversion methods. Priors derived from magnitude data improved the metric scores. Algorithms based on iterative approaches and total variation (and its derivatives) produced the best overall results. The reported results and analysis pipelines have been made public to allow researchers to compare new methods to the current state of the art.

CONCLUSION

The synthetic data provide a consistent framework to test the accuracy and robustness of QSM algorithms in the presence of noise, calcifications and minor voxel dephasing effects. Total Variation-based algorithms produced the best results among all metrics. Future QSM challenges should assess whether this good performance with synthetic datasets translates to more realistic scenarios, where background fields and dipole-incompatible phase contributions are included.

Method for fast lipid reconstruction and removal processing in 1 H MRSI of the brain

PURPOSE

To develop a new rapid spatial filtering method for lipid removal, fast lipid reconstruction and removal processing (FLIP), which selectively isolates and removes interfering lipid signals from outside the brain in a full-FOV 2D MRSI and whole-brain 3D echo planar spectroscopic imaging (EPSI).

THEORY AND METHODS

FLIP uses regularized least-squares regression based on spatial prior information from MRI to selectively remove lipid signals originating from the scalp and measure the brain metabolite signals with minimum cross contamination. FLIP is a noniterative approach, thus allowing a rapid processing speed, and uses only spatial information without any spectral priors. The performance of FLIP was compared with the Papoulis-Gerchberg algorithm (PGA), Hankel singular value decomposition (HSVD), and fast image reconstruction with L2 regularization (L2).

RESULTS

FLIP in both 2D and 3D MRSI resulted in consistent metabolite quantification in a greater number of voxels with less concentration variation than other algorithms, demonstrating effective and robust lipid-removal performance. The percentage of voxels that met quality criteria with FLIP, PGA, HSVD, and L2 processing were 90%, 57%, 29%, and 42% in 2D MRSI, and 80%, 75%, 76%, and 74% in 3D EPSI, respectively. The quantification results of full-FOV MRSI using FLIP were comparable to those of volume-localized MRSI, while allowing significantly increased spatial coverage. FLIP performed the fastest in 2D MRSI.

CONCLUSION

FLIP is a new lipid-removal algorithm that promises fast and effective lipid removal with improved volume coverage in MRSI.

QSM reconstruction challenge 2.0: A realistic in silico head phantom for MRI data simulation and evaluation of susceptibility mapping procedures

PURPOSE

To create a realistic in silico head phantom for the second QSM reconstruction challenge and for future evaluations of processing algorithms for QSM.

METHODS

We created a digital whole-head tissue property phantom by segmenting and postprocessing high-resolution (0.64 mm isotropic), multiparametric MRI data acquired at 7 T from a healthy volunteer. We simulated the steady-state magnetization at 7 T using a Bloch simulator and mimicked a Cartesian sampling scheme through Fourier-based processing. Computer code for generating the phantom and performing the MR simulation was designed to facilitate flexible modifications of the phantom in the future, such as the inclusion of pathologies as well as the simulation of a wide range of acquisition protocols. Specifically, the following parameters and effects were implemented: TR and TE, voxel size, background fields, and RF phase biases. Diffusion-weighted imaging phantom data are provided, allowing future investigations of tissue-microstructure effects in phase and QSM algorithms.

RESULTS

The brain part of the phantom featured realistic morphology with spatial variations in relaxation and susceptibility values similar to the in vivo setting. We demonstrated some of the phantom's properties, including the possibility of generating phase data with nonlinear evolution over TE due to partial-volume effects or complex distributions of frequency shifts within the voxel.

CONCLUSION

The presented phantom and computer programs are publicly available and may serve as a ground truth in future assessments of the faithfulness of quantitative susceptibility reconstruction algorithms.

k-Space-based coil combination via geometric deep learning for reconstruction of non-Cartesian MRSI data

PURPOSE

State-of-the-art whole-brain MRSI with spatial-spectral encoding and multichannel acquisition generates huge amounts of data, which must be efficiently processed to stay within reasonable reconstruction times. Although coil combination significantly reduces the amount of data, currently it is performed in image space at the end of the reconstruction. This prolongs reconstruction times and increases RAM requirements. We propose an alternative k-space-based coil combination that uses geometric deep learning to combine MRSI data already in native non-Cartesian k-space.

METHODS

Twelve volunteers were scanned at a 3T MR scanner with a 20-channel head coil at 10 different positions with water-unsuppressed MRSI. At the eleventh position, water-suppressed MRSI data were acquired. Data of 7 volunteers were used to estimate sensitivity maps and form a base for simulating training data. A neural network was designed and trained to remove the effect of sensitivity profiles of the coil elements from the MRSI data. The water-suppressed MRSI data of the remaining volunteers were used to evaluate the performance of the new k-space-based coil combination relative to that of a conventional image-based alternative.

RESULTS

For both approaches, the resulting metabolic ratio maps were similar. The SNR of the k-space-based approach was comparable to the conventional approach in low SNR regions, but underperformed for high SNR. The Cramér-Rao lower bounds show the same trend. The analysis of the FWHM showed no difference between the two methods.

CONCLUSION

k-Space-based coil combination of MRSI data is feasible and reduces the amount of raw data immediately after their sampling.

Advanced magnetic resonance spectroscopic neuroimaging: Experts' consensus recommendations

Magnetic resonance spectroscopic imaging (MRSI) offers considerable promise for monitoring metabolic alterations associated with disease or injury; however, to date, these methods have not had a significant impact on clinical care, and their use remains largely confined to the research community and a limited number of clinical sites. The MRSI methods currently implemented on clinical MRI instruments have remained essentially unchanged for two decades, with only incremental improvements in sequence implementation. During this time, a number of technological developments have taken place that have already greatly benefited the quality of MRSI measurements within the research community and which promise to bring advanced MRSI studies to the point where the technique becomes a true imaging modality, while making the traditional review of individual spectra a secondary requirement. Furthermore, the increasing use of biomedical MR spectroscopy studies has indicated clinical areas where advanced MRSI methods can provide valuable information for clinical care. In light of this rapidly changing technological environment and growing understanding of the value of MRSI studies for biomedical studies, this article presents a consensus from a group of experts in the field that reviews the state-of-the-art for clinical proton MRSI studies of the human brain, recommends minimal standards for further development of vendor-provided MRSI implementations, and identifies areas which need further technical development.

Water and lipid suppression techniques for advanced 1 H MRS and MRSI of the human brain: Experts' consensus recommendations

The neurochemical information provided by proton magnetic resonance spectroscopy (MRS) or MR spectroscopic imaging (MRSI) can be severely compromised if strong signals originating from brain water and extracranial lipids are not properly suppressed. The authors of this paper present an overview of advanced water/lipid-suppression techniques and describe their advantages and disadvantages. Moreover, they provide recommendations for choosing the most appropriate techniques for proper use. Methods of water signal handling are primarily focused on the VAPOR technique and on MRS without water suppression (metabolite cycling). The section on lipid-suppression methods in MRSI is divided into three parts. First, lipid-suppression techniques that can be implemented on most clinical MR scanners (volume preselection, outer-volume suppression, selective lipid suppression) are described. Second, lipid-suppression techniques utilizing the combination of k-space filtering, high spatial resolutions and lipid regularization are presented. Finally, three promising new lipid-suppression techniques, which require special hardware (a multi-channel transmit system for dynamic B₁⁺ shimming, a dedicated second-order gradient system or an outer volume crusher coil) are introduced.

Accelerated radial echo-planar spectroscopic imaging using golden angle view-ordering and compressed-sensing reconstruction with total variation regularization

PURPOSE

To implement a novel, accelerated, 2D radial echo-planar spectroscopic imaging (REPSI) sequence using undersampled radial k-space trajectories and compressed-sensing reconstruction, and to compare results with those from an undersampled Cartesian spectroscopic sequence.

METHODS

The REPSI sequence was implemented using golden-angle view-ordering on a 3T MRI scanner. Radial and Cartesian echo-planar spectroscopic imaging (EPSI) data were acquired at six acceleration factors, each with time-equivalent scan durations, and reconstructed using compressed sensing with total variation regularization. Results from prospectively and retrospectively undersampled phantom and in vivo brain data were compared over estimated concentrations and Cramer-Rao lower-bound values, normalized RMS errors of reconstructed metabolite maps, and percent absolute differences between fully sampled and reconstructed spectroscopic images.

RESULTS

The REPSI method with compressed sensing is able to tolerate greater reductions in scan time compared with EPSI. The reconstruction and quantitation metrics (i.e., spectral normalized RMS error maps, metabolite map normalized RMS error values [e.g., for total N-acetyl asparate, REPSI = 9.4% vs EPSI = 16.3%; acceleration factor = 2.5], percent absolute difference maps, and concentration and Cramer-Rao lower-bound estimates) showed that accelerated REPSI can reduce the scan time by a factor of 2.5 while retaining image and quantitation quality.

CONCLUSION

Accelerated MRSI using undersampled radial echo-planar acquisitions provides greater reconstruction accuracy and more reliable quantitation for a range of acceleration factors compared with time-equivalent compressed-sensing reconstructions of undersampled Cartesian EPSI. Compared to the Cartesian approach, radial undersampling with compressed sensing could help reduce 2D spectroscopic imaging acquisition time, and offers a better trade-off between imaging speed and quality.

QSMART: Quantitative susceptibility mapping artifact reduction technique

PURPOSE

Quantitative susceptibility mapping (QSM) is a novel MR technique that allows mapping of tissue susceptibility values from MR phase images. QSM is an ill-conditioned inverse problem, and although several methods have been proposed in the field, in the presence of a wide range of susceptibility sources, streaking artifacts appear around high susceptibility regions and contaminate the whole QSM map. QSMART is a post-processing pipeline that uses two-stage parallel inversion to reduce the streaking artifacts and remove banding artifact at the cortical surface and around the vasculature.

METHOD

Tissue and vein susceptibility values were separately estimated by generating a mask of vasculature driven from the magnitude data using a Frangi filter. Spatially dependent filtering was used for the background field removal step and the two susceptibility estimates were combined in the final QSM map. QSMART was compared to RESHARP/iLSQR and V-SHARP/iLSQR inversion in a numerical phantom, 7T in vivo single and multiple-orientation scans, 9.4T ex vivo mouse data, and 4.7T in vivo rat brain with induced focal ischemia.

RESULTS

Spatially dependent filtering showed better suppression of phase artifacts near cortex compared to RESHARP and V-SHARP, while preserving voxels located within regions of interest without brain edge erosion. QSMART showed successful reduction of streaking artifacts as well as improved contrast between different brain tissues compared to the QSM maps obtained by RESHARP/iLSQR and V-SHARP/iLSQR.

CONCLUSION

QSMART can reduce QSM artifacts to enable more robust estimation of susceptibility values in vivo and ex vivo.

ECLIPSE utilizing gradient-modulated offset-independent adiabaticity (GOIA) pulses for highly selective human brain proton MRSI

A multitude of extracranial lipid suppression methods exist for proton MRSI acquisitions. Popular and emerging lipid suppression methods each have their inherent set of advantages and disadvantages related to the achievable level of lipid suppression, RF power deposition, insensitivity to B₁⁺ field and lipid T₁ heterogeneity, brain coverage, spatial selectivity, chemical shift displacement (CSD) errors and the reliability of spectroscopic data spanning the observed 0.9-4.7 ppm band. The utility of elliptical localization with pulsed second order fields (ECLIPSE) was previously demonstrated with a greater than 100-fold in extracranial lipid suppression and low power requirements utilizing 3 kHz bandwidth AFP pulses. Like all gradient-based localization methods, ECLIPSE is sensitive to CSD errors, resulting in a modified metabolic profile in edge-of-ROI voxels. In this work, ECLIPSE is extended with 15 kHz bandwidth second order gradient-modulated RF pulses based on the gradient offset-independent adiabaticity (GOIA) algorithm to greatly reduce CSD and improve spatial selectivity. An adiabatic double spin-echo ECLIPSE inner volume selection (TE = 45 ms) MRSI method and an ECLIPSE outer volume suppression (TE = 3.2 ms) FID-MRSI method were implemented. Both GOIA-ECLIPSE MRSI sequences provided artifact-free metabolite spectra in vivo, with a greater than 100-fold in lipid suppression and less than 2.6 mm in-plane CSD and less than 3.3 mm transition width for edge-of-ROI voxels, representing an ~5-fold improvement compared with the parent, nongradient-modulated method. Despite the 5-fold larger bandwidth, GOIA-ECLIPSE only required a 1.9-fold increase in RF power. The highly robust lipid suppression combined with low CSD and sharp ROI edge transitions make GOIA-ECLIPSE an attractive alternative to commonly employed lipid suppression methods. Furthermore, the low RF power deposition demonstrates that GOIA-ECLIPSE is very well suited for high field (≥3 T) MRSI applications.

SEPIA-Susceptibility mapping pipeline tool for phase images

Quantitative susceptibility mapping (QSM) is a physics-driven computational technique that has a high sensitivity in quantifying iron deposition based on MRI phase images. Furthermore, it has a unique ability to distinguish paramagnetic and diamagnetic contributions such as haemorrhage and calcification based on image contrast. These properties have contributed to a growing interest to use QSM not only in research but also in clinical applications. However, it is challenging to obtain high quality susceptibility map because of its ill-posed nature, especially for researchers who have less experience with QSM and the optimisation of its pipeline. In this paper, we present an open-source processing pipeline tool called SuscEptibility mapping PIpeline tool for phAse images (SEPIA) dedicated to the post-processing of MRI phase images and QSM. SEPIA connects various QSM toolboxes freely available in the field to offer greater flexibility in QSM processing. It also provides an interactive graphical user interface to construct and execute a QSM processing pipeline, simplifying the workflow in QSM research. The extendable design of SEPIA also allows developers to deploy their methods in the framework, providing a platform for developers and researchers to share and utilise the state-of-the-art methods in QSM.

High-resolution metabolic mapping of the cerebellum using 2D zoom magnetic resonance spectroscopic imaging

PURPOSE

The human cerebellum plays an important role in the functional activity of the cerebrum, ranging from motor to cognitive systems given its relaying role between the spinal cord and cerebrum. The cerebellum poses many challenges to Magnetic Resonance Spectroscopic Imaging (MRSI) due to its caudal location, susceptibility to physiological artifacts, and partial volume artifacts resulting from its complex anatomical structure. Thus, in the present study, we propose a high-resolution MRSI acquisition scheme for the cerebellum.

METHODS

A zoom or reduced field of view (rFOV) metabolite-cycled MRSI acquisition at 3 Tesla, with a grid of 48 × 48, was developed to achieve a nominal resolution of 62.5 μL. Single-slice rFOV MRSI data were acquired from the cerebellum of 5 healthy subjects with a nominal resolution of 2.5 × 2.5 × 10 mm³ in 9.6 min. Spectra were quantified using the LCModel package. A spatially unbiased atlas template of the cerebellum was used to analyze metabolite distributions in the cerebellum.

RESULTS

The superior quality of the achieved spectra-enabled generation of high-resolution metabolic maps of total N-acetylaspartate, total Creatine (tCr), total Choline (tCho), glutamate+glutamine, and myo-inositol, with Cramér-Rao lower bounds below 50%. A template-based regions of interest (ROI) analysis resulted in spatially dependent metabolite distributions in 9 ROIs. The group-averaged high-resolution metabolite maps across subjects increased the contrast-to-noise ratio between cerebellum regions.

CONCLUSION

These findings indicate that very high-resolution metabolite probing of the cerebellum is feasible using rFOV or zoomed MRSI at 3 Tesla.

Lipid-suppressed and tissue-fraction corrected metabolic distributions in human central brain structures using 2D 1 H magnetic resonance spectroscopic imaging at 7 T

INTRODUCTION

Magnetic resonance spectroscopic imaging (MRSI) has the potential to add a layer of understanding of the neurobiological mechanisms underlying brain diseases, disease progression, and treatment efficacy. Limitations related to metabolite fitting of low signal-to-noise ratios data, signal variations due to partial-volume effects, acquisition and extracranial lipid artifacts, along with clinically relevant aspects such as scan time constraints, are among the challenges associated with in vivo MRSI.

METHODS

The aim of this work was to address some of these factors and to develop an acquisition, reconstruction, and postprocessing pipeline to derive lipid-suppressed metabolite values of central brain structures based on free-induction decay measurements made using a 7 T MR scanner. Anatomical images were used to perform high-resolution (1 mm³ ) partial-volume correction to account for gray matter, white matter (WM), and cerebral-spinal fluid signal contributions. Implementation of automatic quality control thresholds and normalization of metabolic maps from 23 subjects to the Montreal Neurological Institute (MNI) standard atlas facilitated the creation of high-resolution average metabolite maps of several clinically relevant metabolites in central brain regions, while accounting for macromolecular distributions. Partial-volume correction improved the delineation of deep brain nuclei. We report average metabolite values including glutamate + glutamine (Glx), glycerophosphocholine, choline and phosphocholine (tCho), (phospo)creatine, myo-inositol and glycine (mI-Gly), glutathione, N-acetyl-aspartyl glutamate(and glutamine), and N-acetyl-aspartate in the basal ganglia, central WM (thalamic radiation, corpus callosum) as well as insular cortex and intracalcarine sulcus.

CONCLUSION

MNI-registered average metabolite maps facilitate group-based analysis, thus offering the possibility to mitigate uncertainty in variable MRSI data.

Improved body quantitative susceptibility mapping by using a variable-layer single-min-cut graph-cut for field-mapping

PURPOSE

To develop a robust algorithm for field-mapping in the presence of water-fat components, large B 0 field inhomogeneities and MR signal voids and to apply the developed method in body applications of quantitative susceptibility mapping (QSM).

METHODS

A framework solving the cost-function of the water-fat separation problem in a single-min-cut graph-cut based on the variable-layer graph construction concept was developed. The developed framework was applied to a numerical phantom enclosing an MR signal void, an air bubble experimental phantom, 14 large field of view (FOV) head/neck region in vivo scans and to 6 lumbar spine in vivo scans. Field-mapping and subsequent QSM results using the proposed algorithm were compared to results using an iterative graph-cut algorithm and a formerly proposed single-min-cut graph-cut.

RESULTS

The proposed method was shown to yield accurate field-map and susceptibility values in all simulation and in vivo datasets when compared to reference values (simulation) or literature values (in vivo). The proposed method showed improved field-map and susceptibility results compared to iterative graph-cut field-mapping especially in regions with low SNR, strong field-map variations and high R 2 ∗ values.

CONCLUSIONS

A single-min-cut graph-cut field-mapping method with a variable-layer construction was developed for field-mapping in body water-fat regions, improving quantitative susceptibility mapping particularly in areas close to MR signal voids.

Stromal microenvironment promoted infiltration in esophageal adenocarcinoma and squamous cell carcinoma: a multi-cohort gene-based analysis

The stromal microenvironment has been shown to affect the infiltration of esophageal carcinoma (ESCA), which is linked to prognosis. However, the complicated mechanism of how infiltration is influenced by the stromal microenvironment is not well-defined. In this study, a stromal activation classifier was established with ridge cox regression to calculate stroma scores for training (n = 182) and validation cohorts (n = 227) based on the stroma-related 32 hub genes identified by sequential bioinformatics algorithms. Patients with high stromal activation were associated with high T stage and poor prognosis in both esophagus adenocarcinoma and esophagus squamous cell carcinoma. Besides, comprehensive multi-omics analysis was used to outline stromal characterizations of 2 distinct stromal groups. Patients with activated tumor stoma showed high stromal cell infiltration (fibroblasts, endothelial cells, and monocyte macrophages), epithelial-mesenchymal transition, tumor angiogenesis and M2 macrophage polarization (CD163 and CD206). Tumor mutation burden of differential stromal groups was also depicted. In addition, a total of 6 stromal activation markers in ESCA were defined and involved in the function of carcinoma-associated fibroblasts that were crucial in the differentiation of distinct stromal characterizations. Based on these studies, a practical classifier for the stromal microenvironment was successfully proposed to predict the prognosis of ESCA patients.

High-resolution metabolic imaging of high-grade gliomas using 7T-CRT-FID-MRSI

OBJECTIVES

Successful neurosurgical intervention in gliomas depends on the precision of the preoperative definition of the tumor and its margins since a safe maximum resection translates into a better patient outcome. Metabolic high-resolution imaging might result in improved presurgical tumor characterization, and thus optimized glioma resection. To this end, we validated the performance of a fast high-resolution whole-brain 3D-magnetic resonance spectroscopic imaging (MRSI) method at 7T in a patient cohort of 23 high-grade gliomas (HGG).

MATERIALS AND METHODS

We preoperatively measured 23 patients with histologically verified HGGs (17 male, 8 female, age 53 ± 15) with an MRSI sequence based on concentric ring trajectories with a 64 × 64 × 39 measurement matrix, and a 3.4 × 3.4 × 3.4 mm3 nominal voxel volume in 15 min. Quantification used a basis-set of 17 components including N-acetyl-aspartate (NAA), total choline (tCho), total creatine (tCr), glutamate (Glu), glutamine (Gln), glycine (Gly) and 2-hydroxyglutarate (2HG). The resultant metabolic images were evaluated for their reliability as well as their quality and compared to spatially segmented tumor regions-of-interest (necrosis, contrast-enhanced, non-contrast enhanced + edema, peritumoral) based on clinical data and also compared to histopathology (e.g., grade, IDH-status).

RESULTS

Eighteen of the patient measurements were considered usable. In these patients, ten metabolites were quantified with acceptable quality. Gln, Gly, and tCho were increased and NAA and tCr decreased in nearly all tumor regions, with other metabolites such as serine, showing mixed trends. Overall, there was a reliable characterization of metabolic tumor areas. We also found heterogeneity in the metabolic images often continued into the peritumoral region. While 2HG could not be satisfyingly quantified, we found an increase of Glu in the contrast-enhancing region of IDH-wildtype HGGs and a decrease of Glu in IDH1-mutant HGGs.

CONCLUSIONS

We successfully demonstrated high-resolution 7T 3D-MRSI in HGG patients, showing metabolic differences between tumor regions and peritumoral tissue for multiple metabolites. Increases of tCho, Gln (related to tumor metabolism), Gly (related to tumor proliferation), as well as decreases in NAA, tCr, and others, corresponded very well to clinical tumor segmentation, but were more heterogeneous and often extended into the peritumoral region.

Terminology and concepts for the characterization of in vivo MR spectroscopy methods and MR spectra: Background and experts' consensus recommendations

With a 40-year history of use for in vivo studies, the terminology used to describe the methodology and results of magnetic resonance spectroscopy (MRS) has grown substantially and is not consistent in many aspects. Given the platform offered by this special issue on advanced MRS methodology, the authors decided to describe many of the implicated terms, to pinpoint differences in their meanings and to suggest specific uses or definitions. This work covers terms used to describe all aspects of MRS, starting from the description of the MR signal and its theoretical basis to acquisition methods, processing and to quantification procedures, as well as terms involved in describing results, for example, those used with regard to aspects of quality, reproducibility or indications of error. The descriptions of the meanings of such terms emerge from the descriptions of the basic concepts involved in MRS methods and examinations. This paper also includes specific suggestions for future use of terms where multiple conventions have emerged or coexisted in the past.

Ultra-High-Resolution in vitro MRI Study of Age-Related Brain Subcortical Susceptibility Alteration in Rhesus Monkeys at 9.4 T

Iron concentration in the brain has been suggested as a biomarker of pathologic neurodegeneration. However, the iron concentration changes in healthy aging as well. This study aimed to quantify the age-related changes in iron concentration in the gray matter of healthy rhesus monkeys using quantitative susceptibility mapping (QSM). Three-dimensional gradient-echo images of 16 female rhesus monkey brains aged between 2 and 26 years were acquired in vitro. The susceptibilities in the brain regions of the caudate nucleus (Cd), putamen (Pt), globus pallidus (Gp), and substantia nigra (Sn) were analyzed. The susceptibility varied across different brain regions, with higher levels in the Gp and Sn. Susceptibilities in all analyzed brain regions were linearly correlated with age, yet the plateau period as observed in human brains was absent. This is the first in vitro report of the age-related variability of susceptibility in the deep gray matter of rhesus monkey brains at 9.4 T, with an isotropic resolution of 150 μm. Awareness of age-related changes in susceptibility is vital for the establishment of a baseline to facilitate the differentiation of pathologic neurodegeneration from healthy aging in non-human primate studies.