|Year : 2010 | Volume
| Issue : 3 | Page : 154-163
Quantitation of normal metabolite concentrations in six brain regions by in-vivo 1 H-MR spectroscopy
Ludovico Minati1, Domenico Aquino2, Maria Grazia Bruzzone2, Alessandra Erbetta2
1 Scientific Department and Neuroradiology Unit, Fondazione IRCCS Istituto Neurologico Carlo Besta, Milano, Italy
2 Neuroradiology Unit, Fondazione IRCCS Istituto Neurologico Carlo Besta, Milano, Italy
|Date of Submission||30-Oct-2009|
|Date of Decision||18-Dec-2009|
|Date of Acceptance||24-Jan-2010|
|Date of Web Publication||26-Jul-2010|
Scientific Department and Neuroradiology Unit, Fondazione IRCCS Istituto Neurologico Carlo Besta, via Celoria 11, 20133 Milano
Source of Support: None, Conflict of Interest: None
| Abstract|| |
This study examined the concentrations of brain metabolites visible to in-vivo 1 H-Magnetic Resonance Spectroscopy ( 1 H-MRS) at 1.5 T in a sample of 28 normal subjects. Quantitation was attempted for inositol compounds, choline units, total creatine and N-acetyl moieties, using open-source software. Six brain regions were considered: frontal and parietal white matter, medial temporal lobe, thalamus, pons and cerebellum. Absolute concentrations were derived using tissue water as an internal reference and using an external reference; metabolite signal intensity ratios with respect to creatine were also calculated. The inter-individual variability was smaller for absolute concentrations (internal reference) as compared to that for signal intensity ratios. Significant regional variability in concentration was found for all metabolites, indicating that separate normative values are needed for different brain regions. The values obtained in this study can be used as reference in future studies, provided the same methodology is followed; it is confirmed that despite unsuccessful attempts in the past, smaller coefficients of variation can indeed be obtained through absolute quantification.
Keywords: Brain metabolites, magnetic resonance spectroscopy, normal concentrations, quantitation
|How to cite this article:|
Minati L, Aquino D, Bruzzone MG, Erbetta A. Quantitation of normal metabolite concentrations in six brain regions by in-vivo 1 H-MR spectroscopy. J Med Phys 2010;35:154-63
|How to cite this URL:|
Minati L, Aquino D, Bruzzone MG, Erbetta A. Quantitation of normal metabolite concentrations in six brain regions by in-vivo 1 H-MR spectroscopy. J Med Phys [serial online] 2010 [cited 2020 Mar 31];35:154-63. Available from: http://www.jmp.org.in/text.asp?2010/35/3/154/62128
| Introduction|| |
In proton magnetic resonance spectroscopy ( 1 H-MRS), the signal intensity i s not only proportional to metabolite concentration but is also affected by a large number of variables, including metabolite relaxation rates, pulse sequence parameters and radio-frequency coil sensitivity. As a consequence of the inherent difficulties in obtaining reliable absolute concentrations, the majority of in-vivo MR spectroscopy studies of the brain report relative measurements, calculated assuming that the concentration of creatine is constant and taking its signal as reference. However, the interpretation of relative data can be ambiguous when the concentration of creatine changes, such as in early development and in brain tumors. , In addition to removing interpretation ambiguities, the provision of absolute measurements considerably simplifies the comparison of data recorded under different experimental conditions and the comparison with in-vitro biochemical measurements. 
Absolute metabolite concentrations may be derived using an internal reference, usually unsuppressed tissue water, or an external reference, generally water or sodium acetate contained in a tube positioned next to the head of the patient. Both methods have important limitations: the former can suffer from bias due to altered tissue water content, for example caused by edema, and the latter is associated with large inter-subject variability due to coil sensitivity inhomogeneities. Other methods based on replacement phantoms and coil-loading measurements also exist but are less frequently used. ,,,,
In order to obtain absolute concentrations, it is also necessary to determine the relaxation rates of the metabolites of interest in order to correct for transverse and longitudinal magnetization effects. Direct measurement with relaxometry is the most accurate option but is often impractical due to scan time limitations. As a consequence, metabolite relaxation rates are frequently assumed to be constant, and standard values from literature are used; this can be a source of error due to pathology-related changes and regional differences. ,,,
Despite these difficulties, in recent years quantitative 1 H-MRS of the brain has gained increased acceptance in the clinical domain and has been applied to the study of aging, senile dementia, epilepsy, multiple sclerosis and neuropsychiatric disorders, among others; it provides a tool for quantitative monitoring of disease progression and treatment response and can support differential diagnosis between conditions characterized by similar imaging findings. ,,,,
There remains, however, a relative paucity of normative studies covering multiple cortical and subcortical brain regions, as well as a large inter-study variability in reported metabolite concentrations (for the cerebral hemispheres, in the range 4-8 mM for inositol compounds, 1-5 mM for choline units, 6-14 mM for total creatine, and 10-25 mM for N-acetyl moieties) and inter-subject variation coefficients (between about 10% and more than 50%). ,,,, It has been shown that differences in data processing are a dominant source of inter-study variability; for example, LCModel and AMARES, two widely used commercially available spectroscopy toolkits, may perform differently in terms of sensitivity to noise, linewidth and baseline. ,
Recently, an open-source tool for quantitation of short echo-time spectra, known as AQSES, has become available. In AQSES, the quantitation problem is formulated as a separable nonlinear least-squares fitting problem, and facilities for baseline fitting, removal of residual water and filtering are provided. 
In order to corroborate and extend the existing literature, in this study we determined the normal absolute and relative metabolite concentrations in an extended set of brain regions: frontal and parietal white matter, medial temporal lobe, thalamus, pons and cerebellum. Importantly, another aim was to compare the inter-subject variability of absolute concentrations calculated using the internal and external references with that of relative concentrations. This was done with the purpose of determining whether performing absolute quantification is truly justified, in spite of the need to apply multiple corrections involving parameters which inject measurement uncertainty.
| Materials and Methods|| |
Twenty-eight right-handed healthy volunteers, 14 female and 14 male, aged 41.2±11.9 years (mean±SD), were enrolled. All subjects, unpaid, were informed about the purpose and clinical relevance of the study, and written informed consent was obtained from them according to institutionally approved procedures and regulations. None of them had a positive anamnesis for neurological or psychiatric disorders, and all had normal magnetic resonance imaging findings (see below).
Data acquisition was performed with a Siemens Magnetom Avanto 1.5 T scanner (Siemens AG, Erlangen, Germany) equipped with a transceiver birdcage head coil. For the purpose of positioning the voxels for spectroscopy and to exclude pathology, all subjects were imaged with axial, coronal and sagittal T 2 -weighted turbo spin-echo sequences, with TR=6270 ms, TE=113 ms, field-of-view 260Χ190 mm, 25 slices, slice thickness 4 mm and no gap. A coronal fluid-attenuated inversion-recovery (FLAIR) sequence was also used, with TR=8700 ms, TE=121 ms, TI=2400 ms, field-of-view 240Χ190 mm, 20 slices, slice thickness 4 mm and no gap. High-resolution T 1 -weighted volumetric images were acquired by means of an inversion-recovery rapid gradient-echo sequence with TR=1160 ms, TE=4.1 ms, TI=600 ms, field-of-view 270Χ203 mm, 192 slices, slice thickness 1.1 mm and no gap.
Spectroscopy of the 1 H nucleus was performed using a single-voxel (SVS) point-resolved spectroscopy (PRESS) sequence, with TR=1500 ms, TE=30 ms and voxel size 20Χ20Χ20 mm [with the exception of the medial temporal lobe region (see below), for which it was 40Χ10Χ20 mm]; for each voxel, 1024 data points were acquired with a dwell time of 1 ms and 128 averages. Metabolite spectra were obtained suppressing the water signal by means of a chemical shift-s elective saturation (CHESS) pulse with bandwidth 35 Hz. For each spectroscopy voxel, the signal of unsuppressed water was also acquired, with 16 averages. First- (X, Y and Z), as well as second-order (XY, ZX, XY, Z 2 and X 2 -Y 2 ), shimming was performed automatically by means of a field map-based algorithm and then refined manually by an experienced operator; the width at half height of the water peak was generally below 10 Hz (7±2.2 Hz). A 50-mL polystyrene tube containing distilled water, positioned in contact with the head next to the left earlobe, served as external reference; the corresponding water signal was acquired using a 10Χ10Χ10-mm voxel.
The static field homogeneity, coil sensitivity maps, tuning and signal-to-noise ratio, radio-frequency chain linearity and gradient performance were routinely checked using the standard software and procedures provided by the scanner manufacturer, with the prescribed periodicity.
Six spectroscopy voxels were positioned in the following regions, shown in [Figure 1], together with the corresponding stereotactic coordinates given in Montreal Neurological Institute (MNI) format: frontal white matter (FWM), parietal white matter (PWM), medial temporal lobe (MTL), thalamus (THAL), pons (PONS) and cerebellar hemisphere (CEREB). Notably, as a consequence of choosing equal volume of 8 mL for all voxels, it was not possible to position the voxel entirely within the pons and thalamus; as visible in [Figure 1], a contribution of surrounding structures was also present. The Niewenhuys Atlas More Details  served as anatomical reference. The spectroscopy voxels were normally positioned in the left hemisphere; however, for 9 subjects the cerebellar voxel was placed in the right hemisphere because magnetic susceptibility effects resulted in a water peak width > 10 Hz in the left hemisphere.  For 19 subjects, spectra from the right medial temporal lobe were also available; they were not acquired for all subjects due to variable scan-time limitations.
Preprocessing and quantitation of 1 H spectra
Spectra were processed by means of the freely available AQSES software (Katholieke Universiteit Leuven, Leuven, Belgium) running under MatLab 7 (The MathWorks Inc., Natick, MA, USA) on a Sun Ultra 80 workstation (Sun Microsystems Inc., Palo Alto, CA, USA).  The signal was zero-padded to 4096 points, apodized with a Gaussian function and Fourier-transformed. Zero-order rephasing and frequency correction were performed maximizing the symmetry of the metabolite peaks and centering the N-acetyl moieties peak (see below) on 2.02 ppm. The residual water signal (4.7 ppm) was removed in the range 4.3-5.1 ppm by means of the Hankel-Lanczos singular-value decomposition (HLSVDpro).  Remaining baseline fluctuations were removed by fitting with a fifth-degree polynomial. As shown in [Figure 2], the metabolite spectrum was fit in the range 1.2-4.3 ppm with the linear superposition of five components, corresponding to inositol compounds (MI), choline units (CHO), creatine and phosphocreatine (CR), glutamate and glutamine (GLX), and N-acetyl (NA) moieties. The intensity of the metabolite signals was determined integrating the following resonance peaks: 3.50 ppm for MI, 3.20 ppm for CHO, 3.00 ppm for CR, 2.70-1.70 (range) ppm for GLX, and 2.02 ppm for NA. ,, The basis set, provided with the software, was derived from phantom measurements. 
Volumetric images were normalized to the MNI space and segmented using the SPM5 program (Wellcome Neuroimaging Department, London, UK), and the relative content of white matter (WM), gray matter (GM) and cerebrospinal fluid (CSF) in each spectroscopy voxel was determined.
Absolute concentrations, expressed in millimolarities (mmol/L, mM), were calculated for MI, CHO, CR and NA (not for GLX, due to the inability to separate glutamate and glutamine at a field strength of 1.5 T), using the intensity of unsuppressed water signal in the spectroscopy voxel (internal reference) or in the polystyrene tube (external reference) as reference. Metabolite concentrations cmet were determined using
where kmet and KH2O are correction factors for T 1 recovery and T 2 relaxation of metabolite and water (see below), Smet and SH2O are the signal intensities, nmet and nH2O are the numbers of protons (4 for MI, 9 for CHO, 3 for CR and NA and 2 for water), fCSF is the fraction of cerebrospinal fluid in the spectroscopy voxel, Vmet is the volume of the spectroscopy voxel (8 mL), VH2O is the volume of the voxel used for water signal acquisition (8 mL for internal reference, 1 mL for external reference), and cH2O is the reference water concentration.
The correction factor for metabolite, kmet , was determined using
where fGM and fWM are the fractions of gray and white matter in the spectroscopy voxel, and , T1,met,GM, T2,met,GM, T1,met,WM, and T2,met,WM are the relaxation times for the metabolite, obtained from  and provided in [Table 1]. Notably, the factor 1/(1-fCSF) appears in both eq. 1 and eq. 2; whereas in eq. 1 it accounts for the fact that metabolite signal is not received from CSF, in eq. 2 it accounts for the fact that the weighing factors for the relaxation terms, fGM and fWM , may not sum to 1.
For quantitation with the internal reference, the correction factor for water, kH2O , was determined using
where T1,H2O,GM, T2,H2O,GM, T1,H2O,WM, T2,H2O,WM, T1,H2O,CSF, and T2,H2O,CSF are the relaxation times for water, provided in [Table 1]. For quantitation with the external reference, in the calculation of kH2O the relaxation times of cerebrospinal fluid were used, with fGM=fWM=0 and fCSF=1, and kH2O was multiplied by a factor of 0.80, determined from coil sensitivity maps acquired during quality assessment checks, to account for the difference in coil sensitivity between the two positions: center of the head and coil and average position of the polystyrene tube.
For quantitation with the internal reference, the concentration of water in the voxel, , was determined using
where cH2O,GM, cH2O,WM, and cH2O,CSF are the concentrations of water in gray matter, white matter and cerebrospinal fluid, provided in [Table 1].  For quantitation with the external reference, in the calculation of cH2O we set
In addition to the absolute concentrations, the signal intensity ratios MI/CR, CHO/CR, GLX/CR and NA/CR were also obtained. As these were not normally distributed, the ratios were logarithm-transformed.
For each parameter of interest, mean and standard deviation w ere computed. To ensure that there were no significant deviations from Gaussian distribution, Kolmogorov-Smirnov normality tests were performed.
Inter-subject coefficients of variation (CVs) were calculated for the absolute concentrations and for the signal intensity ratios using the formula CV = SD / mean, and averaged for each metabolite across the six regions. To test for systematic differences between the metabolite concentrations calculated with the internal and external references, paired t-tests were performed.
For those subjects for whom spectra were available for both the left and the right MTL, lateralization indices (LIs) for the concentrations of MI, CHO, CR and NA, determined using the internal reference, were computed with the formula
where cL and cR refer to the concentrations in the left and right voxels; for each metabolite, the statistical significance of lateralization was evaluated by means of paired t-tests.
To test for differences in metabolite concentrations (internal reference) and signal intensity ratios among brain regions, analysis of variance (ANOVA) was performed, followed by Bonferroni-corrected post-hoc t-tests.
| Results|| |
The absolute metabolite concentrations and intensity ratios are given in [Table 2], and the corresponding bar charts are shown in [Figure 3]. All concentrations and logarithm-transformed intensity ratios were normally distributed. Across the spectra of the participants, the following average metabolite signal-to-noise ratios (SNRs) were obtained: 32.1±5.1 for CEREB, 28.7±6.6 for FWM, 26.8±6.7 for MTL, 27.3±7.5 for PONS, 28.6±6.9 for PWM and 26.2±4.8 for THAL.
There were no significant differences between the average metabolite concentrations calculated with the internal and external references. The average CVs for concentrations determined with the internal reference were 32% for MI, 23% for CHO, 20% for CR and 15% for NA. In comparison, the CVs obtained with the external reference were larger: 50% for MI, 44% for CHO, 41% for CR and 41% for NA. The CVs for the corresponding signal intensity ratios were intermediate: 41% for MI/CR, 27% for CHO/CR, 28% for NA/CR and 48% for GLX/CR.
The lateralization indices for the MTL were 0.01±0.18 for MI, 0.05±0.18 for CHO, - 0 .01±0.18 for CR and 0.01±0.10 for NA; there was no significant lateralization.
The difference in metabolite concentrations (internal reference) among regions was largest for NA (F5,27 =16.1, P< .001), intermediate for CR (F5,27 =13.5, P< .001) and CHO (F5,27 =11.0, P< .001), and smallest for MI (F5,27 =5.0, P< .001). The concentration of MI ranged between 5.8 mM (PWM) and 8.8 mM (PONS), that of CHO ranged between 2.9 mM (PWM) and 4.7 mM (PONS), that of CR ranged between 10.5 mM (PONS) and 15.7 mM (CEREB), and that of NA ranged between 14.0 mM (PWM) and 18.4 mM (PONS).
The difference in transformed intensity ratios among regions was largest for CHO/CR (F5,27 =19.3, P< .001), intermediate for NA/CR (F5,27 =13.9, P< .001) and MI/CR (F5,27 =7.4, P< .001), and smallest for GLX/CR (F5,27 =5.2, P< .001). The transformed intensity ratio ln(MI/CR) ranged between -0 .35 (CEREB) and 0.15 (PONS), ln(CHO/CR) ranged between −0.29 (CEREB) and 0.33 (PONS), ln(GLX/CR) ranged between 1.07 (THAL) and 1.56 (PONS), and ln(NA/CR) ranged between 0.10 (CEREB) and 0.66 (PONS).
The results of post-hoc comparisons between regions are given in [Table 3].
| Discussion|| |
While there were no systematic differences between the metabolite concentrations estimated with the internal and external references, the external reference was associated with considerably larger inter-individual variability, reflected in coefficients of variation about twice as large. This finding is in line with the results obtained at the majority of centers that participated in a previous multi-centric study.  Imperfect repeatability in the positioning of the reference tube is the main source of random error, and its effect is amplified by changes in the coil sensitivity profile between subjects, caused by varying shape and size of the head. It therefore appears generally preferable to use the internal reference unless alterations in tissue water content are expected; more successful results with external standards have, however, been obtained by some groups. ,,
The inter-subject coefficients of variation were larger (by about 40%) for the signal intensity ratios than for the corresponding absolute concentrations estimated with the internal reference. This likely results from the combination of two factors: 1) the actual inter-individual variability in concentration is lower for water (as confirmed by proton density-mapping studies) than for creatine, and 2) the water peak is much larger than the creatine peak and therefore associated with a smaller relative measurement error. 
One common criticism to absolute quantification is that it involves the use of multiple correction parameters, namely, cH2O, fGM, fWM, fCSF, T1,met,GM T2,met,GM, T1,met,WM, T2,met,WM which carry uncertainties that are difficult to estimate and which get propagated in the final measurement result. Although this is an important shortcoming, our findings on the inter-subject coefficients of variation (representing random error) demonstrate that performing absolute quantification is nevertheless motivated, potentially leading to increased sensitivity to pathological change when compared to metabolite ratios.
Several studies reported metabolite concentrations in the cerebral hemispheres: the average concentrations obtained in the present work (7 mM for MI, 3 mM for CHO, 11 mM for CR and 14 mM for NA) and the corresponding inter-individual CVs are within the spread of values found in literature (4-8 mM for MI, 1-5 mM for CHO, 6-14 mM for CR and 10-25 mM for NA). ,,,,
The NA peak is frequently assumed to correspond to N-acetyl aspartate, a compound of interest as a putative neuronal marker. Ex-vivo studies of the human and rat brain have reported concentrations in the range 5-8 mM, significantly below the estimates obtained in this study and in most other in-vivo MRS studies. , It is hypothesized that this discrepancy arises because of unresolved compounds contributing to the 2.02 ppm NA resonance observed at 1.5 T, and possibly due to N-acetyl aspartate loss during ex-vivo sample preparation. , The CR peak at 3.00 ppm corresponds to the superposition of creatine and phosphocreatine, which have well-known roles in cell metabolism, in the concentration ratio determined by phosphokinase equilibrium.  Ex-vivo studies of the human and canine brain have reported concentrations in the range 9-12 mM, in line with the results of this study and other in-vivo MRS studies. ,,, The CHO peak at 3.20 ppm is the most complex, receiving contributions from a range of choline-containing compounds, including phosphocholine, glycerophosphocholine, free choline, acetylcholine, phosphatidylcholine and choline-plasmalogen; its intensity is frequently taken as an empirical marker of the density and turnover of cell membranes. , In vitro, the concentration of phosphocholine, glycerophosphocholine plus free choline was determined to be in the range 1-2 mM, in line only with some in-vivo studies but below the estimates of others (including the present one); the determinants of this discrepancy remain unclear. ,,,, The MI peak at 3.50 ppm corresponds to a range of compounds, including phosphatidylinositol, inositol polyphosphatide, inositol monophosphate, myo-inositol and, to a smaller extent, glycine; as inositol is elevated within astrocytes, the intensity of the peak is often taken as an empirical marker of glial density and proliferation. Studies on human brain samples have reported concentrations in the range 5-7 mM, in line with the results of this study and other in-vivo MRS studies. ,,, The GLX complex in the range 2.70-1.70 ppm includes multiple overlapping resonances from glutamate and its precursor glutamine; as these cannot be resolved at 1.5 T, quantitation was not attempted. 
In agreement with the existing literature, significant regional differences were found for all metabolites of interest. ,, These were largest for NA, intermediate for CR and CHO, and smallest for MI. Given that N-acetyl aspartate is present in the soma of neurons, in dendrites and in axons, its regional variability is likely related to differences in neural architecture, population and density; a simple linear relationship with the density of neurons is, however, unlikely given that it also reflects reversible metabolic changes.  Total creatine and choline are less specific as they include contributions from both neurons and glia, and their regional variability is probably related to differences in the density of the cellular matrix; and in the case of choline, also in the level of myelination.  The variability in the concentration of inositol compounds could be more specifically related to the glial population, but it is likely also influenced by regional metabolic differences. 
The concentration of N-acetyl aspartate was higher in the thalamus, pons and cerebellum than in the cerebral hemispheric regions; this effect could be related to higher density of neural soma, axons and dendritic trees in these regions, and is found for the thalamus also in  but not in.  The concentration of choline units was markedly higher in the pons than in other regions; this effect could be related to the high myelination of the dense rostro-caudal pathways in this region  and is also found in  and . In agreement with the same studies, the concentration of choline units was lowest in the parietal lobe; the same trend was observed for inositol, and the interpretation is unclear. The concentration of creatine was markedly higher in the cerebellum than in other regions; this effect could be related to the very high density of neural cells in the cerebellum  and is also found in  and  . The presence of significant regional differences signals the need to have separate reference values for each region of interest, especially for subcortical structures.
The concentration of metabolites in the hippocampus was found to be symmetrical, in line with  significant asymmetries emerge in epilepsy, even in absence of frequent seizures. 
The present study has three main strengths. First, in contrast with several previous ones which considered specific regions only, metabolite concentrations were measured, for all subjects, in an extended set of cortical and subcortical brain regions. Second, high-resolution segmentation was performed to correct for cerebrospinal fluid partial voluming and to report the white/ gray matter contents of each voxel. Third, a freely available fitting program was used, enabling centers that do not have access to commercial spectroscopy software to reproduce our findings.
There are, however, also several important limitations that need to be taken into consideration. First, the metabolite relaxation rates were not measured, due to scan-time limitations related to the large number of voxels under study, and the values from a previous study were used. While this is a potential source of error, it should be noted that the inter-study variability in relaxation times is considerably smaller (CV< 20%) in comparison with that in concentrations.  There are known differences in relaxivities among brain regions, which are less marked at 1.5 T than at higher field strengths; these were indirectly taken into account, determining "effective" T 1 and T 2 for each metabolite on the basis of the content of each voxel and of the average T 1 and T 2 values for white and gray matter. ,,, Although this approach cannot reveal pathology-induced changes in the metabolite T 2 s as found, for example, in stroke, the potential confounding effect was minimized because the T 2 -weighting was very weak, given that the MRS echo time (30 ms) was considerably shorter than the metabolite T 2 s (200-500 ms). , Second, the uncertainty associated with each correction parameter (cH2O, fGM, fWM, fCSF, T1,met,GM, T2,met,GM, T1,met,WM, T2,met,WM) remains unknown, and the finding of reduced inter-subject coefficients of variation does not exclude the presence of significant systematic error, even though there were no systematic differences between internal and external references. This is a general limitation of the technique at the present stage, which needs addressing through dedicated relaxometry and segmentation studies. Our findings of reduced inter-subject variability motivate such studies. Third, with the exception of the medial temporal lobe, voxels were positioned in one hemisphere only, preventing the determination of potential metabolic asymmetries related, for example, to the known difference in axonal density in some regions between the dominant and the non-dominant hemispheres; this limitation is in common with many other similar studies. ,,, Fourth, the number of subjects and the age range were too small to investigate age-related changes. The age range under consideration is, however, relevant to a wide range of pathologies, and previous studies indicate that the age-related changes in metabolite concentrations are negligible until the seventh decade.  Another limitation is that the study was conducted at 1.5 T, while the use of a 3-T scanner would have provided better SNR and peak separation, reducing quantitation uncertainty; however, in the majority of centers, clinical spectroscopy is still routinely practiced at 1.5 T.  Due to the inability to resolve glutamine and glutamate resonances at 1.5 T, the data on the GLX/CR ratio should be interpreted with caution. Finally, the voxel size (8 mL) was comparatively large, especially considering the pons and thalamus regions. Although we have corrected for partial voluming with fluid-containing spaces, inclusion of gray and white matter from surrounding areas affected the anatomical specificity of the measurements. Reducing voxel volume, e.g., to 4 mL, would increase anatomical specificity, as well as limit susceptibility effects. In this study, the choice of 8 mL was made in order to maximize SNR, and all voxels had equal volume to avoid introducing SNR differences across the regional spectra. It has been previously shown that reducing the SNR increases random error and may also result in systematic concentration biases. , While comprehensive evaluations of the effects of SNR are available for LCModel and AMARES, at present these remain lacking for AQSES.  Additional work is therefore necessary both to determine how SNR changes affect the accuracy of concentrations calculated using the proposed method, and to derive further reference values using smaller voxels.
More generally, a limitation of the present study, which is in common with all other similar normative studies, ,,,, is that the reference values are valid only so long as the same methodology and conditions are applied. Replication of the quantification procedure employed in the present study appears unproblematic as a full description has been g iven and the AQSES software is freely available. The exact position of the voxels in normalized space was provided [Figure 1], and the acquisition parameters are standard for clinical MRS at 1.5 T. Nevertheless, there are a number of potential sources of error, such as differences in signal-to-noise ratio, shimming quality and pulse sequence implementation, as a consequence of which our findings cannot immediately be used as normal reference even in absence of explicit methodological differences; rather, they should initially be compared with locally acquired control values and included as reference only after confirming the absence of statistically significant differences. For situations in which the same quantification methodology cannot be followed, for example when a different value of field strength is used, the present values cannot be considered and the relevance of the present work is purely methodological.
In conclusion, we have reported the normal in-vivo concentrations of inositol compounds, choline units, total creatine and N-acetyl moieties in frontal and parietal regions, medial temporal lobe, thalamus, pons and cerebellum. These values can be used as reference for future studies provided the same methodology, based on open-source software, is followed. The inter-individual coefficients of variation were largest for quantitation based on an external reference, discouraging its general use; there was, however, no systematic bias with respect to quantitation based on tissue water. The inter-individual coefficients of variation were larger for metabolite ratios than for absolute concentrations, confirming the potential usefulness of quantitative MRS in spite of the complexity of the correction process. Significant differences in metabolite concentrations were found among brain structures. While their exact determinants cannot be determined by in-vivo studies such as the present one, they are hypothesized to be related to regional variability in neural and glial population, and in myelination. Separate reference values are needed for different brain regions.
| Acknowledgments|| |
We are grateful to Emma Lisa Stanley for her operational contribution in the preprocessing and fitting of the spectra. We would also like to express our gratitude to Dr. Fausto Cevolani, who supported the preparation of some phantoms that were used for preliminary tests. This work is dedicated to his memory. We thank two anonymous reviewers for providing insightful feedback on an earlier version of the manuscript.
| References|| |
|1.||Tong ZY, Toshiaki Y, Wang YJ. Proton magnetic resonance spectroscopy of normal human brain and glioma: A quantitative in vivo study. Chin Med J 2005;118:1251-7. |
|2.||Dezortova M, Hajek M. (1)H MR spectroscopy in pediatrics. Eur J Radiol 2008;67:240-9. |
|3.||Barker PB, Soher BJ, Blackband SJ, Chatham JC, Mathews VP, Bryan RN. Quantitation of proton NMR spectra of the human brain using tissue water as an internal concentration reference. NMR Biomed 1993; 6:89-94. |
|4.||Christiansen P, Henriksen O, Stubgaard M, Gideon P, Larsson HB. In vivo quantification of brain metabolites by 1H-MRS using water as an internal standard. Magn Reson Imag 1993;11:107-18. |
|5.||Kreis R, Ernst T, Ross BD. Absolute quantitation of water and metabolites in the human brain. II: Metabolite concentrations. J Magn Reson 1993;102:9-19. |
|6.||Michaelis T, Merboldt KD, Bruhn H, Hδnicke W, Frahm J. Absolute concentrations of metabolites in the adult human brain in vivo: Quantification of localized proton MR spectra. Radiology 1993;187:219-27. |
|7.||Keevil SF, Barbiroli B, Brooks JC, Cady EB, Canese R, Carlier P, et al. Absolute metabolite quantification by in vivo NMR spectroscopy: II, A multicentre trial of protocols for in vivo localised proton studies of human brain. Magn Reson Imaging 1998;16:1093-106. |
|8.||Tong Z, Yamaki T, Harada K, Houkin K. In vivo quantification of the metabolites in normal brain and brain tumors by proton MR spectroscopy using water as an internal standard. Magn Reson Imag 2004;22:1017-24. |
|9.||Rutgers DR, van der Grond J. Relaxation times of choline, creatine and N-acetyl aspartate in human cerebral white matter at 1.5 T. NMR Biomed 2002;15:215-21. |
|10.||Brief EE, Whittall KP, Li DK, MacKay A. Proton T1 relaxation times of cerebral metabolites differ within and between regions of normal human brain. NMR Biomed 2003;16:503-9. |
|11.||Trδber F, Block W, Lamerichs R, Gieseke J, Schild HH. 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. J Magn Reson Imaging 2004;19:537-45. |
|12.||Walker PM, Ben Salem D, Lalande A, Giroud M, Brunotte F. Time course of NAA T2 and ADCw in ischaemic stroke patients: 1H MRS imaging and diffusion-weighted MRI. J Neurol Sci 2004;220:23-8. |
|13.||Brex PA, Parker GJ, Leary SM, Molyneux PD, Barker GJ, Davie CA, et al. Lesion heterogeneity in multiple sclerosis: A study of the relations between appearances on T1 weighted images, T1 relaxation times, and metabolite concentrations. J Neurol Neurosurg Psychiatry 2000;68:627-32. |
|14.||Brooks JC, Roberts N, Kemp GJ, Gosney MA, Lye M, Whitehouse GH. A proton magnetic resonance spectroscopy study of age-related changes in frontal lobe metabolite concentrations. Cereb Cortex 2001;11:598-605. |
|15.||Varho T, Komu M, Sonninen P, Lδhdetie J, Holopainen IE. Quantitative HMRS and MRI volumetry indicate neuronal damage in the hippocampus of children with focal epilepsy and infrequent seizures. Epilepsia 2005;46:696-703. |
|16.||Zhu X, Schuff N, Kornak J, Soher B, Yaffe K, Kramer JH, et al. Effects of Alzheimer disease on fronto-parietal brain N-acetyl aspartate and myo-Inositol using magnetic resonance spectroscopic imaging. Alzheimer Dis Assoc Disord 2006;20:77-85. |
|17.||Aydin K, Ucok A, Cakir S. Quantitative proton MR spectroscopy findings in the corpus callosum of patients with schizophrenia suggest callosal disconnection. AJNR Am J Neuroradiol 2007;28:1968-74. |
|18.||De Beer R, Barbiroli B, Gobbi G, Knijn A, Kόgel H, Langenberger KW, et al. Absolute metabolite quantification by in vivo NMR spectroscopy: III, Multicentre 1H MRS of the human brain addressed by one and the same data-analysis protocol. Magn Reson Imaging 1998;16:1107-11. |
|19.||Kanowski M, Kaufmann J, Braun J, Bernarding J, Tempelmann C. Quantitation of simulated short echo time 1H human brain spectra by LCMODEL and AMARES. Magn Reson Med 2004;51:904-12. |
|20.||Poullet JB, Sima DM, Simonetti AW, De Neuter B, Vanhamme L, Lemmerling P, et al. An automated quantitation of short echo time MRS spectra in an open source software environment: AQSES. NMR Biomed 2007;20:493-504. |
|21.||Niewenhuys R. The Human Central Nervous System: A Synopsis and Atlas. New York: Springer-Verlag; 1996. |
|22.||Jacobs MA, Horskα A, van Zijl PC, Barker PB. Quantitative proton MR spectroscopic imaging of normal human cerebellum and brain stem. Magn Reson Med 2001;46:699-705. |
|23.||Henriksen O. In vivo quantitation of metabolite concentrations in the brain by means of proton MRS. NMR Biomed 1995;8:139-48. |
|24.||Warntjes JB, Dahlqvist O, Lundberg P. Novel method for rapid, simultaneous T1, T*2, and proton density quantification. Magn Reson Med 2007;57:528-37. |
|25.||Petroff OA, Spencer DD, Alger JR, Prichard JW. High-field proton magnetic resonance spectroscopy of human cerebrum obtained during surgery for epilepsy. Neurology 1989;39:1197-202. |
|26.||Burri R, Bigler P, Straehl P, Posse S, Colombo JP, Herschkowitz N. Brain development: 1H magnetic resonance spectroscopy of rat brain extracts compared with chromatographic methods. Neurochem Res 1990;15:1009-16. |
|27.||Lowry OH, Berger SJ, Chi MM, Carter JG, Blackshaw A, Outlaw W. Diversity of metabolic patterns in human brain tumors-I: High energy phosphate compounds and basic composition. J Neurochem 1977;29:959-77. |
|28.||Nioka S, Chance B, Hilberman M, Subramanian HV, Leigh JS Jr, Veech RL, et al. Relationship between intracellular pH and energy metabolism in dog brain as measured by 31P-NMR. J Appl Physiol 1987;62:2094-102. |
|29.||Dross K, Kewitz H. Concentration and origin of choline in the rat brain. Naunyn Schmiedebergs Arch Pharmacol 1972;274:91-106. |
|30.||Stokes CE, Gillon KR, Hawthorne JN. Free and total lipid myo-inositol concentrations decrease with age in human brain. Biochim Biophys Acta 1983;753:136-8. |
|31.||Gasparovic C, Arfai N, Smid N, Feeney DM. Decrease and recovery of N-acetylaspartate/creatine in rat brain remote from focal injury. J Neurotrauma 2001;18:241-6. |
|32.||Minati L, Grisoli M, Bruzzone MG. MR spectroscopy, functional MRI, and diffusion-tensor imaging in the aging brain: A conceptual review. J Geriatr Psychiatry Neurol 2007;20:3-21. |
|33.||Choi CG, Frahm J. Localized proton MRS of the human hippocampus: Metabolite concentrations and relaxation times. Magn Reson Med 1999;41:204-7. |
|34.||Macrμ MA, Garreffa G, Giove F, Guardati M, Ambrosini A, Colonnese C, et al. In vivo quantitative 1H MRS of cerebellum and evaluation of quantitation reproducibility by simulation of different levels of noise and spectral resolution. Magn Reson Imaging 2004;22:1385-93. |
[Figure 1], [Figure 2], [Figure 3]
[Table 1], [Table 2], [Table 3]
|This article has been cited by|
||Brain Region Mapping Using Global Metabolomics
| ||Julijana Ivanisevic,Adrian A. Epstein,Michael E. Kurczy,Paul H. Benton,Winnie Uritboonthai,Howard S. Fox,Michael D. Boska,Howard E. Gendelman,Gary Siuzdak |
| ||Chemistry & Biology. 2014; 21(11): 1575 |
|[Pubmed] | [DOI]|
||Análisis de espectros de resonancia magnética
| ||Alfonso Lozano,Nathaly Barbosa,Luis Agulles |
| ||Revista Colombiana de Cancerología. 2013; 17(1): 11 |
|[Pubmed] | [DOI]|
||Motor system alterations in retired former athletes: the role of aging and concussion history
| ||Louis De Beaumont,Sébastien Tremblay,Luke C Henry,Judes Poirier,Maryse Lassonde,Hugo Théoret |
| ||BMC Neurology. 2013; 13(1): 109 |
|[Pubmed] | [DOI]|
||Rat brain myo-inositol 3-phosphate synthase is a phosphoprotein
| ||R. N. Parthasarathy,J. Lakshmanan,M. Thangavel,R. S. Seelan,J. I. Stagner,A. J. Janckila,R. E. Vadnal,M. F. Casanova,L. K. Parthasarathy |
| ||Molecular and Cellular Biochemistry. 2013; 378(1-2): 83 |
|[Pubmed] | [DOI]|
||Is the macromolecule signal tissue-specific in healthy human brain? A1H MRS study at 7 tesla in the occipital lobe
| ||Benoît Schaller,Lijing Xin,Rolf Gruetter |
| ||Magnetic Resonance in Medicine. 2013; : n/a |
|[Pubmed] | [DOI]|
||A Guide to the Metabolic Pathways and Function of Metabolites Observed in Human Brain 1H Magnetic Resonance Spectra
| ||Caroline D. Rae |
| ||Neurochemical Research. 2013; |
|[Pubmed] | [DOI]|
||Age-related differences in metabolites in the posterior cingulate cortex and hippocampus of normal ageing brain: A 1H-MRS study
| ||Harmen Reyngoudt,Tom Claeys,Leslie Vlerick,Stijn Verleden,Marjan Acou,Karel Deblaere,Yves De Deene,Kurt Audenaert,Ingeborg Goethals,Eric Achten |
| ||European Journal of Radiology. 2012; 81(3): e223 |
|[Pubmed] | [DOI]|
||Short- and long-term quantitation reproducibility of brain metabolites in the medial wall using proton echo planar spectroscopic imaging
| ||Shang-Yueh Tsai,Yi-Ru Lin,Woan-Chyi Wang,David M. Niddam |
| ||NeuroImage. 2012; 63(3): 1020 |
|[Pubmed] | [DOI]|
||Proton MR Spectroscopy in Metabolic Assessment of Musculoskeletal Lesions
| ||Ty K. Subhawong,Xin Wang,Daniel J. Durand,Michael A. Jacobs,John A. Carrino,Antonio J. Machado,Laura M. Fayad |
| ||American Journal of Roentgenology. 2012; 198(1): 162 |
|[Pubmed] | [DOI]|
||Age-related differences in metabolites in the posterior cingulate cortex and hippocampus of normal ageing brain: A 1H-MRS study
| ||Reyngoudt, H. and Claeys, T. and Vlerick, L. and Verleden, S. and Acou, M. and Deblaere, K. and De Deene, Y. and Audenaert, K. and Goethals, I. and Achten, E. |
| ||European Journal of Radiology. 2012; 81(3): e223-e231 |
||Proton MR spectroscopy in metabolic assessment of musculoskeletal lesions
| ||Subhawong, T.K. and Wang, X. and Durand, D.J. and Jacobs, M.A. and Carrino, J.A. and Machado, A.J. and Fayad, L.M. |
| ||American Journal of Roentgenology. 2012; 198(1): 162-172 |
||Short- and long-term quantitation reproducibility of brain metabolites in the medial wall using proton echo planar spectroscopic imaging
| ||Tsai, S.-Y. and Lin, Y.-R. and Wang, W.-C. and Niddam, D.M. |
| ||NeuroImage. 2012; 63(3): 1020-1029 |
||Applications of chemical shift imaging for AD
| ||Scott, C.J.M., Szilagyi, G.M., Chavez, S., Ganda, A., Black, S.E. |
| ||Current Medical Imaging Reviews. 2011; 7(1): 88-95 |
||In vivo mapping of brain myo-inositol
| ||Mohammad Haris,Kejia Cai,Anup Singh,Hari Hariharan,Ravinder Reddy |
| ||NeuroImage. 2011; 54(3): 2079 |
|[Pubmed] | [DOI]|