F.-Xavier Alario
email: francois-xavier.alario@univ-amu.fr

A.-Sophie Dubarry

Anais Llorens

Agnes Trébuchon

Romain Carron

Catherine Liégeois-Chauvel

Christian-G Bénar

[ Aq

Parallel Processing in a Language Task research-article2017

Keywords: mental chronometry, cascaded processing, word production, intracranial EEG, single trial

We provide a quantitative assessment of the parallel-processing hypothesis included in various language-processing models. First, we highlight the importance of reasoning about cognitive processing at the level of single trials rather than using averages. Then, we report the results of an experiment in which the hypothesis was tested at an unprecedented level of granularity with intracerebral data recorded during a picture-naming task. We extracted patterns of significant high-gamma activity from multiple patients and combined them into a single analysis framework that identified consistent patterns. Average signals from different brain regions, presumably indexing distinct cognitive processes, revealed a large degree of concurrent activity. In comparison, at the level of single trials, the temporal overlap of detected significant activity was unexpectedly low, with the exception of activity in sensory cortices. Our novel methodology reveals some limits on the degree to which word production involves parallel processing.

Research Article

The stages involved in mental information processing [START_REF] Taylor | Stage analysis of reaction time[END_REF]) can be distinguished on the basis of the types of knowledge they involve (e.g., semantic, lexical, phonological) or by the timing with which they operate. Regarding timing, a primary distinction has been made between sequential and parallel architectures. Donders's (1868[START_REF] Donders | On the speed of mental processes (W. G. Koster[END_REF]) influential hypothesis and its variants postulate that, in certain tasks, cognitive stages operate serially, that is, without temporal overlap. An alternative view postulates parallel-that is, concurrent-processing across cognitive stages [START_REF] Mcclelland | On the time relations of mental processes: An examination of systems of processes in cascade[END_REF]. As a case in point, word production is thought to involve multiple types of representations and levels of processing. Their processing architecture has been extensively debated on the basis of various types of evidence [START_REF] Hickok | Computational neuroanatomy of speech production[END_REF][START_REF] Levelt | Spoken word production: A theory of lexical access[END_REF][START_REF] Munding | On the cortical dynamics of word production: A review of the MEG evidence[END_REF][START_REF] Rapp | Discreteness and interactivity in spoken word production[END_REF]. Most recently, intracranial brain activity has been marshaled to address this issue, and this work has led to contrasting conclusions [START_REF] Edwards | Spatiotemporal imaging of cortical activation during verb generation and picture naming[END_REF][START_REF] Flinker | Redefining the role of Broca's area in speech[END_REF][START_REF] Leuthardt | Temporal evolution of gamma activity in human cortex during an overt and covert word repetition task[END_REF][START_REF] Llorens | Intra-cranial recordings of brain activity during language production[END_REF][START_REF] Sahin | Sequential processing of lexical, grammatical, and phonological information within Broca's area[END_REF].

The coordination of mental operations has been framed differently in neurophysiological models than in cognitive models. Neurophysiological models often frame this coordination in terms of communication between brain regions through neuronal oscillations [START_REF] Fries | Rhythms for cognition: Communication through coherence[END_REF]. The timing and spectral characteristics of these oscillations provide the grounds for proposing and contrasting functional-connectivity metrics [START_REF] Wang | A systematic framework for functional connectivity measures[END_REF] that are claimed to characterize various cognitive processes. A wide gap remains between the cognitive and neurophysiological formulations (cf. [START_REF] Fries | Rhythms for cognition: Communication through coherence[END_REF][START_REF] Sigman | Brain mechanisms of serial and parallel processing during dual-task performance[END_REF][START_REF] Sigman | Brain mechanisms of serial and parallel processing during dual-task performance[END_REF], and many reports on highlevel cognitive processing, such as language processing, continue to make explicit reference to the original cognitive distinction between parallel and serial processing (e.g., [START_REF] Edwards | Spatiotemporal imaging of cortical activation during verb generation and picture naming[END_REF][START_REF] Sahin | Sequential processing of lexical, grammatical, and phonological information within Broca's area[END_REF]; see also Lachaux, Axmacher, Mormann, Halgren, & Crone, 2012, p. 287; for integrative views, see [START_REF] Friederici | Grounding language processing on basic neurophysiological principles[END_REF][START_REF] Strijkers | The cortical dynamics of speaking: Present shortcomings and future avenues[END_REF].

Empirically discriminating between serial and parallel cognitive architectures is notoriously difficult. Contrasting models with different underlying architectures often yield indistinguishable predictions, which is referred to as model mimicry, the ability of a model to account for data generated by a competing model (e.g., [START_REF] Logan | Parallel and serial processing[END_REF][START_REF] Meyer | Modern mental chronometry[END_REF][START_REF] Miller | Discrete and continuous models of human information processing: Theoretical distinctions and empirical results[END_REF]. Measures of behavioral performance (e.g., response times, or RTs) record the final product of cognitive processing. Therefore, one must rely on indirect reasoning and complementary qualitative or quantitative assumptions to infer the underlying temporal architecture on the basis of behavioral performance. Also, in both behavioral and neurophysiological studies, data are usually aggregated (averaged) before comparisons are made, which contributes to model mimicry [START_REF] Burle | Spatiotemporal dynamics of information processing in the Brain: Recent advances, current limitations and future challenges[END_REF]. Electroencephalography (EEG) and magnetoencephalography (MEG) can record human cortical dynamics at millisecond resolution. To compensate for relatively low signal-to-noise ratio (SNR), though, the vast majority of studies have focused on trial averages (e.g., typical event-related potential, or ERP, studies). Such averaging is critically based on the assumption that the response is identical (both cognitively and time-wise) across trials. When there is variability across trials (e.g., in cognitive potentials; [START_REF] Jung | Analyzing and visualizing single-trial event-related potentials[END_REF], the average response provides only a blurred view of the underlying processes [START_REF] Holm | Relationship of P300 single-trial responses with reaction time and preceding stimulus sequence[END_REF]. In short, the gain in SNR obtained by averaging is at the expense of spatiotemporal resolution, and therefore compromises what could be termed cognitive resolution (Miller, 1988, Section 2.4.4).

The main distinction between serial and parallel processing in language models is formulated at the level of a single prototypical trial, but averaged data have typically been used to test these models. To illustrate why this is problematic, we simulated an imaginary cognitive task with three processing stages (e.g., Donders, 1868[START_REF] Donders | On the speed of mental processes (W. G. Koster[END_REF]) that varied stochastically across trials; the duration of each processing stage was conceived as a random variable drawn from a uniform distribution. To form each trial, we combined durations of the three processing stages. To illustrate serial architecture, we combined the durations without temporal overlap, and to illustrate parallel architecture, we combined them with a fixed 30% temporal overlap. For each architecture, multiple trials were formed in this way. The resulting trial durations, depicted in Figure 1a, were quite different for the two architectures. Figure 1b shows the proportion of trials in which each stage took place as a function of time elapsed since trial onset. As this figure illustrates, the amount of parallel processing cannot be evaluated properly on the basis of averages; rather, cognitive activity would be better assessed at the level of single trials (for similar considerations, see [START_REF] Rey | Single trial analysis of field potentials in perception, learning and memory[END_REF].

The Current Study

In the current study, we measured the degree to which picture naming involves parallel processing. To do so, we used intracerebral activity in patients with epilepsy while they performed a picture-naming task. Intracerebral activity is recorded directly from cortex with a high SNR and can thus be fruitfully analyzed at the level of single trials. Within the complex recorded intracerebral signals, we focused on significant increases in activity detected in the high-gamma frequency band (80-150 Hz), as this is one of the indices thought to reliably reflect cognitive processing [START_REF] Fries | Rhythms for cognition: Communication through coherence[END_REF][START_REF] Lachaux | High-frequency neural activity and human cognition: Past, present and possible future of intracranial EEG research[END_REF]. Recording sites, located with spatial accuracy within a few millimeters, were grouped across patients in standardized anatomical regions of interest [START_REF] Mai | Atlas of the human brain[END_REF] on the basis of statistical consistency of the functional responses. This novel procedure was used to map all patients' signals onto a common parcellated brain template, thereby compensating for the partial coverage of each patient's brain and allowing a more integrated view (for a different approach, see [START_REF] Kadipasaoglu | Surface-based mixed effects multilevel analysis of grouped human electrocorticography[END_REF]. The resulting regions were tentatively linked to their respective cognitive functions on the basis of previous functional-imaging evidence (mostly from Price's, 2012, meta-analysis; see Table 1 for details). Finally, to quantify parallel processing (i.e., temporal overlap), we systematically computed the proportion of trials in which the activity of two regions coincided in time and assessed the statistical significance of these proportions. This ratio can be thought of as a variant of measures of discrete mutual information; our goal was not to add yet another metric of functional connectivity to the available list of methods [START_REF] Wang | A systematic framework for functional connectivity measures[END_REF], but to formulate and implement a metric that was as directly related as possible to the original parallel-processing cognitive hypothesis.

Method

Patients

The experimental protocol was approved by the relevant institutional review board (i.e., the Comité de Protection des Personnes). The patients (or their legal representatives in the case of minors) provided written informed consent. All patients were undergoing presurgical evaluation for pharmacologically intractable epilepsy at Hôpital de La Timone (a hospital in Marseille, France). Functional stereotactic exploration was performed using 5 to 11 depth electrodes (0.8 mm; Alcis, Besançon, France) each containing 10 to 15 recording sites (also referred to as contacts). Contacts were 2 mm long and separated from each other by 1.5 mm. The locations of the electrode implantations were strictly guided by clinical indications. Anticonvulsant therapy was reduced or withdrawn for the clinical exploration; however, the experiment was conducted only when a patient had been seizure-free for at least the preceding 12 hr.

A total of 23 epileptic patients for whom we had access to all anatomical and functional data volunteered to participate in the experimental protocol. Patients were excluded if they were younger than 10 years old (n = 2), had previously undergone brain surgery (n = 2), or had a pathology that could alter their language abilities (e.g., dyslexia; n = 2).

Procedure

The task consisted of naming 108 objects depicted as black-on-white line drawings, an activity completed independently of the procedures for diagnosing the patients' medical condition. The patients were seated in an electrically shielded room facing a display monitor at a distance such that the pictures, at the center of the screen, subtended an angle of 6° by 6°. An experimenter was present in the room in order to monitor the patients' performance and take written note of any response that was erroneous. For each patient, the pictures were pseudorandomly ordered with the constraints that show the proportion of trials in which each of the three stages is taking place as a function of time. Note that for any given point in time, the proportions for a given stage are similar for the serial model and the parallel model. Thus, if a dependent variable based on averages were used, strong qualitative differences in the temporal organization within single trials (a) might be evident only as quantitative, perhaps subtle, differences in the overlap of stages across trials (b). a.u. = arbitrary units. consecutive trials did not involve items from the same semantic category or words beginning with the same phoneme. On each trial, a single picture target was presented for 1,000 ms; a fixation cross was presented for 1,750 ms ± 350 ms (random jitter) between trials. The patients were instructed to name each pictured object aloud as fast as possible while avoiding errors; they were asked to remain silent if they did not recognize an object or could not come up with an answer.

Responses were recorded with a microphone (Audio-Technica ATR20, Tokyo, Japan) placed 13 cm in front of the patients. The experiment was controlled by E-Prime Version 2.0.1 (Psychology Software Tools, Pittsburgh, PA); the software automatically computed and logged the naming latencies relative to picture onset. We rejected trials on which the response was presumably anticipated (RT < 300 ms), delayed (RT > 1,600 ms), missing, or incorrect. Two patients who had less than 50% of the total number of trials remaining after these exclusions (19% and 46%, respectively) were dropped from analyses (see Table S1 in the Supplemental Material available online for details on the patients included in analyses).

Anatomical data

For all patients, both a structural preoperative MRI scan and an intraoperative computed tomography (CT) scan were acquired as part of the clinical routine. For each patient, we used the coregistration between MRI and CT obtained within the Leksell SurgiPlan software (Elekta, Stockholm, Sweden). Visual inspection of the fused images allowed us to precisely locate every contact [START_REF] Price | A review and synthesis of the first 20 years of PET and fMRI studies of heard speech, spoken language and reading[END_REF] Note: The theoretical framework used was a standard cognitive model of word production in picture naming involving multiple processing stages: visual, semantic, lexical, phonological, and articulatory (e.g., [START_REF] Indefrey | The spatial and temporal signatures of word production components: A critical update[END_REF][START_REF] Levelt | Spoken word production: A theory of lexical access[END_REF][START_REF] Rapp | Discreteness and interactivity in spoken word production[END_REF]. a The regions with consistently significant activity were linked to specific cognitive processes on the basis of a broad review of functionalimaging studies [START_REF] Price | A review and synthesis of the first 20 years of PET and fMRI studies of heard speech, spoken language and reading[END_REF]; when a region was not included or not discussed in detail in [START_REF] Price | A review and synthesis of the first 20 years of PET and fMRI studies of heard speech, spoken language and reading[END_REF], other studies were used. See the General Discussion for comments on this approach.

within each patient's anatomy. These locations were then visually classified by the neurologist on the basis of the human brain atlas [START_REF] Mai | Atlas of the human brain[END_REF], with only minor modifications to its parcellation. Left pars opercularis and left pars triangularis were merged into a single region, left inferior frontal gyrus. In addition, lateral basal and temporal regions, as well as mesial cingular areas, were divided into anterior, middle and posterior subregions. The latter subdivisions were intended to capture relevant functional distinctions made in Price's (2012) review of functional-imaging studies of language. Finally, all patients' contacts were mapped onto a common parcellated brain template.

All subsequent analysis was performed on activity in bipolar channels, calculated by subtracting activity recorded at one site from activity recorded at a neighboring site within the same electrode. When two regions were involved in a bipolar channel (i.e., the two contacts of the channel were classified as located in two different regions), we adopted a conservative approach of classifying the channel as belonging to both regions. All bipolar channels classified as outside the brain (by at least one of their contacts) were rejected from the analysis.

Functional data

Intracerebral EEG signals were recorded at a sampling rate of 1000 Hz using a 256-channel BrainAmp amplifier system (Brain Products GmbH, Munich, Germany). An acquisition band-pass filter was used to limit the bandwidth of the output signal to between 0.16 and 200 Hz. A scalp electrode placed in Fz was used as the recording reference.

Off-line preprocessing was performed with BrainVision Analyzer software (Brain Products GmbH, Munich, Germany). Data were then imported into MATLAB 2012a (The MathWorks, Natick, MA) for signal processing and statistical analysis. Cortex representations were created in Brainstorm [START_REF] Tadel | Brainstorm: A user-friendly application for MEG/ EEG analysis[END_REF], which is freely available for download under the GNU general public license (http://neuroimage.usc.edu/brainstorm).

Signal processing and group analysis

Signals were first processed for individual patients. Epochs were extracted from 1 s before to 2 s after the onset of each picture in the trials that were to be included in the analysis. Epochs containing epileptic spikes were removed by visual inspection (mean number of trials left for analysis = 87, range = 68-108), as were contacts with flat and noisy recorded activity. Bipolar channels were then computed as described earlier. Our analysis targeted activity in the high-gamma range, which has been repeatedly linked to cognitive processing [START_REF] Lachaux | High-frequency neural activity and human cognition: Past, present and possible future of intracranial EEG research[END_REF]. Timefrequency power was computed on consecutive 10-Hz bands between 80 and 150 Hz using a 7-cycle Morlet wavelet transform corresponding to a wavelet duration of approximately 50 ms at 80 Hz (a width that corresponds to 95% of the area under the Gaussian envelope). A baseline correction was applied on each trial at each 10-Hz band by computing a z score relative to activity during the baseline from 600 ms to 50 ms before picture onset (so as to exclude edge effects). Low frequencies have much higher power than high frequencies, and this correction procedure was adopted in order to compensate for this distribution of power across frequencies [START_REF] Buzsáki | The origin of extracellular fields and currents -EEG, ECoG, LFP and spikes[END_REF]. A series of power estimates were then obtained for each trial by summing the z-scored power across these bands.

The statistical significance of broadband power was assessed by computing one-sample Student's t tests (α = .05) across trials at each time stamp. To account for the multiple-comparison problem in the time domain, we estimated a minimum duration threshold for each patient using a bootstrap procedure. This procedure consisted of randomly selecting the same number of trials as in the original data set but with repetitions allowed and identifying periods with significant activity within the baseline window from -600 ms to -50 ms. The procedure was repeated 1,000 times, and the maximum numbers of contiguous points passing the significance threshold (corresponding to uncorrected p = .05) were pooled into a bootstrap distribution. A significance threshold corresponding to 95% of the histogram of maximum values was then obtained.

To compensate for the inevitable variability of electrode implantation sites across patients and the potential idiosyncrasies of the patients, we then combined signals from individual participants in a group analysis. The goal was to identify functional consistency within each region across patients, so as to take "a more global approach" (i.e., more global than using multiple single-case studies or focusing on a small set of brain structures) to examining "global brain dynamics" (Lachaux et al., 2012, pp. 291-292; Lachaux discusses strengths and weaknesses of this approach on these pages). Only regions that were sampled in at least 2 patients (see Fig. S1a in the Supplemental Material) and showed significant signal (see Fig. S1b in the Supplemental Material) were considered further. This procedure allowed us to integrate multiple patients' data (see also [START_REF] Kadipasaoglu | Surface-based mixed effects multilevel analysis of grouped human electrocorticography[END_REF]. To assess functional similarity across patients, we computed an average for each patient within each region within a time window spanning 400 ms before to 1,600 ms after picture onset. Similarity between the different patients' time courses of z-scored power within each region was measured by mean pairwise Pearson correlation coefficients. A heuristically determined criterion of r > .3 was used as a threshold for functional similarity.

For all regions that survived both the significance and the correlation criteria (Fig. S1c in the Supplemental Material, Table 1), the average z scores across patients were computed. For visualization purposes, we colorcoded the regions according to the temporal order in which the averaged z score first exceeded a threshold of 3 (colormap jet in MATLAB). The brain template used for this visualization (and for Fig. S1 in the Supplemental Material) was the ICBM152 as implemented in Brainstorm, and the regions were either adapted from the predefined Destrieux and Desikan-Killiany atlases in Brainstorm or created from Brainstorm's user interface ("scouts" manipulation).

Overlap analysis

The extent of overlapping activity between pairs of regions at the level of single trials was estimated on the basis of thresholded single-trial signals. We transformed the single-trial signals (z scores) into binary significance masks (on/off) by keeping only points that were above a threshold obtained using the local false-discovery rate (corresponding to q < .2; [START_REF] Efron | Local false discovery rates[END_REF]; we note that the reported pattern of results was not strongly modified when more liberal thresholds (e.g., q < .6 were used). When several contacts sampled the same region within a patient, they were grouped together by taking the union of their significance masks in order to form one mask per region.

Then, for each pair of regions at each time sample, the overlap was defined as the ratio between the number of trials that were significantly "on" in both regions (Index 1) and the number of trials that were significantly "on" in at least one of the regions (Index 2). Index 1 effectively captured the simultaneous significant activity present in both regions. Including Index 2 ensured that trials in which both regions were not active were not considered, and that a sufficient number of trials were taken into account in the overlap measure (a minimum of 20 trials was adopted by examining the data histogram).

In order to smooth possibly excessive variability while preserving the temporal resolution of interest, we averaged both indices within 50-ms sliding windows in the poststimulus period (0-1,600 ms). For each 50-ms window, the result was summarized in a matrix showing the maximal overlap strength of all the pairs of regions for which data were available and the 20-trial threshold was met. To assess the significance of overlap, we performed a nonparametric permutation analysis. For each pair of regions, we created a surrogate data set for one of the regions by shuffling the "on" and "off" segments across time and across trials. Overlap in the surrogate data was measured, and the maximum value was kept. This procedure was repeated 1,000 times, and a significance threshold corresponding to 95% of the histogram of maximum values was obtained.

Results

For the 15 patients who remained after the exclusion criteria were applied (see Table S1 in the Supplemental Material), we analyzed intracerebral data from a total of 146 electrodes (mean per patient = 9.8, minimum = 6, maximum = 13) that recorded activity in 1,347 cortical sites, 720 located in the left hemisphere and 627 located in the right hemisphere.

Of the 99 regions that emerged from the classification of recording sites, 49 were sampled in at least 2 patients (Fig. S1a in the Supplemental Material). Significance thresholding revealed 38 regions showing significant high-gamma activity in at least 2 patients (Fig. S1b). A total of 23 regions (13 in the left hemisphere) showed significantly consistent activity across patients, r > .3 (Fig. S1c). The activity in each of these regions was tentatively associated with a particular cognitive process on the basis of previous evidence (mostly from [START_REF] Price | A review and synthesis of the first 20 years of PET and fMRI studies of heard speech, spoken language and reading[END_REF] see Table 1 for full details). This procedure pointed to 13 different cognitive processes.

Average patterns

On average, high-gamma neural activity elicited by the stimuli started around 80 ms after picture onset in posterior regions associated with visual perceptual processing (Figs. 2a and2c). The activity then spread to basal temporal regions bilaterally, up to mid fusiform areas and the left hippocampus. Regions in the dorsal pathway, such as the left inferior frontal gyrus, left precentral gyrus, and right supramarginal gyrus, were then engaged. The activity finally reached the mid part of the superior temporal gyrus, which presumably is involved in auditory perception of the spoken word.

We sorted the time series by increasing order of onset time, which allowed us to identify three subsets of regions that roughly followed different temporal patterns (see also [START_REF] Edwards | Spatiotemporal imaging of cortical activation during verb generation and picture naming[END_REF]: Subset 1 contained 11 bilateral basal-posterior regions (from left occipital gyrus to left posterior fusiform gyrus; see Fig. 2a) in which activity reached a z-score amplitude of 4 to 6 between 50 and 200 ms and slowly decreased back to zero at around 800 ms. The activity present around 1,100 ms in Figure 2a corresponds to the brain response to the disappearance of the visual stimulus (which was presented for 1,000 ms). Subset 2 involved five left-hemisphere regions (from left parieto-occipital transition zone to left hippocampus; see Fig. 2a) in which activity reached a z-score amplitude of 2 to 4 between 200 and 800 ms after picture onset. Subset 3 contained seven frontal and temporal regions all patients. For an explanation of the region labels, see Table 1. R = right; L = left.

from both hemispheres (left precentral gyrus to left anterior middle temporal gyrus; see Fig. 2a) where z-score amplitude started to increase at around 800 ms and remained significant until the end of the time window. Activity in Subsets 1 and 2 largely overlapped in time, whereas activity in Subset 3 overlapped much less with activity in the two other subsets. Within each subset, the regions follow very similar temporal patterns, and therefore activity largely overlapped in time. However, this overlap may have stemmed in part from combining data from different patients and different trials with variable durations (see Fig. 2b for the distribution of response times), so we next analyzed overlap within single trials.

Overlap within single trials

To assess the temporal overlap between cognitive functions more specifically, we computed activity overlap within pairs of regions at the level of single trials. This required that activity for each pair of regions was measured within patients. We included all pairs of regions that were sampled in at least 1 patient (not 2, as in the preceding analyses), bearing in mind that the activity in each region was known to be consistent in at least 2 patients. Of the 107 pairs of regions for which we had data, 81 had enough trials to be further analyzed.

We tested the parallel-processing hypothesis with a metric that quantified, over time, the proportion of trials with overlapping activity for each pair of regions (see Method). This metric showed the strongest overlap in occipital and temporal regions, mainly between 80 ms and 400 ms after stimulus onset. There was also overlapping activity around 1,200 ms, which was maximal for superior temporal regions bilaterally (see Fig. S2 in the Supplemental Material). A histogram of the distribution of overlap across regions and the entire poststimulus period revealed a maximum overlap value of .46 (median = .075; Fig. 3b). The highest overlap values were observed in posterior regions (right posterior fusiform gyrus-right occipital gyrus; left striate area-left occipital gyrus) and in superior temporal regions (left median superior temporal gyrus-left anterior transverse temporal gyrus; right anterior transverse temporal gyrus postriright posterior superior temporal gyrus; Fig. 3c). Eight pairs of regions survived a strict permutation test of the temporal overlap, exceeding the 95% threshold in the permutation distribution of the maximum overlap statistic (black squares in Fig. 3c).

Substantially less overlap was detected in many other pairs of regions (Fig. 3b). This was the case for some pairs of regions located within the ventral stream. For example, there was very limited overlap between the left lingual and posterior fusiform gyri, on the order of .10 (Fig. 3c). Limited, subthreshold, concurrent activation in the ventral stream was also found for the following pairs of regions: left occipital gyrus and left hippocampus, left striate area and left hippocampus, and left occipital gyrus and posterior fusiform gyrus (Fig. S2); in all three cases, the overlap was less than .10. The dorsal stream has been hypothesized to be strongly left-hemisphere dominant. However, our sampling revealed right supramarginal gyrus activity that was moderately (subthreshold) concurrent with activity in both the left (overlap = .17) and right (overlap = .5) lingual gyri, as well as with activity in the right posterior superior temporal gyrus (overlap = ~.12) and right anterior transverse temporal gyrus (overlap = ~.08). In short, our analysis revealed moderate temporal overlap in two groups of regions roughly belonging to the dorsal and ventral pathways, and the strongest overlap in sensory cortices (see Fig. 3d for a circular diagram showing which pairs of regions had overlap values greater than .10 and brain images showing the anatomical location of each region, along with color coding indicating the onset of activity in each region).

Discussion

To quantify parallel processing in word production, we investigated the temporal overlap between cognitive processes indexed by neurophysiological activity in different brain regions. Given the importance of addressing this issue at the level of single trials, we decided to measure activity using intracerebral EEG signal (which has high spatiotemporal resolution and very high SNR) recorded during a picture-naming task known to involve multiple processing levels. We jointly analyzed the data from 15 patients with epilepsy, targeting signal fluctuations in the high-gamma band that were significant and consistent across at least 2 patients.

Averaged and single-trial data offer different views of temporal overlap

Signal averages revealed 23 brain regions with significant activity. The spatiotemporal patterns (Fig. 2) and our cognitive interpretation of the data (Table 1) are in agreement with many previous findings [START_REF] Edwards | Spatiotemporal imaging of cortical activation during verb generation and picture naming[END_REF][START_REF] Indefrey | The spatial and temporal signatures of word production components: A critical update[END_REF][START_REF] Leuthardt | Temporal evolution of gamma activity in human cortex during an overt and covert word repetition task[END_REF][START_REF] Llorens | Intra-cranial recordings of brain activity during language production[END_REF]; see also [START_REF] Munding | On the cortical dynamics of word production: A review of the MEG evidence[END_REF][START_REF] Price | A review and synthesis of the first 20 years of PET and fMRI studies of heard speech, spoken language and reading[END_REF]. Our examination of the onset time of detected activity led us to distinguish three groups of concurrently activated regions: early, middle, and late, roughly corresponding to, respectively, stimulus processing, semantic association and selection processes, and motor and auditory processes (see Table 1). Within each subset, there was substantial temporal overlap among regions and, therefore, among the (presumed) corresponding cognitive processes. These findings might prompt interpretations in favor of , from longest at the top to shortest at the bottom. Overlap was considered to be computable when the total number of trials showing significant activity (third row) was above 20 (red horizontal line). Superimposing the data from the two regions (fourth row) reveals visually the times at which the two regions had overlapping activity (in red). Such temporal overlap between regions was measured by dividing the number of trials with concurrent significant activity by the total number of trials with significant activity. The plot in the fifth row shows the temporal overlap (red) between the left occipital gyrus and the left fusiform gyrus as a function of time in this patient. White indicates periods for which there were (continued on next page) parallel processing within each subset, and perhaps serial processing across the subsets (especially Subset 3 relative to Subsets 1 and 2; [START_REF] Edwards | Spatiotemporal imaging of cortical activation during verb generation and picture naming[END_REF][START_REF] Leuthardt | Temporal evolution of gamma activity in human cortex during an overt and covert word repetition task[END_REF][START_REF] Sahin | Sequential processing of lexical, grammatical, and phonological information within Broca's area[END_REF].

At the level of single trials, the temporal overlap measured was rather low. There was no obvious indication of serial processing (i.e., activity in one region systematically preceding activity in another region). Instead, activity was distributed over the periods of interest (for an example, see Fig. 3a and Fig. S2 in the Supplemental Material). Temporal overlap in activity between a single pair of regions was never found on more than 50% of the trials. Furthermore, in most cases (i.e., in most 50-ms time windows), overlap occurred on about 5% of trials with significant activity (i.e., temporal overlap of .05; Fig. 3b). The highest temporal overlap was observed within sensory cortices: first in the regions presumably capturing progressive visual decoding of the pictures and then in the regions likely capturing the decoding of external auditory feedback from the patients' own voices [START_REF] Gourévitch | Temporal envelope processing in the human auditory cortex: Response and interconnections of auditory cortical areas[END_REF]. Substantially less overlap was detected in pairs of regions located within the ventral stream (Fig. 3c), which is known to be involved in visual stimulus processing, object recognition [START_REF] Malach | The topography of high-order human object areas[END_REF], and lexico-semantic access [START_REF] Fonseca | Electrophysiological study of the basal temporal language area: A convergence zone between language perception and production networks[END_REF]. The dorsal stream, which is mainly involved in mapping phonological onto articulatory representations [START_REF] Hickok | Computational neuroanatomy of speech production[END_REF], has been hypothesized to be strongly left-hemisphere dominant; however, our sampling revealed that right-hemisphere regions were moderately coactivated (subthreshold level).

The observed values of within-trial coactivation (all under .50 and generally under .10) might be taken at face value to reflect the actual levels of simultaneous processing occurring within the pairs of regions. Alternatively, these observed values might reflect the maximum degree of processing overlap that our analysis could detect. 1 In particular, the distribution of single-trial activity might have been influenced by the noise present in the data. It is important to note that the brain activities we report survived a strict nonparametric permutation threshold, such that each region and its level of SNR provided its own control. Furthermore, the estimated levels of overlap were largely unaffected by lowering the statistical threshold, which suggests that moderate overlap was not due only to excess false negatives (i.e., undetected activity that prevented the detection of overlap). With these considerations in mind, we conclude that the single-trial results match the predictions of the parallel-processing model in sensory cortices but not elsewhere, despite the "coactivations" that were clearly present in the (more blurred) averaged data.

Comparison with other connectivity metrics

Our goal to quantify concurrent activity in parallel cognitive architectures led us to define a single-trial, temporal-overlap metric. This metric captured simultaneous (i.e., zero-lag) activation of two regions (for which enough trials were available) at a given point in time. Although the metric was defined on the basis of single time samples of single trials, it still integrated signal over time (and frequency). Our estimates smeared the signal over a temporal window of about 50 ms. Therefore, the window of concurrent activation that was tested was in the range of a few 10s of milliseconds. This is in the range posited by most recent hypotheses about communication through phase or amplitude correlation; interareal communication has been described as occurring in the range of the gamma band [START_REF] Fries | Rhythms for cognition: Communication through coherence[END_REF]. Using a cognitive task very similar to ours (word reading and repetition), [START_REF] Flinker | Redefining the role of Broca's area in speech[END_REF] recently observed interareal communication at longer delays. In that study, Granger causality estimates indicated the direction of the respective influences of Broca's area, superior temporal gyrus, and motor regions. Note that these temporal estimates characterize how information flowed between pairs of regions (with causal influence present within 200-ms delays), rather than the degree of concurrent activation. Indeed, the main findings of Flinker et al. are compatible with both serial and parallel architectures. On the other hand, the estimates of concurrent activation reported here do not address whether there is any causal relationship between concurrent processes ("contingent stages" ;[START_REF] Miller | Discrete and continuous models of human information processing: Theoretical distinctions and empirical results[END_REF]).

Limitations and perspectives

Our study has many of the well-known limitations that are inherent to studies of intracerebral activity in epileptic patients [START_REF] Lachaux | High-frequency neural activity and human cognition: Past, present and possible future of intracranial EEG research[END_REF]. Our group analysis compensated, to some extent, for strong interindividual variability. It yielded a fair spatial sampling (see Figs. S1a more than 20 trials, and gray indicates periods for which there were fewer than 20 trials. The histogram in (b) shows the distribution of 50-ms sliding windows (step = 1 ms) across regions according to the amount of temporal overlap they exhibited in the poststimulus period (0-1,600 ms). The color coding in (c) indicates the maximum temporal overlap observed between all pairs of regions in the poststimulus period. The black squares highlight pairs for which the overlap value exceeded the 95% threshold in the permutation distribution of the maximum overlap statistic. Plus signs indicate that there were not enough trials to compute the temporal overlap. The circular diagram of the 23 regions in (d) shows which pairs had overlap values greater than .10. The color coding in the brain images corresponds to the colors used in the diagram and indicates each region's onset of activity (see also Fig. 2). For an explanation of the region labels, see Table 1. Fig. 3. (continued) and S1b in the Supplemental Material), but at the expense of anatomical details provided by specific stereotactic implantations. In addition, potentially important variability in the signal was excluded from analyses as a result of the cross-patient consistency constraint. Single-trial concurrent coactivation could be assessed for 81 pairs of regions, but often only within 1 patient (48 pairs of regions).

Our focus on high-gamma activity was strongly motivated by current thinking that it is involved in cognitive processing and interregional communication, and we did not explore other potentially meaningful signal components [START_REF] Canolty | Spatiotemporal dynamics of word processing in the human brain[END_REF][START_REF] Gaona | Nonuniform high-gamma (60-500 Hz) power changes dissociate cognitive task and anatomy in human cortex[END_REF]; future research is needed to broaden the search space (e.g., to other frequency components). Data from single trials were reduced to a binary significance mask using a local false-discovery rate threshold, with the consequent risk of losing information. The motivations behind our approach were both statistical and theoretical: to establish periods of significant activity within every trial and to provide the most literal test of the parallel-processing hypothesis.

Finally, to test a hypothesis about cognitive processing with neurophysiological data, we used an essentially anatomical approach to associate cognitive processes with brain regions [START_REF] Poldrack | Can cognitive processes be inferred from neuroimaging data[END_REF][START_REF] Price | A review and synthesis of the first 20 years of PET and fMRI studies of heard speech, spoken language and reading[END_REF]. This common practice has been identified as a case of reverse inference [START_REF] Poldrack | Can cognitive processes be inferred from neuroimaging data[END_REF], and is particularly detrimental when inferences are made across cognitively different tasks. However, the cognition-brain associations we used were quite specific, came from a meta-analysis [START_REF] Price | A review and synthesis of the first 20 years of PET and fMRI studies of heard speech, spoken language and reading[END_REF], and involved the same task we used (picture naming) or very similar tasks (e.g., word reading). Additional links between cognitive processes and brain regions might be established through cognitive contrasts and analysis of activity patterns, as it is possible that cognitive processes have implementations distributed across regions, such that advanced multivariate methods [START_REF] King | Characterizing the dynamics of mental representations: The temporal generalization method[END_REF] would have to be considered as means to detect these links.

Conclusion

Our study assessed the parallel-processing hypothesis at an unprecedented level of granularity in a language task. The group analysis of averaged data provided a detailed functional map of word production. At the level of single trials, significant activity was significantly concurrent in visual and auditory regions, but rarely concurrent among many other task-relevant pairs of regions. Within the boundaries of our exploration, our novel methodology reveals some limits to the degree to which word production involves the parallel-processing architecture postulated in major models of language processing.

Fig. 1 .

 1 Fig.1. Schematic representation of the consequences of averaging for distinguishing between serial processing (left) and parallel processing (right). The graphs in (a) present a raster view of the trials in our simulation ordered by response time. Each line shows the time course of the three theoretical processing stages during one trial. The graphs in (b) show the proportion of trials in which each of the three stages is taking place as a function of time. Note that for any given point in time, the proportions for a given stage are similar for the serial model and the parallel model. Thus, if a dependent variable based on averages were used, strong qualitative differences in the temporal organization within single trials (a) might be evident only as quantitative, perhaps subtle, differences in the overlap of stages across trials (b). a.u. = arbitrary units.

Fig. 2 .

 2 Fig. 2. Average brain activity during picture naming. The raster plot in (a) shows the time course of activity in the high-gamma band, averaged across patients, for all regions with at least 2 patients with similar time courses (r > .3); each line shows the temporal course for a single region. A threshold z score of1.96 was used to identify regions with significant activity. The regions are ordered according to the first time bin in which the z score exceeded 3; this order is color coded at the left and illustrated in the brain maps in (c). The number of contacts and patients that were used to compute each averaged signal is provided at the right (contacts/patients). The histogram in (b) shows the distribution of response times across

Fig. 3 .

 3 Fig.3. Temporal overlap between consistently active regions. The graphs in (a) illustrate an example of overlap between a pair of regions within a single patient. The two raster plots at the top show the temporal course of significant neural activity (color) in single trials for contacts located in the left occipital gyrus (L.OcG; first row) and the left fusiform gyrus (L.pFuG; second row). The trials are ordered by responses time (indicated by the black dots), from longest at the top to shortest at the bottom. Overlap was considered to be computable when the total number of trials showing significant activity (third row) was above 20 (red horizontal line). Superimposing the data from the two regions (fourth row) reveals visually the times at which the two regions had overlapping activity (in red). Such temporal overlap between regions was measured by dividing the number of trials with concurrent significant activity by the total number of trials with significant activity. The plot in the fifth row shows the temporal overlap (red) between the left occipital gyrus and the left fusiform gyrus as a function of time in this patient. White indicates periods for which there were

Table 1 .

 1 The 23 Regions With Consistently Significant Activity

	Subset and
	abbreviation	Region	Associated cognitive process a
	Subset 1
	L.OcG	Left occipital gyrus	Visual association (Price, 2012)
	L.17	Left striate area	Primary visual processes (Price, 2012)
	R.OcG	Right occipital gyrus	Visual association (Price, 2012)
	R.SMG	Right supramarginal gyrus	Articulatory loop (Price, 2012)
			(typically reported for the left hemisphere)
	R.pFuG	Right posterior fusiform gyrus	No straightforward correspondence
	L.LgG	Left lingual gyrus	Visual association (Price, 2012)
	R.LgG	Right lingual gyrus	Visual association (Price, 2012)
	R.mPHG	Right medial parahippocampal gyrus	Object recognition (Malach, Levy, & Hasson, 2002)
	R.mFuG	Right medial fusiform gyrus	No straightforward correspondence
	R.pITG	Right posterior inferior temporal gyrus	No straightforward correspondence
	L.pFuG	Left posterior fusiform gyrus	Transition from visual form to semantic processing
	Subset 2
	L.POTZ	Left parieto-occipital transition zone	No straightforward correspondence
	L.IFG	Left inferior frontal gyrus	Semantic decision and articulation (Price, 2012)
	L.mFuG	Left medial fusiform gyrus	Transition from visual form to semantic processing (Price, 2012)
	L.PCL	Left paracentral lobule	Sequencing of motor movements (Price, 2012)
	L.Hi	Left hippocampus	Associative links (Llorens et al., 2016)
	Subset 3
	L.PrG	Left precentral gyrus	Orofacial motor activity (Price, 2012)
	L.TTG1	Left anterior transverse temporal gyrus	Primary auditory processes (Price, 2012)
	R.SFGL	Right superior frontal gyrus	Semantic word selection (Price, 2012)
			(typically reported for the left hemisphere)
	L.mSTG	Left medial superior temporal gyrus	Auditory processing of complex sounds (Price, 2012)
	R.pSTG	Right posterior superior temporal gyrus	Auditory association (Schuhmann, Schiller, Goebel, & Sack, 2012)
			(linked to self-monitoring in the left hemisphere)
	R.TTG1	Right anterior transverse temporal gyrus	Primary auditory processes (Price, 2012)
	L.aMTG	Left anterior middle temporal gyrus	Semantic association

Funding

We acknowledge the support of European Research Council FP7/2007-2013 Grant 263575, as well as Agence Nationale de la Recherche Grant Multimodel 2010 BLAN 0309 01.

Action Editor

Matthew A. Goldrick served as action editor for this article.

Author Contributions

Declaration of Conflicting Interests

The authors declared that they had no conflicts of interest with respect to their authorship or the publication of this article.

Supplemental Material

Additional supporting information can be found at http:// journals.sagepub.com/doi/suppl/10.1177/0956797616681296 Note 1. External evidence could be useful in assessing this alternative. For example, sensory cortices could be expected to have a very high degree of concurrent activation given intrinsic connectivity patterns reported during resting-state recordings [START_REF] Mantini | Electrophysiological signatures of resting state networks in the human brain[END_REF].