Skip to main content

Reproducibility of brain metabolite concentration measurements in lesion free white matter at 1.5 T



Post processing for brain spectra has a great influence on the fit quality of individual spectra, as well as on the reproducibility of results from comparable spectra. This investigation used pairs of spectra, identical in system parameters, position and time assumed to differ only in noise. The metabolite amplitudes of fitted time domain spectroscopic data were tested on reproducibility for the main brain metabolites.


Proton spectra of white matter brain tissue were acquired with a short spin echo time of 30 ms and a moderate repetition time of 1500 ms at 1.5 T. The pairs were investigated with one time domain post-processing algorithm using different parameters. The number of metabolites, the use of prior knowledge, base line parameters and common or individual damping were varied to evaluate the best reproducibility.


The protocols with most reproducible amplitudes for N-acetylaspartate, creatine, choline, myo-inositol and the combined Glx line of glutamate and glutamine in lesion free white matter have the following common features: common damping of the main metabolites, a baseline using only the points of the first 10 ms, no additional lipid/macromolecule lines and Glx is taken as the sum of separately fitted glutamate and glutamine. This parameter set is different to the one delivering the best individual fit results.


All spectra were acquired in “lesion free” (no lesion signs found in MR imaging) white matter. Spectra of brain lesions, for example tumors, can be drastically different. Thus the results are limited to lesion free brain tissue. Nevertheless the application to studies is broad, because small alterations in brain biochemistry of lesion free areas had been detected nearby tumors, in patients with multiple sclerosis, drug abuse or psychiatric disorders.


Main metabolite amplitudes inside healthy brain can be quantified with a normalized root mean square deviation around 5 % using CH3 of creatine as reference. Only the reproducibility of myo-inositol is roughly twice as bad. The reproducibility should be similar using other references like internal or external water for an absolute concentration evaluation and are not influenced by relaxation corrections with literature values.

Peer Review reports


Nuclear magnetic resonance spectroscopy has developed over the last decades as one of the exciting fields in medical physics for the non-invasive measurement of molecular concentrations inside the body [17]. In principle the method is widely available at magnetic resonance imaging (MRI) units with field strengths above 1 T but its application is restricted for two main reasons. Although data acquisition can be performed automatically for optimal results an experienced operator is necessary. Additionally no standard post processing of the acquired data exists [8]. This investigation tries to enhance the post processing towards reproducible results by varying the post processing parameters for identically acquired spectra pairs.

Biochemical concentrations inside specific body regions can be evaluated by fitting magnetic resonance spectroscopic data using model signals of the metabolites to be quantified. The fit algorithm detects metabolites with a fair signal to noise ratio and sufficient long spin-spin relaxation time T2 in relation to the echo time TE. Metabolites with short T2’s (broad lines) are hardly separable and combined as baseline.

Metabolite model signals in time domain have the following parameters adjustable by the fit algorithm: amplitude (corresponding to the area under a line in frequency domain data), frequency (ppm position), damping (line width) and phase (influencing the line shape in frequency domain) [911]. The adoption of “prior knowledge” can reduce the number of adjustable parameters and/or restrict them to a certain region. For example: the metabolite frequency is fixed to another metabolite or limited to a region around the theoretical value; the metabolites completely or partly have a common damping or phase.

Fit procedures for spectral data adjust the parameters for the specified metabolites to minimize χ2 respecting all prior knowledge restrictions. Usually, χ2 in time domain is defined as the sum Σii|2, where Δi is the difference between measured value and fit value for the data point sampled at time ti. If the assumed errors σi of the measured values are not identical for all i, χ2 is taken as Σiii|2 to weight data points according to the error σi in measurement i. For the ideal case, fitting of all metabolites with appropriate parameters should result in a difference between fit and the acquired signal around zero containing the pure noise signal. This ideal case is not achievable, because signals can be masked by noise, they can overlay each other or signals with low T2 (broad lines) are not completely represented by the base line. But it is evident that a flat difference line with a low χ2 can be easily achieved using a large parameter set (i. e. using many metabolites with all four parameters of each metabolite freely adjustable).

Until today no commonly accepted post processing exists to quantify concentrations. Firstly, frequency [12] or time domain models [10, 11] can be used for the appropriate data. Secondly, a decision is necessary on the metabolites taken into account and on the restrictions of each metabolite’s post processing parameters by the use of prior knowledge.

For a series of acquisitions the fit quality can be rated by the mean of all χ2. Using different metabolite sets with varying post-processing parameters can detect a minimal mean χ2. A large set of metabolites with unrestricted parameters decreases each χ2 and hence also the mean of all χ2’s. But such a parameter setting not necessarily is the one evaluating metabolite concentrations with the best reproducibility. A post processing with a good reproducibility should compute nearly identical metabolite concentrations for data known to have the same metabolite basis.

As approach for the preparation of in vivo data, differing only in noise, ideally a complete acquisition with 2 N additions is divided in two acquisitions with N additions. The two spectra for comparison build a pair consisting of the sum of all even and the sum of all odd acquired signals. The spectra of this procedure are as reproducible as possible for in vivo acquisitions. The MR-system software did not allow adding even and odd excitations to summed signals. Instead the two spectra of a pair were acquired directly consecutively with the same system parameters at an identical position and time point (identical in time with respect to physiological concentration changes). As long as the investigated subject did not move within the two acquisitions this should match the best achievable reproducibility for in vivo spectroscopy. For the comparison it is assumed to have a pair of two spectra formed by an identical metabolite basis overlaid with different noise. The post processing should evaluate similar results for the main (clearly visible) metabolites.

This investigation concentrates on one post processing software package metabolite report [9, 13, 14], mostly identical to the new product software “syngo via” from Siemens Healthcare (Erlangen, Germany) neglecting all aspects and discussions on acquisition schemes or comparisons to other post processing algorithms.

The purpose of this paper is to achieve a good reproducibility of the major metabolite amplitudes from the acquired pairs of in vivo spectra by adjusting the post processing fit parameters. These amplitudes are the basis for all relative or absolute concentration calculations.



This study was retrospective with the data of a master thesis in Biology at the Ruhr University of Bochum (RUB). All participants had given written informed consent. The local ethics committee of the RUB approved the study and all data evaluations.

Fifty-four spectra pairs were acquired, 28 pairs from six females and 26 pairs from five males (Table 1). None of the healthy subjects showed any evidence of a brain lesion according to the acquired T2-weighted sequences. All subjects were of similar age. Gender and age were similarly distributed on the number of subjects as well as on the number of spectra acquisitions (Table 1).

Table 1 Subject demographics

Spectra acquisition

Subjects were positioned head first, supine with their head slightly fixed inside the coil. All measurements were carried out on a Magnetom Espree (horizontal 1.5 T field, wide bore, Siemens Healthcare, Erlangen, Germany). The signal was received with a 12-channel head matrix coil and the circularly polarized body coil was used for excitation. Firstly, for exclusion of any brain lesion as well as for spectroscopy voxel placement, T2 weighted axial “Turbo Spin Echo” images of the entire brain were acquired followed by two (coronal and sagittal) “True Fisp” (true fast imaging with steady state precession) sequences or a sagittal “True Fisp” and a coronal “Tirm” (turbo inversion recovery magnitude) sequence for voxel placement (sequence parameters in Table 2).

Table 2 Imaging sequence parameters

Two pairs of “point resolved single voxel spectra” (PRESS) were acquired with identical system parameters and localized at an identical position on each brain hemisphere for each subject. Always cubic voxels of 2 cm × 2 cm × 2 cm were selected and placed on the left and right brain hemisphere central in anterior posterior direction above the ventricle (marked in Fig. 1 as VOI). Position and orientation were adjusted to contain white matter as much as possible. Voxels were placed within ±5 cm relative to the central z-axis position (center of the magnet in head feet direction). Six saturation pulses were positioned along each plane of the cube (Fig. 1) for outer volume suppression (marked in Fig. 1 as numbered bars surrounding the VOI). The most effective one was oriented towards the ventricle (bar 1 in Fig. 1) and next effective one towards the midline of the cerebrum (bar 2 in Fig. 1). Prior to paired spectra acquisitions on each side, system parameters were automatically adjusted and the magnetic field inside the voxel was homogenized by an automatic shim. Additionally an experienced operator (MB or SM) optimized shim and water suppression.

Fig. 1

Example of voxel position (VOI) in axial (a), coronal (b) and sagittal (c) images, as well as the saturation regions (numbered from 1 to 6)

One thousand twenty four data points were sampled in 512 ms in a two-fold oversampling technique. In the thesis of MS the short sampling window of 512 ms enables faster repetition times TR for measuring spin–lattice relaxation times T1. The echo time TE was short (30 ms) and the repetition time TR moderate (1500 ms). The two spectra of a pair were acquired directly consecutive using 128 additions for each spectrum. For signal evaluation the signal of all 12 channels of the head coil was used. Total acquisition time on one side was 6:36 min (2 × 3:18 min) preceded by the system parameter adjustment time of 7 to 10 min. Total MR time for each subject was about 30 to 40 min for one session.


The post-processing software “metabolite report” was a work in progress package from Siemens Healthcare. It is the precursor of the currently available commercial spectroscopy software within the new product “syngo via” and nearly identical to this package (personal communication with Elisabeth Weiland [15], see acknowledgement). “Automated post processing of this software package consists of the following steps: 1] identification of the prominent metabolites by cross-correlation to a database (i.e., N-acetylaspartate (Naa), choline (Cho), creatine (Cre), myoinositol (Myo) for short echo brain data), 2] determination of the B 0 shift and starting values for the fit parameters, 3] residual water removal, and 4] fit of the spectroscopic data. The fit is based on a basis set of metabolic model signals, which are simulated with the use of literature values for chemical shifts and coupling constants. In addition, truncating and remodeling of the first data points handle baseline artifacts” [9, 15].

For the analysis of the short TE data of this study, the following metabolites were included in the fit (major metabolites are marked bold):

  • Creatine with separately fitted CH2 and CH3 groups (Cre2 and Cre3)

  • N-Acetylaspartat (Naa)

  • Choline (Cho)

  • Myo-inositol (Myo)

  • Glx as a combined Glu/Gln model or separately fitted Glutamate (Glu) and Glutamine (Gln) with the sum of both amplitudes as Glx value (in case of a not found Gln, Glu alone was used as Glx value)

  • Two lipid/macromolecule lines around 0.9 and 1.3 ppm in some cases (Lipid_1, Lipid_2)

  • GABA (Gamma Amino Butter Acid) in a few cases

Metabolites were modeled as Gauss signals. Only the lipids were parameterized as Lorentz singlets with individual damping (line width) as recommended for short echo time acquisitions at 1.5 T [16, 17]. All model signals were filtered by a pass band from 0 to 4.2 ppm to exclude frequencies outside the area of interest (see model files in Fig. 2). This pass band filter influences the amplitude for relative or absolute concentration calculations by a constant factor respected within the software metabolite record. For pure reproducibility investigations a constant factor is not relevant.

Fig. 2

Model metabolites in frequency domain for a 30 ms PRESS sequence calculated with NMR-Scope of jMRUI using a damping of 2 Hz [10, 11]. (Chemical structures build with free Chemscetch from ACD Labs, Toronto, Canada)

For each metabolite the amplitude, the frequency and the damping (line width) are parameters obtained from the fit procedure. It is possible to use an individual phase of each model signal as fit parameter, but as commonly accepted a joint phase is used for all metabolites. The fitted metabolites had to be above the noise level in frequency domain by default factor of eight. Otherwise the metabolite is omitted and the fit is restarted without this metabolite (personal communication with EW). In this study the default factor was never changed. The software marks metabolites chosen for the fit but omitted, because it is too small in signal to noise. For reproducibility evaluation a decision is necessary how to respect pairs with a one signal above and the other below the noise level. In our investigation the not found signal was set to zero. In case of both signals below the noise level this pair was not respected.

A biochemical concentration is calculated with the fitted amplitude of a metabolite. Therefore, the post processing parameters should be optimized to deliver reliable and most reproducible amplitudes (proportional to metabolite concentration) for acquisitions assumed to differ only in noise. The frequency of a metabolite line was adjustable by the fit algorithm within boundaries to the theoretical value or it could be referenced to a fixed distance to another metabolite, which frequency can be more easily detected. If the frequency was not fixed to that of another metabolite the bounds were set to the theoretical value ± 0.1 ppm for all metabolites except the lipids (±0.2 ppm). The damping could be fitted individually for each metabolite or one or more groups of metabolites could have a common damping. The damping for single metabolites or a group of metabolites was restricted to values between 1 and 10 Hz except for lipids (1 to 25 Hz, see below).

The separate fit of the CH3 and CH2 creatine lines was necessary, because the residual water suppression diminishes the amplitude of the CH2 group disturbing the theoretical ratio of 3/2. The water removal consists of an iterative low pass filter around the frequency at 4.7 ppm with a fixed number of 22 iterations. At 1.5 T an accurate separation of Glu and Gln is difficult. Therefore Glx as a combined model with a fixed concentration relation between both or the sum of separately fitted Glu and Gln was evaluated as one of the main metabolites. The baseline was varied using 20 or 50 data points corresponding to 10 or 25 ms, respectively. The two “dummy lipids” were introduced as additional lines to model the base line solely below 1.8 ppm. These lipid lines incorporate macromolecules.

For detecting the most reproducible results using pairs of consecutive acquisitions with identical system parameters, no relaxation corrections had to be applied to the fitted amplitudes. It is sufficient to examine the unchanged amplitudes from which relative or absolute concentrations can be calculated. The error most likely increases (reproducibility decreases) with corrections and/or evaluation of absolute concentrations [mmol/l] or relative concentrations referenced to another metabolite. Both relaxation corrections as well as evaluation of concentrations involve measured or assumed values increasing the concentration errors.

Reproducibility calculations

Reproducibility is calculated with two directly consecutive acquired spectra (a pair) on one side of the brain. They are identical concerning the biochemical metabolite basis including the ones with a short T2 simulated as base line and they coincide in system parameter settings. For concentration independent quantification of the reproducibility the normalized squared difference of a pair is calculated. The possible range of the normalized difference is between −2 and 2. Case of “±2” consists of a zero and a non-zero value, which means the non-zero value is 200 % of the mean value. Typical values are much smaller than 1 for each main metabolite. The root of the mean value of all acquired pairs defines the normalized root mean square deviation, which simply is designated as reproducibility number rn (Eq. 1).

$$ rn(met)=\sqrt{\frac{1}{N}{\displaystyle {\sum}_{i=1}^N{\left(\frac{\left( am{p}_{i_2}(met)- am{p}_{i_1}(met)\right)}{\left( am{p}_{i_2}(met)+ am{p}_{i_1}(met)/2\right)}\right)}^2}} $$


reproducibility number; normalized root mean square deviation for metabolite “met”


number of acquisition pairs


amplitude of metabolite met for acquisition pair with index i. The index within a measurement pair is j = 1, 2.

The “reproducibility number” rn(met) describes the mean relative (percentage) deviation between N measurement pairs for metabolite met. So a low rn(met) value means a good reproducibility of metabolite met.

For the calculation of rn no reference is necessary to exclude different system settings, because the compared values of a pair were acquired with identical parameters. Normalization to the mean value was used to ease the reproducibility comparison between different metabolites.

Firstly the reproducibility is tested with only one spectrum pair of the eleven subjects. The first acquisition pair on the left-brain hemisphere is evaluated using the parameters specified in table 3 for 25 different protocols. As measure of reproducibility for all metabolites the added rn percentage value of Naa, Cre3, Cho, Myo and Glx was used (Table 3). Other sums of rn with less metabolites were also calculated for overall reproducibility rating. Post processings with a low sum (marked underlined and italic Table 3) were evaluated for all 54 pairs (Table 4).

Table 3 Description of tested protocols
Table 4 Reproducibility results using 11 left brain pairs of spectra from eleven subjects

Additionally the metabolite amplitudes for the protocols with the best reproducibility were referenced to the amplitude of Cre3 (ampref(met) = amp(met)/amp(Cre3)) to check the influence on this step to the reproducibility.

In principal an unlimited number of different post-processing protocols are possible. Hence in this investigation the protocols were varied to answer the following questions.

Is the reproducibility increased or decreased by the use of:

  1. 1.

    additional low signal to noise metabolites beside the major metabolites?

    To answer this question as example GABA was used with peaks in between 2 and 3 ppm (Fig. 3a). N-Acetyl Aspartyl Glutamate (NaaG) was also investigated, but in this case the Naa evaluation was distorted too much even using the sum of Naa and NaaG (Fig. 3c) and therefor was not further respected.

    Fig. 3

    Fit (red) of the same data (white) calculated with protocol 5 resulting in a “flat” (a) and with protocol 18 with a “less flat” difference line (green) (b). The baseline is orange. χ2 calculated by the post processing software was 3348 (a) and 5122 (b), respectively and is a little less than the mean values of Table 3. Image (c) shows the distorted Naa evaluation with NaaG as additional metabolite. A:, P: and W: are the amplitude of a metabolite, its frequency in ppm and its line width in Hz. Note: the lipid amplitudes inside healthy brain tissue are very small compared to that of tumors. In most cases, depending on the protocol, no or only one lipid line was detected. This is also documented by both protocols with the best reproducibility (protocol 14 and 18), which did not use lipid models

  2. 2.

    separately fitted glutamate and glutamine values instead of a predefined Glx model?

  3. 3.

    a flat or a more curved baseline using the first 10 or 25 ms for the base line calculation?

  4. 4.

    an independent damping for each metabolite or a common damping for the major metabolites?

  5. 5.

    common damping for major metabolites except for Myo, because its frequency is near the water frequency?

  6. 6.

    two additional lipid/macromolecule signals below 1.8 ppm as additional support for the baseline?

  7. 7.

    the frequencies of glutamate and glutamine as independent fit parameters or the use of fixed frequencies in relation to the most prominent metabolite line of Naa?

  8. 8.

    relative amplitudes (referencing the metabolites to Cre3)?


The results using 11 pairs of acquisitions are summarized in Table 3. The protocol with the lowest χ2 is, as expected, No 5, the one using the most metabolites with freely adjustable parameters. In contrast to the lowest χ2 protocol 5 is not in the top category concerning the reproducibility.

The flow chart of Fig. 4 shows the strategy for the evaluation of the best protocols using eleven pairs, one from each participant. For each question protocols with a significantly lower reproducibility performance were ruled out from further questions. Compared protocols with a difference in the sum of all rn numbers lower than 3 % were defined as similar and both were used for further investigations. Protocols 6, 9a, 13a, 14, 18, 19a are left after this decision process. These protocols except 9a were investigated with all 54 pairs for the detection of the best reproducibility. 9a was replaced for total evaluation by the more reproducible protocols 10 and 12. These seven protocols have the best reproducibility sum (< 31 %) of all main metabolite.

Fig. 4

Decision flow chart showing the protocols tested on eleven pairs towards an increasing reproducibility

The protocols with the best results using all 54 pairs were protocols 14 and 18 (Table 4 and Fig. 3b).

Protocols to compare for each question are listed in the flow chart of Fig. 4. Concerning the prior defined questions the results are as follows.

  1. 1.

    Is the reproducibility increased or decreased by the use of additional metabolites beside the clearly visible metabolites at 1.5 T?

    Question 1 is investigated by comparing identical protocols differing only in the use of GABA. Protocols 1, 3, 5, 7 used GABA as an additional metabolite, whereas the protocols 1a, 3a, 5a, 7a excluded GABA. The protocols without GABA gave a better reproducibility as indicated by a smaller sum of the reproducibility numbers of the five major metabolites (see Fig. 4). Also the mean χ2 is, except for protocol five, less in the protocols without GABA.

  2. 2.

    Is the reproducibility increased or decreased by the use of the sum of separately fitted glutamate and glutamine values instead of a predefined Glx model?

    The comparison shows, that the use of the sum of separately fitted glutamate and glutamine values yields in all cases a much better reproducibility than in the six cases of a predefined Glx model.

  3. 3.

    Is the reproducibility increased or decreased by the use of a “flat” or a more “curved” baseline using the first 10 ms or 25 ms for the base line calculation, respectively?

    Using only the protocols not ruled out by question 1 and 2 the use of a more flat base line is generally better (protocols 5a, 6, 9a, 10, 13a, 14, 18 and 19a are better than protocols 7a, 4, 11a, 8, 15a, 12, 16, 17a).

  4. 4.

    Is the reproducibility increased by the use of an independent damping for each metabolite or a common damping for the major metabolites?

    On the basis protocol 6, 10, 14, 18 no decision is possible (protocols 6 and 18 are slightly better than protocols 14 and 10).

  5. 5.

    Is the reproducibility increased by the use of a common damping for the major metabolites (see 4) except for Myo?

    The use of a freely adjustable damping for Myo and a common damping for Naa, Cre3, Cho and Glx gave less reproducible results than the inclusion of Myo in the common damping (compare Table 3 protocol 20 and 19a).

  6. 6.

    Is the reproducibility increased by the use of two additional lipid/macromolecule signals below 1.8 ppm as additional support for the baseline?

    On the basis protocol 5a, 6, 18, 19a no decision is possible (protocol 6 is much better than 5a and protocol 19a is not significantly better than 18).

  7. 7.

    Is the reproducibility increased or decreased by the use of the frequencies of glutamate and glutamine as independent fit parameters or the use of fixed frequencies in relation to the most prominent metabolite line of Naa?

    On the basis protocol 6, 10, 14, 18 no definite decision possible. Protocol 14 is slightly better than 18 and 6 is better than 10.

  8. 8.

    Is the reproducibility changed by the calculation of relative amplitudes (normalizing the metabolites using Cre3 as reference)?

    The normalization is a division of two amplitudes with errors. So it is expected that the reproducibility is rather decreased by normalization. The largest decrease is found for Naa. For Cho and Glx the reproducibility is slightly decreased, whereas for Myo no relevant increase is found (see Table 4; 14, 18 vs. 14 ref, 18 ref). The reproducibility for normalized Naa, Cho and Glx was around 5 % whereas the reproducibility for Myo was close but above 9 % (Table 4).

Protocols 6, 9a, 13a, 14, 18, 19a are left after this decision process.

These remaining protocols except 9a were investigated with all 54 pairs for the detection of the best reproducibility. 9a was replaced for total evaluation by the more reproducible protocols 10 and 12. Summarizing the results of Table 4 using all 54 pairs, the protocols with most reproducible results for Naa, Cre3, Cho, Myo and Glx in lesion free white matter (protocols 14 and 18) have the following common features: no additional metabolites, Glx is taken as the sum of separately fitted Glu and Gln, a baseline modeled with the points of the first 10 ms, no additional lipid lines, common damping of all metabolites.


This investigation uses in vivo data for reproducibility assessment. An alternative would be the use of a set of theoretical metabolite bases with added statistical noise. The advantage of in vivo data from eleven individuals at different brain sides and different time points lays in the realistic concentration variation of the metabolite basis and the immanent use of non “visible” metabolites with a poor signal to noise ratio.

In contrast to other studies on reproducibility or variability [1824] this investigation did not calculate the coefficient of variation for intra or inter subject reproducibility or variation depending on the location inside the brain. Instead two acquisitions with different noise and an identical metabolite basis are investigated for 54 pairs on their coincidence regarding the amplitudes of the main metabolites.

The protocols with the best results using all 54 pairs were protocols 14 and 18 (Table 4). Protocol 18 is favorable, mainly because protocol 18 achieves the same reproducibility as 14 with two parameters less: the frequencies of Glu and Glu were fixed in relation to the frequency of Naa. Additionally, in the reference version of protocol 18, 18ref, the reproducibility number rn of the pairs individually referenced to Cre3 as well as the sum of all rn’s were less for all cases except Cho in comparison to 14 ref.

The amplitudes of the metabolites strongly depend on the fitted baseline as also described by Hofmann et al. [25]. Depending on the curvature the values can differ significantly. In Fig. 3a for example the baseline had a broad peak between the lines of Cre2 and Cre3 and therefore the amplitude value (f. e.) for Cho is low in contrast to Fig. 3b with a more flat and lower baseline in this region.

This study did not vary all possible parameters. Additionally some parameters should have more steps. For example, using only the data of the first 10 or 25 ms for the base line were used; further investigations with additional times may improve the reproducibility. The default minimal signal to noise ratio for the amplitude of a metabolite was eight and never changed (see “Methods”). But every boundary in the signal to noise ratio–independent of the value–can impede the detection of a metabolite in only one of the two spectra building a pair. In such a case the normalized squared difference is large (four) assuming the not detected signal amplitude as zero (Eq. 1). This explains the large rn values for Glx in protocols 0 to 5a. Using protocols 0 to 3a (Table 5) in some cases the combined Glx was found only in one of the two spectra, whereas in case of protocol 4 to 5a this happened only for Gln with a lower signal to noise ratio compared to Glu, which was found in every case. A missing Gln in one of a pair increases the reproducibility number rn much less, than in the case of Glx.

Table 5 Reproducibility calculated as rn in percent for 54 pairs of consecutive spectra

The residual water suppression in this study was implemented as an iterative filter with 22 stages. The experience of the user (MB) preferred this value instead of the default value of 10. No other stage numbers have been tested, because recent post processing software algorithms (jMRUI [11, 2629] or Syngo Via (Siemens Healthcare)) prefer instead of an iterative filter a singular value decomposition (SVD) algorithm for residual water suppression. This was not possible with the work in progress software metabolite report, but both methods are available in the final product Syngo Via. The comparison of SVD and iterative filter within Syngo Via is not done yet.

This study did not examine absolute metabolite concentrations. Only the amplitudes or relative amplitudes with Cre3 as reference (ampmet/ampCre3) of directly comparable pairs were evaluated. The reason on one hand is the intention to search for an optimized reproducibility not needing a reference with a known concentration. On the other hand the reason is the lack of consensus for acquisition and post-processing methods to evaluate absolute metabolite concentration. Published results show “considerable variations in the approaches” [8] and therefore a large span for absolute concentrations. For the metabolite Naa with the best reproducibility, the concentration inside white matter in healthy subjects can differ from 9.7 to 18.4 mmol/l [3034]. The values depend on sequence design, acquisition system, post-processing algorithm as well as post-processing parameters, relaxation corrections, internal or external references or subject gender, subject age and the position inside the brain [24].

All spectra of this investigation were acquired in “lesion free” (no lesion signs found in MR imaging) white matter brain tissue. Spectra of brain lesions, for example tumors, can be drastically different. For example the lipid components in/around high-grade brain tumors normally are extremely large compared to the small (often zero) values found in healthy tissue. In our examination the detected lipids were very small and superimposed or dominated by macromolecules. In most of our pairs no or only one lipid line was found as expected for healthy subjects without a sign for any brain disease. For that reason the application of the post processing parameters of this investigation is limited to lesion free tissue. But despite of this limitation the application to studies is broad, because alterations in brain biochemistry of lesion free areas had been detected nearby tumors [3537], in patients with multiple sclerosis [38], drug abuse [3941], AIDS [42, 43] or psychiatric disorders [44].


The prerequisites of this study were optimal for in vivo comparisons. The brain is one of the organs that can be easily fixed inside the scanner. Spectrum pairs were acquired within a few minutes guaranteeing a stable metabolite basis for comparisons. All pairs were similar, because only lesion free white matter of healthy subjects at an identical brain position was investigated.

Metabolite amplitudes inside lesion free healthy brain can be quantified with a calculated accuracy (rn(met)) around 5 % for Naa, Cho and Glx using Cre3 as reference. The reproducibility of Myo is roughly twice as bad. Assuming these rn, the real metabolite amplitudes calculated with these post-processing parameters are most likely within a twofold interval of ±10 or ±20 % respectively. The reproducibility (rn) should be similar using other references like internal or external water for an absolute concentration evaluation and are not influenced by relaxation corrections with literature values, which is only a multiplication by a constant factor.


  1. 1.

    Barany M, Langer BG, Glick RP, Venkatasubramanian PN, Wilbur AC, Spigos DG. In vivo H-1 spectroscopy in humans at 1.5 T. Radiology. 1988;167:839–44.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Barany M, Spigos DG, Mok E, Venkatasubramanian PN, Wilbur AC, Langer BG. High resolution proton magnetic resonance spectroscopy of human brain and liver. Magn Reson Imaging. 1987;5:393–8.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Bruhn H, Frahm J, Gyngell ML, Merboldt KD, Hanicke W, Sauter R, et al. Noninvasive differentiation of tumors with use of localized H-1 MR spectroscopy in vivo: initial experience in patients with cerebral tumors. Radiology. 1989;172:541–8.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Frahm J, Bruhn H, Gyngell ML, Merboldt KD, Hanicke W, Sauter R. Localized proton NMR spectroscopy in different regions of the human brain in vivo. Relaxation times and concentrations of cerebral metabolites. Magn Reson Med. 1989;11:47–63.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Bruhn H, Michaelis T, Merboldt KD, Hanicke W, Gyngell ML, Hamburger C, et al. On the interpretation of proton NMR spectra from brain tumours in vivo and in vitro. NMR Biomed. 1992;5:253–8.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Frahm J, Michaelis T, Merboldt KD, Hanicke W, Gyngell ML, Bruhn H. On the N-acetyl methyl resonance in localized 1H NMR spectra of human brain in vivo. NMR Biomed. 1991;4:201–4.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Hennig J, Pfister H, Ernst T, Ott D. Direct absolute quantification of metabolites in the human brain with in vivo localized proton spectroscopy. NMR Biomed. 1992;5:193–9.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Watanabe M, Liao JH, Jara H, Sakai O. Multispectral quantitative MR imaging of the human brain: lifetime age-related effects. Radiographics. 2013;33:1305–19.

    Article  PubMed  Google Scholar 

  9. 9.

    Scheenen TW, Gambarota G, Weiland E, Klomp DW, Futterer JJ, Barentsz JO, et al. Optimal timing for in vivo 1H-MR spectroscopic imaging of the human prostate at 3T. Magn Reson Med. 2005;53:1268–74.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Naressi A, Couturier C, Castang I, de Beer R, Graveron-Demilly D. Java-based graphical user interface for MRUI, a software package for quantitation of in vivo/medical magnetic resonance spectroscopy signals. Comput Biol Med. 2001;31:269–86.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Naressi A, Couturier C, Devos JM, Janssen M, Mangeat C, de Beer R, et al. Java-based graphical user interface for the MRUI quantitation package. MAGMA. 2001;12:141–52.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Provencher SW. Estimation of metabolite concentrations from localized in vivo proton NMR spectra. Magn Reson Med. 1993;30:672–9.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Kobus T, Wright AJ, Weiland E, Heerschap A, Scheenen TW. Metabolite ratios in H MR spectroscopic imaging of the prostate. Magn Reson Med. 2014.

  14. 14.

    Scheenen TW, Futterer J, Weiland E, van Hecke P, Lemort M, Zechmann C, et al. Discriminating cancer from noncancer tissue in the prostate by 3-dimensional proton magnetic resonance spectroscopic imaging: a prospective multicenter validation study. Invest Radiol. 2011;46:25–33.

    Article  PubMed  Google Scholar 

  15. 15.

    E W: Automatische Quantifizierung von Metabolitenkonzentrationen in in vivo Spektren [PhD Thesis], Bremen 2004. 2004.

  16. 16.

    McLean MA, Simister RJ, Barker GJ, Duncan JS. Discrimination between neurochemical and macromolecular signals in human frontal lobes using short echo time proton magnetic resonance spectroscopy. Faraday Discuss. 2004;126:93–102. discussion 169-183.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Seeger U, Klose U, Mader I, Grodd W, Nagele T. Parameterized evaluation of macromolecules and lipids in proton MR spectroscopy of brain diseases. Magn Reson Med. 2003;49:19–28.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Wiebenga OT, Klauser AM, Nagtegaal GJ, Schoonheim MM, Barkhof F, Geurts JJ, et al. Longitudinal absolute metabolite quantification of white and gray matter regions in healthy controls using proton MR spectroscopic imaging. NMR Biomed. 2014;27:304–11.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Maudsley AA, Domenig C, Sheriff S. Reproducibility of serial whole-brain MR spectroscopic imaging. NMR Biomed. 2010;23:251–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Mullins PG, Rowland L, Bustillo J, Bedrick EJ, Lauriello J, Brooks WM. Reproducibility of 1H-MRS measurements in schizophrenic patients. Magn Reson Med. 2003;50:704–7.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Li BS, Babb JS, Soher BJ, Maudsley AA, Gonen O. Reproducibility of 3D proton spectroscopy in the human brain. Magn Reson Med. 2002;47:439–46.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Gasparovic C, Bedrick EJ, Mayer AR, Yeo RA, Chen H, Damaraju E, et al. Test-retest reliability and reproducibility of short-echo-time spectroscopic imaging of human brain at 3T. Magn Reson Med. 2011;66:324–32.

    Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Brooks WM, Friedman SD, Stidley CA. Reproducibility of 1H-MRS in vivo. Magn Reson Med. 1999;41:193–7.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Maudsley AA, Domenig C, Govind V, Darkazanli A, Studholme C, Arheart K, et al. Mapping of brain metabolite distributions by volumetric proton MR spectroscopic imaging (MRSI). Magn Reson Med. 2009;61:548–59.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Hofmann L, Slotboom J, Jung B, Maloca P, Boesch C, Kreis R. Quantitative 1H-magnetic resonance spectroscopy of human brain: Influence of composition and parameterization of the basis set in linear combination model-fitting. Magn Reson Med. 2002;48:440–53.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Gottschalk M, Lamalle L, Segebarth C. Short-TE localised 1H MRS of the human brain at 3 T: quantification of the metabolite signals using two approaches to account for macromolecular signal contributions. NMR Biomed. 2008;21:507–17.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Mazzoni LN, Belli G, Ginestroni A, Pratesi A, Agnoloni S, Diciotti S, et al. Computation of brain metabolite ratios in single-voxel proton MR spectroscopy: comparison between semiautomatic and automatic software. La Radiologia medica. 2010;115:125–32.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Pascual-Lozano AM, Martinez-Bisbal MC, Bosca-Blasco I, Valero-Merino C, Coret-Ferrer F, Marti-Bonmati L, et al. Casanova-Estruch B: [Total brain T2-hyperintense lesion-volume and the axonal damage in the normal-appearing white matter of brainstem in early lapsing-remitting multiple sclerosis]. Rev Neurol. 2007;45:468–73.

    CAS  PubMed  Google Scholar 

  29. 29.

    Pascual AM, Martinez-Bisbal MC, Bosca I, Valero C, Coret F, Martinez-Granados B, et al. Axonal loss is progressive and partly dissociated from lesion load in early multiple sclerosis. Neurology. 2007;69:63–7.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Bonekamp D, Smith MA, Zhu H, Barker PB. Quantitative SENSE-MRSI of the human brain. Magn Reson Imaging. 2010;28:305–13.

    Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Helms G, Stawiarz L, Kivisakk P, Link H. Regression analysis of metabolite concentrations estimated from localized proton MR spectra of active and chronic multiple sclerosis lesions. Magn Reson Med. 2000;43:102–10.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Gasparovic C, Yeo R, Mannell M, Ling J, Elgie R, Phillips J, et al. Neurometabolite concentrations in gray and white matter in mild traumatic brain injury: an 1H-magnetic resonance spectroscopy study. J Neurotrauma. 2009;26:1635–43.

    Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Tisell A, Leinhard OD, Warntjes JB, Lundberg P. Procedure for quantitative (1)H magnetic resonance spectroscopy and tissue characterization of human brain tissue based on the use of quantitative magnetic resonance imaging. Magn Reson Med. 2013;70:905–15.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Minati L, Aquino D, Bruzzone MG, Erbetta A. Quantitation of normal metabolite concentrations in six brain regions by in-vivoH-MR spectroscopy. J Med Phys. 2010;35:154–63.

    Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Busch M, Liebenrodt K, Gottfried S, Weiland E, Vollmann W, Mateiescu S, et al. Influence of brain tumors on the MR spectra of healthy brain tissue. Magn Reson Med. 2011;65:18–27.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Goebell E, Fiehler J, Ding XQ, Paustenbach S, Nietz S, Heese O, et al. Disarrangement of fiber tracts and decline of neuronal density correlate in glioma patients--a combined diffusion tensor imaging and 1H-MR spectroscopy study. AJNR Am J Neuroradiol. 2006;27:1426–31.

    CAS  PubMed  Google Scholar 

  37. 37.

    Inglese M, Brown S, Johnson G, Law M, Knopp E, Gonen O. Whole-brain N-acetylaspartate spectroscopy and diffusion tensor imaging in patients with newly diagnosed gliomas: a preliminary study. AJNR Am J Neuroradiol. 2006;27:2137–40.

    CAS  PubMed  Google Scholar 

  38. 38.

    Tisell A, Leinhard OD, Warntjes JB, Aalto A, Smedby O, Landtblom AM, et al. Increased concentrations of glutamate and glutamine in normal-appearing white matter of patients with multiple sclerosis and normal MR imaging brain scans. PLoS ONE. 2013;8:e61817.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Cowan RL, Joers JM, Dietrich MS. N-acetylaspartate (NAA) correlates inversely with cannabis use in a frontal language processing region of neocortex in MDMA (Ecstasy) polydrug users: a 3 T magnetic resonance spectroscopy study. Pharmacol Biochem Behav. 2009;92:105–10.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Gazdzinski S, Durazzo TC, Mon A, Yeh PH, Meyerhoff DJ. Cerebral white matter recovery in abstinent alcoholics–a multimodality magnetic resonance study. Brain. 2010;133:1043–53.

    Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Licata SC, Renshaw PF. Neurochemistry of drug action: insights from proton magnetic resonance spectroscopic imaging and their relevance to addiction. Ann N Y Acad Sci. 2010;1187:148–71.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Lee PL, Yiannoutsos CT, Ernst T, Chang L, Marra CM, Jarvik JG, et al. A multi-center 1H MRS study of the AIDS dementia complex: validation and preliminary analysis. J Magn Reson Imaging. 2003;17:625–33.

    Article  PubMed  Google Scholar 

  43. 43.

    Paley M, Cozzone PJ, Alonso J, Vion-Dury J, Confort-Gouny S, Wilkinson ID, et al. A multicenter proton magnetic resonance spectroscopy study of neurological complications of AIDS. AIDS Res Hum Retrovir. 1996;12:213–22.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Bernier D, Cookey J, McAllindon D, Bartha R, Hanstock CC, Newman AJ, et al. Multimodal neuroimaging of frontal white matter microstructure in early phase schizophrenia: the impact of early adolescent cannabis use. BMC Psychiatry. 2013;13:264.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


The knowledge of the software “metabolite report” is from Elisabeth Weiland (Siemens Healthcare, Erlangen, Germany), one of the developers of the software, through personal communications. Thanks for the short description of the software in “Methods” and many interesting personal discussions.

Author information



Corresponding author

Correspondence to Martin H J Busch.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MHJB: study design, spectra analysis, statistical analysis, manuscript drafting and revising. WV: spectra analysis, statistical analysis, manuscript revising. SM: spectra acquisition, manuscript revising. MS: spectra acquisition, spectra analysis. MD: spectra acquisition, manuscript revising. MG: spectra analysis, manuscript revising. DHWG: manuscript drafting and revising. All authors read and approved the final manuscript.

Authors’ information

Not applicable.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Busch, M.H.J., Vollmann, W., Mateiescu, S. et al. Reproducibility of brain metabolite concentration measurements in lesion free white matter at 1.5 T. BMC Med Imaging 15, 40 (2015).

Download citation


  • Post Processing
  • Spectrum Pair
  • Metabolite Basis
  • Iterative Filter
  • Absolute Metabolite Concentration