Skip to main content
  • Research article
  • Open access
  • Published:

The neurochemical basis of human cortical auditory processing: combining proton magnetic resonance spectroscopy and magnetoencephalography



A combination of magnetoencephalography and proton magnetic resonance spectroscopy was used to correlate the electrophysiology of rapid auditory processing and the neurochemistry of the auditory cortex in 15 healthy adults. To assess rapid auditory processing in the left auditory cortex, the amplitude and decrement of the N1m peak, the major component of the late auditory evoked response, were measured during rapidly successive presentation of acoustic stimuli. We tested the hypothesis that: (i) the amplitude of the N1m response and (ii) its decrement during rapid stimulation are associated with the cortical neurochemistry as determined by proton magnetic resonance spectroscopy.


Our results demonstrated a significant association between the concentrations of N-acetylaspartate, a marker of neuronal integrity, and the amplitudes of individual N1m responses. In addition, the concentrations of choline-containing compounds, representing the functional integrity of membranes, were significantly associated with N1m amplitudes. No significant association was found between the concentrations of the glutamate/glutamine pool and the amplitudes of the first N1m. No significant associations were seen between the decrement of the N1m (the relative amplitude of the second N1m peak) and the concentrations of N-acetylaspartate, choline-containing compounds, or the glutamate/glutamine pool. However, there was a trend for higher glutamate/glutamine concentrations in individuals with higher relative N1m amplitude.


These results suggest that neuronal and membrane functions are important for rapid auditory processing. This investigation provides a first link between the electrophysiology, as recorded by magnetoencephalography, and the neurochemistry, as assessed by proton magnetic resonance spectroscopy, of the auditory cortex.


Hearing is one of the fundamental and most important abilities of the vertebrate nervous system. The neurophysiology of the auditory system is increasingly well understood, however, little is known about the associations between the electrophysiology and the neurochemistry of the auditory cortex. One way of assessing auditory function is by measuring the amplitude and the decrement of the N1m wave, the major peak of the late auditory response, during rapidly successive presentation of acoustic stimuli. Recurrent acoustic stimulation is associated with characteristic short- and long-term decreases in the response of auditory neurons as seen by electrophysiologic recordings [13] and by functional magnetic resonance imaging [4, 5]. The decrement of auditory evoked responses with rapid stimulation, also known as habituation or sensory gating, is believed to represent cortical filtering of irrelevant input [6].

Recent studies have investigated the effects of various psychopharmacological agents on cortical auditory processing [7]. A single dose of haloperidol, a dopamin D2-receptor antagonist, affected the mismatch negativity, an event-related potential to deviating acoustic stimuli, but had no effect on the amplitude and latency of the N1, the electric counterpart of the N1m [8]. Administration of scopolamine, a centrally acting cholinergic antagonist, increased the latency of the N1m peak in young [9] and old [10] individuals. Alcohol consumption similarly increased the N1m latency and decreased the N1m amplitude [11]. Infusion of ketamine, an N-methyl-D-aspartate (NMDA) antagonist, resulted in significant increase of the N1 peak amplitude without affecting the N1 latency [12]. In addition, serotonergic neurotransmission plays an important role in the regulation of auditory processing [13].

In the present study, magnetoencephalography (MEG) and proton magnetic resonance spectroscopy (1H-MRS) were performed to assess the relationship between the electrophysiology and the neurochemistry of the left auditory cortex. As MEG allows the recording of cortical magnetic fields with supreme temporal resolution [14], this technique is ideally suited for the investigation of the exact timing of neural functions [15]. The recording of auditory evoked fields (AEF) has a high reliability (test-retest reproducibility) [14] and a high validity (consistency with intracranial recordings) [16]. 1H-MRS provides an opportunity to assess the regional neurochemistry of the brain in vivo [17]. We tested the hypotheses that: (i) the amplitude of the N1m wave, the major component of the long-latency auditory evoked responses and (ii) the decrement of this response during rapid stimulation are associated with brain metabolites measurable by 1H-MRS.



Analysis of variance (ANOVA) did not reveal significant differences between N1m baseline-to-peak amplitudes after stimulation with the vowel /a/ and N1m amplitudes after stimulation with a matched sine tone (Figure 1C and 1D, Tukey's multiple comparison test, p > 0.05 for the comparisons N1m1(vowel) vs. N1m1(tone), N1m2(vowel) vs. N1m2(tone), N1m3(vowel) vs. N1m3(tone) and N1m4(vowel) vs. N1m4(tone)). Thus, a grand average of all responses to vowels and tones was calculated for each individual subject. The combined analysis of both conditions yielded a better signal-to-noise ratio, resulting in a more robust source localisation and a more accurate calculation of the amplitude over time. All further analyses were performed with the grand average across individuals and conditions.

Figure 1
figure 1

N1m as a function of stimulus type. (A, B) Acoustic waveforms (upper graphs) and spectrograms (lower graphs) of the vowel (A) and sine tone stimulus (B). In this spectrogram, the acoustic energy in a certain frequency band and at a certain time point is proportional to the darkness of the grey shading. (C, D) Mean ± SD amplitudes of the four consecutive N1m peaks for the vowel (C) and the tone (D) conditions are presented (n = 15). The N1m peak amplitudes in the grand average of vowel and tone stimulation are displayed in Figure 2A.

A repeated-measures ANOVA showed that the N1m amplitude decreased with the number of responses (Figure 2A, stimulation with four successive stimuli in a train. F = 40.14, p < 0.0001). Post hoc tests using Tukey's multiple comparison test indicated that the amplitudes of the second, third and fourth N1m were significantly smaller than the amplitude of the first N1m (p < 0.001 for each comparison). The relative amplitudes of the second, third, and fourth N1m waves (each normalised to the peak amplitude of the first N1m) are depicted in Figure 2C.

Figure 2
figure 2

N1m amplitudes and latencies. (A, B) Mean amplitudes ± SD (A) and latencies (B) of the four N1m peaks (grand average of vowel and tone stimulation). * p < 0.05, * p < 0.01, *** p < 0.001. (C, D) Mean ± SD relative amplitudes (C) and relative latencies (D) of the four N1m peaks. Amplitudes and latencies are relative to the first N1m peak.

The N1m peak latency significantly increased with the number of responses (Figure 2B, Friedman test, a non-parametric repeated measures ANOVA, p < 0.0001). Post hoc tests using Dunn's multiple comparison test demonstrated that the latency of the fourth N1m was significantly longer than the latency of the first (p < 0.05) and the second N1m (p < 0.001). The relative latencies of the second, third, and fourth N1m waves are shown in Figure 2D.

Magnetic resonance spectroscopy

The concentrations of N-acetylaspartate (NAA) (median 12 institutional units (IU), range 6–20 IU), choline-containing compounds (Cho) (median 1.4 IU, range 0.9–3.0 IU) and the glutamine/glutamate pool (Glx) (median 22, range 4–40 IU) demonstrated considerable intersubject variability.

Associations between magnetic resonance spectroscopy and magnetoencephalography

The amplitude of the first N1m peak was significantly associated with the concentration of NAA (Figure 3, linear regression, R 2 = 0.369, F(1,13), p = 0.016) and Cho (Figure 4, R 2 = 0.287, F(1,13), p = 0.040). No significant association was found between Glx and the amplitude of the first N1m. Linear regression analysis revealed no significant associations between the relative amplitude of the second N1m response and NAA, Cho or Glx. However, there was a trend for higher Glx concentrations in individuals with higher relative N1m amplitude (R 2 = 0.26, F(1,13), p = 0.054) (Figure 5).

Figure 3
figure 3

N1m amplitude as a function of N-acetylaspartate. Linear regression plot of the amplitude of the first N1m vs. N-acetylaspartate (NAA) (R 2 = 0.369, F(1,13), p = 0.016). The 95% confidence intervals are represented by dashed lines. The insert shows the normal probability plot of the residuals, indicating that the errors are normally distributed. IU, institutional unit.

Figure 4
figure 4

N1m amplitude as a function of choline. Linear regression plot of the amplitude of the first N1m vs. choline-containing compounds (Cho) (R 2 = 0.287, F(1,13), p = 0.040). The 95% confidence intervals are represented by dashed lines. The insert shows the normal probability plot of the residuals, indicating that the errors are normally distributed. IU, institutional unit.

As this study included participants between 18 and 68 years, we addressed the question of whether age might influence the measured variables and the calculated statistical associations. No significant associations were found between age and the tested variables (NAA, Glx, Cho, first N1m amplitude, relative second N1m amplitude; data not shown). In addition, we reanalysed the data after the exclusion of the two oldest subjects (60 and 68 years old). For this subgroup (n = 13, mean age = 27 years), the amplitude of the first N1m peak was associated with the concentration of NAA (R 2 = 0.324, F(1,11), p = 0.042) and with the concentration of Cho (R 2 = 0.437, F(1,11), p = 0.014). However, no significant association was found between the relative amplitude of the second N1m peak and the concentration of Glx (R 2 = 0.259, F(1,11), p = 0.076).


Magnetoencephalographic and MR spectroscopic investigations of 15 adults revealed: (i) a significant association between the amplitudes of individual N1m responses and the concentrations of NAA and Cho and (ii) a trend for higher Glx concentrations in individuals with higher relative N1m amplitude. The present study is, to our knowledge, the first to combine magnetoencephalographic assessment of cortical auditory processing and MR spectroscopic estimation of neurochemicals in the auditory cortex. In an earlier study from our group, auditory MEG and MRS were used independently to investigate electrophysiological and neurochemical changes in patients with depression [18]. In the present investigation, in contrast, we examined the interactions between electrophysiological and neurochemical variables. Our data suggest that the interindividual differences of cortical auditory processing are generated, at least in part, by the variability of neurochemicals within the auditory cortex.

Auditory evoked magnetic fields

In the present study, no significant differences between the vowel and the sine-tone condition were found regarding the absolute or relative amplitudes of the N1m peak. This result corroborates a previous report in which no significant differences between the amplitudes of the N1m elicited by a sine tone or the vowel /a/ were seen in a passive listening condition (watching a silent movie, as in the present study) [19]. In contrast, vowels evoked a significantly stronger N1m response than tones in a study involving a stimulus detection task [20].

For the present study, we investigated associations between N1m amplitudes, but not N1m latencies, and the neurochemistry of the auditory cortex. We did not assess the influence of MRS variables on N1m latency because this latency depends, in part, on the subcortical auditory pathways, which could not be investigated with 1H-MRS, used here.

Variability of auditory evoked responses

The macroscopic and microscopic variability of primary and secondary auditory cortices has been carefully investigated [21, 22]. In particular, the anatomy of the major generator site of the human N1m response, the planum temporale [23], is characterised by extensive interindividual variability [22]. On the individual level, our AEF recordings displayed a pronounced interindividual variability of the peak amplitude and decrement (Figure 2A and 2C). This variability was expected in the light of previous experiments on auditory evoked potentials [24, 25].

Associations between N1m amplitude and neurochemicals

The amino-acid NAA is found almost exclusively in neurons and axons, and can be used as an indirect measure of neuronal integrity and density. Loss of NAA has been reported in a variety of disorders, e.g., in individuals with epilepsy [26], Alzheimer's disease [27] or cognitive impairment [28]. Cho, the choline-containing compounds measurable with 1H-MRS, are dominated by phosphocholine, glycerophosphocholine, and free choline [29]. As precursors or degradation products of membrane phospholipids make a major contribution to the Cho signal, Cho is regarded as a membrane marker [30]. In addition, choline is the rate-limiting substrate for the synthesis of the neurotransmitter acetylcholine via the enzyme choline acetyltransferase.

The magnetic dipole moment detected by MEG represents the postsynaptic currents of approximately 1 million synapses activated synchronously. Differences in brain magnetic fields are thus assumed to reflect differences in the number of activated synapses and/or in the coherence of activated neurons [31]. Our data suggest that the number and integrity of neurons (NAA) and the density and the functional integrity of cell membranes (Cho) are associated with the amplitude of auditory evoked responses.

The metabolic variables reported here, however, explain only in part the interindividual differences in cortical auditory processing. The association between NAA and the N1m amplitude accounts for 37% of the measured variance, while the association between Cho and the N1m amplitude accounts for 29% of the variance.

Associations between N1m decrement and neurochemicals

This study revealed no significant associations between the concentrations of NAA, Cho or Glx and the reduction of the N1m amplitude after a preceding stimulus (the N1m decrement). However, there was a trend for higher Glx concentrations in individuals with higher relative N1m amplitude (i.e., smaller decrement) (Figure 5, R 2 = 0.26, F(1,13), p = 0.054).

Figure 5
figure 5

Relative N1m amplitude as a function of glutamine. Linear regression plot of the relative amplitude of the second N1m peak vs. the glutamine/glutamate pool (Glx) (R 2 = 0.26, F(1,13), p = 0.054). The 95% confidence intervals are represented by dashed lines. The insert shows the normal probability plot of the residuals, indicating that the errors are normally distributed. IU, institutional unit.

Previous work on cortical processing of rapidly recurring stimuli indicated that inhibitory interneuronal networks modulate the amplitude and the decrement of cortical evoked responses [3, 32]. Several classes of interneurons with distinct functional properties regulate the inhibitory activity in the neocortex [33]. The activity of interneurons appears to be mainly regulated by the major inhibitory neurotransmitter gamma-aminobutyric acid (GABA) [34] and the major excitatory neurotransmitter, glutamate [35].

The concentration of GABA could not be measured by the MRS technique used here. The spectroscopic measurement of GABA in vivo is difficult, due to the overlapping peaks of NAA, Glx, and creatine. Although MR protocols have been developed for the assessment of GABA, these techniques require a measurement volume of more than 35 cm3 (compared with 3.375 cm3 here), much larger than the auditory cortex [25].

Glx represents the concentrations of the amino-acid glutamate and its precursor glutamine. In the peripheral and central auditory system, glutamate receptors of the alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) type play an important role in the rapid synaptic transmission [36]. Moreover, AMPA receptors regulate the processing of rapidly successive stimuli in the peripheral auditory system [37].

Methodological considerations

The aim of the present study was to correlate electrophysiological parameters of auditory processing with the neurochemistry of the auditory cortex. This study relies on the following assumptions.

Reliability of MEG and MRS results

As the MEG and MRS measurements could not be performed simultaneously, the present study relies on the reproducibility of the measured electrophysiologic and neurochemical variables over time. In several investigations using EEG and MEG, a high reproducibility of the N1/N1m amplitude has been found [38, 39]. Using the experimental paradigm of this study, replication measurements with two participants were performed. The amplitude of the first N1m response and the response decrement from the first to the second N1m were highly consistent within 3 months. This result corroborates a study by Lütkenhöner et al, who used the same MEG system as used for the present study [14]. Similarly, the intraindividual reproducibility of NAA, Cho and Glx concentrations was demonstrated by our group [40] and by others [41, 42].

Stability of MRS results under acoustic stimulation

MRS examinations, similar to conventional MRI scans, are associated with intense scanner noise caused by rapidly switching the gradient coils on and off. During the 1H-MRS scans here, short (100 ms) periods of scanner noise, repeated every 1.5 seconds, were produced. As the acoustic stimulus during the MEG recording was (necessarily) different, our study relies on the stability of neurochemicals under different kinds of acoustic stimulation. In a functional MRS study comparing scanner noise alone or in combination with siren noise or music, the concentration of lactate changed in relation to acoustic stimulation, while the other spectroscopic parameters remained stable [43]. Similarly, visual stimulation [44] or a silent word-generation task increased brain lactate concentration, but did not change the levels of the other metabolites as measurable by 1H-MRS. Thus, it appears unlikely that the noise-induced activation of the auditory cortex during MRS significantly changed the metabolites measured here.

Directions for future research

Since the advent of non-invasive techniques to investigate human brain function, such as positron emission tomography, functional magnetic resonance imaging, MRS and MEG, numerous studies have been performed using either modality. Recently, an increasing number of researchers have combined neuroimaging modalities in an attempt to overcome the limitations of each of these methods [45, 46].

The present study correlated the electrophysiology of auditory processing as recorded by MEG with the neurochemistry of the auditory cortex as assessed by MRS. Our study, however, has several limitations. First, the small number of participants limited the power of the statistical analyses. In addition, we acquired data from only the left hemisphere. The 37-channel MEG system used here could be positioned close to the superior temporal plane, resulting in a good signal-to-noise ratio, but could only record the neural activity of one hemisphere at a time. The processing of acoustic stimuli can differ between hemispheres, especially when speech stimuli are used [47]. Similarly, the concentrations of neurochemicals, as assessed by MRS, may display interhemispheric differences [48]. Therefore, we cannot rule out that MEG-MRS associations for the right auditory cortex differ from those reported here for the left auditory cortex. For future studies, a larger sample size and data collection from both hemispheres (using a whole-head MEG system) are desirable. Moreover, we performed MEG recordings while participants were watching a silent movie (i.e., in a passive listening condition). The associations between MEG and MRS parameters reported here might change when MEG is recorded under top-down attentional modulation, as attention is known to influence auditory processing [49].

In addition to the concentration of NAA and Cho, other factors have been shown to influence auditory cortical processing, such as serotonergic neurotransmission (see Background) [7]. In a study combining electrophysiology and genetics, the amplitude increase of the N1/P2 component in response to increasing stimulus intensities (termed the loudness dependence) was significantly different between individuals with different variants of the serotonin transporter gene. We expect that genetic studies will considerably contribute to future research on individual auditory processing.


The combined use of MEG and MRS indicated that the processing of rapid auditory stimuli seems to depend, among other factors, on the functional integrity of neurons and cell membranes in the auditory cortex, as reflected by the concentrations of NAA and Cho. This approach provides evidence that individual differences in the neurochemical composition of the cortex account, in part, for individual electrophysiological differences in auditory processing.



Fifteen healthy volunteers (six women, nine men) with a median age of 24 years (range: 18–68 years) participated in this study. Two subjects were 60 years or older, the remaining 13 were 54 years or younger. Fourteen participants were right-handed, one was left-handed [50]. All participants had a normal audiological status and were without a history of neurological or otological disorders. All individuals gave their informed consent to participate in the study. This project was reviewed and approved by the Research Ethics Board, Medical Faculty, University of Münster. The participants of this study had served as healthy controls in a previous study on metabolic and electrophysiological changes in patients with depression [18]. MEG recordings of subgroups of the participants in this study were also published separately [51, 52].

Magnetic resonance spectroscopy

To investigate individual brain anatomy and to exclude possible brain lesions, structural MR images were obtained from every subject. T1-weighted, three-dimensional, spoiled gradient echo MRI of the whole brain and T2- and proton density-weighted fast spin echo images in axial and coronal orientation were performed on a 1.5-Tesla scanner (Magnetom SP, Siemens, Erlangen, Germany). To assess the concentrations of brain metabolites, single voxel stimulated echo acquisition mode (STEAM) spectroscopy was used (echo time, TE = 20 ms; repetition time, TR = 2.5 s; number of scans = 128) [53]. Using the individual structural MRI, the MRS voxel was centred at the left transverse gyrus of Heschl (Brodmann area 41) on the dorsal surface of the superior temporal plane containing the primary and secondary human auditory cortex. The MRS voxel volume was 1.5 × 1.5 × 1.5 cm3 (Figure 6).

Figure 6
figure 6

MR spectroscopy of the auditory cortex. (A) Anatomical transverse MR image of the left auditory cortex with the superimposed borders of the MRS voxel (1.5 × 1.5 × 1.5 cm3). The image shows the anterior borders of the voxel, which extends 1.5 cm in the posterior direction. (B) 1H-MRS spectrum from the auditory cortex of a single subject. NAA denotes N-acetylaspartate; Glx, glutamate/glutamine; Cr, creatine/phosphocreatine; Cho, choline-containing compounds. The original and the fitted spectrum (using the linear combination model, black line) are displayed.

MRS postprocessing

The postprocessing of MRS spectra was performed as described previously [40, 54, 55]. Concentrations of N-acetylaspartate (NAA), choline-containing compounds (Cho) and the glutamine/glutamate pool (Glx) were quantified using a time domain fitting program, based on the linear combination (LC) model [56]. Owing to the overlapping resonances of glutamate, glutamine and GABA at a field strength of 1.5 T, a combined fitting of the glutamine/glutamate pool was performed. Metabolite concentrations, acquired by the LC model and scaled to water, were corrected for coil loading [57] and partial volume effects [58]. To account for variations in voxel composition (cerebrospinal fluid, grey matter and white matter), brain segmentation was performed with a semiautomatic interactive algorithm (M. Fiebich, University of Applied Sciences, Giessen, Germany) [54]. Metabolites were then normalised to the grey-matter fraction in the voxels, and expressed in institutional units (IU). For robustness of results, only data with a fitting error <20% of the standard deviation was included in the final analysis.

Magnetoencephalographic measurement

Auditory stimulation

Short sequences of rapidly recurring speech and non-speech sounds served as stimuli. The German vowel /a/ (duration 260 ms, fundamental frequency f0 234 Hz), spoken by a female speech-language pathologist, and a matched sine tone were chosen (frequency 234 Hz, duration 260 ms). The envelope of the sine tone was adjusted to the envelope of the vowel /a/. Sequences of four successive stimuli (either four vowels or four tones) were presented with an onset-to-onset interstimulus interval of 450 ms (Figure 7). As electromagnetic activation of the auditory cortex lasts for approximately 400 ms after the onset of a single transient stimulus, an interstimulus interval of 450 ms (or 2.2 stimuli/second) was chosen to avoid overlap between successive brain responses. The onset-to-onset interval between sequences was 4.5 seconds (randomised between 4 and 5 seconds). This interval is long enough to allow a substantial recovery of the N1m component before the onset of the following train [59]. Participants listened passively to 160 trials of the vowel stimulation and to 160 trials of the tone stimulation in randomised order.

Figure 7
figure 7

Auditory evoked magnetic fields. (A) Magnetic waveforms of an individual subject, superimposition of 37 MEG channels, covering the left hemisphere. (B) Amplitude of the four N1m peaks in the same subject. The grey bars indicate the onset and duration of the acoustic stimulation.

MEG data acquisition

AEFs were recorded with a 37-channel axial gradiometer system (Magnes I; BTi, San Diego, USA) in a magnetically shielded room and sampled at rate of 512.4 Hz. The participants were in a right lateral position with the body supported by a vacuum cushion to minimise head and body movements during the measurement. The sensor array was positioned over the auditory cortex as closely as possible to the subject's head. To ensure a stable passive listening condition, subjects watched a self-selected silent video that attracted their attention. Participants were instructed not to move their head, to stay awake, and to keep their eyes open. Immediately before each MEG measurement, the individual hearing thresholds were determined for the vowel and the tone stimuli separately. The stimuli were delivered to a silicon earpiece in the right ear via speakers outside the shielded room and a plastic tube of 6.3 m length. The stimulation system was able to transmit frequencies up to approximately 4500 Hz [60]. All stimuli were presented with an intensity of 60 dB above the individual hearing threshold.

MEG data analysis

Artefact-contaminated epochs were excluded if the magnetic field value exceeded the baseline value (calculated in the time window of -250 to 0 ms before stimulus onset) by 2 pT. Between 4% and 13% of the data had to be discarded. After exclusion of artefacts, three averaged datasets were created per subject: (i) vowels alone, (ii) sine tones alone, and (iii) vowels and sine tones combined. All datasets were baseline corrected (-250 to 0 ms before the first stimulus) and filtered with a band-pass filter of 0.01–40 Hz. Figure 7A shows the overlay of the recordings of all 37 channels in an individual participant. To identify the successive N1m peaks, the root mean square of amplitudes was calculated for a time window of 0 to 2000 ms, based on the amplitudes in each of the 37 MEG channels. To obtain a dipole that was representative of the N1m, we first calculated a single equivalent current dipole for each sampling point and then averaged the dipole parameters for a 30-ms time window around the activity peak (response to the first stimulus of series, as determined by root mean square of amplitude curves). This N1m dipole was finally used to calculate the dipole moment over the entire epoch. Figure 7B displays the dipole moment as a function of time in an individual participant. For this calculation, the location and the direction of the N1m dipole were assumed to be constant (fixed-dipole approach). The strongest deflections in the latency range of the N1m were identified as N1m responses. Based on the dipole moment over time, the peak amplitudes and peak latencies of the N1m responses were determined. In most participants, a baseline shift was detected from the beginning of the first response to the fourth response. To ensure accurate quantification of amplitudes, dipole moments were measured relative to a baseline that was defined as the mean value before the onset of each stimulus (-50 to 0 ms prestimulus). The amplitude ratio of the second and first N1m (termed the relative amplitude of the second N1m) was then calculated.

Assessment of test-retest reliability

To test for reliability of MRS data, we measured six healthy individuals twice, at a mean ± SD of 5.2 ± 0.8 months apart, in a separate study [40]. Repeated-measures ANOVA yielded no significant difference between Glx concentrations of the left dorsolateral prefrontal cortex in the first and second measurement (F = 0.39, p = 0.56). Glx concentrations are usually responsible for the largest variance in MRS data (B. Pfleiderer, unpublished observation). The coefficient of variation for Glx in our previous study [40] was 12%. These values are similar to the coefficients of variation, ranging from 13% to 15%, which were previously reported for 1H chemical shift imaging experiments [61]. These results indicate that the metabolite concentrations measured here are reasonably stable over time.

To assess the reliability of AEF recordings, measurements were repeated in two participants after 4 months (Figure 8). The overall correlation between amplitudes over time in the first and the second recording were high (participant 1: r = 0.93, participant 2: r = 0.87). Most important for this study, the first N1m responses were almost identical. The reliability of later responses, in particular the second N1m, was smaller. These measurements corroborate an earlier study on the reliability of N1m recordings [14].

Figure 8
figure 8

Reliability of auditory evoked field recordings. (A, B) Amplitudes of two AEF measurements, 4 months apart, in two different subjects.

Statistical analysis

To test if the collected data were normally distributed, the Shapiro-Wilk normality test was performed [62]. This test is widely used to test for normality in smaller samples because of its high statistical power compared with alternative tests [63]. The assumption of a normal distribution was rejected if p < 0.01. According to this test, all electrophysiological and neurochemical variables were normally distributed, except the latencies of the four N1m responses. To compare the absolute and relative amplitudes of the four successive N1m waves, a repeated-measures ANOVA with Tukey's multiple-comparison post hoc test was calculated (Figure 2A). A nonparametric repeated-measures ANOVA (Friedman test with Dunn's multiple comparison post hoc test) was performed to compare the absolute and relative latencies of the four N1m waves (Figure 2B). To assess the relationship between absolute and relative N1m amplitudes and neurochemical parameters, a linear regression analysis was performed (Figures 3, 4, 5). Linear-regression models assume that the unpredicted variation (the error term ε) is normally distributed. To test for the normality of errors, a normal-probability plot of the residuals was performed (insets in Figures 3, 4, 5). If the error distribution is normal, the points on the normal-probability plot should fall close to a diagonal line. The linear-regression model was used to fit a regression line through the data and to calculate the 95% confidence intervals (Figures 3, 4, 5). Statistical analyses were performed with the statistical package R for Mac OS X [64].


  1. Fruhstorfer H, Soveri P, Järvilehto T: Short-term habituation of the auditory evoked response in man. Electroencephalogr Clin Neurophysiol. 1970, 28 (2): 153-161. 10.1016/0013-4694(70)90183-5.

    Article  CAS  PubMed  Google Scholar 

  2. Hari R, Kaila K, Katila T, Tuomisto T, Varpula T: Interstimulus interval dependence of the auditory vertex response and its magnetic counterpart: implications for their neural generation. Electroencephalogr Clin Neurophysiol. 1982, 54 (5): 561-569. 10.1016/0013-4694(82)90041-4.

    Article  CAS  PubMed  Google Scholar 

  3. Okamoto H, Kakigi R, Gunji A, Kubo T, Pantev C: The dependence of the auditory evoked N1m decrement on the bandwidth of preceding notch-filtered noise. Eur J Neurosci. 2005, 21 (7): 1957-1961.

    Article  CAS  PubMed  Google Scholar 

  4. Michael N, Ostermann J, Sörös P, Schwindt W, Pfleiderer B: Altered habituation in the auditory cortex in a subgroup of depressed patients by functional magnetic resonance imaging. Neuropsychobiology. 2004, 49: 5-9. 10.1159/000075331.

    Article  PubMed  Google Scholar 

  5. Pfleiderer B, Ostermann J, Michael N, Heindel W: Visualization of auditory habituation by fMRI. Neuroimage. 2002, 17 (4): 1705-1710. 10.1006/nimg.2002.1308.

    Article  CAS  PubMed  Google Scholar 

  6. Grunwald T, Boutros N, Pezer N, von Oertzen J, Fernández G, Schaller C, Elger C: Neuronal substrates of sensory gating within the human brain. Biol Psychiatry. 2003, 53 (6): 511-519. 10.1016/S0006-3223(02)01673-6.

    Article  PubMed  Google Scholar 

  7. Kähkönen S: Magnetoencephalography (MEG): a non-invasive tool for studying cortical effects in psychopharmacology. Int J Neuropsychopharmacol. 2006, 9 (3): 367-372. 10.1017/S1461145705005894.

    Article  PubMed  Google Scholar 

  8. Kähkönen S, Ahveninen J, Jääskeläinen I, Kaakkola S, Näätänen R, Huttunen J, Pekkonen E: Effects of haloperidol on selective attention: a combined whole-head MEG and high-resolution EEG study. Neuropsychopharmacology. 2001, 25 (4): 498-504. 10.1016/S0893-133X(01)00255-X.

    Article  PubMed  Google Scholar 

  9. Pekkonen E, Hirvonen J, Jääskeläinen I, Kaakkola S, Huttunen J: Auditory sensory memory and the cholinergic system: implications for Alzheimer's disease. Neuroimage. 2001, 14 (2): 376-382. 10.1006/nimg.2001.0805.

    Article  CAS  PubMed  Google Scholar 

  10. Pekkonen E, Jääskeläinen I, Kaakkola S, Ahveninen J: Cholinergic modulation of preattentive auditory processing in aging. Neuroimage. 2005, 27 (2): 387-392. 10.1016/j.neuroimage.2005.04.018.

    Article  PubMed  Google Scholar 

  11. Kähkönen S, Marttinen Rossi E, Yamashita H: Alcohol impairs auditory processing of frequency changes and novel sounds: a combined MEG and EEG study. Psychopharmacology (Berl). 2005, 177 (4): 366-372. 10.1007/s00213-004-1960-1.

    Article  Google Scholar 

  12. Umbricht D, Schmid L, Koller R, Vollenweider F, Hell D, Javitt D: Ketamine-induced deficits in auditory and visual context-dependent processing in healthy volunteers: implications for models of cognitive deficits in schizophrenia. Arch Gen Psychiatry. 2000, 57 (12): 1139-1147. 10.1001/archpsyc.57.12.1139.

    Article  CAS  PubMed  Google Scholar 

  13. Manjarrez G, Hernandez E, Robles A, Hernandez J: N1/P2 component of auditory evoked potential reflect changes of the brain serotonin biosynthesis in rats. Nutr Neurosci. 2005, 8 (4): 213-218. 10.1080/10284150500170971.

    Article  CAS  PubMed  Google Scholar 

  14. Lütkenhöner B: Dipole source localization by means of maximum likelihood estimation. II. Experimental evaluation. Electroencephalogr Clin Neurophysiol. 1998, 106 (4): 322-329. 10.1016/S0013-4694(97)00139-9.

    Article  PubMed  Google Scholar 

  15. Hari R, Levänen S, Raij T: Timing of human cortical functions during cognition: role of MEG. Trends Cogn Sci. 2000, 4 (12): 455-462. 10.1016/S1364-6613(00)01549-7.

    Article  PubMed  Google Scholar 

  16. Godey B, Schwartz D, deGraaf J, Chauvel P, Liégeois-Chauvel C: Neuromagnetic source localization of auditory evoked fields and intracerebral evoked potentials: a comparison of data in the same patients. Clin Neurophysiol. 2001, 112 (10): 1850-1859. 10.1016/S1388-2457(01)00636-8.

    Article  CAS  PubMed  Google Scholar 

  17. Rothman D, Behar K, Hyder F, Shulman R: In vivo NMR studies of the glutamate neurotransmitter flux and neuroenergetics: implications for brain function. Annu Rev Physiol. 2003, 65: 401-427. 10.1146/annurev.physiol.65.092101.142131.

    Article  CAS  PubMed  Google Scholar 

  18. Tollkötter M, Pfleiderer B, Sörös P, Michael N: Effects of antidepressive therapy on auditory processing in severely depressed patients: A combined MRS and MEG study. J Psychiatr Res. 2006, 40 (4): 293-306. 10.1016/j.jpsychires.2005.09.003.

    Article  PubMed  Google Scholar 

  19. Tiitinen H, Sivonen P, Alku P, Virtanen J, Näätänen R: Electromagnetic recordings reveal latency differences in speech and tone processing in humans. Brain Res Cogn Brain Res. 1999, 8 (3): 355-363. 10.1016/S0926-6410(99)00028-2.

    Article  CAS  PubMed  Google Scholar 

  20. Gootjes L, Raij T, Salmelin R, Hari R: Left-hemisphere dominance for processing of vowels: a whole-scalp neuromagnetic study. Neuroreport. 1999, 10 (14): 2987-2991.

    Article  CAS  PubMed  Google Scholar 

  21. Rademacher J, Morosan P, Schormann T, Schleicher A, Werner C, Freund H, Zilles K: Probabilistic mapping and volume measurement of human primary auditory cortex. Neuroimage. 2001, 13 (4): 669-683. 10.1006/nimg.2000.0714.

    Article  CAS  PubMed  Google Scholar 

  22. Westbury C, Zatorre R, Evans A: Quantifying variability in the planum temporale: a probability map. Cereb Cortex. 1999, 9 (4): 392-405. 10.1093/cercor/9.4.392.

    Article  CAS  PubMed  Google Scholar 

  23. Lütkenhöner B, Steinsträter O: High-precision neuromagnetic study of the functional organization of the human auditory cortex. Audiol Neurootol. 1998, 3 (2–3): 191-213.

    PubMed  Google Scholar 

  24. Lauter J, Karzon R: Individual differences in auditory electric responses: comparisons of between-subject and within-subject variability. V. Amplitude-variability comparisons in early, middle, and late responses. Scand Audiol. 1990, 19 (4): 201-206.

    Article  CAS  PubMed  Google Scholar 

  25. Lütkenhöner B, Krumbholz K, Seither-Preisler A: Studies of tonotopy based on wave N100 of the auditory evoked field are problematic. Neuroimage. 2003, 19 (3): 935-949. 10.1016/S1053-8119(03)00172-1.

    Article  PubMed  Google Scholar 

  26. Mueller S, Suhy J, Laxer K, Flenniken D, Axelrad J, Capizzano A, Weiner M: Reduced extrahippocampal NAA in mesial temporal lobe epilepsy. Epilepsia. 2002, 43 (10): 1210-1216. 10.1046/j.1528-1157.2002.10802.x.

    Article  PubMed Central  PubMed  Google Scholar 

  27. Schuff N, Capizzano A, Du A, Amend D, O'Neill J, Norman D, Kramer J, Jagust W, Miller B, Wolkowitz O, Yaffe K, Weiner M: Selective reduction of N-acetylaspartate in medial temporal and parietal lobes in AD. Neurology. 2002, 58 (6): 928-935.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Chao L, Schuff N, Kramer J, Du A, Capizzano A, O'Neill J, Wolkowitz O, Jagust W, Chui H, Miller B, Yaffe K, Weiner M: Reduced medial temporal lobe N-acetylaspartate in cognitively impaired but nondemented patients. Neurology. 2005, 64 (2): 282-289.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Miller B, Chang L, Booth R, Ernst T, Cornford M, Nikas D, McBride D, Jenden D: In vivo 1H MRS choline: correlation with in vitro chemistry/histology. Life Sci. 1996, 58 (22): 1929-1935. 10.1016/0024-3205(96)00182-8.

    Article  CAS  PubMed  Google Scholar 

  30. Boulanger Y, Labelle M, Khiat A: Role of phospholipase A(2) on the variations of the choline signal intensity observed by 1H magnetic resonance spectroscopy in brain diseases. Brain Res Brain Res Rev. 2000, 33 (2–3): 380-389. 10.1016/S0165-0173(00)00037-0.

    Article  CAS  PubMed  Google Scholar 

  31. Lü Z, Williamson S: Spatial extent of coherent sensory-evoked cortical activity. Exp Brain Res. 1991, 84 (2): 411-416. 10.1007/BF00231463.

    Article  PubMed  Google Scholar 

  32. Miller C, Freedman R: The activity of hippocampal interneurons and pyramidal cells during the response of the hippocampus to repeated auditory stimuli. Neuroscience. 1995, 69 (2): 371-381. 10.1016/0306-4522(95)00249-I.

    Article  CAS  PubMed  Google Scholar 

  33. McLean M, Busza A, Wald L, Simister R, Barker G, Williams S: In vivo GABA+ measurement at 1.5T using a PRESS-localized double quantum filter. Magn Reson Med. 2002, 48 (2): 233-241. 10.1002/mrm.10208.

    Article  CAS  PubMed  Google Scholar 

  34. Blatow M, Caputi A, Monyer H: Molecular diversity of neocortical GABAergic interneurones. J Physiol. 2005, 562 (Pt 1): 99-105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Virtanen J, Ahveninen J, Ilmoniemi R, Näätänen R, Pekkonen E: Replicability of MEG and EEG measures of the auditory N1/N1m-response. Electroencephalogr Clin Neurophysiol. 1998, 108 (3): 291-298. 10.1016/S0168-5597(98)00006-9.

    Article  CAS  PubMed  Google Scholar 

  36. Parks T: The AMPA receptors of auditory neurons. Hear Res. 2000, 147 (1–2): 77-91. 10.1016/S0378-5955(00)00122-2.

    Article  CAS  PubMed  Google Scholar 

  37. Duan M, Canlon B: Short-term adaptation in the peripheral auditory system is related to the AMPA receptor. Acta Otolaryngol. 2001, 121 (1): 21-27. 10.1080/000164801300006227.

    Article  CAS  PubMed  Google Scholar 

  38. Rotondo E, Bruschetta G, Saccà A, Bramanti P, Di Pasquale M: Straightforward relative quantitation and age-related human standards of N-acetylaspartate at the centrum semiovale level by CSI (1)H-MRS. Magn Reson Imaging. 2003, 21 (9): 1055-1060. 10.1016/S0730-725X(03)00211-X.

    Article  CAS  PubMed  Google Scholar 

  39. Tedeschi G, Bertolino A, Campbell G, Barnett A, Duyn J, Jacob P, Moonen C, Alger J, Di Chiro G: Reproducibility of proton MR spectroscopic imaging findings. AJNR Am J Neuroradiol. 1996, 17 (10): 1871-1879.

    CAS  PubMed  Google Scholar 

  40. Michael N, Gösling M, Reutemann M, Kersting A, Heindel W, Arolt V, Pfleiderer B: Metabolic changes after repetitive transcranial magnetic stimulation (rTMS) of the left prefrontal cortex: a sham-controlled proton magnetic resonance spectroscopy (1H MRS) study of healthy brain. Eur J Neurosci. 2003, 17 (11): 2462-2468. 10.1046/j.1460-9568.2003.02683.x.

    Article  PubMed  Google Scholar 

  41. Prichard J, Rothman D, Novotny E, Petroff O, Kuwabara T, Avison M, Howseman A, Hanstock C, Shulman R: Lactate rise detected by 1H NMR in human visual cortex during physiologic stimulation. Proc Natl Acad Sci USA. 1991, 88 (13): 5829-5831. 10.1073/pnas.88.13.5829.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Richards T, Gates G, Gardner J, Merrill T, Hayes C, Panagiotides H, Serafini S, Rubel E: Functional MR spectroscopy of the auditory cortex in healthy subjects and patients with sudden hearing loss. AJNR Am J Neuroradiol. 1997, 18 (4): 611-620.

    CAS  PubMed  Google Scholar 

  43. Urrila A, Hakkarainen A, Heikkinen S, Vuori K, Stenberg D, Häkkinen A, Lundbom N, Porkka-Heiskanen T: Metabolic imaging of human cognition: an fMRI/1H-MRS study of brain lactate response to silent word generation. J Cereb Blood Flow Metab. 2003, 23 (8): 942-948. 10.1097/01.WCB.0000080652.64357.1D.

    Article  CAS  PubMed  Google Scholar 

  44. Zaehle T, Wüstenberg T, Meyer M, Jäncke L: Evidence for rapid auditory perception as the foundation of speech processing: a sparse temporal sampling fMRI study. Eur J Neurosci. 2004, 20 (9): 2447-2456. 10.1111/j.1460-9568.2004.03687.x.

    Article  CAS  PubMed  Google Scholar 

  45. Kamada K, Houkin K, Takeuchi F, Ishii N, Ikeda J, Sawamura Y, Kuriki S, Kawaguchi H, Iwasaki Y: Visualization of the eloquent motor system by integration of MEG, functional, and anisotropic diffusion-weighted MRI in functional neuronavigation. Surg Neurol. 2003, 59 (5): 352-361. 10.1016/S0090-3019(03)00018-1.

    Article  PubMed  Google Scholar 

  46. Mathiak K, Fallgatter A: Combining magnetoencephalography and functional magnetic resonance imaging. Int Rev Neurobiol. 2005, 68: 121-148.

    Article  PubMed  Google Scholar 

  47. Tiitinen H, Mäkelä A, Mäkinen V, May P, Alku P: Disentangling the effects of phonation and articulation: hemispheric asymmetries in the auditory N1m response of the human brain. BMC Neurosci. 2005, 6: 62-10.1186/1471-2202-6-62.

    Article  PubMed Central  PubMed  Google Scholar 

  48. Jayasundar R, Raghunathan P: Evidence for left-right asymmetries in the proton MRS of brain in normal volunteers. Magn Reson Imaging. 1997, 15 (2): 223-234. 10.1016/S0730-725X(96)00342-6.

    Article  CAS  PubMed  Google Scholar 

  49. Obleser J, Elbert T, Eulitz C: Attentional influences on functional mapping of speech sounds in human auditory cortex. BMC Neurosci. 2004, 5: 24-10.1186/1471-2202-5-24.

    Article  PubMed Central  PubMed  Google Scholar 

  50. Oldfield R: The assessment and analysis of handedness: the Edinburgh inventory. Neuropsychologia. 1971, 9 (1): 97-113. 10.1016/0028-3932(71)90067-4.

    Article  CAS  PubMed  Google Scholar 

  51. Teismann I, Sörös P, Manemann E, Ross B, Pantev C, Knecht S: Responsiveness to repeated speech stimuli persists in left but not right auditory cortex. Neuroreport. 2004, 15 (8): 1267-1270.

    Article  PubMed  Google Scholar 

  52. Sörös P, Dziewas R, Manemann E, Teismann IK, Lütkenhöner B: No indication of brain reorganization after unilateral ischemic lesions of the auditory cortex. Neurology. 2006,

    Google Scholar 

  53. Haase A, Frahm J, Matthaei D, Hänicke W, Bomsdorf H, Kunz D, Tischler R: MR imaging using stimulated echoes (STEAM). Radiology. 1986, 160 (3): 787-790.

    Article  CAS  PubMed  Google Scholar 

  54. Pfleiderer B, Michael N, Gösling M, Erfurth A, Wolgast M, Heindel W: Altered glutamate/glutamine levels within the left dorsolateral prefrontal cortex in aute mania – impact of tissue segmentation and related correlations. Proeceedings of the International Society for Magnetic Resonance in Medicine, 11th Scientific Meeting, Toronto, Canada. 2003

    Google Scholar 

  55. Pfleiderer B, Ohrmann P, Suslow T, Wolgast M, Gerlach A, Heindel W, Michael N: N-acetylaspartate levels of left frontal cortex are associated with verbal intelligence in women but not in men: a proton magnetic resonance spectroscopy study. Neuroscience. 2004, 123 (4): 1053-1058. 10.1016/j.neuroscience.2003.11.008.

    Article  CAS  PubMed  Google Scholar 

  56. Provencher S: Automatic quantitation of localized in vivo 1H spectra with LCModel. NMR Biomed. 2001, 14 (4): 260-264. 10.1002/nbm.698.

    Article  CAS  PubMed  Google Scholar 

  57. Michaelis T, Merboldt K, 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 (1): 219-227.

    Article  CAS  PubMed  Google Scholar 

  58. Weber-Fahr W, Ende G, Braus D, Bachert P, Soher B, Henn F, Büchel C: A fully automated method for tissue segmentation and CSF-correction of proton MRSI metabolites corroborates abnormal hippocampal NAA in schizophrenia. Neuroimage. 2002, 16 (1): 49-60. 10.1006/nimg.2002.1057.

    Article  CAS  PubMed  Google Scholar 

  59. Woods D, Elmasian R: The habituation of event-related potentials to speech sounds and tones. Electroencephalogr Clin Neurophysiol. 1986, 65 (6): 447-459. 10.1016/0168-5597(86)90024-9.

    Article  CAS  PubMed  Google Scholar 

  60. Lütkenhöner B, Seither-Preisler A, Seither S: Piano tones evoke stronger magnetic fields than pure tones or noise, both in musicians and non-musicians. Neuroimage. 2006, 30 (3): 927-937. 10.1016/j.neuroimage.2005.10.034.

    Article  PubMed  Google Scholar 

  61. Chard D, McLean M, Parker G, MacManus D, Miller D: Reproducibility of in vivo metabolite quantification with proton magnetic resonance spectroscopic imaging. J Magn Reson Imaging. 2002, 15 (2): 219-225. 10.1002/jmri.10043.

    Article  PubMed  Google Scholar 

  62. Shapiro S, Wilk M: An analysis of variance test for normality (complete samples). Biometrika. 1965, 52: 591-611. 10.2307/2333709.

    Article  Google Scholar 

  63. Shapiro S, Wilk M, Chen HJ: A comparative study of various tests of normality. J Am Statistical Assoc. 1968, 63: 1343-1372. 10.2307/2285889.

    Article  Google Scholar 

  64. R Development Core Team: R: A language and environment for statistical computing. 2005, Vienna, Austria: R Foundation for Statistical Computing, []

    Google Scholar 

Download references


The authors thank Elisabeth Manemann and Inga K. Teismann who performed MEG recordings in a subgroup of participants. Peter Sörös is funded by the Heart and Stroke Foundation of Ontario Centre for Stroke Recovery.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Peter Sörös.

Additional information

Authors' contributions

PS and NM contributed equally to this work. PS, NM and BP designed the study and drafted the manuscript. PS and MT recorded and analysed MEG data. BP and MT performed and analysed MRS scans. All authors read and approved the final manuscript.

Peter Sörös, Nikolaus Michael contributed equally to this work.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Sörös, P., Michael, N., Tollkötter, M. et al. The neurochemical basis of human cortical auditory processing: combining proton magnetic resonance spectroscopy and magnetoencephalography . BMC Biol 4, 25 (2006).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: