Multimodal brain features at 3 years of age and their relationship with pre-reading measures 1 year later
Kathryn Y. Manning1,2,3*† 
Jess E. Reynolds1,2,3,4† Xiangyu Long1,2,3 Alberto Llera5,6 
Deborah Dewey2,3,7,8 
Catherine Lebel1,2,3- 1Department of Radiology, University of Calgary, Calgary, AB, Canada
- 2Alberta Children’s Hospital Research Institute, Calgary, AB, Canada
- 3Hotchkiss Brain Institute, University of Calgary, Calgary, AB, Canada
- 4Telethon Kids Institute, The University of Western Australia, Perth, WA, Australia
- 5Centre for Cognitive Neuroimaging, Donders Institute for Brain, Cognition and Behavior, Nijmegen, Netherlands
- 6Department of Cognitive Neuroscience, Radboud University Medical Centre, Nijmegen, Netherlands
- 7Department of Pediatrics, University of Calgary, Calgary, AB, Canada
- 8Community Health Sciences, University of Calgary, Calgary, AB, Canada
Pre-reading language skills develop rapidly in early childhood and are related to brain structure and functional architecture in young children prior to formal education. However, the early neurobiological development that supports these skills is not well understood. Here we acquired anatomical, diffusion tensor imaging (DTI) and resting state functional MRI (rs-fMRI) from 35 children at 3.5 years of age. Children were assessed for pre-reading abilities using the NEPSY-II subtests 1 year later (4.5 years). We applied a data-driven linked independent component analysis (ICA) to explore the shared co-variation of gray and white matter measures. Two sources of structural variation at 3.5 years of age demonstrated relationships with Speeded Naming scores at 4.5 years of age. The first imaging component involved volumetric variability in reading-related cortical regions alongside microstructural features of the superior longitudinal fasciculus (SLF). The second component was dominated by cortical volumetric variations within the cerebellum and visual association area. In a subset of children with rs-fMRI data, we evaluated the inter-network functional connectivity of the left-lateralized fronto-parietal language network (FPL) and its relationship with pre-reading measures. Higher functional connectivity between the FPL and the default mode and visual networks at 3.5 years significantly predicted better Phonological Processing scores at 4.5 years. Together, these results suggest that the integration of functional networks, as well as the co-development of white and gray matter brain structures in early childhood, support the emergence of pre-reading measures in preschool children.
Introduction
Reading is an essential skill that plays a fundamental role in academic achievement, social engagement with the world and peers, and mental health (Dehaene-Lambertz et al., 2006a; Torppa et al., 2006; May et al., 2011; Lonigan et al., 2013). Language skills develop dramatically throughout early childhood and lay the foundation for later reading skills. Phonological awareness and speeded naming in particular are strong predictors of future reading success (Lonigan et al., 2009; Snowling and Hulme, 2011). These skills are supported by neurobiological processes that are primed in utero (May et al., 2011) and develop rapidly during early life. A thorough understanding of the early regionally specific structural and functional brain features associated with later reading abilities is essential to determine how reading difficulties emerge, and to help appropriately target early language and reading interventions (Zijlstra et al., 2020).
Reading involves a left lateralized network of white matter connections that support communication between cortical regions (Torppa et al., 2006; Schlaggar and McCandliss, 2007; Zijlstra et al., 2020). Frontal and temporal-parietal cortical regions are connected via dorsal white matter pathways including the arcuate fasciculus and superior longitudinal fasciculus (SLF). Frontal, temporal, and occipital regions are connected by ventral white matter pathways [inferior longitudinal (ILF), inferior fronto-occipital (IFO), and uncinate fasciculi (UF)] (Klingberg et al., 2000; Brown et al., 2001; Hickok and Poeppel, 2004; Welcome et al., 2011). Throughout childhood, both cortical structure (Szaflarski et al., 2006; Holland et al., 2007) and brain functional patterns (Dehaene-Lambertz et al., 2002, 2006b; Peña et al., 2003) become more left-lateralized. White matter structures, such as the arcuate fasciculus, also demonstrate left lateralization in early life (O’Muircheartaigh et al., 2013; Dubois et al., 2016), and are linked to language skill development in young children (Reynolds et al., 2019b).
Relationships between reading performance and cortical brain structure and development have been observed in adults (Klingberg et al., 2000; Brown et al., 2001; Welcome et al., 2011), adolescents (Kronbichler et al., 2008; Lebel et al., 2013), and school-aged children (Beaulieu et al., 2005; Deutsch et al., 2005; Eckert et al., 2005; Niogi and McCandliss, 2006; Lebel and Beaulieu, 2009). Cortical structure in children and adolescents is associated with phonological processing skills (Lu et al., 2007), with higher baseline reading skills associated with faster changes in gray matter volume in typically developing children aged 5–15 years (Houston et al., 2014; Linkersdörfer et al., 2014). Furthermore, children aged 6–7 years who subsequently received a diagnosis of dyslexia had reduced cortical thickness in left hemisphere reading-related regions compared to children who did not eventually receive a diagnosis (Clark et al., 2014). There is evidence that these alternate development patterns emerge even earlier, with slower proportional growth of left cortical regions was observed in preschool aged children with a family history of reading disorders (Ostertag et al., 2021).
White matter microstructural properties like myelin and axonal packing (Geeraert et al., 2020) support pre-reading skills in preschool children (Walton et al., 2018), and reading skills in adolescents (Ben-Shachar et al., 2007) and adults (Welcome and Joanisse, 2014). Longitudinal studies have also demonstrated that white matter development supports reading gains. Faster changes in white matter fractional anisotropy (FA) and mean diffusivity (MD), that represent faster maturation, are associated with larger gains in reading skills in children with and without developmental disorders (Yeatman et al., 2012a; Treit et al., 2013; Wang et al., 2017). Furthermore, similar changes in white matter microstructure have been demonstrated following intensive reading interventions (Keller and Just, 2009; Huber et al., 2018). While these studies suggest associations between early white matter maturation and reading development in children, it is unclear what early white matter characteristics support pre-reading skills later in childhood. Emerging evidence suggests that arcuate fasciculus and corticospinal tract microstructural properties even in infancy can predict phonological processing and vocabulary at 5 years of age (Zuk et al., 2021). However, the coinciding cortical structures and white matter neurobiological properties that support pre-reading skills in preschool aged children, before formal reading education, are not well-understood.
Activation and functional connectivity with reading related regions are related to language skill gains (Xiao et al., 2016), later reading outcomes (Jasińska et al., 2020), and brain development during childhood (Saygin et al., 2016). The process of learning to read is supported by refined functional architecture (Dehaene et al., 2010; Cross et al., 2021). In general, children with higher reading proficiency show faster development of reading-related brain regions than children with poor reading skills (Yeatman et al., 2012b; Wang et al., 2017; Lebel et al., 2019; Reynolds et al., 2019b). A recent cross-sectional study (Benischek et al., 2020) identified higher functional connectivity within the reading network (Wernicke’s and Broca’s temporal-parietal areas), but more negative connectivity between reading areas and the default mode network (DMN), in young children with better pre-reading skills. The brain’s functional architecture during the preschool period has also been shown to relate to pre-reading skills (Raschle et al., 2012; Saygin et al., 2013; Walton et al., 2018; Reynolds et al., 2019b). Furthermore, pre-readers with a family history of dyslexia show reduced temporal-parietal and temporal-occipital activation during sound matching compared to pre-readers without a family history of dyslexia (Raschle et al., 2012).
White matter structure, cortical volume, and functional connectivity have all been studied separately with respect to reading and pre-reading, but their development is inherently linked. In this study, we aimed to determine how early childhood brain morphometry, white matter microstructure, and functional brain network communication predict measures of pre-reading skills later in childhood. Investigating all measures simultaneously can be challenging, as it increases the number of comparisons and typical univariate analyses do not model the shared and spatially specific development between features. Our data-driven linked independent component analysis (ICA) technique (Llera et al., 2019) combined multiple parameters simultaneously to measure co-variation of white matter and cortical structures across participants as they relate to pre-reading measures. This approach has been used to study neurobiological development and aging as well as childhood disorders like autism spectrum disorder and ADHD (Itahashi et al., 2015; Wolfers et al., 2017). In a subset of participants with high-quality functional MRI data, we quantified how the fronto-parietal language network (FPL) communicates with other networks in the brain that may play a role in pre-reading ability. These included the cerebellar (Alvarez and Fiez, 2018), visual, and DMNs which support memory, visual function, and attention while reading (Smallwood et al., 2013), respectively. We hypothesized that multimodal brain structural properties and inter-network connectivity at 3 years of age would relate to pre-reading measures assessed 1 year later.
Materials and methods
Participants
The University of Calgary Conjoint Health Research Ethics Board (CHREB) approved this study (REB13-0020). Informed written consent was obtained from each participant’s legal guardian prior to the commencement of the study, and ongoing verbal assent was obtained from the participants. Participants were 35 children (19 boys/16 girls) selected from the ongoing Calgary Preschool MRI Study (Reynolds et al., 2020) based on the criteria of having high-quality T1-weighted and diffusion weighted scans at 3.5 years (3.49 ± 0.14; range 3.25–3.75 years) and a pre-reading assessment at age 4.5 years (4.50 ± 0.16; range 4.25–4.75 years). The University of Calgary CHREB approved this study (REB13-0020). Informed written consent was obtained from each participant’s legal guardian prior to the commencement of the study, and ongoing verbal assent was obtained from the participants.
Language assessments
Children’s pre-reading skills were assessed using the NEPSY-II Speeded Naming and Phonological Processing subtests (∼20 min) at 3.5 and 4.5 years of age. The Speeded Naming subtest assesses rapid semantic access to and production of names of colors and shapes, and the Phonological Processing subtest assesses phonemic awareness (Korkman et al., 2007). Age standardized Speeded Naming Combined Scaled Scores (accounts for both speed and accuracy) and Phonological Processing Scaled Scores were calculated and used in the analysis. On these measures, higher scores are indicative of better performance.
Magnetic resonance imaging acquisition
All imaging was conducted using the same General Electric 3T MR750w system and a 32-channel head coil (GE, Waukesha, WI) at the Alberta Children’s Hospital in Calgary, Canada. Children were scanned either while awake and watching a movie of their choice, or while sleeping without sedation. fMRI scans during which the child was asleep were excluded from analyses. This approach has been used in our prior work and is effective for scanning young children while awake (Long et al., 2017). Prior to scanning, parents were provided with detailed information on MRI procedures and given the option to complete a practice MRI session in a training scanner to familiarize the child with the scanning environment, or to make use of a take home pack with this information (e.g., noise recordings; Thieba et al., 2018). Families were also provided with a book that incorporates our scanning procedures into an engaging story (Frayne, 2015) and we encouraged the parents/guardians to review the materials with the child.
T1-weighted images were acquired using a FSPGR BRAVO sequence, 210 axial slices; 0.9 × 0.9 × 0.9 mm resolution, TR = 8.23 ms, TE = 3.76 ms, flip angle = 12°, matrix size = 512 × 512, inversion time = 540 ms. Whole-brain diffusion weighted images were acquired using single shot spin echo echo-planar imaging sequence: 1.6 × 1.6 × 2.2 mm resolution (resampled on scanner to 0.78 × 0.78 × 2.2 mm), full brain coverage, TR = 6,750 ms; TE = 79 ms (set to minimum for first year), 30 gradient encoding directions at b = 750 s/mm2, and five interleaved images without gradient encoding at b = 0 s/mm2 for a total acquisition time of approximately 4 min. Passive viewing fMRI data were acquired while children were watching a movie of their choice, with a gradient-echo echo-planar imaging (EPI) sequence: total sequence time = 8 min and 10 s, 36 axial slices, 3.59 × 3.59 × 3.6 mm resolution, TR = 2,000 ms, TE = 30 ms, flip angle = 60°, matrix size = 64 × 64, 250 volumes.
Anatomical imaging
Voxel-based morphometry (VBM) processing was undertaken on the T1-weighted images using FSL [FMRIB (Functional Magnetic Resonance Imaging of the Brain) Software Library freely available at fsl.fmrib.ox.ac.uk] to create gray matter maps. N4-bias correction (Tustison et al., 2010) was performed using ANTs (Avants et al., 2009), then images were transformed to radiological orientation. Brain extraction was undertaken using Brain Extraction Tool (BET) via the standard FSL VBM protocol; when brain extraction using the default settings was unsuccessful (n = 7), BET was performed manually to achieve the best results. Next, images were segmented into gray matter, white matter and CSF. The gray matter images were then affine registered to the NIHPD asymmetrical pediatric brain template (4.5–8.5 years template; this template was used to be consistent with our prior work in the larger Calgary Preschool MRI sample and to permit future analysis spanning 2–8 years) in Montreal Neurological Institute (MNI) standard space (Fonov et al., 2011), resampled to 2 mm isotropic voxels, concatenated and averaged. Using the standard FSL pipeline, these average gray matter images were used to create a study-specific non-linear gray matter template (2 mm isotropic voxels), following which gray matter images were non-linearly registered to the study-specific template, modulated, smoothed (2 mm), and concatenated into one 4D file.
Diffusion imaging
Raw diffusion images were visually quality checked and all motion-corrupted volumes, or volumes with artifacts were removed prior to processing. Datasets passed final quality assurance checks (n = 35) if they had at least 18 high quality diffusion weighted volumes, and two high quality b0 volumes remaining following volume removal. Data was then pipelined through ExploreDTI V4.8.6 (Leemans et al., 2009) to correct for signal drift, Gibbs ringing (non-DWIs), subject motion, and eddy current distortions. FA, axial (AD), and radial diffusivity (RD) image maps were extracted. All measures were included to understand the specific white matter microstructural properties that contribute to pre-reading development. For images where brain extraction during diffusion tensor imaging (DTI) preprocessing did not remove all non-brain material, FSL BET was run on the extracted FA map, and the resulting binary mask was used to mask the FA, AD, and RD images. These DTI maps were then non-linearly warped using ANTs (Avants et al., 2009) to the NIHPD asymmetrical pediatric brain template (ages 4.5–8.5 years) in MNI standard space (Fonov et al., 2011). All of the registered diffusion data was then merged into one four-dimensional image to create a mean FA mask and image for all subjects. The mean FA image was skeletonized with a threshold of FA > 0.2 to create a mean FA skeleton mask. All participants’ FA, AD, and RD images were non-linearly projected onto that skeleton. Analysis was conducted on skeletonized maps of diffusion measures within tracts that have been associated with reading in children (Wandell and Yeatman, 2013): the bilateral uncinate fasciculus (UF), inferior longitudinal fasciculus (ILF), IFO, and SLF. Tracts were identified using the JHU white matter tractography atlas.
Linked independent component analysis
Gray matter morphometry and skeletonized FA, AD and RD maps were aligned to the same space, concatenated, and then used as input for a linked ICA using FSL tools (Llera et al., 2019). This is a data-driven approach that uncovers the neurobiological variations across multiple MRI parameter maps (Groves et al., 2011) so by including multiple diffusion measures we aim to explain the neurobiology more specifically. Briefly, the concatenated imaging data from all subjects are decomposed into a series of linked independent components. Linked components involve a spatial map for each MRI parameter (Figure 1D) that displays where that combination of MRI metrics co-vary across participants according to a common component weighting (Figure 1C). Each MRI parameter has a relative modality loading in each component (Groves et al., 2012) that contributes a fraction of the total variance explained by that component. The component weight describes the combined relative variability of multiple MRI metrics across participants and decreases the number of multiple comparisons while enabling examination of multiple parameters in conjunction (Figure 1). Functional data could not be incorporated in our linked models because only a subset of participants’ functional data passed quality check procedures.
Figure 1. Linked independent component analysis pipeline. (A) Preprocessed T1 and diffusion weighted images (DWI) were processed using voxel-based morphometry (VBM) and tract-based spatial statistics (TBSS) pipelines (with the reading-related tracts mask shown in yellow) to create (B) individual VBM, fractional anisotropy (FA), axial diffusivity (AD), and radial diffusivity (RD) spatial maps for each subject as input for the linked independent component analysis (ICA). The data is decomposed into a series of linked components comprised of (C) component weights that reflect the inter-subject variation shared across (D) all input measures (VBM, FA, AD, and RD) in specific brain regions (for example, shown in red).
Resting state functional magnetic resonance imaging
Resting state functional MRI preprocessing is described in detail elsewhere (Long et al., 2019). Briefly, processing was undertaken using FSL (slice timing, head motion correction, T1 image segmentation, head motion outlier detection, co-registration, and spatial normalization and smoothing) (Jenkinson et al., 2002) and AFNI v17.3.03 (regression of the nuisance signals, band-pass filtering and linear trend removal) (Cox, 1996). fMRI data was registered to a pediatric template in MNI standard space (Fonov et al., 2011). Participants with fMRI data that had < 4 min of low-motion data, or excessive motion (relative frame-wise displacement < 0.3 mm) at any time (Satterthwaite et al., 2013) were excluded from the analysis, resulting in a total of 21 fMRI datasets.
Clean fMRI data (n = 21) was temporally concatenated and analyzed using ICA to decompose the data into 20 components or networks. This model order was chosen to adequately model known neural networks, physiological signals, and noise (Wang and Li, 2015). To identify resting state networks (RSNs), the components were compared (using cross correlation) to the 10 most common RSNs (Smith et al., 2009). We focused on the left-lateralized FPL network because it involves regions, including Broca’s and Wernicke’s areas, that are functionally active during cognitive-language tasks (Smith et al., 2009; Reynolds et al., 2019b). Dual regression algorithms were used to back-reconstruct subject-specific RSNs composed of data-driven voxel-wise clusters. Reading is a complex skill that involves communication between distinct networks in the brain that are involved with visual, memory, attention, and language processing. Therefore, we focused on how the FPL network communicates with other relevant brain networks, as recent findings in preschool aged children suggest that both within-network communication as well as integration with other networks support pre-reading skills (Benischek et al., 2020). We investigated the inter-network communication between the FPL and networks responsible for visual processing, memory, and attention. The average timeseries from each of these networks (FPL, visual, DMN, and the cerebellar RSN) was extracted. Inter-network functional connectivity between the FPL network and the three RSNs were calculated for each participant.
Statistical analysis
This exploratory study used partial correlations controlling for sex. We performed bootstrapping in SPSS (version 26) with 1,000 samples to determine confidence interval of any relationships between the linked ICA component weights at 3.5 years and pre-reading measures at 4.5 years. Sex was included as a covariate because of known sex differences in brain development in this age range (Reynolds et al., 2019a). The inter-network functional connectivity between the FPL and the three RSNs was also tested for correlations with pre-reading measures 1 year later using a partial correlation analysis controlling for sex and age at the time of the scan. As this was an exploratory study, the nominal p-values for the relationships between component weights, inter-network functional connectivity and pre-reading measure are reported. We also report if these relationships survive false-discovery rate (FDR) correction for the linked ICA and functional RSNs analyses separately.
Results
Pre-reading measures
At 4.5-years, the mean Phonological Processing standard score was 11.9 ± 2.1 and ranged from 7 to 16. The mean Speeded Naming standard score was 12.4 ± 2.3 and ranged from 7 to 18. Speeded Naming and Phonological Processing scores at 3.5 years of age were not correlated with scores obtained at 4.5 years (Speeded Naming: r = 0.02, p = 0.91; Phonological Processing: r = 0.26, p = 0.14).
Linked brain structure features that relate to pre-reading
The masked (bilateral UF, IFO/ILF, and SLF) and skeletonized DTI maps (FA, AD, RD) and the gray matter morphometry maps were input into the linked ICA. The linked ICA generated 14 components. Of these, three components were driven by an individual subject; therefore, these components were not included in the correlation analysis that investigated the relationships between the brain imaging component weights and pre-reading measures. Of the remaining 11 components, two correlated with pre-reading measures. To confirm the validity of components of interest, linked independent component models with 13 and 15 components (model order) were also assessed. The two components described below were reproducible using all three model orders (p < 0.05). Intracranial volume was not related to either reading score (p > 0.05).
Component 8 weights were positively associated with Speeded Naming scores (r = 0.37, p = 0.03, CI: [0.04, 0.64]) at age 4.5 years (Figure 2), but this relationship did not survive FDR correction (corrected p = 0.3). This component involved volumetric variability in the angular gyrus, thalamus, fusiform gyrus, middle temporal regions, and the dorsolateral prefrontal cortex (Table 1) and was linked with bilateral changes in FA and AD along the SLF. Specifically, relatively higher component weights amongst participants represent larger volumes in the sensorimotor regions, the angular gyrus, thalamus, fusiform gyrus, and middle temporal area, smaller volumes in the posterior cingulate and visual association areas and the prefrontal cortex, combined with higher FA and AD, and lower RD (though less pronounced than FA and AD) in the bilateral SLF were related to higher pre-reading scores.
Figure 2. Component 8 was approximately evenly weighed by the voxel-based morphometry (VBM) gray matter, fractional anisotropy (FA), radial diffusivity (RD), and axial diffusivity (AD) contributions. This component was positively correlated with Speeded Naming at 4.5 years. For regions in red (positive clusters), children with larger volumes/higher DTI metrics performed better on Speeded Naming, and for regions in blue (negative clusters), children with smaller volumes/DTI metrics performed better on Speeded Naming. For reference, the superior longitudinal fasciculus is also shown.
Table 1. Characterization of component 8 voxel-based morphometry gray matter clusters.
A second component (Component 10) was dominated by volumetric variability in the cerebellum, precuneus, angular gyri, and areas in the occipital cortex (visual regions), and was positively associated with Speeded Naming (r = 0.35, p = 0.04, CI: [0.06, 0.61]) at age 4.5 years (Figure 3). This association did not survive FDR correction (corrected p = 0.2). Smaller cortical volumes in the cerebellum, and supramarginal and angular gyri, and relatively larger volumes in visual areas in the occipital lobe and the fusiform in participants at 3.5 years (Table 2) were associated with higher pre-reading scores at 4.5 years.
Figure 3. Component 10 was dominated by gray matter contributions. This component was positively correlated with Speeded Naming performance at 4.5 years. For regions in red (positive clusters), children with larger volumes performed better on Speeded Naming, and for regions in blue (negative clusters), children with smaller volumes performed better on Speeded Naming.
Table 2. Characterization of component 10 voxel-based morphometry gray matter clusters.
Resting state networks and pre-reading
We examined average inter-network functional connectivity between the FPL RSN (Figure 4) and three other networks including the DMN, the occipital pole visual and the cerebellar RSNs. The inter-network functional connectivity (correlation between average RSN time series) between the FPL and the DMN assessed at 3.5 years of age was correlated with Phonological Processing scores 1year later at 4.5 years of age (r = 0.62, p = 0.006, CI: [0.25, 0.86]) and this relationship survived FDR correction (corrected p = 0.02). Inter-network functional connectivity between the occipital pole visual RSN and the FPL at 3.5 years of age was associated with Phonological Processing scores one year later (r = 0.51, p = 0.03, CI: [0.001, 0.85]), but this did not survive FDR correction (corrected p = 0.05). Inter-network functional connectivity was not significantly related to either structural linked component weights.
Figure 4. The group average left-lateralized fronto-parietal resting state network (left) and the linear relationships between inter-network functional connectivity at 3.5 years predicting Phonological Processing scores at 4.5 years of age.
Discussion
In this exploratory study of typically developing young children, we identified multiparametric brain imaging features at 3.5 years of age that were associated with pre-reading measures 1 year later. Variations in brain volume (angular gyrus, lingual gyrus, thalamus, fusiform gyrus, middle temporal regions, cerebellum, precuneus, and the dorsolateral prefrontal cortex), SLF microstructure, and inter-network functional connectivity of the FPL network at 3.5 years were associated with better pre-reading skills at 4.5 years. More mature SLF microstructure coincided with cortical morphometry variations in portions of Wernicke’s area and Broca’s area. Both are considered critical language comprehension and production regions and were related to better pre-reading scores 1 year later. Also, higher functional connectivity between the FPL and the DMN and visual networks significantly predicted better pre-reading measures 1 year later. Together these findings suggest that gray and white matter structure, as well as inter-network functional communication, contribute to the development of pre-reading skills in young children.
Our multimodal approach was able to capture co-occurring neurobiological variations (Wolfers et al., 2017; Llera et al., 2019) in gray and white matter structures. Shared inter-subject variation (Component 8) of gray matter volume in reading-related regions (e.g., supramarginal and angular gyrus and dorsolateral prefrontal cortex) co-varied with microstructural measures (FA, AD, and a smaller contribution of RD) of the white matter that connects them (SLF). Reading-related brain areas such as the SLF and the cortical regions identified continue to develop throughout childhood and even into adulthood (Lebel et al., 2008) as the brain refines and integrates inter-network communications. The SLF is a key dorsal pathway involved in auditory-motor integration and the phonological aspects of both reading and speech (Vigneau et al., 2006). The associations between the SLF and pre-reading measures is consistent with prior longitudinal diffusion MRI research that has reported relationships between arcuate fasciculus (a component of the SLF) microstructure and reading skills in older children (Wandell and Yeatman, 2013). Furthermore, arcuate fasciculus microstructure during infancy is associated with phonological awareness skills and vocabulary knowledge in kindergarten (Zuk et al., 2021). Our findings are also consistent with prior research that has reported cross-section and longitudinal associations between cortical thickness and volume in these areas and reading abilities in older children. Higher FA and lower MD suggest that participants with increased myelin and/or more tightly packed axons (Beaulieu, 2002; Song et al., 2002, 2003). Because FA generally increases and MD decreases during typical development (Reynolds et al., 2019b), the profiles observed here suggest a more mature SLF in children with better pre-reading skills at age 4.5 years. Congruent cerebral structural development involves both thinning and thickening of the cortex depending on the region and stage of development (Remer et al., 2017). These results suggest that co-occurring white and gray matter structure in early childhood may lay the foundation for more advanced reading abilities throughout childhood and adolescence.
The second component associated with later pre-reading performance was dominated by gray matter volume variations. Smaller volumes in the cerebellum, frontal and sensory regions, and larger angular and visual gyri cortical volumes at 3.5 years of age were related to better Speeded Naming scores 1 year later. The cerebellum plays a critical role in reading development through both the dorsal and ventral circuits that support phonological and semantic processes (Alvarez and Fiez, 2018; Benischek et al., 2020). The inclusion of dorsal and ventral cerebellar gray matter regions in this component suggests that the cerebro-cerebellar pathway plays a key role in the development of pre-reading skills. The angular gyrus also plays a role in language and reading, including orthographic processing, with increased functional interactions with Broca’s and Wernicke’s areas and the visual word form area during reading (Segal and Petrides, 2013; Benischek et al., 2020). Better language ability is related to increased activation in this region (Van Ettinger-Veenstra et al., 2016), and disruptions to functional connectivity with temporal and occipital language regions have been reported in adults with dyslexia compared to controls (Horwitz et al., 1997). The visual system is also important for reading, with sensitivity of the visual cortex for word visibility increasing throughout childhood (Ben-Shachar et al., 2011). Connectivity differences between visual areas and other language and reading regions have been observed in poor compared to good readers (Wandell et al., 2012) suggesting that earlier refinement of the cerebellum and somatosensory cortical regions, combined with growth of the visual cortex, relate to better pre-reading skills in the preschool period.
Higher inter-network functional connectivity between the left-lateralized FPL network and both the DMN (including bilateral angular gyri) and visual RSNs at 3.5 years of age predicted better Phonological Processing scores 1 year later. Prior functional connectivity studies of reading have reported mixed results, but have generally found increased activation and connectivity amongst reading-related, sensory and motor regions (Weiss-Croft and Baldeweg, 2015). However, across development regions that are typically active during reading tasks do not exclusively communicate amongst each other, but instead tend to integrate with other networks including the DMN (Smallwood et al., 2013; Vogel et al., 2013). Our results reflect a similar pattern at a much earlier stage of development, where higher functional connectivity between the FPL and DMN was significantly predictive of Phonological Processing scores in early childhood. This suggests that even before formal reading education, better integration of functional networks support pre-reading ability and may lay the foundation for future reading capability as well.
In most individuals, language and reading are left-lateralized in the brain, with resting state functional connectivity patterns demonstrating increased lateralization throughout early development (Reynolds et al., 2019b; Benischek et al., 2020). Functional connectivity lateralization have been shown to be related to the development of language skills task-based fMRI studies (Szaflarski et al., 2006; Holland et al., 2007; Perani et al., 2011; Yamada et al., 2011; Xiao et al., 2016), where lateralization patterns tend to increase in the dominant hemisphere until a plateau in early adulthood that then gradually decreases in laterality over the life span (Szaflarski et al., 2006). The structural brain properties associated with better speeded naming scores were bilateral in our study, which may reflect that while lateralization strengthens over time, a more broad network of regions support pre-reading skills at this early stage of development (Qiu et al., 2011). Our findings suggest that more inter-connected functional network architecture at this stage of development may lay the foundation for later reading abilities, which require seamless communication between brain networks. Functional brain signals are shared among specific cortical regions and networks, and maturation and refinement of the underlying microstructure of pathways connecting those areas may directly support these functional pathways. Children who utilize those particular pathways through frequent exposure to language and reading in early childhood may establish both the structural and functional brain foundations to support reading abilities later in life.
Limitations
This longitudinal study involved a group of typically developing preschool children with average to high pre-reading skills. Future studies including children with lower scores may assist in determining whether similar patterns hold in children more likely to develop reading problems. We had a tight age range in this study, but the sample size in this exploratory study was relatively small and the statistical tests were under-powered. While VBM has shown similar volumetric results compared to manual region of interest (Asami et al., 2012), other factors like gyrification, surface area and cortical thickness can influence results. Functional connectivity between the FPL and DMN significantly predicted pre-reading measures, specifically phonological processing; however, the structural brain measures associated with pre-reading did not survive a statistical correction. Further studies with larger sample sizes are necessary to confirm these relationships between early brain measures and their influence on a child’s pre-reading abilities later in life.
Conclusion
In this study we found that linked development of brain white and gray matter structure in early childhood was associated with pre-reading measures 1 year later. In particular, the SLF and the cortical regions that it connects, as well as cerebellar-cerebral reading-related circuits, appear to be important for the development of children’s pre-reading skills. Furthermore, a more functionally integrated FPL network predicted better pre-reading skills 1 year later. This analysis approach demonstrates that co-development of white and gray matter brain structures in early life, as well as the integration of functional networks before formal reading education, are associated with pre-reading abilities in preschool children.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: The datasets analyzed for this study can be found in the Calgary Preschool MRI study available through the Open Science Framework (https://osf.io/axz5r/).
Ethics statement
The studies involving human participants were reviewed and approved by the University of Calgary Conjoint Health Research Ethics Board (CHREB). Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author contributions
KM, JR, and XL contributed to data acquisition, analysis, and interpretation. AL contributed to data analysis and manuscript revisions. DD and CL were responsible for study conception, funding, data interpretation, and manuscript revisions. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Canadian Institutes of Health Research (CIHR) (funding reference numbers: IHD-134090 and MOP-136797, New Investigator Award to CL). CL received funding from the Canada Research Chair Program. KM was supported by the NSERC Postdoctoral Fellowship and the T. Chen Fong Postdoctoral Fellowship in Medical Imaging Science. JR was supported by an Eyes High University of Calgary Postdoctoral Scholarship, the T. Chen Fong Postdoctoral Fellowship in Medical Imaging Science, and a CIHR Postdoctoral Fellowship (MFE-164703). AL was supported by the Horizon 2020 Programme CANDY (Grant no. 847818).
Acknowledgments
We thank members of the APrON study for assistance with participant recruitment. A version of this manuscript is available on bioRxiv as a preprint (Manning et al., 2021).
Conflict of interest
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Alvarez, T. A., and Fiez, J. A. (2018). Current perspectives on the cerebellum and reading development. Neurosci. Biobehav. Rev. 92, 55–66. doi: 10.1016/j.neubiorev.2018.05.006
Asami, T., Bouix, S., Whitford, T. J., Shenton, M. E., Salisbury, D. F., and McCarley, R. W. (2012). Longitudinal loss of gray matter volume in patients with first-episode schizophrenia: DARTEL automated analysis and ROI validation. Neuroimage 59, 986–996. doi: 10.1016/J.NEUROIMAGE.2011.08.066
Avants, B., Tustison, N., and Johnson, H. (2009). Advanced Normalization Tools (ANTS). Insight J. 1, 1–35.
Beaulieu, C. (2002). The basis of anisotropic water diffusion in the nervous system - A technical review. NMR Biomed. 15, 435–455. doi: 10.1002/nbm.782
Beaulieu, C., Plewes, C., Paulson, L. A., Roy, D., Snook, L., Concha, L., et al. (2005). Imaging brain connectivity in children with diverse reading ability. Neuroimage 25, 1266–1271. doi: 10.1016/j.neuroimage.2004.12.053
Benischek, A., Long, X., Rohr, C. S., Bray, S., Dewey, D., and Lebel, C. (2020). Pre-reading language abilities and the brain’s functional reading network in young children. Neuroimage 217:116903. doi: 10.1016/j.neuroimage.2020.116903
Ben-Shachar, M., Dougherty, R. F., and Wandell, B. A. (2007). White matter pathways in reading. Curr. Opin. Neurobiol. 17, 258–270. doi: 10.1016/j.conb.2007.03.006
Ben-Shachar, M., Dougherty, R. F., Deutsch, G. K., and Wandell, B. A. (2011). The development of cortical sensitivity to visual word forms. J. Cogn. Neurosci. 23, 2387–2399. doi: 10.1162/jocn.2011.21615
Brown, W. E., Eliez, S., Menon, V., Rumsey, J. M., White, C. D., and Reiss, A. L. (2001). Preliminary evidence of widespread morphological variations of the brain in dyslexia. Neurology 56, 781–783. doi: 10.1212/WNL.56.6.781
Clark, K. A., Helland, T., Specht, K., Narr, K. L., Manis, F. R., Toga, A. W., et al. (2014). Neuroanatomical precursors of dyslexia identified from pre-reading through to age 11. Brain 137, 3136–3141. doi: 10.1093/brain/awu229
Cox, R. W. (1996). AFNI: Software for analysis and visualization of functional magnetic resonance neuroimages. Comput. Biomed. Res. 29, 162–173. doi: 10.1006/cbmr.1996.0014
Cross, A. M., Ramdajal, R., Peters, L., Vandermeer, M. R. J., Hayden, E. P., Frijters, J. C., et al. (2021). Resting-state functional connectivity and reading subskills in children. Neuroimage 243, 118529. doi: 10.1016/j.neuroimage.2021.118529
Dehaene, S., Pegado, F., Braga, L. W., Ventura, P., Nunes Filho, G., Jobert, A., et al. (2010). How learning to read changes the cortical networks for vision and language. Science 330, 1359–1364. doi: 10.1126/science.1194140
Dehaene-Lambertz, G., Dehaene, S., and Hertz-Pannier, L. (2002). Functional neuroimaging of speech perception in infants. Science 298, 2013–2015. doi: 10.1126/science.1077066
Dehaene-Lambertz, G., Hertz-Pannier, L., and Dubois, J. (2006a). Nature and nurture in language acquisition: Anatomical and functional brain-imaging studies in infants. Trends Neurosci. 29, 367–373. doi: 10.1016/j.tins.2006.05.011
Dehaene-Lambertz, G., Hertz-Pannier, L., Dubois, J., Mériaux, S., Roche, A., Sigman, M., et al. (2006b). Functional organization of perisylvian activation during presentation of sentences in preverbal infants. Proc. Natl. Acad. Sci. U.S.A. 103, 14240–14245. doi: 10.1073/pnas.0606302103
Deutsch, G. K., Dougherty, R. F., Bammer, R., Siok, W. T., Gabrieli, J. D. E., and Wandell, B. (2005). Children’s reading performance is correlated with white matter structure measured by diffusion tensor imaging. Cortex 41, 354–363. doi: 10.1016/S0010-9452(08)70272-7
Dubois, J., Poupon, C., Thirion, B., Simonnet, H., Kulikova, S., Leroy, F., et al. (2016). Exploring the Early Organization and Maturation of Linguistic Pathways in the Human Infant Brain. Cereb. Cortex 26, 2283–2298. doi: 10.1093/cercor/bhv082
Eckert, M. A., Leonard, C. M., Wilke, M., Eckert, M., Richards, T., Richards, A., et al. (2005). Anatomical signatures of dyslexia in children: Unique information from manual and voxel based morphometry brain measures. Cortex 41, 304–315. doi: 10.1016/S0010-9452(08)70268-5
Fonov, V., Evans, A. C., Botteron, K., Almli, C. R., McKinstry, R. C., and Collins, D. L. (2011). Unbiased average age-appropriate atlases for pediatric studies. Neuroimage 54, 313–327. doi: 10.1016/j.neuroimage.2010.07.033
Frayne, A. (2015). Pluto and the MRI Rocket Ship Adventure. Available online at: https://www.scribd.com/document/351490812/Pluto-and-the-Mri-Rocket-Ship-Adventure (accessed March 15, 2021).
Geeraert, B. L., Chamberland, M., Marc Lebel, R., and Lebel, C. (2020). Multimodal principal component analysis to identify major features of white matter structure and links to reading. PLoS One 15:e0233244. doi: 10.1371/journal.pone.0233244
Groves, A. R., Beckmann, C. F., Smith, S. M., and Woolrich, M. W. (2011). Linked independent component analysis for multimodal data fusion. Neuroimage 54, 2198–2217. doi: 10.1016/j.neuroimage.2010.09.073
Groves, A. R., Smith, S. M., Fjell, A. M., Tamnes, C. K., Walhovd, K. B., Douaud, G., et al. (2012). Benefits of multi-modal fusion analysis on a large-scale dataset: Life-span patterns of inter-subject variability in cortical morphometry and white matter microstructure. Neuroimage 63, 365–380. doi: 10.1016/j.neuroimage.2012.06.038
Hickok, G., and Poeppel, D. (2004). Dorsal and ventral streams: A framework for understanding aspects of the functional anatomy of language. Cognition 92, 67–99. doi: 10.1016/j.cognition.2003.10.011
Holland, S. K., Vannest, J., Mecoli, M., Jacola, L. M., Tillema, J. M., Karunanayaka, P. R., et al. (2007). Functional MRI of language lateralization during development in children. Int. J. Audiol. 46, 533–551. doi: 10.1080/14992020701448994
Horwitz, B., Rumsey, J. M., and Donohue, B. C. (1997). Functional connectivity of the angular gyrus during single word reading in normal and dyslexic men. Neuroimage 5, 8939–8944.
Houston, S. M., Lebel, C., Katzir, T., Manis, F. R., Kan, E., Rodriguez, G. G., et al. (2014). Reading skill and structural brain development. Neuroreport 25, 347–352. doi: 10.1097/WNR.0000000000000121
Huber, E., Donnelly, P. M., Rokem, A., and Yeatman, J. D. (2018). Rapid and widespread white matter plasticity during an intensive reading intervention. Nat. Commun. 9, 1–13. doi: 10.1038/s41467-018-04627-5
Itahashi, T., Yamada, T., Nakamura, M., Watanabe, H., Yamagata, B., Jimbo, D., et al. (2015). Linked alterations in gray and white matter morphology in adults with high-functioning autism spectrum disorder: A multimodal brain imaging study. NeuroImage Clin. 7, 155–169. doi: 10.1016/j.nicl.2014.11.019
Jasińska, K. K., Shuai, L., Lau, A. N. L., Frost, S., Landi, N., and Pugh, K. R. (2020). Functional connectivity in the developing language network in 4-year-old children predicts future reading ability. Dev. Sci. 24:e13041. doi: 10.1111/desc.13041
Jenkinson, M., Bannister, P., Brady, M., and Smith, S. (2002). Improved optimization for the robust and accurate linear registration and motion correction of brain images. Neuroimage 17, 825–841. doi: 10.1016/S1053-8119(02)91132-8
Keller, T. A., and Just, M. A. (2009). Altering Cortical Connectivity: Remediation-Induced Changes in the White Matter of Poor Readers. Neuron 64, 624–631. doi: 10.1016/j.neuron.2009.10.018
Klingberg, T., Hedehus, M., Temple, E., Salz, T., Gabrieli, J. D. E., Moseley, M. E., et al. (2000). Microstructure of temporo-parietal white matter as a basis for reading ability: Evidence from diffusion tensor magnetic resonance imaging. Neuron 25, 493–500. doi: 10.1016/S0896-6273(00)80911-3
Korkman, M., Kirk, U., and Kemp, S. (2007). NEPSY-Second Edition (NEPSY-II). J. Psychoeduc. Assess. 28, 175–182.
Kronbichler, M., Wimmer, H., Staffen, W., Hutzier, F., Mair, A., and Ladurner, G. (2008). Developmental dyslexia: Gray matter abnormalities in the occipitotemporal cortex. Hum. Brain Mapp. 29, 613–625. doi: 10.1002/hbm.20425
Lebel, C., and Beaulieu, C. (2009). Lateralization of the arcuate fasciculus from childhood to adulthood and its relation to cognitive abilities in children. Hum. Brain Mapp. 30, 3563–3573. doi: 10.1002/hbm.20779
Lebel, C., Benischek, A., Geeraert, B., Holahan, J., Shaywitz, S., Bakhshi, K., et al. (2019). Developmental trajectories of white matter structure in children with and without reading impairments. Dev. Cogn. Neurosci. 36:100633. doi: 10.1016/j.dcn.2019.100633
Lebel, C., Shaywitz, B., Holahan, J., Shaywitz, S., Marchione, K., and Beaulieu, C. (2013). Diffusion tensor imaging correlates of reading ability in dysfluent and non-impaired readers. Brain Lang. 125, 215–222. doi: 10.1016/j.bandl.2012.10.009
Lebel, C., Walker, L., Leemans, A., Phillips, L., and Beaulieu, C. (2008). Microstructural maturation of the human brain from childhood to adulthood. Neuroimage 40, 1044–1055. doi: 10.1016/j.neuroimage.2007.12.053
Leemans, A., Jeurissen, B., Sijbers, J., and Jones, D. K. (2009). ExploreDTI: A graphical toolbox for processing, analyzing, and visualizing diffusion MR data. Proc. Int. Soc. Magn. Reson. Med. 17, 3537.
Linkersdörfer, J., Jurcoane, A., Lindberg, S., Kaiser, J., Hasselhorn, M., Fiebach, C. J., et al. (2014). The association between gray matter volume and reading proficiency: A longitudinal study of beginning readers. J. Cogn. Neurosci. 27, 308–318. doi: 10.1162/jocn_a_00710
Llera, A., Wolfers, T., Mulders, P., and Beckmann, C. F. (2019). Inter-individual differences in human brain structure and morphology link to variation in demographics and behavior. eLife 8:e44443. doi: 10.7554/eLife.44443
Long, X., Benischek, A., Dewey, D., and Lebel, C. (2017). Age-related functional brain changes in young children. Neuroimage 155, 322–330. doi: 10.1016/j.neuroimage.2017.04.059
Long, X., Kar, P., Gibbard, B., Tortorelli, C., and Lebel, C. (2019). The brain’s functional connectome in young children with prenatal alcohol exposure. NeuroImage Clin. 24:102082. doi: 10.1016/j.nicl.2019.102082
Lonigan, C. J., Anthony, J. L., Phillips, B. M., Purpura, D. J., Wilson, S. B., and McQueen, J. D. (2009). The Nature of Preschool Phonological Processing Abilities and Their Relations to Vocabulary. General Cognitive Abilities, and Print Knowledge. J. Educ. Psychol. 101, 345–358. doi: 10.1037/a0013837
Lonigan, C. J., Farver, J. A. M., Nakamoto, J., and Eppe, S. (2013). Developmental trajectories of preschool early literacy skills: A comparison of language-minority and monolingual-english children. Dev. Psychol. 49, 1943–1957. doi: 10.1037/a0031408
Lu, L. H., Leonard, C. M., Thompson, P. M., Kan, E., Jolley, J., Welcome, S. E., et al. (2007). Normal developmental changes in inferior frontal gray matter are associated with improvement in phonological processing: A longitudinal MRI analysis. Cereb. Cortex 17, 1092–1099. doi: 10.1093/cercor/bhl019
Manning, K. Y., Reynolds, J. E., Long, X., Llera, A., Dewey, D., and Lebel, C. (2021). Multimodal brain features at preschool age and the relationship with pre-reading measures one year later: An exploratory study. bioRxiv [Preprint]. doi: 10.1101/2021.04.30.442198
May, L., Byers-Heinlein, K., Gervain, J., and Werker, J. F. (2011). Language and the newborn brain: Does prenatal language experience shape the neonate neural response to speech? Front. Psychol. 2:222. doi: 10.3389/fpsyg.2011.00222
Niogi, S. N., and McCandliss, B. D. (2006). Left lateralized white matter microstructure accounts for individual differences in reading ability and disability. Neuropsychologia 44, 2178–2188. doi: 10.1016/j.neuropsychologia.2006.01.011
O’Muircheartaigh, J., Dean, D. C., Dirks, H., Waskiewicz, N., Lehman, K., Jerskey, B. A., et al. (2013). Interactions between white matter asymmetry and language during neurodevelopment. J. Neurosci. 33, 16170–16177. doi: 10.1523/JNEUROSCI.1463-13.2013
Ostertag, C., Reynolds, J. E., Dewey, D., Landman, B., Huo, Y., and Lebel, C. (2021). Altered gray matter development in pre-reading children with a family history of reading disorder. Dev. Sci. 25:e13160. doi: 10.1111/desc.13160
Peña, M., Maki, A., Kovačić, D., Dehaene-Lambertz, G., Koizumit, H., Bouquet, F., et al. (2003). Sounds and silence: An optical topography study of language recognition at birth. Proc. Natl. Acad. Sci. U.S.A. 100, 11702–11705. doi: 10.1073/pnas.1934290100
Perani, D., Saccuman, M. C., Scifo, P., Awander, A., Spada, D., Baldoli, C., et al. (2011). Neural language networks at birth. Proc. Natl. Acad. Sci.U.S.A. 108, 16056–16061. doi: 10.1073/pnas.1102991108
Qiu, D., Tan, L. H., Siok, W. T., Zhou, K., and Khong, P. L. (2011). Lateralization of the arcuate fasciculus and its differential correlation with reading ability between young learners and experienced readers: A diffusion tensor tractography study in a chinese cohort. Hum. Brain Mapp. 32, 2054–2063. doi: 10.1002/hbm.21168
Raschle, N. M., Zuk, J., and Gaab, N. (2012). Functional characteristics of developmental dyslexia in left-hemispheric posterior brain regions predate reading onset. Proc. Natl. Acad. Sci.U.S.A. 109, 2156–2161. doi: 10.1073/pnas.1107721109
Remer, J., Croteau-Chonka, E., Dean, D. C., D’Arpino, S., Dirks, H., Whiley, D., et al. (2017). Quantifying cortical development in typically developing toddlers and young children, 1–6 years of age. Neuroimage 153, 246–261. doi: 10.1016/j.neuroimage.2017.04.010
Reynolds, J. E., Long, X., Grohs, M. N., Dewey, D., and Lebel, C. (2019b). Structural and functional asymmetry of the language network emerge in early childhood. Dev. Cogn. Neurosci. 39:100682. doi: 10.1016/j.dcn.2019.100682
Reynolds, J. E., Grohs, M. N., Dewey, D., and Lebel, C. (2019a). Global and regional white matter development in early childhood. Neuroimage 196, 49–58. doi: 10.1016/j.neuroimage.2019.04.004
Reynolds, J. E., Long, X., Paniukov, D., Bagshawe, M., and Lebel, C. (2020). Calgary Preschool magnetic resonance imaging (MRI) dataset. Data Br. 29:105224. doi: 10.1016/j.dib.2020.105224
Satterthwaite, T. D., Elliott, M. A., Gerraty, R. T., Ruparel, K., Loughead, J., Calkins, M. E., et al. (2013). An improved framework for confound regression and filtering for control of motion artifact in the preprocessing of resting-state functional connectivity data. Neuroimage 64, 240–256. doi: 10.1016/j.neuroimage.2012.08.052
Saygin, Z. M., Norton, E. S., Osher, D. E., Beach, S. D., Cyr, A. B., Ozernov-Palchik, O., et al. (2013). Tracking the roots of reading ability: White matter volume and integrity correlate with phonological awareness in prereading and early-reading kindergarten children. J. Neurosci. 33, 13251–13258. doi: 10.1523/JNEUROSCI.4383-12.2013
Saygin, Z. M., Osher, D. E., Norton, E. S., Youssoufian, D. A., Beach, S. D., Feather, J., et al. (2016). Connectivity precedes function in the development of the visual word form area. Nat. Neurosci. 19, 1250–1255. doi: 10.1038/nn.4354
Schlaggar, B. L., and McCandliss, B. D. (2007). Development of neural systems for reading. Annu. Rev. Neurosci. 30, 475–503. doi: 10.1146/annurev.neuro.28.061604.135645
Segal, E., and Petrides, M. (2013). Functional activation during reading in relation to the sulci of the angular gyrus region. Eur. J. Neurosci. 38, 2793–2801. doi: 10.1111/ejn.12277
Smallwood, J., Gorgolewski, K. J., Golchert, J., Ruby, F. J. M., Engen, H., Baird, B., et al. (2013). The default modes of reading: Modulation of posterior cingulate and medial prefrontal cortex connectivity associated with comprehension and task focus while reading. Front. Hum. Neurosci. 7:734. doi: 10.3389/fnhum.2013.00734
Smith, S. M., Fox, P. T., Miller, K. L., Glahn, D. C., Fox, P. M., Mackay, C. E., et al. (2009). Correspondence of the brain’s functional architecture during activation and rest. Proc. Natl. Acad. Sci. U.S.A. 106, 13040–13045. doi: 10.1073/pnas.0905267106
Snowling, M. J., and Hulme, C. (2011). Evidence-based interventions for reading and language difficulties: Creating a virtuous circle. Br. J. Educ. Psychol. 81, 1–23. doi: 10.1111/j.2044-8279.2010.02014.x
Song, S. K., Sun, S. W., Ju, W. K., Lin, S. J., Cross, A. H., and Neufeld, A. H. (2003). Diffusion tensor imaging detects and differentiates axon and myelin degeneration in mouse optic nerve after retinal ischemia. Neuroimage 20, 1714–1722. doi: 10.1016/j.neuroimage.2003.07.005
Song, S. K., Sun, S. W., Ramsbottom, M. J., Chang, C., Russell, J., and Cross, A. H. (2002). Dysmyelination revealed through MRI as increased radial (but unchanged axial) diffusion of water. Neuroimage 17, 1429–1436. doi: 10.1006/nimg.2002.1267
Szaflarski, J. P., Holland, S. K., Schmithorst, V. J., and Byars, A. W. (2006). fMRI study of language lateralization in children and adults. Hum. Brain Mapp. 27, 202–212. doi: 10.1002/hbm.20177
Thieba, C., Frayne, A., Walton, M., Mah, A., Benischek, A., Dewey, D., et al. (2018). Factors associated with successful MRI scanning in unsedated young children. Front. Pediatr. 6:146. doi: 10.3389/fped.2018.00146
Torppa, M., Poikkeus, A. M., Laakso, M. L., Eklund, K., and Lyytinen, H. (2006). Predicting delayed letter knowledge development and its relation to grade 1 reading achievement among children with and without familial risk for dyslexia. Dev. Psychol. 42, 1128–1142. doi: 10.1037/0012-1649.42.6.1128
Treit, S., Lebel, C., Baugh, L., Rasmussen, C., Andrew, G., and Beaulieu, C. (2013). Longitudinal MRI reveals altered trajectory of brain development during childhood and adolescence in fetal alcohol spectrum disorders. J. Neurosci. 33, 10098–10109. doi: 10.1523/JNEUROSCI.5004-12.2013
Tustison, N. J., Avants, B. B., Cook, P. A., Zheng, Y., Egan, A., Yushkevich, P. A., et al. (2010). N4ITK: Improved N3 Bias Correction. IEEE Trans. Med. Imaging 29:1310. doi: 10.1109/TMI.2010.2046908
Van Ettinger-Veenstra, H., McAllister, A., Lundberg, P., Karlsson, T., and Engström, M. (2016). Higher language ability is related to angular gyrus activation increase during semantic processing, independent of sentence incongruency. Front. Hum. Neurosci. 10:110. doi: 10.3389/fnhum.2016.00110
Vigneau, M., Beaucousin, V., Hervé, P. Y., Duffau, H., Crivello, F., Houdé, O., et al. (2006). Meta-analyzing left hemisphere language areas: Phonology, semantics, and sentence processing. Neuroimage 30, 1414–1432. doi: 10.1016/j.neuroimage.2005.11.002
Vogel, A. C., Church, J. A., Power, J. D., Miezin, F. M., Petersen, S. E., and Schlaggar, B. L. (2013). Functional network architecture of reading-related regions across development. Brain Lang. 125, 231–243. doi: 10.1016/j.bandl.2012.12.016
Walton, M., Dewey, D., and Lebel, C. (2018). Brain white matter structure and language ability in preschool-aged children. Brain Lang. 176, 19–25. doi: 10.1016/j.bandl.2017.10.008
Wandell, B. A., and Yeatman, J. D. (2013). Biological development of reading circuits. Curr. Opin. Neurobiol. 23, 261–268. doi: 10.1016/j.conb.2012.12.005
Wandell, B. A., Rauschecker, A. M., and Yeatman, J. D. (2012). Learning to see words. Annu. Rev. Psychol. 63, 31–53. doi: 10.1146/annurev-psych-120710-100434
Wang, Y., and Li, T. Q. (2015). Dimensionality of ICA in resting-state fMRI investigated by feature optimized classification of independent components with SVM. Front. Hum. Neurosci. 9:259. doi: 10.3389/fnhum.2015.00259
Wang, Y., Mauer, M. V., Raney, T., Peysakhovich, B., Becker, B. L. C., Sliva, D. D., et al. (2017). Development of tract-specific white matter pathways during early reading development in at-risk children and typical controls. Cereb. Cortex 27, 2469–2485. doi: 10.1093/cercor/bhw095
Weiss-Croft, L. J., and Baldeweg, T. (2015). Maturation of language networks in children: A systematic review of 22years of functional MRI. Neuroimage 123, 269–281. doi: 10.1016/j.neuroimage.2015.07.046
Welcome, S. E., and Joanisse, M. F. (2014). Individual differences in white matter anatomy predict dissociable components of reading skill in adults. Neuroimage 96, 261–275. doi: 10.1016/j.neuroimage.2014.03.069
Welcome, S. E., Chiarello, C., Thompson, P. M., and Sowell, E. R. (2011). Reading skill is related to individual differences in brain structure in college students. Hum. Brain Mapp. 32, 1194–1205. doi: 10.1002/hbm.21101
Wolfers, T., Arenas, A. L., Onnink, A. M. H., Dammers, J., Hoogman, M., Zwiers, M. P., et al. (2017). Refinement by integration: Aggregated effects of multimodal imaging markers on adult ADHD. J. Psychiatry Neurosci. 42, 386–394.
Xiao, Y., Friederici, A. D., Margulies, D. S., and Brauer, J. (2016). Longitudinal changes in resting-state fMRI from age 5 to age 6 years covary with language development. Neuroimage 128, 116–124. doi: 10.1016/j.neuroimage.2015.12.008
Yamada, Y., Stevens, C., Dow, M., Harn, B. A., Chard, D. J., and Neville, H. J. (2011). Emergence of the neural network for reading in five-year-old beginning readers of different levels of pre-literacy abilities: An fMRI study. Neuroimage 57, 704–713. doi: 10.1016/j.neuroimage.2010.10.057
Yeatman, J. D., Dougherty, R. F., Ben-Shachar, M., and Wandell, B. A. (2012a). Development of white matter and reading skills. Proc. Natl. Acad. Sci. U.S.A. 109, 17756. doi: 10.1073/pnas.1206792109
Yeatman, J. D., Dougherty, R. F., Myall, N. J., Wandell, B. A., and Feldman, H. M. (2012b). Tract Profiles of White Matter Properties: Automating Fiber-Tract Quantification. PLoS One 7:e49790. doi: 10.1371/journal.pone.0049790
Zijlstra, H., van Bergen, E., Regtvoort, A., de Jong, P. F., and van der Leij, A. (2020). Prevention of Reading Difficulties in Children With and Without Familial Risk: Short- and Long-Term Effects of an Early Intervention. J. Educ. Psychol. 113, 248–267. doi: 10.1037/edu0000489
Keywords: neurodevelopment, MRI, reading, pediatric, diffusion imaging, resting state fMRI, brain, APrON
Citation: Manning KY, Reynolds JE, Long X, Llera A, Dewey D and Lebel C (2022) Multimodal brain features at 3 years of age and their relationship with pre-reading measures 1 year later. Front. Hum. Neurosci. 16:965602. doi: 10.3389/fnhum.2022.965602
Received: 14 June 2022; Accepted: 03 August 2022;
Published: 22 August 2022.
Edited by:
Matthew Walenski, East Carolina University, United StatesReviewed by:
John D. E. Gabrieli, Massachusetts Institute of Technology, United StatesDouglas Dean, University of Wisconsin-Madison, United States
Copyright © 2022 Manning, Reynolds, Long, Llera, Dewey and Lebel. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kathryn Y. Manning, kathryn.manning@ucalgary.ca
†These authors share first authorship