Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE

User menu

  • Log out
  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log out
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
PreviousNext
Featured ArticleResearch Articles, Behavioral/Cognitive

Toward Robust Functional Neuroimaging Genetics of Cognition

Julia Uddén, Annika Hultén, Katarina Bendtz, Zachary Mineroff, Katerina S. Kucera, Arianna Vino, Evelina Fedorenko, Peter Hagoort and Simon E. Fisher
Journal of Neuroscience 30 October 2019, 39 (44) 8778-8787; DOI: https://doi.org/10.1523/JNEUROSCI.0888-19.2019
Julia Uddén
1Max Planck Institute for Psycholinguistics, Nijmegen, the Netherlands, 6525 XD,
2Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, the Netherlands, 6500 HE,
3Department of Linguistics,
4Department of Psychology, Stockholm University, Sweden, SE-106 91,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Annika Hultén
1Max Planck Institute for Psycholinguistics, Nijmegen, the Netherlands, 6525 XD,
2Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, the Netherlands, 6500 HE,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Katarina Bendtz
4Department of Psychology, Stockholm University, Sweden, SE-106 91,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zachary Mineroff
5Brain and Cognitive Sciences Department, Massachusetts Institute of Technology, Cambridge, Massachusetts, MA 02139-4307,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Katerina S. Kucera
1Max Planck Institute for Psycholinguistics, Nijmegen, the Netherlands, 6525 XD,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Arianna Vino
1Max Planck Institute for Psycholinguistics, Nijmegen, the Netherlands, 6525 XD,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Evelina Fedorenko
5Brain and Cognitive Sciences Department, Massachusetts Institute of Technology, Cambridge, Massachusetts, MA 02139-4307,
6McGovern Institute for Brain Research, Massachusetts Institute of Technology, Cambridge, Massachusetts, MA 02139, and
7Psychiatry Department, Massachusetts General Hospital, Charlestown, Massachusetts MA 02144
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Peter Hagoort
1Max Planck Institute for Psycholinguistics, Nijmegen, the Netherlands, 6525 XD,
2Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, the Netherlands, 6500 HE,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Simon E. Fisher
1Max Planck Institute for Psycholinguistics, Nijmegen, the Netherlands, 6525 XD,
2Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, the Netherlands, 6500 HE,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Simon E. Fisher
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

A commonly held assumption in cognitive neuroscience is that, because measures of human brain function are closer to underlying biology than distal indices of behavior/cognition, they hold more promise for uncovering genetic pathways. Supporting this view is an influential fMRI-based study of sentence reading/listening by Pinel et al. (2012), who reported that common DNA variants in specific candidate genes were associated with altered neural activation in language-related regions of healthy individuals that carried them. In particular, different single-nucleotide polymorphisms (SNPs) of FOXP2 correlated with variation in task-based activation in left inferior frontal and precentral gyri, whereas a SNP at the KIAA0319/TTRAP/THEM2 locus was associated with variable functional asymmetry of the superior temporal sulcus. Here, we directly test each claim using a closely matched neuroimaging genetics approach in independent cohorts comprising 427 participants, four times larger than the original study of 94 participants. Despite demonstrating power to detect associations with substantially smaller effect sizes than those of the original report, we do not replicate any of the reported associations. Moreover, formal Bayesian analyses reveal substantial to strong evidence in support of the null hypothesis (no effect). We highlight key aspects of the original investigation, common to functional neuroimaging genetics studies, which could have yielded elevated false-positive rates. Genetic accounts of individual differences in cognitive functional neuroimaging are likely to be as complex as behavioral/cognitive tests, involving many common genetic variants, each of tiny effect. Reliable identification of true biological signals requires large sample sizes, power calculations, and validation in independent cohorts with equivalent paradigms.

SIGNIFICANCE STATEMENT A pervasive idea in neuroscience is that neuroimaging-based measures of brain function, being closer to underlying neurobiology, are more amenable for uncovering links to genetics. This is a core assumption of prominent studies that associate common DNA variants with altered activations in task-based fMRI, despite using samples (10–100 people) that lack power for detecting the tiny effect sizes typical of genetically complex traits. Here, we test central findings from one of the most influential prior studies. Using matching paradigms and substantially larger samples, coupled to power calculations and formal Bayesian statistics, our data strongly refute the original findings. We demonstrate that neuroimaging genetics with task-based fMRI should be subject to the same rigorous standards as studies of other complex traits.

  • fMRI
  • FOXP2
  • individual differences
  • KIAA0319/TTRAP/THEM2
  • language
  • neuroimaging genetics

Introduction

Advances in genomics are helping to identify genes contributing to key aspects of human cognition, including language (Deriziotis and Fisher, 2017). Insights come partly from rare high-impact mutations disrupting brain development (Fisher and Scharff, 2009). Nonetheless, even for highly heritable traits, the genetic architecture of interindividual variability in cognition and behavior mainly involves effects of many common variants (polymorphisms) at different genomic loci. Typically, one common polymorphism has, at best, tiny effects on the relevant behavioral/cognitive phenotypes, so large cohorts of participants are needed to robustly detect genetic associations.

Neuroimaging enables noninvasive investigation of brain structure and function in living humans. A widely held view is that neuroimaging-derived measures, being closer to the underlying biology, are more amenable for tracing links to genetic variation than behavioral/cognitive assessments (Bigos et al., 2016; but see Flint and Munafò, 2007 for an opposing view). Twin studies using MRI-derived data (Thompson et al., 2001) suggest that a substantive proportion of interindividual variability in brain structure is explained by genetic variation in aggregate. However, systematic genome-wide studies indicate that genetic accounts of neuroanatomical variability are as complex as distal measures of behavior, such that effect sizes of individual polymorphisms remain small. For example, in screening of subcortical volumes in >30,000 individuals, the strongest association signal accounted for <0.52% of variance in putamen volume (Hibar et al., 2015).

Perhaps task-based functional MRI (fMRI), indexing neural activation related to particular cognitive processes, holds greater promise for bridging to genetics? Twin studies again support significant genetic contributions, although heritability estimates vary by task and/or brain region (Koten et al., 2009; Blokland et al., 2011). Whereas large genetic screens are feasible for structural MRI, scaling up task-based fMRI is difficult, with most studies limited to tens of individuals. Lacking power for genome-wide scans, task-based fMRI genetics efforts target candidate genes, yielding claims of positive associations with different polymorphisms for various paradigms (Grabitz et al., 2018). However, small cohorts in studies of complex traits have not only reduced power but also elevated false-positive rates (Button et al., 2013). Moreover, functional neuroimaging genetics involves multidimensional datasets, with multiple flexible parameters when processing/analyzing primary data, dramatically increasing the degrees of freedom (Simmons et al., 2011). This raises risks of p-hacking and HARKing (hypothesizing after results are known) (Grabitz et al., 2018). Few positive genetic associations with task-based fMRI have been formally replicated using matching designs; there are no direct replications in the language processing literature.

Here, we address this gap by focusing on the most highly cited neuroimaging genetics study of language function (Pinel et al., 2012) (179 Google Scholar citations at time of writing, corresponding to ∼26 citations per year). Pinel et al. (2012) correlated single-nucleotide polymorphisms (SNPs) with task-based activation in language-related brain areas in 94 healthy participants, targeting two candidate loci from genetic investigations of language-related disorders: FOXP2 and KIAA0319/TTRAP/THEM2. Rare FOXP2 mutations cause a monogenic disorder, involving speech apraxia, expressive/receptive language impairments (Lai et al., 2001; Fisher and Scharff, 2009), and distributed alterations in brain structure and function (Watkins et al., 2002; Liégeois et al., 2003). The KIAA0319/TTRAP/THEM2 locus was selected from genetic studies of dyslexia; clusters of common SNPs in this region have been associated with variation in reading skills (reviewed by Carrion-Castillo et al., 2013). On analyzing 39 SNPs from FOXP2 and KIAA0319/TTRAP/THEM2, Pinel et al. (2012) found that three were associated with altered activation during covert sentence reading, each with a distinct pattern (Fig. 1). rs6980093 and rs7784315 (FOXP2 SNPs) were correlated with variable activation in left inferior frontal and precentral gyri respectively. rs17243157 (KIAA0319/TTRAP/THEM2) was associated with variation in functional asymmetry of the superior temporal sulcus.

The present study attempted a hypothesis-driven replication of these findings. With 427 healthy participants, a sample four times that of Pinel et al. (2012), we believe this to be the largest fMRI genetics investigation of language-related activation to date. We demonstrate nonreplication of all three findings of the prior report.

Materials and Methods

Cohorts

The 427 participants in the present study came from two independently collected cohorts in which language task-based fMRI testing was coupled with DNA collection: referred to here as MOUS (Mother of Unification Studies) and EvLab. The MOUS cohort comprised 217 participants, 107 of whom performed a sentence reading (visual language) task and 110 performed a sentence listening (auditory language) task. The EvLab cohort comprised 210 participants, all of whom performed a visual language task. Table 1 summarizes the fMRI design and comparison to the Pinel et al. (2012) study, detailed further below.

View this table:
  • View inline
  • View popup
Table 1.

Summary of the key design features for the original study and the two cohorts investigated in the present report

Participants

MOUS cohort.

A total of 242 native Dutch-speaking participants volunteered to participate in MOUS, an in-depth multimethod neuroimaging study of language processing, performed at the Donders Centre for Cognitive Neuroimaging (DCCN) in Nijmegen, the Netherlands. All participants completed an fMRI and a MEG session. Of those, 25 participants were excluded for a range of reasons, including: technical problems during data collection; poor data quality due to excessive blinking during initial MEG measurement, leading to exclusion from subsequent fMRI sessions; failure to complete the study; noncompliance with task instructions assessed via a behavioral threshold (see task description below). Approximately half of the remaining 217 participants read sentences and word lists, presented word-by-word (visual group: 107 participants; 51% male; mean age of 22.3 years, age range 18 to 33 years). The other half listened to auditory versions of the same materials (auditory group: 110 participants; 46% male; mean age of 22.3 years, age range 18 to 30 years). We will refer to the “cohorts” (to distinguish the MOUS and the EvLab sample) as well as the “subsamples” (to distinguish the visual MOUS sample, the auditory MOUS sample and the EvLab cohort). Key design features of the cohorts are noted in Table 1. All participants were right-handed as assessed by the Bever handedness questionnaire (building on Oldfield, 1971), had normal or corrected-to-normal vision, and reported no history of neurological, developmental or language deficits. Participants were instructed to not use medication, alcohol or drugs on the day of measurement. The study was approved by the local ethics committee (CMO, the local “Committee on Research Involving Human Subjects” in the Arnhem-Nijmegen region) and followed the guidelines of the Helsinki declaration.

EvLab cohort.

The EvLab cohort was collected by Evelina Fedorenko's laboratory at the Massachusetts Institute of Technology (MIT) in Boston. This cohort consisted of 210 native English-speaking participants (36% male; mean age of 26.9 years, age range 20 to 59 years) who were students at MIT and members of the larger Boston community. They participated for payment in different neuroimaging experiments that all included a language localizer task (Fedorenko et al., 2010). As in the visual subset of the MOUS cohort, the EvLab participants read sentences presented word-by-word, as well as one or more control conditions, like lists of words and/or nonwords (Table 1). All participants were right-handed as assessed by the Oldfield (1971) handedness questionnaire or by self-report. All participants gave informed consent in accordance with the requirements of the MIT's Committee on the Use of Humans as Experimental Subjects (COUHES).

Language tasks

Both the MOUS and the EvLab cohort used a blocked design, in contrast to the Pinel et al. (2012) study, which relied on an event-related design (Table 1). In a blocked design, several stimuli from the same condition appear in a row (forming a block), and the response is estimated to the entire block. Given the additive nature of the BOLD signal, blocked designs have greater sensitivity (Birn et al., 2002) and are generally recommended over event-related designs except in cases where local predictability of experimental events could be problematic or where it is critical to examine responses to individual experimental events. Given that the primary goal of both the current study and of Pinel et al. (2012) was to reliably identify language-responsive cortical areas, so that properties of these areas could be related to genetic variation, a blocked design provides additional sensitivity.

Specifics of the paradigms have been reported in detail previously (Lam et al., 2016; Mahowald and Fedorenko, 2016; Schoffelen et al., 2017; Hultén et al., 2019; Uddén et al., 2019). Further specifics of the MOUS cohort are present in a peer reviewed data publication (Schoffelen et al., 2019). To check for compliance, participants in the MOUS cohort were presented with yes/no comprehension questions on 10% of the sentence trials. In the EvLab cohort, participants were asked to press a button at the end of each sentence, or to respond to a memory probe task (Fedorenko et al., 2010; Scott et al., 2017; it has been previously established that the nature of the task does not affect the resulting activation maps).

MRI acquisition

For both cohorts, structural and functional data were acquired with a SIEMENS Trio 3 T scanner (the MOUS cohort at the DCCN of the Donders Institute in Nijmegen and the EvLab cohort at the Athinoula A. Martinos Imaging Center at the McGovern Institute for Brain Research at MIT). The functional scans used a TR of 2000 ms in both cases. Further specifics of the acquisition paradigms have been reported in detail previously (Mahowald and Fedorenko, 2016; Uddén et al., 2019).

Preprocessing and first-level analysis

Preprocessing and first-level analysis procedures were matched as closely as possible to those used in Pinel et al. (2012), to maximize the comparability of the results. We used the statistical parametric mapping software (SPM8; Wellcome Trust Centre for Neuroimaging, London; www.fil.ion.ucl.ac.uk/spm), including realignment to correct for individual subject head movement. Thus, one of the most critical sources of noise affecting the BOLD signal (i.e., subject motion) was controlled for in a similar manner across cohorts. A second important noise source: bodily physiology of the participant (e.g., heart rate), was not controlled for in any of the cohorts (Dubois and Adolphs, 2016). A correction for differences in slice acquisition time was performed for the MOUS cohort (ascending slice acquisition) but not for the EvLab cohort (interleaved slice acquisition). Structural images were spatially normalized to the structural image (T1) template provided by SPM8, using affine regularization. In this step, to match the original Pinel et al. (2012) study, the voxel size was resampled to 3.0 × 3.0 × 3.0 mm3. The transformation matrices generated by the normalization algorithm were applied to the corresponding functional EPI-BOLD volumes, after coregistration. All structural and functional images were smoothed with an isotropic 3D spatial Gaussian kernel (FWHM = 5 mm) matching Pinel et al. (2012). For the single-subject fixed effect analyses, we modeled the six realignment parameters from the movement correction. As in the original study, we used the canonical HRF and its temporal derivative. For the MOUS cohort, comprehension questions were modeled separately. As an example, in a block where a question appeared after sentence #1, we modeled two events with sentences: sentence 1 and sentences 2–5. As in Pinel et al. (2012), we computed individual asymmetry maps from the activation maps by subtracting the right from the left hemisphere, voxel by voxel, after normalizing the right hemisphere onto the left hemisphere.

Second-level region of interest (ROI) analysis

Because we sought to directly replicate the key findings reported by Pinel et al. (2012), the ROIs for this study were based on the locations of the peak associations reported for each of the significant SNPs in that study (Fig. 1). Specifically, we used a 10 mm spherical ROI around each of the three peak coordinates described in Pinel et al.: in the left inferior frontal gyrus (rs6980093), left precentral gyrus (rs7784315), and superior temporal gyrus/sulcus (rs17243157). Our analytical approach matched that used for each SNP in the original study. Specifically, for rs6980093 and rs7784315, we assessed left hemisphere activation for the sentence condition (relative to the fixation baseline, as in Pinel et al., 2012), extracted the average (across voxels) BOLD response in each ROI of each participant using MarsBaR (http://marsbar.sourceforge.net/) and analyzed the resulting averages in SPSS. For rs17243157, we analyzed asymmetry maps for the sentences > fixation contrast, averaged the difference values (left minus right, as described above; extracted using the spm_read_vols and Imcalc functions of SPM) across voxels in the ROI of each participant, and analyzed the resulting averages in SPSS. Participants with a deviation larger than four times the interquartile range from the mean within a subsample (and ROI), were considered outliers. Thus, if a participant was considered an outlier with respects to one (or several) ROIs, they were removed from all tests (this concerned two MOUS auditory subsample participants).

To test for the presence of reliable above-baseline activity in these three ROIs at the group level, regardless of the potential modulation by genotype group, we performed one-sample t tests (against 0) at the second level, for each subsample. If there were significant voxels at p < 0.05 (using FWE correction) within the spherical ROI, we report those voxels as indicators of group-level activity within a ROI.

Power analysis

We used G*Power (http://www.gpower.hhu.de/en.html) to calculate the sample size we would need to have 80% power to detect the effect sizes (Cohen's f2 for F tests; Cohen's d for t tests) reported by Pinel et al. (2012). t tests (see “Analysis structure” section for further information) were taken to be one-tailed, since we aimed to replicate a finding in the specific direction. In addition, again assuming 80% power, we estimated the minimal effect size that could be detected in the present investigation, both for the total combined sample and each separate cohort. Power depends not only on the total number of participants but also on the relative sizes of the different genotype groups, in this case, determined by the genotype frequencies for each SNP. We used the established allele frequencies (Table 2) from the European subset of the 1000 Genomes data recorded in the dbSNP database (https://www.ncbi.nlm.nih.gov/projects/SNP/) and assumed Hardy–Weinberg equilibrium.

View this table:
  • View inline
  • View popup
Table 2.

Allele and genotype frequencies used in the power analysis

Genotyping

The participants' saliva was collected with Oragene DNA collection OG500 kits (DNA Genotek) and DNA was extracted according to manufacturer's instructions. This investigation was a hypothesis-driven study in which we selected a priori those SNPs from Pinel et al. (2012) that best represented the major conclusions of that earlier report (rs6980093, rs7784315, and rs17243157, as shown in Fig. 2 of Pinel et al., 2012, and in Fig. 1 of the present manuscript). Because we were focusing on three specific SNPs, we did not perform genome-wide genotyping of hundreds of thousands of markers on DNA-chips for this study. Rather, we directly genotyped only the three SNPs of interest, using a specialized technique that has high sensitivity and specificity, known as KASP (Kompetitive Allele-Specific PCR; developed by LGC Genomics). Targeted genotyping of each SNP of interest was performed using customised KASP assays and a Bio-Rad CFX96 real-time PCR thermocycler. Each 10 μl reaction contained 0.14 μl of 72× KASP primer mix, 5 μl of 2× KASP master mix with standard ROX concentration, and 2 μl of genomic DNA (diluted at 20 ng/μl). Thermocycling was performed as follows: after initial 15 min denaturation at 94°C, 10 cycles were run with 20 s denaturation at 94°C followed by 45 s of annealing/elongation starting at 61°C and decreasing by 0.6°C per cycle, followed by 30 cycles with 10 s denaturation at 94°C followed by 45 s of annealing/elongation at 55°C. When necessary, five further cycles were performed using the same parameters as the final 30 cycles from the initial run. Genotypes were called using the CFX96 Manager software (Bio-Rad). For each SNP assay, every sample was genotyped in duplicate, and we also included negative controls and positive controls (DNA samples for which genotypes were already known) for each allelic combination, allowing us to verify that genotype calls were reproducible and accurate. In the few cases where samples yielded genotype calls that were unclear with KASP, we used direct Sanger sequencing to resolve ambiguities. Moreover, we performed additional validations by Sanger sequencing of randomly selected samples of each genotype. Through these efforts, we successfully determined genotypes for all three SNP markers of interest in all 427 participants.

Statistical tests

The effects of genotype groups on the BOLD response during sentence processing were analyzed within each cohort separately and in a combined analysis across cohorts, with F tests, as follows. For rs6980093, the analysis was run with an F test comparing all three allelic groups: AA homozygotes, AG heterozygotes, GG homozygotes. For rs7784315 and rs17243157, the frequency of the minor allele was only ∼10% (consistent with population data, as shown in Table 2), such that the homozygous minor allele groups had sample sizes that would be much too small for meaningful analyses (see “Results” for further information). Thus, for each of those SNPs, we used an F test comparing two independent means, corresponding to two allelic groups: the major allele homozygotes and the heterozygotes. In these models, we included main effects and two-way interactions with genotype group and the following two potential confounding factors. Sex was included as a potential confounding factor for both cohorts. As the EvLab cohort contained slight variations in the experimental paradigm used across participants, we also included experiment version (“experiment” for short) as a possible confounding factor in addition to sex, for this cohort. Experiments with ≤5 participants were excluded from the EvLab cohort in this analysis to ensure robustness of the statistical tests. Due to the inclusion of potential confounding factors in the models, any results reported on the main effect of genotype, which is our focus, cannot be explained by confounding factors (regardless of their potential significance). In reporting our results, we provide information on the direction of effects observed in the earlier Pinel et al. (2012) study. Note, however, that directionality was not constrained in any test in this analysis.

We performed an additional set of analyses in JASP (https://jasp-stats.org/, JASP Team, 2017, version 0.8.4) to assess the Bayesian evidence for or against the hypothesis presented in the Pinel et al. (2012) report. The effects of genotype group on the BOLD response during sentence processing were analyzed either with a Bayesian F test or t test, as follows. For rs6980093, the analysis was run with a Bayesian F test with three genotype groups, corresponding to the three allelic groups (AA, AG, GG). The H0 corresponds to no effect of genotype group and the H1 to any effect of genotype group. We report the Bayes Factor (BF01) comparing these two hypotheses. Following standard guidelines for interpreting Bayes Factors, BF01 > 1 is considered support for H0 and BF01 < 1 is considered support for H1, although 0.3 < BF01 < 3 is regarded as weak evidence. For rs7784315 and rs17243157, each involving two allelic groups as explained above, we ran Bayesian t tests comparing two independent means. Based on the effects reported by Pinel et al. (2012) the H1 was TC>TT for rs7784315, and CC>TC for rs17243157. In these cases, the H0 corresponds to no genotypic effect, or an effect in the opposite direction. We used the default Cauchy prior distributions in JASP. Bayes factors are reported and labeled as inconclusive/no evidence, weak, substantial, or strong according to the guidelines set out by Jeffreys (1961) and JASP. Note that the Bayesian analyses, as currently implemented in JASP, did not allow for inclusion of potential confounding factors in these models. Bayes factors can only be estimated for complete models—that is, all main effects and interactions must be considered together as the H1. Thus, we used the initial SPSS analysis to first exclude effects of possible confounding factors and then used a model with allelic group as the only factor for the subsequent Bayesian analysis.

Analysis structure

The analytic approach we have described thus had three parts: a power analysis, a classical statistical analysis, and a Bayesian analysis. All analyses were performed with comparable sets of tests that were matched as closely as possible to each other. In this paragraph, we comment on the correspondence of tests across these parts. Because the confounding factors (sex and experiment version; see “Statistical tests” section) were not reported in the original study, the confounding factors were not included in the power analysis. The power analysis was performed assuming an F test for rs6980093 and t tests for rs7784315 and rs17243157, as for rs6980093 there are three allele groups and for the other two SNPs there are two. The same division between F- and t tests across SNPs would also have been relevant for the classical statistical analysis, unless confounding factors had been included. However, because we sought to control for confounding factors in the classical analysis, F tests were in that case used for all SNPs. In the Bayesian statistical analysis, we did not include possible confounding factors, mainly because Bayesian model selection procedure compares the fit of total models and does not allow for single factors to be evaluated. In addition, because possible confounding factors had been included in the classical analysis, we were able to verify that they did not confound the results in this first step and then leave them out of the Bayesian analysis.

Results

The effect sizes of the findings reported by Pinel et al. (2012) were relatively large (reporting Cohen's f2 for F tests; Cohen's d for t tests): Cohen's f2 = 0.29 for rs6980093, Cohen's d = 1.04 for rs7784315, and Cohen's d = 1.00 for rs17243157. Our power calculations indicated that a cohort of 39 (rs6980093), 44 (rs7784315), or 40 (rs17243157) participants would give 80% power to detect effects of this magnitude, assuming an α level of 0.05. In our current study, we used cohorts that are substantially larger than this. Moreover, to account for the potential inflation of effect sizes in the Pinel et al. (2012) study (e.g., due to the winner's curse, Button et al., 2013), we also estimated the minimum effect size that we could detect in our cohorts, both for the separate cohorts and for the entire sample. For rs6980093, the minimum effect size was a Cohen's f2 of 0.07 (separate cohorts) or 0.03 (entire sample). For rs7784315 and rs17243157, the minimum effect sizes (Cohen's d) for separate cohorts were 0.48 and 0.43, respectively, and for the entire sample, 0.32 and 0.29, respectively.

We used direct targeted genotyping to successfully determine the allelic status of rs6980093, rs7784315, and rs17243157 in all 427 participants. The minor allele frequencies were consistent with population data and there was no deviation from Hardy–Weinberg equilibrium for these SNPs in our cohorts or any of the subsamples (Table 3). For rs6980093, there were sufficient numbers of participants in each of the three allelic groups. For rs7784315 and rs17243157, the low minor allele frequency meant that homozygotes for this allele were rare (consistent with prior studies); between 0 and 4 participants in each subsample. Hence, and in line with the analyses performed by Pinel et al. (2012), the minor allele homozygotes were excluded from the imaging genetics analyses.

View this table:
  • View inline
  • View popup
Table 3.

Genotyping results

For all brain ROIs in all subsamples, the presence of significant voxels in the ROI established that there was activity at the group level for the sentence reading/listening task, as expected from prior studies. For all three cohorts, group level activation maps (Pinel, MOUS) or probabilistic activation overlap maps (Evlab), for the sentence versus low level baseline contrast, have been published previously (Pinel et al., 2007; Uddén et al., 2019, Mahowald and Fedorenko, 2016). Visual inspection of these maps showed reasonably high consistency across subsamples, in particular in the areas of interest. Peak voxels for the ROIs had the following locations in MNI space: MOUS Visual, [−57 17 10], left inferior frontal gyrus, T = 7.5, p < 0.001, [−57 −16 40], left precentral gyrus, T = 6.2, p < 0.001, [−57 −37 4], left superior temporal gyrus/sulcus, T = 15.9, p < 0.001; MOUS Auditory, [−57 −1 4], left inferior frontal gyrus, T = 6.49, p < 0.001, [−45 −4 40], left precentral gyrus, T = 4.2, p = 0.002, [51 −40 1], left superior temporal gyrus/sulcus, T = 11.8, p < 0.001; EvLab, [−54 14 10], left inferior frontal gyrus, T = 12.5, p < 0.001, [−48 −4 37], left precentral gyrus, T = 13.4, p < 0.001, [−54 −40 4], left superior temporal gyrus/sulcus, T = 23.6, p < 0.001. Note that small differences in the precise locations of peak group-level voxels are to be expected across cohorts given the well established interindividual anatomical and functional variability (Tomaiuolo et al., 1999; Fedorenko et al., 2010; Frost and Goebel, 2012). Pinel et al. (2012) reported significant effects of the FOXP2 SNP rs6980093 on left inferior frontal gyrus activity for both auditory and visual modalities, with activity increasing with the number of A alleles (i.e., AA>AG>GG) (Fig. 1). We did not replicate this finding in any of the subsamples (Fig. 2, Table 4). For rs7784315, another FOXP2 SNP in a different part of the gene, Pinel et al. (2012) described effects in the left precentral gyrus, with TC participants showing significantly higher activity than TT participants, but the association was only seen for the visual modality (Fig. 1). We did not replicate this finding in either of the two visual subsamples (Fig. 2, Table 4). As an additional analysis that must be considered exploratory, since no effect was observed in the Pinel et al. (2012) study, we tested rs7784315 in the MOUS auditory subsample. This auditory analysis generated a nominally significant uncorrected P-value (TC>TT group; p = 0.02), but the finding did not survive the appropriate Bonferroni correction for multiple comparisons. Finally, for rs17243157, a SNP in the KIAA0319/TTRAP/THEM2 locus, Pinel et al. (2012) reported effects on asymmetry of activation in the superior temporal sulcus, with CC participants showing higher asymmetry than TC participants, for both auditory and visual modalities (Fig. 1). Again, we did not replicate these associations in any of the subsamples (Fig. 2, Table 4).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Hypotheses in the present study were directly based on prior findings of Pinel et al. (2012), displayed here (adapted from Fig. 2 of the original paper). FOXP2 spans ∼610 kilobases of chromosome 7q31.1. First, Pinel et al. (2012) associated rs6980093 and rs7784315 allelic status with reading-related activation of regions within the left inferior frontal gyrus and left precentral gyrus, respectively, as highlighted in the sagittal anatomical slices of the left panels. The right panels show the reported percentage changes in BOLD signal at cortical peaks of the left hemisphere (L.H) and the symmetrical position in the right hemisphere (R.H), plotted separately for the different allelic groups for the reading task. Following up on this association, Pinel et al. (2012) analyzed the corresponding speech-listening data, with results displayed to the very right. KIAA0319, TTRAP, and THEM2 are neighboring genes spanning ∼160 kilobases of chromosome 6p22.3. Pinel et al. (2012) reported that rs17243157 allelic status was associated with altered functional asymmetry indices in the superior temporal gyrus/sulcus during the reading task, for the region highlighted in the sagittal anatomical slice of the left panel. The right panel shows reported percentage changes in peak BOLD signal of each hemisphere for the allelic groups, for the reading and speech-listening tasks. ROI locations of the present study are indicated with red outlines on each brain slice.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Single-subject contrast estimates by genotype group. The data points are horizontally jittered for visibility. The median as well as the first and third quartiles are shown as horizontal bars. For FOXP2 rs6980093 in the original study, there was a negative relationship between the number of G alleles and the BOLD response (AA>AG>GG) to reading or listening to sentences in left inferior frontal gyrus. This finding was not replicated in either modality in the MOUS cohort, nor in the EvLab cohort. For FOXP2 rs7784315 in the original study, the TC group had higher BOLD response in the left precentral gyrus when reading sentences, which is not replicated in any subsample in the current study. The corresponding results for the MOUS auditory subsample are reported for completeness, although we stress that there was no significant effect in the original study. For rs1724315 at the KIAA0319/TTRAP/THEM2 locus, the original results showed greater functional asymmetry of the posterior STS for the CC compared with the TC group in the visual subsample, which is also not replicated in any subsample in the current study. The corresponding results for the MOUS auditory subsample are reported for completeness.

View this table:
  • View inline
  • View popup
Table 4.

Nonreplication of each of the main findings from Pinel et al. (2012) in multiple independent samples and in a combined analysis using F tests to assess genotype associations with ROI activations

We excluded the possibility that the sex of the participants might confound our results. In the MOUS cohort, sex*genotype interactions as well as sex as a main effect were nonsignificant for all SNPs and for both visual (p ≥ 0.14 for rs6980093, p ≥ 0.35 for rs7784315, p ≥ 0.27 for rs17243157) and auditory subsamples (p ≥ 0.54 for rs6980093, p ≥ 0.30 for rs7784315, p ≥ 0.29 for rs17243157). In the EvLab cohort, for the analysis of rs6980093, sex was a marginally significant factor (F(1,201) = 3.96, p = 0.05), and all other factors, including sex*genotype interactions, were nonsignificant (p ≥ 0.48). For the EvLab analysis of rs7784315, no confounding factors were significant (p ≥ 0.18), while for the analysis of rs17243157 “experiment” was significant (F(1,197) = 6.63, p = 0.01) and all other factors were nonsignificant (p ≥ 0.09).

In a joint analysis of the MOUS and the EvLab cohorts, combining all the available data from >400 people, none of the findings from Pinel et al. (2012) were replicated (Table 4). In these analyses, we included sex and “experiment” (see “Statistical tests” section for definitions) in the models, with MOUS auditory and MOUS visual paradigms considered as two different experiment versions. Note, however, that “experiment” in the joint analysis reflects a combination of factors (such as the exact set of participants, slight differences in the materials/procedure, any subtle variations in the scanning environment that are time-dependent, etc.), rather than an effect of sensory input modality per se. For combined analysis of rs6980093 the effect of “experiment” was significant (F(3,416) = 5.61, p = 10−3), as well as the effect of sex (F(1,416) = 5.01, p = 0.03), but the interactions sex*genotype and “experiment”*genotype were all nonsignificant (p ≥ 0.60). For combined analysis of rs7784315 the effect of “experiment” was significant (F(2,304) = 5.78, p = 3 × 10−3), but the effect of sex as well as the interactions sex*genotype and “experiment”* genotype were all nonsignificant (p ≥ 0.58). For combined analysis of rs17243157, the effect of “experiment” was significant (F(3,411) = 3.50, p = 0.02), but the effect of sex as well as the interactions sex*genotype and “experiment”*genotype were all nonsignificant (p ≥ 0.45). Thus, we excluded the possibility that sex and experiment version were confounding factors that would contribute to assessment of genotypic effects.

Finally, we used Bayesian analyses to formally evaluate the robustness of our findings (Table 5). Almost all of the nonreplication results reported above were robustly informative in the Bayesian sense, providing overall substantial to strong evidence for nonreplication (Table 5). The exceptions were in the MOUS visual subsample, with Bayes Factors of 2.96 (i.e., just below the threshold value of three) for rs6980093 and 2.42 for rs17243157. As noted above, rs7784315 showed no effect on left precentral gyrus activation in the auditory modality in the original Pinel et al. (2012) study; their positive (TC>TT) findings were restricted to the visual task. In an exploratory Bayesian analysis of the MOUS auditory subsample we found no evidence for a TC>TT effect of rs7784315 on activation of this region (BF01 = 0.30).

View this table:
  • View inline
  • View popup
Table 5.

Bayesian analysis provides evidence of nonreplication

Discussion

Here, we attempted to replicate an influential fMRI-genetics investigation of language processing; a field with no direct replications in the literature. Pinel et al. (2012) reported association between rs6980093 and rs7784315 (FOXP2 SNPs) and strength of the BOLD response during sentence reading/listening, in the left inferior frontal gyrus and precentral gyrus, respectively. They also reported association between rs17243157 (in the KIAA0319/TTRAP/THEM2 locus) and asymmetry of the BOLD response during sentence reading/listening, in the posterior superior temporal sulcus. Despite including a substantially larger sample, and a more powerful blocked design (Birn et al., 2002), we did not replicate the Pinel et al. (2012) associations. Formal Bayesian analyses yielded substantial-to-strong evidence supporting the null hypotheses for each SNP. The lack of association was consistent across two independent cohorts (MOUS and EvLab), each itself larger than the original cohort. Power calculations indicate that combined investigations of MOUS and EvLab had 80% power to detect effect sizes of considerably lower magnitude than those observed in Pinel et al. (2012). Nevertheless, in combined analyses (>400 participants), association tests were all nonsignificant, indicating unambiguous lack of replication.

Our results contribute to a growing appreciation of the small magnitude of effects in neuroimaging genetics of common DNA variants, demonstrating that this limitation extends to task-based fMRI of higher cognition. It is worth considering factors that might have contributed to observations of positive signals in Pinel et al. (2012).

First, the earlier study began with exploratory screening of 39 SNPs distributed across two chromosomal regions (FOXP2 on 7q31, KIAA0319/TTRAP/THEM2 on 6p22.3), testing every marker for association with three different brain ROIs (based on fMRI-studies of different language-related disorders), each investigated using multiple approaches to phenotype definition, including hemisphere-specific activations and fMRI-based asymmetry indices. With multiple-testing adjustment for the effective number of independent SNPs, the significant p-values from this screen were Pcorr(Meff) = 0.0286 for rs6980093 (left activation, inferior frontal gyrus), Pcorr(Meff) = 0.0182 for rs7784315 (left activation, precentral gyrus), and Pcorr(Meff) = 0.0234 for rs17243157 (functional asymmetry, posterior superior temporal sulcus) (Fig. 1 in Pinel et al. (2012)). These findings cluster around p = 0.05; none would survive further multiple-testing adjustment to account for investigating several brain ROIs, and exploring multiple phenotypic definitions.

Second, to validate effects suggested by their exploratory screen, instead of testing top SNPs and brain ROIs in independently phenotyped/genotyped samples, Pinel et al. (2012) performed further SPM analyses of the same marker data, focusing on the cortical locations that had shown peak association with those SNPs. Since this involved reanalyzing the same input data used in the initial screen, testing hypotheses that had emerged from analyzing those data, it is a foregone conclusion that significant SNP effects were observed for the designated cortical locations (Kriegeskorte et al., 2009).

Another often neglected issue for designing and interpreting neuroimaging genetics studies concerns the question of whether/how SNPs impact gene function. While there is considerable interindividual variation in DNA sequence in humans, spread throughout the genome, the majority of common SNPs have no effects on expression, regulation or activity of genes and their encoded proteins. Since much of the interindividual variation observed at a gene locus is functionally neutral, most SNPs used in association screening have no mechanistic relevance. Experimental efforts are underway in genomics, using cellular models and tissue samples, to pinpoint functional DNA variants against the huge background of nonrelevant noise, and to understand how such functional variants alter expression, regulation and activity of gene products (Gasperini et al., 2016; GTEx Consortium, 2017). For association studies in cohorts with low power, when multiple-testing becomes a limiting factor, evidence from functional experiments should ideally constrain SNP choice, before association testing. Pinel et al. (2012) only considered biological relevance as a post hoc discussion point, after observing which SNPs showed putative genotype-phenotype correlations.

The FOXP2 locus is large (>610 kilobases) with many common SNPs in its introns, noncoding parts that are spliced out of mRNA transcripts before translation into protein (Hoogman et al., 2014). Unlike our mechanistic knowledge concerning rare disruptions (Sollis et al., 2017), we have limited understanding of which common SNPs in FOXP2 have biological impact on its functions (Becker et al., 2018). Pinel et al. (2012) acknowledged the absence of experimental evidence that rs6980093 and/or rs7784315 affect FOXP2 functions; the authors speculated based on their approximate locations within the gene locus, but such arguments could be made post hoc for many different SNPs. Even now, potential molecular mechanisms, along with an explanation of why the two SNPs should affect distinct brain ROIs, remain elusive. For rs17243157 in the KIAA0319/TTRAP/THEM2 locus, Pinel et al. (2012) noted that this marker tends to cosegregate with another SNP, rs9461045, that has been correlated with altered KIAA0319 expression (Dennis et al., 2009). However, rs9461045 was one of the SNPs from Pinel et al.'s association screen that failed to show significant association with functional asymmetry.

Highly penetrant mutations disrupting FOXP2 lead to disordered speech/language development, associated with subtle but detectable alterations in brain structure and function (Lai et al., 2001; Watkins et al., 2002; Liégeois et al., 2003; Fisher and Scharff, 2009). However, the largest targeted studies of common SNPs spanning the gene, in hundreds of people from the general population, report no association with interindividual differences in either neuroanatomy (Hoogman et al., 2014) or measures of language performance (Mueller et al., 2016). Beyond language, FOXP2 was one of 12 significant loci in a large-scale genome-wide scan of susceptibility to attention deficit/hyperactivity disorder (20,183 cases; 35,191 controls), but the associated SNPs differ from those studied here (Demontis et al., 2019).

Multiple investigations have taken Pinel et al. (2012) as the primary motivation for investigating the highlighted SNPs for associations with questionnaire-based data or behavioral/cognitive tests, but results have not been compelling. In a study of 882 healthy undergraduates who completed schizotypal personality questionnaires, Crespi et al. (2017) reported that rs7799109 (a SNP that cosegregates with rs7784315) was associated with self-report items related to inner speech (p = 0.048) and speech fluency (p = 0.049), and with strength of handedness, independent of direction (p = 0.027). A behavioral study of rs6980093 in ∼200 people reported association with individual differences in accuracy and rate of speech-category learning (Chandrasekaran et al., 2015). However, the distribution of rs6980093 diploid genotypes in that study (proportions of people designated as AA, AG or GG) deviated markedly from the expected Hardy–Weinberg equilibrium (Chandrasekaran et al., 2015). This deviation contradicts all other studies involving rs6980093, as well as the information in public databases of worldwide genotypes, casting doubt on genotyping accuracy and undermining their association results.

In another targeted study, Mozzi et al. (2017) assessed rs6980093 in two Italian samples, reporting associations with semantic fluency in a population-based cohort (699 children; 3–11 years old), and with single-word-reading accuracy in 317 families ascertained for dyslexia (572 children; 6–18 years old). Putative allelic effects on these traits contradicted the additive findings of Pinel et al. (2012) (AA>AG>GG; Fig. 1) and were inconsistent between the Italian samples. In the population-based cohort, the G allele seemed dominant (AA<[AG=GG]), whereas the family cohort appeared to show heterozygous advantage (AG>[AA=GG]), although there is no mechanistic reason to favor this unusual mode of action for rs6980093 (Mozzi et al., 2017). Most recently, Zhang et al. (2018) tested 133 Chinese adults for behavioral and event-related potential responses (N1 and P2) to −50 and −200 cents pitch perturbations during vocal production, reporting that rs6980093 GG carriers performed differently from AA and AG carriers, but only for aspects of the −200 condition. Notably, until now, no study of the SNPs of interest attempted direct replication of Pinel et al. (2012); even when SNP choice was the same, different phenotypes were targeted. Moreover, none have reported fMRI associations.

Although unlikely, we acknowledge that lack of replication in the present study could partly reflect specific aspects of study materials, tasks, and/or cohorts that differed from Pinel et al. (2012). Most participants in both studies were young educated Caucasians, with similar age and sex distributions. One difference between the original design and that of our MOUS cohort, was use of a comprehension task in the latter for ensuring compliance. However, we also find strong nonreplication in EvLab, where memory probes or simple button-press tasks were used. Furthermore, prior work has established that language activations in high-level areas are robust to changes in materials, presentation modality, language, and task (Fedorenko et al., 2010; Scott et al., 2017). We note that if putative genetic effects are real, but so sensitive to subtle aspects of study design that they fail to replicate, then their relevance for understanding biological pathways that are fundamental for language function is minimal.

Future investigations in functional neuroimaging genetics should take into account the typically tiny effects of common polymorphisms when justifying sample size, include appropriate power calculations, address multiple testing robustly, and carefully consider functional relevance of targeted SNPs, given that most variants have no biological impact. Independent replications using directly matching study designs must be encouraged, so that the literature of this emerging field can be built on solid foundations.

Footnotes

  • This work was supported by the Max Planck Society and the Donders Centre for Cognitive Neuroimaging. E.F. was supported by the National Institutes of Health (Grants R00-HD057522, R01-DC016607, and R01-DC016950) and by a grant from the Simons Foundation to the Simons Center for the Social Brain at MIT. E.F. would also like to acknowledge the Athinoula A. Martinos Imaging Center at the McGovern Institute for Brain Research at MIT and its support team (Steve Shannon and Atsushi Takahashi). We thank Matt Siegelman at EvLab and MOUS team members Karl Magnus Petersson, Jan-Mathijs Schoffelen, and Nietzsche Lam for contributions.

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Julia Uddén at Julia.Udden{at}mpi.nl or Simon E. Fisher at simon.fisher{at}mpi.nl.

References

  1. ↵
    1. Becker M,
    2. Devanna P,
    3. Fisher SE,
    4. Vernes SC
    (2018) Mapping of human FOXP2 enhancers reveals complex regulation. Front Mol Neurosci 11:47. doi:10.3389/fnmol.2018.00047 pmid:29515369
    OpenUrlCrossRefPubMed
  2. ↵
    1. Bigos KL,
    2. Hariri AR,
    3. Weinberger DR
    (2016) Neuroimaging genetics: principles and practices. New York: OUP.
  3. ↵
    1. Birn RM,
    2. Cox RW,
    3. Bandettini PA
    (2002) Detection versus estimation in event-related fMRI: choosing the optimal stimulus timing. Neuroimage 15:252–264. doi:10.1006/nimg.2001.0964 pmid:11771993
    OpenUrlCrossRefPubMed
  4. ↵
    1. Blokland GA,
    2. McMahon KL,
    3. Thompson PM,
    4. Martin NG,
    5. de Zubicaray GI,
    6. Wright MJ
    (2011) Heritability of working memory brain activation. J Neurosci 31:10882–10890. doi:10.1523/JNEUROSCI.5334-10.2011 pmid:21795540
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Button KS,
    2. Ioannidis JP,
    3. Mokrysz C,
    4. Nosek BA,
    5. Flint J,
    6. Robinson ES,
    7. Munafò MR
    (2013) Power failure: why small sample size undermines the reliability of neuroscience. Nat Rev Neurosci 14:365–376. doi:10.1038/nrn3475 pmid:23571845
    OpenUrlCrossRefPubMed
  6. ↵
    1. Carrion-Castillo A,
    2. Franke B,
    3. Fisher SE
    (2013) Molecular genetics of dyslexia: an overview. Dyslexia 19:214–240. doi:10.1002/dys.1464 pmid:24133036
    OpenUrlCrossRefPubMed
  7. ↵
    1. Chandrasekaran B,
    2. Yi HG,
    3. Blanco NJ,
    4. McGeary JE,
    5. Maddox WT
    (2015) Enhanced procedural learning of speech sound categories in a genetic variant of FOXP2. J Neurosci 35:7808–7812. doi:10.1523/JNEUROSCI.4706-14.2015 pmid:25995468
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Crespi B,
    2. Read S,
    3. Hurd P
    (2017) Segregating polymorphisms of FOXP2 are associated with measures of inner speech, speech fluency and strength of handedness in a healthy population. Brain Lang 173:33–40. doi:10.1016/j.bandl.2017.06.002 pmid:28609679
    OpenUrlCrossRefPubMed
  9. ↵
    1. Demontis D,
    2. Walters RK,
    3. Martin J,
    4. Mattheisen M,
    5. Als TD,
    6. Agerbo E,
    7. Baldursson G,
    8. Belliveau R,
    9. Bybjerg-Grauholm J,
    10. Baekvad-Hansen M,
    11. Cerrato F,
    12. Chambert K,
    13. Churchhouse C,
    14. Dumont A,
    15. Eriksson N,
    16. Gandal M,
    17. Goldstein JI,
    18. Grasby KL,
    19. Grove J,
    20. Gudmundsson OO, et al
    . (2019) Discovery of the first genome-wide significant risk loci for attention deficit/hyperactivity disorder. Nat Genet 51:63–75. doi:10.1038/s41588-018-0269-7 pmid:30478444
    OpenUrlCrossRefPubMed
  10. ↵
    1. Dennis MY,
    2. Paracchini S,
    3. Scerri TS,
    4. Prokunina-Olsson L,
    5. Knight JC,
    6. Wade-Martins R,
    7. Coggill P,
    8. Beck S,
    9. Green ED,
    10. Monaco AP
    (2009) A common variant associated with dyslexia reduces expression of the KIAA0319 gene. PLoS Genet 5:e1000436. doi:10.1371/journal.pgen.1000436 pmid:19325871
    OpenUrlCrossRefPubMed
  11. ↵
    1. Deriziotis P,
    2. Fisher SE
    (2017) Speech and language: translating the genome. Trends Genet 33:642–656. doi:10.1016/j.tig.2017.07.002 pmid:28781152
    OpenUrlCrossRefPubMed
  12. ↵
    1. Dubois J,
    2. Adolphs R
    (2016) Building a science of individual differences from fMRI. Trends Cogn Sci 20:425–443. doi:10.1016/j.tics.2016.03.014 pmid:27138646
    OpenUrlCrossRefPubMed
  13. ↵
    1. Fedorenko E,
    2. Hsieh PJ,
    3. Nieto-Castañón A,
    4. Whitfield-Gabrieli S,
    5. Kanwisher N
    (2010) New method for fMRI investigations of language: defining ROIs functionally in individual subjects. J Neurophysiol 104:1177–1194. doi:10.1152/jn.00032.2010 pmid:20410363
    OpenUrlCrossRefPubMed
  14. ↵
    1. Fisher SE,
    2. Scharff C
    (2009) FOXP2 as a molecular window into speech and language. Trends Genet 25:166–177. doi:10.1016/j.tig.2009.03.002 pmid:19304338
    OpenUrlCrossRefPubMed
  15. ↵
    1. Flint J,
    2. Munafò MR
    (2007) The endophenotype concept in psychiatric genetics. Psychol Med 37:163–180. doi:10.1017/S0033291706008750 pmid:16978446
    OpenUrlCrossRefPubMed
  16. ↵
    1. Frost MA,
    2. Goebel R
    (2012) Measuring structural–functional correspondence: spatial variability of specialised brain regions after macro-anatomical alignment. Neuroimage 59:1369–1381. doi:10.1016/j.neuroimage.2011.08.035 pmid:21875671
    OpenUrlCrossRefPubMed
  17. ↵
    1. Gasperini M,
    2. Starita L,
    3. Shendure J
    (2016) The power of multiplexed functional analysis of genetic variants. Nat Protoc 11:1782–1787. doi:10.1038/nprot.2016.135 pmid:27583640
    OpenUrlCrossRefPubMed
  18. ↵
    1. Grabitz CR,
    2. Button KS,
    3. Munafò MR,
    4. Newbury DF,
    5. Pernet CR,
    6. Thompson PA,
    7. Bishop DVM
    (2018) Logical and methodological issues affecting genetic studies of humans reported in top neuroscience journals. J Cogn Neurosci 30:25–41. doi:10.1162/jocn_a_01192 pmid:28949821
    OpenUrlCrossRefPubMed
  19. ↵
    GTEx Consortium (2017) Genetic effects on gene expression across human tissues. Nature 550:204–213. doi:10.1038/nature24277 pmid:29022597
    OpenUrlCrossRefPubMed
  20. ↵
    1. Hibar DP,
    2. Stein JL,
    3. Renteria ME,
    4. Arias-Vasquez A,
    5. Desrivières S,
    6. Jahanshad N,
    7. Toro R,
    8. Wittfeld K,
    9. Abramovic L,
    10. Andersson M,
    11. Aribisala BS,
    12. Armstrong NJ,
    13. Bernard M,
    14. Bohlken MM,
    15. Boks MP,
    16. Bralten J,
    17. Brown AA,
    18. Chakravarty MM,
    19. Chen Q,
    20. Ching CR, et al
    . (2015) Common genetic variants influence human subcortical brain structures. Nature 520:224–229. doi:10.1038/nature14101 pmid:25607358
    OpenUrlCrossRefPubMed
  21. ↵
    1. Hoogman M,
    2. Guadalupe T,
    3. Zwiers MP,
    4. Klarenbeek P,
    5. Francks C,
    6. Fisher SE
    (2014) Assessing the effects of common variation in the FOXP2 gene on human brain structure. Front Hum Neurosci 8:473. doi:10.3389/fnhum.2014.00473 pmid:25013396
    OpenUrlCrossRefPubMed
  22. ↵
    1. Hultén A,
    2. Schoffelen JM,
    3. Uddén J,
    4. Lam NHL,
    5. Hagoort P
    (2019) How the brain makes sense beyond the processing of single words: an MEG study. Neuroimage 186:586–594. doi:10.1016/j.neuroimage.2018.11.035 pmid:30481591
    OpenUrlCrossRefPubMed
  23. ↵
    1. Jeffreys H
    (1961) Theory of probability, Ed 3. Oxford: Clarendon.
  24. ↵
    1. Koten JW Jr.,
    2. Wood G,
    3. Hagoort P,
    4. Goebel R,
    5. Propping P,
    6. Willmes K,
    7. Boomsma DI
    (2009) Genetic contribution to variation in cognitive function: an FMRI study in twins. Science 323:1737–1740. doi:10.1126/science.1167371 pmid:19325117
    OpenUrlAbstract/FREE Full Text
  25. ↵
    1. Kriegeskorte N,
    2. Simmons WK,
    3. Bellgowan PS,
    4. Baker CI
    (2009) Circular analysis in systems neuroscience: the dangers of double dipping. Nat Neurosci 12:535–540. doi:10.1038/nn.2303 pmid:19396166
    OpenUrlCrossRefPubMed
  26. ↵
    1. Lai CS,
    2. Fisher SE,
    3. Hurst JA,
    4. Vargha-Khadem F,
    5. Monaco AP
    (2001) A forkhead-domain gene is mutated in a severe speech and language disorder. Nature 413:519–523. doi:10.1038/35097076 pmid:11586359
    OpenUrlCrossRefPubMed
  27. ↵
    1. Lam NHL,
    2. Schoffelen JM,
    3. Uddén J,
    4. Hultén A,
    5. Hagoort P
    (2016) Neural activity during sentence processing as reflected in theta, alpha, beta, and gamma oscillations. Neuroimage 142:43–54. doi:10.1016/j.neuroimage.2016.03.007 pmid:26970187
    OpenUrlCrossRefPubMed
  28. ↵
    1. Liégeois F,
    2. Baldeweg T,
    3. Connelly A,
    4. Gadian DG,
    5. Mishkin M,
    6. Vargha-Khadem F
    (2003) Language fMRI abnormalities associated with FOXP2 gene mutation. Nat Neurosci 6:1230–1237. doi:10.1038/nn1138 pmid:14555953
    OpenUrlCrossRefPubMed
  29. ↵
    1. Mahowald K,
    2. Fedorenko E
    (2016) Reliable individual-level neural markers of high-level language processing: a necessary precursor for relating neural variability to behavioral and genetic variability. Neuroimage 139:74–93. doi:10.1016/j.neuroimage.2016.05.073 pmid:27261158
    OpenUrlCrossRefPubMed
  30. ↵
    1. Mozzi A,
    2. Riva V,
    3. Forni D,
    4. Sironi M,
    5. Marino C,
    6. Molteni M,
    7. Riva S,
    8. Guerini FR,
    9. Clerici M,
    10. Cagliani R,
    11. Mascheretti S
    (2017) A common genetic variant in FOXP2 is associated with language-based learning (dis)abilities: evidence from two Italian independent samples. Am J Med Genet B Neuropsychiatr Genet 174:578–586. doi:10.1002/ajmg.b.32546 pmid:28436202
    OpenUrlCrossRefPubMed
  31. ↵
    1. Mueller KL,
    2. Murray JC,
    3. Michaelson JJ,
    4. Christiansen MH,
    5. Reilly S,
    6. Tomblin JB
    (2016) Common genetic variants in FOXP2 are not associated with individual differences in language development. PLoS One 11:e0152576. doi:10.1371/journal.pone.0152576 pmid:27064276
    OpenUrlCrossRefPubMed
  32. ↵
    1. Oldfield RC
    (1971) The assessment and analysis of handedness: the Edinburgh inventory. Neuropsychologia 9:97–113. doi:10.1016/0028-3932(71)90067-4 pmid:5146491
    OpenUrlCrossRefPubMed
  33. ↵
    1. Pinel P,
    2. Thirion B,
    3. Meriaux S,
    4. Jobert A,
    5. Serres J,
    6. Le Bihan D,
    7. Poline JB,
    8. Dehaene S
    (2007) Fast reproducible identification and large-scale databasing of individual functional cognitive networks. BMC Neurosci 8:91. doi:10.1186/1471-2202-8-91 pmid:17973998
    OpenUrlCrossRefPubMed
  34. ↵
    1. Pinel P,
    2. Fauchereau F,
    3. Moreno A,
    4. Barbot A,
    5. Lathrop M,
    6. Zelenika D,
    7. Le Bihan D,
    8. Poline JB,
    9. Bourgeron T,
    10. Dehaene S
    (2012) Genetic variants of FOXP2 and KIAA0319/TTRAP/THEM2 locus are associated with altered brain activation in distinct language-related regions. J Neurosci 32:817–825. doi:10.1523/JNEUROSCI.5996-10.2012 pmid:22262880
    OpenUrlAbstract/FREE Full Text
  35. ↵
    1. Schoffelen JM,
    2. Hultén A,
    3. Lam N,
    4. Marquand AF,
    5. Uddén J,
    6. Hagoort P
    (2017) Frequency-specific directed interactions in the human brain network for language. Proc Natl Acad Sci U S A 114:8083–8088. doi:10.1073/pnas.1703155114 pmid:28698376
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. Schoffelen JM,
    2. Oostenveld R,
    3. Lam NHL,
    4. Uddén J,
    5. Hultén A,
    6. Hagoort P
    (2019) A 204-subject multimodal neuroimaging dataset to study language processing. Sci Data 6:17. doi:10.1038/s41597-019-0020-y pmid:30944338
    OpenUrlCrossRefPubMed
  37. ↵
    1. Scott TL,
    2. Gallée J,
    3. Fedorenko E
    (2017) A new fun and robust version of an fMRI localizer for the frontotemporal language system. Cogn Neurosci 8:167–176. doi:10.1080/17588928.2016.1201466 pmid:27386919
    OpenUrlCrossRefPubMed
  38. ↵
    1. Simmons JP,
    2. Nelson LD,
    3. Simonsohn U
    (2011) False-positive psychology: undisclosed flexibility in data collection and analysis allows presenting anything as significant. Psychol Sci 22:1359–1366. doi:10.1177/0956797611417632 pmid:22006061
    OpenUrlCrossRefPubMed
  39. ↵
    1. Sollis E,
    2. Deriziotis P,
    3. Saitsu H,
    4. Miyake N,
    5. Matsumoto N,
    6. Hoffer MJV,
    7. Ruivenkamp CAL,
    8. Alders M,
    9. Okamoto N,
    10. Bijlsma EK,
    11. Plomp AS,
    12. Fisher SE
    (2017) Equivalent missense variant in the FOXP2 and FOXP1 transcription factors causes distinct neurodevelopmental disorders. Hum Mutat 38:1542–1554. doi:10.1002/humu.23303 pmid:28741757
    OpenUrlCrossRefPubMed
  40. ↵
    1. Thompson PM,
    2. Cannon TD,
    3. Narr KL,
    4. van Erp T,
    5. Poutanen VP,
    6. Huttunen M,
    7. Lönnqvist J,
    8. Standertskjöld-Nordenstam CG,
    9. Kaprio J,
    10. Khaledy M,
    11. Dail R,
    12. Zoumalan CI,
    13. Toga AW
    (2001) Genetic influences on brain structure. Nat Neurosci 4:1253–1258. doi:10.1038/nn758 pmid:11694885
    OpenUrlCrossRefPubMed
  41. ↵
    1. Tomaiuolo F,
    2. MacDonald JD,
    3. Caramanos Z,
    4. Posner G,
    5. Chiavaras M,
    6. Evans AC,
    7. Petrides M
    (1999) Morphology, morphometry and probability mapping of the pars opercularis of the inferior frontal gyrus: an in vivo MRI analysis. Eur J Neurosci 11:3033–3046. doi:10.1046/j.1460-9568.1999.00718.x pmid:10510168
    OpenUrlCrossRefPubMed
  42. ↵
    1. Uddén J,
    2. Hultén A,
    3. Schoffelen JM,
    4. Lam N,
    5. Harbusch K,
    6. van den Bosch A,
    7. Kempen G,
    8. Petersson KM,
    9. Hagoort P
    (2019) Supramodal sentence processing in the human brain: fMRI evidence for the influence of syntactic complexity in more than 200 participants. bioRxiv:576769.
  43. ↵
    1. Watkins KE,
    2. Vargha-Khadem F,
    3. Ashburner J,
    4. Passingham RE,
    5. Connelly A,
    6. Friston KJ,
    7. Frackowiak RS,
    8. Mishkin M,
    9. Gadian DG
    (2002) MRI analysis of an inherited speech and language disorder: structural brain abnormalities. Brain 125:465–478. doi:10.1093/brain/awf057 pmid:11872605
    OpenUrlCrossRefPubMed
  44. ↵
    1. Zhang S,
    2. Zhao J,
    3. Guo Z,
    4. Jones JA,
    5. Liu P,
    6. Liu H
    (2018) The association between genetic variation in FOXP2 and sensorimotor control of speech production. Front Neurosci 12:666. doi:10.3389/fnins.2018.00666 pmid:30294257
    OpenUrlCrossRefPubMed
Back to top

In this issue

The Journal of Neuroscience: 39 (44)
Journal of Neuroscience
Vol. 39, Issue 44
30 Oct 2019
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Advertising (PDF)
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Toward Robust Functional Neuroimaging Genetics of Cognition
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Toward Robust Functional Neuroimaging Genetics of Cognition
Julia Uddén, Annika Hultén, Katarina Bendtz, Zachary Mineroff, Katerina S. Kucera, Arianna Vino, Evelina Fedorenko, Peter Hagoort, Simon E. Fisher
Journal of Neuroscience 30 October 2019, 39 (44) 8778-8787; DOI: 10.1523/JNEUROSCI.0888-19.2019

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Toward Robust Functional Neuroimaging Genetics of Cognition
Julia Uddén, Annika Hultén, Katarina Bendtz, Zachary Mineroff, Katerina S. Kucera, Arianna Vino, Evelina Fedorenko, Peter Hagoort, Simon E. Fisher
Journal of Neuroscience 30 October 2019, 39 (44) 8778-8787; DOI: 10.1523/JNEUROSCI.0888-19.2019
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • fMRI
  • FoxP2
  • individual differences
  • KIAA0319/TTRAP/THEM2
  • language
  • neuroimaging genetics

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Articles

  • Identification of the Acid-Sensitive Site Critical for Chloral Hydrate (CH) Activation of the Proton-Activated Chloride Channel
  • Disentangling Object Category Representations Driven by Dynamic and Static Visual Input
  • Irrelevant Threats Linger and Affect Behavior in High Anxiety
Show more Research Articles

Behavioral/Cognitive

  • Disentangling Object Category Representations Driven by Dynamic and Static Visual Input
  • Irrelevant Threats Linger and Affect Behavior in High Anxiety
  • Multisession Anodal Transcranial Direct Current Stimulation Enhances Adult Hippocampal Neurogenesis and Context Discrimination in Mice
Show more Behavioral/Cognitive
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
(JNeurosci logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.