Comparison of Multivendor Single-Voxel MR Spectroscopy Data Acquired in Healthy Brain at 26 Sites

Model-based frequency-and-phase correction of 1H MRS data with 2D linear-combination modeling

Purpose

Retrospective frequency-and-phase correction (FPC) methods attempt to remove frequency-and-phase variations between transients to improve the quality of the averaged MR spectrum. However, traditional FPC methods like spectral registration struggle at low SNR. Here, we propose a method that directly integrates FPC into a two-dimensional linear-combination model (2D-LCM) of individual transients ('model-based FPC'). We investigated how model-based FPC performs compared to the traditional approach, i.e., spectral registration followed by 1D-LCM in estimating frequency-and-phase drifts and, consequentially, metabolite level estimates.

Methods

We created synthetic in-vivo-like 64-transient short-TE sLASER datasets with 100 noise realizations at 5 SNR levels and added randomly sampled frequency and phase variations. We then used this synthetic dataset to compare the performance of 2D-LCM with the traditional approach (spectral registration, averaging, then 1D-LCM). Outcome measures were the frequency/phase/amplitude errors, the standard deviation of those ground-truth errors, and amplitude Cramér Rao Lower Bounds (CRLBs). We further tested the proposed method on publicly available in-vivo short-TE PRESS data.

Results

2D-LCM estimates (and accounts for) frequency-and-phase variations directly from uncorrected data with equivalent or better fidelity than the conventional approach. Furthermore, 2D-LCM metabolite amplitude estimates were at least as accurate, precise, and certain as the conventionally derived estimates. 2D-LCM estimation of frequency and phase correction and amplitudes performed substantially better at low-to-very-low SNR.

Conclusion

Model-based FPC with 2D linear-combination modeling is feasible and has great potential to improve metabolite level estimation for conventional and dynamic MRS data, especially for low-SNR conditions, e.g., long TEs or strong diffusion weighting.

Linewidth-related bias in modelled concentration estimates from GABA-edited 1H-MRS

J-difference-edited MRS is widely used to study GABA in the human brain. Editing for low-concentration target molecules (such as GABA) typically exhibits lower signal-to-noise ratio (SNR) than conventional non-edited MRS, varying with acquisition region, volume and duration. Moreover, spectral lineshape may be influenced by age-, pathology-, or brain-region-specific effects of metabolite T2, or by task-related blood-oxygen level dependent (BOLD) changes in functional MRS contexts. Differences in both SNR and lineshape may have systematic effects on concentration estimates derived from spectral modelling. The present study characterises the impact of lineshape and SNR on GABA+ estimates from different modelling algorithms: FSL-MRS, Gannet, LCModel, Osprey, spant and Tarquin. Publicly available multi-site GABA-edited data (222 healthy subjects from 20 sites; conventional MEGA-PRESS editing; TE = 68 ms) were pre-processed with a standardised pipeline, then filtered to apply controlled levels of Lorentzian and Gaussian linebroadening and SNR reduction. Increased Lorentzian linewidth was associated with a 2-5% decrease in GABA+ estimates per Hz, observed consistently (albeit to varying degrees) across datasets and most algorithms. Weaker, often opposing effects were observed for Gaussian linebroadening. Variations are likely caused by differing baseline parametrization and lineshape constraints between models. Effects of linewidth on other metabolites (e.g., Glx and tCr) varied, suggesting that a linewidth confound may persist after scaling to an internal reference. These findings indicate a potentially significant confound for studies where linewidth may differ systematically between groups or experimental conditions, e.g. due to T2 differences between brain regions, age, or pathology, or varying T2* due to BOLD-related changes. We conclude that linewidth effects need to be rigorously considered during experimental design and data processing, for example by incorporating linewidth into statistical analysis of modelling outcomes or development of appropriate lineshape matching algorithms.

Simultaneous multi-region detection of GABA+ and Glx using 3D spatially resolved SLOW-editing and EPSI-readout at 7T

GABA+ and Glx (glutamate and glutamine) are widely studied metabolites, yet the commonly used magnetic resonance spectroscopy (MRS) techniques have significant limitations, including sensitivity to B0 and B1+-inhomogeneities, limited bandwidth of MEGA-pulses, high SAR which is accentuated at 7T. To address these limitations, we propose SLOW-EPSI method, employing a large 3D MRSI coverage and achieving a high resolution down to 0.26 ml. Simulation results demonstrate the robustness of SLOW-editing for both GABA+ and Glx against B0 and B1+-inhomogeneities within the range of [-0.3, +0.3] ppm and [40 %, 250 %], respectively. Two protocols, both utilizing a 70 mm thick FOV slab, were employed to target distinct brain regions in vivo, differentiated by their orientation: transverse and tilted. Protocol 1 (n = 11) encompassed 5 locations (cortical gray matter, white matter, frontal lobe, parietal lobe, and cingulate gyrus). Protocol 2 (n = 5) involved 9 locations (cortical gray matter, white matter, frontal lobe, occipital lobe, cingulate gyrus, caudate nucleus, hippocampus, putamen, and inferior thalamus). Quantitative analysis of GABA+ and Glx was conducted in a stepwise manner. First, B1+/B1--inhomogeneities were corrected using water reference data. Next, GABA+ and Glx values were calculated employing spectral fitting. Finally, the GABA+ level for each selected region was compared to the global Glx within the same subject, generating the GABA+/Glx_global ratio. Our findings from two protocols indicate that the GABA+/Glx_global level in cortical gray matter was approximately 16 % higher than in white matter. Elevated GABA+/Glx_global levels acquired with protocol 2 were observed in specific regions such as the caudate nucleus (0.118±0.067), putamen (0.108±0.023), thalamus (0.092±0.036), and occipital cortex (0.091±0.010), when compared to the cortical gray matter (0.079±0.012). Overall, our results highlight the effectiveness of SLOW-EPSI as a robust and efficient technique for accurate measurements of GABA+ and Glx at 7T. In contrast to previous SVS and 2D-MRSI based editing sequences with which only one or a limited number of brain regions can be measured simultaneously, the method presented here measures GABA+ and Glx from any brain area and any arbitrarily shaped volume that can be flexibly selected after the examination. Quantification of GABA+ and Glx across multiple brain regions through spectral fitting is achievable with a 9-minute acquisition. Additionally, acquisition times of 18-27 min (GABA+) and 9-18 min (Glx) are required to generate 3D maps, which are constructed using Gaussian fitting and peak integration.

X-linked creatine transporter (SLC6A8) deficiency in females: Difficult to recognize, but a potentially treatable disease

Creatine transporter deficiency (CTD), caused by pathogenic variants in SLC6A8, is the second most common cause of X-linked intellectual disability. Symptoms include intellectual disability, epilepsy, and behavioral disorders and are caused by reduced cerebral creatine levels. Targeted treatment with oral supplementation is available, however the treatment efficacy is still being investigated. There are clinical and theoretical indications that heterozygous females with CTD respond better to supplementation treatment than hemizygous males. Unfortunately, heterozygous females with CTD often have more subtle and uncharacteristic clinical and biochemical phenotypes, rendering diagnosis more difficult. We report a new female case who presented with learning disabilities and seizures. After determining the diagnosis with molecular genetic testing confirmed by proton magnetic resonance spectroscopy (1H-MRS), the patient was treated with supplementation treatment including creatine, arginine, and glycine. After 28 months of treatment, the patient showed prominent clinical improvement and increased creatine levels in the brain. Furthermore, we provide a review of the 32 female cases reported in the current literature including a description of phenotypes, genotypes, diagnostic approaches, and effects of supplementation treatment. Based on this, we find that supplementation treatment should be tested in heterozygous female patients with CTD, and a prospective treatment underlines the importance of diagnosing these patients. The diagnosis should be suspected in a broad clinical spectrum of female patients and can only be made by molecular genetic testing. 1H-MRS of cerebral creatine levels is essential for establishing the diagnosis in females, and especially valuable when assessing variants of unknown significance.

Quantification of short echo time MRS signals with improved version of QUantitation based on quantum ESTimation algorithm

Magnetic resonance spectroscopy offers information about metabolite changes in the organism, which can be used in diagnosis. While short echo time proton spectra exhibit more distinguishable metabolites compared with proton spectra acquired with long echo times, their quantification (and providing estimates of metabolite concentrations) is more challenging. They are hampered by a background signal, which originates mainly from macromolecules (MM) and mobile lipids. An improved version of the quantification algorithm QUantitation based on quantum ESTimation (QUEST), with MM prior knowledge (QUEST-MM), dedicated to proton signals and invoking appropriate prior knowledge on MM, is proposed and tested. From a single acquisition, it enables better metabolite quantification, automatic estimation of the background, and additional automatic quantification of MM components, thus improving its applicability in the clinic. The proposed algorithm may facilitate studies that involve patients with pathological MM in the brain. QUEST-MM and three QUEST-based strategies for quantifying short echo time signals are compared in terms of bias-variance trade-off and Cramér-Rao lower bound estimates. The performances of the methods are evaluated through extensive Monte Carlo studies. In particular, the histograms of the metabolite and MM amplitude distributions demonstrate the performances of the estimators. They showed that QUEST-MM works better than QUEST (Subtract approach) and is a good alternative to QUEST when measured MM signal is unavailable or unsuitable. Quantification with QUEST-MM is shown for 1 H in vivo rat brain signals obtained with the SPECIAL pulse sequence at 9.4 T, and human brain signals obtained, respectively, with STEAM at 4 T and PRESS at 3 T. QUEST-MM is implemented in jMRUI and will be available for public use from version 7.1.

Continuous Automated Analysis Workflow for MRS Studies

Magnetic resonance spectroscopy (MRS) can non-invasively measure levels of endogenous metabolites in living tissue and is of great interest to neuroscience and clinical research. To this day, MRS data analysis workflows differ substantially between groups, frequently requiring many manual steps to be performed on individual datasets, e.g., data renaming/sorting, manual execution of analysis scripts, and manual assessment of success/failure. Manual analysis practices are a substantial barrier to wider uptake of MRS. They also increase the likelihood of human error and prevent deployment of MRS at large scale. Here, we demonstrate an end-to-end workflow for fully automated data uptake, processing, and quality review.The proposed continuous automated MRS analysis workflow integrates several recent innovations in MRS data and file storage conventions. They are efficiently deployed by a directory monitoring service that automatically triggers the following steps upon arrival of a new raw MRS dataset in a project folder: (1) conversion from proprietary manufacturer file formats into the universal format NIfTI-MRS; (2) consistent file system organization according to the data accumulation logic standard BIDS-MRS; (3) executing a command-line executable of our open-source end-to-end analysis software Osprey; (4) e-mail delivery of a quality control summary report for all analysis steps.The automated architecture successfully completed for a demonstration dataset. The only manual step required was to copy a raw data folder into a monitored directory.Continuous automated analysis of MRS data can reduce the burden of manual data analysis and quality control, particularly for non-expert users and multi-center or large-scale studies and offers considerable economic advantages.

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.

Brain metabolites measured with magnetic resonance spectroscopy in pediatric concussion and orthopedic injury: An Advancing Concussion Assessment in Pediatrics (A-CAP) study

Millions of children sustain a concussion annually. Concussion disrupts cellular signaling and neural pathways within the brain but the resulting metabolic disruptions are not well characterized. Magnetic resonance spectroscopy (MRS) can examine key brain metabolites (e.g., N-acetyl Aspartate (tNAA), glutamate (Glx), creatine (tCr), choline (tCho), and myo-Inositol (mI)) to better understand these disruptions. In this study, we used MRS to examine differences in brain metabolites between children and adolescents with concussion versus orthopedic injury. Children and adolescents with concussion (n = 361) or orthopedic injury (OI) (n = 184) aged 8 to 17 years were recruited from five emergency departments across Canada. MRS data were collected from the left dorsolateral prefrontal cortex (L-DLPFC) using point resolved spectroscopy (PRESS) at 3 T at a mean of 12 days post-injury (median 10 days post-injury, range 2-33 days). Univariate analyses for each metabolite found no statistically significant metabolite differences between groups. Within each analysis, several covariates were statistically significant. Follow-up analyses designed to account for possible confounding factors including age, site, scanner, vendor, time since injury, and tissue type (and interactions as appropriate) did not find any metabolite group differences. In the largest sample of pediatric concussion studied with MRS to date, we found no metabolite differences between concussion and OI groups in the L-DLPFC. We suggest that at 2 weeks post-injury in a general pediatric concussion population, brain metabolites in the L-DLPFC are not specifically affected by brain injury.

Feasibility and implications of using subject-specific macromolecular spectra to model short echo time magnetic resonance spectroscopy data

Expert consensus recommends linear-combination modeling (LCM) of 1 H MR spectra with sequence-specific simulated metabolite basis function and experimentally derived macromolecular (MM) basis functions. Measured MM basis functions are usually derived from metabolite-nulled spectra averaged across a small cohort. The use of subject-specific instead of cohort-averaged measured MM basis functions has not been studied widely. Furthermore, measured MM basis functions are not widely available to non-expert users, who commonly rely on parameterized MM signals internally simulated by LCM software. To investigate the impact of the choice of MM modeling, this study, therefore, compares metabolite level estimates between different MM modeling strategies (cohort-mean measured; subject-specific measured; parameterized) in a lifespan cohort and characterizes its impact on metabolite-age associations. 100 conventional (TE = 30 ms) and metabolite-nulled (TI = 650 ms) PRESS datasets, acquired from the medial parietal lobe in a lifespan cohort (20-70 years of age), were analyzed in Osprey. Short-TE spectra were modeled in Osprey using six different strategies to consider the MM baseline. Fully tissue- and relaxation-corrected metabolite levels were compared between MM strategies. Model performance was evaluated by model residuals, the Akaike information criterion (AIC), and the impact on metabolite-age associations. The choice of MM strategy had a significant impact on the mean metabolite level estimates and no major impact on variance. Correlation analysis revealed moderate-to-strong agreement between different MM strategies (r > 0.6). The lowest relative model residuals and AIC values were found for the cohort-mean measured MM. Metabolite-age associations were consistently found for two major singlet signals (total creatine (tCr])and total choline (tCho)) for all MM strategies; however, findings for metabolites that are less distinguishable from the background signals associations depended on the MM strategy. A variance partition analysis indicated that up to 44% of the total variance was related to the choice of MM strategy. Additionally, the variance partition analysis reproduced the metabolite-age association for tCr and tCho found in the simpler correlation analysis. In summary, the inclusion of a single high signal-to-noise ratio MM basis function (cohort-mean) in the short-TE LCM leads to more lower model residuals and AIC values compared with MM strategies with more degrees of freedom (Gaussian parametrization) or subject-specific MM information. Integration of multiple LCM analyses into a single statistical model potentially allows to identify the robustness in the detection of underlying effects (e.g., metabolite vs. age), reduces algorithm-based bias, and estimates algorithm-related variance.

COVID-19 in people with neurofibromatosis 1, neurofibromatosis 2, or schwannomatosis

PURPOSE

People with pre-existing conditions may be more susceptible to severe COVID-19 when infected by SARS-CoV-2. The relative risk and severity of SARS-CoV-2 infection in people with rare diseases such as neurofibromatosis type 1 (NF1), neurofibromatosis type 2 (NF2), or schwannomatosis (SWN) is unknown.

METHODS

We investigated the proportions of people with NF1, NF2, or SWN in the National COVID Cohort Collaborative (N3C) electronic health record data set who had a positive test result for SARS-CoV-2 or COVID-19.

RESULTS

The cohort sizes in N3C were 2501 (NF1), 665 (NF2), and 762 (SWN). We compared these with N3C cohorts of patients with other rare diseases (98-9844 individuals) and the general non-NF population of 5.6 million. The site- and age-adjusted proportion of people with NF1, NF2, or SWN who had a positive test result for SARS-CoV-2 or COVID-19 (collectively termed positive cases) was not significantly higher than in individuals without NF or other selected rare diseases. There were no severe outcomes reported in the NF2 or SWN cohorts. The proportion of patients experiencing severe outcomes was no greater for people with NF1 than in cohorts with other rare diseases or the general population.

CONCLUSION

Having NF1, NF2, or SWN does not appear to increase the risk of being SARS-CoV-2 positive or of being a patient with COVID-19 or of developing severe complications from SARS-CoV-2.

More than one-half of the variance in in vivo proton MR spectroscopy metabolite estimates is common to all metabolites

The present study characterized associations among brain metabolite levels, applying bivariate and multivariate (i.e., factor analysis) statistical methods to total creatine (tCr)-referenced estimates of the major Point RESolved Spectroscopy (PRESS) proton MR spectroscopy (¹ H-MRS) metabolites (i.e., total NAA/tCr, total choline/tCr, myo-inositol/tCr, glutamate + glutamine/tCr) acquired at 3 T from medial parietal lobe in a large (n = 299), well-characterized international cohort of healthy volunteers. Results supported the hypothesis that ¹ H-MRS-measured metabolite estimates are moderately intercorrelated (M_r  = 0.42, SD_r  = 0.11, ps < 0.001), with more than one-half (i.e., 57%) of the total variability in metabolite estimates explained by a single common factor. Older age was significantly associated with lower levels of the identified common metabolite variance (CMV) factor (β = -0.09, p = 0.048), despite not being associated with levels of any individual metabolite. Holding CMV factor levels constant, females had significantly lower levels of total choline (i.e., unique metabolite variance; β = -0.19, p < 0.001), mirroring significant bivariate correlations between sex and total choline reported previously. Supplementary analysis of water-referenced metabolite estimates (i.e., including tCr/water) demonstrated lower, although still substantial, intercorrelations among metabolites, with 37% of total metabolite variance explained by a single common factor. If replicated, these results would suggest that applied ¹ H-MRS researchers shift their analytical framework from examining bivariate associations between individual metabolites and specialty-dependent (e.g., clinical, research) variables of interest (e.g., using t-tests) to examining multivariable (i.e., covariate) associations between multiple metabolites and specialty-dependent variables of interest (e.g., using multiple regression).

Harmonization of multi-scanner in vivo magnetic resonance spectroscopy: ENIGMA consortium task group considerations

Magnetic resonance spectroscopy is a powerful, non-invasive, quantitative imaging technique that allows for the measurement of brain metabolites that has demonstrated utility in diagnosing and characterizing a broad range of neurological diseases. Its impact, however, has been limited due to small sample sizes and methodological variability in addition to intrinsic limitations of the method itself such as its sensitivity to motion. The lack of standardization from a data acquisition and data processing perspective makes it difficult to pool multiple studies and/or conduct multisite studies that are necessary for supporting clinically relevant findings. Based on the experience of the ENIGMA MRS work group and a review of the literature, this manuscript provides an overview of the current state of MRS data harmonization. Key factors that need to be taken into consideration when conducting both retrospective and prospective studies are described. These include (1) MRS acquisition issues such as pulse sequence, RF and B0 calibrations, echo time, and SNR; (2) data processing issues such as pre-processing steps, modeling, and quantitation; and (3) biological factors such as voxel location, age, sex, and pathology. Various approaches to MRS data harmonization are then described including meta-analysis, mega-analysis, linear modeling, ComBat and artificial intelligence approaches. The goal is to provide both novice and experienced readers with the necessary knowledge for conducting MRS data harmonization studies.

Comparison of different approaches to manage multi-site magnetic resonance spectroscopy clinical data analysis

Introduction

The effects caused by differences in data acquisition can be substantial and may impact data interpretation in multi-site/scanner studies using magnetic resonance spectroscopy (MRS). Given the increasing use of multi-site studies, a better understanding of how to account for different scanners is needed. Using data from a concussion population, we compare ComBat harmonization with different statistical methods in controlling for site, vendor, and scanner as covariates to determine how to best control for multi-site data.

Methods

The data for the current study included 545 MRS datasets to measure tNAA, tCr, tCho, Glx, and mI to study the pediatric concussion acquired across five sites, six scanners, and two different MRI vendors. For each metabolite, the site and vendor were accounted for in seven different models of general linear models (GLM) or mixed-effects models while testing for group differences between the concussion and orthopedic injury. Models 1 and 2 controlled for vendor and site. Models 3 and 4 controlled for scanner. Models 5 and 6 controlled for site applied to data harmonized by vendor using ComBat. Model 7 controlled for scanner applied to data harmonized by scanner using ComBat. All the models controlled for age and sex as covariates.

Results

Models 1 and 2, controlling for site and vendor, showed no significant group effect in any metabolites, but the vendor and site were significant factors in the GLM. Model 3, which included a scanner, showed a significant group effect for tNAA and tCho, and the scanner was a significant factor. Model 4, controlling for the scanner, did not show a group effect in the mixed model. The data harmonized by the vendor using ComBat (Models 5 and 6) had no significant group effect in both the GLM and mixed models. Lastly, the data harmonized by the scanner using ComBat (Model 7) showed no significant group effect. The individual site data suggest there were no group differences.

Conclusion

Using data from a large clinical concussion population, different analysis techniques to control for site, vendor, and scanner in MRS data yielded different results. The findings support the use of ComBat harmonization for clinical MRS data, as it removes the site and vendor effects.

Revisiting the use of proton magnetic resonance spectroscopy in distinguishing between primary and secondary malignant tumors of the central nervous system

BACKGROUND AND PURPOSE

Conventional magnetic resonance images (MRI) has limitations in distinguishing primary from secondary brain tumors. Proton magnetic resonance spectroscopy (1H-MRS) allows evaluation of the concentration of metabolites in a brain lesion and, hence, better characterization of the tumor. Considering that an accurate diagnosis determines the choice of treatment, our purpose was to assess the usefulness of spectroscopy data for differentiating between primary and secondary brain neoplasms.

MATERIALS AND METHODS

We undertook a retrospective analysis of 61 MRI and 1H-MRS images of patients with histologically confirmed tumors (30 primary tumors and 31 metastatic tumors). The metabolite ratios of Cho/Cr and NAA/Cr at short TE were determined from spectroscopic curves, with a single voxel positioned in the enhancing tumor. Additional variables analyzed along with the metabolites, like as age and gender, allowed the construction of a logistic regression model to predict the tumor's nature. The statistical analysis was done using the R software (version 4.0.3 R Core Team, 2020).

RESULTS

The mean NAA/Cr and Cho/Cr ratios were higher in secondary tumors, with a good correlation between NAA/Cr and Cho/Cr (r = 0.61). The mean age of patients with primary tumors was lower than for secondary tumors (43.9 vs 55.9, respectively). Receiver operating characteristic analysis yielded a cut-off value of 0.4 for the NAA/Cr ratio with an accuracy of 73.8%, a sensitivity of 73.3% and a specificity of 74.2% in predicting metastatic tumors.

CONCLUSION

The model was reasonable in predicting the nature of the tumor and provides an additional tool for analyzing brain tumors.

Comparison of seven modelling algorithms for γ-aminobutyric acid-edited proton magnetic resonance spectroscopy

Edited MRS sequences are widely used for studying γ-aminobutyric acid (GABA) in the human brain. Several algorithms are available for modelling these data, deriving metabolite concentration estimates through peak fitting or a linear combination of basis spectra. The present study compares seven such algorithms, using data obtained in a large multisite study. GABA-edited (GABA+, TE = 68 ms MEGA-PRESS) data from 222 subjects at 20 sites were processed via a standardised pipeline, before modelling with FSL-MRS, Gannet, AMARES, QUEST, LCModel, Osprey and Tarquin, using standardised vendor-specific basis sets (for GE, Philips and Siemens) where appropriate. After referencing metabolite estimates (to water or creatine), systematic differences in scale were observed between datasets acquired on different vendors' hardware, presenting across algorithms. Scale differences across algorithms were also observed. Using the correlation between metabolite estimates and voxel tissue fraction as a benchmark, most algorithms were found to be similarly effective in detecting differences in GABA+. An interclass correlation across all algorithms showed single-rater consistency for GABA+ estimates of around 0.38, indicating moderate agreement. Upon inclusion of a basis set component explicitly modelling the macromolecule signal underlying the observed 3.0 ppm GABA peaks, single-rater consistency improved to 0.44. Correlation between discrete pairs of algorithms varied, and was concerningly weak in some cases. Our findings highlight the need for consensus on appropriate modelling parameters across different algorithms, and for detailed reporting of the parameters adopted in individual studies to ensure reproducibility and meaningful comparison of outcomes between different studies.

Harmonization of multi-site MRS data with ComBat

Magnetic resonance spectroscopy (MRS) is a non-invasive neuroimaging technique used to measure brain chemistry in vivo and has been used to study the healthy brain as well as neuropathology in numerous neurological disorders. The number of multi-site studies using MRS are increasing; however, non-biological variability introduced during data collection across multiple sites, such as differences in scanner vendors and site-specific acquisition implementations for MRS, can obscure detection of biological effects of interest. ComBat is a data harmonization technique that can remove non-biological sources of variance in multisite studies. It has been validated for use with structural and functional MRI metrics but not for MRS metabolites. This study investigated the validity of using ComBat to harmonize MRS metabolites for vendor and site differences. Analyses were performed using data acquired across 20 sites and included edited MRS for GABA+ (N=218) and macromolecule-suppressed GABA data (N=209), as well as standard PRESS data to quantify NAA, creatine, choline, and glutamate (N=190). ComBat harmonization successfully mitigated vendor and site differences for all metabolites of interest. Moreover, significant associations were detected between sex and choline levels and between age and glutamate and GABA+ levels that were not detectable prior to harmonization, confirming the importance of removing site and vendor effects in multi-site data. In conclusion, ComBat harmonization can be successfully applied to MRS data in multi-site MRS studies.

High-Grade Glioma Treatment Response Monitoring Biomarkers: A Position Statement on the Evidence Supporting the Use of Advanced MRI Techniques in the Clinic, and the Latest Bench-to-Bedside Developments. Part 2: Spectroscopy, Chemical Exchange Saturation, Multiparametric Imaging, and Radiomics

Objective

To summarize evidence for use of advanced MRI techniques as monitoring biomarkers in the clinic, and to highlight the latest bench-to-bedside developments.

Methods

The current evidence regarding the potential for monitoring biomarkers was reviewed and individual modalities of metabolism and/or chemical composition imaging discussed. Perfusion, permeability, and microstructure imaging were similarly analyzed in Part 1 of this two-part review article and are valuable reading as background to this article. We appraise the clinic readiness of all the individual modalities and consider methodologies involving machine learning (radiomics) and the combination of MRI approaches (multiparametric imaging).

Results

The biochemical composition of high-grade gliomas is markedly different from healthy brain tissue. Magnetic resonance spectroscopy allows the simultaneous acquisition of an array of metabolic alterations, with choline-based ratios appearing to be consistently discriminatory in treatment response assessment, although challenges remain despite this being a mature technique. Promising directions relate to ultra-high field strengths, 2-hydroxyglutarate analysis, and the use of non-proton nuclei. Labile protons on endogenous proteins can be selectively targeted with chemical exchange saturation transfer to give high resolution images. The body of evidence for clinical application of amide proton transfer imaging has been building for a decade, but more evidence is required to confirm chemical exchange saturation transfer use as a monitoring biomarker. Multiparametric methodologies, including the incorporation of nuclear medicine techniques, combine probes measuring different tumor properties. Although potentially synergistic, the limitations of each individual modality also can be compounded, particularly in the absence of standardization. Machine learning requires large datasets with high-quality annotation; there is currently low-level evidence for monitoring biomarker clinical application.

Conclusion

Advanced MRI techniques show huge promise in treatment response assessment. The clinical readiness analysis highlights that most monitoring biomarkers require standardized international consensus guidelines, with more facilitation regarding technique implementation and reporting in the clinic.

Reproducibility of 7-T brain spectroscopy using an ultrashort echo time STimulated Echo Acquisition Mode sequence and automated voxel repositioning

Establishing the reproducibility of brain MRS is important for clinical studies so that researchers can evaluate changes in metabolites due to treatment or the course of a disease and better understand the brain in healthy and disordered states. Prior 7-T MRS reproducibility studies using the stimulated echo acquisition mode (STEAM) sequence have focused on the anterior cingulate cortex or posterior cingulate cortex and precuneus. The purpose of this study was to evaluate the reproducibility of metabolite measurements in the dorsolateral prefrontal cortex (DLPFC) using an ultrashort echo time (TE) STEAM sequence and automated voxel repositioning. Spectra were acquired during two scan sessions from nine subjects using the AutoAlign method for voxel repositioning. Reproducibility was evaluated with coefficients of variation (CVs) and percentage differences. The mean intrasubject CVs were less than 6% for the major metabolites glutamate, N-acetylaspartate, total creatine, total choline, and myo-inositol. The mean CVs were less than 20% for the smaller signals of GABA, glutamine, glutathione, and taurine. These results indicate that 7-T MRS using a STEAM sequence with ultrashort TE and automated voxel repositioning provides excellent reproducibility of metabolites in the DLPFC.

Connectivity-Guided Theta Burst Transcranial Magnetic Stimulation Versus Repetitive Transcranial Magnetic Stimulation for Treatment-Resistant Moderate to Severe Depression: Magnetic Resonance Imaging Protocol and SARS-CoV-2-Induced Changes for a Randomized Double-blind Controlled Trial

Background

Depression is a substantial health and economic burden. In approximately one-third of patients, depression is resistant to first-line treatment; therefore, it is essential to find alternative treatments. Transcranial magnetic stimulation (TMS) is a neuromodulatory treatment involving the application of magnetic pulses to the brain that is approved in the United Kingdom and the United States in treatment-resistant depression. This trial aims to compare the clinical effectiveness, cost-effectiveness, and mechanism of action of standard treatment repetitive TMS (rTMS) targeted at the F3 electroencephalogram site with a newer treatment—a type of TMS called theta burst stimulation (TBS) targeted based on measures of functional brain connectivity. This protocol outlines brain imaging acquisition and analysis for the Brain Imaging Guided Transcranial Magnetic Stimulation in Depression (BRIGhTMIND) study trial that is used to create personalized TMS targets and answer the proposed mechanistic hypotheses.

Objective

The aims of the imaging arm of the BRIGhTMIND study are to identify functional and neurochemical brain signatures indexing the treatment mechanisms of rTMS and connectivity-guided intermittent theta burst TMS and to identify imaging-based markers predicting response to treatment.

Methods

The study is a randomized double-blind controlled trial with 1:1 allocation to either 20 sessions of TBS or standard rTMS. Multimodal magnetic resonance imaging (MRI) is acquired for each participant at baseline (before TMS treatment) with T1-weighted and task-free functional MRI during rest used to estimate TMS targets. For participants enrolled in the mechanistic substudy, additional diffusion-weighted sequences are acquired at baseline and at posttreatment follow-up 16 weeks after treatment randomization. Core data sets of T1-weighted and task-free functional MRI during rest are acquired for all participants and are used to estimate TMS targets. Additional sequences of arterial spin labeling, magnetic resonance spectroscopy, and diffusion-weighted images are acquired depending on the recruitment site for mechanistic evaluation. Standard rTMS treatment is targeted at the F3 electrode site over the left dorsolateral prefrontal cortex, whereas TBS treatment is guided using the coordinate of peak effective connectivity from the right anterior insula to the left dorsolateral prefrontal cortex. Both treatment targets benefit from the level of MRI guidance, but only TBS is provided with precision targeting based on functional brain connectivity.

Results

Recruitment began in January 2019 and is ongoing. Data collection is expected to continue until January 2023.

Conclusions

This trial will determine the impact of precision MRI guidance on rTMS treatment and assess the neural mechanisms underlying this treatment in treatment-resistant depressed patients.

Trial Registration

ISRCTN Registry ISRCTN19674644; https://www.isrctn.com/ISRCTN19674644

International Registered Report Identifier (IRRID)

DERR1-10.2196/31925

Comparison of different linear-combination modeling algorithms for short-TE proton spectra

Short-TE proton MRS is used to study metabolism in the human brain. Common analysis methods model the data as a linear combination of metabolite basis spectra. This large-scale multi-site study compares the levels of the four major metabolite complexes in short-TE spectra estimated by three linear-combination modeling (LCM) algorithms. 277 medial parietal lobe short-TE PRESS spectra (TE = 35 ms) from a recent 3 T multi-site study were preprocessed with the Osprey software. The resulting spectra were modeled with Osprey, Tarquin and LCModel, using the same three vendor-specific basis sets (GE, Philips and Siemens) for each algorithm. Levels of total N-acetylaspartate (tNAA), total choline (tCho), myo-inositol (mI) and glutamate + glutamine (Glx) were quantified with respect to total creatine (tCr). Group means and coefficient of variations of metabolite estimates agreed well for tNAA and tCho across vendors and algorithms, but substantially less so for Glx and mI, with mI systematically estimated as lower by Tarquin. The cohort mean coefficient of determination for all pairs of LCM algorithms across all datasets and metabolites was R 2 ¯ = 0.39, indicating generally only moderate agreement of individual metabolite estimates between algorithms. There was a significant correlation between local baseline amplitude and metabolite estimates (cohort mean R 2 ¯ = 0.10). While mean estimates of major metabolite complexes broadly agree between linear-combination modeling algorithms at group level, correlations between algorithms are only weak-to-moderate, despite standardized preprocessing, a large sample of young, healthy and cooperative subjects, and high spectral quality. These findings raise concerns about the comparability of MRS studies, which typically use one LCM software and much smaller sample sizes.

INSPECTOR: free software for magnetic resonance spectroscopy data inspection, processing, simulation and analysis

In vivo magnetic resonance spectroscopy (MRS) is a powerful tool for biomedical research and clinical diagnostics, allowing for non-invasive measurement and analysis of small molecules from living tissues. However, currently available MRS processing and analytical software tools are limited in their potential for in-depth quality management, access to details of the processing stream, and user friendliness. Moreover, available MRS software focuses on selected aspects of MRS such as simulation, signal processing or analysis, necessitating the use of multiple packages and interfacing among them for biomedical applications. The freeware INSPECTOR comprises enhanced MRS data processing, simulation and analytical capabilities in a one-stop-shop solution for a wide range of biomedical research and diagnostic applications. Extensive data handling, quality management and visualization options are built in, enabling the assessment of every step of the processing chain with maximum transparency. The parameters of the processing can be flexibly chosen and tailored for the specific research problem, and extended confidence information is provided with the analysis. The INSPECTOR software stands out in its user-friendly workflow and potential for automation. In addition to convenience, the functionalities of INSPECTOR ensure rigorous and consistent data processing throughout multi-experiment and multi-center studies.

Reduced striatal GABA in unmedicated children with ADHD at 7T

Attention-deficit hyperactive disorder (ADHD) is characterized by inattention and increased impulsive and hypermotoric behaviors.Despite the high prevalence and impact of ADHD, little is known about the underlying neurophysiology of ADHD. The main inhibitory and excitatory neurotransmitters γ-aminobutyric acid (GABA) and glutamate are receiving increased attention in ADHD and can be measured using Magnetic Resonance Spectroscopy (MRS). However, MRS studies in ADHD are limited. We measured GABA and glutamate in young unmedicated participants, utilizing high magnetic field strength. Fifty unmedicated children (26 with ADHD, 24 controls) aged 5-9 years completed MRS at 7T and behavioral testing. GABA and glutamate were measured in dorsolateral prefrontal cortex (DLPFC), anterior cingulate cortex (ACC), premotor cortex (PMC), and striatum, and estimated using LCModel. Children with ADHD showed poorer inhibitory control and significantly reduced GABA/Cr in the striatum, but not in ACC, DLPFC, or PMC regions. There were no significant group differences for Glu/Cr levels, or correlations with behavioral manifestations of ADHD. The primary finding of this study is a reduction of striatal GABA levels in unmedicated children with ADHD at 7T. These findings provide guidance for future studies or interventions. Reduced striatal GABA may be a marker for specific GABA-related treatment for ADHD.