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 in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • 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
Research Articles, Behavioral/Cognitive

Individual Differences in Reading Skill Are Related to Trial-by-Trial Neural Activation Variability in the Reading Network

Jeffrey G. Malins, Kenneth R. Pugh, Bonnie Buis, Stephen J. Frost, Fumiko Hoeft, Nicole Landi, W. Einar Mencl, Anish Kurian, Ryan Staples, Peter J. Molfese, Rose Sevcik and Robin Morris
Journal of Neuroscience 21 March 2018, 38 (12) 2981-2989; DOI: https://doi.org/10.1523/JNEUROSCI.0907-17.2018
Jeffrey G. Malins
1Haskins Laboratories, New Haven, Connecticut 06511,
2Department of Pediatrics, Yale University School of Medicine, New Haven, Connecticut 06520,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Jeffrey G. Malins
Kenneth R. Pugh
1Haskins Laboratories, New Haven, Connecticut 06511,
3Department of Psychological Sciences, University of Connecticut, Storrs, Connecticut 06269,
4Department of Linguistics, Yale University, New Haven, Connecticut 06511,
5Department of Diagnostic Radiology, Yale University School of Medicine, New Haven, Connecticut 06520,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Bonnie Buis
1Haskins Laboratories, New Haven, Connecticut 06511,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Stephen J. Frost
1Haskins Laboratories, New Haven, Connecticut 06511,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Fumiko Hoeft
1Haskins Laboratories, New Haven, Connecticut 06511,
6Department of Psychiatry and Weill Institute for Neurosciences, University of California, San Francisco, California 94143,
7Dyslexia Center, University of California, San Francisco, California 94143,
8UC6-Stanford Precision Learning Center Neurosciences, San Francisco, California 94143,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Fumiko Hoeft
Nicole Landi
1Haskins Laboratories, New Haven, Connecticut 06511,
3Department of Psychological Sciences, University of Connecticut, Storrs, Connecticut 06269,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Nicole Landi
W. Einar Mencl
1Haskins Laboratories, New Haven, Connecticut 06511,
4Department of Linguistics, Yale University, New Haven, Connecticut 06511,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Anish Kurian
1Haskins Laboratories, New Haven, Connecticut 06511,
3Department of Psychological Sciences, University of Connecticut, Storrs, Connecticut 06269,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Anish Kurian
Ryan Staples
1Haskins Laboratories, New Haven, Connecticut 06511,
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Peter J. Molfese
1Haskins Laboratories, New Haven, Connecticut 06511,
3Department of Psychological Sciences, University of Connecticut, Storrs, Connecticut 06269,
9Section on Functional Imaging Methods, Laboratory of Brain and Cognition, National Institutes of Mental Health, National Institutes of Health, Bethesda, Maryland 20892, and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Peter J. Molfese
Rose Sevcik
10Department of Psychology, Georgia State University, Atlanta, Georgia 30303
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Robin Morris
10Department of Psychology, Georgia State University, Atlanta, Georgia 30303
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Recent work has suggested that variability in levels of neural activation may be related to behavioral and cognitive performance across a number of domains and may offer information that is not captured by more traditional measures that use the average level of brain activation. We examined the relationship between reading skill in school-aged children and neural activation variability during a functional MRI reading task after taking into account average levels of activity. The reading task involved matching printed and spoken words to pictures of items. Single trial activation estimates were used to calculate the mean and standard deviation of children's responses to print and speech stimuli; multiple regression analyses evaluated the relationship between reading skill and trial-by-trial activation variability. The reliability of observed findings from the discovery sample (n = 44; ages 8–11; 18 female) was then confirmed in an independent sample of children (n = 32; ages 8–11; 14 female). Across the two samples, reading skill was positively related to trial-by-trial variability in the activation response to print in the left inferior frontal gyrus pars triangularis. This relationship held even when accounting for mean levels of activation. This finding suggests that intrasubject variability in trial-by-trial fMRI activation responses to printed words accounts for individual differences in human reading ability that are not fully captured by traditional mean levels of brain activity. Furthermore, this positive relationship between trial-by-trial activation variability and reading skill may provide evidence that neural variability plays a beneficial role during early reading development.

SIGNIFICANCE STATEMENT Recent work has suggested that neural activation variability, or moment-to-moment changes in the engagement of brain regions, is related to individual differences in behavioral and cognitive performance across multiple domains. However, differences in neural activation variability have not yet been evaluated in relation to reading skill. In the current study, we analyzed data from two independent groups of children who performed an fMRI task involving reading and listening to words. Across both samples, reading skill was positively related to trial-by-trial variability in activation to print stimuli in the left inferior frontal gyrus pars triangularis, even when accounting for the more conventional measure of mean levels of brain activity. This finding suggests that neural variability could be beneficial in developing readers.

  • β series
  • BOLD variability
  • event-related fMRI
  • individual differences
  • reading disability
  • trial-by-trial variability

Introduction

A growing body of neuroimaging research has linked reading skill to variation in structural and functional circuitry in the brain (Norton et al., 2015). Investigations concerning the functional neuroanatomy of reading have focused on mean levels of activation across trials while children are engaged in different reading tasks. However, an emerging literature suggests that mean differences in activation and connectivity reflect only part of the complex neural foundation of reading ability. Recent studies have linked reading skill to the stability of neural responses to speech sounds (Hornickel and Kraus, 2013). In addition, animal work has shown that expression of the rat homolog of the dyslexia susceptibility gene KIAA0319 is linked to increased trial-by-trial variability in speech sound responses (Centanni et al., 2014a,b). Together, these studies have helped to support the neural noise hypothesis of reading disability, which postulates that levels of neural noise can influence timing mechanisms that affect signal variability and thereby affect reading performance (Hancock et al., 2017).

This previous work leads to an expectation that reading skill in children is related to within-subject measures of neural activation variability. However, to date, neural activation variability has been evaluated with respect to reading skill in children only by examining brainstem electrophysiological responses to speech sounds (Hornickel and Kraus, 2013). In the current study, we instead asked children to perform a task involving word reading and examined trialwise variability in cortical activation using fMRI, a technique that has been used successfully to examine the relationship between neural activation variability and behavioral performance in multiple domains outside of reading (Garrett et al., 2013).

Across these other domains, there exists some debate concerning whether increased variability in the blood oxygen level-dependent (BOLD) signal confers a positive or negative impact on behavior. The directionality of the effect appears in part to be related to the extent to which a task entails cognitive versus sensory processing. For example, increased BOLD signal variability has been associated with faster and more consistent reaction times (RTs) in younger versus older adults during cognitive tasks including attentional cueing and delayed match-to-sample (Garrett et al., 2011), whereas in a study examining cognitive flexibility and stability (Armbruster-Genc et al., 2016), the direction of the relationship between BOLD signal variability and cognitive performance has been characterized as positive or negative depending on the task. Conversely, for sensory processing, increased BOLD signal variability has been associated with increased behavioral variability in older compared with younger adults during audiovisual speech perception (Baum and Beauchamp, 2014) and has been considered maladaptive in adults with autism, who showed greater trial-to-trial variability compared with matched controls in primary sensory areas during a low-level sensory task (Haigh et al., 2016). Given that reading involves both sensory and cognitive components, the direction of the relationship between BOLD signal variability and reading skill therefore remains an open question.

In the current study, we addressed the following novel questions: (1) Does trial-by-trial neural activation variability account for variance in reading skill in children above and beyond differences in mean activation? (2) If so, what is the direction of the relationship between neural activation variability and reading skill? To address these questions, we first conducted analyses on fMRI data from a discovery sample of children who performed a task in which they judged whether printed or spoken words matched pictures of items. We then confirmed whether observed effects held in a separate, independent sample of children. Analyses focused on using single trial β estimates to quantify mean activation across trials, as well as trial-by-trial variability, in the evoked response to print and speech within regions of the reading network. These mean and variability measures were then entered into multiple regression models characterizing the manner in which trial-by-trial activation variability is associated with reading skill after accounting for mean task-related activation as well as predictors of noninterest such as subject age.

Materials and Methods

Discovery sample

Participants.

Children were selected from a larger study examining response to intervention for reading disability; the data presented here correspond to baseline scans before the onset of any intervention. Of this larger sample of 82 children, 44 were selected who met the following inclusion criteria: in third or fourth grade (71/82; the other 11 participants belonged to a cohort of seventh and eighth graders who participated in the larger study), an average Euclidean movement of 0.25 mm or less (58/71), and at least 70% accuracy in each of the auditory and visual mismatch conditions (44/58). Euclidean movement was calculated per volume by first computing point-to-point change for each the six motion parameters (i.e., three translation and three rotation) and then taking the square root of the sum of squares of these measures; average Euclidean movement was calculated by taking the mean value of this measure across all volumes of data collection. The accuracy cutoffs were selected to have a sufficient number of correct trials per participant to calculate dependable standard deviation (SD) measurements. The motion cutoff was selected to increase power to detect effects related to differences in trial-by-trial variability that are unrelated to motion because we expected that intrasubject variability would be affected by movement in the scanner (Lund et al., 2005). We acknowledge that the percentage of data lost is larger than comparable fMRI studies with pediatric populations; however, data quality criteria were particularly stringent for this investigation given our concerns regarding participant motion as well as the requirement of having a sufficient number of trials to calculate valid SD measurements.

All children completed a battery of standardized cognitive assessments (Table 1). These included assessments of single word reading, pseudoword decoding, and passage comprehension from the Woodcock-Johnson III Tests of Achievement (WJ-III; Woodcock et al., 2001); the Peabody Picture Vocabulary Test (PPVT-4; Dunn and Dunn, 2007), which measures receptive vocabulary; the Comprehensive Test of Phonological Awareness (CTOPP-2; Wagner et al., 2013), which measures metalinguistic knowledge of the structure of speech, or phonological awareness, by assessing skills including phoneme elision, blending, and isolation; and the Wechsler Abbreviated Scale of Intelligence (WASI-II; Wechsler, 2011), which measures verbal and nonverbal intelligence. As can be observed in Table 1, the range of reading scores was very broad and some children in the sample would be considered typically developing, whereas others would be classified as having reading disability using traditional diagnostic criteria. However, we treated reading skill as a continuous dimension, consistent with recent views concerning the multifactorial nature of reading skill, as well as the pitfalls of grouping children into diagnostic categories using cutoff scores (Pennington et al., 2012; Branum-Martin et al., 2013).

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

Descriptive information concerning the two groups of children who performed the fMRI experiment

fMRI task.

Functional volumes were acquired while participants completed a task in which they judged whether picture cues matched auditory and visual target words (Frost et al., 2009; Landi et al., 2013; Jasińska et al., 2016; Preston et al., 2016). In this task, participants were presented with pictures of common items (e.g., “cake”) that remained on the screen for 40–65 s, corresponding to between seven and eight trials, before being replaced by another picture. This procedure encouraged participants to generate strong expectations of target items, thereby maximizing responses to mismatches, and also obviated the need to associate targets with a new picture on every single trial, which could have been overly taxing. Whereas each picture remained on the screen, participants were presented with target items in an event-related fashion; specifically, printed words appeared in a box below the picture (presented for 3000 ms in 40-point Arial font), or auditory words were presented via headphones. Importantly, in 1/6 of trials, the printed or spoken word matched the picture, whereas in the other 5/6 of trials, the printed or spoken word mismatched the picture. Participants were asked to indicate via button press whether the printed or spoken word matched the picture. In total, participants completed 25 trials in each of the auditory (spoken) and visual (print) mismatch conditions. A sample trial sequence is illustrated in Figure 1.

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

A sample trial sequence for the fMRI task.

Acquisition of MRI data.

Images were acquired using a 3 T Siemens Trio scanner with a 12-channel head coil located at the GSU/GaTech Center for Advanced Brain Imaging in Atlanta, Georgia. T2*-weighted images were acquired in an axial-oblique orientation parallel to the intercommissural line (32 slices; 4 mm slice thickness; no gap) using single-shot echo planar imaging (matrix size = 64 × 64; voxel size = 3.438 × 3.438 × 4 mm; FoV = 220 mm; TR = 2000 ms; TE = 30 ms; flip angle = 80°). To allow for stabilization of the magnetic field, the first four volumes within each run were discarded. Anatomical scans were collected in the same orientation as the functional volumes (MPRAGE; matrix size = 256 × 256; voxel size = 1 × 1 × 1 mm; FoV = 256 mm; TR = 2530 ms; TE = 2.77 ms; flip angle = 7°); these were acquired either after or between the functional runs. In total, participants completed two runs of the functional task, which had a combined duration of 7 min 32 s (226 volumes). Across all trials in the experiment, the time between trial onsets was jittered between 4 and 13 s; trial order and intertrial intervals were optimized by an in-house MATLAB program that balanced intertrial intervals and null trials across conditions and minimized the variability of the measured response in Monte Carlo simulations.

Preprocessing.

Data were analyzed using AFNI (Cox, 1996; RRID:SCR_005927). Functional images were preprocessed by first correcting for slice acquisition time (3dTshift). After this, functional images were aligned with anatomical images, corrected for motion using a six-parameter rigid-body transform (3dvolreg), and normalized to the Colin27 brain in Talairach space using an affine transform (@auto_tlrc). These three steps were combined into a single transform that also forced a 3 mm isotropic voxel size on the data. All images were then smoothed (3dmerge) using a Gaussian kernel with a full width at half maximum of 8 mm (i.e., twice the between-plane distance of 4 mm; Skudlarski et al., 1999) and data were scaled (3dcalc) so that each voxel's time series had a mean of 100 for each run. During this scaling step, values in excess of 200 were clipped; this is the default value for scaling in AFNI and was selected to retain the precision of scaled short values.

We elected to use the Talairach atlas for normalization because Burgund et al. (2002) have shown that, relative to the resolution of fMRI data, there are minimal anatomical differences between children ages 7 and 8 compared with adults. Given that the children in the current study were even older than the children in the Burgund et al. (2002) study (i.e., between 8 and 12 years of age), our view is that use of the Talairach atlas should allow for broader comparability between our study and others, including developmental investigations with adult samples.

GLM analysis.

Single trial β estimates were obtained using a single General Linear Model (GLM) including nuisance regressors for the six motion parameters as well as a separate regressor for each trial (Rissman et al., 2004; Mumford et al., 2012; least-squares all, or LS-A). This model was specified using the -stim_times_IM flag for 3dDeconvolve in AFNI. The hemodynamic response function was approximated using a gamma function. Because we were interested in intrinsic neural variability as opposed to variability related to individual differences in behavioral performance on the task, we included reaction times (RTs) for each trial as duration modulators in the GLM (Grinband et al., 2008; Yarkoni et al., 2009). For trials in which a participant either did not respond, responded with a RT <200 ms (i.e., invalid anticipation), or responded with an RT >1.5 times the interquartile range above the third quartile for a participant's distribution of RTs, overall mean RT for that participant was used as a duration modulator; however, these trials were not considered in further analyses.

β estimates corresponded to the amplitude assigned to each regressor in the GLM, and the set of β estimates across trials for a given voxel constituted that voxel's β series. When performing the GLM, any volume that exceeded the thresholds of 0.3 mm Euclidean movement and/or 10% outliers were censored from further analysis, resulting in an average loss of less than one trial in each of the auditory and visual mismatch conditions. It should be noted that this approach gave rise to some extreme outlier β values due to rare spikes that were still present in the data even after censoring these volumes. To handle these, outlier β values were identified for each participant using the program 3dToutcount in AFNI, which flags outliers using an algorithm based on median absolute deviation. Trials with outliers in >10% of voxels in the brain were censored from analysis; this occurred for an average of two trials in each of the auditory and visual mismatch conditions. In all other trials, outlier values were replaced with zeroes and ignored when calculating average β values within regions of interest (the mean number of voxels with outlier values across all trials and participants was less than one in both the auditory and visual mismatch conditions in each of the ROIs detailed in the next section).

ROI selection.

Given that mean activation was one of the predictors we aimed to include in the multiple regression models evaluating relationships with reading skill, we elected not to analyze mean activation at the whole brain level because this would have biased selection of ROIs. Moreover, we did not predict perfect concordance between areas in which the reading task resulted in overall levels of activation and areas in which the task resulted in increased levels of variance in activation. Therefore, we instead defined ROIs using a recent meta-analysis that took the results of 20 different imaging studies of reading in children and combined them to identify a set of coordinates which showed convergence across studies (Martin et al., 2015). Because this meta-analysis combined results across tasks examining different aspects of reading, our view was that, by using these coordinates, we were more likely to include regions that may show a relationship between neural activation variability and reading skill even if these regions do not appear in a map of mean activation for the current task; with that said, we acknowledge the limitation that the meta-analysis also used mean activation to define ROIs.

We created spheres with a radius of 6 mm (two voxels) centered on the Martin et al. (2015) coordinates for the following regions: left inferior frontal gyrus (IFG) pars opercularis, left IFG pars triangularis, left middle temporal gyrus, left superior temporal gyrus (STG), left superior parietal lobule, and left inferior temporal gyrus. In addition, we also included an ROI for left thalamus, given extant findings indicating that the thalamus contributes to the reading network (Galaburda et al., 1985; Pugh et al., 2013). The full set of ROIs selected for analysis are detailed in Table 3 and displayed in Figure 2. For reference, we overlaid these ROIs on a conjunction map that shows the extent of overlap in task-related activation between the discovery sample and the confirmation sample (further details below). This map was created by running a standard GLM for each participant with a single regressor per condition; groupwise evoked response maps across all task conditions were then generated for both the discovery and confirmation samples using the program 3dANOVA2 (corrected at a false discovery rate (FDR) < 0.01). The conjunction map was created using step functions (3dcalc) and adding together resultant maps.

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

Lateral and medial views of the left hemisphere showing the location of each of the regions of interest (ROIs) tested in the multiple regression analyses. These ROIs, shown in white, are peak coordinates from the Martin et al. (2015) meta-analysis and are overlaid on a conjunction map that shows the overlap in task-related activation between the discovery and confirmation samples (each corrected at a false discovery rate (FDR) < 0.01; for more details on how this map was constructed, refer to the Materials and Methods section). IFGtr, Inferior frontal gyrus pars triangularis; IFGop, inferior frontal gyrus pars opercularis; SPL, superior parietal lobule; MTG, middle temporal gyrus; ITG, inferior temporal gyrus.

Analysis of trialwise variability.

For analysis of trialwise variability, we considered correct trials in the auditory and visual word mismatch conditions with RTs within the acceptable range. Our rationale for analyzing both print and speech trials was that, even though we were specifically interested in responses to printed words, the paradigm included spoken words and analyzing neural responses in this condition afforded us the ability to determine whether any potential relationships between neural activation variability and reading skill were print specific or were instead more general for language. Our rationale for analyzing only the mismatch conditions was that these were the predominant conditions in the experiment in terms of the overall number of trials; related to this, there were too few trials in the match conditions to calculate valid SD measurements. For each trial and each ROI, we calculated the average β weight across the voxels in the ROI (3dROIstats), ignoring outlier voxels that had been replaced with zero.

Next, we calculated the mean and SD of the β series in each ROI. Intrasubject SD measures were calculated by using leave-one-out jackknife estimation in version 2015.2 of the package “bootstrap” (Tibshirani and Leisch, 2015) in the R Project for Statistical Computing (RRID:SCR_001905) and taking the mean across estimates. Jackknife estimation was used to mitigate bias of SD estimates, especially given the relatively small number of measurements from which these SDs were derived (Efron, 1981). Then, in each ROI, we ran separate multiple regression models for the auditory and visual mismatch conditions with reading skill as the dependent variable, which was quantified using raw scores for Letter–Word Identification (LWID) from the WJ-III Tests of Achievement. We started with a full model that included the mean and SD of the β series in either the auditory or visual mismatch condition, as well as the following predictors of noninterest: age in months (McIntosh et al., 2008; Garrett et al., 2011); amount of subject motion, defined as the average point-to-point Euclidean movement across all volumes of data collection (Power et al., 2012); and the number of trials used to calculate the mean and SD of the β series (i.e., the number of correct trials following removal of trials that exceeded motion, outlier, or RT thresholds). Using the program dropterm in version 7.43–45 of the “MASS” package in R (Venables and Ripley, 2002), we removed, in a stepwise fashion, any of the three predictors of noninterest that did not account for significant variance in reading skill (an α criterion of 0.05 was used for backward selection; at each step, the predictor with the largest associated p-value was removed). Then, for the resulting models, change in the Akaike Information Criterion (AIC), change in the Bayesian Information Criterion (BIC), and change in adjusted R2 were quantified for both the mean and SD of the β series by comparing final models with models in which each of these respective terms were removed.

Code Accessibility.

Analysis scripts can be accessed via the Open Science Framework at: https://osf.io/6zrak/.

Confirmation sample

Participants.

Children were selected for this analysis from a large dataset that has been the subject of other reports (Frost et al., 2009; Landi et al., 2013; Jasińska et al., 2016; Preston et al., 2016). From this larger sample of 122 children, we first selected participants whose average Euclidean movement was 0.25 mm or less (81/122). Next, because the distribution of reading ability in the larger sample differed from the discovery sample, which was weighted toward the lower end of the reading skill distribution, we selected a subset of children who were matched to the discovery sample in age and raw single word reading scores (WJ-III LWID) using version 3.0.1 of the R package “MatchIt” (Ho et al., 2011). From this subset of children, we then selected those who attained at least 70% accuracy in each of the auditory and visual mismatch conditions, which resulted in 32 children in the confirmation sample (14 female). Assessment scores for the confirmation sample are listed in Table 1. As can be noted from the table, mean raw reading scores and mean age were not significantly different across the two samples (WJ-III LWID raw scores: t(55) = −1.51, p = 0.14; age: t(47) = −0.30, p = 0.77), although standard single word reading scores, phonological awareness, vocabulary, and IQ were lower in the discovery sample compared with the confirmation sample (WJ-III LWID standard scores: t(64) = −2.26, p = 0.03; CTOPP phonological awareness composite standard scores: t(57) = −5.51, p < 0.001; PPVT standard scores: t(72) = −3.39, p = 0.001; WASI FSIQ-2 standard scores: t(68) = −3.15, p < 0.01), with the caveat that children in the confirmation sample were tested using CTOPP (Wagner, Torgesen, & Rashotte, 1999) instead of CTOPP-2, PPVT-III (Dunn & Dunn, 1997) instead of PPVT-4, and WASI-I (Wechsler, 1999) instead of WASI-II. This caveat notwithstanding, we would argue that such sample differences provide for increased generalizability of results.

fMRI task.

Functional volumes were acquired while participants completed the same picture cue–target word identification task as the children in the discovery sample. However, the task in this sample included a larger number of conditions. More specifically, for both the auditory and visual modalities, mismatches were either real words or pseudowords; in addition, for the visual modality, some mismatches were either semantically related words or meaningless consonant strings. As a result of this different design, the match to mismatch ratio was 1:4 instead of 1:5; in addition, printed words were presented for a duration of 2000 ms and in 18-point Verdana font. To keep the confirmatory analyses as comparable as possible to the analyses used for the discovery sample, only the real word conditions in both modalities were considered.

Acquisition of MRI data.

Images were acquired using a 1.5 T Siemens Sonata scanner with a one-channel head coil located at the Yale Magnetic Resonance Research Center in New Haven, Connecticut. T2*-weighted images were acquired in an axial-oblique orientation parallel to the intercommissural line (20 slices; 6 mm slice thickness; no gap) using single-shot echoplanar imaging (matrix size = 64 × 64; voxel size = 3.125 × 3.125 × 6 mm; FoV = 200 mm; TR = 2000 ms; TE = 50 ms; flip angle = 80°). To allow for stabilization of the magnetic field, the first four volumes within each run were discarded. Anatomical scans were collected in a sagittal orientation (MPRAGE; matrix size = 256 × 256; voxel size = 1 × 1 × 1 mm; FoV = 256 mm; TR = 2000 ms; TE = 3.65 ms; flip angle = 8°); these were acquired either following or between the functional runs. Participants completed between 7 and 10 functional runs each 3:46 (113 volumes) in length, which corresponded to up to 40 trials in each of the auditory and visual mismatch conditions (i.e., 4 trials in each condition in each run). Across all trials in the experiment, the time between trial onsets was jittered between 4 and 13 s.

Analysis pipeline.

All analyses were conducted in the exact same fashion as they were for the discovery sample. Removal of volumes that exceeded the thresholds of 0.3 mm Euclidean movement and/or 10% outliers before the GLM resulted in an average loss of less than one trial in each of the auditory and visual mismatch conditions. After the GLM, removal of trials with outliers in >10% of voxels in the brain resulted in a further loss of three trials on average in the auditory mismatch condition and two trials on average in the visual mismatch condition; subsequently, the mean number of voxels with outlier values across all trials and participants was less than one in both the auditory and visual mismatch conditions in each ROI. Analyses focused solely on the ROIs and experimental conditions for which we observed effects in the discovery sample for either the mean or SD of the β series. When performing these confirmatory analyses, we opted to use a Bonferroni-corrected α threshold for significance of 0.0125, which was calculated by dividing 0.05 by 4, the total number of models tested.

Results

Behavioral performance for the in-scanner task for both samples is summarized in Table 2, along with data concerning average amount of movement in the scanner and the number of trials in the β series in each condition. Average values across participants for the mean and SD of the β series within each ROI are presented in Table 3; results for the multiple regression analysis are detailed in Table 4 for the discovery sample and Table 5 for the confirmation sample.

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

MRI quality control parameters and performance for the in-scanner picture cue–target word identification task

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

Location of the regions of interest selected for the multiple regression analysis and average mean and standard deviation of the beta series across participants

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

Multiple regressions for the discovery sample. Each model quantifies the contribution of mean activation and trial-by-trial activation variability to WJ-III Letter Word Identification raw scores

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

Multiple regressions for the confirmation sample. Each model quantifies the contribution of mean activation and trial-by-trial activation variability to WJ-III Letter Word Identification raw scores

Across the two samples, we observed a positive relationship between reading skill and trial-by-trial neural activation variability for printed words in the left IFG pars triangularis (Fig. 3). In this region, trial-by-trial neural activation variability for printed words not only accounted for significant variance in reading skill above and beyond mean levels of activation, but actually accounted for a greater proportion of variance in reading skill than did mean activation. The relationship between reading skill and neural activation variability in this region appears to be fairly selective for print because we did not observe a significant relationship between reading skill and trial-by-trial variability in neural activation for spoken words in this region. For spoken words, the only relationship that we observed between reading skill and trial-by-trial neural activation variability was a negative association in the left STG in the discovery sample; however, this finding did not hold in the confirmation sample. These results also appear to be fairly selective for reading ability because we ran a secondary analysis with performance IQ (WASI Matrix Reasoning raw scores) as the dependent variable to assess whether neural activation variability was related to general cognitive ability. These analyses did not reveal any significant relationship between performance IQ and neural activation variability for print or speech in any region of interest. In addition, we also reran the analysis with single word reading skill as the dependent variable, this time including performance IQ as a covariate in the multiple regression models. The positive relationship between reading skill and neural activation variability for print in the left IFG pars triangularis was still marginally significant in both the discovery and confirmation samples, even when accounting for individual differences in performance IQ (discovery sample: p = 0.055; confirmation sample: p = 0.017).

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

Partial correlation results for the multiple regression analyses performed for the left inferior frontal gyrus pars triangularis region of interest. The x-axis specifies adjusted trial-by-trial neural activation variability, which corresponds to residuals from a model where the standard deviation of the β series is the dependent variable and the mean of the β series, as well as any significant predictors of noninterest (i.e., age, average Euclidean motion, and/or the number of trials in the β series) are the regressors; the y-axis specifies adjusted reading scores, which are the size of the residuals from a model where letter–word identification raw scores are the dependent variable and the regressors are the mean of the β series as well as the same predictors of noninterest. The regression line is superimposed on the plot; the shaded region represents the 95% confidence interval.

The above results stand in contrast to those for mean activation, which only showed a significant relationship between reading skill and activation for printed and spoken words in the left thalamus in the discovery sample but not the confirmation sample. We should also note that the alternative predictors did not show systematic patterns across both samples; however, we did observe that average Euclidean motion and the number of trials within the β series accounted for significant variance in reading skill in a number of regions. To further test for the influence of subject motion and the number of trials in the β series on the observed results, we ran a secondary analysis in which we relaxed the subject inclusion criteria to a maximum of 0.40 mm average Euclidean motion and a minimum of 50% accuracy for both the auditory and visual mismatch conditions. This resulted in a sample size of 50 for the discovery sample and 43 for the confirmation sample. The relationship between reading skill and neural activation variability for print in the left IFG pars triangularis was still significant in the same direction in both samples (discovery sample: p = 0.012; confirmation sample: p = 0.002).

Discussion

Our aim was to assess the relationship between reading skill in school-aged children and trial-by-trial variability in fMRI activation for print or speech. This stemmed from recent advances concerning individual differences in neural response stability in relation to reading skill (Hornickel and Kraus, 2013), as well as the potential impact of neural noise on the timing and systematic variability of processes important for reading (Hancock et al., 2017). Based on this previous work, we hypothesized that reading skill would be related to trial-by-trial variability in neural activation even after accounting for intrasubject differences in mean levels of task-related activation. However, the direction of this relationship remained an open question because domains outside of reading have shown different relationships between behavioral performance and variability in the fMRI BOLD response, whether measured from moment-to-moment within blocks (Garrett et al., 2011) or from trial-to-trial in experiments using event-related designs (Baum and Beauchamp, 2014; Armbruster-Genc et al., 2016; Haigh et al., 2016).

Trial-by-trial activation variability versus mean activation

For each of two samples of children who performed an fMRI picture-word matching task, we entered intrasubject means and SDs of single trial β estimates for print and speech trials into multiple regression models predicting reading skill. We observed that in the left inferior frontal gyrus pars triangularis, the SD of the β series for printed words not only accounted for additional variance in reading skill that was not captured by the mean of the β series, but actually accounted for a greater proportion of variance in reading skill than did mean levels of activation. This effect held across the two samples despite differences in participants, scanners, and slight differences in trial context. Moreover, these effects were fairly selective for print because we did not observe a relationship between reading skill and activation variability for spoken words in this region; the only effect we observed for spoken words was a negative relationship between reading skill and activation variability in the left STG in the discovery sample that was not observed in the confirmation sample. In contrast to the findings for activation variability, for mean activation, the only relationships that we observed were positive correlations between reading skill and mean activation for both printed and spoken words in the left thalamus in the discovery sample, supporting previous studies documenting the important contributions of the thalamus to reading (Galaburda et al., 1985; Pugh et al., 2013); however, these relationships were not observed in the confirmation sample.

The left IFG has long been implicated as a critical part of the skilled reading network in adults, with more anterior and lateral subregions of IFG such as pars triangularis thought to be involved in semantic processing (Poldrack et al., 1999; Bookheimer, 2002). Moreover, activation of the left IFG, as well as connectivity between the left IFG and certain components of the reading network, has been associated with age-related increases over the course of reading development and this region has been linked to processes such as phonological segmentation and covert articulation (Schlaggar et al., 2002; Turkeltaub et al., 2003; Bitan et al., 2007). However, the current study is the first time that individual differences in reading ability in children have been associated with variability in fMRI activation for printed words in this region.

Direction of the relationship between reading skill and BOLD signal variability

In the left IFG pars triangularis, trial-by-trial variability in neural activation for print was positively related to reading skill. A positive relationship between BOLD signal variability and behavioral performance has been previously observed by Garrett et al. (2011), who found that increased levels of BOLD signal variability were associated with faster and more consistent behavioral performance in younger versus older adults across a range of cognitive tasks, as well as by Armbruster-Genc et al. (2016), who found that increased levels of BOLD signal variability were associated with greater cognitive flexibility in adults, which manifested as reduced behavioral switching costs in a task-switching paradigm.

Based on EEG measures in children, McIntosh et al. (2008) suggest that increased neural variability reflects a greater dynamic range of cognitive states as well as a greater ability to transition between them, which perhaps translates to a greater ability to adapt to the environment. It is possible that the increased neural variability observed in the better readers in the current study could reflect greater neural adaptability; however, we interpret this with caution given that we assessed trial-by-trial variability as opposed to moment-to-moment variability at a finer within-trial timescale. Furthermore, increased neural variability may not always have a positive effect on behavioral performance. Armbruster-Genc et al. (2016) observed that one of the brain regions that showed a positive association between neural variability and cognitive flexibility—that is, the left inferior frontal junction—actually showed a negative association with cognitive stability, which manifested as more extensive behavioral costs for distractor inhibition. This dissociation in terms of directionality may depend on the level of hierarchical organization in which a brain region is situated as well as the extent to which a task is weighted toward sensory versus cognitive processing. For example, individuals with autism, a neurodevelopmental disorder that can co-occur with reading disability, have shown increased fMRI BOLD signal variability in primary sensory areas in response to low-level sensory stimulation and this finding has been used to explain why individuals with autism may experience difficulties in highly sensory environments (Haigh et al., 2016). This distinction between sensory and cognitive processing may help to reconcile the current results with the observation that low-level neural responses to speech sounds show greater variability in children with reading disability compared with typically developing children (Hornickel and Kraus, 2013). In addition, these findings, as well as future experiments that more directly tease apart sensory versus cognitive processing, may inform the neural noise hypothesis put forth by Hancock et al. (2017) by elucidating the conditions that promote greater versus lesser neural variability in developing readers as well as how these relationships pattern across different brain regions as a function of reading experience. The differentiation of the role of random neural “noise” versus systemic components that drive “dynamic range” or “adaptability” within such greater neural variability indices may be a critical conceptual and analytic challenge.

In the electrophysiological literature, He and Zempel (2013) have asserted that a certain amount of neural variability is beneficial, but if the level of variability is too high, brain activity could be scattered across too wide a range, which could be detrimental. Therefore, it is possible that the relationship between neural variability and behavioral performance is non-monotonic, and that we are only observing the ascending portion of an inverted U-shaped curve. Future studies could address this possibility by including children who are more severely impaired than the poorest readers in the present sample.

Possible mechanisms of neural variability

The increased neural variability observed in the more skilled readers in the current study could be the result of spontaneous fluctuations in the BOLD signal that are intrinsically generated in the brain and not attributable to specific inputs or outputs (Fox et al., 2005, 2007; Fox and Raichle, 2007). These spontaneous fluctuations may serve to coordinate neuronal activity between distal brain regions and may be the product of changes in the power of high-frequency electrical activity such as the gamma band (Leopold et al., 2003). These changes in gamma oscillation frequency may in turn be associated with differences in GABA concentrations and their resulting influence on the balance of neural excitation and inhibition (Muthukumaraswamy et al., 2009). It is also possible that the balance of excitation and inhibition could have been influenced by glutamatergic inputs, especially given recent findings documenting an association between glutamate concentrations and reading skill (Pugh et al., 2014), as well as similar support from animal models (Che et al., 2016). Based on these findings, future investigations should target the neural mechanisms of BOLD signal variability and their links to other neurobiological indices, including neural oscillations, neurochemistry, indices of neural noise, and neuroanatomical measures (Becker et al., 2011).

Conclusions and future directions

Overall, this investigation lends support to work advocating for the added value of evaluating intrasubject variability in brain signals compared with solely evaluating mean levels of neural activity (Faisal et al., 2008; Garrett et al., 2010, 2011, 2013; Pernet et al., 2011) and highlights the importance of considering individual difference dimensions beyond subject age as contributors to, or reflections of, individual differences in neural activation variability (McIntosh et al., 2008; Grady and Garrett, 2014). Furthermore, the current findings constitute a critical first step in considering the role of adaptability in developing brain systems involved in reading and motivate future investigations concerning the mechanistic link between neural activation variability, neural noise, and reading skill. Finally, from an applied standpoint, these results beg the tantalizing question of whether trial-by-trial activation variability in reading-related brain areas could serve as a useful biomarker for clinically relevant phenotypes such as response to intervention for reading disability.

Footnotes

  • This work was supported by the Eunice Kennedy Shriver National Institute of Child Health and Human Development of the National Institutes of Health (Grant P01HD070837 to Georgia State University, Grants P01HD001994, R01HD086168, and R37HD090153 to Haskins Laboratories). F.H. was supported by the National Institutes of Health (Grants R01HD078351, R01HD067254, P50HD052120, and R01HD044073), University of California Office of the President (Grant MRP-17-454926), the Oak Foundation (Grant ORIO-16-012), and the National Science Foundation (Grant 1540854). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. We thank the children and their parents, as well as their schools and teachers, for participation in this study and Candice Goerger for her extensive role in recruitment and data collection at the Atlanta site.

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Jeffrey G. Malins, Haskins Laboratories, 300 George Street Suite 900, New Haven, CT, 06511. jeffrey.malins{at}yale.edu

References

  1. ↵
    1. Armbruster-Genç DJ,
    2. Ueltzhöffer K,
    3. Fiebach CJ
    (2016) Brain signal variability differentially affects cognitive flexibility and cognitive stability. J Neurosci 36:3978–3987. doi:10.1523/JNEUROSCI.2517-14.2016 pmid:27053205
    OpenUrlAbstract/FREE Full Text
  2. ↵
    1. Baum SH,
    2. Beauchamp MS
    (2014) Greater BOLD variability in older compared with younger adults during audiovisual speech perception. PLoS One 9:e111121. doi:10.1371/journal.pone.0111121 pmid:25337918
    OpenUrlCrossRefPubMed
  3. ↵
    1. Becker R,
    2. Reinacher M,
    3. Freyer F,
    4. Villringer A,
    5. Ritter P
    (2011) How ongoing neuronal oscillations account for evoked fMRI variability. J Neurosci 31:11016–11027. doi:10.1523/JNEUROSCI.0210-11.2011 pmid:21795550
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Bitan T,
    2. Cheon J,
    3. Lu D,
    4. Burman DD,
    5. Gitelman DR,
    6. Mesulam MM,
    7. Booth JR
    (2007) Developmental changes in activation and effective connectivity in phonological processing. Neuroimage 38:564–575. doi:10.1016/j.neuroimage.2007.07.048 pmid:17884585
    OpenUrlCrossRefPubMed
  5. ↵
    1. Bookheimer S
    (2002) Functional MRI of language: new approaches to understanding the cortical organization of semantic processing. Annu Rev Neurosci 25:151–188. doi:10.1146/annurev.neuro.25.112701.142946 pmid:12052907
    OpenUrlCrossRefPubMed
  6. ↵
    1. Branum-Martin L,
    2. Fletcher JM,
    3. Stuebing KK
    (2013) Classification and identification of reading and math disabilities: the special case of comorbidity. Journal of Learning Disabilities 46:490–499. doi:10.1177/0022219412468767 pmid:23232442
    OpenUrlCrossRefPubMed
  7. ↵
    1. Burgund ED,
    2. Kang HC,
    3. Kelly JE,
    4. Buckner RL,
    5. Snyder AZ,
    6. Petersen SE,
    7. Schlaggar BL
    (2002) The feasibility of a common stereotactic space for children and adults in fMRI studies of development. Neuroimage 17:184–200. doi:10.1006/nimg.2002.1174 pmid:12482076
    OpenUrlCrossRefPubMed
  8. ↵
    1. Centanni TM,
    2. Booker AB,
    3. Sloan AM,
    4. Chen F,
    5. Maher BJ,
    6. Carraway RS,
    7. Khodaparast N,
    8. Rennaker R,
    9. LoTurco JJ,
    10. Kilgard MP
    (2014) Knockdown of the dyslexia-associated gene Kiaa0319 impairs temporal responses to speech stimuli in rat primary auditory cortex. Cereb Cortex 24:1753–1766. doi:10.1093/cercor/bht028 pmid:23395846
    OpenUrlCrossRefPubMed
  9. ↵
    1. Centanni TM,
    2. Chen F,
    3. Booker AM,
    4. Engineer CT,
    5. Sloan AM,
    6. Rennaker RL,
    7. LoTurco JJ,
    8. Kilgard MP
    (2014) Speech sound processing deficits and training-induced neural plasticity in rats with dyslexia gene knockdown. PLoS One 9:e98439. doi:10.1371/journal.pone.0098439 pmid:24871331
    OpenUrlCrossRefPubMed
  10. ↵
    1. Che A,
    2. Truong DT,
    3. Fitch RH,
    4. LoTurco JJ
    (2016) Mutation of the dyslexia-associated gene Dcdc2 enhances glutamatergic synaptic transmission between layer 4 neurons in mouse neocortex. Cereb Cortex 26:3705–3718. doi:10.1093/cercor/bhv168 pmid:26250775
    OpenUrlCrossRefPubMed
  11. ↵
    1. Cox RW
    (1996) AFNI: software for analysis and visualization of functional magnetic resonance neuroimages. Computers and Biomedical Research 29:162–173. doi:10.1006/cbmr.1996.0014 pmid:8812068
    OpenUrlCrossRefPubMed
  12. ↵
    1. Dunn LM,
    2. Dunn LM
    (1997) PPVT-III: Peabody Picture Vocabulary Test. Bloomington, MN: Pearson Assessments.
  13. ↵
    1. Dunn LM,
    2. Dunn DM
    (2007) PPVT-4: Peabody Picture Vocabulary Test. Bloomington, MN: Pearson Assessments.
  14. ↵
    1. Efron B
    (1981) Nonparametric estimates of standard error: the jackknife, the bootstrap and other methods. Biometrika 68:589–599. doi:10.1093/biomet/68.3.589
    OpenUrlCrossRef
  15. ↵
    1. Faisal AA,
    2. Selen LP,
    3. Wolpert DM
    (2008) Noise in the nervous system. Nat Rev Neurosci 9:292–303. doi:10.1038/nrn2258 pmid:18319728
    OpenUrlCrossRefPubMed
  16. ↵
    1. Fox MD,
    2. Raichle ME
    (2007) Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging. Nat Rev Neurosci 8:700–711. doi:10.1038/nrn2201 pmid:17704812
    OpenUrlCrossRefPubMed
  17. ↵
    1. Fox MD,
    2. Snyder AZ,
    3. Vincent JL,
    4. Corbetta M,
    5. Van Essen DC,
    6. Raichle ME
    (2005) The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc Natl Acad Sci U S A 102:9673–9678. doi:10.1073/pnas.0504136102 pmid:15976020
    OpenUrlAbstract/FREE Full Text
  18. ↵
    1. Fox MD,
    2. Snyder AZ,
    3. Vincent JL,
    4. Raichle ME
    (2007) Intrinsic fluctuations within cortical systems account for intertrial variability in human behavior. Neuron 56:171–184. doi:10.1016/j.neuron.2007.08.023 pmid:17920023
    OpenUrlCrossRefPubMed
  19. ↵
    1. Frost SJ,
    2. Landi N,
    3. Mencl WE,
    4. Sandak R,
    5. Fulbright RK,
    6. Tejada ET,
    7. Jacobsen L,
    8. Grigorenko EL,
    9. Constable RT,
    10. Pugh KR
    (2009) Phonological awareness predicts activation patterns for print and speech. Annals of Dyslexia 59:78–97. doi:10.1007/s11881-009-0024-y pmid:19306061
    OpenUrlCrossRefPubMed
  20. ↵
    1. Galaburda AM,
    2. Sherman GF,
    3. Rosen GD,
    4. Aboitiz F,
    5. Geschwind N
    (1985) Developmental dyslexia: four consecutive cases with cortical anomalies. Annals of Neurology 18:222–233. doi:10.1002/ana.410180210
    OpenUrlCrossRefPubMed
  21. ↵
    1. Garrett DD,
    2. Kovacevic N,
    3. McIntosh AR,
    4. Grady CL
    (2010) Blood oxygen level-dependent signal variability is more than just noise. J Neurosci 30:4914–4921. doi:10.1523/JNEUROSCI.5166-09.2010 pmid:20371811
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Garrett DD,
    2. Kovacevic N,
    3. McIntosh AR,
    4. Grady CL
    (2011) The importance of being variable. J Neurosci 31:4496–4503. doi:10.1523/JNEUROSCI.5641-10.2011 pmid:21430150
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Garrett DD,
    2. Samanez-Larkin GR,
    3. MacDonald SW,
    4. Lindenberger U,
    5. McIntosh AR,
    6. Grady CL
    (2013) Moment-to-moment brain signal variability: a next frontier in human brain mapping? Neurosci Biobehav Rev 37:610–624. doi:10.1016/j.neubiorev.2013.02.015 pmid:23458776
    OpenUrlCrossRefPubMed
  24. ↵
    1. Grady CL,
    2. Garrett DD
    (2014) Understanding variability in the BOLD signal and why it matters for aging. Brain Imaging and Behavior 8:274–283. doi:10.1007/s11682-013-9253-0 pmid:24008589
    OpenUrlCrossRefPubMed
  25. ↵
    1. Grinband J,
    2. Wager TD,
    3. Lindquist M,
    4. Ferrera VP,
    5. Hirsch J
    (2008) Detection of time-varying signals in event-related fMRI designs. Neuroimage 43:509–520. doi:10.1016/j.neuroimage.2008.07.065 pmid:18775784
    OpenUrlCrossRefPubMed
  26. ↵
    1. Haigh SM,
    2. Gupta A,
    3. Barb SM,
    4. Glass SAF,
    5. Minshew NJ,
    6. Dinstein I,
    7. Heeger DJ,
    8. Eack SM,
    9. Behrmann M
    (2016) Differential sensory fMRI signatures in autism and schizophrenia: analysis of amplitude and trial-to-trial variability. Schizophrenia Research 175:12–19. doi:10.1016/j.schres.2016.03.036 pmid:27083780
    OpenUrlCrossRefPubMed
  27. ↵
    1. Hancock R,
    2. Pugh KR,
    3. Hoeft F
    (2017) Neural noise hypothesis of developmental dyslexia. Trends Cogn Sci 21:434–448. doi:10.1016/j.tics.2017.03.008 pmid:28400089
    OpenUrlCrossRefPubMed
  28. ↵
    1. He BJ,
    2. Zempel JM
    (2013) Average is optimal: an inverted-U relationship between trial-to-trial brain activity and behavioral performance. PLoS Computational Biology 9:e1003348.
    OpenUrl
  29. ↵
    1. Ho DE,
    2. Imai K,
    3. King G,
    4. Stuart EA
    (2011) MatchIt: nonparametric preprocessing for parametric causal inference. Journal of Statistical Software 42:1–28.
    OpenUrl
  30. ↵
    1. Hornickel J,
    2. Kraus N
    (2013) Unstable representation of sound: a biological marker of dyslexia. J Neurosci 33:3500–3504. doi:10.1523/JNEUROSCI.4205-12.2013 pmid:23426677
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Jasińska KK,
    2. Molfese PJ,
    3. Kornilov SA,
    4. Mencl WE,
    5. Frost SJ,
    6. Lee M,
    7. Pugh KR,
    8. Grigorenko EL,
    9. Landi N
    (2016) The BDNF Val66Met polymorphism influences reading ability and patterns of neural activation in children. PLoS One 11:e0157449. doi:10.1371/journal.pone.0157449 pmid:27551971
    OpenUrlCrossRefPubMed
  32. ↵
    1. Landi N,
    2. Frost SJ,
    3. Mencl WE,
    4. Preston JL,
    5. Jacobsen LK,
    6. Lee M,
    7. Yrigollen C,
    8. Pugh KR,
    9. Grigorenko EL
    (2013) The COMT Val/Met polymorphism is associated with reading-related skills and consistent patterns of functional neural activation. Dev Sci 16:13–23. doi:10.1111/j.1467-7687.2012.01180.x pmid:23278923
    OpenUrlCrossRefPubMed
  33. ↵
    1. Leopold DA,
    2. Murayama Y,
    3. Logothetis NK
    (2003) Very slow activity fluctuations in monkey visual cortex: implications for functional brain imaging. Cereb Cortex 13:422–433. doi:10.1093/cercor/13.4.422 pmid:12631571
    OpenUrlCrossRefPubMed
  34. ↵
    1. Lund TE,
    2. Nørgaard MD,
    3. Rostrup E,
    4. Rowe JB,
    5. Paulson OB
    (2005) Motion or activity: their role in intra- and inter-subject variation in fMRI. Neuroimage 26:960–964. doi:10.1016/j.neuroimage.2005.02.021 pmid:15955506
    OpenUrlCrossRefPubMed
  35. ↵
    1. Martin A,
    2. Schurz M,
    3. Kronbichler M,
    4. Richlan F
    (2015) Reading in the brain of children and adults: a meta-analysis of 40 functional magnetic resonance imaging studies. Hum Brain Mapp 36:1963–1981. doi:10.1002/hbm.22749 pmid:25628041
    OpenUrlCrossRefPubMed
  36. ↵
    1. McIntosh AR,
    2. Kovacevic N,
    3. Itier RJ
    (2008) Increased brain signal variability accompanies lower behavioral variability in development. PLoS Comput Biol 4:e1000106. doi:10.1371/journal.pcbi.1000106 pmid:18604265
    OpenUrlCrossRefPubMed
  37. ↵
    1. Mumford JA,
    2. Turner BO,
    3. Ashby FG,
    4. Poldrack RA
    (2012) Deconvolving BOLD activation in event-related designs for multivoxel pattern classification analyses. Neuroimage 59:2636–2643. doi:10.1016/j.neuroimage.2011.08.076 pmid:21924359
    OpenUrlCrossRefPubMed
  38. ↵
    1. Muthukumaraswamy SD,
    2. Edden RA,
    3. Jones DK,
    4. Swettenham JB,
    5. Singh KD
    (2009) Resting GABA concentration predicts peak gamma frequency and fMRI amplitude in response to visual stimulation in humans. Proc Natl Acad Sci U S A 106:8356–8361. doi:10.1073/pnas.0900728106 pmid:19416820
    OpenUrlAbstract/FREE Full Text
  39. ↵
    1. Norton ES,
    2. Beach SD,
    3. Gabrieli JD
    (2015) Neurobiology of dyslexia. Curr Opin Neurobiol 30:73–78. doi:10.1016/j.conb.2014.09.007 pmid:25290881
    OpenUrlCrossRefPubMed
  40. ↵
    1. Pennington BF,
    2. Santerre-Lemmon L,
    3. Rosenberg J,
    4. MacDonald B,
    5. Boada R,
    6. Friend A,
    7. Leopold DR,
    8. Samuelsson S,
    9. Byrne B,
    10. Willcutt EG,
    11. Olson RK
    (2012) Individual prediction of dyslexia by single versus multiple deficit models. Journal of Abnormal Psychology 121:212–224. doi:10.1037/a0025823 pmid:22022952
    OpenUrlCrossRefPubMed
  41. ↵
    1. Pernet CR,
    2. Sajda P,
    3. Rousselet GA
    (2011) Single-trial analyses: why bother? Front Psychol 2:322. pmid:22073038
    OpenUrlCrossRefPubMed
  42. ↵
    1. Poldrack RA,
    2. Wagner AD,
    3. Prull MW,
    4. Desmond JE,
    5. Glover GH,
    6. Gabrieli JD
    (1999) Functional specialization for semantic and phonological processing in the left inferior prefrontal cortex. Neuroimage 10:15–35.
    OpenUrlCrossRefPubMed
  43. ↵
    1. Power JD,
    2. Barnes KA,
    3. Snyder AZ,
    4. Schlaggar BL,
    5. Petersen SE
    (2012) Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage 59:2142–2154. doi:10.1016/j.neuroimage.2011.10.018 pmid:22019881
    OpenUrlCrossRefPubMed
  44. ↵
    1. Preston JL,
    2. Molfese PJ,
    3. Frost SJ,
    4. Mencl WE,
    5. Fulbright RK,
    6. Hoeft F,
    7. Pugh KR
    (2016) Print-speech convergence predicts future reading outcomes in early readers. Psychol Sci 1–1.
  45. ↵
    1. Pugh KR,
    2. Landi N,
    3. Preston JL,
    4. Mencl WE,
    5. Austin AC,
    6. Sibley D,
    7. Fulbright RK,
    8. Seidenberg MS,
    9. Grigorenko EL,
    10. Constable RT,
    11. Molfese P,
    12. Frost SJ
    (2013) The relationship between phonological and auditory processing and brain organization in beginning readers. Brain Lang 125:173–183. doi:10.1016/j.bandl.2012.04.004 pmid:22572517
    OpenUrlCrossRefPubMed
  46. ↵
    1. Pugh KR,
    2. Frost SJ,
    3. Rothman DL,
    4. Hoeft F,
    5. Del Tufo SN,
    6. Mason GF,
    7. Molfese PJ,
    8. Mencl WE,
    9. Grigorenko EL,
    10. Landi N,
    11. Preston JL,
    12. Jacobsen L,
    13. Seidenberg MS,
    14. Fulbright RK
    (2014) Glutamate and choline levels predict individual differences in reading ability in emergent readers. J Neurosci 34:4082–4089. doi:10.1523/JNEUROSCI.3907-13.2014 pmid:24623786
    OpenUrlAbstract/FREE Full Text
  47. ↵
    1. Rissman J,
    2. Gazzaley A,
    3. D'Esposito M
    (2004) Measuring functional connectivity during distinct stages of a cognitive task. Neuroimage 23:752–763. doi:10.1016/j.neuroimage.2004.06.035 pmid:15488425
    OpenUrlCrossRefPubMed
  48. ↵
    1. Schlaggar BL,
    2. Brown TT,
    3. Lugar HM,
    4. Visscher KM,
    5. Miezin FM,
    6. Petersen SE
    (2002) Functional neuroanatomical differences between adults and school-age children in the processing of single words. Science 296:1476–1479. doi:10.1126/science.1069464
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Skudlarski P,
    2. Constable RT,
    3. Gore JC
    (1999) ROC analysis of statistical methods used in functional MRI: individual subjects. Neuroimage 9:311–329. doi:10.1006/nimg.1999.0402 pmid:10075901
    OpenUrlCrossRefPubMed
  50. ↵
    1. Tibshirani R,
    2. Leisch F
    (2015) Bootstrap: functions for the book “an introduction to the bootstrap”. R package version 2015.2.
  51. ↵
    1. Turkeltaub PE,
    2. Gareau L,
    3. Flowers DL,
    4. Zeffiro TA,
    5. Eden GF
    (2003) Development of neural mechanisms for reading. Nat Neurosci 6:767–773. doi:10.1038/nn1065 pmid:12754516
    OpenUrlCrossRefPubMed
  52. ↵
    1. Venables WN,
    2. Ripley BD
    (2002) Modern applied statistics with S, Ed 4. New York, NY: Springer.
    1. Wagner RK,
    2. Torgesen JK,
    3. Rashotte CA
    (1999) CTOPP: Comprehensive Test of Phonological Processing. Austin, TX: Pro-Ed.
  53. ↵
    1. Wagner RK,
    2. Torgesen JK,
    3. Rashotte CA,
    4. Pearson NA
    (2013) Comprehensive Test of Phonological Processing: CTOPP2. Austin, TX: Pro-Ed.
  54. ↵
    1. Wechsler D
    (1999) Wechsler Abbreviated Scale of Intelligence (WASI-I). San Antonio, TX: The Psychological Corporation.
  55. ↵
    1. Wechsler D,
    2. Hsiao-pin C
    (2011) WASI-II: Wechsler Abbreviated Scale of Intelligence-Second Edition (WASI-II). San Antonio, TX: Pearson Assessments.
  56. ↵
    1. Woodcock RW,
    2. McGrew KS,
    3. Mather N
    (2001) Woodcock-Johnson III Tests of Achievement. Itasca, IL: Riverside Publishing.
  57. ↵
    1. Yarkoni T,
    2. Barch DM,
    3. Gray JR,
    4. Conturo TE,
    5. Braver TS
    (2009) BOLD correlates of trial-by-trial reaction time variability in gray and white matter: a multi-study fMRI analysis. PLoS One 4:e4257. doi:10.1371/journal.pone.0004257 pmid:19165335
    OpenUrlCrossRefPubMed
Back to top

In this issue

The Journal of Neuroscience: 38 (12)
Journal of Neuroscience
Vol. 38, Issue 12
21 Mar 2018
  • 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.
Individual Differences in Reading Skill Are Related to Trial-by-Trial Neural Activation Variability in the Reading Network
(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
Individual Differences in Reading Skill Are Related to Trial-by-Trial Neural Activation Variability in the Reading Network
Jeffrey G. Malins, Kenneth R. Pugh, Bonnie Buis, Stephen J. Frost, Fumiko Hoeft, Nicole Landi, W. Einar Mencl, Anish Kurian, Ryan Staples, Peter J. Molfese, Rose Sevcik, Robin Morris
Journal of Neuroscience 21 March 2018, 38 (12) 2981-2989; DOI: 10.1523/JNEUROSCI.0907-17.2018

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
Individual Differences in Reading Skill Are Related to Trial-by-Trial Neural Activation Variability in the Reading Network
Jeffrey G. Malins, Kenneth R. Pugh, Bonnie Buis, Stephen J. Frost, Fumiko Hoeft, Nicole Landi, W. Einar Mencl, Anish Kurian, Ryan Staples, Peter J. Molfese, Rose Sevcik, Robin Morris
Journal of Neuroscience 21 March 2018, 38 (12) 2981-2989; DOI: 10.1523/JNEUROSCI.0907-17.2018
Reddit logo Twitter logo Facebook 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

  • β series
  • BOLD variability
  • event-related fMRI
  • individual differences
  • reading disability
  • trial-by-trial variability

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

  • The human centromedial amygdala contributes to negative prediction error signaling during appetitive and aversive Pavlovian gustatory learning
  • Perceptual Difficulty Regulates Attentional Gain Modulations in Human Visual Cortex
  • Concurrent- and after-effects of medial temporal lobe stimulation on directed information flow to and from prefrontal and parietal cortices during memory formation
Show more Research Articles

Behavioral/Cognitive

  • The human centromedial amygdala contributes to negative prediction error signaling during appetitive and aversive Pavlovian gustatory learning
  • Perceptual Difficulty Regulates Attentional Gain Modulations in Human Visual Cortex
  • Concurrent- and after-effects of medial temporal lobe stimulation on directed information flow to and from prefrontal and parietal cortices during memory formation
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.