Resting-state functional connectivity (RSFC) approaches offer a novel tool to delineate distinct functional networks in the brain. In the present functional magnetic resonance imaging (fMRI) study, we elucidated patterns of RSFC associated with 6 regions of interest selected primarily from a meta-analysis on word reading (Bolger DJ, Perfetti CA, Schneider W. 2005. Cross-cultural effect on the brain revisited: universal structures plus writing system variation. Hum Brain Mapp. 25: 92–104). In 25 native adult readers of English, patterns of positive RSFC were consistent with patterns of task-based activity and functional connectivity associated with word reading. Moreover, conjunction analyses highlighted the posterior left inferior frontal gyrus and the posterior left middle temporal gyrus (post-LMTG) as potentially important loci of functional interaction among 5 of the 6 reading networks. The significance of the post-LMTG has typically been unappreciated in task-based studies on unimpaired readers but is frequently reported to be a locus of hypoactivity in dyslexic readers and exhibits intervention-induced changes of activity in dyslexic children. Finally, patterns of negative RSFC included not only regions of the so-called default mode network but also regions involved in effortful controlled processes, which may not be required once reading becomes automatized. In conclusion, the current study supports the utility of resting-state fMRI for investigating reading networks and has direct relevance for the understanding of reading disorders such as dyslexia.
Reading requires dynamic interactions among multiple brain systems. As early as the mid-20th century, Geschwind (1965) claimed that language and reading disorders arise from diminished connectivity between distinct functional systems in the brain. The recent advent of neuroimaging techniques such as positron emission tomography (PET) and functional magnetic resonance imaging (fMRI) have provided researchers with a means of examining connectionist frameworks for reading disorders in the living human brain. For example, using PET, Horwitz et al. (1998) demonstrated that during single-word reading, dyslexic readers, when compared with unimpaired readers, exhibited decreased functional connectivity of the left angular gyrus with the left frontal, temporal, and occipito–temporal areas. This finding was replicated in part by a later fMRI study, in which dyslexic readers showed diminished functional connectivity between the left angular gyrus and the left occipito–temporal areas during tasks that required intensive phonological processing (Pugh et al. 2000).
Despite these initial successes, researchers continue to struggle with the challenge of developing effective probes for examining reading networks, in both nonclinical and clinical populations. In other words, the field lacks consensus regarding the “ideal task” or set of tasks for probing reading networks, yielding inconsistent findings across studies of reading. For example, the specific task (e.g., naming words and semantic categorization) employed in fMRI studies of word reading can produce marked variation in the locations of peak activations (Booth et al. 2006; Nakamura et al. 2007) and in patterns of functional connectivity (Mechelli et al. 2005; Nakamura et al. 2007). Additionally, given that reading abilities are dependent on both age and educational level, researchers face the challenge of devising probe tasks appropriate for a broad range of ages and abilities. Such difficulties are reflected in inconsistent findings from studies of unimpaired (Turkeltaub et al. 2002; Bolger et al. 2005) and dyslexic (Richlan et al. 2009) readers. For example, the importance of the left angular gyrus in reading is suggested by reduced task-related activation in these areas in dyslexic readers (Shaywitz et al. 1998; Richlan et al. 2009), but this region is not reported to be strongly activated during word-reading tasks in unimpaired readers (Turkeltaub et al. 2002; Bolger et al. 2005). As such, task selection can confound the examination of neural mechanisms of reading and reading disorders.
A potential solution to this problem is offered by resting-state functional connectivity (RSFC) approaches. In recent years, the application of these approaches to fMRI data has emerged as a powerful and highly efficient method for mapping and assessing the functional architecture of the brain, without task constraints (Biswal et al. 1995). When individuals are scanned at rest, low-frequency (<0.1 Hz) spontaneous fluctuations in the blood oxygen level–dependent (BOLD) signal are temporally synchronized between functionally related brain areas (Biswal et al. 1997; Lowe et al. 1998; Greicius et al. 2003; Damoiseaux et al. 2006; Margulies et al. 2007). Such RSFC patterns are considered to represent intrinsically organized functional networks in the brain (see review by Fox and Raichle 2007). Among identified intrinsic connectivity networks, the so-called default mode network has been most intensively investigated (Raichle et al. 2001; Raichle and Snyder 2007). This network includes the medial prefrontal cortex, posterior cingulate cortex, and precuneus and is typically deactivated during goal-oriented tasks. Spontaneous BOLD fluctuations within the default mode network exhibit negative correlations with fluctuations in brain regions that are typically activated during the performance of attention-demanding tasks (so-called task-positive regions, Fox et al. 2005). Interestingly, alterations in the RSFC of the default mode network have been observed in several clinical populations, such as individuals with Alzheimer's disease (Li et al. 2002; Greicius et al. 2004), attention-deficit hyperactivity disorder (Castellanos et al. 2008; Uddin et al. 2008), and schizophrenia (Liang et al. 2006; Zhou et al. 2007).
Previous studies have successfully applied RSFC approaches to the study of language and reading networks. Hampson et al. (2006) showed that patterns of functional connectivity among reading-related areas are strikingly consistent across reading and rest conditions. In particular, Broca's area (i.e., the left inferior frontal gyrus [IFG]) was significantly correlated with the left angular gyrus (BA39) as well as with the left occipitotemporal cortex (i.e., fusiform gyrus [FFG]). Intriguingly, the strength of functional connectivity between Broca's area and the left angular gyrus during both reading and rest differentiated better readers from worse readers (although none were impaired readers). Most recently, Lohmann et al. (2009) demonstrated that functional connectivity between language areas observed during 4 different language tasks remained even after regressing out specifics of the task paradigms and low-pass filtering (<0.01 Hz). This finding supports the suggestion that patterns of RSFC reflect an intrinsic functional organization underlying cognitive processes. Furthermore, these initial results support the utility of RSFC methods for further characterization of the neural networks underlying reading.
In order to build on the work of Hampson et al. (2006), we examined patterns of RSFC associated with 6 regions of interest (ROIs) that were found to be consistently activated by word reading tasks in 2 recently published meta-analyses (Turkeltaub et al. 2002; Bolger et al. 2005). Given that even subtle differences in the specific paradigms and cognitive constructs examined can produce marked variation in findings across task-based studies, we selected our ROIs on the basis of a meta-analysis rather than on any single study so as to best sample the regions activated by word-reading tasks. We hypothesized that patterns of positive RSFC (positive correlations) associated with the 6 reading-related ROIs would be consistent with previously shown patterns of task-related activation and functional connectivity during word reading. Second, we tested whether patterns of negative RSFC (negative correlations) associated with the 6 reading-related ROIs would coincide with the default mode network. Finally, we aimed to identify potentially important “loci of functional interaction,” defined as regions of overlap among the 6 reading networks, by performing a conjunction analysis. The networks and loci of interaction identified using this approach have clinical relevance in that they provide targets for investigation in reading disorders such as dyslexia.
Materials and Methods
Twenty-five right-handed healthy adults (mean age: 32.0 ± 8.4 years, 19 males) participated in the current fMRI study. Data from these participants have been included in several previously published studies from our laboratory (e.g., Margulies et al. 2007; Di Martino et al. 2008; Shehzad et al. 2009). All participants were college educated and native speakers (and readers) of English, with no history of neurological, psychiatric, or learning disorders, including dyslexia, as confirmed by both clinical interview and questionnaires. The study was approved by the institutional review boards of the New York University School of Medicine and New York University. Prior written informed consent was obtained from all participants.
fMRI Data Acquisition
MRI data were collected on a Siemens Allegra 3.0-T scanner at the NYU Center for Brain Imaging. Scans were brief (6 min 38 s) and comprised of 197 continuous echo planar imaging functional volumes (time repetition [TR] = 2000 ms; time echo [TE] = 25 ms; flip angle = 90; 39 slices, matrix = 64 × 64; field of view [FOV] = 192 mm; acquisition voxel size = 3 × 3 × 3 mm). During the scan, participants were instructed to relax with their eyes open while the word “Relax” was projected at the center of the display. A high-resolution T1-weighted anatomical image was also acquired using a magnetization prepared gradient echo sequence (TR = 2500 ms; TE = 3.93 ms; inversion time = 900 ms; flip angle = 8; 176 slices; and FOV = 256 mm).
Preprocessing steps for slice timing correction, motion correction, and despiking (detection and reduction of extreme time series outliers) were performed using analysis of functional neuroimages (AFNI; Cox 1996). Further preprocessing steps were performed using FMRIB Software Library (FSL) (Smith et al. 2004) and comprised spatial smoothing (using a Gaussian kernel of full width half maximum 6 mm), mean-based intensity normalization of all volumes by the same factor, temporal bandpass filtering (0.01–0.1 Hz), and prewhitening (correction for time series autocorrelation). Finally, each individual's time series was spatially normalized by a 12-degrees-of-freedom affine registration to the MNI152 standard brain template (Montreal Neurological Institute), with 2-mm3 resolution.
Nuisance Signal Regression
To control for the effects of physiological processes (such as fluctuations related to cardiac and respiratory cycles) and motion, we removed signal associated with several nuisance covariates. Specifically, we regressed each subject's preprocessed 4-D volume on 9 predictors that modeled nuisance signals from white matter, cerebrospinal fluid, the global signal, and 6 motion parameters, as detailed elsewhere (Kelly, de Zubicaray, et al. 2009; Kelly, Di Martino, et al. 2009).
ROI Selection and Seed Generation
We selected a total of 6 ROIs on the basis of 2 meta-analyses of single word reading in alphabetic languages (Turkeltaub et al. 2002; Bolger et al. 2005). These meta-analyses used an activation likelihood estimate approach to identify several foci reliably activated across studies of single word reading. First, we identified the 5 brain areas that were identified as being colocated in regions defined by the Harvard–Oxford Cortical Structural Probabilistic Atlas (Kennedy et al. 1998; Makris et al. 1999) in both meta-analyses (the left FFG, the left superior temporal gyrus (STG), the left temporoparietal junction (TPJ), the left precentral gyrus (PCG), and the left IFG). Because the peak coordinates between the 2 meta-analyses differed by 17 mm for the left FFG and 14 mm for the left IFG, we used the coordinates reported by the more recent Bolger et al. (2005) meta-analysis, which included more contrasts (i.e., 38 contrasts from 25 studies) than Turkeltaub et al. (2002), which contained 16 contrasts from 11 studies. Limiting coordinate selection to 1 meta-analysis also avoided potentially over representing regions. Of note, the 5 regions selected have also been reported to be hypoactivated in dyslexic readers compared with unimpaired readers (Richlan et al. 2009).
Second, we added a sixth ROI located in the posterior part of the left inferior occipital gyrus. This cortical area was reported only in Bolger et al.'s (2005) meta-analysis but was included in light of the clinical observation that damage to the posterior occipital area, particularly in the left hemisphere, can result in pure alexia characterized by letter-by-letter reading (Behrmann et al. 1998; Sakurai et al. 2001). We verified that all ROIs were at least 12 mm apart in Euclidean distance to ensure that the regions were spatially distinct. For each ROI, we converted the reported Talairach coordinates into MNI coordinates (Brett et al. 2001) and created a spherical seed centered on these coordinates in 2 × 2 × 2-mm space, with a radius of 6 mm.
The 6 seed ROIs were spatially distributed in the left hemisphere (Fig. 1) and labeled according to the Harvard–Oxford Cortical Structural Atlas (Kennedy et al. 1998; Makris et al. 1999); 1) inferior occipital gyrus (IOG) (x = −25, y = −87, z = −10), in the medial posterior part of the left inferior occipital gyrus, part of the visual ventral pathway; 2) FFG (x = −48, y = −57, z = −20), in the posterior part of the left FFG, at a location coinciding with the “Visual Word Form Area” that is responsive to word and word-like orthography (Cohen et al. 2000; Dehaene et al. 2002); 3) STG (x = −53, y = −31, z = 9), in the posterior part of the left STG, which is included in Wernicke's area and involved in speech perception (Scott and Wise 2004); 4) TPJ (x = −59, y = −45, z = 15), in the left TPJ, where the ventral posterior part of the left supramarginal gyrus meets the posterior part of STG. This region's activation is often observed in association with grapheme-phoneme conversion (phonological decoding, Fiez and Petersen 1998; Turkeltaub et al. 2003); 5) PCG (x = −48, y = −12, z = 45), in the dorsal part of the left PCG, part of the primary motor cortex; and 6) IFG (x = −51, y = 10, z = 10), in the pars opercularis of the left IFG, which is included in Broca's area and involved in speech articulation (Wise et al. 1999; Nixon et al. 2004; Owen et al. 2004). The location of each seed is shown in MNI space in Figure 1. For each participant, we calculated the mean time series of each seed by averaging across all voxels within the seed.
We note that our set of ROIs did not include BA39, the activation level and connectivity strength of which are diminished in dyslexic readers (Horwitz et al. 1998; Shaywitz et al. 1998; Pugh et al. 2000). This was because neither of the recent meta-analyses of unimpaired reading (Turkeltaub et al. 2002; Bolger et al. 2005) highlighted the left angular gyrus as a locus of consistent activation for word reading. The exclusion of the left angular gyrus can be further rationalized given that this region's functional significance in reading appears to be specific to less experienced readers (children; Church et al. 2008) and individuals with reading disorders (Dejerine 1891).
Subject- and Group-Level RSFC Maps
For each participant, we performed a multiple regression analysis (using the general linear model implemented in the FSL program FEAT) for each seed. The mean time series of each seed was calculated by averaging across all voxels within the seed for each subject, using the AFNI program 3dmaskave. This analysis produced individual subject-level maps of all voxels that were positively and negatively correlated with the seed-time series. Group-level analyses were carried out using a mixed-effects model as implemented in the FSL program FLAME. Corrections for multiple comparisons were conducted at the cluster level for each RSFC map (Z > 2.3; P < 0.05, corrected). The group analysis created 6 positive and 6 negative thresholded Z-statistic maps (RSFC maps).
To identify loci of functional interactions among intrinsic reading networks, the 6 group-level thresholded maps of positive RSFC were each binarized and then summed to create a conjunction map. The resultant map was then thresholded to identify areas that are common to all 6 RSFC maps and to 5, 4, 3, and 2 of the 6 RSFC maps. Similarly, we conducted a conjunction analyses for the 6 maps of negative RSFC. To reduce the rate of false positives, areas of overlap were required to comprise >10 voxels.
For each of 6 brain regions previously implicated in reading (Bolger et al. 2005), seed-based correlation analyses were employed to characterize their associated functional systems during rest. Figure 2 depicts positive RSFC for each of the 6 reading-related seeds (second to seventh columns). Detailed descriptions of each RSFC pattern and peak coordinates of significant positive correlations are found in Supplementary Document 1 and Supplementary Table 1, respectively.
Overall, the patterns of RSFC associated with each of our 6 seed ROIs were consistent with reading networks previously identified using task-based approaches. Most notably, RSFC between the TPJ seed and left frontal and temporal areas showed remarkable similarities to previously shown functional connectivity patterns of the left angular gyrus (in the vicinity of the TPJ seed used in the present study) during reading tasks (Horwitz et al. 1998; Pugh et al. 2000) as well as during rest (Hampson et al. 2006). In addition, the FFG seed exhibited long-range RSFC with the left IFG extending into the PCG. This was exactly mirrored by long-range RSFC of the IFG seed with the left FFG. Such long-range connectivity patterns between occipito-temporal and frontal areas have been robustly demonstrated by task-based fMRI studies on reading (Mechelli et al. 2005; Bitan et al. 2006). Of note, RSFC patterns for the IFG seed also were comparable with those demonstrated by Hampson et al. (2006). In sum, even though no reading tasks or word stimuli were included, our findings were highly consistent with those previously reported in the task-based reading literature.
As shown in the red-framed column (“overlap”) of Figure 2, a broad array of brain areas covering much of the cerebral cortex was found to be common to one or more of the reading-related seeds. In order to identify “loci of functional interaction” among the 6 RSFC networks, we carried out a series of conjunction analyses. Although no voxels exhibited significant positive correlations with all 6 seeds, 2 loci in the left hemisphere were common to 5 of the 6 RSFC maps (Fig. 3, shown in red). The first was located within the posterior part of the left middle temporal gyrus (post-LMTG; center of gravity in MNI space: x = −48, y = −62, z = 4; cluster extent = 82 voxels). The second region of overlap was located in the posterior boundary of the left IFG (post-LIFG; center of gravity in MNI space: x = −54, y = 6, z = 18; cluster extent = 47 voxels). Of note, inspection of the data revealed that the IOG seed was the only seed that did not exhibit connectivity with these 2 loci. Finally, Figure 3 shows that overlap between 3 (in green) and 4 (in yellow) of the 6 correlation maps was evident in multiple brain areas, distributed across frontal, motor, insula, temporal, and parietal regions, bilaterally.
Brain areas that exhibited significant negative connectivity with one or more of the 6 seeds included the medial prefrontal cortex, dorsolateral prefrontal cortex, anterior cingulate cortex, superior lateral parietal cortex, lateral temporal cortex, posterior cingulate cortex, and adjacent precuneus (Fig. 3). Many, but not all, of these negatively connected brain areas are topographically consistent with the default mode network, which is typically deactivated during tasks and negatively correlated with task-active networks (Raichle et al. 2001; Fox et al. 2005; Toro et al. 2008). Similar to the positive connectivity maps, although no voxels exhibited significant negative correlations with all 6 seeds, there was overlap among 5 of the 6 negative RSFC maps (again, excluding the map associated with the IOG seed) in several cortical areas. These included the posterior cingulate cortex extending into the precuneus, the dorsolateral prefrontal cortex (i.e., the right middle frontal gyrus), and the superior lateral parietal cortex (i.e., the right angular gyrus; Fig. 3 in dark blue). In addition, Figure 3 illustrates loci of overlap among 3 (in light blue) and 4 (in purple) of the 6 negative RSFC maps, which extended outward from the aforementioned loci of overlap to comprise much of the classically defined default mode network.
Secondary Analysis: Post-LMTG
In order to confirm the importance of the post-LMTG as a locus of putative functional interaction for reading, we used the post-LMTG region that exhibited overlap among 5 of the 6 maps of positive RSFC as a seed in a secondary RSFC analysis (Fig. 4). This analysis allowed us to determine the extent to which the pattern of RSFC associated with this region overlapped with those associated with the 6 seed regions. The post-LMTG seed was positively connected with a range of cortical areas in both hemispheres, including the lateral occipital cortex, the lingual gyrus, the temporal cortices, the parietal cortices, the insula, and the precentral/postcentral gyri extending into the posterior boundary of the IFG. This map of positive RSFC comprised the locations of 5 of the 6 reading seeds, with the one exception being the IOG (Fig. 4). This confirms the observed result of the conjunction analysis, which suggested that the post-LMTG region is a locus of functional interaction among the 5 intrinsic reading networks. Notably, the post-LMTG seed was not positively connected with areas of the prefrontal cortex nor the posterior–medial part of the occipital gyrus, where the IOG seed was located.
In short, our resting-state fMRI results successfully replicated previous patterns of task-based activity and functional connectivity associated with word reading, even in the absence of reading tasks or word stimuli. Moreover, the results highlighted the post-LIFG and post-LMTG as loci of functional interaction among 5 of the 6 reading networks we examined.
We examined patterns of RSFC associated with each of 6 reading-related ROIs. We selected these ROIs on the basis of recent meta-analyses of single word reading (Turkeltaub et al. 2002; Bolger et al. 2005), rather than deriving them from a single study or on the basis of prior theoretical knowledge about language networks. There were 3 main findings. First, the observed patterns of positive RSFC largely replicated patterns of task-based activity and functional connectivity previously demonstrated during word reading, despite the absence of reading tasks or word stimuli. Second, patterns of negative RSFC included both the default mode network as well as several executive regions involved in effortful or controlled processing (e.g., dorsolateral prefrontal and superior parietal cortices; Cazalis et al. 2003; MacDonald 2008). Third, and most importantly, conjunction analyses highlighted the post-LIFG and the post-LMTG as potentially important loci of functional interaction among the reading networks. The post-LMTG was not highlighted in the recently published meta-analyses on word reading in unimpaired readers (Turkeltaub et al. 2002; Bolger et al. 2005) but has been reported as a region of dysfunction in dyslexic readers (Richlan et al. 2009). Our findings suggest the utility of RSFC approaches for studying the functional organization of the brain underlying reading and reading disorders.
Reading Networks at Rest
The first finding of the current study is that reading networks are intrinsically represented in spontaneous BOLD signal fluctuations, at least for experienced adult readers. For example, the observed RSFC between the left temporoparietal area (i.e., the TPJ seed) with temporal areas as well as with frontal areas in the left hemisphere has been robustly demonstrated during nonword reading in unimpaired readers (Horwitz et al. 1998; Pugh et al. 2000; Hampson et al. 2006). Functional connectivity in these circuits was diminished in individuals with dyslexia, suggesting the importance of its integrity to reading ability (Horwitz et al. 1998). Similarly, long-range patterns of left-hemisphere RSFC between frontal (i.e., the IFG seed) and posterior parietal areas as well as those between occipito–temporal (i.e., the FFG seed) and frontal areas are consistent with reading networks identified in task-based fMRI studies (Mechelli et al. 2005; Nakamura et al. 2007) and with reading networks identified in a previous resting-state fMRI study (Hampson et al. 2006). The observed long-range RSFC between the IFG and the FFG in the left hemisphere is also consistent with a recent magnetoencephalography study, in which these 2 regions showed temporal coactivation at the early stage (130–140 ms) of visual word recognition (Cornelissen et al. 2009).
The observed intrinsic representations for reading are likely to result from long-term experience of reading in our participants. RSFC patterns undergo developmental changes (Fair et al. 2007; 2008; Kelly, de Zubicaray, et al. 2009; Kelly, Di Martino, et al. 2009; Supekar et al. 2009), which may reflect experience-dependent (e.g., sensory and cognitive learning) changes in regional coactivation patterns (Dosenbach et al. 2007; Pinsk and Kastner 2007). In other words, experience likely modifies functional connectivity throughout development. Such experiential effects on brain networks are not limited to the highly plastic developing brain (Johnston 2009). Even in adults, learning experiences (e.g., motor learning) can alter RSFC patterns (Albert et al. 2009; Xiong et al. 2009). To test experience-dependent formation of RSFC associated with reading, future research should investigate developmental trajectories of intrinsic representations of reading in children and adolescents in relation to proficiency.
Bilaterality of RSFC
The observed RSFC patterns for reading were largely bilateral, rather than left lateralized, as expected for regions supporting language functions (including reading) in the human brain (Frost et al. 1999; Pujol et al. 1999; Springer et al. 1999; Powell et al. 2006). Findings of largely bilateral patterns of RSFC are not surprising, as the robust homotopic RSFC for the majority of brain regions has been documented (Salvador et al. 2005; Stark et al. 2008). With that said, Stark et al. (2008) noted that primary sensory/motor regions exhibit the highest degree of homotopic connectivity and that higher-order heteromodal association areas exhibit weaker homotopic connectivity. This observation of weaker homotopic connectivity among higher-order heteromodal association areas may reflect hemispheric specialization for higher cognitive functions, such as language. In this regard, we did observe some left-lateralized RSFC patterns, for example, between the left IFG and the left FFG (see Fig. 2).
Alternatively, task-free RSFC may not be sensitive enough to detect lateralization for reading, given that the degree of left lateralization is highly dependent on the nature of tasks (Spironelli and Angrilli 2006; Yang et al. 2009). Although the present study did not focus on the laterality of reading, investigating relationships between asymmetric RSFC patterns and reading performance remains of interest in future studies.
Interpretations of Negative Connectivity
We found evidence of negative connectivity between reading-related regions and 2 key sets of regions. First, we observed negative RSFC between reading-related regions and the posterior cingulate cortex extending into the precuneus, areas commonly included in the default mode network. This result is consistent with previous studies that note strong negative connectivity between task-positive networks and the default mode network (Raichle et al. 2001; Fox et al. 2005; Fransson 2005). The strength of such negative connectivity has been related to behavioral variability (e.g., reaction time in an attentional task; Kelly et al. 2008). Based on these findings, we propose that 1) negative connectivity may represent functional segregation or differentiation between distinct functional systems and 2) the strength of negative connectivity between such systems is a potential index of optimally balanced competition, which may be one of the keys to successful behavior. These theoretical accounts are consistent with the findings of Weissman et al. (2006), who demonstrated that less deactivation in the default mode network during an attention task was significantly related to inefficient behavior (i.e., longer response times). However, such a dichotomous view has been challenged by Sadaghiani et al. (2009), who demonstrated beneficial effects of increased activity in the default mode network on perceptual performance.
Second, we also found negative RSFC in the dorsolateral prefrontal and superior parietal cortices, which are strongly activated by effortful control tasks (Cazalis et al. 2003; MacDonald 2008) and tasks placing high demand on working memory (Marklund et al. 2007; Wendelken et al. 2008). This result is consistent with a recent report that negative connectivity between 2 task-positive networks (i.e., between the visual cortex and the visuospatial attention network) was induced by intensive visual perceptual training (Lewis et al. 2009). In most adults with extensive reading experience, word reading is a relatively automatized process and may not require (or may even require suppression of) high-level effortful control. This is best exemplified by the Stroop task, where word reading occurs obligatorily, despite controlled efforts to ignore the word (MacLeod 1991; Carter et al. 1995; Milham et al. 2001). Thus, the observed negative connectivity between the 2 task-positive regions and our reading seeds may reflect segregation of these functional organizations as a result of maturation of reading fluency/automatization. This assumption needs to be tested in future studies, particularly by comparing adults and children (i.e., less experienced readers).
Negative relationships, however, should be interpreted with some caution. Recent analyses have shown that negative RSFC is, at least in part, a consequence of preprocessing strategies such as global correction (Chang and Glover 2009; Fox et al. 2009; Murphy et al. 2009). On the other hand, intrinsic negative relationships are clearly present in direct neurophysiological recordings in animals (Walters et al. 2007; Hayden et al. 2009; Popa et al. 2009). Additionally, resting-state fMRI studies in humans have linked behavioral measures to the magnitude of intrinsic negative RSFC (Kelly et al. 2008; Di Martino et al. 2009), further supporting the potential functional significance of this measure.
Loci of Functional Interaction
The conjunction analysis demonstrated that RSFC for 5 of the 6 ROIs overlapped in the post-LMTG and the post-LIFG. These results suggest that these cortical areas represent important loci of functional interaction among intrinsic reading networks. Notably, the post-LMTG region was not included in the original set of reading seed regions nor was it reported in the aforementioned meta-analyses of normal word reading.
One may question why the meta-analyses failed to report the post-LMTG as a region reliably activated by word reading. This may be due to differences in the nature of the contrasts and tasks employed in the meta-analyses. One of the meta-analyses (Turkeltaub et al. 2002) strictly employed loose contrasts (e.g., reading aloud vs. rest or fixation), which tend to induce broad networks of reading-related activity that are often oversimplified in reporting (i.e., large clusters comprising several regions are represented by a single peak coordinate). Consequently, post-LMTG activations may be missed in meta-analyses performed on the basis of such peak coordinates. The other meta-analysis (Bolger et al. 2005) used both loose and tight contrasts (e.g., words vs. pseudowords). Both words and pseudowords involve visual–phonological recoding, subserved by the post-LMTG (Fiebach et al. 2002), which may consequently be activated by both reading conditions. Thus, tighter contrasts may underestimate this region's activation associated with word reading. However, when individuals with dyslexia have been studied, both abnormal (i.e., reduced) activations (Richlan et al. 2009) and normalized (i.e., increased) activations in response to intervention training (Simos, Fletcher, Sarkari, Billingsley, et al. 2007; Simos, Fletcher, Sarkari, Billingsley-Marshall, et al. 2007) have been detected in the post-LMTG. Recently, significant increases in gray-matter volumes were identified in several regions, including the post-LMTG, in a sample of late-literates who acquired reading skills in adulthood, relative to illiterate adults (Carreiras et al. 2009). These results further support the use of RSFC for investigating functional brain organization, particularly in clinical populations.
Considering that the post-LMTG is strongly activated both for reading real words (Fiebach et al. 2002) and for naming objects (Price et al. 2005), automatized phonological (and semantic) retrieval of visually presented items may be the cognitive process responsible for neural responses in the post-LMTG. Such automaticity for visually presented words is often impaired in dyslexic readers (Shaywitz and Shaywitz 2008). The post-LMTG region identified in our conjunction analysis is adjacent to the middle temporal area (MT), known to be the motion-sensitive cortex (Maunsell and Newsome 1987; Livingstone and Hubel 1988). Some, but not all, individuals with dyslexia are impaired in tasks that require sensitivity to visual motion (Stein and Walsh 1997; Talcott et al. 2000; Stein 2001) and show reduced activation in the MT during the perception of moving stimuli (Eden et al. 1996). Although the present study does not permit conclusive interpretation regarding the role of the post-LMTG in reading, it does suggest the potential merits of investigating RSFC patterns of the post-LMTG and their relationship with behavioral performance (e.g., reading and visual processing) in individuals with dyslexia.
Another locus of functional interaction, the post-IFG, is located immediately posterior to the pars opercularis of the IFG. This area has been a focus of literature investigating the functional organization of language (Bookheimer 2002; Davis et al. 2008; Keller et al. 2009) and reading (Paulesu et al. 1993; Cornelissen et al. 2009). Saur et al. (2008) combined diffusion tensor imaging with fMRI and reported that phonological processing (sound articulation) is subserved by structural connectivity between the posterior–frontal (including the post-LIFG) and superior–temporal regions. This is consistent with the classic Wernicke–Geschwind model of language (Geschwind 1965; Anderson et al. 1999) grounded in postmortem studies. From this finding, we can infer that the reported hyperactivation of the post-IFG in dyslexics (Shaywitz et al. 1998; Hoeft et al. 2007) reflects increased dependence on articulation mechanisms to compensate for impaired phonological processing. However, a recent subtraction meta-analysis of dyslexia arrived at the opposite result, that is, the post-IFG is hypoactivated in dyslexic readers (Richlan et al. 2009). Such an inconsistency may be ascribable to methodological differences such as the nature of tasks used or participant ages. Our results suggest the value of examining the functional connectivity of the IFG during rest, without task demands, in dyslexic participants.
Inferior Occipital RSFC
Notably, the IOG seed was not significantly correlated with either of the 2 identified loci of functional interaction. This suggests that the functional network associated with the IOG seed contributes to reading differently than the other 5 reading networks. Given that the lateral occipital cortex, including the IOG, is robustly activated during visual object recognition tasks (Grill-Spector et al. 2001; Kourtzi and Kanwisher 2000, 2001), the IOG network may be associated with visual perception of letters and words. This is consistent with patterns of deficits shown in alexic patients with damage to the IOG and adjacent posterior occipital areas; their deficits are not restricted to letter processing per se but are rather characterized by general impairments in visual perceptual processing including size identification (Behrmann et al. 1998; Sakurai et al. 2001). Moreover, the IOG activation pattern does not discriminate dyslexics from unimpaired readers, at least in alphabetic languages (Richlan et al. 2009). Although compelling, our examination of IOG connectivity only during rest leaves open the possibility that specific task demands (e.g., visual word recognition) could increase IOG coordination with reading-related areas. Effective connectivity approaches such as psychophysiological interaction (Friston et al. 1997, 2003; Hasson et al. 2009) are needed to determine whether task and visual demands impact IOG connectivity during reading.
Limitations of the Current Study
Our results suggest that reading-related networks are intrinsically represented in the human brain. Such representations likely reflect the result of learning and experience in adult readers. One way to test this hypothesis would be to map the patterns of RSFC for reading-related regions throughout development, as individuals systematically learn and gain reading experience (e.g., in children, adolescents, and adults). As reading is an acquired skill, RSFC patterns for reading in adults are expected to differ from those in children, who are less experienced with reading skills. Given robust developmental changes in RSFC patterns (Fair et al. 2007; Kelly, de Zubicaray, et al. 2009; Kelly, Di Martino, et al. 2009; Supekar et al. 2009), examining the developmental trajectory of RSFC associated with reading is a high priority in future studies.
Finally, another interesting avenue for research is the investigation of the extent to which RSFC strength predicts individual differences in reading-related performance (Hampson et al. 2006) . Unfortunately, formal measures of reading performance were not available in the present data set. Future work will address these questions.
In summary, our resting-state fMRI study demonstrated large-scale resting-state functional connectivity networks associated with word reading. The patterns of positive RSFC replicated patterns of task-based activity and functional connectivity, previously demonstrated during word reading, despite the absence of reading tasks or word stimuli. Our approach also revealed potentially important loci of functional interaction within reading networks. The significance of one of these loci, the post-LMTG, has been underestimated in task-based approaches aiming to localize neural correlates of reading. Importantly, this cortical area has been previously implicated in dyslexia (Richlan et al. 2009) and has been shown to exhibit intervention-induced changes in cortical activity in dyslexic children (Simos, Fletcher, Sarkari, Billingsley, et al. 2007; Simos, Fletcher, Sarkari, Billingsley-Marshall, et al. 2007). In conclusion, the current study suggests that resting-state fMRI can inform our understanding of the functional organization of the brain underlying reading and reading disorders.
NIMH grant (R01MH081218); Stavros Niarchos Foundation grant (to Dr Castellanos); Leon Levy Foundation grant (to Dr Milham); and gifts from Joseph Healy, Linda and Richard Schaps, and Jill and Bob Smith (Dr Castellanos).
The authors thank Dr Xinian Zuo for assistance with analyses and Mr Jonathan Adelstein for editorial assistance. Conflict of Interest: None declared.