Children With Dyslexia and Familial Risk for Dyslexia Present Atypical Development of the Neuronal Phonological Network

Learning to read changes the brain language system. Phonological processing is the language domain most crucial for reading, but it is still unknown how reading acquisition modifies the neural phonological network in children who either develop dyslexia or are at risk of dyslexia. For the two first years of formal education, we followed 90 beginning readers with (n = 55) and without (n = 35) familial risk of dyslexia who became typical readers (n = 70) or developed dyslexia (n = 20). We used functional magnetic resonance imaging to identify the neural correlates of phonological awareness using an auditory rhyme judgment task. This task was applied when participants were starting formal education, and repeated 2 years later. By applying two alternative group splits, we analyzed the effects of dyslexia and the effects of familial risk of dyslexia separately. We found that the phonological brain network undergoes reorganization during the first 2 years of formal education. This process proceeds differently depending on the presence of a familial history of dyslexia and reading impairment. Typical readers without risk for dyslexia activate structures responsible for phonological processing already at the beginning of literacy. This group shows reduced brain activation over time during phonological processing, perhaps due to automatization of phonological skills. Children who develop reading impairment present a delay in the development of phonological structures such as the bilateral superior temporal gyri, left middle temporal gyrus, right insula and right frontal cortex, where we observed time and group interaction. Finally, typical readers with familial risk of dyslexia also present an atypical development of the neural phonological structures, visible both at the beginning of reading instruction and 2 years later. These children used a presumably efficient neural mechanism of phonological processing, based on the activation of the precentral and postcentral gyri, and achieved a typical level of phonological awareness.

Learning to read changes the brain language system. Phonological processing is the language domain most crucial for reading, but it is still unknown how reading acquisition modifies the neural phonological network in children who either develop dyslexia or are at risk of dyslexia. For the two first years of formal education, we followed 90 beginning readers with (n = 55) and without (n = 35) familial risk of dyslexia who became typical readers (n = 70) or developed dyslexia (n = 20). We used functional magnetic resonance imaging to identify the neural correlates of phonological awareness using an auditory rhyme judgment task. This task was applied when participants were starting formal education, and repeated 2 years later. By applying two alternative group splits, we analyzed the effects of dyslexia and the effects of familial risk of dyslexia separately. We found that the phonological brain network undergoes reorganization during the first 2 years of formal education. This process proceeds differently depending on the presence of a familial history of dyslexia and reading impairment. Typical readers without risk for dyslexia activate structures responsible for phonological processing already at the beginning of literacy. This group shows reduced brain activation over time during phonological processing, perhaps due to automatization of phonological skills. Children who develop reading impairment present a delay in the development of phonological structures such as the bilateral superior temporal gyri, left middle temporal gyrus, right insula and right frontal cortex, where we observed time and group interaction. Finally, typical readers with familial risk of dyslexia also present an atypical development of the neural phonological structures, visible both at the beginning of reading instruction and 2 years later. These children used a presumably efficient neural mechanism of phonological processing, based on the activation of the precentral and postcentral gyri, and achieved a typical level of phonological awareness.

INTRODUCTION
Learning to read is a long-lasting process which requires mastering a number of skills. Phonological skills, i.e., understanding that spoken words are composed of discrete sounds which can be mapped onto letters, are particularly important, and a deficit in phonological skills is a common etiology of developmental dyslexia (Wagner and Torgesen, 1987;Ziegler et al., 2010). The brain network responsible for phonological processing includes mostly left hemisphere areas like the inferior parietal lobule (including supramarginal and angular gyri), inferior frontal cortex, postcentral and precentral gyri, superior and middle temporal gyri, and fusiform and dorsolateral prefrontal cortex (Poldrack et al., 1999;Booth et al., 2006;Hoeft et al., 2006;Bitan et al., 2007;Cao et al., 2008;Kovelman et al., 2012). Patterns of brain activation during phonological processing differ between dyslexic and typical readers and have been observed with functional magnetic resonance imaging (fMRI; Richlan et al., 2011). However, only a few fMRI experiments have explored phonological processing independent of reading abilities using a task based on spoken language (Desroches et al., 2010;Kovelman et al., 2012;Raschle et al., 2012Raschle et al., , 2014Dębska et al., 2016;Yu et al., 2018). Studies on dyslexic and typically reading children have produced inconsistent results: Desroches et al. (2010) reported hypoactivation of the left fusiform gyrus in the dyslexic group while Kovelman et al. (2012) found hypoactivation only in the left dorsolateral prefrontal cortex.
In order to explore patterns of brain activation independent of reading experience and already present before schooling, Raschle et al. (2014) recruited pre-reading children with and without a familial history of dyslexia. Children with familial history of dyslexia constitute a common risk group and tend to have around a 45% chance to develop a reading disorder (Snowling and Melby-Lervåg, 2016). Children with a familial history of dyslexia exhibited underactivation during phonological processing in the bilateral occipito-temporal and left temporoparietal areas, as well as in the cerebellum (Raschle et al., 2014). They also showed disruption in the left frontal brain regions during auditory processing (Raschle et al., 2014). Dębska et al. (2016), who tested Polish-speaking pre-schoolers and first-graders, reported similar results. Children with a familial history of dyslexia showed reduced activation of many brain regions, including the left occipito-temporal cortex. For this particular region, reduced activation was observed only among pre-schoolers in the familial risk group. At the same time, a grade effect of lower activity in first-graders compared to kindergarten pupils was observed exclusively in children without familial history of dyslexia. This cross-sectional study raises the question as to what the connection is between familial history of dyslexia and formal reading instruction. It may be that the phonological network engaged by preschoolers is not essential for phonological processing during the time when both phonological awareness and reading are becoming more fluent. As a result, any differences between children with and without familial risk may only be visible at the pre-reading state. However, it is still not clear how development and literacy acquisition affect differences in phonological processing between dyslexic and typically reading children, and whether familial risk for dyslexia plays a role in this process.
Literacy affects the brain language network (Dehaene et al., 2015) since reading acquisition supports the development of more advanced phonemic skills through experience with an alphabet (Bentin, 1992). Phonological processing is also affected by learning to read. The level of reading skills at the beginning of schooling is a predictor of growth in phonological awareness in the following year (Nation and Hulme, 2011). Reading skills support the development of subvocal rehearsal which, in turn, improves performance in verbal working memory tasks (Demoulin and Kolinsky, 2016). However, it is still unknown how reading development affects neural correlates of auditory phonological processing, even though phonology is a language domain crucial for reading. Two cross-sectional studies provided evidence for age-related increases of brain activation in the phonological processing network in alphabetic languages (Cone et al., 2008;Brennan et al., 2013). Specifically, English-speaking typical readers during rhyme judgment showed increased activation with age in the left dorsal inferior frontal and temporal gyri (Cone et al., 2008;Brennan et al., 2013) and inferior parietal cortex (Brennan et al., 2013). However, in cross-sectional studies, age-related effects may also come from existing between-group differences independent of age or reading experience. So far only one study has used a longitudinal design to explore the development of the phonological network in typical readers from the beginning of formal education until the emergent reading stage (5-7-9 years of age, Yu et al., 2018). They found a developmental decrease in activation in the left inferior parietal cortex and bilateral precuneus. This decrease suggests that experience in reading induced specialization of brain regions responsible for phonological processing. As participants were restricted to typical readers without familial risk of dyslexia, it is unknown how the phonological network develops in atypical reading development where children have a familial history of dyslexia or reading disabilities.
Here, we aim to longitudinally explore changes in the phonological processing network that occur when children with and without familial risk for dyslexia learn to read, and to follow these changes as they either become proficient readers or develop dyslexia. First, we explore the development of the neural basis of phonological processing in children who acquire reading typically and children who develop dyslexia during the first 2 years of their school education. We hypothesize that typical readers will show reduced brain activity in the language processing network in later stages of reading development as compared with the early stage, in-line with experience-induced fine-tuning reported previously (Yu et al., 2018). We expect children with dyslexia to display behavioral and brain activation differences as compared to typical readers at both early and later stages of reading development. In particular, we expect children with dyslexia to present low accuracy in both reading and phonological assessments at all measurement points and to show hypoactivation of left hemisphere structures (as in Desroches et al., 2010;Kovelman et al., 2012).
Second, by applying alternative group splitting, we examine the impact of familial history of dyslexia on the neural correlates of phonological processing during the same time range of the first years of education. We expect that the phonological network of children with familial risk for dyslexia who become typical readers should differ from that of typical readers without familial history of dyslexia at the beginning of literacy acquisition. Those children with a family history of dyslexia that become typical readers should show a vast hypoactivation of the typical phonological network (as in Raschle et al., 2014;Dębska et al., 2016). We hypothesize that a familial history of dyslexia has a persistent impact on brain function, even in typically reading children. Therefore, we expect to find the effects of familial risk not only at the beginning of education, but also in later stages of reading development. Presumably, the effects of familial history of dyslexia will be observed independent of dyslexia itself.

Participants
All participants took part in a longitudinal study approved by the Warsaw University Ethical Committee. The study consisted of three time points, a year apart from each other. The first and the third time points involved both behavioral and fMRI sessions. The second time point was limited to behavioral testing.
At the first time point, 120 native Polish-speaking children were recruited from first grade and kindergarten classes. Due to the educational reform taking place in Poland during our study, children could begin formal school education at age six or seven, depending on their parents' choice. Therefore, the age range of children who attended first grade in September was 5.9 (if a child was born in December and her/his parents decided that she/he would begin schooling as 6-year-old) to 7.8 (if a child was born in January and her/his parents decided that she/he would begin schooling as 7-year-old). This massive variation of age of first-graders is reflected in our sample. Similarly, there was a variation in age among children who attended the last year of kindergarten, who were on average 6 months younger than the first-graders in our sample. Because of the relatively small difference in age between the two groups, we decided to recruit children from both the first grade and kindergarten. Though formal reading training was supposed to start in elementary school, children were already taught letters in kindergarten.
The participants had typical or high IQs (higher than 25. percentile as measured with Raven Matrices), were monolingual, right-handed (as reported by parents and confirmed by experimenters during testing, when children were supposed to use the hand they prefer for writing), and born at term (≥37 weeks of pregnancy). None of them had any history of neurological illness or brain damage and no prior symptoms of ADHD as reported by their parents. From the group of 102 participants in our previous study (Dębska et al., 2016), nine children either left the study due to moving to another city or losing interest in the project. Additionally, data from three children (the oldest one and the two youngest ones) were excluded from analyses in order to clearly separate the age range of participants during the first and the third time points.
The current study included 90 children (53 girls and 37 boys) aged 5.94-7.95 years (6.91 years on average) at the first time point, 6.99-9.43 (7.88 years on average) at the second time point, and 8.05-10.05 years (8.95 years on average) at the third time point. At the first time point, 27 participants were attending the last year of kindergarten and 63 were first graders. At consecutive time points, children progressed in their education until, at the third time point, there were 27 second graders and 63 third graders. The distribution of children who attended school and kindergarten did not differ between children with a familial history of dyslexia and without. There was also no difference in distribution between children diagnosed later with dyslexia and typical readers. Therefore, the results were affected by the school grade to the same extent in all subgroups.
There were 55 children with familial history of dyslexia and 35 children without familial history of dyslexia. Children with a familial history of dyslexia were identified as those who had a first-degree relative with a formal diagnosis of developmental dyslexia or at least one parent who scored more than 0.4 points in the Adult Reading History Questionnaire (ARHQ, Lefly and Pennington, 2000), as specified in Black et al. (2012). These criteria overlapped in 32 children while 23 children fulfilled only the questionnaire criterion. There were no children who had a first-degree relative with a diagnosis of dyslexia and for which both parents scored below 0.4 points in the questionnaire. The two groups of children whose parents scored high in ARHQ and had a first-degree relative with diagnosis of dyslexia (maternal ARHQ: M = 0.35, SD = 0.15, paternal ARHQ: M = 0.45, SD = 0.15) and children whose parents scored high in ARHQ but had no formal diagnosis of dyslexia in family (maternal ARHQ: M = 0.42, SD = 0.12, paternal ARHQ: M = 0.38, SD = 0.16) did not differ in maternal or paternal ARHQ average scores (maternal: t(53) = 1.80, p = 0.078, 95% CI DV = [−0.01; 0.14], d = 0.50; paternal: t(48) = 1.68, p = 0.100, 95% CI DV = [−0.16; 0.015], d = 0.48). We used the questionnaire measure due to the fact that in Poland, reading impairment was not diagnosed when the parents of our participants were school aged. However, according to the norming study of the Polish version of ARHQ, the criterion of scoring more than 0.4 points is fulfilled by 11% of mothers and 16% of fathers (Krasowicz-Kupis et al., 2014).
The third time point included a formal diagnosis of dyslexia and made it possible to select a group of children with dyslexia (n = 20). Children selected as having a reading disorder achieved low scores (at least 1 SD below the population mean, which corresponds to below 16 percentile) in at least two reading subtests (out of four: single-word reading, pseudo-word reading, reading with lexical decision, and text reading). Children who achieved a low score in no more than one reading subtest were assigned to the typically reading group. Children assessed as dyslexic achieved low scores in 3.05 tests (out of four) on average (SD = 0.69), and children assessed as typically reading scored low in 0.23 tests (out of four) on average (SD = 0.43). The dyslexic group included 15 children with a familial history of dyslexia and 5 children without a familial history of dyslexia, whereas the group of typical readers (n = 70) included 40 children with a familial history of dyslexia and 30 children without such risk factor.
Children with dyslexia did not differ from typical readers in terms of age at each time point, sex, grade, or verbal and nonverbal IQ (as measured with WISC-R). Although their parental socioeconomic status, paternal score in the Adult Reading History Questionnaire, and performance IQ (measured with Raven Matrices) was slightly lower than in the typically reading group, these differences did not survive Bonferroni correction for multiple comparisons (see Table 1). Similar effects were observed when the groups were restricted to children with a familial history of dyslexia (see Supplementary Table S1).
Groups of children with and without familial history of dyslexia did not differ in age at each time point, sex, grade, parental socioeconomic status, or IQ. Children with a familial history of dyslexia had significantly higher maternal and paternal scores in the Adult Reading History Questionnaire than children without familial risk (see Table 1), as it was the group division criterion. Similar effects were observed when the groups were restricted to typically reading children (see Supplementary Table S1).

Experimental Design: Behavioral Measures
Participants completed three phases of behavioral tests which included measuring reading skills, letter knowledge, rapid naming, and phonological awareness (all time points), and language and cognitive skills (the first time point). Standard tests designed for dyslexia diagnosis were applied at the third time point. Parents of participants completed questionnaires about FIGURE 1 | Word reading, phoneme analysis, and phoneme deletion of typical readers (TR) and children with dyslexia (DR) across three time points.
Frontiers in Neuroscience | www.frontiersin.org  their reading history (ARHQ, Lefly and Pennington, 2000), as well as about their home literacy environment.
Letter knowledge, word and pseudo-word reading, elision and phoneme analysis were measured at each time point with the same test battery (Szczerbiński and Pelc-Pêkala, 2013). Rapid automatized naming was tested with subtests of object and color naming at the first time point, and object, color, letter, and digit naming at further time points (Fecenec et al., 2013). Language skills were assessed with the Picture Vocabulary Test: Comprehension (vocabulary assessment at the first time point; Haman et al., 2012), and selected subtests of Test of Language Development (vocabulary and grammar assessment at the second time point; Smoczyńska et al., 2015). Intelligence was assessed with Raven's Colored Progressive Matrices (Szustrowa and Jaworowska, 2003) at the first time point, and with Wechsler Intelligence Scale for Children -Revised (WISC-R; Matczak et al., 2008) at the second time point. Syllable and digit span were measured at the first time point with tasks in which series of increasing length were repeated by participants, and the total number of correctly repeated series was used as the outcome measure.
At the third time point, a standardized battery of tests for dyslexia diagnosis was applied (Bogdanowicz et al., 2009). The battery consisted of nine tests: four of them assessing reading, two assessing writing, and three measuring phonological skills. Reading tests included word reading, pseudoword reading, reading with lexical decision, and text reading with comprehension control. Writing was measured with word writing and text writing. Phonological abilities were assessed using phoneme deletion, pseudoword repetition, and a battery of tasks (phoneme analysis and synthesis, syllable analysis and synthesis) based on pseudowords.

Statistical Analyses
We divided the analyses into two parts: (A) concerning the effects of dyslexia and (B) concerning the effects of familial history of dyslexia. For each time point we compared performance in behavioral tests between children who developed typical reading skills and children with dyslexia for the whole sample. Next, we repeated the analyses for the sample restricted to children with familial history of dyslexia to further explore if differences between typical and dyslexic readers were driven by reading ability regardless of familial risk. Similarly, all comparisons were performed between children with and without familial history of dyslexia for the whole sample and then repeated for only typical readers. For comparisons between groups of unequal sizes (typical readers vs. children with dyslexia in the whole sample and in children with familial history of dyslexia; children with familial history of dyslexia vs. children without familial history of dyslexia) non-parametric methods were applied, whereas for comparison of typical readers with and without familial history of dyslexia, parametric tests were used. Each series of comparisons was followed by Bonferroni corrections.
The optimal way to compare the effects of presence of reading impairment and familial risk of dyslexia in one model, would be to run 2 × 2 comparisons. However, the small size of the group of children without familial risk of dyslexia, who developed reading impairment (n = 5) did not allow us to run this type of analyses. Therefore, we applied the same model of analyses as used before for research of white matter volume in a similar group of children with and without familial risk of dyslexia and with and without reading impairment (Vanderauwera et al., 2017). This involved analyzing the effects of dyslexia first in whole sample, then in children with familial history of dyslexia, and the effects of familial history of dyslexia first in whole sample, then in typical readers.

Experimental Design: fMRI Tasks
The same fMRI tasks were applied at the first and the third time points and have been described by Dębska et al. (2016). Twenty noun pairs were presented to both ears of participants via headphones. The noun pairs were accompanied by pictures depicting the words. Subjects viewed the stimuli on a backprojection screen through an angled mirror. The exact list of stimuli is provided as Supplementary Table S2. After each pair, children had to decide whether the words rhymed or not (Rhyme task), similar to Kovelman et al. (2012). The control experiment included exactly the same stimuli, but the participants' task was to assess whether the words were spoken by speakers of the same gender or not (Voice task), similar to Raschle et al. (2012). The yes/no decision was made by pressing the corresponding button.
Both tasks were contrasted with a rest condition. During the rest condition, children looked at a fixation cross for the duration of the block. The accuracy and reaction times were analyzed at both TPs using repeated-measures ANOVA. The procedure was the same for both time points. Children were familiarized with the task in a mock-scanner using items not included in the following scanning session. The experimental scheme consisted of two functional runs: one with the experimental Rhyme task and one with the control Voice task. The timing and duration of tasks were identical. The order of runs was counterbalanced between the children and reversed at the third time point (compared to the first time point). Each word was presented to children via headphones, and at the same time the picture of its meaning appeared on the screen for 2 s. After that, the second word was played and the second picture appeared for 2 s. Then, a question mark appeared for 2 s prompting the child to give an answer. Each run consisted of ten blocks: five blocks with stimulation, with four trials per block, and five with the rest condition, each lasting 24 s. Half of the trials matched regarding rhyme and half of the words were spoken by the same gender voice. Stimuli were presented using Presentation software (Neurobehavioral Systems).

Behavioral Results
Children with dyslexia had lower performance than typical readers in early reading and reading related tests, already at the first time point when children had only just started formal education (see Table 1 and Figure 1). Differences between these two groups were observed in letter knowledge (the first two time points), word and pseudoword reading (at each time point), phoneme analysis (only at the first time point due to ceiling effect at later time points), phoneme deletion (at each time point), rapid automatized naming of colors and objects (the first and the third time points), and rapid automatized naming of letters and digits (the second and the third time points). However, vocabulary and working memory performance measured with digit and syllable span of children with dyslexia did not differ from the typically reading group at the first time point. Additionally, in all tests from the normalized battery for dyslexia diagnosis, children with dyslexia scored significantly lower than typically reading children (see Table 2). This pattern of results was not qualitatively different when the group was restricted to children with a familial history of dyslexia (see Supplementary Table S1 for demographic and reading related tests and Supplementary Table S3 for dyslexia diagnosis tests).

In-Scanner Performance
Children with dyslexia underperformed typical readers in accuracy in the in-scanner Rhyme task, but only at the first time point (see Table 3). This difference was also significant when the group was restricted to children with a familial history of dyslexia (see Supplementary Table S4). The two groups did not differ either in accuracy in the Voice task, or reaction times in the in-scanner tasks.

fMRI Results
Figure 2 depicts brain activation to Rhyme > Voice (see Experimental Design: fMRI tasks) contrast in typical readers FIGURE 3 | Effects of dyslexia (typical readers vs. dyslexic readers) and familial history of dyslexia (group with familial history of dyslexia vs. children without familial history of dyslexia) at the first and the third time points as revealed by two-sample t-tests. and children with dyslexia at each time point (for details see Supplementary Table S5).
At the first time point, typical readers engaged widespread brain networks including the bilateral inferior frontal areas, the inferior, middle, superior, and anterior temporal areas, the left fusiform gyrus and calcarine sulcus, the cingulate cortex and putamen and caudate (subcortically). At the third time point, activation was restricted to the bilateral putamen, caudate, and occipital areas.
At the first time point, children with dyslexia showed modest activation in the right insula and precentral and postcentral gyri. At the third time point, they showed activation in numerous regions including the bilateral middle, superior temporal, and parietal areas, the bilateral inferior frontal areas, the left insula, the cerebellum (IV, V), and subcortical structures such as the putamen, caudate, amygdala, and hippocampus.
At the first time point, the only significant difference in brain activation to phonological processing between typical readers and children with dyslexia was in the left middle and inferior occipital gyri, where children with dyslexia had reduced activation compared to typical readers (Figure 3 and Table 4). However, at the third time point, children with dyslexia had higher brain activation than typical readers in several areas. These differences were found mainly in the bilateral superior temporal gyrus (STG), middle temporal gyrus (MTG), Heschl's gyri (HG), Rolandic operculum and insula, but also in the left supramarginal gyrus (SMG), precentral and postcentral gyri, and subcortically in the right putamen. After 2 years, typically reading children showed a reduction of brain activation in the language regions of the left hemisphere (STG, insula, inferior frontal gyrus: IFG, precentral gyrus: PrCG, superior and inferior parietal lobules: SPL, IPL, hippocampus; Figure 4). In children with dyslexia, brain activation during phonological processing in the right STG and insula increased with time. Diverging developmental trajectories related to literacy acquisition in typical readers and children with dyslexia were confirmed in a time x group interaction present in the bilateral STG, insula, left MTG, and right frontal cortex (see Figure 5 and Table 4).
The pattern of effects seen in typical readers and children with dyslexia was different when the sample was restricted to children with a familial history of dyslexia. This is due to generally reduced brain activation to Rhyme > Voice in typical readers at risk for dyslexia that was restricted to the left fusiform (TP1) and postcentral gyrus (TP3) (for details see Supplementary Table S5). This is why, at the first time point, children with a familial history of dyslexia who later developed reading impairment had higher brain activation in the bilateral middle frontal gyri as well as in the right inferior parietal lobule than children who became typical readers. Again, at the third time point, children with a familial history of dyslexia who developed dyslexia had higher brain activation than children who did not develop any reading impairment. This Group effects were tested with one-way ANOVA, TP effects were tested with paired t-tests, Group and TP interaction was tested with flexible factorial design.
activation was more extensive than at the first time point and included the bilateral STG, MTG, SMG, Heschl's gyri, Rolandic operculum, left precentral and postcentral gyri, bilateral putamen (subcortically), left thalamus, right amygdala, and hippocampus (see Figure 3 and Supplementary Table S7). Time effects in typical readers with a familial history of dyslexia were restricted to higher activation in the left precentral and postcentral gyri at the third compared to the first time point. Increased activation in the right STG was observed (Figure 5) in children with dyslexia from the familial risk group just as was observed for the whole sample of children with dyslexia. However, in the case of children with a familial history of dyslexia, the interaction between time and group (typical readers vs. dyslexic readers) was not significant.

Behavioral Results
No significant differences in early reading and reading related tests were observed between children with and without a familial history of dyslexia (see Table 1 and Figure 6), although the groups differed in parental scores in the Adult Reading History Questionnaire and in home literacy environment. These two groups also had similar performance in all tests belonging to the dyslexia diagnosis battery (see Table 2) run after 2 years of school education. This pattern of lack of between group results remained unchanged when the sample was restricted to typical readers (see Supplementary Tables S2, S3).

In-Scanner Performance
No significant differences in accuracy in in-scanner tasks were observed between children with and without a familial history of dyslexia (see Table 3). However, children with a familial history of dyslexia were slower in the Rhyme task at the first time point. On the other hand, this difference appeared only at the trend level in typical readers (see Supplementary Table S4).

fMRI Results
Figure 7 depicts brain activation to Rhyme > Voice (see Experimental Design: fMRI tasks) for children with and without a familial history of dyslexia (for details see Supplementary Table S6). When the total sample was split into children with and without a familial history of dyslexia, irrespective of a later dyslexia diagnosis, at the first time point, children with a familial history of dyslexia showed activations only in left subcortical areas. Reduced activation persisted at the third time point and included the bilateral occipital areas and putamen.
Children without a familial history of dyslexia at the first time point engaged an extensive network including the bilateral middle, superior temporal, and parietal areas, the bilateral inferior frontal areas, insula, occipital areas, and the fusiform gyrus as well as subcortical structures such as the inter alia putamen, caudate, amygdala, hippocampus, and thalamus. At TP3, activation was similar to TP1 and present mainly in the middle and superior occipital areas, the inferior frontal areas, and the left fusiform and subcortical structures (for details see Supplementary Table S6).
Comparing children with a familial history of dyslexia to those without showed reduced brain activation in widespread brain areas. These areas included the bilateral STG, MTG, Insula, left IFG and superior frontal gyrus (SFG), SPL, IPL, right lingual gyrus (LG) and cerebellum (III, IV, V), left postcentral and precentral gyri, and left thalamus (subcortically). After 2 years, children with a familial history of dyslexia still presented reduced brain activation but this was restricted to the left MTG, caudate and bilateral anterior cingulate (see Figure 3 and Table 5). In children without a familial history of dyslexia, a reduction of brain activation was seen with time, especially in left hemispheric areas. These included the STG, IFG, SMG, insula and bilaterally in the hippocampi. Children with a familial history of dyslexia showed an increase in activity over time in the left postcentral and precentral gyri (see Figure 4). An interaction between group and time was present in left precentral and postcentral gyri as well as in the left STG and MTG (see Figure 5).
When the sample of children was restricted to typical readers, significant differences were still observed between children with and without a familial history of dyslexia. Again, at both time points, children with a familial history of dyslexia had reduced brain activation compared to children without dyslexia in families (see Supplementary Table S6). At the first time point, group differences were found in the bilateral STG, MTG, insula, left precentral and postcentral gyri, IFG, MFG and SFG, right LG and cerebellum (IV, V), and left thalamus (subcortically). At the third time point, the differences were reduced to the right STG, MTG, left caudate and bilateral putamen (see Supplementary  Table S8 and Figure 3). While children without a familial history of dyslexia showed a reduction of brain activation with time in the left hemispheric language areas (STG, IFG, SMG, HG, Insula, Rolandic Operculum) and bilaterally in the hippocampi, children with a familial history of dyslexia increased activation in the left precentral and postcentral gyri (see Figure 4). The interaction between group and time was significant in the left precentral and postcentral gyri (see Figure 5).

Effects of Dyslexia
In the present study, we investigated longitudinally how neural correlates of phonological processing change during the two first years of reading instruction in typical readers and in children with dyslexia. We also examined these effects in children with and without familial history of dyslexia, irrespectively of dyslexia itself, which will be discussed in the following section.
At the behavioral level, children with dyslexia performed lower than typical readers in reading, phonological awareness, and rapid automatized naming tests at each time point, even at the very beginning of education. The behavioral differences were stable when the sample was restricted to children with a familial history of dyslexia. Similar differences between typical and poor readers have also been reported at a very early stage of literacy development (kindergarden) in other orthographies (Czech and Slovak: Moll et al., 2016, Dutch: Dandache et al., 2014, Finnish: Torppa et al., 2010, German: Schneider et al., 2000, English: Gallagher et al., 2000, for meta-analysis see: Snowling and Melby-Lervåg, 2016). This means that in many orthographies, including Polish, behavioral differences between future poor and fluent readers may be observed much earlier than after several years of reading instruction.
In terms of BOLD activations as registered by fMRI, the pattern of differences observed between typical and dyslexic readers at the beginning of education was different depending on whether the analyses included the whole sample or only children with a familial risk for dyslexia. Reduced activation of the left visual cortex in dyslexics was present only when all children were included in the analysis (consistent with previous fMRI results in visual and orthographic processing tasks, e.g., Dehaene et al., 2010;Boros et al., 2016;Cao et al., 2018). It was no longer significant when the group was restricted to children with a familial history of dyslexia. When the group was restricted to children with a familial history of dyslexia, hyperactivation in the dyslexic group was observed, which included mainly the right fronto-parietal regions. The differences between dyslexia effects in the total sample and in the analyses restricted to children with a familial risk seem to be primarily related to a large reduction in activation of the phonological network in typical readers with a familial risk for dyslexia.
At the third time point, 2 years after the first time point, children with dyslexia showed increased activation compared to typical readers in the bilateral temporal cortices including the auditory cortex, as well as in the left supramarginal and precentral and postcentral gyri, and in the putamen (subcortically). These areas are typically associated with the neural phonological network (for e.g., Brennan et al., 2013) and were also employed by typical readers at the first time point. The observed overactivation might support the notion that children with dyslexia present a delay in the development of their phonological network (Raschle et al., 2011;Morken et al., 2017), as after 2 years of education they activate the regions that typical readers activated at an earlier stage of development. There is, however, a debate on the issue as to whether dyslexia is a developmental delay or a deficit with an altered developmental pathway. With respect to phonological skills, a cross-sectional study in dyslexic children using a developmental trajectory method (Thomas et al., 2009) revealed a delayed trajectory for phonological short term memory and RAN, but showed an atypical trajectory for phonological awareness (Kuppen and Goswami, 2016). An atypical rather than delayed phonological brain network in dyslexia was also found in a cross-sectional fMRI study where dyslexic children exhibited reduced activation during a rhyme judgment task in the bilateral temporo-parietal and frontal cortex relative to both age-matched and readingmatched children (Hoeft et al., 2006). However, cross-sectional studies cannot definitively distinguish between atypical and delayed development of the phonological brain network. More Group effects were tested with one-way ANOVA, TP effects were tested with paired t-tests, Group and TP interaction was tested with flexible factorial design.
longitudinal studies are needed to resolve this debate. As our data covers only the first 2 years of education, the current set of data does not allow us to predict what happens with the phonological network after this period. Nor can we say whether or not the activations observed in the dyslexic group would begin to resemble those of typical readers or whether their behavior over time will be the same. At the very first stage of literacy acquisition, typical readers engaged not only structures typically involved in phonological processing (bilateral superior and middle temporal gyri, left IFG) but also those involved in semantic analysis of words (anterior temporal areas) and in movement planning (premotor and motor areas, caudate, putamen). The reduced activation at the third time point in the phonological network, especially in the left perisylvian areas, suggests that with growing reading experience or time, typical readers automatize phonological processing and therefore the neural circuitry becomes more specialized (Pugh et al., 2013;Dębska et al., 2016;Yu et al., 2018). This finding contrasts with the agerelated increases in activation found in cross-sectional studies (Cone et al., 2008;Brennan et al., 2013). These reductions are more pronounced (not restricted to left IPL) than what was observed in a previous longitudinal study of typical readers (Yu et al., 2018), probably because of a wider time period between the two time points and a larger sample size in the current study.
After 2 years, children with dyslexia (irrespective of familial history status) showed increased brain activation, and they engaged the right hemisphere superior temporal cortex. The right hemisphere is commonly employed during reading by typical beginning readers (Waldie and Mosley, 2000) and its activity declines as literacy develops (Shaywitz et al., 2007). Previous studies reported compensatory shifts to right hemisphere in dyslexia in terms of activation increases (Shaywitz et al., 2007;Simos et al., 2007).
The interaction of group and time observed in the bilateral STG, insula, left MTG, and right frontal cortex supports the hypothesis of a delay in the development of phonological structures in dyslexic readers (Raschle et al., 2011;Morken et al., 2017). These brain areas seem to be involved in early phonological processing related to low reading skills, as they were employed by typical readers only at the beginning of literacy, and by dyslexic children that are 2 years older. Lack of interaction when the group is restricted to children with a familial history of dyslexia suggests some additive influence of both factors. As the interaction appears only if both children with and without a familial history of dyslexia are included, perhaps it is driven by activations of children without a familial risk who reshape the phonological network more clearly.

Effects of Familial History of Dyslexia
In contrast to previous studies, no behavioral differences were found between children with and without a familial history of dyslexia at any time point, even though the groups differed in home literacy environment and parental reading history. Previously, in the relatively transparent Norwegian language, behavioral performance of children with and without a familial history of reading impairment started to differ in the second grade (Specht et al., 2009), though there were no differences at earlier stages. In the more opaque English language, the gap was visible even at earlier stages of education (Raschle et al., 2012). The effect of a familial history of dyslexia may depend on the depth of the orthography. As shallow orthographies are more easily acquired, differences between readers with and without a familial history of dyslexia are less visible, whereas more demanding opaque orthographies lead to earlier observable differences. However, this hypothesis was not supported by a recent meta-analysis (Snowling and Melby-Lervåg, 2016), where no effect of language transparency was revealed. In a metaanalysis based on two to nine studies, depending on the analyzed cognitive skill, the authors found that typical readers with a familial history of dyslexia overcome early delays in vocabulary, grammar, and phonological skills by the time of schooling. However, they still present poorer performance in word and nonword decoding and spelling. This was not the case in our study.
It is plausible that we have overestimated the number of children with a familial history of dyslexia as the group selection was based on quite liberal criteria. Likewise, perhaps fewer children were at severe risk of dyslexia. In the analyzed sample, 27.27% of children with familial history of dyslexia received dyslexia diagnosis, in contrast to the 45% reported in a recent meta-analysis (Snowling and Melby-Lervåg, 2016). On the other hand, the relatively high reading level of children with familial risk of dyslexia may have resulted from the presence of several protective factors (Eklund et al., 2013). These children not only scored moderately to highly in cognitive tests, but also came from well-educated families of relatively high socioeconomic status. What is more, their parents were probably concerned about reading development of their children and perhaps provided some additional training, which we did not control for. Also, we have interpreted the risk of dyslexia as a dichotomous variable, although it can also be understood as a continuum (Snowling et al., 2003). Perhaps a stricter split point would result in betweengroup differences. However, the criterion we used in the current study have previously been applied in research on other languages (Lefly and Pennington, 2000) and have been approved for use in Polish (Bogdanowicz et al., 2015).
Consistent hypoactivation was present for children with a familial history of dyslexia compared to children without such a history, and included the bilateral superior and middle temporal cortex, left precentral and postcentral gyri, left inferior frontal areas, right cerebellum and visual cortex, as well as left thalamus. This pattern is in line with previous studies comparing prereaders with and without a familial history of dyslexia (Raschle et al., 2012) and beginning readers (Dębska et al., 2016). We show for the first time that a familial history of dyslexia modifies brain activation during phonological processing even in typical readers, not only children diagnosed later with dyslexia.
The hypoactivation in the group at risk of dyslexia was visible at the third time point. After 2 years of education, children with a familial history of dyslexia compared to children without such history (even when restricted only to typical readers) still show hypoactivations in the temporal cortex and subcortically in the caudate and putamen, areas which have been implicated in phonological processing (Georgiewa et al., 1999;Cone et al., 2008). Even though cognitive and readingrelated skills are still largely developing at school age, having a familial risk for dyslexia can be considered a more stable risk factor also at the neural level. It may also be interpreted as an indirect proxy for genotype-related properties, in this case resulting in hypoactivation of language structures during phonological processing. Having a familial risk rather than dyslexia itself is also related to atypical planum temporale asymmetry important in speech processing (Vanderauwera et al., 2017) as well as basic auditory processing deficits at the brain level (Hakvoort et al., 2015). Though we are still far from seeing the fuller picture, we provide additional evidence that some cerebral anomalies might be present early in development and might be related to familial risk, irrespective of the later reading outcome.
Typical readers with a familial history of dyslexia show a different neurodevelopmental trajectory, as revealed by group and time interaction in typical readers. With literacy development, typical readers with a familial history of dyslexia show increased engagement of the left precentral and postcentral gyri, activated in covert articulation (Price, 2010), whereas children without a familial history of dyslexia present a decrease in activation of these areas. Perhaps these brain regions support phonological processing in children with a familial history of dyslexia who happen not to develop reading impairment by relying more on covert articulation than on orthographic processing.

Limitations
The main limitation of the project is that, despite a relatively high number of recruited participants, we did not find enough children without familial risk of reading impairment who developed dyslexia to run analyses in a full 2 × 2 model. There were only five children without familial risk of dyslexia identified as having a reading impairment. This number was too small to include this group in analyses. Therefore, instead of including effects of dyslexia and familial history of dyslexia in one model, we decided to run the analyses in a model previously applied in a study of white matter dynamics (Vanderauwera et al., 2017) having a similar group of participants.
Moreover, the age range of our sample was quite wide at each time point. As described in the Section "Materials and Method, " this was a result of ongoing educational reform in Poland. The sample reflected the group of children starting formal education in Poland in the years of our study. However, due to the large variability in the ages of children starting schooling, the oldest participants of the first time point (kindergartens or first graders) were almost as old as the youngest participants of the third time point (second or third graders). Thus, we were unable to include age as a factor in the analyses, and our focus was on educational experience.
Additionally, due to the short history of diagnosing dyslexia in Poland, it was impossible to select children with a familial history of dyslexia only on the basis of formal diagnosis of parents and/or siblings. Although some of the parents presented a full spectrum of dyslexia symptoms, they had not been diagnosed with reading impairment due to the fact that when they attended school, dyslexia was not formally diagnosed in Poland. Instead of relying solely on a formal diagnosis, we were forced to apply a questionnaire measure of familial risk of dyslexia, which was perhaps biased by the memories of the parents. Relying on this questionnaire may have had an impact on the lack of behavioral differences observed between children selected as having a familial risk of dyslexia and children without such risk.
Although we tried to include as many participants of various socio-economic status as possible, the variance in socio-economic status was relatively small and children diagnosed with dyslexia had slightly (although statistically insignificantly) lower socioeconomic status than typical readers. The majority of the participants of the current study lived in Warsaw and had parents of high educational and socio-economic status. As the study was longitudinal and quite demanding for participants, it was impossible to include participants from rural areas of Poland.
Finally, for the results reported in the current paper we used p < 0.005 uncorrected and a cluster size of 50 voxels. Though the results were not corrected for multiple comparisons, in pediatric neuroimaging more liberal thresholds are usually used (usually not corrected), especially in children this young. This is because children show reduced signal-to noise ratios (see e.g., Thomason et al., 2005), time shifted hemodynamic response functions, increased movement throughout the scan, and reduced overall compliance. Additionally, less data is usually acquired to keep scanning time as low as possible. Indeed an uncorrected thresholds of p < 0.005 are commonly reported in fMRI studies on children (Raschle et al., 2012(Raschle et al., , 2014Langer et al., 2015;Saygin et al., 2016). We believe that applying standards common for adult participants' studies may be possible upon solving some technical issues related to scanning (e.g., child appropriate head coils, age appropriate hemodynamic response functions in standard analysis packages etc.).

CONCLUSION
We conclude that the phonological brain network undergoes reorganization during the first 2 years of reading acquisition and that this process proceeds differently depending on the presence of a familial history of dyslexia and reading impairment. Typical readers without risk for dyslexia activate structures responsible for phonological processing already at the beginning of literacy. In this group, reduced brain activation over time during phonological processing is plausibly due to automatization of phonological skills. Children who develop reading impairment present a kind of delay in the development of language and, in particular, phonological structures such as the bilateral STG, left MTG, right insula, and right frontal cortex, where we observed time and group interaction. Finally, typical readers with familial risk of dyslexia also present an atypical development of the neural phonological structures, visible both at the beginning of reading instruction and 2 years later. These children used a presumably efficient neural mechanism of phonological processing, based on the activation of the precentral and postcentral gyri, and achieved a typical level of phonological awareness possibly through the use of silent articulation.

DATA AVAILABILITY STATEMENT
The datasets generated for this study are available on request to the corresponding author.

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by the Warsaw University Ethical Committee. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.

AUTHOR CONTRIBUTIONS
MŁ and KJ wrote the manuscript. AG, AD, and KC contributed to the writing process. MŁ, KJ, KC, and AD analyzed and interpreted the data and prepared the tables and figures. MŁ, AB, KC, AD, and AŻ acquired the data. AM and AG additionally contributed to conception and design of the study. All authors read and revised the article.

FUNDING
This work was supported by grants from the Polish Ministry of Science and Higher Education (IP2011 020271) and the National Science Center (DEC-2011/03/D/HS6/05584 and 2014/14/A/HS6/00294). The project was realized with the aid of CePT research infrastructure purchased with funds from the European Regional Development Fund as part of the Innovative Economy Operational Programme 2007-2013.