Fast reproducible identification and large-scale databasing of individual functional cognitive networks
BMC Neurosciencevolume 8, Article number: 91 (2007)
Although cognitive processes such as reading and calculation are associated with reproducible cerebral networks, inter-individual variability is considerable. Understanding the origins of this variability will require the elaboration of large multimodal databases compiling behavioral, anatomical, genetic and functional neuroimaging data over hundreds of subjects. With this goal in mind, we designed a simple and fast acquisition procedure based on a 5-minute functional magnetic resonance imaging (fMRI) sequence that can be run as easily and as systematically as an anatomical scan, and is therefore used in every subject undergoing fMRI in our laboratory. This protocol captures the cerebral bases of auditory and visual perception, motor actions, reading, language comprehension and mental calculation at an individual level.
81 subjects were successfully scanned. Before describing inter-individual variability, we demonstrated in the present study the reliability of individual functional data obtained with this short protocol. Considering the anatomical variability, we then needed to correctly describe individual functional networks in a voxel-free space. We applied then non-voxel based methods that automatically extract main features of individual patterns of activation: group analyses performed on these individual data not only converge to those reported with a more conventional voxel-based random effect analysis, but also keep information concerning variance in location and degrees of activation across subjects.
This collection of individual fMRI data will help to describe the cerebral inter-subject variability of the correlates of some language, calculation and sensorimotor tasks. In association with demographic, anatomical, behavioral and genetic data, this protocol will serve as the cornerstone to establish a hybrid database of hundreds of subjects suitable to study the range and causes of variation in the cerebral bases of numerous mental processes.
Inter-subjects variability is a missing facet of the current neuroimaging literature [1–3], and until recently has been viewed more as a nuisance for brain imaging studies than as a relevant dimension to investigate the mechanisms of human cognition. Indeed, most of the published studies described the cerebral bases of various cognitive processes from voxel-based group analyses performed on the data from 10–15 subjects. Group analysis of a small collection of brains assures that the description of these functional invariants may be extended to other healthy subjects. However we usually do not know if a cerebral network involved in a task is homogenous enough among the healthy population to be analyzed in only one group or if several groups have to be considered, nor how many subjects are required to correctly describe different sub-groups  (This question was also recently addressed in  on the basis of the present database). Consequently, it is plausible that in many cases, especially in those involving associative areas in complex tasks, we just capture the common denominator of each individual cognitive circuit and lose a large amount of information.
Describing more completely the parts of cerebral networks used but not shared by all of our subjects require considering variability of brain activation, which may have various origins: ▪ Intra-subject inter-sessions variability due to movement artifacts, physiological noise, etc... [6, 7] ▪ Spatial variability caused by the shape and location of cortical sulci  even for tasks requiring low-level processing. ▪ Biological factors such as sex [9, 10], genotype [11–13], or protein expression  ▪ Cognitive skills or difficulties, which may reflect heterogeneity of the healthy ('control') population of volunteers . ▪ Cognitive strategies spontaneously chosen by subjects to perform a task [16–18] or constrained by the protocol  ▪ Education and learning, that may locally modulate activation or structural anatomy [15, 20, 21]. Exploring inter-individual variability thus requires investigating various types of co-variation in a multi-dimensional space.
Toward a multidimensional database
Characterizing this functional variability, particularly when considering the genetic level, ideally requires acquiring functional imaging data from hundreds of subjects and organizing these data into a large-scale database, together with genetic, behavioral and biomorphological data. Databasing and analysis of structural magnetic resonance images has already resulted in probabilistic anatomical atlases [22, 23]. However, a similar large scale description of functional networks is still lacking.
Given that we are in the early stages of exploration of the causes of inter-individual variability, it would be desirable for such a functional imaging study to cover a broad variety of cortical territories and to describe cerebral correlates associated with various level of cognitive processing, from simple perceptual processing to higher-level cognitive functions that require explicit learning and education. For instance, recent advances in the genetics of dysphasia, dyslexia, and dyscalculia have provided several candidate genes whose impact on inter-individual variability in the normal population remains unknown. Considering some cognitive tasks that have been extensively described in the neuroimaging literature, we chose to include: a mental calculation task to investigate superior fronto-parietal networks  and a language comprehension task which focuses on the inferior frontal and superior temporal lobes . Using auditory and visual stimulation allows us to isolate cortical pathways associated with perceptual processing (superior temporal sulcus and occipito-temporal cortex [26, 27]) while the use of conjunction analysis across modalities also allows us to isolate correlates of amodal processing (associative cortices, for instance). Finally, evolutionary and developmental models suggest that some primitive mechanisms may be (partially) shared between hand/finger motor representation, speech language areas and correlates of mental arithmetic [28–30], noticeably in frontal and parietal lobes. If such assumptions are verified, crossing analysis among these tasks may then help to dissect the task-related networks into a more subtle functional parcellation, and enlighten developmental issues of human brain organization. In brief, these considerations suggest that it would be particularly valuable to obtain images of the cerebral substrates of speech comprehension, reading, and calculation in a large number of subjects, associated with genetic, anatomical and behavioral data, in a highly standardized manner, and at a low cost.
As shown in Figure 1, we planed to acquire four types of data: functional images and a high-resolution anatomical scan for a fine description of sulci, grey and white matter, as well as (not described here) behavioral and personal data, aimed to create a rough cognitive profile of the subject, and DNA sampling (cheek swab). Recall that this data collection occurs within the constraints of this project: to be added to other running protocol, with a minimal cost of people and time.
A fast brain mapping sequence
In the present research, our goal was to define a simple fMRI test, less than 5 minutes long, that could delineate, in a subject-specific manner, those cerebral circuits. A functional sequence was added to each functional imaging session performed in our lab (Figure 1), taking advantage of the continuous flow of volunteers recruited for various protocols. Because we wanted to capture the maximal amount of functional information in the minimum amount of time, we designed the sequence according the following challenging constraints:
▪ the sequence had to be short, so as to disrupt as little as possible the main protocol. We choose 5 minutes for performing 100 trials.
▪ we aimed to obtain for each subject a description of different levels of functional architecture, from sensori-motor areas (perception and action) to more associative areas involved in reading, language processing and calculation.
▪ we aimed to capture in 5 minutes most of the individual networks related to each task.
▪ individual networks described in 5 min had to be reproducible over sessions and time.
The feasibility of using short stimulation designs (ranging from 10 to 25 min long) to reveal individual functional maps has been previously assessed for language mapping , for visual areas  and recently for a set of functional networks covering sensorimotor processes, working memory, executive functions and emotional processes . Beyond that point, the main goal of the current study was to use the data obtained with this fMRI protocol with individual subjects as the cornerstone of a large-scale hybrid database. Because the individual functional information that can be captured in such a short sequence should be considered with caution, we focus here on the design efficiency and within-subject. We then describe preliminary data obtained from 81 subjects scanned in a 3T scanner and address new methodological issues including statistical methods for analysis and visualization of inter-individual functional variability. Subsequent publications will exploit the potential of this database to focus on characterizing inter-individual variability.
Individual information captured in 5 minutes
Figure 2 shows the evolution of the statistics of four functional contrasts of two subjects with scanning time (T value, number of activated voxels, spatial location), allowing us to determine the minimal number of blocks required to reach a satisfactory description of the individual maps. Interestingly, all of the selected peaks exceeded the p < 0.001 threshold at voxel level after only one block and averaging t-values showed that one block was enough to reach about 60% of the final statistical significance. Locations of these peaks were on average 7 mm away from the final coordinates obtained after completion of 6 blocks. Finally, we reported that mean intra-subject and intra-session activation significance varies with block and across contrasts, suggesting that weaker activation peaks may be missing in some blocks when a fixed statistical threshold is applied. Importantly, this variability does not seem to reflect overall factors such as arousal or attention, which would have been predicted to yield similar profiles of activation for each contrast across sessions.
Reliability of individual maps obtained from one single block is assessed in Figure 3. Functional networks obtained from the first 5-minutes block appeared very similar to those obtained after a 30-minute session, except that the statistical significance increased with the number of blocks performed, thus allowing for a more stringent threshold. A first threshold-dependant quantitative approach, detailed in Table 1, showed that about the three-fourths of maxima reported for the two subjects after the completion of 6 blocks were already present as local or main maxima in the one-block analysis, with an averaged spatial shift of 6 mm (2 voxels), varying between 2 and 10 millimeters depending of subject and task. Conversely, about 40% of maxima isolated in the one block analysis did not appear as a local/main maxima in the six blocks analysis at the selected threshold. More interestingly, a second threshold-free approach (ROC analysis, Figure 3) revealed a high probability of a consistent classification of voxels into active and inactive categories across the 5-minutes and 30-minutes maps of a subject. For right hand action, auditory stimulation, calculation and reading task, respectively, the discriminative power (Dp) values were 1, 0.99, 0.98 and 0.87 for the first subject and 0.99, 0.99, 0.95 and 0.94 for the second one. We also reported a gradient of subject-specificity over contrasts: for right hand action, auditory stimulation, calculation and reading task, respectively, the averaged inter-subject Dp values were 0.95, 0.92, 0.78 and 0.63. This result confirms the sufficiency of the 5-minutes map and suggests high intra-subject consistency but lower cross-subject reliability for reading and calculation tasks.
Interestingly we directly compared individual activation maps obtained from the present protocol with those obtained from a more classic block-design paradigm completed during the same session. In the Additional File 1 we show six examples of within-subject reproducibility across experimental designs. These results suggest that the main individual foci of activation collected in our database truly reflect areas crucial for a cognitive task independently of the design (block vs. task shifting), number of trials and task details.
Examination of individual representative contrast images from subjects who participated in two fMRI sessions often indicated a high reliability of the activations (Figure 4, right column). Inter-subject variability is also illustrated and appears to vary with task: for the motor contrast, subjects displayed a similar pattern of activation around the left central sulcus. For the reading contrast, the inter-subject variability seems much higher, due to various combinations of activation, but mainly restricted to the left hemisphere (due to language left-lateralization). For calculation, spatial variability across individuals extended both in lateralization and in the relative amounts of frontal versus parietal activity. For instance, subject 5 strongly activated bilateral fronto-parietal regions, but only a subpart of this pattern was isolated for subject 2 and 3.
To quantify these observations, an inter-scans distance based on the calculation of a similarity coefficient were computed between activation patterns of the 13 subjects and 2 sessions and plotted in a reduced 3-dimensional space for three contrasts (Figure 4). A high degree of within-subject reliability was observed for reading and calculation tasks, as the two scans obtained in a given subject (dots with same color) were mostly tightly grouped together in this summary space. Considering the small size of the motor activations, we suspected that the distance computation was strongly influenced by voxels of non-interest. Indeed, when we masked the distance computation by the right-left motor activation from the RFX group analysis, all pairs of intra-subject scans were grouped together and the variance explained increased up to 66%. Overall, this suggests a high degree of similarity between intra-subject scans both in location and time course of activations. This observation was statistically assessed by a non-parametric Wilcoxon signed rank test computed on the 13 differences between inter-scan intra-subject distance (scan 1 – scan 2 distance) and averaged inter-scan inter-subject distances (scan 1 – 24 other subject's scans). P values for the difference between intra and inter-subject distances were 0.0215, 0.0398 and 0.0012 for motor, reading and calculation task, respectively.
Group analysis of main cerebral networks
A voxel-based random effect analysis (RFX) performed on the 81 subjects' contrast images allowed us to compare the efficiency of our fast protocol with results already described in the neuroimaging literature (see Figure 5 and Table 2 for a detailed listing of areas – maxima were labeled using the Anatomical Automatic Labeling (AAL) program package  except for visual areas where a functional labeling was preferred). Contrasts based on checkerboards, auditory and visual stimuli activate cortical areas related to perception and stimulus encoding: occipital lobe, several sites of the occipito-temporal pathway, auditory temporal cortex. Left and right hand action contrasts demonstrated activation of contralateral sensorimotor cortex, SMA and ipsilateral cerebellum, while the conjunction of audio and visual motor – sentences contrasts showed a large bilateral set of areas included SMA, rolandic operculum, cerebellum, thalami, postcentral and precentral areas. The main correlates of reading were a set of middle and superior temporal areas with a strong trend toward left lateralization, left frontal sites and inferotemporal areas and a set of subcortical sites. The conjunction of auditory and visual sentences mainly restricted the previous network to bilateral middle and superior temporal sites that surround auditory cortex and left frontal areas. The conjunction of audio and visual calculation minus sentences showed increased activation in bilateral intraparietal areas, putamen and left precentral gyrus. Because it is supposed that the cerebral mechanisms of mental numerical manipulation may partially overlap with visuospatial and language processing, we isolated areas common to calculation and reading: the conjunction of auditory, visual calculation and reading versus checkerboards isolated a small left lateralized network comprising SMA, precentral and temporal area. The conjunction of auditory, visual calculation and reading versus rest revealed additional left frontal and superior occipital sites, right cerebellum and right calcarine activation.
Preliminary database exploitation
We compared different approaches for the functional group analysis to test methods that should both improve the statistical power of the group analysis and automatically collect individual functional information of a large fMRI database. A set of highly significant temporal, frontal and parietal activations related to motor, reading and calculation processing were isolated with the voxel-based RFX and plotted on a left sagital slice (column 6a). However, the overlapping of individual functional maps thresholded at p < 10-2 at the voxel level shows that the activated voxels are rarely common to more than a third of the subjects (Figure 6b). Even in left motor cortex, the highest degree of overlap was about 66%. Conversely, some areas, like the left inferior temporal gyrus in the calculation contrast, were present in overlap maps but did not pass the RFX threshold. We enhanced the level of cross-subjects replication when using non-voxel based group analyses. The first one aimed to extract automatically for each contrast a list of individual maxima (called functional landmarks, or BFL) reliable across subjects but not located at exactly identical coordinates that conveys subject-specific information of activation magnitude, statistical significance and spatial location. Detection of single-subject peaks then reached 90% in left motor cortex and 93% in left intraparietal area (column 6c), that are supposed to be crucial for motor and calculation tasks, respectively. Interestingly, the average coordinates of the BFLs were close to the corresponding peak locations found in the overlap and RFX maps. The variances of BFLs location were comparable across all contrasts: the standard deviation error for Talairach coordinates was about 3.5–4 mm for most of the functional landmarks located in both sensori-motor and associative areas. A second procedure was performed to identify similar brain areas (parcel) across subjects mainly on the basis of their functional profile across conditions (Figure 6d). This functional parcellation isolated networks with a very similar topology to those observed with RFX analysis, but with an increase of sensitivity and more extended regions of activation functionally dissociated in subsets of areas: in Figure 7, we report in detail the functional profiles for some of these parcels. It appeared that the parietal activations reported in the voxel-based RFX correspond to a functional gradient, with an anterior site (area P1) strongly involved in all motor trials, the horizontal part of the IPS (P2) equally involved in motor action and calculation, while motor activation was essentially absent in the posterior intraparietal sulcus (P3) which seems more specifically activated by calculation trials, and shows a significant deactivation for language processing. Similar mosaics can be described in frontal lobe: descending from upper part of the precentral gyrus to inferior frontal sites, we observe an area strongly activated for all tasks except checkerboard flashing (parcel F1), then an area which shows preferential activation for motor and calculation tasks (F2). Below this (F3), preference for auditory sentence increases to reach its maximum specificity in (F4) close to Broca's area, but also in the anterior part of the superior temporal gyrus (T5). The posterior superior temporal areas show a high sensitivity to auditory material compared to visual ones (T2) whereas this tendency is inverted in a very close middle temporal area (T1) and reaches its maximal preference for visual word perception in the fusiform gyrus (T3).
Designing a functional database of hundreds of healthy brains while minimizing time and resources costs was the original question addressed by the present work. Here, we present an approach based on the accumulation of homogeneous fMRI data (in combination with behavioral and genetic data) from a sample of a hundred of subjects in a very brief standardized cognitive protocol.
Efficiency of the fast multi-functional localizer
Protocol duration was the main constraint in the current project. Adding a 5-minute functional scan to each fMRI acquisition performed in the lab could be done without difficulty. Our results indicate that this was sufficient to describe in a subject-specific and reliable way individual topologies of brain activations, with an average spatial accuracy of 2 voxels.
Furthermore, when several subjects were tested twice at an interval of several months, the topology of activations was largely reproducible, even for cognitive tasks such as mental calculation which might be expected to be subject to fluctuations due to learning or attention.
It is worth mentioning that for some subjects, similar patterns were obtained for the two sessions only when lowering the threshold of one session. As shown for subjects who performed 6 blocks in one session, the significance level of a given contrast may be subject to fluctuation even during a single session, and considering all subjects' contrasts images at a fixed threshold may obscure similarities and increases the risk of type I or type II errors in maxima detection. This observation is supported by the high discriminative power of individual statistical maps obtained after a 5 min session and strongly argues for a structural description of individual activations [35, 36], less dependent on anatomical location and statistical threshold but more consistent over time and sessions.
We emphasize that reproducibility should be assessed relative to our study's goals: as shown in Figures 3 and 4, some activation peaks can be missed by a given 5-minute scan, presumably due to statistical noise, but possibly also due to changes in strategies, attention, or task-related adaptation. In addition, experimental and/or analysis procedures are not infallible as illustrated by a few subjects for whom no BFL was extracted for motor tasks. However, in the present context, which is to create a database intended to perform behavioral-brain correlations at the scale of hundreds of individuals, this level of reliability is likely to suffice.
Cerebral networks covered by the protocol and cognitive issues
Group analysis performed on a first database of 81 subjects first allowed us to report several functional networks which fit with previous fMRI results in the fields of speech perception, reading, motor execution and number manipulation. Visual stimuli activated the ventral occipito-temporal pathway from basic retinotopic organization to a more anterior site of the left fusiform gyrus involved in visual word processing . Another gradient of functional specialization was observed in superior temporal gyri, from bilateral primary auditory cortices, easily detected with auditory stimulation, up to temporal areas located around the superior temporal sulcus (with a trend toward left hemisphere dominance) and recruited during multimodal sentence comprehension. They closely mirror the speech-processing areas detailed by Price et al. , with an anterior temporal activation associated to the processing of sentence structure, and a posterior middle temporal activation close to sites previously associated with sentence-level semantics . We also observed left frontal areas and SMA, which have been already reported for syntactic and semantic processing of sentences [40, 41]. The calculation task involved the bilateral intraparietal and fronto-cingular network classically reported as active during simple number manipulation [42–46], together with the bilateral putamen. These subcortical areas are not usually considered as a part of the core numerical system, but have been tentatively linked to sequential processing of multi-stage calculations or to the retrieval of arithmetical facts .
In addition to this coarse functional description of brain areas, due to the simplicity of tasks, the co-existence in our protocol of motor, language and calculation tasks allows examination of their respective neural correlates and gives new insights into some debated cognitive issues. For instance, a restricted subset of three left-lateralized cortical areas (precentral gyri, SMA, and posterior middle temporal gyrus) was isolated by examining the conjunction of sentence comprehension and calculation in both visual and auditory modalities. Although studies of patients have suggested a relative semantic and syntactic independence between language and arithmetic in the adult brain, such functional overlap may represent a core system on which language and calculation are articulated regardless of modality. In particular the left posterior temporal gyrus that has been described as a multimodal integration region [47, 48] could be a candidate for convergence of visual, verbal and non-verbal magnitude codes. The left superior occipital gyrus and right visual areas, involved in visual stimuli processing (reading task and checkerboard perception), were also equally recruited during calculation trials performed from visual or auditory inputs. This observation is compatible with the hypothesis that processes sustaining mental calculation may involve a top-down activation of a symbolic digital code  and may also share some cortical territories with visuo-spatial areas . It was possible to perform a complementary and improved functional dissection of cortical maps using a parcellation technique over experimental conditions and designed to compensate for inter-subjects anatomical differences. For instance, the two bilateral intraparietal clusters isolated in the calculation RFX, while very close to those reported literature, appeared spatially restricted compared to the individual activation maps where more anterior or posterior additional sites can be seen along the intraparietal sulcus. After individual parcellation of functional areas, the parietal areas related to calculation appeared more extended along the intraparietal sulcus, with an antero-posterior functional gradient that corroborates the geometrical layout reported by Simon et al. , as well as the hypothesis of distributed overlapping parietal representations proposed by . A detailed examination of the anterior parietal area indicates that a co-location of motor- and calculation-related activations exists at the single subject level, perhaps illustrating an extension of the numerical system to a sensori-motor representation of hands . Interestingly, a recent study of inter-individual variability of the infero-parietal cytoarchitecture showed a reproducible topography of areas that however vary in size and extent . This biological evidence supports the framework of our individual functional parcellation algorithm which assumed, as a methodological constraint, that parcels are connected in a similar spatial organization across brains. However, further investigations have to be done to see if functional parcellation matches is some respects cytoarchitectonic boundaries. This observation is also true for the mosaic of inferior frontal areas, which were underestimated, both in terms of spatial extent and statistical significance by the voxel-based RFX analysis; this lobe appears to have a complex functional topography, with a widespread intermingling of areas involved in motor response, language comprehension and mental calculation. Interestingly, the left inferior temporal area was absent at the selected threshold (Z > 4) for the calculation task voxel-based RFX analysis. Because the functional profile of this parcel shows a consistent involvement in both auditory and visual calculation tasks and because the overlap of activations across different individuals is only about 27%, an elevated degree of anatomical variability may be suspected in this area. In conclusion, the combination of tasks and modalities allows drawing a detailed functional continuum for middle and inferior frontal lobes, middle and inferior temporal lobes and superior parietal lobe, from more modality-dependant areas up to more abstract cortical regions.
Dealing with anatomical and functional variability at the group level
Even if cerebral correlates of sensori-motor, reading, language comprehension and calculation processes are well identified and reproducible at the group level, subject level analyses revealed that individual activations were distributed in more complex patterns. Here, one of the most variable patterns across subjects, though associated with the highest intra-subject reliability, was the substrate of calculation. All subjects presented several sites of activation in the intraparietal sulci and frontal lobes but their extents, locations and combination varied. Because this individual patterns were highly reproducible over sessions, not only for large but also for small clusters, inter-subject variability should not be considered as noise but as a physiological signature of subject's brain activation during the numerical task. Whether such individual organizations depend on anatomical or functional constraints is a fundamental question. While we still presented our results in the MNI coordinate system and do not take into account the underlying anatomical structures, we believe that the automatic extraction of individual functional landmarks presented here was a first step to estimate for each functionally defined area both realistic spatial variability and frequency in the population and provides a promising context to address this issue. For instance, spatial variability of the left frontal functional landmarks described in our results fits well with the standard deviation associated to normalized left frontal anatomical landmarks location , which ranges from 2.5 to 5.7 mm. This underlies how structural changes may directly affect the location of functional sites and matches well with our estimate of their spatial variability. A recent paper reported a similar conclusion considering only visual cortical areas, where spatial variability of sulcal features was found to reflect those of the functional topography . As suggested by Juch et al. , automatically extracting and labeling activation and cerebral structures (sulci) jointly at the individual level will help to describe reliable cerebral organizations while separating anatomical and functional sources of inter-individual variability. Interestingly, a preliminary analysis of the spatial distribution of the functional landmarks in the left hemisphere did not show strong differences between tasks, suggesting that higher cognitive circuits are not associated with a greater spatial variability than sensori-motor ones. These results must be viewed with great caution, as this might be due to the algorithm of BFL extraction, which is based on a recursive method primarily derived from the RFX group analysis. However, the range of variability reported here is similar to that reported in another fast functional localizer study , at least for left-lateralized language areas (Broadmann area 44/45 and 22/42).
Concerning the functional features that should constrain individual circuits, it is particularly tempting to describe some frontal or parietal associative areas as functional 'nodes' because of their implication in various of tasks. We suspect that these sites are highly relevant to define individual functional reference frames for brain organization, similar to 'sulcal roots' in global brain gyrification  that appear early in the fetal brain and seem to constraint the adult global cortical folding. In a similar vein, detection of systematic spatial co-variation/co-lateralization of activation sites with these functionally defined 'nodes' (e.g. between calculation- and language-related areas) may be informative about the links that exist between different cognitive functions, and may help to specify the developmental constraints or evolutionary roots of the functional cerebral networks in the adult brain .
Perspectives for database exploitation
Further analysis and extension of the database to a greater number of subjects will be needed to disentangle the various sources of inter-subject variability listed in the introduction. We estimate that a minimum of 150–200 subjects will be necessary to begin to describe the variety of activation patterns in the population as well as to reach the minimal statistical power required to correctly isolate sub-groups that can be characterized by a combination of behavioral, anatomical, physiological and genetic features. We have no a priori assumptions about the rate of these possible 'variants' in the healthy population. For example, we may recall that 5% of right-handed subjects present an absence of usual left-hemispheric dominance in language correlates. Atypical organization of functional networks may also be suspected in subjects with self-reported developmental difficulties in arithmetic (~5% of our 'normal subjects' sample) or in adults who reported suffering of verbal difficulties during childhood (~10% of our sample). Identifying and characterizing some of these variants in various cognitive domains (verbal, visuo-motor, visuo-spatial, numerical...) may require hundreds of subjects. As demonstrated by a few studies, mostly from clinical populations, genetic variation may also contribute to the phenotype of cerebral activation in verbal  or numerical cognition [14, 56]. A recent work has emphasized that genetic studies of human cognition might greatly benefit from considering a continuum of cognitive abilities in healthy populations . In the further, we plan to explore the impact of genotype on the structural and physiological features of non pathological functional networks, considering the normal range of variation in subjects' verbal, visuo-spatial and numerical abilities.
We designed a fast multi-functional localizer paradigm that is now routinely applied in our laboratory to isolate, in five minutes, individual cerebral correlates of visual, auditory and sensorimotor processes, reading, language comprehension and mental calculation. We showed here that isolated activation patterns vary noticeably between subjects, especially for more complex cognitive tasks, but are nonetheless replicable between sessions within each subject. Presenting cognitive tasks in different modalities helped to separate distinct areas, and parcellation tools enhanced this functional dissection. In association with the acquisition of brain anatomies, cognitive profile and genotype of subjects, this protocol will serve to establish a hybrid database of hundreds of subjects suitable to study the range and causes of variation in the cerebral bases of numerous mental processes.
94 functional scans were acquired: 66 subjects performed one session, 13 subjects were recruited twice with an interval of several weeks, and 2 subjects performed 6 times the protocol in a thirty minutes session. In total, we obtained a database of 81 different functional scans and anatomy of French healthy subjects (35 males and 46 females; 24 ± 4.1 years old; 80 right handed and one left handed according to the Edinburgh inventory), 13 test/re-test situations and 2 multi-sessions to illustrate the evolution of results precision according to the number of scans acquired. All subjects gave their written informed consent. The study was approved by the regional ethical committee (Hopital de Bicêtre, France).
fMRI experimental design
Ten types of trials were mixed together and presented successively in a fixed order, which appeared random to the subject (Figure 1), with the E-Prime software (Psychology Software Tool, Inc.): 1-passive viewing of flashing horizontal checkerboards (10 trials), 2-passive viewing of flashing vertical checkerboards (10 trials), 3-pressing three times the left button with the left thumb button according to visual instructions (5 trials), 4-pressing three times the right button according to visual instruction (5 trials), 5-pressing three times the left button according to auditory instruction (5 trials), 6-pressing three times the right button according to auditory instruction (5 trials), 7-silently reading short visual sentences (10 trials), 8-listening to short sentences (10 trials), 9-solving silently visual subtraction problems (10 trials), 10-solving silently auditory subtraction problems (10 trials). 20 rest periods (black screen) were inserted in the sequence and served as null event for a better hemodynamic deconvolution. Subjects were briefly instructed before starting acquisition with a visual sequence projected in the scanner during a non-functional acquisition. A post-acquisition debriefing served to ensure that the subject correctly understood and performed the tasks.
Visual stimuli were displayed as four successive screens (250 ms) separated by 100 ms interval, and each composed of group of one to three words, resulting in 1.3 sec of visual stimulation. Auditory stimuli were digitally recorded by a male speaker (resolution of 16 bits and sampling frequency of 22.05 kHz) and had a similar duration (about 1.6 sec for motor instruction, 1.2–1.7 sec for sentences and 1.2–1.3 sec for subtraction). Calculation consisted of a one-digit number (range from 4 to 9) subtracted from a two-digit number (range from 10 to 19). (See Additional File 2 for a detailed description of the stimuli).
The experimental protocol was organized as a fast event-related paradigm. The 100 trials were presented in a fixed sequence with a stochastic SOA (2400 ms, 2700 ms, 3000 ms, 3300 ms or 3600 ms; mean SOA = 3 sec). This sequence was optimized for both statistical detection and hemodynamic response estimation using a Matlab script developed in the lab and inspired by the genetic algorithm of Wager and Nichols 2003 . One original feature of our procedure was to optimize the sequence according to more than one contrast. Eight contrasts of interest were selected: right vs. left hand action, vertical vs. horizontal checkerboards, auditory stimuli vs. rest, visual stimuli vs. rest, auditory calculations vs. auditory non-numerical stimuli, visual calculations vs. visual non-numerical stimuli, auditory stimuli vs. visual stimuli and visual stimuli vs. checkerboards.
Data acquisition and processing
Functional images were acquired on a 3T Brucker scanner using an EPI sequence (TR = 2400 ms, TE = 30 ms, matrix size = 64 × 64, FOV = 24 cm × 24 cm). Each volume consisted of 34 slices of 4 mm thickness. Anatomical T1 images were acquired with a spatial resolution of 1 × 1 × 1.2 mm. Data were pre-processed using SPM2 software  in Matlab7 environment according to the following procedure: slice timing, subject motion estimation and correction by realignment, coregistration of the anatomical image to the MNI template, spatial normalization of functional images (resampled voxel size = 3 × 3 × 3 mm) and smoothing (5 mm FWHM). Each voxel time series was fitted with a linear combination of the canonical hemodynamic response function and its temporal derivative. A temporal high pass filter was applied (cutoff 128 sec. and AR(1) whitening).
Individual contrast images were generated using SPM2. Individual conjunction maps were constructed to isolate voxels that were activated over a fixed threshold for multiple tasks: these maps were created by taking the min value of each of the contrasts considered for each voxel.
▪ Multi-sessions analysis
Two subjects performed the 5 min short protocol six times in a single session (all blocks used the same temporal design but four of the six included new stimuli to partially avoid habituation effects). We performed six SPM2 models, using a third of one block (i.e the first 30% of the trials for each condition of the first block), half of one block, one block, one and half blocks, two blocks, three blocks, four blocks, five blocks and six blocks, respectively. Statistics were collected from various maxima of interest listed from each subject's analysis: two frontal, four parietal and two temporal sites for calculation, four frontal, two inferior temporal, one fusiform and one parietal site for reading, the two left superior temporal maxima for auditory stimuli, and the two left motor area maxima for the right hand action.
For four contrasts, we conducted two SPM2 analyses for each subject using the first experimental block and all six blocks, respectively. Main maxima were listed for both analyses (p < 10-3, uncorrected, 10 voxels cluster extent for one block analysis; p < 10-3, corrected, 10 voxels cluster extent for the 6 blocks analysis reflecting the increased statistical power with number of trials). For each subject and contrast we calculated the proportion of peaks detected or missed when using only the first block compared to the overall session results, and the precision of their spatial location using Euclidian distance from final coordinates. Because these proportions may be task- and threshold-dependant, we first isolated local maxima as all the peaks listed with the SPM2 interface at the threshold noted above, separated by a minimal distance of 8 mm. We then selected as main maxima of a functional network peaks greater than two thirds of the highest t-value for the respective contrast.
The threshold-free receiver operating characteristic (ROC) approach compare the power to discriminate overall activated voxels identified by an SPM t-value map generated from trials of a 5 min session with the SPM thresholded t-value map obtained after a 30 min session (p < 10-3, corrected, 10 voxels cluster extent). ROC curves represent the plot of the true positive rate (sensitivity) against the false positive rate (1 – specificity) for different thresholds of the 5 min t-value map. The discriminative power (Dp) was estimated by calculating the area under the curve and may be interpreted as the probability that voxels classified as activated (or non-activated) in the 5 min map were also classified as activated (non-activated) in the 30 min sequence. In addition we calculated an inter-subject Dp by comparing the 5-min t-map of each subject with the 30-min t-map for the other subject. Similar patterns of activation across subjects would lead to comparable intra- and inter-subject Dpvalues.
▪ Intra-/inter-subjects variability
Inter-scan distance were calculated using the SPM2 Distance toolbox . It provides a global and objective measure of the relative 'distance' in time course and pattern of activation between all fMRI sessions for any contrast. A Multi-Dimensional Scaling (MDS) procedure allows a visual representation of the inter-session distance in a reduced 3-dimensional space which captures the greatest proportion of variance.
▪ Group analysis
Various methods of group analysis were applied to the first scan of our group of 81 subjects:
Random Effect Analysis (RFX)
Voxel-based RFX analysis were performed on the individual smoothed contrast images and conjunction images (5 mm FWHM) at a voxelwise threshold of p < 0.05 corrected for multiple comparison across the brain volume, with a minimal cluster extent of 5 voxels.
Brain Functional Landmark detection
We used an automatic identification method described by Thirion et al. . The list of landmarks is primarily derived from the RFX group analysis thresholded at p < 10-3, uncorrected. Reliability of the BFL is based on a recursive leave-one-out procedure and an agglomerative clustering of candidate individual local maxima.
Parcellation of fMRI datasets
An automatic parcellation method [62, 63] was used allowing up to 10 mm of variability in anatomical location of each parcel. It is based on a spectral clustering algorithm that delineates functionally homogeneous and spatially connected regions over the entire brain. A recent report of Thyreau et al., 2006  suggested 500 parcels per brain as spatial definition. Parcel-based RFX analysis was then performed on the set of individual parcellation, resulting in Z-value maps for each functional contrast.
Seghier ML, Lazeyras F, Pegna AJ, Annoni JM, Zimine I, Mayer E, Michel CM, Khateb A: Variability of fMRI activation during a phonological and semantic language task in healthy subjects. Hum Brain Mapp. 2004, 23: 140-155. 10.1002/hbm.20053.
Mechelli A, Penny WD, Price CJ, Gitelman D, Friston KJ: Effective connectivity and intersubject variability: Using a multisubject network to test differences and commonalities. NeuroImage. 2002, 17: 1459-1469. 10.1006/nimg.2002.1231.
Brett M, Johnsrude IS, Owen AM: The problem of functional localization in the human brain. Nat Rev Neurosci. 2002, 3: 243-249. 10.1038/nrn756.
Friston KJ, Holmes AP, Worsley KJ: How many constitute a study?. NeuroImage. 1999, 10: 1-5. 10.1006/nimg.1999.0439.
Thirion B, Pinel P, Meriaux S, Roche A, Dehaene S, Poline JB: Analysis of a large fMRI cohort: Statistical and methodological issues for group analysis. NeuroImage. 2007, 35 (1): 105-120. 10.1016/j.neuroimage.2006.11.054.
McGonigle DJ, Howseman AM, Athwal BS, Friston KF, Frackowiak RSJ, Holmes AP: Variability in fMRI: An examination of intersession differences. NeuroImage. 2000, 11: 708-734. 10.1006/nimg.2000.0562.
Smith SM, Beckmann CF, Ramnani N, Woolrich MW, Bannister PR, Jenkinson M, Matthews : Variability in fMRI: a re-examination of inter-session differences. Hum Brain Mapp. 2005, 24: 248-57. 10.1002/hbm.20080.
Juch H, Zimine I, Seghier ML, Lazeyras F, Fasel JHD: Anatomical variability of the lateral frontal lobe surface: implication for intersubject variability in language neuroimaging. NeuroImage. 2005, 24: 504-514. 10.1016/j.neuroimage.2004.08.037.
Shaywitz BA, Shaywitz SE, Pugh KR, Constable RT, Skudlarski P, Fulbright RK, Bronen RT, Fletcher JM, Shankweiler DP, Katz L, Gore JC: Sex differences in the functional organization of the brain for language. Nature. 1995, 373: 607-609. 10.1038/373607a0.
Grön G, Wunderlich AP, Spitzer M, Tomczak R, Riepe MW: Brain activation during human navigation: gender-different neural networks as substrate of performance. Nat Neurosci. 2000, 3: 404-408. 10.1038/73980.
Hariri AR, Weinberger DR: Imaging genomics. Brit Med Bull. 2003, 65: 259-270. 10.1093/bmb/65.1.259.
Fan J, Fossella J, Sommer T, Wu Y, Posner MI: Mapping the genetic variation of executive attention onto brain activity. Proc Natl Acad Sci USA. 2003, 100: 7406-7411. 10.1073/pnas.0732088100.
Egan MF, Goldberg TE, Kolachana BS, Callicott JH, Mazzanti CM, Straub RE, Goldman D, Weinberger DR: Effect of COMT Val108/158 Met genotype on frontal lobe function and risk for schizophrenia. Proc Natl Acad Sci USA. 2001, 98: 6917-6922. 10.1073/pnas.111134598.
Menon V, Kwon H, Eliez S, Taylor AK, Reiss AL: Functional brain activation during cognition is related to FMR1 gene expression. Brain Res. 2000, 8777: 367-370. 10.1016/S0006-8993(00)02617-2.
Menon V, Rivera SM, White CD, Eliez GH, Glover GH, Reiss AL: Functional optimization of arithmetic processing in perfect performers. Cogn Brain Res. 2000, 9: 343-345. 10.1016/S0926-6410(00)00010-0.
Burbaud P, Camus O, Guehl D, Bioulac B, Caillé JM, Allard M: Influence of cognitive strategies on the pattern of cortical activation during mental subtraction A functional imaging study in human subjects. Neurosci Letters. 2000, 287: 76-80. 10.1016/S0304-3940(00)01099-5.
Reichle ER, Carpenter PA, Just MA: The neural bases of strategy and skill in sentence-picture verification. Cog Psychol. 2000, 40: 261-295. 10.1006/cogp.2000.0733.
Nadeau SE, Williamson DJ, Crosson B, Gonzalez-Rothi LJ, Heilman KM: Functional imaging: heterogeneity in task strategy and functional anatomy and the case for individual analysis. Neuropsychiatry Neuropsychol Behav Neurol. 1998, 11: 83-96.
Sohn M, Goode A, Koedinger KR, Stenger VA, Fissel K, Carter CA, Anderson JR: Sex differences in the functional organization of the brain for language. Nat Neurosci. 2004, 7: 1193-1194. 10.1038/nn1337.
Delazer M, Domahs F, Bartha L, Brenneis C, Lochy A, Trieb T, Benke T: Learning complex arithmetic – an fMRI study. Cog Brain Res. 2003, 18: 76-88. 10.1016/j.cogbrainres.2003.09.005.
Dragansky B, Gaser C, Busch V, Schuierer G, Bogdahn U, May A: Changes in grey matter induced by training. Nature. 2004, 427: 311-312. 10.1038/427311a.
Toga AW, Thompson PM, Mega MS, Narr KL, Blanton RE: Probabilistic approaches for atlasing normal and disease-specific brain variability. Anat Embryol. 2001, 204: 267-282. 10.1007/s004290100198.
Van Essen DC: Windows on the brain: the emerging role of atlases and databases in neuroscience. Curr Opin Neurobiol. 2002, 12: 574-579. 10.1016/S0959-4388(02)00361-6.
Simon O, Kherif F, Flandin G, Poline JB, Rivière D, Mangin JF, Le Bihan D, Dehaene S: Automatized clustering and functional geometry of human parietofrontal networks for language, space, and number. NeuroImage. 2004, 23: 1192-1202. 10.1016/j.neuroimage.2004.09.023.
Price CJ: The anatomy of language: contributions from functional neuroimaging. J Anat. 2000, 197: 335-359. 10.1046/j.1469-7580.2000.19730335.x.
Binder JR, Frost JA, Hammeke TA, Bellgowan PSF, Springer JA, Kaufman JN, Possing ET: Human temporal lobe activation by speech and nonspeech sounds. Cereb Cortex. 2000, 10: 512-528. 10.1093/cercor/10.5.512.
Dehaene S, Cohen L, Sigman M, Vinckier F: The neural code for written words: a proposal. Trends Cogn Sci. 2005, 9: 335-341. 10.1016/j.tics.2005.05.004.
Rizzolatti G, Arbib MA: Language within our grasp. Trends Neurosci. 1998, 21: 188-194. 10.1016/S0166-2236(98)01260-0.
Butterworth B: The Mathematical Brain. 1999, Macmillan London
Gerstmann J: Fingeragnosie und isolierte agraphie; ein neues syndrom. Z Ges Neurol Psychiat. 1927, 108: 152-77.
Fernandez G, de Greiff A, von Oertzen J, Reuber M, Lun S, Klaver P, Ruhlmann J, Reul J, Elger CE: Language mapping in less than 15 minutes: Real-time functional MRI during routine clinical investigation. NeuroImage. 2001, 14: 585-594. 10.1006/nimg.2001.0854.
Stiers P, Peeters R, Lagae L, Van Hecke P, Sunaert S: Mapping multiple visual areas in the human brain with a short fMRI sequence. NeuroImage. 2005, 29: 74-89.
Drobyshevsky A, Baumann SB, Schneider W: A rapid fMRI task battery for mapping of visual motor cognitive and emotional function. NeuroImage. 2006, 31: 732-44. 10.1016/j.neuroimage.2005.12.016.
Tzourio-Mazoyer N, Landeau B, Papathanassiou D, Crivello F, Etard O, Delcroix N, Mazoyer B, Joliot M: Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. NeuroImage. 2002, 15: 273-289. 10.1006/nimg.2001.0978. [http://www.cyceron.fr/freeware]
Coulon O, Mangin JF, Poline JB, Zilbovicius M, Roumenov D, Samson Y, Frouin V, Bloch I: Structural Group Analysis of Functional Activation Maps. NeuroImage. 2000, 11: 767-782. 10.1006/nimg.2000.0580.
Mangin JF, Rivière D, Coulon O, Poupon C, Cachia A, Cointepas Y, Poline JB, Le Bihan D, Régis J, Papadopoulos-Orfanos D: Coordinate-based versus structural approaches to brain image analysis. Artif Intell Med. 2004, 30: 177-197. 10.1016/S0933-3657(03)00064-2.
Cohen L, Lehéricy S, Chochon F, Lemer C, Rivaud S, Dehaene S: Language-specific tuning of visual cortex? Functional properties of the Visual Word Form area. Brain. 2002, 125: 1054-1069. 10.1093/brain/awf094.
Price CJ, Thierry G, Griffiths T: Speech-specific auditory processing: where is it?. Trends Cogn Sci. 2005, 9: 272-276. 10.1016/j.tics.2005.03.009.
Vanderberg R, Nobre AC, Price CJ: The response of left temporal cortex to sentences. J Cogn Neurosci. 2002, 14: 550-560. 10.1162/08989290260045800.
Friederichi AG, Rüschemeyer SA, Hahne A, Fiebach CJ: The role of left inferior frontal and superior temporal cortex in sentence comprehension: Localizing syntactic and semantic processes. Cereb Cortex. 2003, 13: 170-177. 10.1093/cercor/13.2.170.
Just MA, Carpenter PA, Keller TA, Eddy WF, Thulborn KR: Brain activation modulated by sentence comprehension. Science. 1996, 274: 114-116. 10.1126/science.274.5284.114.
Chochon F, Cohen L, van de Moortele PF, Dehaene S: Differential contributions of the left and right inferior parietal lobules to number processing. J Cogn Neurosci. 1999, 11: 617-30. 10.1162/089892999563689.
Dehaene S, Piazza M, Pinel P, Cohen L: Three parietal circuits for number processing. Cogn Neuropsychol. 2003, 20: 487-506. 10.1080/02643290244000239.
Dehaene S, Cohen L: Towards an anatomical and functional model of number processing. Math Cogn. 1995, 1: 83-120.
Dehaene S: Varieties of numerical abilities. Cognition. 1992, 44: 1-42. 10.1016/0010-0277(92)90049-N.
Hubbard EM, Piazza M, Pinel P, Dehaene S: Interactions between number and space in parietal cortex. Nat Rev Neurosci. 2005, 6: 435-448. 10.1038/nrn1684.
Calvert GA: Crossmodal processing in the human brain: Insight from functional neuroimaging studies. Cereb Cortex. 2001, 11: 1110-1123. 10.1093/cercor/11.12.1110.
Beauchamp MS, Argall BD, Bodurka J, Duyn JH, Martin A: Unraveling multisensory integration: patchy organization within human STS multisensory cortex. Nat Neurosci. 2004, 7: 1190-1191. 10.1038/nn1333.
Simon O, Mangin JF, Cohen L, Le Bihan D, Dehaene S: Topographical layout of hand, eye, calculation, and language-related areas in the human parietal lobe. Neuron. 2002, 33: 475-487. 10.1016/S0896-6273(02)00575-5.
Pinel P, Piazza M, Le Bihan D, Dehaene S: Distributed and overlapping cerebral representations of number size and luminance during comparative judgments. Neuron. 2004, 41: 983-993. 10.1016/S0896-6273(04)00107-2.
Caspers S, Geyer S, Schleicher A, Mohlberg H, Amunts K, Zilles K: The human inferior parietal cortex: Cytoarchitectonic parcellation and interindividual variability. NeuroImage. 2006, 33 (2): 430-448. 10.1016/j.neuroimage.2006.06.054.
Hasnain MK, Fox PT, Woldorff MG: Hemispheric Asymmetry of sulcus-function correspondence: quantization and developmental implications. Hum Brain Mapp. 2006, 27: 277-287. 10.1002/hbm.20183.
Cachia A, Mangin JF, Rivière D, Kherif F, Boddaert N, Andrade A, Papadopoulos-Orfanos D, Poline JB, Bloch I, Zilbovicius M, Sonigo P, Brunelle F, Régis J: A primal sketch of the cortex mean curvature: a morphogenesis based approach to study the variability of the folding patterns. IEEE Trans Med Imag. 2003, 22 (6): 754-765. 10.1109/TMI.2003.814781.
Hasson U, Harel M, Levy I, Malach R: Large-Scale Mirror-Symmetry Organization of Human Occipito-Temporal Object Areas. Neuron. 2003, 37: 1027-1041. 10.1016/S0896-6273(03)00144-2.
Liegeois F, Baldeweg T, Connelly A, Gadian DG, Mishkin M, Vargha-Khadem F: Language fMRI abnormalities associated with FOXP2 gene mutation. Nat Neurosci. 2003, 6: 1230-1237. 10.1038/nn1138.
Molko N, Cachia A, Rivière D, Mangin JF, Bruandet M, Le Bihan D, Cohen L, Dehaene S: Functional and structural alterations of the intraparietal sulcus in a developmental dyscalculia of genetic origin. Neuron. 2003, 40: 847-858. 10.1016/S0896-6273(03)00670-6.
Leonard CM, Eckert MA, Kuldau JM: Exploiting human anatomical variability as a link between genome and cognomen. Genes Brain Behav. 2006, 5: 64-77. 10.1111/j.1601-183X.2006.00196.x.
Wager TD, Nichols TE: Optimization of Experimental Design in fMRI: A General Framework Using a Genetic Algorithm. NeuroImage. 2003, 18: 293-209. 10.1016/S1053-8119(02)00046-0.
SPM2 Software. [http://www.fil.ion.ucl.ac.uk/spm/software/spm2/]
Kherif F, Poline JB, Mériaux S, Banali H, Flandin G, Brett M: Group analysis in functional neuroimaging: selecting subjects using similarity measures. NeuroImage. 2003, 20: 2197-2208. 10.1016/j.neuroimage.2003.08.018. [http://groups.google.com/group/distance-toolbox]
Thirion B, Pinel P, Poline JB: Finding landmarks in the functional brain: detection and use for group characterization. Med Image Comput Comput Assist Interv Int Conf Med Image Comput Comput Assist Interv. 2005, 8 (Pt 2): 476-483.
Thirion B, Pinel P, Tucholka A, Roche A, Ciuciu P, Mangin JF, Poline JB: Improving sensitivity and reliability of fMRI group studies. IEEE TMI. 2007, 1256-69.
Thirion B, Flandin G, Pinel P, Roche A, Ciuciu P, Poline JB: Dealing with the shortcomings of spatial normalization: Multi-subject parcellation of fMRI datasets. Hum Brain Mapp. 2005, 27: 678-693. 10.1002/hbm.20210.
Thyreau B, Thirion G, Flandin G, Poline JB: Anatomo-functional description of the brain: a probabilistic approach. Proc 31th Proc IEEE ICASSP vol 5 Toulouse France. 2006, 1109-1112.
We thank the following investigators for sharing fMRI scanning time and subjects with us: Christophe Pallier, Sid Kouider, Manuela Piazza, Edward Hubbard and Anne-Dominique Devauchelle. These study also benefits to the helpful presence of Anne-Dominique Devauchelle and Mariagrazia Ranzini for subjects questioning. We also thank Catherine Oppenheim and Edward Hubbard for correction and discussion of the manuscript, and Jeremy Mattout for additional help in image analysis.
The project was conceived by PP, JBP and SD. JS was involved in optimization of the paradigm. Data were collected, processed and organized by PP and AJ. Inclusion of volunteers in the protocol was under the responsibility of DLB. Non-voxel based statistical methods were developed by BT, SM and JBP. All authors read and approved the manuscript.