Page Last Updated: October 18, 2025
Magnetic Resonance Spectroscopy (MRS)🔗
Magnetic Resonance Spectroscopy (MRS) provides measures of biochemicals involved in neuronal metabolism, neurotransmission, and oxidative stress. These measures enable researchers to interrogate biochemical mechanisms underlying the structural, functional, and behavioral trajectories. HBCD is the first study of this magnitude to include MRS in a comprehensive pediatric neuroimaging protocol through the development of Integrated Short-TE and Hadamard Multi-Sequence, or ISTHMUS (Hui et al., 2024).
Release Data🔗
MRS release data include both file-based (raw and processed data files in modality-specific formats) and tabulated (instrument and derived data in a standardized table format) data.
See the Data Structure Overview for a full explanation of these data types.
- Raw BIDS stored under subject- and session-specific
mrs/folders - Derivatives generated by the OSPREY-BIDS pipeline
- Tabulated data derived from Osprey pipeline outputs — see the full list of tables here
MRS files include metabolite (*_svs.nii.gz) and water reference (*_ref.nii.gz) data acquired via short-echo-time (TE = 35 ms; acq-shortTE) and HERCULES (spectral-edited, TE = 80 ms; acq-hercules). The JSON sidecar files include the dimensions of the NIfTI-MRS data array, holding different coil elements in dimension 5 and different transients in dimension 6. See BIDS Conversion Procedures.
hbcd/
|__ rawdata/
|__ sub-{ID}/
|__ ses-{V0X}/
|__ mrs/
|__ sub-{ID}_ses-{V0X}_acq-shortTE_run-{X}_svs.nii.gz (+JSON)
|__ sub-{ID}_ses-{V0X}_acq-shortTE_run-{X}_ref.nii.gz (+JSON)
|__ sub-{ID}_ses-{V0X}_acq-hercules_run-{X}_svs.nii.gz (+JSON)
|__ sub-{ID}_ses-{V0X}_acq-hercules_run-{X}_ref.nii.gz (+JSON)
Only the HERCULES/ file tree is displayed below. The unedited/ files generally follow similar naming conventions, with some exceptions (e.g., the BIDS field acq-shortTE is used instead of acq-hercules). See Osprey MRS Output User Guide on this page for usage guidance.
hbcd/
|__ derivatives/
|__ osprey/
|__ sub-{ID}/
|__ ses-{V0X}/
|__ HERCULES/
| |__ PreOspreyLocalizerReg/
| | |__ figures/
| | |__ <aal|c1|c2|c3>reference_seg_aligned_to_localizer.nii.gz
| | |__ reference_<img|seg>_aligned_to_localizer.nii.gz
| | |__ registration_summary.json
| | |__ transform_mat.npy
| | |__ readme.txt
| |
| |__ QuantifyResults/ (ALL +JSON)
| | |__ <diff1|diff2|sum>_AlphaCorrWaterScaledGroupNormed_Voxel_1_Basis_1.tsv
| | |__ <diff1|diff2|sum>_AlphaCorrWaterScaled_Voxel_1_Basis_1.tsv
| | |__ <diff1|diff2|sum>_CSFWaterScaled_Voxel_1_Basis_1.tsv
| | |__ <diff1|diff2|sum>_TissCorrWaterScaled_Voxel_1_Basis_1.tsv
| | |__ <diff1|diff2|sum>_amplMets_Voxel_1_Basis_1.tsv
| | |__ <diff1|diff2|sum>_rawWaterScaled_Voxel_1_Basis_1.tsv
| | |__ <diff1|diff2|sum>_tCr_Voxel_1_Basis_1.tsv
| |
| |__ Reports/
| | |__ reportFigures/
| | |__ sub-{ID}-report.html
| |
| |__ SegMaps/
| | |__ TissueFractions_Voxel_1.tsv (+JSON)
| | |__ sub-{ID}_ses-{V0X}_acq-hercules_svs.nii_space-scanner_Voxel-1_label-CSF.nii.gz
| | |__ sub-{ID}_ses-{V0X}_acq-hercules_svs.nii_space-scanner_Voxel-1_label-GM.nii.gz
| | |__ sub-{ID}_ses-{V0X}_acq-hercules_svs.nii_space-scanner_Voxel-1_label-Tha.nii.gz
| | |__ sub-{ID}_ses-{V0X}_acq-hercules_svs.nii_space-scanner_Voxel-1_label-WM.nii.gz
| |
| |__ VoxelMasks/
| | |__ sub-{ID}_ses-{V0X}_acq-hercules_svs_space-scanner_mask.nii.gz
| |
| |__ LogFile.txt
| |__ QM_processed_spectra.tsv (+JSON)
| |__ SummaryMRSinMRS.md
| |__ subject_names_and_excluded.tsv (+JSON)
| |__ wrapper_settings.mat (+JSON)
|
|__ unedited/
Data Acquisition🔗
The MRS acquisition protocol was optimized to maximize signal-to-noise across multiple low-concentration metabolites while maintaining an acquisition time (TA) under 9 minutes. The MRS acquisition centers on a single voxel Point-RESolved Spectroscopy (PRESS) (Bottomley, 1987) localization (30×23×23 mm^3) in the bilateral thalamus, with SVS localizer acquisitions to define the ROI. The ISTHMUS sequence incorporates a short TE (35 ms) unedited block at the beginning of the sequence for optimal measurement of high concentration metabolites, including N-acetylasparte, followed by an advanced Hadamard encoded spectral editing scheme (Oeltzschner et al., 2019) to derive reliable and reproducible measures of the low-concentration metabolites.
The primary metabolites measured for HBCD include:
- N-acetylaspartate (NAA): a marker of neuronal mitochondrial metabolism
- Glutamate and γ-aminobutyric acid (GABA): the principal excitatory and inhibitory neurotransmitters
- Glutathione (GSH): the most abundant antioxidant involved in protection against reactive oxygen species in the human brain
Additional metabolites measured include NAA, lactate, ascorbate, creatine, myo-inositol, glutamine, and total choline (Oeltzschner et al., 2019).
One limitation to the incorporation of MRS into human connectome studies is proper control for scanner drift. Because MRS relies on the frequency of the measured signals, uncorrected frequency drift during data acquisition is very detrimental to data quality, as it changes the contribution of coedited signals as well as editing efficiency (Harris et al., 2014). To mitigate drift, an innovative approach was taken to incorporate interleaved water referencing (Edden et al., 2016) for real-time frequency correction, implemented on the Philips platform at the outset of HBCD, and in Y3 for Siemens and GE.
Data Processing🔗
HBCD MRS data are processed with a customized automated pipeline based on OSPREY (Oeltzschner et al., 2020; Zöllner et al., 2023) - see the OSPREY-BIDS documentation. Derivatives include quality control metrics and metabolite estimates (ratios to creatine and water, with and without tissue correction). See the list of OSPREY derivatives included in the release above.
Osprey MRS Output User Guide🔗
The primary outcome variables from MRS data processed through the Osprey pipeline are metabolite concentrations. The ISTHMUS acquisition used in HBCD generates four spectra, each modeled separately using linear combination modeling with inclusive basis sets (since most metabolites may contribute, at least minimally, to all spectra).
Each spectrum provides optimal quantification for a different subset of metabolites. The following recommended files are all available as tabulated derivatives that summarize participant data across Osprey derivatives.
| Spectrum | Best-Quantified Metabolites | Table Name |
|---|---|---|
| Short-TE Unedited | tNAA, tCr, tCho, mI, Glx, Scyllo | img_osprey_unedited_A_TissCorrWaterScaled_Voxel_1_Basis_1 |
| HERCULES Sum | NAA, Glu, Gln | img_osprey_HERCULES_sum_TissCorrWaterScaled_Voxel_1_Basis_1 |
| HERCULES Diff1 | GABA+ | img_osprey_HERCULES_diff1_TissCorrWaterScaled_Voxel_1_Basis_1 |
| HERCULES Diff2 | GSH, Lac, NAAG, PE | img_osprey_HERCULES_diff2_TissCorrWaterScaled_Voxel_1_Basis_1 |
Quantification Approaches🔗
While MRS is a quantitative approach, the values are relative, relying on an internal reference signal for quantification. Because the optimal reference for HBCD data has not yet been established, several quantification methods are provided (see Harris et al., 2015):
| Method | Recommended Use |
|---|---|
| Unscaled model amplitudes | Do not use for analyses |
| Metabolite-to-creatine ratios | Normalized to total creatine (tCr), which is only found in tissue, not CSF (so not suitable for all cases). tCr is relatively stable across the lifespan |
| Raw water-scaled concentrations | |
| Tissue-Corrected Water-Scaled Concentrations | Recommended default (TissCorrWaterScaled TSV files) |
| Alpha-corrected variants | Use only for GABA+ |
References🔗
Bottomley, P. A. (1987). Spatial localization in NMR spectroscopy in vivo. Annals of the New York Academy of Sciences, 508(1), 333–348. https://doi.org/10.1111/j.1749-6632.1987.tb32915.x
Edden, R. A. E., Oeltzschner, G., Harris, A. D., Puts, N. A. J., Chan, K. L., Boer, V. O., Schär, M., & Barker, P. B. (2016). Prospective frequency correction for macromolecule-suppressed GABA editing at 3T. Journal of Magnetic Resonance Imaging, 44(6), 1474–1482.. https://doi.org/10.1002/jmri.25304
Harris, A. D., Glaubitz, B., Near, J., John Evans, C., Puts, N. A. J., Schmidt-Wilcke, T., Tegenthoff, M., Barker, P. B., & Edden, R. A. E. (2014). Impact of frequency drift on gamma-aminobutyric acid-edited MR spectroscopy. Magnetic Resonance in Medicine, 72(4), 941–948. https://doi.org/10.1002/mrm.25009
Harris, A. D., Puts, N. A. J., & Edden, R. A. E. (2015). Tissue correction for GABA-edited MRS: Considerations of voxel composition, tissue segmentation, and tissue relaxations.Journal of Magnetic Resonance Imaging, 42(5), 1431–1440. https://doi.org/10.1002/jmri.24903
Hui, S. C. N., Murali-Manohar, S., Zöllner, H. J., Hupfeld, K. E., Davies-Jenkins, C. W., Gudmundson, A. T., Song, Y., Yedavalli, V., Wisnowski, J. L., Gagoski, B., Oeltzschner, G., & Edden, R. A. E. (2024). Integrated Short-TE and Hadamard-edited Multi-Sequence (ISTHMUS) for advanced MRS. Journal of Neuroscience Methods, 409(110206), 110206. https://doi.org/10.1016/j.jneumeth.2024.110206
Oeltzschner, G., Saleh, M. G., Rimbault, D., Mikkelsen, M., Chan, K. L., Puts, N. A. J., & Edden, R. A. E. (2019). Advanced Hadamard-encoded editing of seven low-concentration brain metabolites: Principles of HERCULES. NeuroImage, 185, 181–190. https://doi.org/10.1016/j.neuroimage.2018.10.002
Oeltzschner, G., Zöllner, H. J., Hui, S. C. N., Mikkelsen, M., Saleh, M. G., Tapper, S., & Edden, R. A. E. (2020). Osprey: Open-source processing, reconstruction & estimation of magnetic resonance spectroscopy data. Journal of Neuroscience Methods, 343(108827), 108827. https://doi.org/10.1016/j.jneumeth.2020.108827
Träber, F., Block, W., Lamerichs, R., Gieseke, J., & Schild, H. H. (2004). 1H metabolite relaxation times at 3.0 tesla: Measurements of T1 and T2 values in normal brain and determination of regional differences in transverse relaxation. Journal of Magnetic Resonance Imaging, 19(5), 537–545. https://doi.org/10.1002/jmri.20053
Zöllner, H. J., Davies-Jenkins, C. W., Lee, E. G., Hendrickson, T. J., Clarke, W. T., Edden, R. A. E., Wisnowski, J. L., Gudmundson, A. T., & Oeltzschner, G. (2023). Continuous automated analysis workflow for MRS studies. Journal of Medical Systems, 47(1), 69. https://doi.org/10.1007/s10916-023-01969-6