The timing and strength of regional brain activation associated with word recognition in children with reading difficulties.

The study investigates the relative degree and timing of cortical activation across parietal, temporal, and frontal regions during performance of a continuous visual-word recognition task in children who experience reading difficulties (N = 44, RD) and typical readers (N = 40, NI). Minimum norm estimates of regional neurophysiological activity were obtained from magnetoencephalographic recordings. Children with RD showed bilaterally reduced neurophysiological activity in the superior and middle temporal gyri, and increased activity in rostral middle frontal and ventral occipitotemporal cortices, bilaterally. The temporal profile of activity in the RD group, featured near-simultaneous activity peaks in temporal, inferior parietal, and prefrontal regions, in contrast to a clear temporal progression of activity among these areas in the NI group. These results replicate and extend previous MEG and fMRI results demonstrating atypical, latency-dependent attributes of the brain circuit involved in word reading in children with reading difficulties.

Introduction

Among several explanations of developmental reading difficulties (RD), the presence of a core deficit in the ability to mentally represent and manipulate the sounds of the language (phonological awareness) appears to account for a significant amount of variance in word-level reading scores, especially during the earlier stages of reading acquisition. The ultimate goal of reading, however, is automatic recognition of words and comprehension of their meaning (Gough and Tumner, 1986; Share and Stanovich, 1995; Berninger et al., 2006), which involves a variety of other skills such as vocabulary and working memory (Kaplan et al., 1998; O'Shaughnessy and Swanson, 1998; Willcutt et al., 2001; Storch and Whitehurst, 2002; Ouellette and Sénéchal, 2008). Although in each individual with RD several alternative or complementary factors may account for such difficulties, the observation that RD is often found in children with intact or even superior intelligence and oral (non-phonological) language skills indicates that RD is associated with functional deficits specific to the brain circuitry that supports reading.

With the advent of functional neuroimaging methods, it has become possible to detect differences in the cerebral activation patterns associated with processes such as phonological decoding and word recognition between RD and non-impaired readers (NI). In particular, evidence from positron emission tomography (PET) and functional magnetic resonance imaging (fMRI) studies suggested that performance deficits in children with RD is paralleled by hypoactivation in several distinct regions critical for reading ability, including the superior (BA 22) and middle (BA 21) temporal gyri, the inferior parietal lobule [supramarginal and angular gyri (BA 40/39)] and the ventral portion of occipitotemporal cortex (mainly in the left hemisphere; BA 37; Temple et al., 2001; Shaywitz et al., 2002; Cao et al., 2006; Hoeft et al., 2007; for a recent review and meta-analysis see Maisog et al., 2008). An associated finding of either increased or decreased activation in prefrontal (mainly left inferior frontal, IFG) regions has also been reported in these studies.

In addition to PET and fMRI, magnetoencephalography, also known as magnetic source imaging (MSI), is an alternative neuroimaging technique that offers information both on spatial localization and on the timing (latency) of neurophysiological processes. The capacity of MSI to generate spatial maps of language-specific cortex has been validated through direct comparisons with invasive cortical mapping procedures in large-scale clinical populations (e.g., Papanicolaou et al., 2004; Kamada et al., 2007).

Over the course of several studies we have used MSI to differentiate between children with RD and NI readers during execution of both explicit phonological decoding (Simos et al., 2000b, 2002a, 2007), and continuous visual-word recognition tasks (silent real word reading; Simos et al., 2000a). In the latter study the magnitude of neurophysiological activity in left hemisphere posterior temporal (superior and middle temporal gyri) and inferior parietal regions (supramarginal gyrus) was significantly reduced among children with severe RD as compared to children with average to above average reading achievement. Children with RD also showed increased activity in homologous right hemisphere regions. Another reported feature of aberrant profiles in children with RD on activation tasks involving pseudoword reading involves the relative timing of regional activation, such that the early occipitotemporal (visual association cortex) activity was immediately followed by early engagement of prefrontal activity and only later by, activation of temporoparietal cortex.

Earlier MEG studies of RD suffered from a number of potential limitations. One potential limitation was that spatiotemporal activation profiles were obtained by iterative application of the single equivalent dipole model (Sarvas, 1987). While this method has been routinely applied to estimate source locations for sensory- and language-related magnetic fields, and has been subjected to external validation against invasive diagnostic techniques (Szymanski et al., 2001; Papanicolaou et al., 2004; Kamada et al., 2007), a single-dipole model may not fully represent the observed field distributions, a particular problem for modeling brain activation in complex cognitive tasks involving multiple brain areas. A means of overcoming this potential problem is the application of distributed source models and, in particular, the technique known as minimum-norm estimate (MNE; Hämäläinen and Ilmoniemi, 1994; Moran and Tepley, 2000) chosen in the present study. Furthermore, previous MEG studies (as well as some PET and fMRI studies) were limited by their reliance on participant samples that lacked adequate diversity in terms of socioeconomic status and IQ. Finally, another potential source of experimental confounds in early MEG studies of word recognition was that brain activation profiles were obtained with tasks that placed heavy demands on both word recognition and memory processes, despite evidence that task difficulty may significantly impact the degree of activation in key cortical areas (Cao et al., 2006; Pugh et al., 2008).

The objectives of the current investigation were to (a) replicate and extend previous MEG results on the key anatomical and temporal features of the brain activation profiles associated with silent word reading; and (b) assess the relative importance of time-dependent engagement of posterior temporal, inferior parietal, and frontal regions for word recognition in children who are either typical or struggling readers. In order to address some of the limitations of previous studies outlined above: (i) a larger, more representative sample of the general school population was selected; (ii) brain activation profiles were reconstructed in the context of much easier version of the continuous word recognition task used in our previous study (Simos et al., 2000a); and (iii) task-related neurophysiological activity was estimated by application of the MNE technique, mentioned above.

Materials and Methods

Participants

The primary target group consisted of 44 children with reading difficulties (RD group) as indicated by scores at or below the 16th percentile (standard score of 85) on the basic reading composite (average of word attack and letter–word identification subtest score of the Woodcock–Johnson III; Woodcock et al., 2001). A second group of 40 children, who had never experienced difficulties in reading (NI group), served as the comparison group, obtaining standard scores ≥ 88 on the reading composite index. All 84 participants had full scale IQ scores > 80 on the Wechsler (1999) abbreviated scale of intelligence. None of the participants had a history or current diagnosis of ADHD. In addition, all children were screened for the presence of attention and/or hyperactivity/impulsivity problems using one of two parental-report symptom checklists (child behavior checklist, CBCL; Achenbach, 1991; SNAP-IV, Swanson, 1992). All children had T scores on the attention problems CBCL scale below 55 or a score lower than 1.6 on the inattention and hyperactivity–impulsivity scales of the SNAP-IV indicating low risk for ADHD (Chen et al., 1994). As shown in Table 1, the two reading groups were comparable in age, PIQ, gender, ethnicity, and handedness. The differences in VIQ (and FSIQ) are not surprising and found in most behavioral and imaging studies; nonetheless, the average VIQ (and FSIQ) of the group with RD was clearly in the average range.

Table 1

Demographic data, educational history, and performance on standardized tests for students with reading difficulties (RD) and non-reading impaired students (NI; range, mean and SD in parentheses).

	NI	RD
Gender (boys/girls)	25/15	27/17
Age (months)	7.0–15.0 (11.0 ± 2.3)	7.7–14.3 (11.4 ± 1.7)
Ethnicity (Caucasian, H, A-AM2)	16/13/11	18/14/12
Handedness (R/L)	35/5	37/7
Number of students retained in school	–	10
Number of students receiving special education services	–	29
VIQ*	86–147 (108 ± 17)	80–128 (97 ± 12)
PIQ	80–117 (100 ± 10)	80–129 (97 ± 13)
FSIQ*	81–130 (104 ± 12)	83–119 (97 ± 10)
WJ-III Reading Composite†	88–126 (107 ± 10)	53–85 (78 ± 8)
WJ-III Word Attack†	87–131 (107 ± 12)	59–91 (80 ± 7)
WJ-III Letter-Word Identification†	87–126 (107 ± 11)	35–88 (75 ± 12)
WJ-III Spelling†	83–136 (112 ± 14)	20–91 (73 ± 14)

Pairwise group differences: *p < 0.001, .

The potential mediating effects of age on reading-group differences were assessed by treating age as a dichotomous, between-subjects variable. To ensure adequate cell size, each reading ability group was broken down to two subgroups: younger participants (21 NI and 24 RD) aged 84–130 months, and older participants (19 NI and 20 RD) aged 131–180 months. All participants were recruited through ads placed in area schools and presentations of the study by the researchers in parent–teacher association meetings. Parent who expressed an interest in participating in the study were subsequently contacted to schedule a meeting for administering the screening instruments. Following approval of the study by the university of Texas-Health Science Center IRB, all parents were asked to sign an informed consent form while underage participants signed an assent form.

Procedures

Tasks

Each participant was tested on a printed word recognition task, adapted from Simos et al. (2000a). Changes in the selection of the stimuli and the number of targets were made to reduce overall task difficulty and memory load. There were only five target stimuli which were repeated in a different random order, mixed with a different set of 30 distractors (non-repeating words) in each of four blocks of stimuli (for a total of 140 stimuli). Target words (age, enjoy, hope, jar, and road) included four monosyllabic and one disyllabic word, and had a mean frequency in the Zeno et al. (1995) G6–7 corpus of 158 occurrences per million (range: 32–194 occurrences). A slightly higher proportion of distractors were disyllabic (40%) and the remaining monosyllabic, with a mean frequency of occurrence of 150 words per million in the same corpus (range: 18–820). Stimuli were presented for 1 s, one at a time (with a randomly varied interstimulus interval of 3–4 s), through a Sony LCD projector (Model VPL-PX21) on a back-projection screen located approximately 60 cm in front of the participant, and subtended 1.0–3.0° and 0.5° of horizontal and vertical visual angle, respectively. Participants were presented the target words prior to the scan and instructed to read them aloud (to ensure that they could) and to “try to remember the words.” During the scan they were presented with the target stimuli mixed with the distractors and asked to raise their index finger “as soon as they saw one of the words they had read before the scan.” Responding hand was counterbalanced across blocks of trials. Responses were manually recorded by the experimenter.

Imaging procedures

Magnetoencephalography recordings were obtained with a whole-head neuromagnetometer array (4-D Neuroimaging, Magnes WH3600) equipped with 248 first-order axial gradiometer coils and housed in a magnetically shielded chamber. The magnetic flux measurements were digitized at 250 Hz, filtered off-line with a bandpass filter between 0.1 and 20 Hz and subjected to a noise reduction algorithm that is part of the 4D-Neuroimaging software. The single-trial event-related field segments (ERFs) in response to 85–115 stimulus presentations (out of a total of 140 recorded), were averaged after excluding those containing eye movement or other myogenic or mechanical artifacts. Only correct responses were analyzed for targets and non-targets.

To identify the intracranial origin of ERFs the magnetic flux distribution recorded simultaneously over the entire head surface at successive time points (4 ms apart) was analyzed using a minimum norm model to obtain estimates of the time-varying strength of intracranial currents (MNE Software, v. 2.5, http://www.nmr.mgh.harvard.edu/martinos/userInfo/data/sofMNE.php). This method permits detection of simultaneous magnetic sources distributed along the entire cortical surface. The model assumes a continuous distribution of current along the cortical surface which has some minimum norm (Hämäläinen and Ilmoniemi, 1994). Estimated current sources were anatomically constrained by an MRI-derived surface model of each participant's brain (T1-weighted: TR 13.6 ms; TE 4.8 ms; recording matrix 256 × 256 pixels, one excitation, 240 mm field of view, and 1.4-mm slice thickness). This surface model was generated by a fully automated cortical surface reconstruction procedure using FreeSurfer software (Dale et al., 1999) for producing a detailed geometric description (regular tessellation of the cortical surface consisting of equilateral triangles known as vertices) of the gray-white matter boundary of the neocortical mantle and the mesial temporal lobe. Each hemisphere consisted of approximately 150,000 vertices (depending on each subject's cortical surface area). For estimating current sources, the MNE software requires the Freesurfer-derived cortical surface reconstruction for defining the solution source space. A grid-spacing of 7 mm was used to construct icosahedrons to decimate the number of vertices from 150,000 to approximately 3,000 per hemisphere. Additionally, the MNE software was used to construct a single compartment boundary element model using triangular tessellations to compute the forward model for a current dipole located at each vertex perpendicular to the cortical surface. The inverse solution was subsequently reduced to obtaining an estimate of the scalar distribution of dipole strength across current sources within orientation-specific cortical patches of vertices (Dale and Sereno, 1993). Co-registration of each MEG dataset with its corresponding MRI dataset was performed using an automated co-registration routine within MNE which aligns digitization points in the MEG headshape file with the fiducial points demarcated on the outer skin surface reconstruction of the MRI.

Based on previous reports on the cortical areas that constitute major components of the brain mechanism for reading in children, the following ROIs were examined (in both hemispheres): superior (BA22) and middle temporal gyri (BA 21), supramarginal gyrus (BA 40), angular gyrus (BA 39), inferior (BA 44/45), rostral middle frontal cortex (BA 46), fusiform gyrus (BA 37), lateral occipitotemporal cortex (BA 19), and mesial temporal cortex (entorhinal and parahippocampal cortices; MTL). The anatomical extent of ROIs was determined automatically using the cortical map provided by Freesurfer (Dale et al., 1999). The program outputs a current estimate measure for each voxel at each 4-ms time point. This measure is then used to compute the two dependent measures used in the analyses outlined in the following paragraphs. First, the average current across all voxels defining each of the ROIs listed above and across all of the 4-ms time points comprising 14 successive 50-ms time bins (100–150, 150–200, etc up to 800 ms). Second, the latency (in ms after stimulus onset) when averaged current reached peak amplitude within the entire recording epoch. This approach aimed at (a) determining if reading-group differences in the degree of regional activity were age-, time- and hemisphere-dependent (i.e., more systematic for particular time windows and restricted to one hemisphere), (b) establishing the temporal progression of regional activity in each group of participants, and (c) identifying ROIs and time bins where activity correlated with reading ability (measured through standardized achievement tests and during performance of the activation task in the MEG scanner).

Analyses

Preliminary analyses indicated that group differences (RD > NI) were present as early as between 0 and 50 ms after stimulus onset [as indicated by a reading group by ROI interaction, F(8,672) = 59.35, p < 0.0001] in three ROIs: BA46, BA44/45, and BA37 (in all cases p < 0.0001 for the simple main effects of reading group). In order to correct for potential group differences in peristimulus activity, current estimates in all time bins between 100 and 800 ms were adjusted by subtracting the average current during 0–50 ms (to ensure that no stimulus-evoked activity was included) in the corresponding ROI for each participant. The average, adjusted current for each 50 ms time bin and each ROI was then submitted to an ANOVA with ROI (9), Hemisphere (2), and Time bin (14) as the within subjects variables. Reading group (2) and age group (2) served as the between subjects variables. Significant interactions involving reading group were further evaluated by examining four-way (e.g., time by ROI by hemisphere by reading group) or three way interactions (e.g., Time by ROI by reading group) which, if significant, were explored by testing two-way interactions (time by reading group). Significant hemisphere asymmetries in the degree of adjusted regional activity will be mentioned (in the form of Hemisphere or Time by Hemisphere effects), even in the absence of reading group-specific hemisphere effects, in order to facilitate comparisons between the present and findings from previous studies using similar tasks. All ANOVA results were evaluated using the Huynh–Feldt method as a precaution against inhomogeneity of variance problems.

The second set of analyses established the outline of the spatiotemporal profiles of activity for participants in each reading group – essentially the temporal progression of regional activity for each group. The dependent variable in these analyses was the median time point of the 50-ms time bin when average, adjusted MNE current reached peak value. ROIs from both hemispheres were ranked according to peak latency and a series of dependent-sample t-tests were computed between ROIs (earliest peaking ROI with each subsequent ROI, second earliest ROI with each subsequent ROI, etc), in order to test the hypothesis that regional onset latency differences were statistically significant. In order to control for Type I error for multiple comparisons, all tests were evaluated at α = 0.001. This procedure was performed separately for each reading group. Results were used to identify sets of ROIs that peaked at approximately the same latency after stimulus onset. When activity did not peak at exactly the same time, statistical tests indicated that peak latencies were not systematically different across regions that formed each set. Activity in the majority of ROIs belonging to a particular set of ROIs peaked significantly earlier than activity in the majority of ROIs in the set of subsequently active regions. Across reading groups two or three sets of successively active ROIs were identified.

The final set of analyses consisted of partial correlations between reading measures and peak latency or average, adjusted current in ROIs (and time bins) where significant reading group differences were revealed by previous analyses. Age and VIQ were entered in these analyses as covariates. Notably, the metric of the degree of regional activation provided by MNE is void of the deviation-from-normality problems permitting the application of linear regression algorithms. Count data are generally more susceptive to deviations from normality and alternatives to linear models are more appropriate (such as the general estimable function procedure). In the case of ECD count data the range of values in a particular region is often very large (e.g., between 0 and 70) with greatly uneven distributions. This problem is especially apparent for regression models.

Results

In-scanner task performance

Percent hits (percent correct identification of targets) was slightly higher for the NI (80.3 ± 14%) compared to the RD group (70.5 ± 17%), F(1,82) = 3.10, p < 0.06). A similar trend was found for false alarms (14.2 ± 5% and 17.4 ± 4%, respectively, p < 0.08). Although, the main effect of Age Group did not reach significance on either measure (p > 0.2), there was a reading by age group interaction on the percentage of hits, F(1,82) = 4.30, p < 0.03, suggesting that age-related differences in response accuracy were larger among RD (younger: 60.4 ± 14%, older: 78.3 ± 16%) than NI students (younger: 77.6 ± 9%, older: 83.6 ± 10%).

MEG data

Spatiotemporal profiles of activity

Sample snapshots of the activation profiles of two participants, one from the NI group and a second from the RD group in Figure 1, demonstrate notable group differences in the relative timing of regional activity. Activity in both groups peaked initially in ventral and lateral occipitotemporal (LOC, ROIs bilaterally (NI readers: between 175 ± 54 and 188 ± 46 ms, RD readers: between 168 ± 64 and 185 ± 74 ms). For the NI group, activity in the left MTG, STG, and SMG peaked after a significant temporal delay (at 217 ± 66, 223 ± 70, and 248 ± 62 ms, respectively). Activity in left hemisphere inferior and middle frontal ROIs peaked significantly later (at 284 ± 82 and 289 ± 79 ms, respectively) than activity in left temporoparietal ROIs. The temporal progression of activity for the RD group was distinct from that of the NI group in that (a) on average, activity in occipitotemporal ROIs and the left MTG (at 206 ± 69 ms) peaked nearly at the same time and (b) although activity in the left STG and SMG peaked significantly later than activity in occipitotemporal ROIs (at 214 ± 89 and 259 ± 107 ms, respectively, as for the NI group), activity peaks in left inferior and middle frontal cortices occurred near simultaneously with temporoparietal ROIs (at 233 ± 73 and 249 ± 87 ms, respectively).

Figure 1

Brain activation map snapshots computed at different time points during the visual-word recognition task from two representative participants [a student from the NI group (left pair of columns) and a student with reading difficulties (right pair of columns)], emphasizing activity in the lateral occipitotemporal (LOC), ventral occipitotemporal (VOC), superior/middle temporal (STG/MTG), inferior frontal (IFG), and middle frontal regions of the left (LH) and right (RH) hemispheres. The latency, post-stimulus onset that each map was obtained is noted next to the ROI label. Values represent “raw” estimated current (without subtracting the peristimulus current for each ROI).

Reading group differences in degree of activity

The omnibus ANOVA for the degree of activity data revealed a significant Time by ROI by reading group interaction, F(104,8320) = 2.51, p < 0.0001, which was further explored via two-way ANOVAs (time by reading group) separately for each ROI [evaluated at alpha = 0.05/9 (ROIs) = 0.006 to correct for family-wise Type I error]. There were no interactions of reading with age group.

Group differences as a function of latency are displayed in Figure 2. Significant time by reading group interactions were found for activity (adjusted for peristimulus estimated current) in MTG, F(13,1066) = 4.32, p < 0.0001, and STG, F(13,1066) = 4.99, p < 0.0001. Follow-up one-way ANOVAs performed at each successive time bin across hemispheres [evaluated at alpha = 0.05/14 (time bins) = 0.004] revealed greater degree of activity for NI than RD participants during the late portion of the source density waveform [between 200 and 650 ms: F(1,82) = 30.42, p < 0.0001 and F(1,82) = 18.87, p < 0.0001, respectively]. Group effects in the opposite direction (RD > NI) were found in the fusiform and middle frontal ROIs as indicated by time by reading group interactions, F(13,1066) = 4.32, p < 0.0001 and F(13,1066) = 4.00, p < 0.0001, respectively. Follow-up one-way ANOVAs revealed that group differences were significant in the early portion of the source density waveform (150–450 ms) in the fusiform, F(1,82) = 20.76, p < 0.0001, and middle frontal gyri, F(1,82) = 21.69, p < 0.0001.

Figure 2

Time course (in ms after stimulus onset at 0) of estimated neurophysiological current associated with word reading and recognition in four ROIs [in the left (LH) and right hemisphere (RH)]. Values on the vertical axis represent the difference between the peristimulus current in each ROI and estimated current in each subsequent time bin. Time windows of significant group differences are marked by boxes. BA22, superior temporal gyrus; BA21, middle temporal gyrus; BA46, rostral middle frontal cortex; BA37, fusiform gyrus.

Given previous reports of hemisphere-dependent reading-group effects in these areas, marginally significant time by hemisphere by reading group interactions were explored further. One such interaction was found in FUS, F(13,1066) = 2.01, p < 0.015. Inspection of the source current density waveforms in Figure 2 indicates that increased degree of activity for the RD group persisted longer in the RH as compared to the LH. Indeed the simple main effect of reading group was significant in both hemispheres between 150 and 300 ms (p < 0.0001), but only in the RH between 300 and 450 ms [F(1,82) = 13.11, p < 0.001; p > 0.05 in the LH]. The waveforms for each reading group were essentially identical between hemispheres in MTG, STG, and middle frontal gyrus (interactions involving hemisphere and reading group were all associated with p values > 0.1).

Greater degree of left than right hemisphere activity across reading groups was found in the LOC cortex [as indicated by a hemisphere main effect F(1,82) = 7.41, p < 0.008]. Latency-dependent hemispheric asymmetries for both groups, as indicated by significant time by hemisphere effects, were found in MTG [200–300 ms, F(1,82) = 9.01, p < 0.003], STG [400–650 ms, F(1,82) = 10.55, p < 0.001], SMG [300–700 ms, F(1,82) = 12.77, p < 0.0001], fusiform gyrus [200–350 ms, F(1,82) = 8.32, p < 0.004], and middle frontal gyrus [350–550 ms, F(1,82) = 18.54, p < 0.0001]. There were no significant effects involving reading group or hemisphere on peak latency.

Correlations between degree of activity and achievement

Although the two groups were comparable in age, moderate yet significant negative correlations were found between age and achievement scores (WJ-III word attack, letter–word identification, and spelling), as well as between age and degree of activity in the majority of ROIs (both within and across reading groups). Accordingly, the relation between degree of activity and achievement scores was explored by controlling for the effects of age on both sets of variables. Initially, and in order to preserve the temporal information that is inherent in the MNE time series data, partial correlation coefficients were computed between estimates of neurophysiological activity (adjusted for peristimulus estimated current) for each time bin and ROIs and each of the achievement measures. These analyses were performed separately for NI and RD children and evaluated at p < 0.005 in order to control for family-wise Type I error.

For NI readers significant positive correlations were obtained between WJ-III Spelling scores and degree of activity in the following left hemisphere regions: MTG (between 500 and 550 ms; r = 0.47), STG (between 400 and 500 ms: r = 0.48), and fusiform gyrus (between 200 and 250 ms: r = 0.49). Significant positive correlations were also noted between two peak latency measures and spelling (left MTL: r = 0.42 and left fusiform gyrus: r = 0.41). Correlations with letter–word identification scores were generally smaller (ranging between r = 0.18 and 0.30) and failed to reach significance. Correlations with WJ-III word attack scores were negligible (r < 0.13). For the RD group, only negative correlations were found between spelling scores and two right hemisphere ROIs, MTL (between 250 and 450 ms: r = −0.52) and fusiform gyrus (between 350 and 600 ms: r = −0.48). There were no significant correlates of achievement measures among peak latency measures in this group.

Each of the MNE estimates of activity listed in the previous paragraph were entered as independent variables in a series of hierarchical multiple regression analyses. The dependent variable in each of these analyses was WJ-III spelling or letter–word identification scores. Age was always entered in the first step of the analysis, degree of activity measures in the second step, and peak latency measures in the third step of each model. For the NI group degree of activity in left posterior temporal cortices (MTG, STG, fusiform gyrus) accounted for 30% of the variance in spelling scores [Adj R2 = 0.30, F(4,36) = 4.30, p < 0.006]. Adding peak latency measures significantly improved the model [Adj R2 = 0.41, F(6,34) = 5.30, p < 0.001]. A similar pattern emerged for letter–word identification scores with degree of activity measures as predictors [Adj R2 = 0.16, F(4,36) = 2.30, p < 0.04], and after addition of peak latency measures into the model [Adj R2 = 0.39, F(6,34) = 4.86, p < 0.001]. Figure 3 displays the degree of data fit (concordance between measured spelling performance and standardized predicted spelling scores based on the linear combination of the independent variables in the regression equation). Entering VIQ in the first step along with age did not affect the overall regression results or the individual beta weights of the remaining independent variables. The prediction models for Word Attack scores were not significant. For the RD group prediction models for each of the three achievement measures did not reach significance (p > 0.2 in all cases).

Figure 3

Regression plot of measured Woodcock–Johnson III (WJ-III) standard spelling scores (vertical axis) over predicted standard spelling scores (horizontal axis) based on a linear combination of the following MEG variables: adjusted degree of activity in left posterior ROIs (MTG [500–550 ms], STG [400–450 ms], and VOC [fusiform gyrus; 250–300 ms]) and peak latency in the left MTL and VOC among NI students (Adj. . The regression model is described by the following equation: Spelling = 3.065 × DMTG + 6.195 × DSTG + 5.164 × DVOC + 0.030 × LVOC + 0.040 × LMTL + 88.436, where D and L denote degree of activity and peak latency, respectively.

Discussion

In the present study, distributed source modeling of reading-related MEG data derived from a large, representative sample of typical and struggling readers partially replicated and extended previous reports (Simos et al., 2000a; Shaywitz et al., 2002; Booth et al., 2007; Pugh et al., 2008) regarding the nature of the functional disruption of the brain mechanism for word reading in RD. Sample selection ensured that the distribution of key demographic variables (age, gender, ethnic background), in addition to handedness and general cognitive ability (measured by PIQ), did not differ significantly between groups. Spatiotemporal brain activation profiles were obtained in the context of a continuous recognition task, the level of difficulty of which was titrated to reduce performance differences between groups compared to previous studies (between 70 and 80% correct identification of printed words from the training set). This was achieved by using shorter, higher frequency words and a smaller target set compared to our earlier study (Simos et al., 2000a). In this manner demands for both decoding and word recognition as well as memory encoding and retrieval were kept relatively low for both groups. The task was, nevertheless, successful in activating all key reading-related brain areas with clear left hemisphere dominance in their degree of activation as has been our experience with previous language studies (Simos et al., 2001; Papanicolaou et al., 2004).

The main findings of the present study can be summarized as follows. Firstly, children with RD showed reduced degree of neurophysiological activity in posterior temporal regions (MTG and STG) in both hemispheres. This activity was associated with relatively late stages of stimulus processing (200–650 ms after stimulus onset). Secondly, children with RD showed significant overactivation of the fusiform gyrus and middle frontal areas (between 150 and 450 ms), bilaterally, independent of age. Thirdly, individual differences in the degree and latency of activity in posterior temporal, ventral occipitotemporal (VOC), and mesial temporal regions mapped onto print-related achievement measures.

Although reading group differences in MTG and STG were found in both hemispheres, two additional lines of evidence suggest that disrupted functional brain organization for reading was restricted to the left hemisphere in RD children. Thus, while the degree of activity in left posterior temporal regions was a significant predictor of reading and spelling skill among NI readers, such an association was not noted among RD children. In the latter group, enhanced activity in right occipitotemporal and mesial temporal regions was associated with lower achievement scores. Importantly, reading-skill related differences in the degree of activity in these regions were evident across the age span covered by the present study (7–15 years). Moreover, for NI children, a clear temporal separation was noted between peaks of neurophysiological activity in left occipitotemporal (lateral and ventral visual association cortices), posterior temporal (middle and superior temporal gyri), and inferior/middle frontal regions. In contrast, the regular temporal progression of activity peaks between left hemisphere occipitotemporal, posterior temporal, and inferior/middle frontal regions was disrupted in the RD group.

Temporoparietal cortex and reading ability

There is ample evidence that neurophysiological processes that take place in the posterior portion of MTG (BA 21) reflect lexical/semantic processing (Damasio and Damasio, 1983; Pugh et al., 1996; Gaillard et al., 2001; Booth et al., 2002; Fiebach et al., 2002; Halgren et al., 2002; Simos et al., 2002b; Turkeltaub et al., 2002; McCandliss et al., 2003; Tyler et al., 2005; Wehner et al., 2007a). These operations may include storage or gaining access to lexical/semantic information for printed stimuli (Pugh et al., 1996; Fiebach et al., 2002). Damage to this region has been linked to deficits in a variety of tasks tapping into lexical/semantic processing, such as semantic priming (Tyler et al., 2005). Moreover, the onset latency of magnetic activity in this region significantly correlated with pronunciation latency of orthographically exceptional words while this was clearly not the case for regular words and pseudowords (Simos et al., 2002b).

With respect to cortex lying on the lateral and superior aspects of the STG (BA 22), extending posteriorly into the supramarginal gyrus, the bulk of empirical evidence implicates neurophysiological processes taking place in this region in phonological processing of spoken and written language (Beauvois and Derouesne, 1979; Caplan et al., 1995; Papanicolaou et al., 2003; Specht et al., 2003; Majerus et al., 2005). These processes include, but are not limited to, storage of sublexical phonological representations, automatic retrieval and short-term maintenance of these representations in consciousness (Hickok and Poeppel, 2000; Wise et al., 2001; Simos et al., 2002b; Mustovic et al., 2003). It should be noted, however, that both imaging and lesion data do not clearly indicate exclusive involvement of the left STG in phonological processing. There is also evidence that this region plays a crucial role for lexical/semantic processing of spoken language. Patients exhibiting the clinical profile of fluent aphasia, which typically results from damage to this region (Benson, 1985; Vignolo, 1988), almost invariably display notable deficits in a variety of experimental tasks tapping into lexical/semantic processing (Hagoort, 1993; Rodd et al., 2005; Janse, 2006). There are also reports that at least in adult experienced readers, the posterior portion of BA 22 is involved in lexical-semantic processing during silent reading tasks (Fiez et al., 1999; Cannestra et al., 2000; Haist et al., 2001; Fiebach et al., 2002; Okada and Hickok, 2006; Simos et al., 2009).

Evidence of increased activation in RD

In contrast to the bilateral reduction in MTG and STG activity, a significant increase in the strength of activation in middle frontal regions was noted for the RD group concurring with similar findings from PET and fMRI (Shaywitz et al., 1998; Brunswick et al., 1999; Grunling et al., 2004), as well as our more recent MEG investigations (Simos et al., 2007). It has been suggested that elevated levels of neurophysiological activity in prefrontal regions indicates engagement of cognitive operations supporting alternative strategies employed by individuals with RD to compensate for a deficiency in mapping the sound structure of printed words through covert articulation (Shaywitz et al., 1998; Shaywitz et al., 2004). However, the putative compensatory increase in prefrontal activity in RD is not a universal finding (e.g., Rumsey et al., 1997; Eden et al., 2004). In a recent meta-analysis of nine fMRI and PET studies Maisog et al. (2008) found no evidence to suggest a universal pattern of bilateral prefrontal hyperactivation in RD, but this meta-analysis involved only adults. In a meta-analysis with no age restrictions, Richlan et al. (2009) reported underactivation of the IFG gyrus and hyperactivity of the primary motor cortex and anterior insula. Several lines of evidence indicate that prefrontal hyperactivation may not reflect a key feature of the brain circuit for reading in RD, but is related to task and age-related processes (Hoeft et al., 2007).

Prefrontal hyperactivation in children with RD in the present study may also reflect additional demands for both encoding and retrieval operations posed by our continuous word recognition task. A comparable finding of increased memory-related prefrontal hemodynamic response among RD children has recently been reported (Vasic et al., 2008). There is ample evidence from lesion and functional imaging studies regarding the critical role played by prefrontal (Knowlton and Squire, 1995; Skinner and Fernandes, 2007; MacPherson et al., 2008) and mesial temporal cortices (Brown and Aggleton, 2001; Turriziani et al., 2004; Wixted and Squire, 2004) for recognition of previously encountered stimuli. The notion that task demands were related to increased prefrontal activity in RD is consistent with our finding that this effect was significantly larger among younger students, given that for this group memory/attention demands are expected to be higher than among older participants, rendering our continuous recognition task even more taxing among young RD students.

Increased strength of activation for children with RD was also found in VOC cortices (anatomically corresponding to the fusiform gyrus). This region is suitable for hosting neurophysiological processes responsible for graphemic processing and also possibly for the integration of orthographic, phonological and morphological information derived from print (McCandliss et al., 2003; Flowers et al., 2004; Pammer et al., 2004; Vigneau et al., 2005). Accumulating evidence suggests that VOC cortices are involved in storing or gaining access to familiar orthographic representations, which may take the shape of “word forms” (Kronbichler et al., 2004; Simos et al., 2009), and play an important role in both word recognition and spelling (Philipose et al., 2007; Mani et al., 2008). In agreement with previous reports, group differences in activity in this region were found during early processing stages at the peaks of the current density waveforms (between approximately 150 and 350 ms; e.g., Tarkiainen et al., 1999; Usui et al., 2009). One possibility is that fusiform activation among struggling readers may reflect increased effort to engage a visual/orthographic strategy for encoding and recognition of the printed word stimuli. This possibility is consistent with behavioral evidence in children and adults with reading disability (Lennox and Siegel, 1996; Milne et al., 2003; McNeil and Johnston, 2004, 2008; Best and Howard, 2005; Zoccolotti et al., 2005; Miller and Kupfermann, 2009). An alternative account of increased VOC activation relies on the putative role of this area in graphemic processing using linguistic constraints (Vigneau et al., 2005; Devlin et al., 2006) and in language functions in general (Schäffler et al., 1994).

Timing of regional activation

Another key finding of the represent study was the lack of temporal differentiation between activity peaks in VOC and the left MTG and also between left temporoparietal and IFG ROIs in the RD group. Similar results have been reported previously in MEG studies with younger struggling readers (Simos et al., 2005, 2007, in press). In contrast, non-impaired readers showed a systematic progression of activity peaks between VOC, middle temporal, temporoparietal, and IFG regions. In the absence of empirical data establishing direct links between the temporal progression of activity between any two cortical regions and anatomical connections between these regions, interpretations of these findings can only be speculative at present. The timing data, however, present an intriguing parallel with findings of reduced white matter integrity in poor readers (Deutsch et al., 2005; Niogi and McCandliss, 2006) and indirect evidence from hemodynamic studies for reduced functional “connectivity” between VOC and temporoparietal regions (Pugh et al., 2000; Cao et al., 2008) and between temporoparietal and frontal regions (Hampson et al., 2006).

Comparison with previous MEG studies of RD

A notable difference between the current findings and those of our earlier MEG study (Simos et al., 2000a) concerns the hemispheric profile of reading-group activation differences during the visual-word recognition task. Significant right hemisphere predominance of temporoparietal activity had been observed earlier for RD children, whereas in the present study both groups of children demonstrated L > R activation of posterior temporal cortices. Moreover, whereas RD children had previously been distinguished by a noticeable reduction in activation of the left temporoparietal cortex and accompanying hyperactivation of the right homotopic region, bilateral group effects were found here, namely reduced activation of the left and right STG and MTG, relative to NI children. Although, bilateral hemodynamic hypoactivation in posterior temporal and inferior parietal areas has sometimes been reported (Rumsey et al., 1997; McCrory et al., 2000; Meyler et al., 2007), these discrepancies merit some consideration.

One possible explanation may concern differences in sample selection: The present investigation employed a more representative sample of children selected from schools as opposed to community recruitment, and a slightly more lenient reading achievement criterion to form the study groups. This account is not supported by additional exploratory analyses comparing two age- and PIQ-comparable groups of children derived from the current data set, one with severe RD (N = 20, with word attack, letter–word identification, and spelling scores below 75 points) and the second with high-average reading achievement scores (over 95 points on each of the three subtests; N = 22). The results were essentially identical to those from the entire sample without notable increases in the effect size of interactions involving hemisphere.

Another potential explanation may relate to differences in task demands. The task used in the present investigation was significantly easier to perform than the version used in Simos et al. (2000a) as evidenced by comparable performance between reading achievement groups. Evidence from neuroimaging studies in individuals with RD have shown the degree of brain of activity to be affected by reading task difficulty (Cao et al., 2006), including lateralization of STG and MTG activation (Pugh et al., 2008). Word reading tasks, especially if they do not pose specific demands for phonological processing are likely to be more susceptible to the effects of overall task difficulty, and result in reduced hemispheric asymmetries in perisylvian activation (e.g., McCrory et al., 2000). A similar effect appears to hold for auditory processing as well (e.g., Breier et al., 2003; Wehner et al., 2007b).

Furthermore, the use of a higher-density gradiometer array in conjunction with distributed source modeling may account for the additional group-related findings in prefrontal and mesial temporal cortices. By affording increased sensitivity to smaller sources (mainly in terms of the spatial extent of their respective magnetic surface maps), the current MEG system was more likely to reveal multiple sources in each hemisphere which were invisible or difficult to model using the single-dipole method. Similarly, the potential for detecting smaller, more extended sources of activation is likely to have been facilitated by the use of distributed source modeling.

Conclusion and Future Directions

Although task demands were intentionally reduced to ensure high levels of performance among students with RD, clear indications of atypical functional organization of the brain mechanism for word recognition were detected. In addition to assessing group differences in average degree of regional neurophysiological activity, our results explored the relation between individual differences in achievement scores and degree of regional activity. Significant positive relations were restricted to the NI group in agreement with recent functional (Hoeft et al. 2007) and anatomical studies (Pernet et al., 2009). These findings further corroborate results from this and previous studies suggesting that atypical functional organization of cortical networks for reading form the core pathophysiology of developmental RD, whereas reliance on operations that take place in prefrontal and right VOC cortices (among other areas) reflects a relatively inefficient mechanism for word recognition. An important caveat of the source mapping approach used in the present study, however, is that ROI selection was restricted by the cortical segmentation procedure available in FreeSurfer, resulting in relatively large ROIs (e.g., the entire extent of MTG). Future studies adopting a voxel-based approach may provide more detailed spatiotemporal data associated with engagement of the brain circuit for reading.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.