Osprey: Open-source processing, reconstruction & estimation of magnetic resonance spectroscopy data

Background Processing and quantitative analysis of magnetic resonance spectroscopy (MRS) data are far from standardized and require interfacing with third-party software. Here, we present Osprey, a fully integrated open-source data analysis pipeline for MRS data, with seamless integration of pre-processing, linear-combination modelling, quantification, and data visualization. New Method Osprey loads multiple common MRS data formats, performs phased-array coil combination, frequency-and phase-correction of individual transients, signal averaging and Fourier transformation. Linear combination modelling of the processed spectrum is carried out using simulated basis sets and a spline baseline. The MRS voxel is coregistered to an anatomical image, which is segmented for tissue correction and quantification is performed based upon modelling parameters and tissue segmentation. The results of each analysis step are visualized in the Osprey GUI. The analysis pipeline is demonstrated in 12 PRESS, 11 MEGA-PRESS, and 8 HERMES datasets acquired in healthy subjects. Results Osprey successfully loads, processes, models, and quantifies MRS data acquired with a variety of conventional and spectral editing techniques. Comparison with Existing Method(s) Osprey is the first MRS software to combine uniform pre-processing, linear-combination modelling, tissue correction and quantification into a coherent ecosystem. Compared to existing compiled, often closed-source modelling software, Osprey’s open-source code philosophy allows researchers to integrate state-of-the-art data processing and modelling routines, and potentially converge towards standardization of analysis. Conclusions Osprey combines robust, peer-reviewed data processing methods into a modular workflow that is easily augmented by community developers, allowing the rapid implementation of new methods.

[1]  Jon-Fredrik Nielsen,et al.  Across‐vendor standardization of semi‐LASER for single‐voxel MRS at 3T , 2019, NMR in biomedicine.

[2]  Kenneth Levenberg A METHOD FOR THE SOLUTION OF CERTAIN NON – LINEAR PROBLEMS IN LEAST SQUARES , 1944 .

[3]  A. Henning,et al.  In vivo estimation of transverse relaxation time constant (T2) of 17 human brain metabolites at 3T , 2018, Magnetic resonance in medicine.

[4]  D. Marquardt An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[5]  R. Edden,et al.  Gannet: A batch‐processing tool for the quantitative analysis of gamma‐aminobutyric acid–edited MR spectroscopy spectra , 2014, Journal of magnetic resonance imaging : JMRI.

[6]  M. Garwood,et al.  Simultaneous in vivo spectral editing and water suppression , 1998, NMR in biomedicine.

[7]  R. Mattson,et al.  Localized 1H NMR measurements of gamma-aminobutyric acid in human brain in vivo. , 1993, Proceedings of the National Academy of Sciences of the United States of America.

[8]  E Moser,et al.  Proton T 1 and T 2 relaxation times of human brain metabolites at 3 Tesla , 2001, NMR in biomedicine.

[9]  S. Holland,et al.  NMR relaxation times in the human brain at 3.0 tesla , 1999, Journal of magnetic resonance imaging : JMRI.

[10]  Arvind Caprihan,et al.  Use of tissue water as a concentration reference for proton spectroscopic imaging , 2006, Magnetic resonance in medicine.

[11]  U. Klose In vivo proton spectroscopy in presence of eddy currents , 1990, Magnetic resonance in medicine.

[12]  Darren Price,et al.  Methodology for improved detection of low concentration metabolites in MRS: Optimised combination of signals from multi-element coil arrays , 2014, NeuroImage.

[13]  Matthew D Robson,et al.  OXSA: An open-source magnetic resonance spectroscopy analysis toolbox in MATLAB , 2017, PloS one.

[14]  S. Provencher Estimation of metabolite concentrations from localized in vivo proton NMR spectra , 1993, Magnetic resonance in medicine.

[15]  Karl J. Friston,et al.  CHAPTER 2 – Statistical parametric mapping , 2007 .

[16]  Richard A. E. Edden,et al.  Simultaneous edited MRS of GABA and glutathione , 2016, NeuroImage.

[17]  Peter B Barker,et al.  Measuring T2 in vivo with J‐difference editing: Application to GABA at 3 tesla , 2012, Journal of magnetic resonance imaging : JMRI.

[18]  Jun Shen,et al.  Fast computation of full density matrix of multispin systems for spatially localized in vivo magnetic resonance spectroscopy , 2017, Medical physics.

[19]  Muhammad G Saleh,et al.  Edited 1H magnetic resonance spectroscopy in vivo: Methods and metabolites , 2017, Magnetic resonance in medicine.

[20]  Michael-Paul Schallmo,et al.  Big GABA: Edited MR spectroscopy at 24 research sites , 2017, NeuroImage.

[21]  Li Pan,et al.  Multi-vendor standardized sequence for edited magnetic resonance spectroscopy , 2019, NeuroImage.

[22]  Kirstie J. Whitaker,et al.  Raincloud plots: a multi-platform tool for robust data visualization , 2018, PeerJ Prepr..

[23]  Theodoros N. Arvanitis,et al.  A constrained least‐squares approach to the automated quantitation of in vivo 1H magnetic resonance spectroscopy data , 2011, Magnetic resonance in medicine.

[24]  Jorge Nocedal,et al.  Algorithm 778: L-BFGS-B: Fortran subroutines for large-scale bound-constrained optimization , 1997, TOMS.

[25]  Karl Landheer,et al.  Theoretical description of modern 1H in Vivo magnetic resonance spectroscopic pulse sequences , 2020, Journal of magnetic resonance imaging : JMRI.

[26]  Blaise deB. Frederick,et al.  Grey and white matter GABA level differences in the human brain using two‐dimensional, J‐resolved spectroscopic imaging , 2005, NMR in biomedicine.

[27]  D. van Ormondt,et al.  Improved algorithm for noniterative time-domain model fitting to exponentially damped magnetic resonance signals , 1987 .

[28]  Sabine Van Huffel,et al.  An automated quantitation of short echo time MRS spectra in an open source software environment: AQSES , 2007, NMR in biomedicine.

[29]  Karl J. Friston,et al.  Statistical parametric mapping , 2013 .

[30]  R. Edden,et al.  Measuring the longitudinal relaxation time of GABA in vivo at 3 tesla , 2013, Journal of magnetic resonance imaging : JMRI.

[31]  S. A. COONS Constrained least-squares , 1978, Comput. Graph..

[32]  M. Terpstra,et al.  MEGA‐PRESS of GABA+: Influences of acquisition parameters , 2019, NMR in biomedicine.

[33]  Ashley D Harris,et al.  Tissue correction for GABA‐edited MRS: Considerations of voxel composition, tissue segmentation, and tissue relaxations , 2015, Journal of magnetic resonance imaging : JMRI.

[34]  Zenon Starčuk,et al.  Quantitation of magnetic resonance spectroscopy signals: the jMRUI software package , 2009 .

[35]  Vanhamme,et al.  Improved method for accurate and efficient quantification of MRS data with use of prior knowledge , 1997, Journal of magnetic resonance.

[36]  Michael Schär,et al.  HERMES: Hadamard encoding and reconstruction of MEGA‐edited spectroscopy , 2016, Magnetic resonance in medicine.

[37]  Julien Cohen-Adad,et al.  SCT: Spinal Cord Toolbox, an open-source software for processing spinal cord MRI data , 2017, NeuroImage.

[38]  Richard A. E. Edden,et al.  Advanced Hadamard-encoded editing of seven low-concentration brain metabolites: Principles of HERCULES , 2019, NeuroImage.

[39]  Richard A. E. Edden,et al.  Big GABA II: Water-referenced edited MR spectroscopy at 25 research sites , 2019, NeuroImage.

[40]  M. Terpstra,et al.  Influence of fitting approaches in LCModel on MRS quantification focusing on age‐specific macromolecules and the spline baseline , 2019, NMR in biomedicine.

[41]  Peter Jezzard,et al.  Frequency and phase drift correction of magnetic resonance spectroscopy data by spectral registration in the time domain , 2015, Magnetic resonance in medicine.

[42]  A. Haley,et al.  Molality as a unit of measure for expressing 1H MRS brain metabolite concentrations in vivo. , 2003, Magnetic resonance imaging.

[43]  Anke Henning,et al.  Methodological consensus on clinical proton MRS of the brain: Review and recommendations , 2019, Magnetic resonance in medicine.

[44]  Peter Jezzard,et al.  Advanced processing and simulation of MRS data using the FID appliance (FID‐A)—An open source, MATLAB‐based toolkit , 2017, Magnetic resonance in medicine.

[45]  Davide Poggiali,et al.  Raincloud plots: a multi-platform tool for robust data visualization , 2019, Wellcome open research.

[46]  Jason C. Crane,et al.  SIVIC: Open-Source, Standards-Based Software for DICOM MR Spectroscopy Workflows , 2013, Int. J. Biomed. Imaging.

[47]  Ralph Noeske,et al.  Advanced single voxel 1H magnetic resonance spectroscopy techniques in humans: Experts' consensus recommendations , 2020, NMR in biomedicine.

[48]  Michael Elad,et al.  ESPIRiT—an eigenvalue approach to autocalibrating parallel MRI: Where SENSE meets GRAPPA , 2014, Magnetic resonance in medicine.

[49]  Dennis W J Klomp,et al.  1H–MRS processing parameters affect metabolite quantification: The urgent need for uniform and transparent standardization , 2017, NMR in biomedicine.

[50]  Jorge Nocedal,et al.  A Limited Memory Algorithm for Bound Constrained Optimization , 1995, SIAM J. Sci. Comput..