# Nurturing the reading brain: home literacy practices are associated with children's neural response to printed words through vocabulary skills

Cléa Girard, Thomas Bastelica, Jessica Léone, Justine Epinat-Duclos, Léa Longo, Jérôme Prado

► To cite this version:

Cléa Girard, Thomas Bastelica, Jessica Léone, Justine Epinat-Duclos, Léa Longo, et al.. Nurturing the reading brain: home literacy practices are associated with children's neural response to printed words through vocabulary skills. npj Science of Learning, 2021, 6 (1), pp.34. 10.1038/s41539-021-00112-9. hal-03914953

HAL Id: hal-03914953

https://hal.science/hal-03914953v1

Submitted on 11 Dec 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ARTICLE **OPEN**

<figure><img src="image_2.png" /></figure>

Check for updates

# Nurturing the reading brain: home literacy practices are associated with children's neural response to printed words through vocabulary skills

Cléa Girard<fnref n="1" />✉, Thomas Bastelica<fnref n="1" />, Jessica Léone<fnref n="1" />, Justine Epinat-Duclos<fnref n="1" />, Léa Longo<fnref n="1" /> and Jérôme Prado<fnref n="1" />✉

Previous studies indicate that children are exposed to different literacy experiences at home. Although these disparities have been shown to affect children's literacy skills, it remains unclear whether and how home literacy practices influence brain activity underlying word-level reading. In the present study, we asked parents of French children from various socioeconomic backgrounds (n = 66; 8.46 ± 0.36 years, range 7.52–9.22; 20 girls) to report the frequency of home literacy practices. Neural adaptation to the repetition of printed words was then measured using functional magnetic resonance imaging (fMRI) in a subset of these children (n = 44; 8.49 ± 0.33 years, range 8.02–9.14; 13 girls), thereby assessing how sensitive was the brain to the repeated presentation of these words. We found that more frequent home literacy practices were associated with enhanced word adaptation in the left posterior inferior frontal sulcus (r = 0.32). We also found that the frequency of home literacy practices was associated with children's vocabulary skill (r = 0.25), which itself influenced the relation between home literacy practices and neural adaptation to words. Finally, none of these effects were observed in a digit adaptation task, highlighting their specificity to word recognition. These findings are consistent with a model positing that home literacy experiences may improve children's vocabulary skill, which in turn may influence the neural mechanisms supporting word-level reading.

npj Science of Learning (2021)6:34; https://doi.org/10.1038/s41539-021-00112-9

## INTRODUCTION

Literacy skills are fundamental to children's personal and professional growth in our society. Yet, children vary greatly in their reading development throughout the school years<fnref n="1" />. Because such individual differences can be observed as early as in the preschool years, it has been argued that the home literacy environment (HLE) may play an important role in scaffolding early literacy skills². The HLE typically refers to the literacy-related interactions, resources, and attitudes that children experience at home³. Though the quality of the HLE may be estimated through various measures (e.g., title recognition, reading exposure, language input), it is most often assessed using parental self-administered questionnaires⁴,⁵. Consistent with a scaffolding role of the HLE for the development of literacy, studies have found that a supportive HLE is associated with both enhanced vocabulary⁴,⁶,⁷ and improved reading skills⁴,⁵,⁸,⁹ in children. This has notably led some to hypothesize that the relation between the HLE and children's reading may be mediated by vocabulary⁵,¹⁰,¹¹. In other words, more frequent home literacy activities may improve children's reading efficiency by fostering vocabulary skills (which are an important foundation for word recognition)¹² (see Fig. 1a). This also prompted various organizations to call for greater parent and caregiver involvement in home literacy education¹³,¹⁴.

What are the neural mechanisms supporting the relation between the HLE and literacy development in children? The first line of evidence comes from studies investigating how differences in children's brain structure and function may relate to disparities in family socioeconomic status (SES), an index of family status and position in society (often obtained using measures of parental income, education, or occupation). Indeed, differences in SES are well known to be associated with differences in language input

and literacy environment¹⁵,¹⁶. Overall, these studies have found that SES is associated with structural and functional differences in brain regions typically associated with literacy. For example, studies have found a positive relation between emergent literacy skills and cortical thickness in left perisylvian regions¹⁷. However, lower SES has been associated with lesser cortical thickness in these regions¹⁸,¹⁹. Children from disadvantaged families have also less cortical surface area²⁰ in left perisylvian regions than more advantaged peers. Finally, these children have been found to have a weaker specialization of the left inferior frontal gyrus (IFG) during a rhyming task²¹. Thus, there is growing evidence that SES is associated with the structure and function of left perisylvian areas in children, a relation that might be mediated by cognitive precursors to reading acquisition, such as vocabulary²² (see Fig. 1a).

The second line of evidence comes from recent studies that have more directly investigated the relation between HLE and brain development, either using parental questionnaires or home recordings. For instance, home reading exposure and maternal shared reading quality have been associated with white matter integrity²³ and left perisylvian activity during story-listening tasks in preschoolers²⁴,²⁵. Using home audio recordings, studies have also showed that the quality of conversational experience with parents at home is associated with greater surface area²⁶ and structural connectivity²⁷ within the left perisylvian system in children, as well as increased activity in the left IFG during a story-listening task²⁸. Another study employing parental questionnaires has found that structural connectivity of left perisylvian areas is associated with both HLE and vocabulary development²⁹, in line with the idea that both dimensions may interact to account for reading acquisition²². Finally, more supportive HLE has also been

<footnote n="1">¹Lyon Neuroscience Research Center (CRNL), INSERM U1028 - CNRS UMR5292, University of Lyon, 69500 Bron, France. ✉email: clea.girard@etu.univ-lyon1.fr; jerome.prado@univ-lyon1.fr</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Model and experimental design. **a** Hypothesized model of the relations between home literacy practices, vocabulary and reading. **b** Experimental task. Participants were sequentially presented with words (top) or digits (bottom) that were identical (adaptation blocks) or different (no-adaptation blocks).</figcaption></figure>

linked to increased activity in the left IFG during a sound matching task in 5-year-olds³⁰. Therefore, in keeping with studies investigating SES, neuroimaging studies that have more explicitly assessed the HLE also suggest that the quality of adult-child literacy interactions at home might affect both structure and function of left perisylvian areas (particularly at the level of the left IFG).

These previous studies, however, only provide indirect evidence that the HLE relates to brain regions associated with word-level reading in children. This is because studies to date have either assessed structural (not functional) properties of the brain of reading children²⁶,²⁹ or have only used functional magnetic resonance imaging (fMRI) to measure the brain activity of pre-reading children using auditory tasks that did not involve word recognition²⁴,²⁵,²⁸,³⁰. Thus, although one previous study has assessed the relation between SES and brain regions supporting pseudo-word processing¹⁶, we do not know if the frequency of home literacy practices is associated with fMRI activity during a word recognition task in children. We also do not know whether this relation may be mediated by pre-literacy skills such as vocabulary, as some have hypothesized²².

In the present study, we aimed to test whether fMRI activity in areas supporting word recognition may be related to the HLE of elementary school children. To test this, parents of 8-year-olds completed an extensive questionnaire evaluating the frequency of literacy practices shared with children at home. Children's literacy skills were also assessed using measures of reading fluency and vocabulary. Finally, fMRI activity associated with word recognition was measured using a word adaptation task developed by Perrachione et al.³¹. In this task, participants are passively presented with blocks of visual (written) words that are either identical (adaptation blocks) or different (no-adaptation blocks) (see Fig. 1b). Such fMRI-adaptation paradigms rely on the idea that presenting a given stimulus repeatedly leads to a decrease of fMRI signal in the region that processes this stimulus because firing rates of the neuronal population decrease³². Thus, how sensitive a region is to word processing can be measured by comparing no-adaptation to adaptation blocks, i.e., the so-called neural adaptation effect³¹. Critically, this neural adaptation effect appears to be functionally relevant because it is related to reading skills in adults (i.e., it is smaller in dyslexic readers than in typical readers³¹). To compare our findings to that of Perrachione et al.³¹ (who tested English-speaking adults) and therefore assess task replicability³³, we first identified the network in which a neural adaptation effect was found across the whole brain in our

sample of French-speaking children. To increase statistical power³⁴ and provide an unbiased estimate of effect sizes, we then employed a hypothesis-driven approach and measured the word adaptation effect in the exact same regions that were found to exhibit a word adaptation effect in Perrachione et al. in typical readers³¹.

We made five predictions. First, we expected that the frequency of home literacy practices would be associated with both vocabulary and reading fluency skills in children, in line with previous behavioral studies⁵. Second, we anticipated replicating and extending Perrachione et al.'s findings by showing word adaptation effects in the regions that were identified in that previous study using the same task in French-speaking children. These were our regions of interest (ROIs). Third, we expected that more frequent home literacy practices would also be associated with greater neural adaptation to words (i.e., a larger difference in fMRI activity between no-adaptation and adaptation blocks) in at least some of the ROIs, reflecting a relation between HLE and fMRI activity associated with word recognition. Fourth, in line with the model shown in Fig. 1a, we hypothesized that the relation between HLE and both reading fluency skills and word adaptation would be mediated by children's vocabulary skills. Fifth, we expected that none of the above effects would be observed in a control adaptation task in which words were replaced by other types of stimuli. Therefore, we also presented participants with a control adaptation task in which words were replaced by Arabic digits (see Fig. 1b). This allowed us to test whether any of our effects were specific to words (as compared to any kind of symbolic information), thereby also controlling for general perceptual adaptation. Indeed, processing symbolic numbers involves parietal mechanisms that support quantity processing to a much greater extent than the brain circuits of reading³⁵.

## RESULTS

### Samples

Demographic information about children and parents can be found in Table 1 and Methods. The behavioral sample consisted of 66 French-speaking children from age 7.52 to 9.22 (mean = 8.46), while the fMRI sample consisted of 44 of these participants (age range = 8.02–9.14; mean = 8.49). SES ranged from relatively low to relatively high in both samples (see Table 1 and Methods for details). Age-normalized scores for reading fluency and vocabulary skills (assessed using the Alouette-R test³⁶ and the "vocabulaire" subtest of the NEMI-2 test³⁷, see Methods) were in the normal to superior range in both samples. IQ was also in the normal to the superior range. Parental reading scores were in the normal range (see Table 1). The subset of children in the fMRI sample did not differ from children in the behavioral sample in terms of demographics and scores (all *p*'s > 0.27).

### Behavioral results

Figure 2 shows the frequency of the different home literacy practices reported by parents across items in our questionnaire (parents indicated the frequency of each practice on a six-point rating scale, see Methods). Practices were categorized according to whether they were informal or formal as well as basic or advanced, in line with the previous literature on the HLE⁵,³⁸. As can be seen on Fig. 2, there was large individual variability in the reported frequencies of home literacy practices. Frequencies were then averaged across items within each category to gather separate scores of informal, formal basic, and formal advanced practices. Between parents, average responses on the rating scale ranged from 1 to 3.75 for practices that were considered informal (M = 2.080, SD = 0.548), from 0.5 to 4 for practices that were considered formal but basic for the child's age (M = 2.140, SD = 0.876), and from 0.670 to 3.840 for practices that were considered
<table><caption>**Table 1.** Demographic information and test scores of children and parents.</caption><thead><tr><th rowspan="3">Measure</th><th colspan="4">Behavioral sample (n = 66)</th><th colspan="4">fMRI sample (n = 44)</th></tr><tr><th colspan="2">Children</th><th colspan="2">Parents</th><th colspan="2">Children</th><th colspan="2">Parents</th></tr><tr><th>Mean (SD) or n</th><th>Range</th><th>Mean (SD) or n</th><th>Range</th><th>Mean (SD) or n</th><th>Range</th><th>Mean (SD) or n</th><th>Range</th></tr></thead><tbody><tr><td>Age (years)<sup>a</sup></td><td>8.46 (0.36)</td><td>7.52-9.22</td><td>39.53 (4.95)</td><td>29-51</td><td>8.49 (0.33)</td><td>8.02-9.14</td><td>39.25 (4.78)</td><td>32-51</td></tr><tr><td>Sex (female)</td><td>n = 20</td><td>-</td><td>n = 59</td><td>-</td><td>n = 13</td><td>-</td><td>n = 39</td><td>-</td></tr><tr><td>Child grade (second grade)<sup>b</sup></td><td>n = 20</td><td>-</td><td>-</td><td>-</td><td>n = 12</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Parental education<sup>c</sup></td><td>-</td><td>-</td><td>3.05 (2.36)</td><td>-7-8</td><td>-</td><td>-</td><td>3.41 (2.17)</td><td>-1-8</td></tr><tr><td>primary</td><td>-</td><td>-</td><td>n = 1</td><td>-</td><td>-</td><td>-</td><td>n = 0</td><td>-</td></tr><tr><td>secondary</td><td>-</td><td>-</td><td>n = 8</td><td>-</td><td>-</td><td>-</td><td>n = 6</td><td>-</td></tr><tr><td>undergraduate</td><td>-</td><td>-</td><td>n = 33</td><td>-</td><td>-</td><td>-</td><td>n = 20</td><td>-</td></tr><tr><td>master or higher</td><td>-</td><td>-</td><td>n = 24</td><td>-</td><td>-</td><td>-</td><td>n = 17</td><td>-</td></tr><tr><td>Monthly income<sup>d</sup></td><td>-</td><td>-</td><td>1727 (1078)</td><td>500-5500</td><td>-</td><td>-</td><td>1659 (1098)</td><td>500-5500</td></tr><tr><td>€0-1000</td><td>-</td><td>-</td><td>n = 18</td><td>-</td><td>-</td><td>-</td><td>n = 14</td><td>-</td></tr><tr><td>€1000-2000</td><td>-</td><td>-</td><td>n = 25</td><td>-</td><td>-</td><td>-</td><td>n = 15</td><td>-</td></tr><tr><td>€2000-3000</td><td>-</td><td>-</td><td>n = 16</td><td>-</td><td>-</td><td>-</td><td>n = 11</td><td>-</td></tr><tr><td>€3000-4000</td><td>-</td><td>-</td><td>n = 5</td><td>-</td><td>-</td><td>-</td><td>n = 3</td><td>-</td></tr><tr><td>€4000-5000</td><td>-</td><td>-</td><td>n = 1</td><td>-</td><td>-</td><td>-</td><td>n = 0</td><td>-</td></tr><tr><td>€5000-6000</td><td>-</td><td>-</td><td>n = 1</td><td>-</td><td>-</td><td>-</td><td>n = 1</td><td>-</td></tr><tr><td>Reading accuracy<sup>e</sup></td><td>105 (8)</td><td>86-120</td><td>96 (15)</td><td>40-115</td><td>106 (8)</td><td>86-120</td><td>95 (16)</td><td>40-115</td></tr><tr><td>Reading speed<sup>e</sup></td><td>112 (16)</td><td>81-155</td><td>107 (13)</td><td>74-157</td><td>114 (17)</td><td>81-155</td><td>106 (13)</td><td>74-147</td></tr><tr><td>Vocabulary<sup>f</sup></td><td>9 (2)</td><td>4-14</td><td>-</td><td>-</td><td>10 (2)</td><td>6-14</td><td>-</td><td>-</td></tr><tr><td>IQ composite<sup>e</sup></td><td>112 (11)</td><td>83-135</td><td>-</td><td>-</td><td>113 (12)</td><td>83-135</td><td>-</td><td>-</td></tr></tbody></table>

There was no difference between the behavioral and the fMRI samples in terms of demographics or performance (all *p*'s > 0.27).

<sup>a</sup>All children were 8 years of age, except for one who was younger (7.51) and three who were slightly older (9.06 to 9.22) at the time of testing.

<sup>b</sup>Children were either in second or third grade.

<sup>c</sup>Education level was measured using the number of years pre- or post- high school graduation, with high school graduation set at 0.

<sup>d</sup>As a reference, the median monthly income in France is about €1,700 (Robin, 2019).

<sup>e</sup>Standard score.

<sup>f</sup>Scaled score.

formal and relatively advanced (M = 2.006, SD = 0.695). Correlation analyses indicated that neither education nor income was associated with home literacy practices in the behavioral sample (all *p*'s > 0.19, two-tailed). However, education tended to be associated with both basic and formal advanced home literacy practices in the fMRI sample (basic: *r*(64) = -0.259, *p* = 0.090, two-tailed; advanced: *r*(64) = -0.280, *p* = 0.066, two-tailed). As is standard in the literature on the home learning environment, all further analyses were controlled for parental SES (parental income and education).

We first investigated the relations between overall home literacy practices (a composite measure calculated by summing the frequency of informal, formal basic, and formal advanced practices) and both reading fluency and vocabulary across this entire behavioral sample. Our first hypothesis was only partly supported. In line with expectations, partial correlations controlling for parental income and education indicated a small positive relation between overall frequency of home literacy practices and child vocabulary (partial *r*(64) = 0.250, *p* = 0.023). However, there was no positive relation between overall frequency of home literacy practices and either child reading accuracy (partial *r*(64) = 0.148, *p* = 0.126) or child reading speed (partial *r*(64) = 0.121, *p* = 0.174) (see Fig. 3). Follow-up analyses revealed that the relation between home literacy practices and vocabulary was relatively consistent across all types and complexity of literacy practices (informal: partial *r*(64) = 0.187, *p* = 0.070, basic formal: partial *r*(64) = 0.182, *p* = 0.075; advanced formal: partial *r*(64) = 0.240, *p* = 0.028). Therefore, neuroimaging analyses focus on the composite measure of the overall frequency of literacy practices.

## fMRI results

We then turned to neuroimaging data and analyzed adaptation effects measured in the word adaptation task (see Fig. 1b) in the fMRI sample. Critically, this task was used in a previous study by Perrachione et al.<sup>31</sup> (see their Experiment 2B), who investigated neural adaptation to words in a sample of typical (n = 23, standardized reading score = 108.0 ± 6.7, see Table 1 in Perrachione et al.<sup>31</sup>) and dyslexic (n = 23, standardized reading score = 84.2 ± 6.6, see Table 1 in Perrachione et al.<sup>31</sup>) adult readers. In typical readers, Perrachione et al.<sup>31</sup> found widespread word adaptation effects in a left-lateralized network, including the temporal (fusiform gyrus, inferior temporal gyrus, superior, and middle temporal gyrus), frontal (inferior frontal gyrus, frontal operculum, premotor cortex, and pre-supplementary motor area), and visual (pericalcarine) cortices (see Fig. 3b and Table S5 in Perrachione et al.<sup>31</sup>). In dyslexic readers, word adaptation effects were weaker and centered around the left inferior frontal cortex (with more limited adaptation in the left middle temporal and inferior temporal cortices).

In our sample of French-speaking 8-year-olds, whole-brain analyses (see Methods) showed enhanced word adaptation effects across all participants in a left-lateralized network that was largely similar to that found in typical readers in Perrachione et al.<sup>31</sup>. In our sample, word adaptation effects were also found in the temporal (superior, middle, and inferior temporal gyrus), frontal (notably the left triangular and orbital parts of the inferior frontal gyrus), and visual (inferior occipital gyrus) cortices (see Fig. 4a and Table 2). There was, however, less extensive activation in the left
<figure><img src="image_4.png" /><figcaption>**Fig. 2** Divergent stacked bar chart showing the frequencies of informal, basic formal, and advanced formal practices among parents (n = 66). Percentages on the right indicate the share of parents who engage in each practice at least once a week. Practices are ordered from the most (top) to the least frequent (bottom) for each level of complexity.</figcaption></figure>

precentral gyrus and left parietal cortex (see Fig. 3b in Perrachione et al.³¹).

Given the similarity between the brain regions showing a word adaptation effect for typical readers in Perrachione et al.³¹ and in our study (see above), we used the activation peaks for typical readers reported in Table S5 of Perrachione et al.³¹ to define regions of interest (ROIs) in the present study (see Methods). Specifically, Perrachione et al. (2016) reported peaks in the left posterior inferior frontal sulcus (pIFs; MNI coordinates: x = -39, y = 11, z = 26), left posterior dorsal superior temporal sulcus (pdSTS; MNI coordinates: x = -56, y = -39, z = 6), left pre-supplementary motor area (pre-SMA; MNI coordinates: x = -6, y = 10, z = 52), left occipital cortex (OC; MNI coordinates: x = -17, y = -76, z = 12), left posterior temporal fusiform gyrus (pTF; MNI coordinates: x = -38, y = -41, z = -19), left Putamen (MNI coordinates: x = -21, y = 7, z = 0), and right Putamen (MNI coordinates: x = 22, y = 12, z = -1) (see Fig. 4b). Note that the pTF corresponds to the anterior part of the visual word form area (VWFA)³⁹,⁴⁰.

Averaging across all children from the fMRI sample (n = 44), word adaptation effects were significant after Bonferroni correction for multiple comparisons in three of these ROIs: the pIFs (t(43) = 2.698, *p*<sub>corr</sub> = 0.035, *d* = 0.407), pdSTS (t(43) = 3.794, *p*<sub>corr</sub> = 0.002, *d* = 0.572), and pTF (t(43) = 4.613, *p*<sub>corr</sub> < 0.001, *d* = 0.695). The word adaptation effect, however, was only partially significant the putamen (left Putamen: t(43) = 2.674, *p*<sub>corr</sub> = 0.037, *d* = 0.403, right Putamen: t(43) = 1.912, *p*<sub>corr</sub> = 0.217, *d* = 0.288) and was not

significant in the pre-SMA (t(43) = 1.335, *p*<sub>corr</sub> = 0.658, *d* = 0.201) and the OC (word: t(43) = -3.103, *p*<sub>corr</sub> > 0.999, *d* = -0.468) (see Fig. 5a).

We then tested our third prediction that the overall frequency of home literacy practices would be related to word adaptation effects in ROIs. We focused here on the ROIs for which we found significant word adaptation effects across all participants (excluding the putamen as a word adaptation effect was not uniformly found in that ROI). After correcting for multiple comparisons, there was a positive relation between frequency of home literacy practices and word adaptation in the pIFs (r(42) = 0.322, *p*<sub>corr</sub> = 0.048), but not in the pdSTS (r(42) = 0.052, *p*<sub>corr</sub> > 0.999) and pTF (r(42) = -0.056, *p*<sub>corr</sub> > 0.999). The relation between frequency of home literacy practices and word adaptation remained significant in the pIFs after controlling for parental income and education (partial r(42) = 0.302, *p* = 0.026) (see Fig. 5b, left).

Finally, our fourth prediction was that the relation between HLE and either reading fluency or word adaptation would be influenced by children's vocabulary skills. Because we did not find a relation between home literacy practices and reading fluency, we focused on word adaptation. First, we tested whether vocabulary was related to word adaptation effects in the pIFs, the only region in which there was a relation between activity and frequency of home literacy practices (see above). Controlling for parental income and education, there was a positive relation between vocabulary and word adaptation in the pIFs (partial r(42) = 0.317, *p* = 0.020) (see Fig. 5b, right). Therefore, the frequency of
<figure><img src="image_5.png" /><figcaption>**Fig. 3 Behavioral results.** Scatterplots of child reading accuracy, reading speed, and vocabulary as a function of the overall frequency of home literacy practices (controlled for parental income and education, *n* = 66).</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 4 Whole-brain neural adaptation to words and location of ROIs.**
**a** Brain regions more activated for no-adaptation than adaptation blocks in the word adaptation task across the whole brain. **b** Location of ROIs on lateral and ventral views of an inflated rendering of the cortical surface and an axial slice of an MNI normalized brain. pIFs left inferior frontal sulcus, pdSTS left posterior dorsal superior temporal sulcus, pre-SMA left pre-supplementary motor area, OC left occipital cortex, pTF left posterior temporal fusiform gyrus.</figcaption></figure>

home literacy practices was related to both vocabulary and word adaptation in the pIFs. Vocabulary was then entered as a potential mediator of the relation between overall frequency of home literacy practices and word adaptation in the pIFs in a mediation analysis (see Methods and Fig. 6). As expected, the total relationship between home literacy practices and word adaptation was significantly positive, even in the fMRI sample (path c, coefficient = 0.063, standard error, or STE = 0.033, 95% CI Lower = 0.011, *p* = 0.024). Also as expected, more frequent home literacy practices were associated with improved vocabulary (path a, coefficient = 0.427, STE = 0.198, 95% CI Lower = 0.063, *p* = 0.024). Improved vocabulary also tended to be associated with larger word adaptation even after controlling for frequency of home literacy practices (path b, coefficient = 0.036, STE = 0.026, 95% CI Lower = -0.003, *p* = 0.062). Critically, the relation between frequency of home literacy practices and word adaptation was significantly mediated by vocabulary (path a*b, coefficient = 0.015, STE = 0.014, 95% CI Lower = 0.001, *p* = 0.046). Finally, the relation between frequency of home literacy practices and word adaptation was no longer significant after accounting for vocabulary (direct effect or path c', coefficient = 0.048, STE = 0.031, 95% CI Lower = -0.001, *p* = 0.052).

## Additional analyses

We performed four sets of additional control analyses. First, to evaluate the specificity of our effects for print words, we measured

the relation between adaptation to digits in the ROIs and both home literacy practices and vocabulary. There was no positive relation between digit adaptation and either frequency of home literacy practices or vocabulary in any of the seven ROIs (home literacy practices: all *r*'s < 0.004, all *p*'s > 0.493; vocabulary: all *r*'s < 0.188, all *p*'s > 0.117), including the pIFs (see Supplementary Fig. 1 and Supplementary Table 1). Figure 7 shows the unthresholded *t*-maps of the relation between either word or digit adaptation and home literacy practices and vocabulary across the whole brain.

Second, adaptation effects were measured using a passive task. Therefore, it is possible that individual differences in task engagement might influence our results. To test for this possibility, we embedded a target detection task among the stimuli (see Methods). On average, children detected 89% of targets (SD = 15) in the word adaptation task and 89% of targets (SD = 19) in the digit adaptation task. There was no difference in target detection rate between tasks (Wilcoxon rank-sum test: *p* = 0.730). Critically, target detection rates were not correlated with home literacy practices or vocabulary in any of the tasks (word adaptation: all rhos < -0.045, all *p*'s > 0.615; digit adaptation: all rhos < 0.114, all *p*'s > 0.231).

Third, although it has been hypothesized that home literacy practices are related to children's reading through improved vocabulary, it is also possible that parents engage in more frequent literacy practices at home because they perceive that their child is skilled in reading (which may result in both improved vocabulary and enhanced brain sensitivity to print words). To tease apart these possibilities, we asked parents to indicate their perception of their child's reading skill. This parental estimate was positively correlated with both reading accuracy (*r*(64) = 0.386, *p* < 0.001) and reading speed (*r*(64) = 0.534, *p* < 0.001), though it was not positively correlated with vocabulary (*r*(64) = 0.124, *p* = 0.161). However, this parental perception of reading skills was not correlated with home literacy practices (controlling for parental income and education, partial *r*(64) = 0.025, *p* = 0.421) or word adaptation in the pIFs (controlling for parental income and education, partial *r*(42) = 0.115, *p* = 0.229). When controlling for this parental estimate of child's reading skill, literacy practices remained related to vocabulary (controlling for parental income and education, partial *r*(64) = 0.249, *p* = 0.022) and word adaptation in the pIFs (controlling for parental income and education, partial *r*(42) = 0.302, *p* = 0.023).

Finally, associations between literacy practices and children's skills may also result (at least to some extent) from genetic predispositions to high academic achievement that parents may pass on to children. For instance, skilled parents might maintain relatively high-quality literacy environments⁴¹. Following Hart, Little, and Bergen⁴², we attempted to control for this effect by measuring parental reading fluency as a genetic proxy for predisposition to literacy skill. Critically, the relation between literacy practices and vocabulary remained significant after controlling for parent reading accuracy and reading speed
<table><caption>**Table 2.** Activation peaks in the whole-brain analysis of the word adaptation task.</caption><thead><tr><th rowspan="2">Anatomical Location</th><th rowspan="2">Cluster size (mm³)</th><th colspan="3">MNI coordinates</th><th rowspan="2">t-score</th></tr><tr><th>X</th><th>Y</th><th>Z</th></tr></thead><tbody><tr><td>Left occipital cortex</td><td>7224</td><td>-22</td><td>-98</td><td>-4</td><td>6.03</td></tr><tr><td>Left inferior temporal gyrus</td><td>-</td><td>-46</td><td>-46</td><td>-15</td><td>5.35</td></tr><tr><td>Left inferior occipital cortex</td><td>-</td><td>-38</td><td>-86</td><td>-8</td><td>5.31</td></tr><tr><td>Left inferior frontal gyrus pars orbitalis</td><td>1246</td><td>-36</td><td>32</td><td>-8</td><td>4.65</td></tr><tr><td>Left inferior frontal gyrus pars orbitalis</td><td>-</td><td>-42</td><td>28</td><td>-15</td><td>4.18</td></tr><tr><td>Left inferior frontal gyrus pars orbitalis</td><td>-</td><td>-26</td><td>30</td><td>-12</td><td>3.44</td></tr><tr><td>Left inferior frontal gyrus pars triangularis</td><td>840</td><td>-52</td><td>26</td><td>16</td><td>4.61</td></tr><tr><td>Left superior temporal gyrus</td><td>938</td><td>-58</td><td>-40</td><td>6</td><td>4.50</td></tr><tr><td>Left middle temporal gyrus</td><td>-</td><td>-62</td><td>-48</td><td>6</td><td>4.43</td></tr></tbody></table>

MNI Montreal neurological institute.

a

<figure><img src="image_7.png" /></figure>

b

<figure><img src="image_8.png" /><figcaption>**Fig. 5** Neural adaptation to words (n = 44). **a** Average neural adaptation to words in the ROIs, ranked by effect size. The dark outline shows the brain regions for which word adaptation was significant after correction for multiple comparisons (excluding the left putamen, see text). Error bars represent the standard error of the mean (SEM). Each dot represents a participant. **b** Scatterplots of neural adaptation to words in the pIFs as a function of the overall frequency of home literacy practices (left) and child vocabulary (right). Analyses were controlled for parental income and education.</figcaption></figure>

(also controlling for parental income and education, partial *r*(64) = 0.248, *p* = 0.022).

DISCUSSION

Studies indicate large variations in the home literacy environment of children. Although these variations have been shown to relate to children's reading skills, the relation between home literacy practices and brain activity underlying word-level reading remains unknown. The main goal of this study was to assess whether home literacy practices are associated with activity in brain

regions supporting word recognition in 8-year-olds. We also aimed to test whether the relation between home literacy practices and activity associated with word recognition may be influenced by children's vocabulary, as suggested by previous studies<sup>22,29</sup>.

First, based on prior literature, we predicted that the frequency of home literacy practices would be positively related to both children's reading fluency and vocabulary skills. This hypothesis was only partially supported. We did find a relation between more frequent parental practices and enhanced children's vocabulary, in keeping with a long line of previous studies on the HLE<sup>4,6,43</sup>.
However, there was no significant relation between home literacy practices and any of our measures of reading fluency (accuracy or speed). Although this may be surprising given multiple reports of such a relation in the literature<sup>5,8</sup>, we could think of at least two reasons why we failed to find this association here. A first possibility is that the reading fluency test we used<sup>36</sup> might have lacked sensitivity to discriminate between the reading fluency of our participants. Indeed, the Alouette test was initially designed to screen for dyslexia among children based on phonological-recoding skills<sup>44</sup>, which are early precursors of decoding skills. However, all of our participants were typical readers and it is possible that this 3-min test might not have been sensitive enough to capture variance within this population. A second possible reason is that the influence of the home environment on reading fluency skills may vary depending on children's age<sup>1,45</sup>. Our participants were 8-year-olds who were already relatively efficient readers. For example, it is possible that individual differences in the reading fluency of these children might have been better captured by a reading comprehension task<sup>46</sup>. Note that we also did not find a relation between the frequency of home literacy practices and SES. Although this is inconsistent with several previous studies<sup>15,16</sup>, it is possible that this lack of relation might be due to the relatively high concentration of children from middle to high SES families in our sample.

Second, we assessed the replicability of our task by evaluating whether word adaptation effects would be found across participants in the same regions reported by Perrachione et al.<sup>31</sup>, who used the same task in a different population (adults) and language (English). Both whole-brain and ROI analyses

showed that we largely replicated Perrachione et al.'s findings<sup>31</sup>. That is, word adaptation effects were found in a left-lateralized brain network that encompassed the inferior frontal, temporal, and occipital cortices. Not only are these regions located in the classic reading brain network<sup>47,48</sup>, these findings suggest that the task robustly assesses neural adaptation to print words in both French and English and generalizes across two different samples. This is significant given recent concerns about the reliability and replicability of fMRI studies<sup>49</sup>.

Third, we predicted that home literacy practices would be associated with activity in brain regions supporting word recognition. Somewhat surprisingly, we did not find any relation between home literacy practices and brain regions associated with visual processing, including the VWFA that has been shown to be specialized for letter strings<sup>50</sup>. Because this region develops during the first years of reading acquisition in children<sup>51</sup>, it is possible that VWFA activity during word recognition might be influenced by home literacy practices in younger children. In any case, our results suggest that home literacy practices may not strongly influence the visual processing of words in 8-year-olds. However, we found that home literacy practices were associated to word adaptation in the left pLFs in our sample of 8-year-olds. In other words, children benefiting from a higher quality home literacy environment had heightened sensitivity to the repetition of print words in an area located in the dorsal and posterior part of the left inferior frontal gyrus (IFG). This area, together with the larger adjacent left IFG pars triangularis, has consistently been associated with reading in both adults and children, as shown in a meta-analysis of studies involving reading or reading-related tasks with visual word, nonword, or letter string stimuli<sup>50</sup>. The specific role of this region - and of the larger anterior reading circuit- during reading remains unclear and may involve both phonological and semantic processing<sup>52,53</sup>. Nonetheless, there is evidence that reading skills are positively correlated with activity associated with print words in the pars triangularis in a sample of children from third and fourth grade<sup>54</sup>. Interestingly, several studies have previously found associations between the home literacy environment and the structure and function of the left IFG more broadly<sup>21,25,28,30</sup>. However, to our knowledge, none of these studies have directly measured activity associated with word-level reading in this region. Therefore, our study provides the first evidence that disparities in children's home literacy environment may relate to functional differences in how the dorsal part of the left IFG processes words.

Fourth, based on prior literature<sup>22,29</sup> we aimed to assess whether vocabulary may mediate the relation between home literacy practices and activity associated with word recognition.

<figure><img src="image_9.png" /><figcaption>**Fig. 6** Mediation analysis (n = 44). The relation between home literacy practices and word adaptation in the pLFs was influenced by vocabulary (see main text). The solid line from home literacy practices to word adaptation in the pLFs represents the total effect (Path c). The dotted line represents the direct effect (Path c'). The analysis was controlled for parental income and education.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 7** Unthresholded maps of *t* values showing the relations between neural adaptation, home literacy practices, and child vocabulary across the whole brain (n = 44). a Lateral views of the left and right hemispheres showing regions in which neural adaptation to words (left) and digits (right) increased with the overall frequency of home literacy practices. b Lateral views of the left and right hemispheres showing regions in which word adaptation (left) and digit adaptation (right) increased with child vocabulary.</figcaption></figure>
Not only did we find an association between children's vocabulary skills and word adaptation in the left pLFS, our results indicate that the relation between home literacy practices and word adaptation in the left pLFS was mediated by children's vocabulary. This is consistent with a model positing that children's vocabulary skills mediate the relation between home literacy environment and mechanisms supporting word-level reading (see Fig. 1a). This is also consistent with previous behavioral studies showing that HLE effects on literacy outcomes may be mediated by earlier precursor skills, most notably vocabulary<sup>3,29,55</sup>. Overall, these findings support the general idea that the relation between the home environment and the brain mechanisms supporting reading in children is influenced by cognitive precursors to reading, such as vocabulary<sup>22</sup> (see Fig. 1a).

Fifth, we anticipated that the relations between home literacy practices and neural adaptation would not be specific to words. Therefore, we also presented participants with an adaptation task in which words were replaced by Arabic numerals. We did not find any significant relation between home literacy practices and digit adaptation in any of our ROIs. There was also no significant relation between vocabulary and digit adaptation in the ROIs. Therefore, our effects cannot be easily explained by low-level perceptual adaptation effects or adaptation effects involving domain-general symbolic processing.

Finally, it is important to consider some limitations of this study. For example, our cross-sectional design does not allow us to formally test causal claims in Fig. 1a. Specifically, the relations that we observed are in line with this model. However, there may be other models that are compatible with our findings (e.g., the frequency of home literacy practices may increase in response to children's reading skills). Our measure of the HLE also relies on a questionnaire. Not only may questionnaires be associated with a recall bias (i.e., parents may inaccurately or incompletely remember events that happened in the past), they might also be associated with a social desirability bias (i.e., parents may have reported more frequent literacy practices in an effort to avoid embarrassment and project a more favorable image to the experimenters)<sup>56</sup>. Finally, our study is restricted in terms of sample size, age range investigated, and reading task. Clearly, reading encompasses much more than word recognition and it is unclear to what extent our findings may apply to other age groups (e.g., home literacy practices may change with age) as well as to the relation between home literacy practices and brain mechanisms underlying reading comprehension or other higher-level reading processes.

We attempted to minimize these limitations in several ways. For instance, although our study is fundamentally correlational, the model we test here is based on extensive prior research on the HLE<sup>5,11,22,29</sup>. We also attempted to control for other plausible interpretations by gathering information about parental skills and perception of children's reading skill (see Additional analyses). Even though recall and social desirability biases may still be present, we attempted to limit these biases by asking parents to report information about a very wide range of home activities and expectations (e.g., math, sport, music). Although our sample size is similar to that of previous studies that explored associations between home literacy practices and the reading network in the brain<sup>26,28</sup>, we adopted ROI analyses to maximize power<sup>34</sup> (though a drawback of this approach is that it leaves open the possibility that other meaningful relations may be observed elsewhere in the brain). Notwithstanding these attempts to address limitations, it is clear that our results need to be substantiated by interventional studies in order to test the causal claims suggested here, as well as by future neuroimaging studies investigating the relation between home literacy practices and the reading brain in younger children and different tasks.

To conclude, we found that home literacy practices are associated with brain mechanisms supporting word recognition

in 8-year-olds, in line with the idea that experiences beyond the classroom may influence brain regions supporting reading and learning<sup>57</sup>. Our results also support a model suggesting that home literacy practices may affect children's reading through precursors skills such as vocabulary<sup>10,11</sup>. Nonetheless, our study was cross-sectional and the effects found were relatively small in size. Therefore, it calls for future longitudinal and interventional studies investigating how home literacy practices in younger children shape the development of the brain networks supporting reading.

# METHODS

## Participants

Seventy-three right-handed children from second and third grade and one of their parents were recruited through flyers sent to schools and advertisements on social media. The experiment involved two sessions. In the first session, parents and children completed tests and questionnaires in the lab. In the second session, children completed the experimental tasks in the scanner. Seven children were excluded from analyses of the first session because they (1) were seeing a speech-language pathologist on a regular basis (n = 3), (2) had an intelligence quotient (IQ) lower than the 25th percentile (n = 2), (3) had a delay in speech and language acquisition (n = 1), and (4) were diagnosed with attention deficit disorder (n = 1). Therefore, 66 children were included in the behavioral sample. Out of the 58 children who participated the second (i.e., fMRI) session, 14 were excluded from the fMRI analyses because of incomplete data acquisition (n = 7) or excessive motion in the scanner (n = 7) (see criteria below). Therefore, our fMRI sample consisted of 44 children who had at least one run of data analyzable in both the word and digit adaptation tasks.

Children had a full-scale IQ comprised between 83 and 135 (mean = 112, standard deviation [SD] = 11), as measured by the NEMI-2 test<sup>37</sup>. A previous study<sup>58</sup> explored the relation between brain activity and home numeracy practices in the same sample of participants. Demographic information about children and parents can be found in Table 2. All children and parents were native French speakers. Eighty-nine percent of the parents were mothers. Parental income ranged from less than €6000 to more than €60,000 per year (the median annual income in France is about €20,400;<sup>59</sup>). Fourteen percent of parents reported only having a secondary degree, 50% reported having an undergraduate degree, and 36% a master's degree or higher. One parent did not go to high school. Therefore, SES ranged from relatively low to relatively high. Parents gave written informed consent and children gave their assent to participate in the study. The study was approved by a French ethics committee (Comité de Protection des Personnes Sud-Est 2). Families were paid 80 euros for their participation.

## Justification of fMRI sample size

The relation between children's literacy skills and home literacy practices typically ranges between $r = 0.40$ to $r = 0.60^{30,60}$. This range is in line with two previous studies that investigated the relation between home literacy practices and language-related (phonological processing and story listening) brain activity in young children. Specifically, Powers et al.<sup>30</sup> and Romeo et al.<sup>28</sup> ($n = 36$) found effect sizes of this relation ranging from $r = 0.40$ and $r = 0.58$ (see Table 5 in Powers et al.,<sup>30</sup> and openly available dataset in Romeo et al.<sup>28</sup>). Considering these estimates, power analyses using G* Power version 3.1<sup>61</sup> indicate that our final fMRI sample size of $n = 44$ would provide at least 88% power to detect correlations between home literacy practices and reading-related activity in the expected direction ($\alpha = 0.05$, one-tailed). Because the present study also aims at investigating whether reading-related brain activity mediates the relation between home literacy practices and literacy skills, we performed a power analysis of such a mediation using the exact parameters reported by Romeo et al.<sup>28</sup> in their mediation analysis involving language-related brain activity, conversational turns, and literacy skills. Monte Carlo simulations implemented in the online tool developed by Schoemann, Boulton, and Short<sup>62</sup> (https://schoemanna.shinyapps.io/mc_power_med/) indicated a power of 86% to detect an influence of the mediator (i.e., brain activity) on the relation between the predictor (i.e., home literacy practices) and the outcome (i.e., literacy skills) based on Romeo et al.'s estimates (1000 iterations, 95% Confidence level). Overall, although it is important to note that this effect size is likely inflated because of publication bias<sup>63</sup>, our final fMRI sample size appears to be adequately powered to detect the hypothesized relations.
## Assessment of child and parent literacy skills

Children's reading and vocabulary skills were assessed using the Alouette-R test³⁶ and the "vocabulaire" subtest of the NEMI-2³⁷. The Alouette-R test is a standardized test that is frequently used to test reading skills in French-speaking countries. This test requires participants to read a 265-word text aloud in 3 min. The number of words read and the number of pronunciation errors are used to calculate indices of reading speed and reading accuracy, respectively. The test has notably been shown to be a sensitive and specific screening tool detecting individuals with reading difficulties⁴⁴. The NEMI-2 is a standardized IQ test that is similar to the Stanford-Binet Intelligence Scales. It uses measures of verbal intelligence (i.e., "general knowledge", "vocabulary", and "comparison" subtests) and matrix reasoning (i.e., "Raven's matrices" subtest) to provide a composite IQ score that is highly correlated ($r = 0.80$) with the WISC score³⁷. The "vocabulaire" subtest of the NEMI-2 requires children to define up to 27 words that decrease in frequency. The subtest is untimed and is stopped after three consecutive errors. Age-normalized reading and vocabulary scores are presented in Table 2. Parental reading skills were also tested using the Alouette-R test.

## Parental questionnaire

The frequency of home literacy practices was assessed using an electronic questionnaire given to parents on a tablet. After providing information about SES, parents were asked how often they engaged in home literacy activities with their children. Questions also concerned practices involving other academic and non-academic domains (e.g., math, music, biology). These questions, which helped reduce the focus on literacy practices, are considered fillers in the present study (items about math practices are analyzed in Girard et al.⁵⁸). Practices were adapted from the questionnaire used by Skwarchuk and colleagues³⁸. Using a six-point rating scale (Did not occur/Activity is not relevant to my child, 1–3 times per month, Once per week, 2–4 times per week, almost daily, daily), parents rated the frequency with which they engaged in that practice with their child at home during the past month. Parents also had the option to indicate whether the practice occurred in the past but no longer at the time of testing (rated 1 point) and whether the child engaged in the practice at home but on her own (also rated 1 point).

Among the 77 questions, 18 concerned home literacy practices. These included a mixture of informal and formal practices (Fig. 2). As defined in Sénéchal and LeFevre (p. 1552⁵), "informal literacy experiences are those where the print is present but is not the focus of the parent–child interaction. In contrast, formal literacy activities are those where the attention is on the print itself⁶⁰. Formal practices were of varying difficulty, encompassing both basic or more advanced activities for an 8-year-old child. For example, listening when children read out loud was considered relatively basic but asking questions when children read was considered more advanced. Because some advanced activities were inherently related to more basic activities (e.g., reading/writing short texts and reading/writing long texts), those advanced items were presented to parents only if they stated that the corresponding more basic activity had occurred at home. This was done to avoid parents having to repeatedly say no to an activity they had never engaged in). Ratings associated with frequencies of home literacy practices are shown in Fig. 2. Frequencies of activities were averaged across items within each category to gather separate scores of informal, formal basic, and formal advanced practices. We also used a composite score, corresponding to the sum of these three subscores. Finally, parents were asked to provide a subjective estimate of their child's skills in reading (and other domains) using a six-point rating scale (i.e., Not sure, Severe difficulty, Difficulty, Average skills, Good skills, Very good skills).

## Experimental tasks

We adapted from Experiment 2B of Perrachione et al.³¹ a task in which a series of words were passively presented in blocks at the center of the screen (see Fig. 1b). As in Perrachione et al.³¹, words were monosyllabic nouns of three to five letters in length (mode = 4). The only difference between our task and that of Perrachione et al.³¹ was that our stimuli were in French instead of English. Words were thematically varied and equated for frequency using the CHACQFAM database in http://www.lexique.org/shiny/openlexicon/. In another version of the task that served as a control in the present study, words were replaced by Arabic digits ranging from 1 to 8. Word and digit adaptation tasks were each composed of two runs. In each of these runs, participants were presented with adaptation and

no-adaptation blocks. Adaptation blocks consisted of the repetition of the same stimulus (word or digit) eight times. No-adaptation blocks consisted of the presentation of eight different stimuli. Another task involving dot arrays was included in the scanning session but is not analyzed here.

## fMRI data acquisition

Images were collected using a Siemens Prisma 3T MRI scanner with a transmit body coil and a 64-channel receiver head-neck coil (Siemens Healthcare, Erlangen, Germany) at the CERMEP Imagerie du vivant in Lyon, France. The BOLD signal was measured with a susceptibility-weighted single-shot EPI sequence. Imaging parameters were as follows: TR = 2000 ms, TE = 24 ms, flip angle = 80°, matrix size = 128 × 120, field of view = 220 × 206 mm, slice thickness = 3 mm (0.48 mm gap), number of slices = 32. A high-resolution T1-weighted whole-brain anatomical volume was also collected for each participant. Parameters were as follows: TR = 3500 ms, TE = 2.24 ms, flip angle = 8°, matrix size = 256 × 256, field of view = 224 × 224 mm, slice thickness = 0.9 mm, number of slices = 192.

## Experimental timeline

The experimental timeline, which was directly taken from Perrachione et al.³¹, was identical in both the word and digit adaptation task. In each block, stimuli remained on the screen for 700 ms, with a 500 ms inter-stimulus interval (for a total block duration of 9.6 s). Ten adaptation blocks and ten no-adaptation blocks were presented along with ten blocks of visual fixation (duration = 9.6 s) in each run. Block presentation was pseudo-randomized such that two blocks of the same type could not follow each other. Finally, ten target stimuli (a picture of a rocket) randomly appeared in each run (outside of blocks). Participants were asked to press a button every time this target appeared.

## fMRI data preprocessing

Images were analyzed with SPM12 (Wellcome Department of Cognitive Neurology, London, UK). The first four images of each run were discarded to allow for T1 equilibration effects. Functional images were corrected for slice acquisition delays and spatially realigned to the first image of the first run to correct for head movements. Realigned images were smoothed with a Gaussian filter (4 mm × 4 mm × 7 mm full-width at half maximum). Using ArtRepair (https://cibsr.stanford.edu/tools/human-brain-project/artrepair-software.html), functional volumes with a global mean intensity greater than three standard deviations from the average of the run or a volume-to-volume motion greater than 2 mm were identified as outliers and substituted by the interpolation of the two nearest non-repaired volumes²⁸. Participants with outliers in more than 20% of volumes were excluded from the analyses ($n = 7$). After outlier exclusion, the movement range was on average 0.01 (SD = 0.04), 0.17 (SD = 0.30), and 0.18 (SD = 0.64) mm in the x, y, and z direction, with 0.27 (SD = 0.87), 0.37 (SD = 1.17), 0.45 (SD = 1.46) degrees of roll, pitch, and yaw.

Finally, functional images were normalized into the same stereotaxic space. Studies have found that anatomical differences between children older than 7–8-year-olds and adults are small enough that they are beyond the resolution of fMRI experiments⁶⁴,⁶⁵. Therefore, considering the age of our participants, the resolution of our data, and the fact that we wanted to be able to compare our results to that from Perrachione et al.³¹ with adult participants, we normalized all individual brains into the standard adult Montreal Neurological Institute (MNI) space. This was done in two steps. First, after coregistration with the functional data, the structural image was segmented into gray matter, white matter, and cerebrospinal fluid by using a unified segmentation algorithm⁶⁶. Second, the functional data were normalized to the MNI space by using the normalization parameters estimated during unified segmentation (normalized voxel size, $2 mm^3 \times 2 mm^3 \times 3.5 mm^3$).

## fMRI data analysis

Statistical analysis of fMRI data was performed according to the GLM. Brain activity associated with periods of adaptation and no-adaptation was modeled as epochs with onsets time-locked to the beginning of each block and a duration of 9.6 s. All epochs were convolved with a canonical hemodynamic response function. The time-series data were high-pass filtered (1/128 Hz), and serial correlations were corrected using an auto-regressive AR(1) model.

For each participant, the word adaptation effect was identified by subtracting activity associated with adaptation blocks from activity associated with no-adaptation blocks. Individual contrasts were then
submitted to one-sample *t*-tests across all participants. First, we used a whole-brain voxelwise approach to evaluate whether the overall brain network involved in word adaptation was similar to that found in Perrachione et al.³¹. An FDR-corrected cluster-level threshold of *p* < 0.05 (defined using voxel-level thresholds of *p* < 0.001 and *p* < 0.002) was applied to the whole-brain statistical map to assess brain activations. Second, we used an a priori region of interest (ROI) approach to measure activity in the seven brain regions that showed a word adaptation effect (i.e., lower activity in adaptation than no-adaptation blocks of the word adaptation task) in normal adults in Perrachione et al.³¹ (see Fig. 3). All ROIs were 6-mm radius spheres centered on coordinates reported in Table S5 for Experiment 2B of Perrachione et al.³¹. We tested for the presence of a word adaptation effect in each ROI using a series of one-sample *t*-tests. *P* values were Bonferroni-corrected for multiple comparison. In ROIs for which there was a significant word adaptation effect, the neural adaptation effect was then correlated with the frequency of home literacy practices. Mediation analyses were then performed to test whether, in regions for which there was a relation between home literacy practices and word adaptation, there was a mediating effect of vocabulary. These analyses were performed using the M3 toolbox in Matlab (https://github.com/canlab/MediationToolbox⁶⁷). Briefly, a mediation estimates the degree to which a variable (e.g., *M*) can explain the relation between two other variables (e.g., *X* and *Y*). Here, *X* was defined as the frequency of home literacy practices. *Y* was defined as the word adaptation effect. *M* was defined as vocabulary. Path *A* reflected the relation between home literacy practices and vocabulary. Path *B* reflected the relation between vocabulary and word adaptation, controlling for home literacy practices. The total relationship between frequency of home literacy practices and word adaptation (including direct and indirect effects) is reflected in path *c*, while the direct effect of the relationship between frequency of home literacy practices and word adaptation (controlling for vocabulary) is reflected in path *c'*. Finally, the significance of the mediator was tested in the product *A*\*b⁶⁷⁻⁶⁹. A bias-corrected bootstrap test over 10,000 iterations was used to measure the statistical significance of the mediator. Mediation analyses systematically included as covariates parental education and income to control for differences in SES between children.

Because our hypotheses were directional (i.e., we only anticipated lesser activity in adaptation than no-adaptation blocks, as well as increases in either literacy skills or neural adaptation effects as a function of home literacy practices), *p* values for behavioral and fMRI results are one-sided (unless otherwise noted). Accordingly, only 95% lower bound confidence intervals (95% CI Lower) are reported.

## Reporting Summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

## DATA AVAILABILITY

The task, electronic version of the questionnaire (in French), behavioral data, and individual MRI data are available from Zenodo (https://doi.org/10.5281/zenodo.5112814). The ROI images and the whole-brain unthresholded p-maps corresponding to Fig. 4a and Fig. 7 are available from NeuroVault (https://identifiers.org/neurovault.collection:9453).

## CODE AVAILABILITY

The Matlab scripts used to analyze fMRI data are available from GitHub (https://github.com/BBL-lab/BBL-batch-system). The M3 toolbox is also available from GitHub (https://github.com/canlab/MediationToolbox).

Received: 25 February 2021; Accepted: 29 October 2021;
Published online: 03 December 2021

## REFERENCES

1. Storch, S. A. & Whitehurst, G. J. Oral language and code-related precursors to reading: evidence from a longitudinal structural model. *Dev. Psychol.* **38**, 934 (2002).

2. Janus, M. & Duku, E. The school entry gap: socioeconomic, family, and health factors associated with children's school readiness to learn. *Early Educ. Dev.* **18**, 375-403 (2007).

3. Hamilton, L. G., Hayiou-Thomas, M. E., Hulme, C. & Snowling, M. J. The home literacy environment as a predictor of the early literacy development of children at family-risk of dyslexia. *Sci. Stud. Read.* **20**, 401–419 (2016).

4. Sénéchal, M. Testing the home literacy model: parent involvement in kindergarten is differentially related to grade 4 reading comprehension, fluency, spelling, and reading for pleasure. *Sci. Stud. Read.* **10**, 59–87 (2006).

5. Sénéchal, M. & LeFevre, J.-A. Continuity and change in the home literacy environment as predictors of growth in vocabulary and reading. *Child Dev.* **85**, 1552–1568 (2014).

6. Pan, B. A., Rowe, M. L., Singer, J. D. & Snow, C. E. Maternal correlates of growth in toddler vocabulary production in low-income families. *Child Dev.* **76**, 763–782 (2005).

7. Weizman, Z. O. & Snow, C. E. Lexical input as related to children's vocabulary acquisition: effects of sophisticated exposure and support for meaning. *Dev. Psychol.* **37**, 265–279 (2001).

8. Molfese, V. J., Modglin, A. & Molfese, D. L. The role of environment in the development of reading skills: a longitudinal study of preschool and school-age measures. *J. Learn. Disabil.* **36**, 59–67 (2003).

9. Sénéchal, M. & Young, L. The effect of family literacy interventions on children's acquisition of reading from kindergarten to grade 3: a meta-analytic review. *Rev. Educ. Res.* **78**, 880–907 (2008).

10. Hood, M., Conlon, E. & Andrews, G. Preschool home literacy practices and children's literacy development: a longitudinal analysis. *J. Educ. Psychol.* **100**, 252–271 (2008).

11. Relyea, J. E., Zhang, J., Liu, Y. & Wui, M. G. L. Contribution of home language and literacy environment to English reading comprehension for emergent bilinguals: sequential mediation model analyses. *Read. Res. Q.* **55**, 473–492 (2020).

12. Ouellette, G. & Beers, A. A not-so-simple view of reading: How oral vocabulary and visual-word recognition complicate the story. *Read. Writ.* **23**, 189–208 (2010).

13. Council on Early Childhood, High, P. C. & Klass, P. Literacy promotion: an essential component of primary care pediatric practice. *Pediatrics* **134**, 404–409 (2014).

14. World Health Organization. Guidelines on Physical Activity, Sedentary Behaviour and Sleep for Children Under 5 Years of Age: Web Annex: Evidence Profiles (WHO, 2019).

15. Hoff, E. How social contexts support and shape language development. *Dev. Rev.* **26**, 55–88 (2006).

16. Noble, K. G., Wolmetz, M. E., Ochs, L. G., Farah, M. J. & McCandliss, B. D. Brain-behavior relationships in reading acquisition are modulated by socioeconomic factors. *Dev. Sci.* **9**, 642–654 (2006).

17. Hutton, J. S. et al. Validation of the reading house and association with cortical thickness. *Pediatrics* **147**, e20201641 (2021).

18. Jednoróg, K. et al. The influence of socioeconomic status on children's brain structure. *PLoS ONE* **7**, e42486 (2012).

19. Mackey, A. P. et al. Neuroanatomical correlates of the income-achievement gap. *Psychol. Sci.* **26**, 925–933 (2015).

20. Noble, K. G. et al. Socioeconomic disparities in neurocognitive development in the first two years of life. *Dev. Psychobiol.* **57**, 535–551 (2015).

21. Raizada, R. D. S., Richards, T. L., Meltzoff, A. & Kuhl, P. K. Socioeconomic status predicts hemispheric specialisation of the left inferior frontal gyrus in young children. *NeuroImage* **40**, 1392–1401 (2008).

22. Noble, K. G., Farah, M. J. & McCandliss, B. D. Socioeconomic background modulates cognition-achievement relationships in reading. *Cogn. Dev.* **21**, 349–368 (2006).

23. Hutton, J. S., Dudley, J., Horowitz-Kraus, T., DeWitt, T. & Holland, S. K. Associations between home literacy environment, brain white matter integrity and cognitive abilities in preschool-age children. *Acta Paediatrica* **109**, 1376–1386 (2020).

24. Hutton, J. S. et al. Home reading environment and brain activation in preschool children listening to stories. *Pediatrics* **136**, 466–478 (2015).

25. Hutton, J. S. et al. Shared reading quality and brain activation during story listening in preschool-age children. *J. Pediatrics* **191**, 204–211. e1 (2017).

26. Merz, E. C., Maskus, E. A., Melvin, S. A., He, X. & Noble, K. G. Socioeconomic disparities in language input are associated with children's language-related brain structure and reading skills. *Child Dev.* **91**, 846–860 (2020).

27. Romeo, R. R. et al. Language exposure relates to structural neural connectivity in childhood. *J. Neurosci.* **38**, 7870–7877 (2018).

28. Romeo, R. R. et al. Beyond the 30-million-word gap: children's conversational exposure is associated with language-related brain function. *Psychol. Sci.* **29**, 700–710 (2018).

29. Su, M. et al. Influences of the early family environment and long-term vocabulary development on the structure of white matter pathways: a longitudinal investigation. *Dev. Cogn. Neurosci.* **42**, 100767 (2020).

30. Powers, S. J., Wang, Y., Beach, S. D., Sideridis, G. D. & Gaab, N. Examining the relationship between home literacy environment and neural correlates of phonological processing in beginning readers with and without a familial risk for dyslexia: an fMRI study. *Ann. Dyslexia* **66**, 337–360 (2016).
31. Perrachione, T. K. et al. Dysfunction of rapid neural adaptation in dyslexia. *Neuron* **92**, 1383-1397 (2016).

32. Grill-Spector, K. & Malach, R. fMR-adaptation: a tool for studying the functional properties of human cortical neurons. *Acta Psychol.* **107**, 293-321 (2001).

33. Bossier, H. et al. The empirical replicability of task-based fMRI as a function of sample size. *Neuroimage* **212**, 116601 (2020).

34. Saxe, R., Brett, M. & Kanwisher, N. Divide and conquer: a defense of functional localizers. *Neuroimage* **30**, 1088-1096 (2006). discussion 1097-1099.

35. Amalric, M. & Dehaene, S. Origins of the brain networks for advanced mathematics in expert mathematicians. *Proc. Natl Acad. Sci. USA* **113**, 4909-4917 (2016).

36. Lefavrais, P. Test de l'Alouette, version révisée. Paris: ECPA (2005).

37. Cognet, G. NEMI-2: Nouvelle Echelle Métrique de l'Intelligence-2. (2006).

38. Skwarchuk, S.-L., Sowinski, C. & LeFevre, J.-A. Formal and informal home learning activities in relation to children's early numeracy and literacy skills: the development of a home numeracy model. *J. Exp. Child Psychol.* **121**, 63-84 (2014).

39. Olulade, O. A., Flowers, D. L., Napoliello, E. M. & Eden, G. F. Dyslexic children lack word selectivity gradients in occipito-temporal and inferior frontal cortex. *Neuroimage Clin.* **7**, 742-754 (2015).

40. Szwed, M., Qiao, E., Jobert, A., Dehaene, S. & Cohen, L. Effects of literacy in early visual and occipitotemporal areas of Chinese and French readers. *J. Cogn. Neurosci.* **26**, 459-475 (2014).

41. van Bergen, E., van Zuijen, T., Bishop, D. & de Jong, P. F. Why are home literacy environment and children's reading skills associated? What parental skills reveal. *Read. Res. Q.* **52**, 147-160 (2017).

42. Hart, S. A., Little, C. & van Bergen, E. Nurture might be nature: cautionary tales and proposed solutions. *npj Sci. Learn.* **6**, 1-12 (2021).

43. Hart, B. & Risley, T. R. *Meaningful Differences in the Everyday Experience of Young American Children* (Paul H Brookes, 1995).

44. Cavalli, E. et al. Screening for dyslexia in French-speaking university students: an evaluation of the detection accuracy of the Alouette test. *J. Learn. Disabil.* **51**, 268-282 (2018).

45. de Jong, P. F. & Leseman, P. P. M. Lasting effects of home literacy on reading achievement in school. *J. Sch. Psychol.* **39**, 389-414 (2001).

46. Perfetti, C. Reading ability: lexical quality to comprehension. *Sci. Stud. Read.* **11**, 357-383 (2007).

47. Price, C. J., Wise, R. J. S. & Frackowiak, R. S. J. Demonstrating the implicit processing of visually presented words and pseudowords. *Cereb. cortex* **6**, 62-70 (1996).

48. Turkeltaub, P. E., Gareau, L., Flowers, D. L., Zeffiro, T. A. & Eden, G. F. Development of neural mechanisms for reading. *Nat. Neurosci.* **6**, 767-773 (2003).

49. Elliott, M. L. et al. What is the test-retest reliability of common task-functional MRI measures? New empirical evidence and a meta-analysis. *Psychol. Sci.* **31**, 792-806 (2020).

50. Martin, A., Schurz, M., Kronbichler, M. & Richlan, F. Reading in the brain of children and adults: a meta-analysis of 40 functional magnetic resonance imaging studies. *Hum. Brain Mapp.* **36**, 1963-1981 (2015).

51. Dehaene-Lambertz, G., Monzalvo, K. & Dehaene, S. The emergence of the visual word form: longitudinal evolution of category-specific ventral visual areas during reading acquisition. *PLoS Biol.* **16**, e2004103 (2018).

52. Taylor, J. S. H., Rastle, K. & Davis, M. H. Can cognitive models explain brain activation during word and pseudoword reading? A meta-analysis of 36 neuroimaging studies. *Psychol. Bull.* **139**, 766 (2013).

53. Binder, J. R. & Desai, R. H. The neurobiology of semantic memory. *Trends Cogn. Sci.* **15**, 527-536 (2011).

54. Malins, J. G. et al. Individual differences in reading skill are related to trial-by-trial neural activation variability in the reading network. *J. Neurosci.* **38**, 2981-2989 (2018).

55. Inoue, T., Georgiou, G. K., Parrila, R. & Kirby, J. R. Examining an extended home literacy model: the mediating roles of emergent literacy skills and reading fluency. *Sci. Stud. Read.* **22**, 273-288 (2018).

56. King, M. F. & Bruner, G. C. Social desirability bias: a neglected aspect of validity testing. *Psychol. Mark.* **17**, 79-103 (2000).

57. Noble, K. G. & Giebler, M. A. The neuroscience of socioeconomic inequality. *Curr. Opin. Behav. Sci.* **36**, 23-28 (2020).

58. Girard, C. et al. Nurturing the mathematical brain: home numeracy practices are associated with children's neural responses to Arabic numerals. *Psychol. Sci.* (in the press).

59. Robin, M. Davantage de personnes vivent dans un ménage médian dans les régions de l'Ouest. *Insee Focus* (2019).

60. Bus, A. G., van IJzendoorn, M. H. & Pellegrini, A. D. Joint book reading makes for success in learning to read: a meta-analysis on intergenerational transmission of literacy. *Rev. Educ. Res.* **65**, 1-21 (1995).

61. Faul, F., Erdfelder, E., Lang, A.-G. & Buchner, A. G* power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences. *Behav. Res. Methods* **39**, 175-191 (2007).

62. Schoemann, A. M., Boulton, A. J. & Short, S. D. Determining power and sample size for simple and complex mediation models. *Soc. Psychol. Personal. Sci.* **8**, 379-386 (2017).

63. Anderson, S. F., Kelley, K. & Maxwell, S. E. Sample-size planning for more accurate statistical power: a method adjusting sample effect sizes for publication bias and uncertainty. *Psychol. Sci.* **28**, 1547-1562 (2017).

64. Burgund, E. D. et al. The feasibility of a common stereotactic space for children and adults in fMRI studies of development. *Neuroimage* **17**, 184-200 (2002).

65. Kang, H. C., Burgund, E. D., Lugar, H. M., Petersen, S. E. & Schlaggar, B. L. Comparison of functional activation foci in children and adults using a common stereotactic space. *Neuroimage* **19**, 16-28 (2003).

66. Ashburner, J. & Friston, K. J. Unified segmentation. *NeuroImage* **26**, 839-851 (2005).

67. Wager, T. D., Davidson, M. L., Hughes, B. L., Lindquist, M. A. & Ochsner, K. N. Prefrontal-subcortical pathways mediating successful emotion regulation. *Neuron* **59**, 1037-1050 (2008).

68. Wager, T. D. et al. Brain mediators of cardiovascular responses to social threat: Part I: reciprocal dorsal and ventral sub-regions of the medial prefrontal cortex and heart-rate reactivity. *NeuroImage* **47**, 821-835 (2009).

69. Shrout, P. E. & Bolger, N. Mediation in experimental and nonexperimental studies: new procedures and recommendations. *Psychol. Methods* **7**, 422-445 (2002).

## ACKNOWLEDGEMENTS

We thank Franck Lamberton and Danielle Ibarrola for their help with data collection. We are also grateful to the children and parents who participated. This study was funded by the Agence Nationale de la Recherche (ANR-14-CE30-0002 and ANR-17-CE28-0014) to J.P.

## AUTHOR CONTRIBUTIONS

C.G. and J.P. developed the study concept, designed the study, and analyzed and interpreted the data. C.G., T.B., J.L., J.E.-D., and L.L. collected the data. C.G. drafted the manuscript. T.B., J.L., J.E.-D., L. L., and J.P. provided critical revisions. All authors approved the final version of the manuscript for submission.

## COMPETING INTERESTS

The authors declare no competing interests.

## ADDITIONAL INFORMATION

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41539-021-00112-9.

**Correspondence** and requests for materials should be addressed to Cléa. Girard or Jérôme Prado.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons

Attribution 4.0 International License, which permits use, sharing,

adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.

org/licenses/by/4.0/.

© The Author(s) 2021