# Shared genetic aetiology between cognitive performance and brain activations in language and math tasks

Yann Le Guen, Marie Amalric, Philippe Pinel, Christophe Pallier, Vincent Frouin

► To cite this version:

Yann Le Guen, Marie Amalric, Philippe Pinel, Christophe Pallier, Vincent Frouin. Shared genetic aetiology between cognitive performance and brain activations in language and math tasks. Scientific Reports, 2018, 8, pp.17624. 10.1038/s41598-018-35665-0. cea-01944702

HAL Id: cea-01944702

https://cea.hal.science/cea-01944702

Submitted on 4 Dec 2018

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>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

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

OPEN

# Shared genetic aetiology between cognitive performance and brain activations in language and math tasks

Yann Le Guen<fnref n="1" />, Marie Amalric², Philippe Pinel², Christophe Pallier² & Vincent Frouin<fnref n="1" />

Received: 5 January 2018

Accepted: 10 November 2018

Published online: 04 December 2018

Cognitive performance is highly heritable. However, little is known about common genetic influences on cognitive ability and brain activation when engaged in a cognitive task. The Human Connectome Project (HCP) offers a unique opportunity to study this shared genetic etiology with an extended pedigree of 785 individuals. To investigate this common genetic origin, we took advantage of the HCP dataset, which includes both language and mathematics activation tasks. Using the HCP multimodal parcellation, we identified areals in which inter-individual functional MRI (fMRI) activation variance was significantly explained by genetics. Then, we performed bivariate genetic analyses between the neural activations and behavioral scores, corresponding to the fMRI task accuracies, fluid intelligence, working memory and language performance. We observed that several parts of the language network along the superior temporal sulcus, as well as the angular gyrus belonging to the math processing network, are significantly genetically correlated with these indicators of cognitive performance. This shared genetic etiology provides insights into the brain areas where the human-specific genetic repertoire is expressed. Studying the association of polygenic risk scores, using variants associated with human cognitive ability and brain activation, would provide an opportunity to better understand where these variants are influential.

The genetic contribution to intellectual ability in human has been studied for almost a century<fnref n="1" />. These first studies on twins and families constitute the starting point of the behavioral genetic field. One of the first goal was to estimate the heritability of behavioral traits, meaning the proportion explained by genetics in the variance of an observed trait. These heritability studies were historically conducted in dizygotic and monozygotic twins for which a sensible model known as Falconer formula² links the variance of a trait to the a priori genetic and environmental shares, conveniently determined in twins. Among the behavioral traits, language and math functions in humans have been extensively studied in fundamental neuroscience as distinctive abilities of human lineage. They are frequently assessed through standard behavioral tests as regular phenotypes but also with neuroimaging to provide brain characteristics, named endophenotypes³,⁴. Those two kinds of phenotype are used jointly as a way to classify the broad behavioral symptoms of language impairments into stable characteristics that in turn are candidates to search for potential associations with either medical treatment responses or genetic profiles⁵,⁶. Structural properties observed using magnetic resonance imaging (MRI) or activations observed with functional MRI (fMRI) have been used to produce such endophenotypes⁷. These can reveal differences between control and disease groups in language-specific regions⁸ or distinguish disorder subtypes such as grammatical-SLI (specific language impairment)⁹. In this study, the endophenotypes correspond to fMRI activations, which are hypothesized to reflect brain spatial display of underlying molecular mechanisms when a subject is tested. This raises the question of their potential common genetic roots with cognitive ability assessed by tests. Imaging-genetics resources, such as the Human Connectome Project (HCP) provide an unprecedented opportunity to study the variability of behavioral endophenotypes along with behavioral tests in control subjects, as well as to determine their potential heritability or association with genetics. This dataset provides a unique opportunity to link in one study issues addressed by the genetics, neuroimaging and psychology communities because it contains: pedigree

<footnote n="1">¹Neurospin, Institut Joliot, CEA, Université Paris-Saclay, Gif-sur-Yvette, France. ²Cognitive Neuroimaging Unit, U992, INSERM, Neurospin, Institut Joliot, CEA, Université Paris-Saclay, Gif-sur-Yvette, France. Correspondence and requests for materials should be addressed to Y.L.G. (email: yann.leguen@cea.fr)</footnote>
information, fMRI images and behavioral tests. There are existing bridges between these fields, but to the best of our knowledge data is rarely available at the same time. The rationale of this study is thus to look into the heritability of fMRI activations and behavioral scores separately and then to investigate the shared genetics roots between the two.

Following up on these ideas, we first proposed to study the genetic influence involved in fMRI activation differences among typically developed individuals. We used the pedigree data from the HCP language comprehension and verbal math fMRI tasks. These tasks recruit regions directly implicated in brain disorders, such as Broca's area in SLI¹⁰, the angular gyrus in developmental dyslexia³ and the intraparietal sulcus in dyscalculia⁴. A few studies have already attempted to estimate the narrow sense heritability of brain activations for various tasks. They notably include digit and n-back working memory¹¹,¹², visual math subtraction¹³, and stimuli such as written words, faces and spoken language¹⁴. However, these studies had relatively small sample sizes for reliably detecting heritability estimates ranging between 25 and 50%. In addition to including a larger sample size, the HCP data were processed using state of the art methods, providing 2-mm isotropic resolution and finer inter-individual registration. In particular, the so-called grayordinate activations are computed on the surface of the cortex¹⁵ for each individual, and inter-subject fMRI alignment is performed using areal-feature-based registration¹⁶. The grayordinate approach refers to fMRI analyses performed on the cortical surface, as opposed to a volume-based approach. We decomposed our analysis onto the HCP's multimodal parcellation of the human cerebral cortex¹⁷, which enabled us to map the genetic influence on fMRI activations on a very fine scale.

Furthermore, it is known that neural activation endophenotypes from MRI may reflect not only impairments in language but also differences in cognitive scores. For example, weaker left-lateralizations have been reported for some developmental language disorders¹⁸, and in normal populations, fMRI activations during simple tasks correlate with various cognitive scores. Notably, single digit calculation fMRI activations are predictive of high school math scores¹⁹, and the fronto-parietal functional connectivity in children performing a task that required them to match Arabic numbers to an array of dots correlated with their score on a standardized math test²⁰. Using HCP data and in line with these approaches, we show in this paper how variations in language related fMRI activations correlate with cognitive abilities assessed by the median reaction time (RT), average accuracy and difficulty level during the HCP language and math tasks. Remarkably, recent studies have shown that, beyond the age-related heritability of general cognitive ability²¹,²² and of various indicators of academic performance²³, these scores are highly pleiotropic²¹,²⁴ [pleiotropy occurs when one gene regulates one or more phenotypic traits].

This raises the question of the potential pleiotropy between neural activations and cognitive abilities. Thus, as a second contribution, we studied the shared genetic variance of fMRI activations and cognitive performance scores measured during the MRI session or behavioral scores acquired independently from the task. We studied behavioral variables measured by the HCP using standardized tests from the National Institute of Health (NIH): fluid intelligence, working memory, and language assessments such as vocabulary comprehension and oral reading decoding. Details of these variables and their heritability estimates can be found in Table S1, and how well they correlate phenotypically and genetically with the behavioral scores measured during the task is reported in Table S2.

Recent genome wide association studies have unveiled new loci and genes influencing human cognitive performance (e.g. human intelligence²⁵, general cognitive function²⁶,²⁷ and educational attainment²⁷) and possibly intelligence as a construct in differential psychology²¹. However, for these human-specific characteristics, little is known about the underlying integration mechanism of molecular functions or the brain areas where they are most influential. The shared genetic etiology investigated in this work provides new perspectives to decipher the basis of cognitive abilities such as language in humans. This study had two major aims: (1) to estimate the heritability of fMRI activations during story comprehension and math tasks; and (2) to determine the shared genetic etiology between these activations and cognitive performance.

## Results

**Task fMRI Activations in MATH and STORY tasks.** Figure 1 shows the activations for MATH (vs the intercept of the general linear model (GLM) being considered as baseline), STORY and the contrast STORY - MATH. The intercept reflects the mean of the residual BOLD time series after removing variance explained by all other regressors. Both tasks show clear activations in the planum temporale and Heschl's gyrus area, reflecting the fact that the stimuli were presented in the auditory modality. The MATH task, in which participants were requested to perform addition and subtraction, activates areas traditionally implicated in mathematical calculations, that is, the intraparietal sulcus, the middle frontal and the inferior temporal regions²⁸,²⁹. The story listening task activates the language understanding network, encompassing bilateral temporal regions and left frontal regions³⁰,³¹. As expected, the group activations for the STORY task are more left lateralized, notably in the left posterior superior temporal and inferior frontal regions, which correspond to Wernicke's and Broca's areas, respectively. Moreover, regions implicated in inhibition networks are also activated by these tasks³²,³³, notably the middle frontal gyrus in the math task and the medial prefrontal cortex, implicated in motivation and execution, and above the anterior cingulate cortex, controlling selective attention³⁴,³⁵. In addition, both tasks activate complementary networks; in particular, the math task deactivates the semantic and episodic memory processes, known as the default mode network, which is also active in resting or passive states³⁶. This last remark makes the STORY-MATH contrast particularly relevant for studying the genetic influence on activation specifically elicited by math and story tasks.

**Univariate Genetic Analyses.** We performed a cortex-wise heritability analysis on the median activation (β-value) in the 360 areals of the HCP multi-modal parcellation. After stringent Bonferroni correction (p < 0.05/360 ≈ 1.4·10⁻⁴), we found 54 regions whose activations during the MATH task are heritable and 46 regions for the STORY task. These results are summarized in Fig. 2 and heritability estimates are included in
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Group average activations for the HCP language tasks, including MATH (a) and STORY (b) blocks, and the STORY-MATH contrast (c). Group maps are shown with a lower threshold of z = ±10 and saturation from z = 30 to introduce the main areas activated by the tasks. Due to the large number of subjects, the associated p-values are significant; we arbitrarily set the thresholds to emphasize the regions that are known to be recruited by these tasks.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Heritability estimates for the activations of the MATH (a) and STORY (c) tasks, and their associated p-values (respectively b and d). Only the estimates significant after correction (p < 0.05/360, with 180 areals in each hemisphere) are displayed. Activations correspond to the median parameter estimate (β) in each areal of the HCP multimodal parcellation.</figcaption></figure>

Tables S3–S6. The details and names of the areals can be found in the Supplementary Information of the paper describing the multimodal parcellation of the human cerebral cortex¹⁷. For the MATH (resp. STORY) task, the heritability estimates range from 0.23 to 0.45, with the maximum in the left “Area PGp” corresponding to the angular gyrus (resp. 0.22 to 0.55, with the maximum in the left “PeriSylvian Language Area”). In addition, we performed heritability analysis using the median z-stat value in each areal instead of the median parameter estimate (β-value) and obtained similar results (Fig. S1).
<table><caption>**Table 1.** Heritability estimates for the behavioral scores associated with the tasks. Language accuracy and reaction time (RT) correspond to the average of the respective MATH and STORY variables. The p-values associated with the covariates related to age and sex, ethnic group and education level are also displayed.</caption><thead><tr><th rowspan="2">Trait</th><th rowspan="2">h<sup>2</sup> ± SE (p)</th><th colspan="2">Age</th><th colspan="2">Sex</th><th colspan="2">Hispanic</th><th rowspan="2">Educ</th><th rowspan="2">h<sup>2</sup>cov%</th></tr><tr><th>p-val</th><th>Age<sup>2</sup></th><th>Age<sup>2</sup></th><th>p-val</th><th>Age<sup>2</sup></th><th>p-val</th></tr></thead><tbody><tr><td>Language Accuracy</td><td>0.34 ± 0.06 (2.3·10<sup>-8</sup>)</td><td>0.49</td><td>0.75</td><td>0.01</td><td>0.96</td><td>0.22</td><td>0.97</td><td>1.2·10<sup>-7</sup></td><td>7.2</td></tr><tr><td>Language RT</td><td>0.22 ± 0.07 (8.4·10<sup>-4</sup>)</td><td>0.44</td><td>0.91</td><td>0.69</td><td>0.85</td><td>0.95</td><td>0.71</td><td>0.04</td><td>0.5</td></tr><tr><td>Math Accuracy</td><td>0.4 ± 0.06 (1.6·10<sup>-10</sup>)</td><td>0.5</td><td>0.96</td><td>9.2·10<sup>-4</sup></td><td>0.81</td><td>0.39</td><td>0.76</td><td>1.9·10<sup>-6</sup></td><td>7.2</td></tr><tr><td>Math Difficulty Level</td><td>0.33 ± 0.07 (1.0·10<sup>-6</sup>)</td><td>0.53</td><td>0.17</td><td>0.02</td><td>0.61</td><td>0.05</td><td>0.42</td><td>8.7·10<sup>-8</sup></td><td>7.9</td></tr><tr><td>Math Median RT</td><td>0.17 ± 0.07 (7.4·10<sup>-3</sup>)</td><td>0.6</td><td>0.41</td><td>0.78</td><td>0.82</td><td>0.69</td><td>0.26</td><td>0.01</td><td>1.1</td></tr><tr><td>Story Accuracy</td><td>0.18 ± 0.06 (1.6·10<sup>-3</sup>)</td><td>0.78</td><td>0.48</td><td>0.86</td><td>0.76</td><td>0.23</td><td>0.61</td><td>2.2·10<sup>-3</sup></td><td>1.4</td></tr><tr><td>Story Difficulty Level</td><td>0.33 ± 0.07 (1.3·10<sup>-6</sup>)</td><td>0.66</td><td>0.43</td><td>0.81</td><td>0.44</td><td>0.28</td><td>0.91</td><td>0.11</td><td>0.0</td></tr><tr><td>Story Median RT</td><td>0.2 ± 0.07 (1.4·10<sup>-3</sup>)</td><td>0.41</td><td>0.63</td><td>0.39</td><td>0.59</td><td>0.68</td><td>0.45</td><td>0.39</td><td>0.0</td></tr></tbody></table>

The univariate genetic analysis of the activations associated with verbal math emphasizes mainly areals spanning the math network, including the intraparietal sulcus, middle frontal, inferior temporal and angular gyri. The analysis of activations associated with story comprehension distinctively underlines regions of the language network as bilaterally heritable. Among these regions are the superior temporal sulcus dorsal and ventral parts, Brodmann area (BA) 47 in Broca's area, and the middle frontal gyrus at the junction with the precentral sulcus. Interestingly, the heritability networks of the MATH and STORY tasks overlap very little except in the auditory cortex, around the planum temporale, in the frontal cortex (BA 8), and in the inferior temporal region.

Table 1 presents the heritability estimates of the behavioral scores gathered during the MRI scans. The global accuracy on the HCP language tasks, averaging the scores in the MATH and STORY tasks, was significantly heritable, with h² = 0.34, close to the traditionally high heritability estimate of cognitive performance³⁷,³⁸. The heritability estimates for the median reaction time (RT) were approximately 0.2. Furthermore, RT Story and RT Math were significantly correlated (phenotypic correlation: ρₚ = 0.34, genetic correlation: ρₗ = 0.45, Table S2). Regarding the accuracy and average difficulty level of the HCP language task, we observed that the MATH task variables have higher heritability estimates than those of the STORY task. This result might indicate a higher genetic influence on performance during simple arithmetic tasks than during language comprehension. However, this result needs to be considered in light of the different distribution patterns of MATH and STORY accuracies. The STORY accuracy reported by the HCP displays discrete values and might not be sufficiently informative (Fig. S2). Table S2 underlines a significant correlation between math and story accuracies (ρₚ = 0.15, ρₗ = 0.27). The discrete distribution of STORY accuracy likely occurs because each story lasted approximately 20–30 s, few story questions were presented to the subjects, and most subjects tended to choose the correct answer in the two-alternative forced-choice question.

**Bivariate Genetic Analyses.** We performed bivariate genetic analyses to quantify the shared genetic influence between intellectual performance, represented by the behavioral measures, and the neural activation in each areal. The genetic correlation estimates are usually subject to substantial sampling errors and therefore inaccurate. The large sample size of the HCP offers the opportunity to reduce the standard errors. The distribution of STORY accuracy values is concentrated on a small number of values (Fig. S2), thus, we chose to use the average of the STORY and MATH accuracies as the behavioral score to characterize the individual performance. Thus, the results presented here concern the relationship between this average score and the activations or deactivations revealed by the contrast STORY-MATH (activation map shown in Fig. 1c). As a first step of our analysis and to filter out the areals for which the neural activation was not significantly correlated with the behavioral score, we computed the phenotypic correlation between these two variables in each areal of the HCP multi-modal parcelation. Figure 3(a,b) summarizes the phenotypic correlation and associated p-values, for areals significant after Bonferroni correction (p < 0.05/360). The language network is clearly encompassed along the left superior temporal sulcus (STS) and Broca's area, as well as in the anterior part of the right STS. The activations in the angular gyrus (area PGp), supporting the manipulation of numbers in verbal form³⁹, were also significantly correlated with the behavioral scores.

Among the 360 areals of the HCP multimodal parcellation, 39 (resp. 38) were significantly phenotypically correlated and were kept for the bivariate analysis in the right (resp. left) hemisphere. The shared genetic variance estimates for these areals are presented Fig. 3(c,d) (p < 0.05 without correction), and detailed values can be found in Tables S7 and S8. With stringent Bonferroni correction the p-value threshold for ρₗ is p < 0.05/(39 + 38) ≈ 6.5·10⁻⁴. Among the areals with significant shared genetic variance, we found the left anterior ventral insular area (AVI, ρₗ = 0.61) and the right angular gyrus (PGp, ρₗ = −0.40). Noticeably, in the left hemisphere areals, parts of the language network in the posterior STS had activations that shared significant genetic variance with language accuracy. These include the posterior ventral (STSv posterior, ρₗ = 0.47) and dorsal (STSd posterior, ρₗ = 0.45) parts of the STS, adjacent to the auditory 5 complex area (A5, ρₗ = 0.54), the perisylvian language area (PSL, ρₗ = 0.47) and the temporo-intraparietal junction (PGi, ρₗ = 0.54). On the left hemisphere internal face, we also found the superior frontal language area (SFL, ρₗ = 0.61) and, adjacent to this areal the Brodmann 8 decomposed into medial (8Bm, ρₗ = 0.75) and lateral (8Bl, ρₗ = 0.73) parts. Additionally, we noted two right hemisphere regions implicated in language processing and significantly genetically correlated with the fMRI task average score: the temporal pole (area TG dorsal, TGd, ρₗ = 0.65) and the lateral part of Brodmann area 47 (471, ρₗ = 0.51).
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Bivariate genetic analysis results between HCP LANGUAGE task accuracy and activation for the STORY-MATH contrast in each areal. (a) After strict Bonferroni correction (p < 0.05/360), significant phenotypic correlations between the language task accuracy (average of story and math accuracies) and the median activation of the contrast STORY-MATH in each areal (b) with their associated p-values. (c) Proportion of variability due to shared genetic effects with (d) their associated uncorrected p-values < 0.05.</figcaption></figure>

The latter is adjacent to Brodmann areas 44 and 45 in the inferior frontal, which are connected through the arcuate fasciculus with the language temporal regions.

Additionally, we extended our analysis to behavioral variables measured by the HCP following a standardized NIH protocol. Among these, we selected the variables that are most likely to reflect cognitive performance. Then, we estimated their heritability (Table S1) and correlations with the behavioral scores measured during the fMRI task (Table S2). In this set of variables, fluid intelligence (heritability: $h^2$ = 0.43, correlations with language accuracy: $\rho_p$ = 0.36, $\rho_g$ = 0.61), working memory ($h^2$ = 0.52, $\rho_p$ = 0.34, $\rho_g$ = 0.50), vocabulary comprehension ($h^2$ = 0.64, $\rho_p$ = 0.40, $\rho_g$ = 0.57) and oral reading decoding ($h^2$ = 0.67, $\rho_p$ = 0.46, $\rho_g$ = 0.67) were the ones with the highest heritability estimates and correlations with the average accuracy of the two fMRI tasks. Thus, we performed a bivariate genetic analysis between the STORY-MATH activations (difference between $\beta_{STORY}$ and $\beta_{MATH}$) and these four variables. Regardless of whether one considers the NIH scores or the ones directly related to the fMRI tasks, the study of shared genetic influence with the median activation yields approximately the same set of regions (Figs 4 and 5). This observation reinforces our claim that these regions have common genetic roots with the parts of general cognitive performance accounted for by the four cognitive variables under scrutiny, namely fluid intelligence, working memory, vocabulary comprehension and reading decoding.

## Discussion

In this paper, we have shown that brain activation pattern in the language and math networks are heritable. Additionally, we highlighted a particular set of regions along the superior temporal sulcus and in the inferior frontal whose activations share a common genetic basis with some aspects of general cognitive ability, assessed through fMRI task accuracy and behavioral scores.

Previous studies that have estimated the heritability of fMRI task activation had low sample sizes for significantly estimating heritably values ranging from 25% to 50%. The previous samples included 30 subjects (10 triplets of male monozygotic (MZ) twins with one additional brother)¹¹, 64 subjects (19 MZ and 13 dizygotic (DZ) pairs)¹³,¹⁴ or 319 subjects (75 MZ and 66 DZ pairs, 37 unpaired)¹². We must emphasize that these results correspond to fMRI activations associated with verbal math and semantic comprehension tasks. Thus, regions not recruited by the tasks cannot be found to be significantly correlated with cognitive ability in our case, because activations in these regions are incoherent across individuals. Notably, the visual word form area, related to literacy, is not activated in our oral tasks because they did not require word reading.

Furthermore, combining data from various cohorts is unfeasible, because neural activations from different tasks are not comparable when estimating inter-individual variance. This highlights the necessity of utilizing large cohorts with standardized fMRI protocols to perform such genetic analyses. To our knowledge, this is the first study to address the heritability of fMRI activation cortex-wise on a multimodal parcellation of the human cerebral cortex. Our results confirm the genetic influence on the formation of neural circuits implicated in language⁴⁰ and math¹³. Using the HCP
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Significant phenotypic correlations between the grayordinate activations of the STORY-MATH contrast and the NIH behavioral scores. (a) Fluid intelligence (PMAT24_A_CR). (b) Working memory (ListSort). (c) Vocabulary comprehension (PicVocab). (d) Reading decoding (ReadEng). Associated p-values (*p* < 0.05/360, Bonferroni correction) can be found Fig. S3.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Shared genetic variance (absolute value) between the grayordinate activations of the STORY-MATH contrast and the NIH behavioral scores. (a) Fluid intelligence (PMAT24_A_CR). (b) Working memory (ListSort). (c) Vocabulary comprehension (PicVocab). (d) Reading decoding (ReadEng). Associated p-values (*p* < 0.05, uncorrected) can be found in Fig. S4. Genetic correlation was investigated only for areals that were significantly phenotypically correlated (Fig. 4).</figcaption></figure>
fine scale parcellation<sup>17</sup> allowed us, for instance, to distinguish the genetic effects on the temporo-parietal junction implicated in language<sup>30</sup> (area PGI) and on the adjacent angular gyrus (area PGp), which is particularly involved in the manipulation of numbers in verbal form<sup>39</sup>. Indeed, these two areas, part of BA 39, present different cytoarchitec- tonic properties, such as a slightly broader layer II for PGI<sup>41</sup>, which might explain their involvement in different tasks. In a previous work, Pinel and Dehaene also found the left angular gyrus and the posterior superior parietal lobule bilaterally to be heritable<sup>13</sup>. Adding to these observations, our results underline a left hemisphere intraparietal sulcus specificity, with more heritable areals and slightly higher heritability compared to the right for the MATH task. This finding is consistent with results reported by Vogel and colleagues demonstrating a correlation of activations in the left intraparietal sulcus modulated by age, which was not observed in the right intraparietal sulcus<sup>42</sup>. Heritability represents the proportion of observed inter-individual phenotypic variance that is explained by genetics. Thus, it might be that inter-individual variance is not sufficiently pronounced in the right hemisphere, whereas activations have evolved over one's lifetime in the left hemisphere. Overall, the heritability maps for the STORY and MATH tasks pinpoint regions known to be disrupted in neurodevelopmental disorders. For instance, the inferior frontal area and the temporo-parietal junction activations are impaired in developmental dyslexia<sup>3,43</sup>, and the intraparietal region activations are less modulated by the numerical distance between two numbers being compared in developmen- tal dyscalculia<sup>4,44,45</sup>. Highlighted areas might provide new insights into brain regions where normal gene expression might be disrupted, leading to brain dysfunction and neurodevelopmental disorders. Frequently replicated genes associated with neurobehavioral disorders, such as developmental dyslexia or SLI, likely play such a role in structural brain maturation by interfering with neuronal migration and neurite growth<sup>6</sup>.

Several studies have already described some phenotypic correlations between cognitive abilities and neural activations in language<sup>46-48</sup> and math<sup>19,20,49,50</sup>. Our study replicates these observations, notably the correlation with language processing regions, including Broca's area and the posterior superior temporal gyrus<sup>46</sup>. Moreover, we estimated the genetic proportion in these phenotypic correlations. Hence, we demonstrated a shared genetic etiology between brain activations and cognitive performance, assessed in our study by the following tests: fluid intelligence, working memory, vocabulary comprehension and reading decoding. Interestingly, in the right hemi- sphere, mainly the anterior STS was found to be genetically correlated with the language task accuracy. This result seems consistent with the hypothesized role of the right anterior STS in the processing of figurative language, likely involved in the Aesop's fable metaphors presented to the subjects<sup>51</sup>.

The observed genetic correlations shed light on the genetic links between cognitive performance and activation level in cognitive task-related fMRI. These links might be related to the development and maturation of myelin, enhancing brain connectivity. In children with difficulties processing syntactically complex sentences, arcuate fasciculus maturation was incomplete compared to adults<sup>52,53</sup>. Thus, we could look for additive genetic effects implicated in the various levels of fiber tract maturation, which improves brain connectivity and efficiency. Indeed, Skeide and colleagues reported an example of such a genetic risk variant for dyslexia. They showed that this variant is related to the functional connectivity of left fronto-temporal phonological processing areas during the resting state<sup>54</sup>. Similarly, children with higher arithmetic scores present a more mature response modulation in their left intraparietal lobe<sup>49</sup>. Our study suggests that a proportion of the observed inter-individual variance in cognitive performance partly results from the same additive genetic effects as those contributing to brain activation variance. The moderate shared genetic basis suggests that a crucial interaction occurs between the environment and gene networks to enable the brain to develop to its full potential.

In this study, we could not accurately model the shared environment between siblings because HCP did not provide household information to model the shared environment. However, we pinpointed in the method section that we used the education level to model environmental differences in family socioeconomic status to mitigate the absence of common environment information. This choice is likely conservative and our heritability estimates and shared genetic are likely underestimated. Indeed, this model is conservative because the number of year of education shares genetic roots with the general cognitive ability<sup>55</sup>. In addition, the genotyping material from HCP was not available at the time of our study, thus we did not study gene × environment interaction.

Recently, with the emergence of large cohorts, such as UK Biobank, new loci and genes influencing human cognitive ability have been discovered<sup>25-27</sup>. However, little is known about how these genes contribute to this human-specific trait. Our study pinpoints brain regions where activations genetically correlate with global cognition scores. These regions might help elucidate the mechanism in which these genes are implicated. When the HCP genotyping data are released, a polygenic score of these newly discovered variants could be used to determine the explained proportion of the neural activation variance in these regions.

## Material and Methods

**Subjects.** This study utilized the dataset of the Human Connectome Project (HCP). The HCP scans and data were released in April 2017 (humanconnectome.org). The details of the release are available in the HCP reference manual. In this project, 1046 subjects aged between 22 and 37 years old (μ ± σ = 28.8 ± 3.7 years) completed the fMRI language task in the HCP S1200 release. To avoid population stratification, we only included the 785 Caucasian individuals (372/413 M/F) that were classified as race = "White" by the HCP; among these, the ethnicity of 69 was "Hispanic/Latino". This subgroup of the HCP contains 178 twin pairs (117 monozygotic twins (MZ) with 103 siblings and 61 dizygotic twins (DZ) with 61 siblings and 1 half sibling), 203 siblings, 1 half sibling and 60 unpaired individuals. The unpaired individuals did not contribute to the genetic parameter estimation but allowed a more accurate estimation of mean and variance effects. Subjects were chosen by the HCP consortium to represent healthy adults beyond the age of major neurodevelopmental changes and before the onset of neurodegenerative changes<sup>56</sup>. They underwent a battery of tests to determine if they met the inclusion/exclusion criteria of the HCP<sup>56</sup>. All subjects provided written informed consent on forms approved by the Institutional Review Board of Washington University. All of the following methods were carried out in accordance with relevant guidelines and regulations. All experimental MRI protocols were approved by the Institutional Review Board of Washington University.
**Image acquisition and processing.** MR images were acquired using a 3T Connectome Scanner, adapted from Siemens Skyra, housed at Washington University in St Louis, using a 32-channel head coil. T1-weighted images with 256 slices per slab were acquired with the three-dimensional magnetization-prepared rapid gradient echo (3D-MPRAGE) sequence: TR = 2400 ms, TE = 2.14 ms, TI = 1000 ms, flip angle = 8°, FOV = 224 × 224 mm, and resolution 0.7 mm isotropic. T2-weighted images, 256 slices per slab, were acquired with a 3D-T2SPACE sequence: TR = 3200 ms, TE = 565 ms, variable flip angle, FOV = 224 × 224 mm, and resolution 0.7 mm isotropic. fMRI data acquisition parameters were as follows: TR = 720 ms, TE = 33.1 ms, flip angle = 52 deg, BW = 2290 Hz/Px, in-plane FOV = 208 × 180 mm, 72 slices, and 2.0 mm isotropic voxels, with a multi-band acceleration factor of 8. Two runs of each task were acquired, one with right-to-left and the other with left-to-right phase encoding 2; each run interleaved 4 blocks of a story task with 4 blocks of a math task. The lengths of the blocks varied (average of approximately 30 seconds), but the task was designed so that the math task blocks matched the length of the story task blocks, with some additional math trials at the end of the task to complete the 3:57 (min:sec) run.

The details of the HCP data analysis pipelines are described elsewhere<sup>15,57</sup>. Briefly, they are primarily built using tools from FSL<sup>58</sup> and Freesurfer<sup>59</sup>. The HCP fMRI Volume pipeline generates "minimally preprocessed" 4D time series that include gradient unwarping, motion correction, fieldmap-based EPI distortion correction, brain-boundary- based registration of EPI to structural T1-weighted scans, non-linear (FNIRT) registration into MNI152 space, and grand-mean intensity normalization<sup>57</sup>. For the S500 release, two smoothing approaches were chosen by the HCP: volume-based smoothing or smoothing constrained to the cortical surface and subcortical gray-matter parcels. For the former, standard FSL tools can be applied for analysis, while for the latter, the HCP adapted these tools to this the 'grayordinate' approach<sup>15,57</sup>. The grayordinate approach refers to fMRI analyses performed on the cortical surface, as opposed to a volume-based approach. This is more accurate spatially because activation occurs in gray, not white, matter. Unconstrained volume-based smoothing causes blurring effects by mixing signals from cortex regions adjacent in volume but not on the surface. For these reasons, our study analyses were carried out on the surface of the cortex.

The HCP fMRI Surface pipeline brings the time series from the volume into the CIFTI grayordinate standard space. This is accomplished by mapping the voxels within the cortical gray matter ribbon onto the native cortical surface, transforming them according to the surface registration onto the 32k Conte69 mesh, and mapping the set of subcortical gray matter voxels from each subcortical parcel in each individual to a standard set of voxels in each atlas parcel. The result is a standard set of grayordinates in every subject (i.e., the same number in each subject, with spatial correspondence) with 2 mm average surface vertex and subcortical volume voxel spacing. These data are smoothed with surface and parcel constrained smoothing of 2 mm FWHM (full width half maximum) to regularize the mapping process<sup>57</sup>.

**The language task.** The HCP language task was developed by Binder and colleagues<sup>36</sup>. In the story blocks, participants were presented with brief auditory stories adapted from Aesop's fables, followed by a 2-alternative forced-choice question to check the participants' understanding of the story topic. The example provided in the original paper is "*For example, after a story about an eagle that saves a man who had done him a favor, participants were asked, "Was that about revenge or reciprocity?"*". In the math blocks, participants were also presented auditory series of addition and subtraction (e.g., "fourteen plus twelve"), followed by "equals" and then two choices (e.g., "twenty-nine or twenty-six"). To ensure similar level of difficulty across participants, math trials automatically adapted to the participants responses. As shown by Binder and colleagues<sup>36</sup>, the story and math trials were well matched in terms of duration, auditory and phonological input, and attention demand. Furthermore, they were likely to elicit distinct brain activation – on the one hand, anterior temporal lobes classically involved in semantic processing, and parietal cortex on the other hand, classically involved in numerical processing, thus spanning a broad set of regions involved in conceptual semantic processing.

**HCP task fMRI analysis.** The analysis of fMRI data was carried out by the HCP consortium and we describe briefly their pipeline<sup>15</sup>. The Story predictor covered the variable duration of a short story, question, and response period (~30 s). The Math predictor covered the duration of a set of math questions designed to roughly match the duration of the story blocks. The grayordinate data for individual task runs were processed in a level 1 analysis. Activity estimates were computed for the preprocessed functional time series from each run using a general linear model (GLM) implemented in FSL's FILM (FMRIB's Improved Linear Model with autocorrelation correction)<sup>60</sup>. Predictors were convolved with a double gamma "canonical" hemodynamic response function<sup>61</sup> to generate the main model regressors. The two runs for each task and subject were then combined in a level 2 fixed-effects analysis<sup>15</sup>, which we used as our phenotype. Fixed-effects analyses were conducted using FEAT (fMRI Expert Analysis Tool) to estimate the average effects across runs within-subjects, and then mixed-effects analyses treating subjects as random effects were conducted using FLAME (FMRIB's Local Analysis of Mixed Effects) to estimate the average effects of interest for the group third-level analysis.

**Phenotype definitions.** To define our phenotypes, we consider separately the regression analyses on STORY and MATH tasks, and the contrast STORY-MATH. We used the beta values (pe1.dtseries.nii files) of the results of the level 2 analysis, which essentially average the level 1, i.e., the individual, runs. The contrasts were defined by the HCP in level 1 and averaged for level 2: thus, the grayordinate values of the beta and contrast values (cope1.dtseries.nii) are identical in this case, as they did not define any "new" contrasts specifically at level 2. Therefore, we could have used the cope1.dtseries.nii.files with no difference in results. We used the MSMAll registered the functional analysis results from HCP and the HCP multimodal parcellation<sup>17</sup>. We analyzed each of the 180 areals separately. We computed the median beta values in each areal for both hemispheres. These phenotypes constitute our proxy to estimate the activation in each part of the brain.
Moreover, we also included in our phenotypes the accuracy, reaction time and average difficulty level for the MATH and STORY tasks. We called these “behavioral scores”, as opposed to the grayordinate activation phenotypes previously defined.

The STORY-MATH contrast allows to cancel out: (i) the mean of the residual BOLD signal, which reflects the mean of the random error, (ii) signal in areas activated by both tasks, for example the primary auditory area. In addition, as emphasized by Binder and colleagues in the paper introducing this contrast³⁶, the story and math tasks activate complementary networks. In particular, the math task deactivates the semantic and episodic memory processes, known as the default mode network, which remains active in resting or passive states. The STORY task activates regions traditionally involved in the semantic comprehension of phrases along the superior temporal sulcus, including Wernicke's area, and in the inferior frontal region, including Broca's area⁶². This set of activation is often left lateralized as observed here (Fig. 1). The MATH task involves brain regions in the intraparietal, middle and inferior frontal regions, as well as the angular gyrus, when the math task is performed orally³⁹.

In the STORY-MATH contrast positive (respectively negative) beta weights reflect brain regions more activated by the STORY task than the MATH task, or deactivated by the MATH task (respectively more activated by the MATH task than the STORY task, or deactivated by the STORY task). Brain regions that are jointly activated by the STORY and MATH tasks have approximately null beta weights.

**Univariate analysis of additive genetic variance.** The variance components method, as implemented in the Sequential Oligogenic Linkage Analysis Routines (SOLAR) software package⁶³, was used for the heritability estimations of the phenotypes under analysis, such as the median activation in each areal¹⁷. The SOLAR algorithms use maximum variance decomposition methods derived from the strategy developed by Amos⁶⁴. The covariance matrix $\Omega$ for a pedigree of individuals is given by:

$$ \Omega = 2 \cdot \Phi \cdot \sigma_g^2 + I \cdot \sigma_e^2, \text{ where } \sigma_g^2 \text{ is the genetic variance due to the additive genetic factors, } \Phi \text{ is the kinship matrix representing the pair-wise kinship coefficients among all individuals, } \sigma_e^2 \text{ is the variance due to individual-specific environmental effects, and } I \text{ is the identity matrix.} $$

Narrow sense heritability is defined as the fraction of the phenotype variance $\sigma_p^2$ attributable to additive genetic factors: $h^2 = \sigma_g^2/\sigma_p^2$.

The significance of the heritability is tested by comparing the likelihood of the model in which $\sigma_g^2$ is constrained to zero with that of the model in which $\sigma_g^2$ is estimated. Before testing for the significance of heritability, phenotype values for each individual within the HCP cohort were adjusted for the following covariates: sex, age, age², age-sex interaction, age²-sex interaction, ethnicity (Hispanic or not) and education level. We used the number of years of education as a proxy for the education level to account for environmental differences in family socioeconomic status. This is a conservative approach because the number of years of education was shown to be associated not only with the family socioeconomic status (7%) but also with the general cognitive ability (3.5%)⁵⁵. Thus, it likely has shared environmental ground with the former and shared genetic origin with the latter. HCP data do not contain the information that would disentangle this issue. Following this last remark, one should note that the heritability estimates and shared genetic variances, described in the next section, were underestimated.

**Bivariate genetic analyses.** To assess the relationship between math dexterity/language comprehension and activation in brain areas, we computed the Pearson correlation between the median activation in each of the 180 areals of both hemispheres, and the behavioral scores.

Furthermore, we assessed the degree of shared genetic variance in the areals for which activation was significantly correlated with the behavioral scores; we performed a genetic correlation analysis using SOLAR, relying on the following model:

$$ \rho_p = \sqrt{h_a^2} \sqrt{h_b^2} \cdot \rho_g + \sqrt{1 - h_a^2} \sqrt{1 - h_b^2} \cdot \rho_e, \text{ where Pearson's phenotypic correlation } \rho_p \text{ is decomposed into } \rho_g \text{ and } \rho_e \cdot \rho_g \text{ is the proportion of variability due to shared genetic effects and } \rho_e \text{ that due to the environment, while } h_a^2 \text{ and } h_b^2 \text{ correspond to the previously defined narrow sense heritability for phenotypes a and b, respectively. In our case, one corresponds to the heritability of fMRI activation in one areal, while the second is the heritability of one of our behavioral scores.} $$

**Data Availability statement**

HCP (https://db.humanconnectome.org) is a publicly available dataset. Investigators need to apply to be granted access to restricted data.

**References**

1. Merriman, C. The intellectual resemblance of twins. *Psychol. Monogr.* **33**, i–57 (1924).

2. Falconer, D. S. *Introduction to quantitative genetics*. (Oliver and Boyd, 1960).

3. Peterson, R. L. & Pennington, B. F. Developmental dyslexia. *Lancet* **379**, 1997–2007 (2012).

4. Price, G. & Ansari, D. Dyscalculia: Characteristics, Causes, and Treatments. *Numeracy* **6** (2013).

5. Docherty, S. J. et al. A genome-wide association study identifies multiple loci associated with mathematics ability and disability. *Genes, Brain Behav.* **9**, 234–247 (2010).

6. Mascheretti, S. et al. Neurogenetics of developmental dyslexia: from genes to behavior through brain neuroimaging and cognitive and sensorial mechanisms. *Transl. Psychiatry* **7**, e987 (2017).

7. Männel, C. et al. Working-memory endophenotype and dyslexia-associated genetic variant predict dyslexia phenotype. *Cortex* **71**, 291–305 (2015).

8. Friederici, A. D. The Neural Basis of Language Development and Its Impairment. *Neuron* **52**, 941–952 (2006).

9. van der Lely, H. K. J. & Pinker, S. The biological basis of language: Insight from developmental grammatical impairments. *Trends Cogn. Sci.* **18**, 586–595 (2014).

10. Leonard, L. B. In *Language competence across populations: Toward a definition of specific language impairment*. (2003).

11. Koten, J. W. et al. Genetic Contribution to Variation in Cognitive Function: An fMRI Study in Twins. *Science* (80-). **323**, 1737–1740 (2009).
12. Blokland, G. A. et al. Heritability of working memory brain activation. *J Neurosci* **31**, 10882–10890 (2011).

13. Pinel, P. & Dehaene, S. Genetic and environmental contributions to brain activation during calculation. *Neuroimage* **81**, 306–316 (2013).

14. Pinel, P. et al. Genetic and environmental influences on the visual word form and fusiform face areas. *Cereb. Cortex* **25**, 2478–2493 (2015).

15. Barch, D. M. et al. Function in the human connectome: Task-fMRI and individual differences in behavior. *Neuroimage* **80**, 169–189 (2013).

16. Robinson, E. C. et al. MSM: A new flexible framework for multimodal surface matching. *Neuroimage* **100**, 414–426 (2014).

17. Glasser, M. F. et al. A multi-modal parcellation of human cerebral cortex. *Nature* **536**, 171–8 (2016).

18. Bishop, D. V. M. Cerebral Asymmetry and Language Development: Cause, Correlate, or Consequence? *Science* (80–). **340**, 1230531–1230531 (2013).

19. Price, G. R., Mazzocco, M. M. M. & Ansari, D. Why Mental Arithmetic Counts: Brain Activation during Single Digit Arithmetic Predicts High School Math Scores. *J. Neurosci.* **33**, 156–163 (2013).

20. Emerson, R. W. & Cantlon, J. F. Early math achievement and functional connectivity in the fronto-parietal network. *Dev. Cogn. Neurosci.* **2**, S139–S151 (2012).

21. Plomin, R. & Deary, I. J. Genetics and intelligence differences: five special findings. *Mol. Psychiatry* **20**, 98–108 (2015).

22. Davies, G. et al. Genome-wide association studies establish that human intelligence is highly heritable and polygenic. *Mol. Psychiatry* **16**, 996–1005 (2011).

23. Haworth, C. M. A., Dale, P. & Plomin, R. A. Twin Study into the Genetic and Environmental Influences on Academic Performance in Science in nine-year-old Boys and Girls. *Int. J. Sci. Educ.* **30**, 1003–1025 (2008).

24. Rimfeld, K., Kovas, Y., Dale, P. S. & Plomin, R. Pleiotropy across academic subjects at the end of compulsory education. *Sci. Rep.* **5**, 11713 (2015).

25. Sniekers, S. et al. Genome-wide association meta-analysis of 78,308 individuals identifies new loci and genes influencing human intelligence. *Nat. Genet.* **49**, 1107–1112 (2017).

26. Trampush, J. W. et al. GWAS meta-analysis reveals novel loci and genetic correlates for general cognitive function: a report from the COGENT consortium. *Mol. Psychiatry* **22**, 336–345 (2017).

27. Davies, G. et al. Genome-wide association study of cognitive functions and educational attainment in UK Biobank (N = 112 151). *Mol. Psychiatry* **21**, 1–10 (2016).

28. Amalric, M. & Dehaene, S. Origins of the brain networks for advanced mathematics in expert mathematicians. *Proc. Natl. Acad. Sci.* **113**, 4909–4917 (2016).

29. Rickard, T. et al. The calculating brain: an fMRI study. *Neuropsychologia* **38**, 325–335 (2000).

30. Pallier, C., Devauchelle, A.-D. & Dehaene, S. Cortical representation of the constituent structure of sentences. *Proc. Natl. Acad. Sci.* **108**, 2522–2527 (2011).

31. Skeide, M. A. & Friederici, A. D. The ontogeny of the cortical language network. *Nat. Rev. Neurosci.* **17**, 323–332 (2016).

32. Simmonds, D. J., Pekar, J. J. & Mostofsky, S. H. Meta-analysis of Go/No-go tasks demonstrating that fMRI activation associated with response inhibition is task-dependent. *Neuropsychologia* **46**, 224–232 (2008).

33. Erika-Florence, M., Leech, R. & Hampshire, A. A functional network perspective on response inhibition and attentional control. *Nat Commun* **5**, 4073 (2014).

34. Leisman, G., Machado, C., Melillo, R. & Mualem, R. Intentionality and “free-will” from a neurodevelopmental perspective. *Front. Integr. Neurosci.* **6**, 1–12 (2012).

35. Etkin, A., Egner, T. & Kalisch, R. Emotional processing in anterior cingulate and medial prefrontal cortex. *Trends Cogn. Sci.* **15**, 85–93 (2011).

36. Binder, J. R. et al. Mapping anterior temporal lobe language areas with fMRI: A multicenter normative study. *Neuroimage* **54**, 1465–1475 (2011).

37. Shakeshaft, N. G. et al. Strong Genetic Influence on a UK Nationwide Test of Educational Achievement at the End of Compulsory Education at Age 16. *PLoS One* **8**, e80341 (2013).

38. Krapohl, E. et al. The high heritability of educational achievement reflects many genetically influenced traits, not just intelligence. *Proc. Natl. Acad. Sci.* **111**, 15273–15278 (2014).

39. Dehaene, S., Piazza, M., Pinel, P. & Cohen, L. Three Parietal Circuits for Number Processing. *Cogn. Neuropsychol.* **20**, 487–506 (2003).

40. Vernes, S. C. & Fisher, S. E. In *Animal Models of Speech and Language Disorders* (ed. Helekar, S. A.) 13–40 (Springer New York, 2013).

41. Caspers, S. et al. The human inferior parietal cortex: Cytoarchitectonic parcellation and interindividual variability. *Neuroimage* **33**, 430–448 (2006).

42. Vogel, S. E., Goffin, C. & Ansari, D. Developmental specialization of the left parietal cortex for the semantic representation of Arabic numerals: An fMR-adaptation study. *Dev. Cogn. Neurosci.* **12**, 61–73 (2015).

43. Eicher, J. D. & Gruen, J. R. Imaging-genetics in dyslexia: Connecting risk genetic variants to brain neuroimaging and ultimately to reading impairments. *Mol. Genet. Metab.* **110**, 201–212 (2013).

44. Price, G. R., Holloway, I., Räsänen, P., Vesterinen, M. & Ansari, D. Impaired parietal magnitude processing in developmental dyscalculia. *Curr. Biol.* **17**, 1042–1043 (2007).

45. Ashkenazi, S., Rosenberg-Lee, M., Tenison, C. & Menon, V. Weak task-related modulation and stimulus representations during arithmetic problem solving in children with developmental dyscalculia. *Dev. Cogn. Neurosci.* **2**, S152–S166 (2012).

46. Skeide, M. A., Brauer, J. & Friederici, A. D. Brain Functional and Structural Predictors of Language Performance. *Cereb. Cortex* **26**, 2127–2139 (2016).

47. van Ettinger-Veenstra, H. M. et al. Right-hemispheric brain activation correlates to language performance. *Neuroimage* **49**, 3481–3488 (2010).

48. Booth, J. R. et al. Relation between brain activation and lexical performance. *Hum. Brain Mapp.* **19**, 155–169 (2003).

49. Bugden, S., Price, G. R., McLean, D. A. & Ansari, D. The role of the left intraparietal sulcus in the relationship between symbolic number processing and children's arithmetic competence. *Dev. Cogn. Neurosci.* **2**, 448–457 (2012).

50. Emerson, R. W. & Cantlon, J. F. Continuity and change in children's longitudinal neural responses to numbers. *Dev. Sci.* **18**, 314–326 (2015).

51. Diaz, M. T. & Hogstrom, L. J. The influence of context on hemispheric recruitment during metaphor processing. *J. Cogn. Neurosci.* **23**, 3586–3597 (2011).

52. Friederici, A. D. & Gierhan, S. M. The language network. *Curr. Opin. Neurobiol.* **23**, 250–254 (2013).

53. Brauer, J., Anwander, A. & Friederici, A. D. Neuroanatomical prerequisites for language functions in the maturing brain. *Cereb. Cortex* **21**, 459–466 (2011).

54. Skeide, M. A. et al. Genetic dyslexia risk variant is related to neural connectivity patterns underlying phonological awareness in children. *Neuroimage* **118**, 414–421 (2015).

55. Selzam, S. et al. Predicting educational achievement from DNA. *Mol. Psychiatry* **22**, 267–272 (2017).

56. Van Essen, D. C. et al. The Human Connectome Project: A data acquisition perspective. *Neuroimage* **62**, 2222–2231 (2012).

57. Glasser, M. F. et al. The minimal preprocessing pipelines for the Human Connectome Project. *Neuroimage* **80**, 105–124 (2013).
58. Jenkinson, M., Beckmann, C. F., Behrens, T. E. J., Woolrich, M. W. & Smith, S. M. Fsl. *Neuroimage* **62**, 782–790 (2012).

59. Fischl, B. FreeSurfer. *Neuroimage* **62**, 774–781 (2012).

60. Woolrich, M. W., Ripley, B. D., Brady, M. & Smith, S. M. Temporal Autocorrelation in Univariate Linear Modeling of FMRI Data. *Neuroimage* **14**, 1370–1386 (2001).

61. Glover, G. H. Deconvolution of Impulse Response in Event-Related BOLD fMRI. *Neuroimage* **9**, 416–429 (1999).

62. Binder, J. R., Desai, R. H., Graves, W. W. & Conant, L. L. Where is the semantic system? A critical review and meta-analysis of 120 functional neuroimaging studies. *Cereb. Cortex* **19**, 2767–2796 (2009).

63. Almasy, L. & Blangero, J. Multipoint Quantitative-Trait Linkage Analysis in General Pedigrees. *Am. J. Hum. Genet.* **62**, 1198–211 (1998).

64. Amos, C. I. Robust variance-components approach for assessing genetic linkage in pedigrees. *Am. J. Hum. Genet.* **54**, 535–543 (1994).

## Acknowledgements

This work was supported by the CEA DRF Impulsion program (SRAVS project 2017, P.I. Vincent Frouin). Data were provided by the Human Connectome Project, WU-Minn Consortium (Principal Investigators: David Van Essen and Kamil Ugurbil; 1U54MH091657) funded by the 16 NIH Institutes and Centers that support the NIH Blueprint for Neuroscience Research; and by the McDonnell Center for Systems Neuroscience at Washington University.

## Author Contributions

Y.L.G. and V.F. designed the analysis and wrote the manuscripts. Y.L.G. prepared the figures and tables. Y.L.G., V.F., M.A., P.P. and C.P. analyzed the results and reviewed the manuscripts.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-018-35665-0.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as 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. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018