# A SENtence Supramodal Areas AtlaS (SENSAAS) based on multiple task-induced activation mapping and graph analysis of intrinsic connectivity in 144 healthy right-handers

Loic Labache, Marc Joliot, Jérôme Saracco, Gaël Jobard, Isabelle Hesling, Laure Zago, E. Mellet, Laurent Petit, Fabrice Crivello, Bernard Mazoyer, et al.

► To cite this version:

Loic Labache, Marc Joliot, Jérôme Saracco, Gaël Jobard, Isabelle Hesling, et al. A SENtence Supramodal Areas AtlaS (SENSAAS) based on multiple task-induced activation mapping and graph analysis of intrinsic connectivity in 144 healthy right-handers. Brain Structure and Function, 2019, 224 (2), pp.859-882. 10.1007/s00429-018-1810-2. hal-01963242

HAL Id: hal-01963242

https://hal.science/hal-01963242

Submitted on 13 Nov 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

# A SENtence Supramodal Areas AtlaS (SENSAAS) based on multiple task-induced activation mapping and graph analysis of intrinsic connectivity in 144 healthy right-handers

L. Labache<fnref n="1" />,<fnref n="2" />,<fnref n="3" />,<fnref n="4" />,<fnref n="5" /> D · M. Joliot<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> D · J. Saracco<fnref n="5" />,<fnref n="6" /> D · G. Jobard<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> D · I. Hesling<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> D · L. Zago<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> D · E. Mellet<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> D ·

L. Petit<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> D · F. Crivello<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> D · B. Mazoyer<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> D · Nathalie Tzourio-Mazoyer<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> D

Received: 5 March 2018 / Accepted: 1 December 2018 / Published online: 7 December 2018

© The Author(s) 2018

## Abstract

We herein propose an atlas of 32 sentence-related areas based on a 3-step method combining the analysis of activation and asymmetry during multiple language tasks with hierarchical clustering of resting-state connectivity and graph analyses. 144 healthy right-handers performed fMRI runs based on language production, reading and listening, both with sentences and lists of over-learned words. Sentence minus word-list BOLD contrast and left-minus-right BOLD asymmetry for each task were computed in pairs of homotopic regions of interest (hROIs) from the AICHA atlas. Thirty-two hROIs were identified that were conjointly activated and leftward asymmetrical in each of the three language contrasts. Analysis of resting-state temporal correlations of BOLD variations between these 32 hROIs allowed the segregation of a core network, SENT_CORE including 18 hROIs. Resting-state graph analysis applied to SENT_CORE hROIs revealed that the pars triangularis of the inferior frontal gyrus and the superior temporal sulcus were hubs based on their degree centrality (DC), betweenness, and participation values corresponding to epicentres of sentence processing. Positive correlations between DC and BOLD activation values for SENT_CORE hROIs were observed across individuals and across regions regardless of the task: the more a SENT_CORE area is connected at rest the stronger it is activated during sentence processing. DC measurements in SENT_CORE may thus be a valuable index for the evaluation of inter-individual variations in language areas functional activity in relation to anatomical or clinical patterns in large populations. SENSAAS (SENtence Supramodal Areas AtlaS), comprising the 32 supramodal sentence areas, including SENT_CORE network, can be downloaded at http://www.gin.cnrs.fr/en/tools/.

**Keywords** fMRI · Left hemisphere · Sentence · Language production · Reading · Speech comprehension · Intrinsic connectivity · Graph analysis · Right-handers · Resting-state

**Electronic supplementary material** The online version of this article (https://doi.org/10.1007/s00429-018-1810-2) contains supplementary material, which is available to authorized users.

<figure><img src="image_3.png" /></figure>

Nathalie Tzourio-Mazoyer
nathalie.tzourio-mazoyer@u-bordeaux.fr

## Abbreviations

<table><tr><td>aHCA</td><td>Agglomerative hierarchical cluster analysis</td></tr><tr><td>AICHA</td><td>Atlas of intrinsic connectivity of homotopic areas</td></tr><tr><td>AU</td><td>Approximately unbiased p value</td></tr><tr><td>BC</td><td>Betweenness centrality</td></tr><tr><td>DC</td><td>Degree centrality</td></tr><tr><td>DMN</td><td>Default mode network</td></tr><tr><td>FIR</td><td>Finite impulse response</td></tr><tr><td>fMRI</td><td>Functional magnetic resonance imaging</td></tr><tr><td>FWHM</td><td>Full width at half maximum</td></tr><tr><td>hROIS</td><td>Homotopical regions of interest</td></tr><tr><td>LISN<sub>SENT-WORD</sub></td><td>Sentence minus word-list contrast for listening task</td></tr></table>

<footnote n="1">¹ Univ. Bordeaux, IMN, UMR 5293, 33000 Bordeaux, France</footnote>

<footnote n="2">² CNRS, IMN, UMR 5293, 33000 Bordeaux, France</footnote>

<footnote n="3">³ CEA, GIN, IMN, UMR 5293, 33000 Bordeaux, France</footnote>

<footnote n="4">⁴ Univ. Bordeaux, IMB, UMR 5251, 33405 Talence, France</footnote>

<footnote n="5">⁵ INRIA Bordeaux Sud-Ouest, CQFD, UMR 5251, 33405 Talence, France</footnote>

<footnote n="6">⁶ Bordeaux INP, IMB, UMR 5251, 33405 Talence, France</footnote>
<table><tr><td>MNI</td><td>Montreal Neurological Institute</td></tr><tr><td>pIndex</td><td>Participation index</td></tr><tr><td>PROD<sub>SENT-WORD</sub></td><td>Sentence minus word-list contrast for production task</td></tr><tr><td>READ<sub>SENT-WORD</sub></td><td>Sentence minus word-list contrast for reading task</td></tr><tr><td>rs-fMRI</td><td>Resting-state fMRI</td></tr><tr><td>SD</td><td>Standard deviation</td></tr><tr><td>SPM</td><td>Statistical parametric mapping</td></tr></table>

## Introduction

Defining language areas is a complex enterprise because of the numerous possible approaches currently available to identify language-related regions. The gold standard is to consider that language areas correspond to regions wherein lesions lead to aphasia. Even when limiting the definition of language areas to that of essential language areas, different identification methods exist that provide various kinds of information. Wada testing allows identification of the hemisphere controlling language but does not provide regional information (Wada and Rasmussen 1960). By contrast, surgical cortical stimulation studies have documented left hemisphere language areas in large samples of patients (Ojemann et al. 1989; Tate et al. 2014), but such mapping of eloquent areas is still limited to the cortical regions available to the neurosurgeon and is conducted in patients having potentially modified language organization. The probabilistic mapping of lesions combined with fine-grained aphasic patient evaluations of language performance have provided the community with very accurate descriptions of essential language areas (Dronkers and Ogar 2004; Dronkers et al. 2004) although this very important approach does not reveal how these cortical areas are organized in networks. Because each multiple cortical area altered by a given pathology is not involved in the language deficit, the comprehensive identification of language areas from lesions is a complex issue [see (Genon et al. 2018a, b) for a review].

Functional neuroimaging provides a way to map multiple areas activated during the completion of various language tasks in a large number of individuals. Furthermore, neuroimaging methodology is very efficient at compiling results obtained in multiple laboratories across the world, thereby allowing meta-analyses across laboratories that provide the location of areas activated at an acceptable spatial resolution within a common normalization space for a variety of language tasks. Similar to the results obtained with cortical stimulation (Ojemann et al. 1989), meta-analyses of neuroimaging data have provided the landscape of the left hemisphere cortical areas involved in language tasks in healthy individuals, which covers nearly the entire hemisphere surface (Price 2000, 2010, 2012; Vigneau et al. 2006).

Despite the vast amount of information obtained from the methods cited above, an atlas of left hemisphere language areas in healthy individuals having a typical left-hemisphere dominance for language is still lacking, and with respect to language areas, the absence of a consensus is clear. The posterior part of the superior temporal gyrus and the supramarginal gyrus (Tomasi and Volkow 2012; Klingbeil et al. 2017) are phonological regions that can be found under the label “Wernicke’s area”, while lesion-based studies (Dronkers and Ogar 2004; Yourganov et al. 2015) as well as lesion studies in association with activation studies (Saur et al. 2006) have shown deep aphasia associated with lesions of the posterior region of the middle temporal gyrus and superior temporal sulcus (Binder 2015, 2017). There is greater consistency concerning the location of frontal language areas under the label of Broca’s area because its original definition was anatomical. Most people define Broca’s area as the pars triangularis of the inferior frontal gyrus (Clos et al. 2013; Friederici and Gierhan 2013; Yourganov et al. 2015). However, the extent of Broca’s area in the left frontal lobe varies, and the anterior insula (Baldo et al. 2011) is sometimes added, as reviewed in Amunts (Amunts and Zilles 2012). Moreover, posterior lesions can also lead to Broca’s aphasia (Richardson et al. 2012), demonstrating that these anterior and posterior language poles work tightly together. This relation enhances the importance of networks in cognitive processing, as defined by Fuster (Fuster and Bressler 2012). An atlas of language areas and networks in healthy individuals would thus be a useful tool, especially when individual task-induced mapping is not available. This atlas would be especially helpful for patients having difficulties completing language tasks and for the exploration of genetic language bases in large cohorts of individuals, in cohorts targeting normal or pathological brains, including those with developmental pathologies, and/or in individuals mapped for their anatomy and/or resting state while not performing a language task (Thompson et al. 2017).

To elaborate such an atlas, increasing the specificity for language areas is important because as uncovered by lesion studies, not all areas revealed by task-induced activation studies are essential language areas. Components of the task, such as monitoring, selecting, and holding the instructions, as well as paralinguistic processing, such as context, emotional and prosodic processing, are responsible for activations that exceed the essential language areas of the left hemisphere. The strong right-hemisphere activations observed with functional imaging during various language tasks have even led some authors to claim that neuroimaging methods are not adequate to map language regions (Sidtis 2007). One way to overcome this issue is using appropriate reference tasks. To discriminate language areas among those involved in the completion of a given task, Binder has suggested using well-designed reference tasks. The idea is

<figure><img src="image_4.png" /></figure>
to remove the non-specific or non-lateralized activations of primary areas and/or executive regions by applying the difference paradigm (Binder 2011). Compared to a non-verbal reference, the use of a verbal reference tasks allows left hemisphere language areas to be specifically highlighted, as shown by Ferstl’s meta-analysis (Ferstl et al. 2008). The use of verbal reference tasks with functional magnetic resonance imaging (fMRI) has proven to successfully measure activation asymmetry, a proxy of language dominance strongly concordant with Wada testing. Note that this is true whether hemispheric or regional asymmetry of activations is used for the evaluation of language hemispheric dominance [review in Dym et al. (2011)].

Thus, asymmetry represents an additional method for increasing the specificity of identifying left hemisphere language areas. Typical language organization, seen in 90% of the healthy population (Mazoyer et al. 2014) and 97% of healthy right-handers (Zago et al. 2017), is characterized by a strong left hemisphere dominance, giving rise to regional leftward asymmetries in fMRI. Adding to the detection of activated areas (by comparison to a high-level verbal reference task), a criterion based on leftward asymmetry would certainly increase the specificity of identifying left hemisphere language areas.

Another difficulty in identifying essential language areas with functional imaging is the fact that different tasks lead to different patterns of activation. One way to overcome this difficulty is to combine several language tasks in the same participant and apply conjunction analyses to unravel the activated and asymmetrical regions independent of the type of task or modality involved (Papathanassiou et al. 2000; Jobard et al. 2007; Dodoo-Schittko et al. 2012).

Finally, the task-induced approach does not provide any information on how the different activated areas are organized. The co-activation of a group of regions does not indicate that they are all strongly functionally connected and thus constitute a network. Resting-state intrinsic connectivity has proven to be capable of identifying the organization of brain networks underpinning cognitive functions. A good illustration of this concept is provided by Turken and Dronkers, who conducted correlation analysis in resting-state images of healthy participants using the posterior middle temporal gyrus region as the seed (Turken and Dronkers 2011). In this work, this seed region was selected because its lesion was associated with strong comprehension deficits in aphasics, and its location was previously identified by probabilistic lesion mapping (Dronkers et al. 2004). Using this approach, Turken et al. revealed a network of areas connected at rest that support speech comprehension in healthy individuals. Investigating intrinsic connectivity would thus be an interesting means to investigate the networks existing at rest among the areas activated during language tasks. Connectivity measures can provide essential information on

how regions are connected and how they are organized in networks. Graph analysis methodology applied on resting-state connectivity also permits the measurement of the connectivity strength of each region with all other regions of a given network to which it belongs, thereby characterizing its role in the network. In particular, identifying the topological roles of the regions is possible, i.e., identifying hubs, regions essential to a given network and therefore essential to the cognitive function(s) they support (Sporns et al. 2007).

To propose an atlas of left hemisphere high-order language areas, we first combined multiple-language fMRI task-induced activation mapping and conjunction analysis to select a set of both activated and leftward asymmetrical areas during sentence processing. Second, we clustered the regions identified in the first step into networks based on their intrinsic connectivities at rest. Third, we applied graph analysis to characterize the roles of the regions in communication within and across networks. To this end, we utilized BIL&GIN, a database dedicated to the study of hemispheric specialization (Mazoyer et al. 2015), and selected 144 right-handers who were mapped during sentence production, reading and listening tasks compared to the production, reading and listening of lists of words, respectively. All but six participants were also mapped during the resting state. Most investigations of the resting-state and task-induced activation networks have relied on whole-brain comparisons between the functional connectivities measured in these two conditions [review in Wig (2017)], although they correspond to very different physiological states (Raichle and Mintun 2006; Raichle 2015). Here, we aimed to find resting-state markers of left hemisphere activation in discrete language areas to provide a comprehensive tool for further research on the inter-individual variability of language areas. Such markers of language activation are likely to be of interest for studies in which no task-induced activations are documented but instead include a resting-state acquisition. Consequently, we used homotopical regions of interest (hROIs) from the AICHA atlas, a functional atlas obtained from intrinsic connectivity analysis (Joliot et al. 2015). We used AICHA hROIs because (1) we needed an atlas suitable for functional imaging and, in that respect, AICHA, which was elaborated from resting-state connectivity, is optimal for analysing functional data; (2) AICHA has been specifically designed to identify functionally homotopic regions of interest, enabling the accurate computation of functional asymmetries since it avoids the potential bias that anatomical and functional areas do not strictly overlap. The different sets of hROIs corresponding to the supramodal sentence processing areas of the proposed atlas (SENSAAS) are available in the Montreal Neurological Institute (MNI) space at http://www.gin.cnrs.fr/en/tools/.
# Materials and methods

## Participants

From the BIL&GIN database, we selected 144 healthy right-handers (72 women) who completed the fMRI battery, including several language tasks (Mazoyer et al. 2015). The sample mean age was 27 years (SD = 6 years), and the women were two years younger than the men (women $26 \pm 5$; men: $28 \pm 7$, $p=0.053$). The mean educational level of the participants was 16 years (SD = 6 years), with no significant difference between the men and women ($p > 0.05$). All participants reported themselves as right-handed; their mean normalized finger tapping test asymmetry [((right number of taps – left number of taps)/(left + right number of taps)) × 100] was 6.25 (SD = 4.3), and their mean Edinburgh score was 93.5 (SD = 11), confirming their right-handedness. There was no difference between gender for the Edinburgh score ($p=0.47$), although there was a slightly stronger rightward manual laterality in women (finger tapping test asymmetry in women: $6.9 \pm 3.8$; men: $5.7 \pm 4.7$, $p=0.08$, controlling for age).

Of these participants, 138 (mean age 27 years (SD = 6 years), 68 women) also completed a resting-state fMRI (rs-fMRI) acquisition lasting 8 min. Note that this resting-state acquisition was performed on average 9 months (SD = 9.6 months) before the language task acquisition in all but five cases. In these five cases the resting-state acquisition occurred approximately 1 year after the language session [range (11.2–13.8) months].

## Image acquisition and processing

### Structural imaging

Structural images were acquired using the same 3T Philips Intera Achieva scanner including high-resolution T1-weighted volumes (sequence parameters: TR, 20 ms; TE, 4.6 ms; flip angle, 10°; inversion time, 800 ms; turbo field echo factor, 65; sense factor, 2; field of view, 256×256×180 mm³; isotropic voxel size 1×1×1 mm³). For each participant, the line between the anterior and posterior commissures was identified on a mid-sagittal section, and the T1-MRI volume was acquired after orienting the brain in the bi-commissural coordinate system. T2*-weighted multi-slice images were also acquired [T2*-weighted fast field echo (T2*-FFE), sequence parameters: TR=3,500 ms; TE=35 ms; flip angle=90°; sense factor=2; 70 axial slices; 2×2×2 mm³ isotropic voxel size].

### Task-induced functional imaging

**Training** To ensure proper task execution, the participants were trained outside the scanner in the hour preceding the

fMRI session. The training used stimuli that were of the same nature but different from those used during the fMRI session.

**Language tasks** Three runs were administered to the participants. They included a sentence task involving phonological, semantic, prosodic and syntactic processing and a word-list reference task, a less complex, albeit high-level, verbal task. To achieve homogeneity in the sentence task material, 51 line drawings illustrating the stories of ‘Le petit Nicolas’ (Little Nicholas), a classic French children’s series, were used. The three tasks consisted of a randomized alternation of event-related trials devoted to sentence processing, with event-related trials devoted to the verbal reference task, i.e., lists of words. The drawings used for the reference task were scrambled versions of the line drawings, and the stimuli presented either orally or visually were lists of months, days and/or seasons. Within each trial, the subject was shown either a line drawing or a scrambled drawing for 1 s, immediately followed by a central fixation cross-hair. While fixating the cross, the subject performed either the sentence task or the word reference task. Once the task was completed, a low-level reference task, detecting the transformation of a centrally displayed cross into a square, was presented. When the subjects detected this change, they were asked to press a button with their index finger of the assigned hand. The square was then displayed until the end of the trial. This second part of the trial, which lasted at least half of the total trial duration, aimed at refocusing the subject’s attention to a non-verbal stimulus and controlling for manual motor response activation, which was also present in the first part of the trial. A 12-s presentation of a fixation cross-hair preceded and followed the first and last trial. Note that except during the drawings display, the subjects were asked to keep fixating the cross, and the star and square were then presented on the centre of the screen.

*Sentence and list of word production tasks* During the production run, after seeing a Little Nicholas line drawing, the subject was instructed to covertly generate a sentence beginning with a subject (The little Nicholas..., The gentleman...) and a complement (with his satchel..., in shorts..., with glasses...), followed by a verb describing the action taking place and ending with an additional complement of a place (in the street..., in the playground..., on the beach...) or a manner (with happiness..., nastily...). When a scrambled drawing was displayed, the subject was asked to covertly generate the list of the months of the year. The production paradigm randomly alternated ten 18-s trials of sentence generation with ten 18-s trials of generating the list of months. The response time limit, indicated by the transformation of the cross in a star, was 9 s, including the 1-s drawing display. The entire experimental run lasted 6 min and 24 s. The mean sentence production time was 5617 ms
(SD = 935 ms), while the mean duration of word-list production was 5249 ms (SD = 1131 ms).

*Sentence and list of word-listening tasks* When a Little Nicholas line drawing was displayed, the subject was instructed to carefully listen to a sentence dealing with the line drawing and click at the end of the sentence. For the LISN, when a scrambled drawing was displayed, he/she was instructed to listen to the list of the months, days of the week and/or seasons and click at the end of the list.

The paradigm consisted of a randomized alternation of thirteen 14-s sentence-listening trials with thirteen 14-s list-listening trials. The mean durations of auditory presentation were 4371 ± 468 ms for the sentences and 4386 ± 484 ms for the lists. The entire experimental run lasted 6 min and 28 s. The reaction times after sentence and list listening were 387 ms (SD = 125 ms) and 478 ms (SD = 97 ms), respectively.

*Sentence and list of word-reading tasks* Like in the other two tasks, when a line drawing was displayed, the subject was instructed to read a sentence based on the line drawing. When a scrambled drawing was displayed, he/she was instructed to read the list of months, days of the week and/or seasons.

The paradigm consisted of a randomized alternation of thirteen 14-s sentence-reading trials with thirteen 14-s list-reading trials. The entire experimental run lasted 6 min and 28 s. The average time for reading sentences was 3729 ms (SD = 567 ms), while reading the lists of words required 4412 ms (SD = 602 ms).

**Debriefing the fMRI tasks** Right after the fMRI sessions, the participants were asked to rate the difficulty of the task on a five-point scale (1-easy to 5-very difficult) and answer some debriefing questions about how they accomplished the task.

The production task had the highest task difficulty score reported by the participants (2.73), while the reading and listening tasks had low scores (1.14 and 1.20, respectively). All participants were able to recollect the sentence they produced when presented with the corresponding drawing for at least 5 of 10 images (mean = 9.43 images, SD = 0.96), with the mean number of words per sentence being 12.4 (SD = 2).

## Functional image acquisition

The functional volumes were acquired as T2*-weighted echo-planar EPI images (TR = 2 s; TE = 35 ms; flip angle = 80°; 31 axial slices with a 240×240 mm² field of view and 3.75×3.75×3.75 mm³ isotropic voxel size). In the three runs, 192, 194 and 194 T2*-weighted volumes were acquired for the sentence production, listening and reading tasks, respectively.

## Resting-state functional imaging (rs-fMRI)

Spontaneous brain activity was monitored for 8 min (240 volumes) using the same imaging sequence (T2*-weighted echo-planar images) as that used for the language tasks. Immediately prior to rs-fMRI scanning, the participants were instructed to “keep their eyes closed, to relax, to refrain from moving, to stay awake and to let their thoughts come and go”.

## Image analysis

**Functional imaging analysis common to task-induced and resting-state acquisitions** For each participant, (1) the T2*-FFE volume was rigidly registered to the T1-MRI; (2) the T1-MRI was segmented into three brain tissue classes: grey matter, white matter and cerebrospinal fluid; and (3) the T1-MRI scans were normalized to the BIL&GIN template including 301 volunteers from the BIL&GIN database (aligned to the MNI space) using the SPM12 “normalise” procedure with otherwise default parameters.

For each of the three fMRI runs, data were corrected for slice timing differences. To correct for subject motion during the runs, all the T2*-weighted volumes were realigned using a six-parameter rigid-body registration. The EPI-BOLD scans were then registered rigidly to the structural T2*-FFE image. The combination of all registration matrices allowed for warping the EPI-BOLD functional scans to the standard space with a single trilinear interpolation.

**Specific task-induced functional imaging analysis** First, a 6-mm full width at half maximum (Gaussian filter) was applied to each run. Global linear modelling (statistical parametric mapping (SPM), http://www.fil.ion.ucl.ac.uk/spm/) was used for processing the task-related fMRI data. For each participant, BOLD variations corresponding to each sentence versus the list belonging to the same run were computed [sentence minus word-list production (PROD<sub>SENT-WORD</sub>), sentence minus word-list reading (READ<sub>SENT-WORD</sub>), and sentence minus word-list listening (LISN<sub>SENT-WORD</sub>)]. Finally, contrast maps (defined at the voxel level) were subjected to hROI analysis. BOLD signal variations were measured in 192 pairs of functionally defined hROIs of the AICHA atlas (Joliot et al. 2015) adapted to SPM12, excluding 7 hROIs pairs belonging to the orbital and inferior temporal parts of the brain in which signals were reduced due to susceptibility artefacts. For each participant, we computed contrast maps of the three language conditions. We then calculated the right and left hROI BOLD signal variations for each of the 185 remaining pairs by averaging the contrast BOLD values of all voxels located within the hROI volume.
**Specific analysis of resting-state functional images** Time series white matter and cerebrospinal fluid (individual average time series of voxels that belonged to each tissue class) and temporal linear trends were removed from the rs-fMRI data series using regression analysis. Additionally, rs-fMRI data were temporally filtered using a least squares linear-phase finite impulse response filter design bandpass (0.01–0.1 Hz).

For each participant and hROI (the same 185 homotopic ROIs as those used in the task-induced analysis), an individual BOLD rs-fMRI time series was computed by averaging the BOLD fMRI time series of all voxels located within the hROI volume.

## Part 1: Identification and characterization of hROIs exhibiting both leftward activation and leftward asymmetrical activation in all three tasks

To complete the identification of high-order language areas, we first searched for hROIs that were both significantly co-activated and significantly leftward asymmetrical on average among the 144 participants during the PROD<sub>SENT-WORD</sub>, READ<sub>SENT-WORD</sub>, and LISN<sub>SENT-WORD</sub> tasks.

### Statistical analysis

#### hROI selection

Using JMP14 (http://www.jmp.com, SAS Institute Inc., 2018), conjunction analysis was conducted to select the left-hemisphere hROIs exhibiting BOLD signal variations that were both significantly positive and significantly larger than that in their right counterparts in all three tasks. An hROI was selected whenever it was significantly activated in each of the three task contrasts using a significance threshold set to *p* < 0.05 per contrast. The significance threshold for the conjunction of activation in the three tasks was thus 0.05 × 0.05 × 0.05 = 1.25 × 10<sup>-4</sup>. The second criterion for hROI selection was the existence of a significant leftward asymmetry in each of the three task contrasts, the threshold of significance

of this second conjunction being again 1.25 × 10<sup>-4</sup>. Finally, since to be selected, a given hROI had to fulfil both criteria, the overall significance threshold for the conjunction of conjunction analyses was 1.5 × 10<sup>-8</sup> = (1.25 × 10<sup>-4</sup>)<sup>2</sup>.

### Results

#### hROI selection

Among the 80 hROIs jointly activated in the 3 contrasts (Table 1), 46 also showed joint asymmetries. In total, 32 hROIs showed both joint activation on the left and joint asymmetry (Fig. 1; Table 2).

On the lateral surface of the left frontal lobe, the regions having both joint activation and leftward asymmetry during the three language tasks covered the left inferior frontal gyrus (pars triangularis: F3t and pars opercularis: F3O1), the adjacent inferior frontal sulcus (f2_2), the junction of the middle frontal gyrus with the precentral sulcus (prec4), and the upper part of the precentral sulcus (prec3) located dorsally to prec4. The medial part of the superior frontal gyrus (F1_2), the upper paracentral gyrus (pCENT4), and the pre-superior motor areas (SMA2 and SMA3) were also part of these areas in the medial frontal lobe. Two hROIs were located within the anterior insula (INSa2 and INSa3), while the INSa1 hROI was located medially and ventrally close to the amygdala. On the lateral surface of the temporal lobe, the hROIs overlapped the entire length of the superior temporal sulcus (STS2, STS3 and STS4), extending to the temporal pole anteriorly (STS1), to the superior temporal gyrus dorsally (T1_4), to the supramarginal (SMG7) and angular gyri (AG2) posteriorly, crossing the middle temporal gyrus (T2_3 and T2_4) and joining the inferior temporal gyrus (T3_4), the inferior occipital gyrus (O3_1), and ventrally the fusiform gyrus (FUS4). Regions located within the hippocampus (HIPP2), parahippocampal gyrus (pHIPP1) and amygdala (AMYG) were also part of the selected areas. In the posterior medial wall, the dorsal part of the precuneus (PRECU6) together with the posterior cingulum (CINGp3) were selected using this approach. Sub-cortical areas jointly activated and leftward asymmetrical during the three tasks

<table><caption>**Table 1** Results of conjunction analyses across each sentence minus word-list contrasts for production (PROD<sub>SENT-WORD</sub>), listening (LISN<sub>SENT-WORD</sub>) and reading (READ<sub>SENT-WORD</sub>) tasks in terms of the number of hROIs</caption><thead><tr><th></th><th>L activation</th><th>L asymmetry</th><th>Conjunction of activation and asymmetry</th></tr></thead><tbody><tr><td>PROD<sub>SENT-WORD</sub></td><td>133</td><td>93</td><td>75</td></tr><tr><td>LISN<sub>SENT-WORD</sub></td><td>116</td><td>73</td><td>64</td></tr><tr><td>READ<sub>SENT-WORD</sub></td><td>97</td><td>60</td><td>43</td></tr><tr><td>Conjunctions of 3 contrasts</td><td>80</td><td>46</td><td>32</td></tr></tbody></table>

Numbers of hROIs with significant left activation, leftward asymmetry or conjunction of activation and asymmetry for the three “sentence minus word” contrasts
<figure><img src="image_5.png" /><figcaption>**Fig. 1** Locations of the 32 hROIs co-leftward activated and co-leftward asymmetrical during the completion of 3 sentence minus word-list tasks by 144 healthy right-handers and corresponding networks after hROI clustering based on resting-state connectivity. **a** Left lateral view of 3D surfaces rendering the 32 hROIs on the BIL&GIN display template in the MNI space with Surf Ice software (https://www.nitrc.org/projects/surface/). **b** Representation of hROIs on left</figcaption></figure>

hemisphere axial slices from the BIL&GIN display template; the hROI numbers correspond to the z-axis in the MNI space. **c** Lateral and medial views of the three identified networks. SENT_CORE network: red, SENT_MEM: light blue and SENT_VISU: green. Correspondences between the abbreviations and the full names of the AICHA atlas can be found in Table 2
<table><caption>**Table 2** Names and abbreviations of the 32 hROIs showing joint left activation and left asymmetry during the three sentences minus word-list contrasts for production (PROD<sub>SENT-WORD</sub>), listening</caption><thead><tr><th></th><th>AICHA hROI name</th><th>Abbreviation</th><th>Cluster</th><th>X (mm)</th><th>Y (mm)</th><th>Z (mm)</th></tr></thead><tbody><tr><td rowspan="8">Frontal and insula</td><td>S_Precentral-3</td><td>prec3</td><td>SENT_MEM</td><td>-18.2</td><td>-8.7</td><td>69.3</td></tr><tr><td>S_Precentral-4</td><td>prec4</td><td>SENT_CORE</td><td>-42.2</td><td>0.7</td><td>49.9</td></tr><tr><td>G_Frontal_Sup-2</td><td>F1_2</td><td>SENT_CORE</td><td>-11.9</td><td>46.5</td><td>41.4</td></tr><tr><td>S_Inf_Frontal-2</td><td>f2_2</td><td>SENT_CORE</td><td>-43.1</td><td>14.8</td><td>29.4</td></tr><tr><td>G_Frontal_Inf_Tri-1</td><td>F3t</td><td>SENT_CORE</td><td>-49.4</td><td>25.6</td><td>4.7</td></tr><tr><td>G_Frontal_Inf_Orb-1</td><td>F3O1</td><td>SENT_CORE</td><td>-42.2</td><td>30.5</td><td>-16.9</td></tr><tr><td>G_Insula-anterior-1</td><td>INSa1</td><td>SENT_VISU</td><td>-20.3</td><td>5.0</td><td>-19.3</td></tr><tr><td>G_Insula-anterior-2</td><td>INSa2</td><td>SENT_CORE</td><td>-33.8</td><td>16.8</td><td>-12.7</td></tr><tr><td rowspan="10">Temporal and parietal</td><td>G_Insula-anterior-3</td><td>INSa3</td><td>SENT_CORE</td><td>-33.7</td><td>23.7</td><td>0.6</td></tr><tr><td>G_Temporal_Sup-4</td><td>T1_4</td><td>SENT_CORE</td><td>-58.7</td><td>-23.3</td><td>3.7</td></tr><tr><td>G_Temporal_Mid-3</td><td>T2_3</td><td>SENT_CORE</td><td>-61.0</td><td>-35.0</td><td>-4.8</td></tr><tr><td>G_Temporal_Mid-4</td><td>T2_4</td><td>SENT_CORE</td><td>-53.1</td><td>-59.4</td><td>7.0</td></tr><tr><td>G_Temporal_Inf-4</td><td>T3_4</td><td>SENT_VISU</td><td>-50.0</td><td>-60.6</td><td>-7.6</td></tr><tr><td>S_Sup_Temporal-1</td><td>STS1</td><td>SENT_CORE</td><td>-49.7</td><td>14.0</td><td>-21.5</td></tr><tr><td>S_Sup_Temporal-2</td><td>STS2</td><td>SENT_CORE</td><td>-54.9</td><td>-7.2</td><td>-12.8</td></tr><tr><td>S_Sup_Temporal-3</td><td>STS3</td><td>SENT_CORE</td><td>-54.7</td><td>-33.0</td><td>-1.7</td></tr><tr><td>S_Sup_Temporal-4</td><td>STS4</td><td>SENT_CORE</td><td>-56.5</td><td>-48.4</td><td>13.4</td></tr><tr><td>G_SupraMarginal-7</td><td>SMG7</td><td>SENT_CORE</td><td>-55.2</td><td>-51.7</td><td>25.5</td></tr><tr><td>G_Angular-2</td><td>AG2</td><td>SENT_CORE</td><td>-37.5</td><td>-70.4</td><td>39.5</td></tr><tr><td rowspan="7">Internal surface</td><td>G_Occipital_Inf-1</td><td>O3_1</td><td>SENT_VISU</td><td>-48.4</td><td>-69.0</td><td>-4.3</td></tr><tr><td>G_Fusiform-4</td><td>FUS4</td><td>SENT_VISU</td><td>-43.1</td><td>-49.8</td><td>-17.4</td></tr><tr><td>G_ParaHippocampal-1</td><td>pHIPP1</td><td>SENT_VISU</td><td>-15.7</td><td>-4.0</td><td>-18.4</td></tr><tr><td>G_Hippocampus-2</td><td>HIPP2</td><td>SENT_MEM</td><td>-24.9</td><td>-32.5</td><td>-2.7</td></tr><tr><td>G_Supp_Motor_Area-2</td><td>SMA2</td><td>SENT_CORE</td><td>-10.6</td><td>18.2</td><td>63.1</td></tr><tr><td>G_Supp_Motor_Area-3</td><td>SMA3</td><td>SENT_CORE</td><td>-7.2</td><td>7.6</td><td>65.6</td></tr><tr><td>G_Paracentral_Lobule-4</td><td>pCENT4</td><td>SENT_MEM</td><td>-6.4</td><td>-29.2</td><td>75.9</td></tr><tr><td rowspan="4">Sub-cortical</td><td>G_Cingulum_Pos-3</td><td>CINGp3</td><td>SENT_MEM</td><td>-5.1</td><td>-42.9</td><td>10.0</td></tr><tr><td>G_Precuneus-6</td><td>PRECU6</td><td>SENT_MEM</td><td>-7.4</td><td>-61.3</td><td>64.1</td></tr><tr><td>N_Amygdala-1</td><td>AMYG</td><td>SENT_VISU</td><td>-21.9</td><td>-0.4</td><td>-11.5</td></tr><tr><td>N_Thalamus-4</td><td>THA4</td><td>SENT_MEM</td><td>-3.2</td><td>-14.4</td><td>8.4</td></tr><tr><td rowspan="2">Sub-cortical</td><td>N_Putamen-2</td><td>PUT2</td><td>SENT_VISU</td><td>-23.3</td><td>6.3</td><td>0.8</td></tr><tr><td>N_Putamen-3</td><td>PUT3</td><td>SENT_VISU</td><td>-28.0</td><td>-6.3</td><td>1.8</td></tr></tbody></table>

(LISN<sub>SENT-WORD</sub>) and reading (READ<sub>SENT-WORD</sub>) tasks; the network label to which they were clustered; and their coordinates in MNI space after SPM12 normalization of the AICHA atlas

covered almost the entire putamen (PUT2 and PUT3) and a thalamic hROI located medially (THA4).

## Part 2: Identification of networks based on the resting-state connectivity matrix of the 32 hROIs co-activated and co-leftward asymmetrical during the 3 sentence minus word-list tasks

In a second step, we investigated the intrinsic functional organization of the 32 hROIs selected in the first step. We computed the intrinsic connectivity matrix between these 32 hROIs for the subsample of 138 right-handed participants

who completed a resting-state acquisition. We completed a hierarchical clustering analysis of this intrinsic connectivity matrix to identify temporally coherent networks within this set of hROIs.

## Methods

### Calculation of the intrinsic connectivity matrix

An intrinsic connectivity matrix was calculated for each of the 138 individuals and for each of the 496 possible pairs of hROIs (N × (N − 1))/2, with N = 32). The intrinsic connectivity matrix off-diagonal elements were the Pearson correlation coefficients between the rs-fMRI time series of
the hROI pairs. The intrinsic connectivity matrix diagonal elements were set to zero because no information on the correlation for a specific hROI with itself exists (Rubinov and Sporns 2010). The individual intrinsic connectivity matrix was then Fisher z-transformed before being averaged over the subsample of 138 individuals, thereby producing a mean intrinsic connectivity matrix.

## Identification and characterization of networks

The agglomerative hierarchical cluster analysis (aHCA) method was applied to extract brain networks from this mean intrinsic connectivity matrix. We first transformed the Pearson correlation ($r_{i,j}$) between hROI i and hROI j into a distance ($d_{i,j}$) using the equation $d_{i,j} = (1 - r_{i,j})/2$, as in Doucet et al. (2011), resulting in a 32×32 dissimilarity matrix. According to Lance and Williams (1979), the previous equation is “unequivocally the best” to transform a correlation into a distance on real data sets. Finally, we used an agglomerative hierarchical clustering algorithm [aHCA, (Sneath and Sokal 1973)] for clustering and the Ward distance (Ward 1963) for aggregating the different hROIs into clusters. The number of clusters (networks) was determined using the R library NbClust (Charrad et al. 2014). This package provides 30 statistical indices for determining the optimal number of clusters and proposes the best clustering scheme from the different results obtained by varying all combinations of number of clusters for the chosen method, in this case, aHCA with Ward’s distance. We chose the number of clusters that fulfilled a maximum of indices.

To characterize each network, we calculated its mean volume activation for each task contrast as the sum of the activations of all hROIs composing the network weighted by their individual volume and then divided by the sum of their volumes. The same computation was performed for the right hemisphere equivalent of each network, which was then used for computing the left-minus-right asymmetry of each network activation. We then compared activation amplitude and asymmetry values across networks and across tasks using a mixed-model ANOVA.

## Reliability of the network identification across individuals

We used multiscale bootstrap resampling (Efron et al. 1996) to assess the reliability of the identification of each cluster. In total, 10,000 multiscale bootstrap resampling datasets, including 50–140% of sample data from the 138 participants, were processed. Applying the R package “pvclust” (Suzuki and Shimodaira 2006) function to the multiscale bootstrap resampling outputs, we measured the approximately unbiased (AU) *p* value for each cluster. The AU *p* value for a network, the probability of this network occurring among the 138 participants, indicates the network’s reliability.

## Robustness of the networks identification with respect to the clustering method

We also assessed the robustness of the clustering method by comparing its output to those of three other clustering methods: aHCA with the average distance method (instead of Ward’s), Gaussian mixture model, and k-means (see supplementary material). Gaussian mixture modelling was conducted with the “Rmixmod” package with Normalized Entropy Criterion to find well-separated clusters and with a Gaussian model with diagonal variance matrices (Lebret et al. 2015).

We then compared the four different partitions through the adjusted Rand index (Hubert and Arabie 1985) allowing to get a similarity measure between two different classifications, an adjusted Rand index of 1 indicating similar partitions.

## Temporal correlation across networks and significance

To compute the mean intrinsic functional correlations between two networks, we used the same methodology as that used to compute the mean intrinsic connectivity matrix (see above).

First, for each individual and for each network, we computed the corresponding rs-fMRI time series by averaging the individual resting time series of all voxels of all hROIs belonging to this network.

Then, for each individual, we computed the Pearson correlation coefficients between all pairs of networks that we further Fisher z-transformed. Finally, each of these z-transformed coefficients was averaged across the sample of 138 individuals, providing a mean intrinsic functional correlation (*r*) for each pair of networks. We assessed the significance of each of these mean intrinsic functional correlations compared to 0 using a non-parametric sign test at the 0.05 significance level (Bonferroni correction for the number of network pairs).

## Results

### Identification and characterization of networks

Hierarchical clustering analysis revealed three networks from the selected set of 32 hROIs (Fig. 2).

**SENT_CORE network** The first network (Fig. 1c, red), termed SENT_CORE, was composed of 18 hROIs and was the most distant from the 2 others in terms of inertia. SENT_CORE included all lateral and medial hROIs of the frontal lobe, apart prec3, pCENT4 and anterior insula INSa2 and INSa3. SENT_CORE also included all temporal and parietal hROIs
of the lateral surface, except T3_4, which was aggregated with the network gathering visual hROIs.

We named this network SENT_CORE because it included essential sentence processing regions, as further described below. SENT_CORE was the largest network in terms of volume, as it was 9.2 times larger than SENT_MEM and 3.4 times larger than SENT_VISU (Table 3).

**SENT_VISU network** This group of clusterized areas included the inferior temporal and occipital gyri laterally (T3_4, O3_1); the mid-fusiform (FUS4) ventrally; the parahippocampal region; the amygdala (AMYG) and INSA1, close to the amygdala medially (Fig. 1c, green); and the two hROIs of the putamen (PUT2 and PUT3). We labelled it

SENT_VISU because it aggregated four hROIs acknowledged as involved in visual processing.

**SENT_MEM network** This network (Fig. 1c, light blue) included three regions of the medial wall—the paracentral gyrus (pCENT4), the precuneus (PRECU6) and the posterior cingulate (CINGp3)—and the posterior part of the hippocampus (HIPP2) as well as one frontal area at the upper end of the precentral sulcus (prec3) and the THA4 hROI located medially in the thalamus. We named it SENT_MEM because these posterior areas belong to both the posterior regions of the DMN involved in the posterior hippocampus in episodic memory, as further discussed below.

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 2** Results of the agglomerative hierarchical cluster analysis method. **a** Dendrogram of aHCA of the mean intrinsic connectivity matrix (SENT_CORE network: red, SENT_MEM: light blue, SENT_VISU: green). Approximately unbiased *p* values are indicated for</figcaption></figure>

each identified network. **b** Scree plot of aHCA of the mean intrinsic connectivity matrix. For both graphs, the red horizontal line corresponds to the threshold applied to select the number of networks

<table><caption>**Table 3** Mean volumetric activation (and standard deviation) of the four language networks in each sentence minus word-list contrast for production (PROD<sub>SENT-WORD</sub>), listening (LISN<sub>SENT-WORD</sub>) and reading (READ<sub>SENT-WORD</sub>) tasks in 144 healthy right-handers</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Volume (mm<sup>3</sup>)</th><th colspan="3">Left activation</th><th colspan="3">Leftward asymmetry</th></tr><tr><th>PROD<sub>SENT-WORD</sub></th><th>LISN<sub>SENT-WORD</sub></th><th>READ<sub>SENT-WORD</sub></th><th>PROD<sub>SENT-WORD</sub></th><th>LISN<sub>SENT-WORD</sub></th><th>READ<sub>SENT-WORD</sub></th></tr></thead><tbody><tr><td>SENT_CORE</td><td>83,232</td><td>0.73±0.31</td><td>0.43±0.20</td><td>0.55±0.26</td><td>0.41±0.22</td><td>0.25±0.14</td><td>0.28±0.19</td></tr><tr><td>SENT_MEM</td><td>9024</td><td>0.40±0.36</td><td>0.28±0.23</td><td>0.21±0.26</td><td>0.11±0.17</td><td>0.08±0.12</td><td>0.07±0.14</td></tr><tr><td>SENT_VISU</td><td>24,368</td><td>0.37±0.27</td><td>0.30±0.17</td><td>0.25±0.21</td><td>0.16±0.13</td><td>0.11±0.08</td><td>0.12±0.10</td></tr></tbody></table>

The mean volumetric activation for a network was calculated from the sum of the activations of the hROIs comprising the network weighted by their individual volumes and then divided by the volume of the network
## Profile comparisons between networks

In terms of mean voluminal activity (Table 3), ANOVA revealed a significant network effect, a task effect and a network × task interaction (all post hoc Tukey’s HSD: $p < 10^{-4}$).

The interaction occurred because while in `SENT_CORE`, there was greater activation in `PROD`, followed by `READ` and then by `LISN` (all $p < 10^{-4}$), `SENT_MEM` and `SENT_VISU` had different profiles. Although the greatest activation values were also observed during `PROD` (all $p < 0.003$), activation values were significantly higher during `LISN` than during `READ` in both networks (`SENT_VISU: $p = 0.01$; `SENT_MEM: $p = 0.0067`).

ANOVA on asymmetries also revealed a significant network effect, a task effect and a network × task interaction (all post hoc Tukey’s HSD: $p < 10^{-4}$). The interaction was due to different profiles of asymmetry in `SENT_CORE` than in the two other networks. In `SENT_CORE`, the profile of asymmetries was the same as that of activation: larger in `PROD`, followed by `READ` and then by `LISN` (all post hoc Tukey’s HSD: $p < 10^{-4}$). In `SENT_VISU`, the asymmetry during `PROD` was slightly larger than that during `READ` (`SENT_VISU: $p = 0.0005$; `SENT_MEM: $p = 0.01$`) and larger than that during `LISN` ($p = 0.0001$). In addition, there was no difference in asymmetry between `PROD` and `LISN` in `SENT_MEM` ($p = 0.075$) and no difference in asymmetry between `READ` and `LISN` in `SENT_VISU` ($p = 0.52$) or `SENT_MEM` ($p = 0.25$). The task main effect was related to larger asymmetries in `PROD` and the network main effect to larger asymmetries in `SENT_CORE` (all $p < 0.0001$).

## Assessing the reliability of the identification of 3 networks across individuals in the 138 participants

The AU $p$ values provided by the multiscale bootstrap resampling method showed that each network of the first partition was reliable at levels of 93%, corresponding to `SENT_CORE` on one side and to the `SENT_VISU` and `SENT_MEM` on the other side (Fig. 2). However, for the second partition, `SENT_VISU` and `SENT_MEM` were reliable at 73% and 61%, respectively, indicating a lower reliability.

## Robustness of the identified networks with respect to the clustering method

The `SENT_CORE` network was identified by all 4 clustering methods, including at least 13 of the 18 hROIs initially found with the aHCA method using Ward’s distance (see supplementary Table 1). The aHCA method using the average distance metric led to an adjusted Rand index of 1, indicating a clustering similar to that achieved using the aHCA method with Ward’s distance. Comparing the Gaussian

Mixture Model and aHCA methods led to an adjusted Rand index of 0.76, indicating two highly similar partitions. Comparing the aHCA and $k$-means methods led to the lowest adjusted Rand index of 0.43.

Only one hROI of `SENT_CORE` (AG2) was segregated in `SENT_MEM` by Gaussian mixture modelling and by $k$-means, while all 17 other hROIs were classified together by all clustering methods but $k$-means (supplementary Table 1), although `PUT3` joined `SENT_CORE` according to Gaussian mixture modelling. $k$-means classified `INSa3`, `T2-4`, `T1-4` and `f2-2` in `SENT_VISU` rather than in `SENT_CORE`.

Only `prec3` was classified in `SENT_VISU` rather than `SENT_MEM` with both Gaussian mixture modelling and $k$-means, while `pHIPP1` shifted from `SENT_VISU` to `SENT_MEM` only with $k$-means.

## Temporal correlation across networks and significance

The chord diagram shown in Fig. 3 describes the average correlations between each pair of hROIs in the three networks. Strong and highly significant negative mean intrinsic correlations were found between `SENT_CORE` and `SENT_MEM` ($R = -0.27$; 92.03% of the participants showed a negative correlation, $p < 10^{-4}$), and a positive correlation was present between `SENT_MEM` and `SENT_VISU` ($R = 0.058$, 62.32% of the participants showed a positive correlation, $p = 0.0024$), while there was no significant correlation between `SENT_CORE` and `SENT_VISU` ($R = -0.037$; 56.52% of the participants showed a positive correlation, $p = 0.074$).

# Part 3: Graph theory analysis of the `SENT_CORE` network

We applied graph theory analysis to the `SENT_CORE` hROI pairwise correlation matrix, including only positive correlations since the inclusion of negative correlations in graph theory analysis remains controversial (Rubinov and Sporns 2010). Note that the graph theory analysis of intra-network communication was completed for only `SENT_CORE`, as the other two networks had too few nodes.

## Statistical analysis

### Identification of hubs using graph analysis metrics of the networks

**Measurements of weighted centrality** We measured the degree centrality (DC) of the hROIs composing the `SENT_`
<figure><img src="image_8.png" /><figcaption>**Fig. 3** Chord diagram of the temporal correlation across each hROI composing the three networks averaged in the whole group. Abbreviations for hROIs of the AICHA atlas can be found in Table 2 (colour scale goes from red for positive correlation to blue for negative correlations, and the line width indicates the strength of the correlation)</figcaption></figure>

CORE network, corresponding to the sum of the strength of the positive correlation of each node (hROI). DC can thus be interpreted as the amount of information that a given hROI receives from the hROIs to which it is directly connected, i.e., the DC measures the importance of a given hROI within its network according to the number and strength of interactions it undergoes with the other hROIs.

The betweenness centrality (BC) was also measured for SENT_CORE as defined by Opsahl et al. (2010). The BC of an hROI can be interpreted as the participation rate of that hROI in the set of shortest paths between any pair of nodes within the network; i.e., BC measures the dependence of the network on a specific hROI for its communication.

**Hub definition and clustering** To discriminate hubs among SENT_CORE, we applied a combination of Sporns et al. (2007) and van den Heuvel et al. (2010) definitions. We considered that an hROI had the properties of a hub when its DC and BC values were larger than the means plus one standard deviation of the DC and BC values of the hROI set in the network.

To assess whether the hubs identified in SENT_CORE participated in communication with the other two networks or whether its communication was only intra-SENT_CORE,

we calculated the participation index (pIndex) criteria as defined by Guimera (Guimera and Amaral 2005). hROIs having the 15% highest pIndex values were considered connector hubs (i.e., between networks) (van den Heuvel and Sporns 2011), while the other hROIs corresponded to provincial hubs, i.e., an hROI communicating with only its own network.

**Investigation of the relationship between intrinsic connectivity and activation measured during the language tasks**

**Relationships between DC and BOLD variation at the hROI level** We investigated whether an hROI exhibiting high intrinsic connectivity with other areas of the SENT_CORE was more activated during language tasks. For this, we performed a MANCOVA with repeated measures with a TASK main effect, a DC main effect and a DC by TASK interaction. Correlations values and corresponding *p* values between the DC and activation were computed for each hROI and each task.

To test the specificity of this relationship, we completed similar MANCOVA with DC measurements obtained for the 185 hROIs of the AICHA atlas covering the entire
left hemisphere. Therefore, the DC values were computed considering the connections of hROIs belonging to SENT_CORE with all other hROIs of the left hemisphere. This analysis made it possible to more deeply characterize whether the relationship between the DC and activation during the language tasks was specific to the essential language network intrinsic connectivity or whether this relationship was held at the hemispheric level.

**Relationships between the DC measured in the SENT_CORE network and BOLD variation upon pooling all hROIs and participants** To test whether the relationship previously identified between DC and BOLD signal variation for each hROI of SENT_CORE was a general property that could be extrapolated to any hROI of any participant, we applied the method proposed by Buckner for evaluating the relationship between DC and beta-amyloid accumulation in Alzheimer's disease (Buckner et al. 2009). Correlation coefficients obtained for the three tasks between DC values and BOLD variations were compared using the R package “cocor” (Diedenhofen and Musch 2015) to determine whether there was any difference across language tasks.

# Results

## Graph analysis of SENT_CORE

**Sample distributions of DC and BC values** DC variation across the hROIs spanned from 2.24 to 6.48 (Table 4), and the DC standard deviation was very consistent across hROIs ranging from 1.17 to 1.79.

By contrast, BC variation across the hROIs spanned from 0.72 to 18.98 (Table 4). Notably, only 3 hROIs had low numbers of BC null values across the sample of 138 participants: F3t, STS3 and STS4 (1%, 6% and 2% null values, respectively, Table 4).

**Hub identification and characterization** Three hROIs corresponded to the hub definition, i.e., BC and DC values above the chosen significance thresholds (mean+SD) of 10.26 and 5.55, respectively (Table 4). The first hROI (F3t) was located in the frontal lobe, and the other two were located in the posterior third of the STS (STS3 and STS4). The BC values of these 3 hubs were over 11, and their DC values were over 5.5, with F3t having the strongest values (Fig. 4; Table 4). Note that no other hROI exhibited a supra-threshold value for any of the centrality indices.

<table><caption>**Table 4** Betweenness and degree centrality of SENT_CORE hROIs</caption><thead><tr><th rowspan="2"></th><th rowspan="2">AICHA hROI</th><th colspan="3">Betweenness centrality</th><th colspan="5">Degree centrality</th></tr><tr><th>Mean</th><th>SD</th><th>% Null values</th><th>Mean</th><th>SD</th><th>Skewness</th><th>Kurtosis</th><th>p norm</th></tr></thead><tbody><tr><td rowspan="4">Frontal and insula</td><td>prec4</td><td>6.88</td><td>7.00</td><td>11</td><td>4.94</td><td>1.50</td><td>-0.07</td><td>-0.24</td><td>0.87</td></tr><tr><td>F1_2</td><td>6.30</td><td>5.42</td><td>17</td><td>4.64</td><td>1.40</td><td>0.28</td><td>-0.32</td><td>0.27</td></tr><tr><td>f2_2</td><td>3.33</td><td>4.45</td><td>37</td><td>3.35</td><td>1.53</td><td>0.60</td><td>-0.42</td><td>0.0002</td></tr><tr><td><b>F3t</b></td><td><b>18.98</b></td><td><b>9.77</b></td><td><b>1</b></td><td><b>6.48</b></td><td><b>1.26</b></td><td><b>-0.06</b></td><td><b>-0.29</b></td><td><b>0.78</b></td></tr><tr><td rowspan="4">Temporal and parietal</td><td>F3O1</td><td>2.18</td><td>4.04</td><td>54</td><td>4.31</td><td>1.50</td><td>0.50</td><td>-0.02</td><td>0.018</td></tr><tr><td>INSa2</td><td>2.49</td><td>3.66</td><td>43</td><td>3.84</td><td>1.46</td><td>0.27</td><td>-0.23</td><td>0.14</td></tr><tr><td>INSa3</td><td>2.54</td><td>3.98</td><td>40</td><td>3.01</td><td>1.29</td><td>0.76</td><td>0.63</td><td>0.0005</td></tr><tr><td>T1_4</td><td>0.75</td><td>2.19</td><td>77</td><td>3.12</td><td>1.47</td><td>0.56</td><td>-0.20</td><td>0.0012</td></tr><tr><td rowspan="4">Internal surface</td><td>T2_3</td><td>5.28</td><td>5.96</td><td>29</td><td>5.44</td><td>1.43</td><td>0.05</td><td>-0.17</td><td>0.71</td></tr><tr><td>T2_4</td><td>2.21</td><td>4.28</td><td>58</td><td>3.46</td><td>1.79</td><td>0.24</td><td>-0.60</td><td>0.074</td></tr><tr><td>STS1</td><td>2.46</td><td>3.57</td><td>41</td><td>4.55</td><td>1.52</td><td>-0.10</td><td>-0.18</td><td>0.44</td></tr><tr><td>STS2</td><td>3.07</td><td>4.85</td><td>41</td><td>4.30</td><td>1.44</td><td>0.14</td><td>-0.58</td><td>0.11</td></tr><tr><td rowspan="2"></td><td><b>STS3</b></td><td><b>13.04</b></td><td><b>9.66</b></td><td><b>6</b></td><td><b>6.32</b></td><td><b>1.36</b></td><td><b>-0.36</b></td><td><b>0.23</b></td><td><b>0.088</b></td></tr><tr><td><b>STS4</b></td><td><b>13.33</b></td><td><b>8.25</b></td><td><b>2</b></td><td><b>5.58</b></td><td><b>1.51</b></td><td><b>0.09</b></td><td><b>-0.22</b></td><td><b>0.95</b></td></tr><tr><td rowspan="2"></td><td>SMG7</td><td>6.23</td><td>6.86</td><td>22</td><td>5.28</td><td>1.60</td><td>-0.21</td><td>-0.64</td><td>0.037</td></tr><tr><td>AG2</td><td>0.72</td><td>1.65</td><td>71</td><td>2.24</td><td>1.17</td><td>0.88</td><td>1.11</td><td>0.0002</td></tr><tr><td rowspan="2"></td><td>SMA2</td><td>2.21</td><td>3.58</td><td>46</td><td>4.35</td><td>1.33</td><td>0.33</td><td>-0.09</td><td>0.32</td></tr><tr><td>SMA3</td><td>2.24</td><td>3.60</td><td>49</td><td>4.22</td><td>1.28</td><td>0.05</td><td>-0.36</td><td>0.45</td></tr></tbody></table>

The means and standard deviations (SD) of the betweenness centrality (BC) and degree centrality (DC) were computed by averaging the BC and DC values of each participant for each SENT_CORE hROI. For BC, the percentage of null values is based on the number of BC values at zero among the 138 subjects for one hROI. For DC, the skewness, kurtosis and Shapiro-Wilk normality test (p norm) correspond to information regarding the normality of the DC distribution for each hROI. A value above 0.05 for the Shapiro-Wilk normality test indicates that the DC was normally distributed. hROIs in bold are those we found to be hubs
Concerning the pIndex, hubs were defined as the top 15% of the highest index (pIndex ranging from 0.587 to 0.989). Five hROIs were thus defined as connector hubs: T2_3 (pIndex = 0.989), F3t (pIndex = 0.987), F1_2 (pIndex = 0.984), STS3 (pIndex = 0.983) and SMA2 (pIndex = 0.983).

Note that the centrality hubs F3t and STS3 were also connector hubs, meaning that they are important for both communication among the three different networks and for communication within the SENT_CORE network.

Note that T2_3 was a connector hub characterized by high DC and BC values (Table 4), although it did not meet the criteria to be labelled as a centrality hub (Fig. 4).

## Relationship between the DC at rest and activations during the language tasks in the SENT_CORE network

**Relationship at the individual hROI level** Using DC values computed from only SENT_CORE hROIs, we observed significant positive correlations between activations during the 3 language tasks and these DC values in 12 hROIs among the 18 constituting SENT_CORE together with a trend for PREC4 and f2_2 (Table 5). Among these hROIs, DC values were positively correlated with activations during the language tasks and the *R* value varied between 0.17 and 0.33. Moreover, in 8 of these 12 hROIs, there was no DC by Task interaction, meaning that the correlation between the DC and activation did not differ between the tasks. In the f2_2, INSa2, INSa3, T1_4 and SMG7 hROIs, a significant DC by Task interaction was observed. In f2_2, the interaction was due to non-significant correlation for the production task contrast, while the correlation was strong and significant for the reading and listening task contrasts. In INSa2, INSa3 and T1_4, the interaction was due to non-significant correlation for the listening task, while there were strong and significant correlations for production and reading. In

SMG7, the interactions were due to a lower correlation during listening (Table 5).

The results obtained using DC values computed from the entire set of 185 left hemisphere hROIs were strikingly different. There was a significant main effect of the DC in only 2 hROIs (F3t and SMG7, see supplementary Table 2), meaning that, except for these two regions, the strength of the correlation when the DC was calculated across the hROIs of the entire hemisphere did not explain the activation variations in SENT_CORE hROIs.

**Relationship at the global level using all participants and hROIs** There was a significant correlation between the DC values and BOLD variations measured in each of the 3 tasks when considering the 18 SENT_CORE hROIs and the 138 participants in a single analysis (Fig. 5) for each task. The coefficient correlation values were 0.158, 0.216 and 0.294 for sentence production, sentence listening and sentence reading, respectively. The correlation for reading was significantly larger than that for both listening (p=0.0025) and production (p=0.0075), and the latter two were not significantly different (p=0.80).

## Summary of results

Conjunction analysis of left-activated and leftward asymmetrical hROIs in 144 right-handed participants performing three language tasks (PROD<sub>SENT-WORD</sub>, READ<sub>SENT-WORD</sub> and LISN<sub>SENT-WORD</sub>) uncovered a set of 32 supramodal regions involved in lexico-syntactic processing. The hierarchical bottom-up clustering of the intrinsic connectivity between these 32 hROIs led to the identification of 3 networks, including a network of essential language areas (SENT_CORE) with strong positive correlations at rest across its 18 hROIs in more than 90% of the participants. The two other identified

<figure><img src="image_9.png" /><figcaption>**Fig. 4** Plot of degree centrality (DC) versus betweenness centrality (BC) in SENT_CORE. The mean plus standard deviation values of DC and BC define the quadrants. hROIs located in the superior right quadrant are hubs. Abbreviations for the hROIs of the AICHA atlas can be found in Table 2</figcaption></figure>

<figure><img src="image_10.png" /></figure>
<table><caption>**Table 5** Correlation analysis between the degree centrality measured in the SENT_CORE network and the mean activation in each of the three language tasks</caption><thead><tr><th rowspan="2"></th><th rowspan="2">AICHA hROI</th><th colspan="2">MANCOVA</th><th colspan="2">PROD<sub>SENT-WORD</sub></th><th colspan="2">LISN<sub>SENT-WORD</sub></th><th colspan="2">READ<sub>SENT-WORD</sub></th></tr><tr><th>DC</th><th>DC*Task</th><th>R</th><th>p</th><th>R</th><th>p</th><th>R</th><th>p</th></tr></thead><tbody><tr><td rowspan="7">Frontal and insula</td><td>prec4</td><td>0.072</td><td>0.84</td><td>0.066</td><td>0.44</td><td>0.163</td><td>0.056</td><td>0.146</td><td>0.088</td></tr><tr><td>F1_2*</td><td>0.0065</td><td>0.24</td><td>0.181</td><td>0.033</td><td>0.116</td><td>0.18</td><td>0.206</td><td>0.016</td></tr><tr><td>f2_2</td><td>0.092</td><td>0.0008</td><td>-0.115</td><td>0.18</td><td>0.229</td><td>0.0068</td><td>0.248</td><td>0.0034</td></tr><tr><td>F3t*</td><td>&lt;0.0001</td><td>0.35</td><td>0.264</td><td>0.0018</td><td>0.255</td><td>0.0025</td><td>0.245</td><td>0.0037</td></tr><tr><td>F3O1*</td><td>0.02</td><td>0.75</td><td>0.161</td><td>0.060</td><td>0.196</td><td>0.022</td><td>0.113</td><td>0.19</td></tr><tr><td>INSa2*</td><td>0.0025</td><td>0.032</td><td>0.200</td><td>0.019</td><td>0.072</td><td>0.40</td><td>0.267</td><td>0.0015</td></tr><tr><td>INSa3*</td><td>0.0006</td><td>0.017</td><td>0.231</td><td>0.0063</td><td>0.037</td><td>0.67</td><td>0.289</td><td>0.0006</td></tr><tr><td rowspan="10">Temporal and parietal</td><td>T1_4*</td><td>0.0006</td><td>0.0008</td><td>0.316</td><td>0.0002</td><td>0.026</td><td>0.77</td><td>0.260</td><td>0.0021</td></tr><tr><td>T2_3*</td><td>0.009</td><td>0.06</td><td>0.266</td><td>0.0016</td><td>0.173</td><td>0.043</td><td>0.210</td><td>0.013</td></tr><tr><td>T2_4*</td><td>0.0008</td><td>0.21</td><td>0.124</td><td>0.15</td><td>0.183</td><td>0.032</td><td>0.261</td><td>0.0020</td></tr><tr><td>STS1*</td><td>0.0002</td><td>0.27</td><td>0.289</td><td>0.00060</td><td>0.228</td><td>0.0073</td><td>0.231</td><td>0.0064</td></tr><tr><td>STS2</td><td>0.67</td><td>0.33</td><td>0.084</td><td>0.33</td><td>-0.028</td><td>0.75</td><td>0.033</td><td>0.70</td></tr><tr><td>STS3*</td><td>0.021</td><td>0.36</td><td>0.198</td><td>0.020</td><td>0.153</td><td>0.073</td><td>0.140</td><td>0.10</td></tr><tr><td>STS4*</td><td>0.0010</td><td>0.08</td><td>0.250</td><td>0.0032</td><td>0.187</td><td>0.028</td><td>0.283</td><td>0.0008</td></tr><tr><td>SMG7*</td><td>&lt;0.0001</td><td>0.007</td><td>0.329</td><td>0.0001</td><td>0.204</td><td>0.016</td><td>0.303</td><td>0.0003</td></tr><tr><td>AG2</td><td>0.41</td><td>0.78</td><td>0.009</td><td>0.92</td><td>0.091</td><td>0.29</td><td>0.065</td><td>0.45</td></tr><tr><td rowspan="2">Internal surface</td><td>SMA2</td><td>0.28</td><td>0.058</td><td>0.095</td><td>0.27</td><td>-0.061</td><td>0.48</td><td>0.149</td><td>0.082</td></tr><tr><td>SMA3</td><td>0.16</td><td>0.66</td><td>0.111</td><td>0.20</td><td>0.074</td><td>0.39</td><td>0.076</td><td>0.37</td></tr></tbody></table>

Correlations (R) were calculated within each hROI of the left hemisphere constituting the SENT_CORE network, and the DC values were calculated in the SENT_CORE network. hROIs with a * are those with significant correlations between activation and DC values (p < 0.05)

<figure><img src="image_11.png" /><figcaption>**Fig. 5** Correlation between DC values and activations in SENT_CORE across participants and across the 18 hROIs during each of the 3 language tasks. Plots of DC values and BOLD variations of the sentence minus word contrasts calculated for sentence production (left),</figcaption></figure>

sentence listening (middle), sentence reading (right) and degree centrality. The positive correlation coefficients (N=138×18=2484) are 0.158 for sentence production, 0.216 for sentence listening, and 0.295 for sentence reading

networks had lower inter-individual consistency, one including visual language areas at the interface between visual and syntactic processing (SENT_VISU) and the other including posterior DMN areas including posterior hippocampus (SENT_MEM). Intrinsic connectivity analysis showed that SENT_CORE was negatively correlated with SENT_MEM but was not correlated with SENT_VISU. Graph analysis metrics obtained for the SENT_CORE network revealed that F3t, STS3, and STS4 were hubs of both degree and

betweenness centrality, and F3t and ST3 were also hubs of participation, meaning that these are key areas for both intra-network communication and inter-network communication between SENT_CORE and the other 2 networks. Importantly, a positive correlation across individuals was observed between the DC measured at rest and the strength of activation in most SENT_CORE regions, meaning that participants with higher DC values in a given region had higher activations than participants with lower DC values.
Moreover, such a positive correlation between the DC and activation was still significant when all regions of all participants in the three tasks were pooled, meaning that this was true regardless of the cortical area considered.

# Discussion

## Methodological issues

In this study, we selected right-handers from the BIL&GIN database because we previously demonstrated that these participants have a left hemisphere dominance for language at both the group level (Tzourio-Mazoyer et al. 2016) and the individual level (Zago et al. 2017), with only 5 (3%) of the 144 participants having a co-dominant right hemisphere. This sample group is optimal for selecting areas specific for sentence processing based on a conjunction of activations and leftward asymmetries. In addition, the inclusion of a fairly considerable number of participants (N = 144) provided us a high sensitivity for detecting supramodal sentence areas while minimizing the risk of overlooking some.

However, we must underline that the present atlas is not all-inclusive. First, we selected map regions involved in only high-order language processing and lexico-syntactic processing. Using the list of familiar words as the reference condition, we removed the dorsal route of language, including the phonological loop, responsible for articulation and sound-to-articulation mapping (Saur et al. 2008; Rauschecker and Scott 2009). In addition, the regions selected herein focused on the left hemisphere and did not account for right hemisphere-specialized aspects of sentence processing, such as emotional prosody (Beaucousin et al. 2007; Hurschler et al. 2012) and context processing (Grindrod and Baum 2003; Ferstl et al. 2005).

Second, the presence of susceptibility artefacts combined with averaging the large number of participants led to incomplete mapping of the inferior part of the temporal lobe, prohibiting us from documenting some areas, such as the basal language area in the anterior part of the fusiform gyrus. This essential language area, first identified using deep electrical recordings (Nobre et al. 1994), has been shown with positron emission tomography (PET) to be activated during both the production and auditory comprehension of language (Papathanassiou et al. 2000).

Third, small size regions may also be lacking in this atlas since we provided data at the hROI scale rather than at the voxel scale.

Concerning the clustering methods of correlation values at rest, a perfect match was observed between the Ward and average clustering methods (see supplementary material Table 1), and a good score was obtained with the Gaussian mixture model for global clustering at the 32 hROI levels.

The weakest score was that of the *k*-means method, and such a difference in clustering observed with *k*-means compared to that in the other 3 methods is consistent with the fact that, as reported by Thirion et al., *k*-means forms clusters spatially close and connected but with poor reproducibility using the sample studied. By contrast, hierarchical clustering using Ward's method, which we selected to segregate the networks, was reported to create connected clusters that are highly reproducible using the studied samples (Thirion et al. 2014).

## A large set of supramodal language areas is involved in sentence processing tasks

We carefully designed each of the language tasks such that joint analyses were possible; the design was identical in the three tasks, and we chose to make them close enough to allow comparisons and conjunctions in terms of the number of words or the complexity of sentences. As mentioned above, the use of a high-level verbal reference task for controlling the involvement of primary areas (auditory, visual and motor) and removing phonological and automatic word processing kept the lexico-syntactic aspects common to all three tasks.

The first set of 32 hROIs provides left hemispheric regions that are dedicated to the monitoring and completion of tasks based on sentence processing. Although not all regions can be considered essential language areas, all were determined to be modulated by the verbal material with which they are associated (left activation and leftward asymmetry) and are thus part of an extended language network functioning during language tasks.

Unsupervised hierarchical clustering based on the resting-state connectivity between hROIs successfully segregated three different networks, including networks hosting core language, visual areas, and posterior areas of the DMN and posterior hippocampus. Within the systems to which they belonged, these networks hosted areas dedicated to the interaction/interface with language systems. For example, the current analysis extracted from among the visual areas involved in picture processing those areas specifically dealing with picture-sentence meaning integration.

## Sentence comprehension essential network (SENT_CORE)

Clustering the resting-state correlation between these 32 hROIs allowed the discrimination of SENT_CORE, a network of 18 strongly and positively correlated hROIs, including frontal and temporo-parietal hROIs located on the lateral surface of the left hemisphere and anterior insula areas. In particular, SENT_CORE included areas of the antero-posterior language networks, named in reference to the
Broca–Wernicke model in aphasia literature and reported with consistency in meta-analyses of healthy individuals mapped during language tasks see Fig. 6 (Vigneau et al. 2006; Price 2010, 2012). Note that SENT_CORE was the largest network in terms of volume (in mm³), as it included more than half of the hROIs (18/32), all of which were strongly activated and leftward asymmetrical.

In the following, we discuss the potential roles of the identified areas in relation to the literature. However, it is now acknowledged that, apart from very specific regions where a lesion can be closely associated with a specific defect, the role of a given area documented with functional imaging must be understood as the combination of its functional properties with those of the regions with which it constitutes a network to complete a given cognitive task. For example, prec4 is not part of the regions commonly labelled as “language areas”. The present work shows that prec4, located at the junction between the precentral and middle frontal gyrus, is both strongly activated and leftward asymmetrical in the three sentence tasks. Language meta-analyses have reported prec4 as part of the language

areas involved in lexico-syntactic processing (named F2p in Vigneau et al. 2006, Fig. 6), and in word selection and hierarchical sequencing (named dPrec in Price 2010, Fig. 6). Applying Neurosynth to prec4 coordinates ($x = -42.2, y=0.7, z=50$, Table 3) with an association test reveals the greatest number of studies with the terms “sentence”, “comprehension”, “language”, and “sentences”, followed by “eye” and “premotor”. Jouen et al. (2018) propose that prec4 is “involved in the understanding of actions during verbal and non-verbal tasks”. In the present protocol, the sentences involved human actions, closely consistent with that role proposed by Jouen et al. (2018). Using a network approach, Saur et al. (2008) found that prec4 belongs to the sentence comprehension functional network. In Glasser’s atlas, prec4 corresponds to “language area 55b” (Glasser et al. 2016) and overlaps the posterior part of the “rostro-ventral module” described by Genon et al. (2018a, b). In this last work based on peaks meta-analysis, this module is connected with the inferior frontal gyrus, orbital frontal and inferior parietal as prec4 in the present work. The present study further demonstrates that prec4 is involved in a supramodal manner during

<figure><img src="image_12.png" /><figcaption>**Fig. 6** Schematic comparison of SENSAAS with the literature. This figure shows the results of neuroimaging meta-analyses and clinical studies superimposed on the hROI of SENSAAS of the lateral surface of the left hemisphere of the BIL&GIN display template. In the first row: left SENSAAS hROIs of the left hemisphere lateral surface; middle: clusters of the meta-analysis of semantics (red) and sentence processing (green) adapted from Vigneau et al. (2006) with their labels: right: schematic representation of the meta-analysis of language-related activation studies (adapted from Price (2012); sentence:</figcaption></figure>

purple; semantics: light and dark pink; visual: red; word retrieval: green; integration: white). In the second row, left: functional connectivity of the left middle temporal gyrus centred on the site where lesion results in deep aphasia (orange, adapted from Turken 2011); middle: functional connectivity from a seed centred on the left inferior frontal gyrus (BA 45, red, adapted from Margulies and Petrides 2013); right: zones of atrophy observed when pooling all types of PPA (orange, adapted from Mesulam et al. 2014)
sentence processing, is strongly leftward asymmetrical, and
is strongly and positively connected at rest with the network
of areas we named SENT_CORE that hosts essential lan-
guage areas. These findings are consistent with the fact that
prec4 is part of the areas that are conjointly atrophied in
patients suffering from nonfluent primary progressive apha-
sia (PPA) (Mesulam et al. 2014) (Fig. 6).

Indeed, in the frontal lobe, SENT_CORE includes the
whole inferior frontal gyrus (F3t and F3O) corresponding
to Broca's area according to most authors, whose lesion is
responsible for conversational deficits, as is also the case for
the anterior insula (INSa2, INSa3) (Borovsky et al. 2007).
The posterior part of the inferior frontal sulcus (f2_2),
also part of SENT_CORE, has been underlined as an area
involved in lexico-syntactic processing (named F3opd in
Vigneau et al. 2006, Fig. 6). Price also targeted f2_2 as
being involved in word selection and hierarchical sequenc-
ing (named mFG in Price 2010, Fig. 6). Another indication
of the important role of f2_2 in language is that together
with F3t and prec4, it is the location of atrophy in nonfluent
PPA [(Mesulam et al. 2014) see Fig. 6].

In the medial part of the frontal lobe, SENT_CORE
includes both preSMA and the superior frontal gyrus (here,
SMA2, SMA3 and F1_2), which have been reported in tasks
involving sentences dealing with characters (Hervé et al.
2012), as in the present paradigm. The activation of these
areas has been attributed to processing the social aspects of
verbal material (Ferstl and von Cramon 2002). As discussed
below, these medial frontal areas are strongly connected with
the posterior part of the middle temporal gyrus (Turken and
Dronkers 2011). Notably, these medial frontal areas are also
the sites of atrophy in PPA (Wilson et al. 2010a, b; Tetzloff
et al. 2018).

In the temporal lobe, SENT_CORE included the STS
together with the posterior part of the middle and superior
temporal gyri, extending to the angular and supramarginal
gyri. This characterization is consistent with the proposal
by Price that these areas are involved in amodal semantic
combinations, a process common to the three sentence
tasks (Price 2010). A recent work on the time course of
sentence processing areas has also shown that the posterior
temporal lobe, here corresponding to STS3 and T2_3, is
involved in lexico-syntactic processing, with the process-
ing of individual words in relation to the syntactic structure
(Matchin et al. 2018). These posterior temporal areas have
been documented to be essential language areas since the
lesion of each results in specific deficits in sentence com-
prehension (Dronkers et al. 2004). In particular, lesion of
the posterior parts of the middle temporal gyrus results in
deep aphasia due to a loss of word meaning (Dronkers et al.
2004), consistent with the specific atrophy observed in logo-
penic PPA (Wilson et al. 2010a, b). STS4, located dorsally
to STS3 and T2_3, is involved in sentence-level high-order

processing, and Matchin, in his work on the time course
of sentence processing investigated with MEG, proposed
that “increased activation at the end of the sentence sug-
gests a response associated with the interpretation of the
sentence” (Matchin et al. 2018). The integration gradient
of sentence meaning from lexical to conceptual progresses
posteriorly with SMG7, which corresponds to the area
labelled STSp in Vigneau et al.’s meta-analysis (Fig. 6).
Vigneau et al. reported the following: “STSp role seems
to process the semantic integration of complex linguistic
material. This statement comes from the observation that
it is recruited when subjects listen to coherent rather than
syntactically or pragmatically incoherent sentences (Kuper-
berg et al. 2000; Luke et al. 2002), and it is involved in con-
text processing and syntactic generation—more activated
when subjects have to choose between two words to end
a sentence or have to generate the final word of a sentence
(Kircher et al. 2001). STSp activity is very likely related to
the linkage of linguistic structure to meaning: it is more acti-
vated when sentences are linked as dialogue (Homae et al.
2002) or syllogisms (Goel et al. 1998) than when they are
unlinked. It is also more activated during text comprehen-
sion, either presented auditory [compared to reverse speech
(Kansaku et al. 2000; Crinion et al. 2003) or words (Jobard
et al. 2007)] or visually [compared to words (Jobard et al.
2007) or pseudo-word reading (Vingerhoets et al. 2003)]”..
Price made an alternative proposal in her meta-analysis in
2000 stating that the equivalent of SMG7 (vSMG and pPT)
is activated by sentence processing, particularly when the
sentences increase in difficulty, and she therefore suggested
that SMG7 is involved in subvocal articulation. The present
results demonstrating that SMG7 is involved and leftward
asymmetrical not only during production and reading but
also during sentence listening, which was cited as a easy
task by the participants, suggests its involvement in meaning
elaboration. The angular gyrus is the final location where
this integration towards conceptual knowledge operates dors-
sally. Known to be involved in lexico-semantic processing,
as identified by language task-induced activation studies
[(Vigneau et al. 2006; Price 2010), see Fig. 6], the angular
gyrus is specifically involved in conceptual knowledge as
shown by neuropsychological studies. Lesion of this area
leads to the inability to associate a sound or image related
to the same concept (Saygin et al. 2004), and inactivation of
the angular gyrus in the left hemisphere induces a deficit in
semantic integration (Price et al. 2016). The angular hROI
identified herein corresponds to Wang’s C5 and C6 parcels,
involved in both language and theory-of-mind tasks (Wang
et al. 2017). This last observation confirms that the role of
the angular gyrus includes the interfacing between sentence
processing and the understanding of human actions, a pro-
cess that is part of the present sentence tasks.

<figure><img src="image_13.png" /></figure>
The similitude between the areas composing the SENT_
CORE network and the regions showing atrophy in all types
of PPA [(Mesulam et al. 2014), see Fig. 6] is striking and a
key element evidencing that SENT_CORE contains essen-
tial language areas, although the anterior part of the inferior
temporal gyrus hosting the basal temporal language area
(Nobre et al. 1994; Papathanassiou et al. 2000) is lacking
for the methodological reasons discussed above. Finally,
the language functional network, which shares a common
molecular basis of interaction highlighted by Zilles et al.
using neurotransmitter receptor fingerprints (Zilles et al.
2015), overlaps the hROIs of SENT_CORE, apart from
AG2. Based on Zilles et al. report, we can hypothesize that
the two other sentence-related networks we have identified
are likely to have different fingerprints from those of SENT_
CORE. However, because the hROIs of these networks were
selected as both activated and leftward asymmetrical in the
three sentence tasks, it would be of interest to investigate
what subset of fingerprints the hROIs of SENT_VISU and
SENT_MEM may share in common with SENT_CORE to
deepen our understanding of network interactions.

## Hubs of SENT_CORE

The computation of betweenness and degree centralities allowed us to identify three hubs within SENT_CORE: one in the inferior frontal gyrus (F3t) and two others along the posterior part of the temporal cortex (STS3 and STS4). These high centrality figures demonstrate that these hROIs play a central role in communications with other parts of SENT_CORE, and the participation indices of F3t and STS3 further indicate their key role in communication with the other networks identified. Such properties are consistent with the definition of epicentres proposed by Mesulam et al. (2014), and these regions can be considered essential for sentence/test comprehension independent of the modality. According to Mesulam, the network epicentre specializes in a specific behavioural component, which is language in this study, and the destruction of transmodal epicentres causes global impairments. In fact, as noted above, these hROIs overlap the middle temporal gyrus targeted in the aphasics investigation by Dronkers et al. (2004) as well as the regions showing the highest atrophy in all types of PPA (Mesulam et al. 2014), confirming that they do correspond to epicentres. In addition, these hubs are distributed in the anterior and posterior cortices, constituting an antero-posterior loop across areas belonging to the same hierarchical level in terms of cortical organization, consistent with Fuster's model for cognition (Fuster and Bressler 2012).

Other studies involving functional connectivity at rest
have shown the strong connection of the left inferior frontal
gyrus at rest with all SENT_CORE areas [(Margulies and
Petrides 2013), Fig. 6]. The similarity between the regions

found by these authors when seeding the posterior part of
the middle temporal gyrus (here, STS3 and STS4) and the
areas constituting SENT_CORE is also clear [(Turken and
Dronkers 2011), Fig. 6]. These two observations are consist-
ent with the present observation that these regions are hubs
for sentence processing and, more generally, the language
comprehension network. Such a framework is also thought
to correspond to the Broca–Wernicke language model, with
F3t serving as Broca’s area and STS3 and STS4 serving as
regions involved in the supramodal integration of meaning,
consistent with the location of posterior areas leading to
comprehension deficits (Pillay et al. 2017). These posterior
STS areas are also considered by Binder as areas supporting
meaning integration during sentence comprehension (Binder
2017). Considering that a left deficit in this area leads to
deficits in language comprehension, we propose to label it
Wernicke’s area, although as reviewed by Binder, such a
definition is different from the location currently proposed
for Wernicke’s area in the superior temporal gyrus. In the
present work, closely adhering to an anatomo-functional
definition in reference to deficits in comprehension associ-
ated with Wernicke’s aphasia (Dronkers et al. 2004; Binder
2015, 2017), we propose that F3t, STS3 and STS4, regions
strongly activated and asymmetrical during sentence pro-
cessing in different modalities and hubs of the SENT_CORE
network, are the epicentres of sentence comprehension.

## Other areas contributing to sentence tasks

Because there was a low inter-individual consistency in
the segregation of the networks labelled SENT_VISU and
SENT_MEM, we will discuss the potential role of the 14
hROIs that were not part of SENT_CORE together.

Since each event of the sentence tasks included task mon-
itoring, such as shifting between the word list and sentence
tasks when a picture was presented or providing a motor
response at the end of sentence processing, the involvement
of executive areas was expected. Conjunction analyses
indeed revealed that, in addition to the anterior insula at play
in the three tasks (see above), putamen hROIs were involved
and are a key neural support of executive functions and task
monitoring (Monchi et al. 2006; Sefcsik et al. 2009). Moreo-
ver, a recent meta-analysis involving connectivity modelling
showed that the left and right putamen areas are different in
terms of their respective co-activations, which are specifi-
cally co-involved in language areas (Viñas-Guasch and Wu
2017).

The posterior cingulate, precuneus, and paracentral lobule together with the posterior hippocampus are part of the DMN, which has been shown to be involved in both episodic thinking and processing of “self”. At rest, these areas constituted a network that had a very strong and negative correlation with SENT_CORE, confirming that, although
they belong to different networks, they are related at rest. Task-induced activation studies have demonstrated that posterior DMN regions are part of mind-reading areas together with the angular gyrus (Spreng et al. 2009). Activated and leftward lateralized during language tasks, this set of areas likely interacts with SENT_CORE during language tasks to process the sentence content dealing with social interactions (Hervé et al. 2012). The angular gyrus is a likely candidate for such an interaction since it segregated with SENT_MEM areas in two clustering methods and is involved in both language and theory-of-mind tasks (Wang et al. 2017).

The fact that these posterior areas were segregated with the posterior hippocampus, involved in both the actual perception and the encoding of scenes (Zeidman and Maguire 2016), suggests that the SENT_MEM network is also involved in the image processing of the drawings, a component common to the three sentence tasks. In fact, all participants performed well in the recall of these sentences (more than nine out of ten images recalled).

Finally, three hROIs were located along the occipito-temporal junction. The posterior part of the inferior temporal lobe (T3_4), the inferior occipital gyrus (O3_1) and the mid-fusiform gyrus (FUS4), with activity and leftward asymmetry, are likely related to the processing of the sentence content in relation to the one-second drawing presentation (minus the scramble version of the word-list reference tasks) common to all tasks. In her meta-analysis, Price considered these areas to be involved in direct visual combinations (Price 2010), consistent with the design of the present paradigm wherein the participants dealt with images related to sentence content regardless of the sentence task modality and the fact that the mean value and asymmetry of this network did not vary with the language tasks. Furthermore, in her 2012 review, Price considered these areas to be involved in sentence processing depending on the task demand (Price 2012). FUS4 is a region of the ventral route that corresponds to the visual word form area (Mellet et al. 2018). It is notable that this region was activated and leftward asymmetrical during the three tasks and segregated with the other visual regions at rest. Such a behaviour is in accordance with Price and Devlin's proposal that the mid-fusiform area is involved not only in word and visual processing but is a multimodal area whose function is determined by the set of interacting areas (Price and Devlin 2003). FUS4 involvement in the present study is likely to be in the association between the meaning conveyed by the pictures and the sentences and appears to be at play whatever the sentences' modality. The fact that FUS4 did segregate with visual areas rather than sentence areas at rest supports this hypothesis (Fig. 2). In fact, visual hROIs were not correlated with SENT_CORE networks, meaning that they do not exchange information at rest. However, these areas involved in visual processing, as well as the other hROIs that did not constitute a network

at rest with sentence core areas, showed a strong leftward activation during the language tasks, suggesting that their involvement during sentence is, at least in part, related to top-down influences from language networks.

**The degree of centrality measured in SENT_CORE explains the activation variability during the three language tasks**

To our knowledge, this study is the first to report a positive correlation between DC and task-induced activation values. Considering that the correlations between DC values and activation were positive for all hROIs, we suggest that the mechanism underlying this result may be that regions more highly connected within the SENT_CORE network are more highly recruited during language tasks. Such a hypothesis is supported by a very recent work conducted with electrocortical stimulation that demonstrated that critical language sites have significantly higher connectivity (Rolston and Chang 2018). It is also important to note that for 8 of the 13 hROIs, no significant effect of the nature of the task on the correlation strength was observed: the positive correlation did not differ regardless of whether the participants produced, read or listened to sentences compared to listing words. Such a relationship was also observed when plotting all SENT_CORE regions of all subjects, as reported by others in a different context (Buckner et al. 2009), meaning that the relationship between resting-state intrinsic connectivity and activation strength remains regardless of the participant and cortical area being considered. In Buckner's study, the authors hypothesized that beta amyloid accumulates in high DC regions because their high metabolism makes them more vulnerable to the disease, consistent with the report that the DC calculated at the voxel level in the entire hemisphere is correlated with cerebral blood flow values (Liang et al. 2013). Similarly, DC values computed at the hROI level in the SENT_CORE network may also indicate the metabolism of language areas in a given individual, which may be of interest in the evaluation of pathological states. Mesulam et al. reported the selective atrophy of right hemisphere areas in a PPA patient with right hemisphere dominance for language, an observation suggesting that language networks are specifically targeted by this illness (Mesulam et al. 2005). DC measurements in SENT_CORE may thus be a valuable index with which to evaluate inter-individual variations in language area activities in relation to anatomical and clinical patterns in such pathologies.

Previous investigations dealing with the relationships between tasks and the resting state have compared the functional connectivity during cognitive tasks with that measured during the resting state (Cole et al. 2014; Gerchen and Kirsch 2017), as reviewed by Wig (2017). Reports more closely related to the approach used herein have compared resting-state

<figure><img src="image_14.png" /></figure>
connectivity and hemispheric activation asymmetries obtained during language production in healthy individuals (Joliot et al. 2016) and epileptic patients (Doucet et al. 2014), or during story listening (Raemaekers et al. 2018). Other investigations have found correlations between the task-induced and intrinsic connectivity asymmetries measured in selected sets of hROIs (Liu et al. 2009), such as the set of regions involved in a semantic decision task (Wang et al. 2014). Seeding specific language areas, a previous study involving epileptic patients demonstrated that the functional connectivity measured at rest was correlated with lateralization indices measured during language tasks. In this last work, the intrinsic connectivity most strongly correlated with asymmetry of task-induced fMRI was between the inferior frontal and temporo-parietal language areas, all regions constituting SENT_CORE (Teghipco et al. 2016). The set of consistent findings across the latter studies show that asymmetries of intrinsic connectivity partially explain the variability in activation asymmetries measured during language tasks. However, to our knowledge, there are no previous reports of a relationship between intrinsic DC and activation strength, probably because such a relationship is lacking when comparing whole-brain intrinsic connectivity to variations in activity triggered by cognitive tasks that are underpinned by specific networks. In fact, in the present study, this relationship was observed when the DC computation considered only the 18 hROIs of the SENT_CORE network and disappeared when the DC computation included all of the left hemisphere hROIs. This is an important observation because it underlines the necessity of exploring the properties of intrinsic connectivity within specific networks rather than at the whole-brain level.

Considering that the rs-fMRI acquisition was completed on average 11 months before the language fMRI session herein, it was surprising that the rs-fMRI-derived DC values explained up to 11% of the variance in the language fMRI-derived activation amplitudes in 12 of the 18 SENT_CORE network regions. Thus, DC values at rest in regions constituting SENT_CORE can be considered as proxies of their potential involvement during sentence processing. However, to generalize this observation, we need to both investigate how the SENT_CORE DC is modified in individuals typical for language and to confirm that such a relationship between the DC value during the resting state and activation strength also exists for networks supporting other cognitive domains, specifically the attentional system, in the same participant.

## Conclusion

Based on the fMRI analysis of 3 language tasks performed by 144 healthy adult right-handers combined with the analysis of intrinsic resting-state connectivity in 138 of the participants,

we propose a SENSAAS atlas of high-order sentence processing areas. This atlas includes 32 regions decomposed into 3 networks, including one (SENT_CORE) specifically composed of essential areas for sentence reading, listening and production. This atlas also contains the features of these three networks, a graph analysis of the intrinsic connectivity of regions that compose SENT_CORE (their degree centrality values that correlate with their strength of activation during the language tasks) as well as the relationships across the different language networks at rest. Such a positive correlation between the DC at rest and the language task-induced activation amplitude in the left hemisphere language network opens the way for investigating participants with language pathologies or population neuroimaging studies searching for the genetic basis of language by analysing only resting-state acquisition. Finally, the methodology we applied, identifying regions from activation studies for selecting the networks at play, advanced the specificity of resting-state graphical analysis and shed light on the relationships between resting-state and task-related networks.

**Acknowledgements** The authors thank Daniel Margulies for having provided images from his work to complete the Fig. 6 and enrich the discussion.

**Funding** This work was supported by a grant from the FLAG-ERA Human Brain Project 2015 (ANR-15-HBPR-0001-03-MULTI-LAT-ERAL) awarded to FC, BM, NTM, and MJ and the Ginesis-Lab project (ANR 16 LCV2 0006 01), a grant awarded to MJ, BM, NTM, EM, CF, and LP.

## Compliance with ethical standards

**Conflict of interest** The authors declare that they have no conflict of interest.

**Ethical approval** The study has been approved by the Basse-Normandie research ethics committee and was performed in accordance with the ethical standards as laid out in the 1964 Declaration of Helsinki and its later amendments or comparable ethical standards.

**Informed consent** Informed written consent was obtained from all individual participants included in the study.

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

## References

Amtums K, Zilles K (2012) Architecture and organizational principles of Broca's region. *Trends Cogn Sci* 16:418–426
Baldo JV, Wilkins DP, Ogar J, Willock S, Dronkers NF (2011) Role of the precentral gyrus of the insula in complex articulation. *Cortex* 47:800–807

Beaucousin V, Lacheret A, Turbelin MR, Morel M, Mazoyer B, Tzourio-Mazoyer N (2007) fMRI study of emotional speech comprehension. *Cereb Cortex* 17:339–352

Binder JR (2011) Functional MRI is a valid noninvasive alternative to Wada testing. *Epilepsy Behav* 20:214–222

Binder JR (2015) The Wernicke area: modern evidence and a reinterpretation. *Neurology* 85:2170–2175

Binder JR, O'M (2017) Current controversies on Wernicke's area and its role in language. *Curr Neurol Neurosci Rep* 17: 58

Borovsky A, Saygin AP, Bates E, Dronkers N (2007) Lesion correlates of conversational speech production deficits. *Neuropsychologia* 45:2525–2533

Buckner RL, Sepulcre J, Talukdar T, Krienen FM, Liu H, Hedden T, Andrews-Hanna JR, Sperling RA, Johnson KA (2009) Cortical hubs revealed by intrinsic functional connectivity: mapping, assessment of stability, and relation to Alzheimer's disease. *J Neurosci* 29:1860–1873

Charrad M, Ghazzali N, Boiteau V, Niknafs A (2014) NbClust: an R package for determining the relevant number of clusters in data set. *J Stat Soft* 61:1–36

Clos M, Amunts K, Laird AR, Fox PT, Eickhoff SB (2013) Tackling the multifunctional nature of Broca's region meta-analytically: co-activation-based parcellation of area 44. *Neuroimage* 83:174–188

Cole MW, Bassett DS, Power JD, Braver TS, Petersen SE (2014) Intrinsic and task-evoked network architectures of the human brain. *Neuron* 83:238–251

Crinion JT, Lambon-Ralph MA, Warburton EA, Howard D, Wise RJ (2003) Temporal lobe regions engaged during normal speech comprehension. *Brain* 126:1193–1201

Diedenhofen B, Musch J (2015) cocor: a comprehensive solution for the statistical comparison of correlations. *PLoS One* 10:e0121945

Dodoo-Schittko F, Rosengarth K, Doenitz C, Greenlee MW (2012) Assessing language dominance with functional MRI: the role of control tasks and statistical analysis. *Neuropsychologia* 50:2684–2691

Doucet G, Naveau M, Petit L, Delcroix N, Zago L, Crivello F, Jobard G, Tzourio-Mazoyer N, Mazoyer B, Mellet E, Joliot M (2011) Brain activity at rest: a multiscale hierarchical functional organization. *J Neurophysiol* 105:2753–2763

Doucet GE, Pustina D, Skidmore C, Sharan A, Sperling MR, Tracy JI (2014) Resting-state functional connectivity predicts the strength of hemispheric lateralization for language processing in temporal lobe epilepsy and normals. *Hum Brain Mapp* 36:288–303

Dronkers N, Ogar J (2004) Brain areas involved in speech production. *Brain* 127:1461–1462

Dronkers NF, Wilkins DP, Van Valin RD, Redfern BB, Jaeger JJ (2004) Lesion analysis of the brain areas involved in language comprehension. *Cognition* 92:145–177

Dym RJ, Burns J, Freeman K, Lipton ML (2011) Is functional MR imaging assessment of hemispheric language dominance as good as the Wada test?: a meta-analysis. *Radiology* 261:446–455

Efron B, Halloran E, Holmes S (1996) Bootstrap confidence levels for phylogenetic trees. *Proc Natl Acad Sci* 92:13429

Ferstl EC, von Cramon DY (2002) What does the frontomedian cortex contribute to language processing: coherence or theory of mind? *Neuroimage* 17:1599–1612

Ferstl EC, Rinck M, von Cramon DY (2005) Emotional and temporal aspects of situation model processing during text comprehension: an event-related fMRI study. *J Cogn Neurosci* 17:724–739

Ferstl EC, Neumann J, Bogler C, Yves von CD (2008) The extended language network: a meta-analysis of neuroimaging studies on text comprehension. *Hum Brain Mapp* 29:581–593

Friederici AD, Gierhan SM (2013) The language network. *Curr Opin Neurobiol* 23:250–254

Fuster JM, Bressler SL (2012) Cognit activation: a mechanism enabling temporal integration in working memory. *Trends Cogn Sci* 16:207–218

Genon S, Reid A, Li H, Fan L, Müller VI, Cieslik EC, Hoffstaedter F, Langner R, Grefkes C, Laird AR (2018a) The heterogeneity of the left dorsal premotor cortex evidenced by multimodal connectivity-based parcellation and functional characterization. *Neuroimage* 170:400–411

Genon S, Reid A, Langner R, Amunts K, Eickhoff SB (2018b) How to characterize the function of a brain region. *Trends Cogn Sci* 22:350–364

Gerchen MF, Kirsch P (2017) Combining task-related activation and connectivity analysis of fMRI data reveals complex modulation of brain networks. *Hum Brain Mapp* 38:5726–5739

Glasser MF, Coalson TS, Robinson EC, Hacker CD, Harwell J, Yacoub E, Ugurbil K, Andersson J, Beckmann CF, Jenkinson M, Smith SM, Van Essen DC (2016) A multi-modal parcellation of human cerebral cortex. *Nature* 536:171–178

Goel V, Gold B, Kapur S, Houle S (1998) Neuroanatomical correlates of human reasoning. *J Cogn Neurosci* 10:293–302

Grindrod CM, Baum SR (2003) Sensitivity to local sentence context information in lexical ambiguity resolution: evidence from left- and right-hemisphere-damaged individuals. *Brain Lang* 85:503–523

Guimera R, Amaral LAN (2005) Functional cartography of complex metabolic networks. *Nature* 433:895–900

Hervé P-Y, Razafimandimby A, Vigneau M, Mazoyer B, Tzourio-Mazoyer N (2012) Disentangling the brain networks supporting affective speech comprehension. *Neuroimage* 61:1255–1267

Homae F, Hashimoto R, Nakajima K, Miyashita Y, Sakai KL (2002) From perception to sentence comprehension: the convergence of auditory and visual information of language in the left inferior frontal cortex. *Neuroimage* 16:883–900

Hubert L, Arabie P (1985) Comparing partitions. *J Classif* 2:193–218

Hurschler MA, Liem F, Jäncke L, Meyer M (2012) Right and left perisylvian cortex and left inferior frontal cortex mediate sentence-level rhyme detection in spoken language as revealed by sparse fMRI. *Hum Brain Mapp* 34:30182–33192

Jobard G, Vigneau M, Mazoyer B, Tzourio-Mazoyer N (2007) Impact of modality and linguistic complexity during reading and listening tasks. *Neuroimage* 34:784–800

Joliot M, Jobard G, Naveau M, Delcroix N, Petit L, Zago L, Crivello F, Mellet E, Mazoyer B, Tzourio-Mazoyer N (2015) AICHA: an atlas of intrinsic connectivity of homotopic areas. *J Neurosci Methods* 254:46–59

Joliot M, Tzourio-Mazoyer N, Mazoyer B (2016) Intra-hemispheric intrinsic connectivity asymmetry and its relationships with handedness and language lateralization. *Neuropsychologia* 93:437–447

Jouen AL, Ellmore TM, Madden-Lombardi CJ, Pallier C, Dominey PF, Ventre-Dominey J (2018) Beyond the word and image: II—structural and functional connectivity of a common semantic system. *Neuroimage* 166:185–197

Kansaku K, Yamaura A, Kitazawa S (2000) Sex differences in lateralization revealed in the posterior language areas. *Cereb Cortex* 10:866–872

Kircher TT, Brammer M, Tous Andreu N, Williams SC, McGuire PK (2001) Engagement of right temporal cortex during processing of linguistic context. *Neuropsychologia* 39:798–809

Klingbeil J, Wawrzyniak M, Stockert A, Saur D (2017) Resting-state functional connectivity: an emerging method for the study of language networks in post-stroke aphasia. *Brain Cogn.* https://doi.org/10.1016/j.bande.2017.08.005
Kuperberg GR, McGuire PK, David AS (2000) Sensitivity to linguistic anomalies in spoken sentences: a case study approach to understanding thought disorder in schizophrenia. *Psychol Med* 30:345–357

Lance GN, Williams WT (1979) INVER: a program for the computation of distance-measures between attributes of mixed types. * Aust Comput J* 11:27–28

Lebreit R, Iovleeff S, Langronnet F, Biernacki C, Celeux G, Govaert G (2015) Rmixmod: the R package of the model-based unsupervised, supervised and semi-supervised classification mixmod library. *J Stat Softw* 67:1–29

Liang X, Zou Q, He Y, Yang Y (2013) Coupling of functional connectivity and regional cerebral blood flow reveals a physiological basis for network hubs of the human brain. *Proc Natl Acad Sci USA* 110:1929–1934

Liu H, Stufflebeam SM, Sepulcre J, Hedden T, Buckner RL (2009) Evidence from intrinsic activity that asymmetry of the human brain is controlled by multiple factors. *Proc Natl Acad Sci USA* 106:20499–20503

Luke KK, Liu HL, Wai YY, Wan YL, Tan LH (2002) Functional anatomy of syntactic and semantic processing in language comprehension. *Hum Brain Mapp* 16:133–145

Margulies DS, Petrides M (2013) Distinct parietal and temporal connectivity profiles of ventrolateral frontal areas involved in language production. *J Neurosci* 33:16846–16852

Matchin W, Brodbeck C, Hammerly C, Lau E (2018) The temporal dynamics of structure and content in sentence comprehension: Evidence from fMRI-constrained MEG. *Hum Brain Mapp.* https://doi.org/10.1002/hbm.24403

Mazoyer B, Zago L, Jobard G, Crivello F, Joliot M, Perchey G, Mellet E, Petit L, Tzourio-Mazoyer N (2014) Gaussian mixture modeling of hemispheric lateralization for language in a large sample of healthy individuals balanced for handedness. *PLoS One* e.e101165

Mazoyer B, Mellet E, Perchey G, Zago L, Crivello F, Jobard G, Delcroix N, Vigneau M, Leroux G, Petit L, Joliot M, Tzourio-Mazoyer N (2015) BIL&GIN: a neuroimaging, cognitive, behavioral, and genetic database for the study of human brain lateralization. *Neuroimage* 124:1225–1231

Mellet E, Salagnon M, Majkic AC, Joliot S, Jobard M, Tzourio Mazoyer GB, d'Errico N F (2018) Neuroimaging supports the representational nature of the earliest human engravings. *BioRxiv.* https://doi.org/10.1101/464784

Mesulam M, Weintraub S, Parrish T, Gitelman D (2005) Primary progressive aphasia: reversed asymmetry of atrophy and right hemisphere language dominance. *Neurology* 64:556–557

Mesulam M-M, Rogalski EJ, Wieneke C, Hurley RS, Geula C, Bigio EH, Thompson CK, Weintraub S (2014) Primary progressive aphasia and the evolving neurology of the language network. *Nat Rev Neurol* 10:554–569

Monchi O, Petrides M, Strafella AP, Worsley KJ, Doyon J (2006) Functional role of the basal ganglia in the planning and execution of actions. *Ann Neurol* 59:257–264

Nobre AC, Allison T, McCarthy G (1994) Word recognition in the human inferior temporal lobe. *Nature* 372:260–263

Ojemann G, Ojemann J, Lettich E, Berger M (1989) Cortical language localization in left, dominant hemisphere. An electrical stimulation mapping investigation in 117 patients. *J Neurosurg* 71:316–326

Opsahl T, Agneessens F, Networks JS-S (2010) Node centrality in weighted networks: generalizing degree and shortest paths. *Soc Netw* 32:245–251

Papathanassiou D, Etard O, Mellet E, Zago L, Mazoyer B, Tzourio-Mazoyer N (2000) A common language network for comprehension and production: a contribution to the definition of language epicenters with PET. *Neuroimage* 11:347–357

Pillay SB, Binder JR, Humphries C, Gross WL, Book DS (2017) Lesion localization of speech comprehension deficits in chronic aphasia. *Neurology* 88:970–975

Price CJ (2000) The anatomy of language: contributions from functional neuroimaging. *J Anat* 197:335–359

Price CJ (2010) The anatomy of language: a review of 100 fMRI studies published in 2009. *Ann N Y Acad Sci* 1191:62–88

Price CJ (2012) A review and synthesis of the first 20 years of PET and fMRI studies of heard speech, spoken language and reading. *Neuroimage* 62:816–847

Price CJ, Devlin JT (2003) The myth of the visual word form area. *Neuroimage* 19:473–481

Price AR, Peelle JE, Bonner MF, Grossman M, Hamilton RH (2016) Causal evidence for a mechanism of semantic integration in the angular gyrus as revealed by high-definition transcranial direct current stimulation. *J Neurosci* 36:3829–3838

Raemaekers M, Schellekens W, Petridou N, Ramsey NF (2018) Knowing left from right: asymmetric functional connectivity during resting state. *Brain Struct Funct* 223:1909–1922

Raichle ME (2015) The brain's default mode network. *Annu Rev Neurosci* 38:433–447

Raichle ME, Mintun MA (2006) Brain work and brain imaging. *Annu Rev Neurosci* 29:449–476

Rauschecker JP, Scott SK (2009) Maps and streams in the auditory cortex: nonhuman primates illuminate human speech processing. *Nat Neurosci* 12:718–724

Richardson JD, Fillmore P, Rorden C, Lapointe LL, Fridriksson J (2012) Re-establishing Broca's initial findings. *Brain Lang* 123:125–130

Rolston JD, Chang EF (2018) Critical language areas show increased functional connectivity in human cortex. *Cereb Cortex* 28:4161–4168

Rubinov M, Sporns O (2010) Complex network measures of brain connectivity: uses and interpretations. *Neuroimage* 52:1059–1069

Saur D, Lange R, Baumgaertner A, Schraknepper V, Willmes K, Rijntjes M, Weiller C (2006) Dynamics of language reorganization after stroke. *Brain* 129:1371–1384

Saur D, Kreher BW, Schnell S, Kümmerer D, Kellmeyer P, Vry MS, Umarova R, Musso M, Glauche V, Abel S (2008) Ventral and dorsal pathways for language. *Proc Natl Acad Sci* 105:18035

Saygin AP, Wilson SM, Dronkers NF, Bates E (2004) Action comprehension in aphasia: linguistic and non-linguistic deficits and their lesion correlates. *Neuropsychologia* 42:1788–1804

Sefcsik T, Nemeth D, Janacsek K, Hoffmann I, Scialabba J, Klivenyi P, Ambrus GG, Haden G, Vecsei L (2009) The role of the putamen in cognitive functions—a case study. *Learn Percept* 1:215–227

Sidtis JJ (2007) Some problems for representations of brain organization based on activation in functional imaging. *Brain Lang* 102:130–140

Sneath PH, Sokal RR (1973) Numerical taxonomy. The principles and practice of numerical classification (a series of books in biology), vol 573. WF Freeman and Co., San Francisco

Sporns O, Honey CJ, Kötter R (2007) Identification and classification of hubs in brain networks. *PLoS One* 2:e1049

Spreng RN, Mar RA, Kim AS (2009) The common neural basis of autobiographical memory, prospection, navigation, theory of mind, and the default mode: a quantitative meta-analysis. *J Cogn Neurosci* 21:489–510

Suzuki R, Shimodaira H (2006) Pvclust: an R package for assessing the uncertainty in hierarchical clustering. *Bioinformatics* 22:1540–1542

Tate MC, Herbet G, Moritz-Gasser S, Tate JE, Duffau H (2014) Probabilistic map of critical functional regions of the human cerebral cortex: Broca's area revisited. *Brain* 137:2773–2782
Teghipco A, Hussain A, Tivarus ME (2016) Disrupted functional connectivity affects resting state based language lateralization. Neuroimage Clin 12:910–927

Tetzloff KA, Duffy JR, Clark HM, Strand EA, Machulda MM, Schwarz CG, Senjem ML, Reid RI, Spychalla AJ, Tosakulwong N, Lowe VJ, Jack CR, Josephs KA, Whitwell JL (2018) Longitudinal structural and molecular neuroimaging in agrammatic primary progressive aphasia. Brain 141:302–317

Thirion B, Varoquaux G, Dohmatob E, Poline JB (2014) Which fMRI clustering gives good brain parcellations. Front Neurosci 8:167

Thompson PM, Andreassen OA, Arias-Vasquez A, Bearden CE, Boedhoe PS, Brouwer RM, Buckner RL, Buitelaar JK, Bulayeva KB, Cannon DM, Cohen RA, Conrod PJ, Dale AM, Deary IJ, Dennis EL, de Reus MA, Desrivieres S, Dima D, Donohoe G, Fisher SE, Fouche JP, Francks C, Frangou S, Franke B, Ganjgahi H, Garaman H, Glahn DC, Grabe HJ, Guadalupe T, Gutman BA, Hashimoto R, Hibar DP, Holland D, Hoogman M, Pol HEH, Hosten N, Jahanshad N, Kelly S, Kochunov P, Kremen WS, Lee PH, Mackey S, Martin NG, Mazoyer B, McDonald C, Medland SE, Morey RA, Nichols TE, Paus T, Pausova Z, Schmaal L, Schumann G, Shen L, Sisodiya SM, Smit DJA, Smoller JW, Stein DJ, Stein JL, Toro R, Turner JA, van den Heuvel MP, van den Heuvel OL, van Erp TGM, van Rooij D, Veltman DJ, Walter H, Wang Y, Wardlaw JM, Whelan CD, Wright MJ, Ye J, ENIGMA C (2017) ENIGMA and the individual: predicting factors that affect the brain in 35 countries worldwide. Neuroimage 145:389–408

Tomasi D, Volkow ND (2012) Resting functional connectivity of language networks: characterization and reproducibility. Mol Psychiatry 17:841–854

Turken AU, Dronkers NF (2011) The neural architecture of the language comprehension network: converging evidence from lesion and connectivity analyses. Front Syst Neurosci 5:1

Tzourio-Mazoyer N, Joliot M, Marie D, Mazoyer B (2016) Variation in homotopic areas’ activity and inter-hemispheric intrinsic connectivity with type of language lateralization: an fMRI study of covert sentence generation in 297 healthy volunteers. Brain Struct Funct 221:2735–2753

van den Heuvel MP, Sporns O (2011) Rich-club organization of the human connectome. J Neurosci 31:15775–15786

van den Heuvel MP, Mandl RC, Stam CJ, Kahn RS, Hulshoff Pol HE (2010) Aberrant frontal and temporal complex network structure in schizophrenia: a graph theoretical analysis. J Neurosci 30:15915–15926

Vigneau M, Beaucousin V, Hervé PY, Duffau H, Crivello F, Houdé O, Mazoyer B, Tzourio-Mazoyer N (2006) Meta-analyzing left

hemisphere language areas: phonology, semantics, and sentence processing. Neuroimage 30:1414–1432

Viñas-Guasch N, Wu YJ (2017) The role of the putamen in language: a meta-analytic connectivity modeling study. Brain Struct Funct 22:3991–4004

Vingerhoets G, Van Borsel J, Tesink C, van den Noort M, Deblaere K, Seurinck R, Vandemaele P, Achten E (2003) Multilingualism: an fMRI study. Neuroimage 20:2181–2196

Wada J, Rasmussen T (1960) Intracarotid injection of sodium amytal for the lateralization of cerebral speech dominance: experimental and clinical observations. J Neurosurg 17:266–282

Wang D, Buckner RL, Liu H (2014) Functional specialization in the human brain estimated by intrinsic hemispheric interaction. J Neurosci 34:12341–12352

Wang J, Xie S, Guo X, Becker B, Fox PT, Eickhoff SB, Jiang T (2017) Correspondent functional topography of the human left inferior parietal lobule at rest and under task revealed using resting-state fMRI and coactivation based parcellation. Hum Brain Mapp

Ward JH Jr (1963) Hierarchical grouping to optimize an objective function. J Am Stat Assoc 58:236–244

Wig GS (2017) Segregated systems of human brain networks. Trends Cogn Sci 21:981–996

Wilson SM, Dronkers NF, Ogar JM, Jang J, Growdon ME, Agosta F, Henry ML, Miller BL, Gorno-Tempini ML (2010a) Neural correlates of syntactic processing in the nonfluent variant of primary progressive aphasia. J Neurosci 30:16845–16854

Wilson SM, Henry ML, Besbris M, Ogar JM, Dronkers NF, Jarrold W, Miller BL, Gorno-Tempini ML (2010b) Connected speech production in three variants of primary progressive aphasia. Brain 133:2069–2088

Yourganov G, Smith KG, Fridriksson J, Rorden C (2015) Predicting aphasia type from brain damage measured with structural MRI. Cortex 73:203–215

Zago L, Hervé P-Y, Genuer R, Laurent A, Mazoyer B, Tzourio-Mazoyer N, Joliot M (2017) Predicting hemispheric dominance for language production in healthy individuals using support vector machine. Hum Brain Mapp 38:5871–5889

Zeidman P, Maguire EA (2016) Anterior hippocampus: the anatomy of perception, imagination and episodic memory. Nat Rev Neurosci 17:173–182

Zilles K, Bacha-Trams M, Palomero-Gallagher N, Amunts K, Friederici AD (2015) Common molecular basis of the sentence comprehension network revealed by neurotransmitter receptor fingerprints. Cortex 63:79–89