Evidence of functional connectivity between auditory cortical areas revealed by amplitude modulation sound processing.

The human auditory cortex includes several interconnected areas. A better understanding of the mechanisms involved in auditory cortical functions requires a detailed knowledge of neuronal connectivity between functional cortical regions. In human, it is difficult to track in vivo neuronal connectivity. We investigated the inter-area connection in vivo in the auditory cortex using a method of directed coherence (DCOH) applied to depth auditory evoked potentials (AEPs). This paper presents simultaneous AEPs recordings from insular gyrus, primary and secondary cortices (Heschl’s gyrus and planum temporale) and associative areas (BA 22) with multilead intracerebral electrodes in response to sinusoidal modulated white noises in four epileptic patients who underwent invasive monitoring with depth electrodes for epilepsy surgery. DCOH allowed estimation of the causality between two signals recorded from different cortical sites. The results showed: (i) a predominant auditory stream within the primary auditory cortex (PAC) from the most medial region to the most lateral one whatever the modulation frequency (MF), (ii) unidirectional functional connection from the primary to secondary auditory cortex, (iii) a major auditory propagation from the posterior areas to the anterior ones, particularly at 8, 16 and 32 Hz, (iv) a particular role of Heschl’s sulcus (HS), dispatching information to the different auditory areas. These findings suggest that cortical processing of auditory information is performed in serial and parallel streams. Our data showed that the auditory propagation could not be associated to an unidirectional traveling wave but to a constant interaction between these areas, that could reflect the large adaptive and plastic capacities of auditory cortex. The role of the insular gyrus is discussed.


INTRODUCTION
The human auditory cortex includes several interconnected areas.The complexity of these interacting structures presents a challenge to researchers to determine how this system analyzes and uses auditory information.Indeed, auditory perception comprises the discrimination of a vast array of sound features such as pitch, rhythm, timbre, loudness, or amplitude modulations.
Neurophysiologically speaking, the discrimination of two successive sounds differing in one parameter reflects the discrimination of two different neural representations evoked in the same or a different group of neurons.Perception of stimulus change is strongly correlated with the amplitude, latency and other parameters of surface as well as depth-recorded evoked potentials (EPs) (for review, Eggermont and Ponton, 2002).Depth EPs, only detected if a sufficient number of neurons are activated synchronously, have allowed physiological distinction of some of the auditory areas defined by cytoarchitectonic criteria (Braak, 1980;Galaburda and Sanides, 1980;Liégeois-Chauvel et al., 1991, 1994;Howard et al., 2000).Eloquent auditory cortical areas have been found within the posterior two-thirds of the superior temporal gyrus including Heschl's gyrus (HG), planum temporale (PT), the most posterior part of superior temporal gyrus (STG) (Post T1), insular gyrus (IG) and parietal operculum (Liégeois-Chauvel et al., 2004).Little is known about their functional organization and their connectivity, and most current beliefs are derived from work on animals in whom physiological and functional data could be interpreted within an anatomical framework (Rauschecker et al. 1997(Rauschecker et al. , 1998;;Kaas andHackett, 1998, 2000).The primate auditory cortex is subdivided into three areas: the core, belt, and parabelt.Most work has shown that auditory information is distributed from the core areas to surrounding belts that relay information to parabelt regions (for review, Kaas et al., 1999;Kaas and Hackett, 2000).
To our knowledge, only a few studies exist on the connectivity between human auditory areas.Most of these have been performed on post-mortem human brains by use of tracer injections (Rivier and

Patients
Four patients (two males, two females, 25-35 years of age), among the 20 patients who participated in the study reported in (Liégeois-Chauvel et al., 2004), participated in the present study.They suffered from drug-resistant partial epilepsy and were implanted with chronic SEEG electrodes in the right (2 patients) or in the left auditory cortex (2 patients).Several additional electrodes were implanted in various cortical structures in order to: (i) determine which structures were involved in the initiation and propagation of seizures, and (ii) accurately delineate the limits of future cortical excision (Talairach et al., 1974).The choice of the anatomical location of electrodes was based on clinical and video-EEG recordings and on magnetic resonance imaging (MRI), and was made independently of the present study.Physiology of the auditory cortex has been studied and considered as normal.Neuropsychological assessment indicated that the patients showed typical language representation (left hemisphere).Recordings of brainstem evoked potentials and pure-tone audiograms carried out before SEEG indicated intact cochlear and brainstem auditory functions.This study did not add any invasive procedure to the depth EEG recordings performed routinely in the neurological evaluation.The patients were informed about the research protocol during SEEG and gave their fully informed consent for participating in this study.

Anatomical definition of depth electrode position
The stereotactic method was based on the co-registration of the patient's MRI with the stereotactic angiogram.This was performed in order to avoid any injury of brain vessels.Moreover, this method allowed orthogonal introduction of multilead electrodes (0.8 mm diameter, 10 or 15 electrode contacts of 2 mm length each with 1.5 mm spacing between contacts) in the stereotactic space (Szikla et al., 1977;Talairach and Tournoux, 1988).The anatomical position of each contact was then identified on the basis of (i) an axial scanner image acquired before the removal of electrodes, and (ii) an MRI scan performed after the removal of electrodes (Liégeois-Chauvel et al., 1991).
Anatomical location and responses elicited by pure tones and amplitude modulated tones from these electrodes have been previously reported, and for legibility in describing these electrodes we have attributed the same case number and indicated this in superscript (Liégeois-Chauvel et al., 2004).
One of the four patients (Case 3) had three electrodes recording from four different auditory regions giving an exceptional and unique sampling.The first is electrode P C3 , whose medial contacts (1-4) recorded activity from the medial part of the right HG (i.e.PAC), whereas the lateral contacts (5-9) recorded activity from PT.The first contacts (1-2) of electrode H C3 (second electrode) were located in the insular gyrus, contacts 3-5 were located in PAC, contact H6 was in the transverse Heschl's sulcus (HS) and contacts 7-9 were in PT.The third electrode is T C3 : contacts 4 and 5 were located in the lateral part of HG, which is in the secondary auditory cortex (SAC), whereas the lateral contacts (6-9) were in BA 22.The three remaining patients had only two electrodes in the auditory cortex.In Case 12, three areas are recorded: the left PAC (H8-H12 C12 ), the left SAC (T12-T15 C12 ) and the left posterior part of Superior temporal gyrus (Post T1 H C12 ).In Case 1, four contacts of electrode H C1 (H6-H9) were in the right PAC, H10 C1 was in the Transverse HS.Eight contacts were located in right postero-lateral part of the superior temporal gyrus (Post T1 H C1 and Post T1 P C1 ).Case 20 had two electrodes, which explored the left PAC (H10-H15 C20 and T10-T15 C20 ).

Stimuli
All stimuli were generated using a 16-bit D/A converter at a sampling frequency of 44.1 kHz.The stimuli were white noises modulated sinusoidally in amplitude at frequencies MF = 4,8,16,32,64 and 128 Hz, with a 100 % depth modulation.The starting phase of the modulation was fixed at 270° (thus, each stimulus started at an amplitude minimum in the modulation waveform).All AM stimuli were 1 s in duration and were shaped by rising and falling 25-ms cosine ramps/damps.They were equated in rms (root mean square) and presented binaurally via Sennheiser headphones at 75 dB SPL (rms).Series of 100 mixed stimuli (4/32 Hz, 8/64 Hz, 16/128 Hz,) were delivered to the listener in random order.

Recordings
The recordings of intracerebral AEPs were monopolar, with each contact of a given depth electrode referenced to an extra-dural lead.All signals were amplified and bandpass filtered between 0.15 and 200 Hz (more precisely, a highpass filter (cutoff frequency: 0.15 Hz, rolloff: 12 dB/oct) and a lowpass filter (cutoff frequency: 200 Hz, rolloff: 24 dB/oct) were used to bandpass filter signals).
Data acquisition started 164 ms before the presentation of the sound and lasted for 1476 ms.53 trials were averaged for each sound.During each recording session, the patient lays comfortably in a chair in a sound-attenuated room and listened passively to the sounds.
The averaged signals presented a frequency peak at the modulation frequency in their power spectral density.The oscillations at the modulation frequency were also visible in the temporal shape of the signal in the time window between 200 to 1000 ms post-stimulus onset.

Method
Based on a bivariate autoregressive (AR) model and the Granger causality, the directed coherence (DCOH) has been introduced in (Saito and Harashima, 1981).Compared to the ordinary coherence defined in (Bendat and Piersol, 1986), the directed coherence provides the direction of propagation between two time series in addition to their degree of linear interaction.
The principle of this method, whose complete description is given in Appendix A (cf. Supplementary Material), is to model two signals X and Y by a bivariate AR model.From this model, two directed coherences DCOH yx (f) and DCOH xy (f), representing the degree of linear interaction from X to Y and from Y to X respectively, are obtained at any frequency f.In this modeling, different sources are introduced, W x (respectively W y ) is the source from which signal X (respectively Y) is generated, and W s is a common source to X and Y. Sources W x , W y , W s are gaussian white noises with a null mean and a variance of one.The directed coherences are normalized in the range [0-1] and represent the ratio between the different contributions of sources to one signal.
A simple example is given in equation (1), in which X is autoregressive and Y is computed thanks to X: where k is the current sample.
In this equation, the samples of Y depend on the preceding sample of X, while the samples of X do as observed in the previous study (Liégeois-Chauvel et al., 2004).

Threshold determination for directed coherence significance
To study the connectivity between auditory cortical areas revealed by AM sound processing, we

Anatomical distribution of AM responses within the different auditory areas
The magnitude of the evoked responses ( 53 III for Case 3 and in Table A1 for Case 20 at MF = 8 Hz in Supplementary Material.These tables follow the reading convention of source (column) to target (row).In Tables III and A1 the highlighted values correspond to DCOH modules upper than 0.6.

***Insert Table III about here***
In Table III, we observe that the auditory flow tends mainly to propagate serially from the medial leads (lead P1 C3 ) to the most lateral ones (lead P4 C3 ).This direction of propagation is found whatever the AM frequency within the posterior part of PAC (electrode P C3 ).In the anterior part of PAC (electrode H C3 ), this unidirectionality seems less clear since it exists for H3 C3 whose activity is propagated to H4 C3 and H5 C3 , but less for H4 C3 .Moreover, some variations are seen according to the AM sound.Interestingly, the oscillatory activity was propagated from the posterior to anterior part of PAC (i.e. from electrode P C3 to H C3 ) except for 4 Hz.
The values in Tables III and A1 for MF = 8 Hz are reported on Fig. 4a and Fig. 4b for Cases 3 and 20 respectively.The arrow width codes interaction strength (0 ≤ no arrow ≤ 0.6 < increasing width solid lines in step of 0.2).For reasons of clarity in Fig. 4b, only the directed coherences modules superior to 0.8 are represented.
***Insert Figure 4 about here*** Fig. 4a shows that the activity propagates from the posterior part to the anterior part of PAC of Case 3, mainly to H4 and H5.In Fig. 4b dealing with Case 20, the location of the electrode H C20 is in the same region as H C3 and the electrode T C20 is more anterior and lateral than H C20 .There is a strong tendency of medio-lateral propagation but which is not strictly serial, for instance auditory stream mainly goes from H10 C20 and H1 C20 to H14 C20 and H15 C20 and from T12 C20 to T15 C20 skipping T13 C20 and T14 C20 .In the anterior-lateral part of PAC C20 there is a constant interaction inside this area, predominantly from T C20 to H C20 .

Between different anatomical areas
Firstly, we will describe the results obtained from Case 3 because the exceptional sampling of the recording sites within the auditory cortex (see Methods) allowed recording of activity in four different anatomical areas of the right auditory cortex (PAC, SAC, PT, BA 22) and in the insular gyrus.In the same manner that we distinguished the posterior part (PAC P: P1-4) from the anterior part of PAC (PAC H: H3-5), we separate evoked activity recorded from posterior PT P (P5 to P9) to anterior PT H (H7-9).We studied the functional connectivity between these seven sub-regions.
The directed coherence algorithm was then applied to each couple constituted of AEPs relative to leads located in the different sub-regions.The corresponding AEPs are taken in the time window between 200 to 1000 ms post-stimulus onset.The means of all the DCOH modules over all the frequencies MF = 4, 8, 16 and 32 Hz from the different couplings were summarized in Table A2 (see Supplementary Material) and displayed in Fig. 5  In Fig. 6b and Fig. 6c, auditory information propagates from primary and secondary auditory areas (PAC and SAC) to the posterior part of STG (Post T1).As previously observed, HS sends information to PAC but it receives information from Post T1 in Case 12.
Fig. 7a shows the information propagation at MF = 16 Hz in Case 3. As in the previous case, PAC P sends information to PT P, PAC H, PT H, SAC T and A22 T. Mainly, information flows from posterior to anterior zones and from medial to lateral zones.We observe the striking influence of HS, which directs information to the PAC, SAC and BA 22 and even to the Insula.
Interestingly, we notice that the connectivity is frequency dependent as shown in Fig. 7b for 32 Hz.
One can readily see that the interaction strengths are very high, particularly from PAC P, Insula H

DISCUSSION
Overall, the present data indicate: (i) a medio-lateral propagation within the PAC (ii) a predominant auditory stream from primary to secondary and associative auditory areas, and (iii) a hierarchical and parallel processing between the cortical areas.Finally, the data showed (iv) the special role played by the HS in the dispatching of auditory stream.These results are discussed below.
The human auditory cortex is currently subdivided into several architectonically defined areas.The number of these areas as well as their extent and exact position vary between authors.Using cytoarchitectonic criteria, Brodmann (Brodmann 1909) distinguished three areas, called 41, 42, and 22, and von Economo and Koskinas (von Economo and Koskinas 1925) four areas, as did Braak (Braak 1978).Combining cyto-and myelo-architecture, Galaburda and Sanides (Galaburda and Sanides 1980) identified in the same region eight areas.

Primary auditory cortex
As reported by Rademacher et al. (Rademacher et al. 2001), the precise location and absolute size of PAC cannot be reliably inferred from the macroanatomic landmarks and could comprise between 16 and 92% of the cortical volume of HG according to subjects.Evidence for identification of PAC came from electrophysiological recordings reported in vivo studies (Célesia, 1976;Liégeois-Chauvel et al., 1991).Comparison between cytoarchitectonic studies and our data suggests that areas PAC P C3 , PAC H C3,C20 and PAC T C20 would correspond to A1 (Rivier and Clarke, 1997;Clarke and Rivier, 1998) or Te1 (Morosan et al., 2001;Rademacher et al., 2001).The striking result is the co-existence of simplicity and complexity within PAC connectivity.In one hand, we observed a medio-lateral auditory propagation observed along electrodes located in different parts of right and left PAC suggesting a hierarchical processing.On the other hand, the differentiation of physiological regions showed distinct interactions between themselves.Interestingly these areas were not differentiated by AEPs, but by their connectivities.The first two recorded from P C3 and H C3 showed an unidirectional propagation of auditory stream from right posterior to anterior PAC.
This postero-anterior propagation did not seem to continue to occur in the antero-lateral part of PAC, at least in the left hemisphere.We observed that the electrode located in the antero-lateral part (T12-15 C20 ) sent information to H C20 , which is located medially and posteriorly to T C20 .Although we do not have data on cytoarchitectony, we may be able to derive some information on a putative distinction between the three primary areas from the connectivity patterns.The first one is located in the tip of HG recorded by P C3 , the second one is located in the medio-anterior part (H C3 and H C20 ) and the third one in the antero-lateral of PAC (T C20 ).This physiological distinction could be related to the subdivisions (Te1.1,Te1.0 and Te2) of PAC reported by Morosan et al. (Morosan et al. 2001) and is close to the organization of auditory core of monkey (Kaas and Hackett, 2000).Our data showed that the auditory propagation could not be associated to an unidirectional traveling wave but to a constant interaction between these areas.That could reflect the large adaptive and plastic capacities of auditory cortex (Engineer et al., 2004;Nelken, 2004).Fritz et al. (Fritz et al. 2003) have recently demonstrated rapid plastic changes in ferret auditory cortex within minutes of the beginning of a tone detection task.Such changes were smaller or absent when the animal was listening passively to same sounds.In human, Schneider et al. (Schneider et al. 2002) demonstrated correlation among the amplitude of N19/P30 primary auditory components, the gray matter volume of PAC and musical aptitude in groups of musicians versus non musicians.

The role of HS
The morphology of AEPs in response to tones in HS was quite different from those recorded from the PAC, the PT or Post T1 that may result in part from differences in thalamic input to these fields.
Interestingly, HS dispatches auditory stream to several regions such as Insula, PAC, SAC, PT, BA 22.It receives information from Post T1.Its place between upstream auditory areas (PAC, SAC) and downstream areas (PT, Post T1) (Rivier and Clarke, 1997) leads to speculation regarding its role.Auditory scene analysis needs a continuous processing of sounds that are spectrally and temporally complex.HS would represent an interface between processing in other auditory cortical regions and could play a modulator role since it sends auditory upstream to PAC only for 8 and 16 Hz.

Role of Insula
Strong connectivity is seen between the posterior part of the insular gyrus (PIC) and SAC, PT and BA 22 especially at highest AM frequency (32 Hz).A functional magnetic resonance imaging (f-MRI) study (Giraud et al., 2000) also investigated the cortical representation of AM using a set of white noises modulated sinusoidally in amplitude at frequencies covering a broad range (4-256 Hz).
Consistent with animal studies, this study revealed that, in humans, (i) temporal resolution degrades from the brainstem to the auditory cortex, and (ii) only a few cortical regions were tuned to high AM frequencies, for instance PIC is known to be an eloquent auditory region and seems to be strongly modulated by high AM frequency rather than low AM.Its functional role remains unclear.
Firstly, there were discrepancies in the interpretation of functional and anatomical data.As a matter of fact, PAC and PIC had a cytochrome oxidase profile that was compatible with a primary sensory area (Rivier and Clarke, 1997), whereas electrophysiological responses recorded from PIC were clearly different and could not be considered as primary responses (Liégeois-Chauvel et al., 2004).

Hierarchical and Parallel Processing within Human Auditory Areas?
Given that, similarly to the visual system (for review, Felleman and Van Essen, 1991), auditory areas might be organized in a hierarchical system (Rauschecker et al., 1997), we may be able to derive some information on a putative hierarchical order from the DCOH values.An auditory stream from the PAC to the secondary cortex (lateral part of HG and PT) and to BA 22 is observed.
The hierarchical organization within the human auditory cortex proposed by previous studies (Mesulam and Geula, 1994;Hutsler and Gazzaniga, 1996) seemed to be confirmed here.However, our results suggest an organization with a certain degree of parallel processing as observed in a recent study using magnetoencephalography (Inui et al., 2006).These authors found six different cortical areas as a function of latency onset of cortical activity and inferred the existence of a serial model of auditory processing along the medio-lateral axis as well as several parallel streams running from the PAC to the belt region and then the Posterior STG and PT and anteriorly from PAC to BA 22.We also observed a direct propagation from the PAC to BA 22 and to Post T1.The projections from PAC to BA 22 have been shown in neurophysiological studies of the temporal auditory fields in the chimpanzee (Bailey et al., 1943) and rhesus monkey (Ward et al., 1946;Sugar et al., 1948).Later, electrophysiological (Bignall, 1969) and anatomical (Hackett et al., 1998;Kaas andHackett, 1998, 2000) results have confirmed these auditory corticocortical connections in the nonhuman primate.The existence of a projection of HG upon Post T1 has been observed with direct cortical stimulation (Howard et al., 2000;Brugge et al., 2003).While applying electrical stimulation to sites in HG, Howard et al. (Howard et al. 2000) recorded short-latency evoked responses in posterior lateral superior temporal area (PLST) which corresponds to the area that we called Post T1.Moreover, Brugge et al. (Brugge et al. 2003) showed widespread convergence and divergence of input from HG to posterior STG with evidence for a reciprocal functional projection, from posterolateral STG to HG.By our study, we specified that this reciprocal projection was done via HS.
We were unable, however, to ascertain whether these were direct connections between two areas or an indirect one, the latter being consistent with the hierarchical connectivity model of Kaas and Hackett (Kaas and Hackett 1998) based on extensive electrophysiological mapping and neuroanatomic tracer studies in the monkey.

Is directed coherence an accurate tool?
The objective of using directed coherence is to find the propagation direction between two given leads located in auditory cortical areas.The problem may be viewed in terms of delay.Various methods based on linear systems have been proposed.Among non-parametric methods, the most classical corresponds to the maximum of the cross-correlation function.Some variants introducing a weighting window in the frequency domain were derived, such as the maximum likelihood.In the field of parametric methods, besides methods based on the fact that an observation is the sum of one signal derived from another observation and an additive noise, a vectorial autoregressive model was proposed in (Hannan and Thomson, 1981) from which a delay estimation is derived using the crosspower spectral density matrix (Hamon and Hannan, 1974).The first non-linear method, the average amount of mutual information based on mutual entropy, was proposed in (Mars and Van Arragon, 1982).In (Moddemeijer, 1989), the estimation of the probability law is carried out using the histogram method leading to the information theoretical delay estimator.Methods based on higher order moments (Nikias and Raghuveer, 1987) containing phase information are increasingly studied (Emile and Comon, 1995).Parametric approaches have also been proposed in the non-linear field such as the method based on the non-linear regression coefficient (Pijn, 1989).
Contrary to these methods that we have briefly summarised, and other methods of directed coherence such as directed transfer function (DTF) (Kaminski and Blinowska, 1991) or partial directed coherence (PDC) (Baccala and Sameshima, 2001), directed coherence (DCOH) uses a different approach.It takes into account the influence of the rest of the brain on the two investigated signals, in modeling the rest of the brain by a common source influencing these two signals, and thus seems more appropriate to study neuronal connections.This characteristic of the DCOH method, although it does not give the delay, provides the propagation information directly from the relations between samples.The DCOH has already been applied to EEG signals recorded by scalp electrodes in healthy subjects (Jing and Takigawa, 2000) and in epileptic subjects (Takigawa et al.,21 HAL author manuscript inserm-00188495, version 1 1996), or by intracranial electrodes in rats (Wang et al., 2002).Auditory event-related potentials have also been investigated by DCOH (Jing et al., 2001a(Jing et al., , 2001b)), but to our knowledge this paper is the first to study intracranial auditory event-related potentials.

Conclusion
The aim of the present work was to study the connectivity which exists between human auditory areas using the method of directed coherence (DCOH), applied to depth auditory evoked potentials.
DCOH allowed the estimation of causality between two signals recorded from different cortical sites.The results showed: (i) a predominant auditory stream within the primary auditory cortex (PAC) from the most medial region to the most lateral one whatever the MF, (ii) unidirectional functional connection from the primary to secondary auditory cortex, (iii) a major auditory propagation from the posterior areas to the anterior ones, particularly at 8, 16 and 32 Hz, (iv) a particular role of Heschl's sulcus (HS), dispatching information to the different auditory areas.
These findings suggest that cortical processing of auditory information is performed in serial and parallel streams.Table I.DCOH modules at the frequency MF = 8 Hz and the corresponding 95% confidence intervals between all the leads located in the PAC of Case 3 for the background EEG activity.Table II.DCOH modules at the frequency MF = 8 Hz and the corresponding 95% confidence intervals determined by bootstrap between all the leads located in the PAC of Case 3 for the auditory activity.HAL author manuscript inserm-00188495, version 1

TABLES
Fig. 1 (top)    shows typical synchronized AEPs recorded from third lead of electrode H C3 (lead H3 located in the right PAC) in response to a 16 Hz-MF stimulus.Periodic activity characterized by phasic responses close to the onset of each modulation cycle is clearly visible in this AEP preceded by initial evoked response (N1/P2).The power spectral density computed in the 200-1000 ms range is shown on Fig.
not depend on the samples of Y. Then X is autogenerated and generates Y.The resulting modules of DCOH yx (f) and DCOH xy (f) are given in Figure 2, for the frequency f in the range [0-500 Hz], with a sampling frequency F s = 1 kHz.***Insert Figure 2 about here*** In Figure 2, we can see that |DCOH yx (f)| (left) is close to 1, and |DCOH xy (f)| (right) is null, which shows the influence of X on Y whatever the frequency.This figure confirms the model given in equation (1).Our aim is to study the propagation between the auditory areas of the oscillatory part of eventrelated potentials, which oscillates at the modulation frequency of the stimuli.Consequently signals were examined in their quasi-stationary oscillatory part (from 200 to 1000 ms post-stimulus onset) and the DCOH method was applied for each pair of leads at each modulation frequency.We particularly computed, for each modulation frequency MF, the modules of DCOH from Y to X and from X to Y taken at the frequency f = MF, denoted |DCOH xy (MF)| and |DCOH yx (MF)| respectively.For the study of DCOH in Case 3, who had the largest sampling of electrodes and auditory regions explored, four different stimuli (amplitude modulated white noise at MF = 4, 8, 16 or 32 Hz) were used and the corresponding AR model orders were 140, 100, 70 and 50 respectively.For the three other patients, only stimulus at MF = 8 Hz (AR model order = 100) was used because it corresponded to the most representative modulation transfer function (MTF) across auditory regions must distinguish the activity induced by AM stimuli from the background EEG activity by determining a threshold for directed coherence.To determine this threshold, we compare the values obtained on EEG signals with those obtained on responses to stimuli.This preliminary study is carried out at 8 Hz for the leads of PAC of Case 3. In this purpose, we sent 53 stimuli at 8 Hz, which provided us with 53 responses to stimuli (53 trials) and 53 interstimuli signals (background EEG activity without auditory activity) for each recorded lead.We first applied the directed coherence algorithm to the 53 interstimuli signals recorded on each lead of the PAC.A standard deviation is then computed and a confidence interval derived.Secondly, we applied the directed coherence algorithm to the 53 responses to stimuli using a bootstrap technique: we built an AEP by averaging a set of 53 trials randomly chosen with replacement and we computed the DCOH measure on this AEP.This step is iterated one hundred times and a confidence interval is derived.Results of the propagation expressed by DCOH modules at frequency f = 8 Hz and the 95% confidence intervals of each DCOH measure are reported in trials) was analyzed as a function of the stimulus modulation frequency(Liégeois-Chauvel et al., 2004).The magnitude values of the spectral components were determined by considering the Fast Fourier Transforms of the responses over 1200 ms.These values are averaged and Fig.3shows the amplitude of the spectral peaks at the modulation frequency within main auditory areas (PAC, SAC, BA 22 and Post T1) as a function of AM frequency (4, 8, 16 and 32 Hz) and patients.***Insert Figure 3 about here*** As reported in the previous study, the distribution of MF varies from one contact to another and predominant response of cortical auditory areas to the low AM frequencies (4-8 Hz) is observed in each region.Best synchronization to 16 Hz could be, however, observed as on the three contacts H13-H15 in the PAC of Case 20.The response to 32 Hz AM is greater within the PAC than in the SAC or BA 22 in Case 3. Propagation within the primary auditory cortex (PAC) To study the propagation between the different leads located within the PAC, we distinguished the posterior part from the anterior part of PAC in Cases 3 and 20.The corresponding AEPs are taken in the time window between 200 to 1000 ms post-stimulus onset.Results of the propagation expressed by DCOH modules at each frequency f = MF are reported in Table for Case 3. ***Insert Figure 5 about here***This graph shows the major propagations of the oscillations between the seven auditory sub-regions of Case 3 independently of the AM frequency.The main source of information is the PAC P, which directs auditory information to the PAC H, confirming the posterior-to-anterior influence within the primary auditory cortex, and to the SAC T showing a primary-to-secondary auditory areas connectivity.In the same way, the PAC H is linked to the PT H (medial-to-lateral connection) and to the SAC T. This figure presents examples of strong unidirectional influence of primary auditory areas (PAC P and PAC H) on secondary ones (SAC T and A22 T) and on planum temporale (PT P and PT H).The insular gyrus (Insula H) also seems to send information to other areas (PT H, SAC T and A22 T) and receives information only from PAC P.For Case 3, the modules of DCOH calculated for 8, 16 and 32 Hz are given in Tables A3, A4 and A5, presented as Supplementary Material.We considered the Heschl's sulcus (lead H6) as a single area.Modules of DCOH ofCases 3, 1 and 12 calculated for 8 Hz are summarized on the graphs shown on Fig. 6 a, b, c respectively.Only modules of directed coherences superior to 0.6 are displayed.***Insert Figure 6 about here***In Fig.6a, the main connections are similar to those given in Fig.5, which correspond to the principal directed connections present between the 8 zones whatever the modulation frequency.But one can readily see that the HS dispatches information to several zones including PAC P, PT P, Insula H, PT H, SAC T and A22 T. The other main source of information is the PAC P which propagates information to the PT P, Insula H, PAC H, PT H, SAC T and A22 T. In general, information flows from posterior to anterior zones and from medial to lateral zones.
Fig. 8 summarizes all the results of Case 3 represented in terms of its anatomy.Two main streams

LEGENDSFigure 1 .
Figure 1.Above: Shape of the AEP recorded on lead H3 of Case 3 in response to a 16 Hz-MF

Figure 4a .
Figure 4a.Graph of propagation between the leads of electrode P and of electrode H located in the PAC, at the frequency MF = 8 Hz for Case 3. The arrow width codes interaction strength (0 ≤ no arrow ≤ 0.6 < increasing width solid lines in step of 0.2).The most medial lead is the first one and

Figure 4b .
Figure 4b.Graph of propagation between the leads of electrode T and of electrode H located in the PAC, at the frequency MF = 8 Hz for Case 20.The arrow width codes interaction strength (0 ≤ no arrow ≤ 0.8 < increasing width solid lines in step of 0.2).

Figure 5 .
Figure 5. Graph of propagation between 7 auditory sub-regions of Case 3, averaged over all the modulation frequencies MF = 4, 8, 16 and 32 Hz.The arrow width codes interaction strength (0 ≤ no arrow ≤ 0.6 < increasing width solid lines in step of 0.2).

Figure 6a .
Figure 6a.Graph of propagation between the 8 sub-regions of Case 3 at the frequency MF = 8 Hz.The arrow width codes interaction strength (0 ≤ no arrow ≤ 0.6 < increasing width solid lines in step

Figure 6b .
Figure 6b.Graph of propagation between the 3 auditory sub-regions of Case 1 at the frequency MF = 8 Hz.The arrow width codes interaction strength (0 ≤ no arrow ≤ 0.6 < increasing width solid

Figure 6c .
Figure 6c.Graph of propagation between the 4 auditory sub-regions of Case 12 at the frequency MF = 8 Hz.The arrow width codes interaction strength (0 ≤ no arrow ≤ 0.6 < increasing width solid

Figure 7a .
Figure 7a.Graph of propagation between the 8 auditory sub-regions of Case 3 at the frequency MF = 16 Hz.The arrow width codes interaction strength (0 ≤ no arrow ≤ 0.6 < increasing width solid

Figure 8 .
Figure 8. Cortico-cortical connections within the auditory areas of Case 3. Arrows indicate the Table I for the background EEG activity and in TableIIfor auditory activity.These tables follow the reading convention of source (column) to target (row).The DCOH obtained on background EEG activity is lower than 0.6 with a low variance.We therefore consider that DCOH measures higher than 0.6 are revealing of flow information for auditory activity.