Anatomical localization of the sEEG electrodes for each patient (N=16).

Paradigm and coordination indexes.

(A) Top : illustration of one trial of the interactive synchronous speech repetition task (orange: virtual partner speech; blue: participant speech; stimulus papi m’a dit repeated 10 times ; only the 10 first seconds are represented). Bottom : the four speech utterances used in the task and the experimental procedure. (B) Speech signals processing stages. The top panel corresponds to the speech envelope, the second to the phase of speech envelope and the third panel to the phase difference between VP and participant speech envelopes, illustrating the coordination dynamics along one trial. (C) Left: distributions of verbal coordination index (phase locking values between VP and participant speech envelopes, for each trial) for all participants (top) and patients . Right: boxplots for control participants (top) and patients showing the trial-averaged verbal coordination index as a function of the virtual partner parameters (in-phase coupling vs coupling with a 180 shift) .

Power spectrum analyses.

Each dot represents a significant channel (permutation test). (A) Global activity (Task versus Rest) for each frequency band. The activity is expressed in % of power change compared to resting. (B) Behaviour-related activity: r values of the Spearman correlation across trials between the iEEG power and the verbal coordination index (VCI). (C) The proportion of significant channels in the task vs rest (orange), in the brain-behaviour correlation (green) and in both (blue). The percentage in the center indicates the overall proportion of significant channels from the three categories (with respect to the total number of channels).

Group analysis by regions of interest.

(A) Regions of interest (ROI) defined according to the cluster analysis (see Figure S2), the delimitation of regions is based on the Brainnetome atlas. (B) For each ROI, boxplots illustrate, in red, channels with significant global power changes (HFa, task vs rest) and, in blue, their corresponding r values (correlation between HFa power and verbal coordination index, VCI). Red and blue stars indicate a significant difference from a null distribution. Dots represent independent iEEG channels. The « n » below each region of interest specifies the number of patients. STG : superior temporal gyrus ; IPL : inferior parietal lobule ; IFG : inferior frontal gyrus ; BA : brodmann area.

Phase-amplitude coupling between virtual partner speech signal or coordination dynamics and high frequency activity (HFa).

A) Representation of the increase in PAC expressed in % compared to surrogates when using the virtual partner speech (left) or the coordination dynamics (phase difference between VP and patient, right).The shaded (slight blue) area corresponds to the location of the IFG BA44. B) PAC values for VP (in red) and phase difference (in blue) by regions of interest. Statistical difference calculated using paired wilcoxon test (STG BA41/42 : p=0.01 ; STG BA22 : p=0.004 ; IPL BA40 : p=0.6 ; IFG BA44 : p=0.02). Y-axis range has been adjusted to better illustrate the contrast between VP speech and coordination dynamics.

Distribution of verbal coordination index across patients.

For each of the sixteen patients, this figure depicts the histogram of the coordination index for all trials (in blue) as well as the null distribution (random phase shift) computed using 500 permutations per trial (in red).

Cluster analysis (silhouette score).

Illustration of the mean silhouette scores according to the number of clusters for global activity (in red) and behaviour-related activity (correlation between power changes and coordination index, in blue). The highest silhouette score was obtained for five clusters in the HFa range for behaviour-related activity (framed in full black square). Bottom right : spatial cluster representation for the highest mean silhouette score value in HFa range.