Spectral decomposition for resolving partial volume effects in MRSI

A comparison of non-negative matrix underapproximation methods for the decomposition of magnetic resonance spectroscopy data from human brain tumors

Magnetic resonance spectroscopy (MRS) is an MR technique that provides information about the biochemistry of tissues in a noninvasive way. MRS has been widely used for the study of brain tumors, both preoperatively and during follow-up. In this study, we investigated the performance of a range of variants of unsupervised matrix factorization methods of the non-negative matrix underapproximation (NMU) family, namely, sparse NMU, global NMU, and recursive NMU, and compared them with convex non-negative matrix factorization (C-NMF), which has previously shown a good performance on brain tumor diagnostic support problems using MRS data. The purpose of the investigation was 2-fold: first, to ascertain the differences among the sources extracted by these methods; and second, to compare the influence of each method in the diagnostic accuracy of the classification of brain tumors, using them as feature extractors. We discovered that, first, NMU variants found meaningful sources in terms of biological interpretability, but representing parts of the spectrum, in contrast to C-NMF; and second, that NMU methods achieved better classification accuracy than C-NMF for the classification tasks when one class was not meningioma.

Spectroscopic MRI-Based Biomarkers Predict Survival for Newly Diagnosed Glioblastoma in a Clinical Trial

Despite aggressive treatment, glioblastoma has a poor prognosis due to its infiltrative nature. Spectroscopic MRI-measured brain metabolites, particularly the choline to N-acetylaspartate ratio (Cho/NAA), better characterizes the extent of tumor infiltration. In a previous pilot trial (NCT03137888), brain regions with Cho/NAA ≥ 2x normal were treated with high-dose radiation for newly diagnosed glioblastoma patients. This report is a secondary analysis of that trial where spectroscopic MRI-based biomarkers are evaluated for how they correlate with progression-free and overall survival (PFS/OS). Subgroups were created within the cohort based on pre-radiation treatment (pre-RT) median cutoff volumes of residual enhancement (2.1 cc) and metabolically abnormal volumes used for treatment (19.2 cc). We generated Kaplan-Meier PFS/OS curves and compared these curves via the log-rank test between subgroups. For the subgroups stratified by metabolic abnormality, statistically significant differences were observed for PFS (p = 0.019) and OS (p = 0.020). Stratification by residual enhancement did not lead to observable differences in the OS (p = 0.373) or PFS (p = 0.286) curves. This retrospective analysis shows that patients with lower post-surgical Cho/NAA volumes had significantly superior survival outcomes, while residual enhancement, which guides high-dose radiation in standard treatment, had little significance in PFS/OS. This suggests that the infiltrating, non-enhancing component of glioblastoma is an important factor in patient outcomes and should be treated accordingly.

Exploratory Multisite MR Spectroscopic Imaging Shows White Matter Neuroaxonal Loss Associated with Complications of Type 1 Diabetes in Children

BACKGROUND AND PURPOSE

Type 1 diabetes affects over 200,000 children in the United States and is associated with an increased risk of cognitive dysfunction. Prior single-site, single-voxel MRS case reports and studies have identified associations between reduced NAA/Cr, a marker of neuroaxonal loss, and type 1 diabetes. However, NAA/Cr differences among children with various disease complications or across different brain tissues remain unclear. To better understand this phenomenon and the role of MRS in characterizing it, we conducted a multisite pilot study.

MATERIALS AND METHODS

In 25 children, 6-14 years of age, with type 1 diabetes across 3 sites, we acquired T1WI and axial 2D MRSI along with phantom studies to calibrate scanner effects. We quantified tissue-weighted NAA/Cr in WM and deep GM and modeled them against study covariates.

RESULTS

We found that MRSI differentiated WM and deep GM by NAA/Cr on the individual level. On the population level, we found significant negative associations of WM NAA/Cr with chronic hyperglycemia quantified by hemoglobin A1c (P < .005) and a history of diabetic ketoacidosis at disease onset (P < .05). We found a statistical interaction (P < .05) between A1c and ketoacidosis, suggesting that neuroaxonal loss from ketoacidosis may outweigh that from poor glucose control. These associations were not present in deep GM.

CONCLUSIONS

Our pilot study suggests that MRSI differentiates GM and WM by NAA/Cr in this population, disease complications may lead to neuroaxonal loss in WM in children, and deeper investigation is warranted to further untangle how diabetic ketoacidosis and chronic hyperglycemia affect brain health and cognition in type 1 diabetes.

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.

Applying a Radiation Therapy Volume Analysis Pipeline to Determine the Utility of Spectroscopic MRI-Guided Adaptive Radiation Therapy for Glioblastoma

Accurate radiation therapy (RT) targeting is crucial for glioblastoma treatment but may be challenging using clinical imaging alone due to the infiltrative nature of glioblastomas. Precise targeting by whole-brain spectroscopic MRI, which maps tumor metabolites including choline (Cho) and N-acetylaspartate (NAA), can quantify early treatment-induced molecular changes that other traditional modalities cannot measure. We developed a pipeline to determine how spectroscopic MRI changes during early RT are associated with patient outcomes to provide insight into the utility of adaptive RT planning. Data were obtained from a study (NCT03137888) where glioblastoma patients received high-dose RT guided by the pre-RT Cho/NAA twice normal (Cho/NAA ≥ 2x) volume, and received spectroscopic MRI scans pre- and mid-RT. Overlap statistics between pre- and mid-RT scans were used to quantify metabolic activity changes after two weeks of RT. Log-rank tests were used to quantify the relationship between imaging metrics and patient overall and progression-free survival (OS/PFS). Patients with lower Jaccard/Dice coefficients had longer PFS (p = 0.045 for both), and patients with lower Jaccard/Dice coefficients had higher OS trending towards significance (p = 0.060 for both). Cho/NAA ≥ 2x volumes changed significantly during early RT, putting healthy tissue at risk of irradiation, and warranting further study into using adaptive RT planning.

Spectroscopic MRI-Guided Proton Therapy in Non-Enhancing Pediatric High-Grade Glioma

Radiation therapy (RT) is a critical part of definitive therapy for pediatric high-grade glioma (pHGG). RT is designed to treat residual tumor defined on conventional MRI (cMRI), though pHGG lesions may be ill-characterized on standard imaging. Spectroscopic MRI (sMRI) measures endogenous metabolite concentrations in the brain, and Choline (Cho)/N-acetylaspartate (NAA) ratio is a highly sensitive biomarker for metabolically active tumor. We provide a preliminary report of our study introducing a novel treatment approach of whole brain sMRI-guided proton therapy for pHGG. An observational cohort (c1 = 10 patients) receives standard of care RT; a therapeutic cohort (c2 = 15 patients) receives sMRI-guided proton RT. All patients undergo cMRI and sMRI, a high-resolution 3D whole-brain echo-planar spectroscopic imaging (EPSI) sequence (interpolated resolution of 12 µL) prior to RT and at several follow-up timepoints integrated into diagnostic scans. Treatment volumes are defined by cMRI for c1 and by cMRI and Cho/NAA ≥ 2x for c2. A longitudinal imaging database is used to quantify changes in lesion and metabolite volumes. Four subjects have been enrolled (c1 = 1/c2 = 3) with sMRI imaging follow-up of 4-18 months. Preliminary data suggest sMRI improves identification of pHGG infiltration based on abnormal metabolic activity, and using proton therapy to target sMRI-defined high-risk regions is safe and feasible.

Replicability of proton MR spectroscopic imaging findings in mild traumatic brain injury: Implications for clinical applications

PURPOSE

Proton magnetic resonance spectroscopy (1H MRS) offers biomarkers of metabolic damage after mild traumatic brain injury (mTBI), but a lack of replicability studies hampers clinical translation. In a conceptual replication study design, the results reported in four previous publications were used as the hypotheses (H1-H7), specifically: abnormalities in patients are diffuse (H1), confined to white matter (WM) (H2), comprise low N-acetyl-aspartate (NAA) levels and normal choline (Cho), creatine (Cr) and myo-inositol (mI) (H3), and correlate with clinical outcome (H4); additionally, a lack of findings in regional subcortical WM (H5) and deep gray matter (GM) structures (H6), except for higher mI in patients' putamen (H7).

METHODS

26 mTBI patients (20 female, age 36.5 ± 12.5 [mean ± standard deviation] years), within two months from injury and 21 age-, sex-, and education-matched healthy controls were scanned at 3 Tesla with 3D echo-planar spectroscopic imaging. To test H1-H3, global analysis using linear regression was used to obtain metabolite levels of GM and WM in each brain lobe. For H4, patients were stratified into non-recovered and recovered subgroups using the Glasgow Outcome Scale Extended. To test H5-H7, regional analysis using spectral averaging estimated metabolite levels in four GM and six WM structures segmented from T1-weighted MRI. The Mann-Whitney U test and weighted least squares analysis of covariance were used to examine mean group differences in metabolite levels between all patients and all controls (H1-H3, H5-H7), and between recovered and non-recovered patients and their respectively matched controls (H4). Replicability was defined as the support or failure to support the null hypotheses in accordance with the content of H1-H7, and was further evaluated using percent differences, coefficients of variation, and effect size (Cohen's d).

RESULTS

Patients' occipital lobe WM Cho and Cr levels were 6.0% and 4.6% higher than controls', respectively (Cho, d = 0.37, p = 0.04; Cr, d = 0.63, p = 0.03). The same findings, i.e., higher patients' occipital lobe WM Cho and Cr (both p = 0.01), but with larger percent differences (Cho, 8.6%; Cr, 6.3%) and effect sizes (Cho, d = 0.52; Cr, d = 0.88) were found in the comparison of non-recovered patients to their matched controls. For the lobar WM Cho and Cr comparisons without statistical significance (frontal, parietal, temporal), unidirectional effect sizes were observed (Cho, d = 0.07 - 0.37; Cr, d = 0.27 - 0.63). No differences were found in any metabolite in any lobe in the comparison between recovered patients and their matched controls. In the regional analyses, no differences in metabolite levels were found in any GM or WM region, but all WM regions (posterior, frontal, corona radiata, and the genu, body, and splenium of the corpus callosum) exhibited unidirectional effect sizes for Cho and Cr (Cho, d = 0.03 - 0.34; Cr, d = 0.16 - 0.51).

CONCLUSIONS

We replicated findings of diffuse WM injury, which correlated with clinical outcome (supporting H1-H2, H4). These findings, however, were among the glial markers Cho and Cr, not the neuronal marker NAA (not supporting H3). No differences were found in regional GM and WM metabolite levels (supporting H5-H6), nor in putaminal mI (not supporting H7). Unidirectional effect sizes of higher patients' Cho and Cr within all WM analyses suggest widespread injury, and are in line with the conclusion from the previous publications, i.e., that detection of WM injury may be more dependent upon sensitivity of the 1H MRS technique than on the selection of specific regions. The findings lend further support to the corollary that clinic-ready 1H MRS biomarkers for mTBI may best be achieved by using high signal-to-noise-ratio single-voxels placed anywhere within WM. The biochemical signature of the injury, however, may differ and therefore absolute levels, rather than ratios may be preferred. Future replication efforts should further test the generalizability of these findings.

A multi-institutional pilot clinical trial of spectroscopic MRI-guided radiation dose escalation for newly diagnosed glioblastoma

Background

Glioblastomas (GBMs) are aggressive brain tumors despite radiation therapy (RT) to 60 Gy and temozolomide (TMZ). Spectroscopic magnetic resonance imaging (sMRI), which measures levels of specific brain metabolites, can delineate regions at high risk for GBM recurrence not visualized on contrast-enhanced (CE) MRI. We conducted a clinical trial to assess the feasibility, safety, and efficacy of sMRI-guided RT dose escalation to 75 Gy for newly diagnosed GBMs.

Methods

Our pilot trial (NCT03137888) enrolled patients at 3 institutions (Emory University, University of Miami, Johns Hopkins University) from September 2017 to June 2019. For RT, standard tumor volumes based on T2-FLAIR and T1w-CE MRIs with margins were treated in 30 fractions to 50.1 and 60 Gy, respectively. An additional high-risk volume based on residual CE tumor and Cho/NAA (on sMRI) ≥2× normal was treated to 75 Gy. Survival curves were generated by the Kaplan-Meier method. Toxicities were assessed according to CTCAE v4.0.

Results

Thirty patients were treated in the study. The median age was 59 years. 30% were MGMT promoter hypermethylated; 7% harbored IDH1 mutation. With a median follow-up of 21.4 months for censored patients, median overall survival (OS) and progression-free survival were 23.0 and 16.6 months, respectively. This regimen appeared well-tolerated with 70% of grade 3 or greater toxicity ascribed to TMZ and 23% occurring at least 1 year after RT.

Conclusion

Dose-escalated RT to 75 Gy guided by sMRI appears feasible and safe for patients with newly diagnosed GBMs. OS outcome is promising and warrants additional testing. Based on these results, a randomized phase II trial is in development.

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.

The clinical utility of proton magnetic resonance spectroscopy in traumatic brain injury: recommendations from the ENIGMA MRS working group

Proton (1H) magnetic resonance spectroscopy provides a non-invasive and quantitative measure of brain metabolites. Traumatic brain injury impacts cerebral metabolism and a number of research groups have successfully used this technique as a biomarker of injury and/or outcome in both pediatric and adult TBI populations. However, this technique is underutilized, with studies being performed primarily at centers with access to MR research support. In this paper we present a technical introduction to the acquisition and analysis of in vivo 1H magnetic resonance spectroscopy and review 1H magnetic resonance spectroscopy findings in different injury populations. In addition, we propose a basic 1H magnetic resonance spectroscopy data acquisition scheme (Supplemental Information) that can be added to any imaging protocol, regardless of clinical magnetic resonance platform. We outline a number of considerations for study design as a way of encouraging the use of 1H magnetic resonance spectroscopy in the study of traumatic brain injury, as well as recommendations to improve data harmonization across groups already using this technique.

Bolus infusion scheme for the adjustment of steady state [11C]Flumazenil levels in the grey matter and in the blood plasma for neuroreceptor imaging

The use of hybrid PET/MR imaging facilitates the simultaneous investigation of challenge-related changes in ligand binding to neuroreceptors using PET, while concurrently measuring neuroactivation or blood flow with MRI. Having attained a steady state of the PET radiotracer using a bolus-infusion protocol, it is possible to observe alterations in ligand neuroreceptor binding through changes in distribution volumes. Here, we present an iterative procedure for establishing an administration scheme to obtain steady state [¹¹C]flumazenil concentrations in grey matter in the human brain. In order to achieve a steady state in the shortest possible time, the bolus infusion ratio from a previous examination was adapted to fit the subsequent examination. 17 male volunteers were included in the study. Boli and infusions with different weightings were given to the subjects and were characterised by kbol values from 74 ​min down to 42 ​min. Metabolite analysis was used to ascertain the value of unmetabolised flumazenil in the plasma, and PET imaging was used to assess its binding in the grey matter. The flumazenil time-activity curves (TACs) in the brain were decomposed into activity contributions from pure grey and white matter and analysed for 12 ​vol of interest (VOIs). The curves highlighted a large variability in metabolic rates between the subjects, with kbol ​= ​54.3 ​min being a reliable value to provide flumazenil equilibrium conditions in the majority of the VOIs and cases. The distribution volume of flumazenil in all 12 VOIs was determined.