Stroke disconnectome decodes reading networks

Cognitive functional neuroimaging has been around for over 30 years and has shed light on the brain areas relevant for reading. However, new methodological developments enable mapping the interaction between functional imaging and the underlying white matter networks. In this study, we used such a novel method, called the disconnectome, to decode the reading circuitry in the brain. We used the resulting disconnection patterns to predict a typical lesion that would lead to reading deficits after brain damage. Our results suggest that white matter connections critical for reading include fronto-parietal U-shaped fibres and the vertical occipital fasciculus (VOF). The lesion most predictive of a reading deficit would impinge on the left temporal, occipital, and inferior parietal gyri. This novel framework can systematically be applied to bridge the gap between the neuropathology of language and cognitive neuroscience.


Introduction
"This discovery of fire [is] probably the greatest ever made by man, excepting language […]" (Darwin 1871).
The invention of the written language is a pivotal hallmark in human history, and mastering reading and writing (i.e. literacy) requires many years. Written language is critical to society and spans from passing information between individuals to written rules that govern society. In today's world, literacy abilities are indispensable to learning and navigating our social and professional environments.
Reading and writing require the interplay of several processes, including abstraction. For alphabetic languages (sound-based, e.g. English), when we write, we graphically encode (grapheme) linguistic utterances (phoneme) with the aim of a reader being able to decode and reconstruct the meaning during reading accurately. For this information transaction to succeed, one must understand the language associated with the symbols to comprehend the text.
Seminal work on the functional brain activations associated with proficiency in literacy has shown that reading selectively activates bilateral temporo-occipital areas, including the left visual word form area and right occipital areas (Dehaene et al. 2010;Murphy et al. 2019). However, reading proficiency modulates a more comprehensive left hemisphere language network, suggesting that literacy acquisition is tapping into the extended language circuitry in the brain. This network includes the left temporal cortex, premotor, and supplementary motor area (SMA) (Dehaene et al. 2010). Literacy-related activations also extended to the planum temporale (PT) suggesting a top-down modulation Michel Thiebaut de Schotten and Isabelle Hesling have contributed equally. of orthography from spoken inputs. Although several studies have offered evidence that the left PT is involved in speech production and perception, some studies have also enlightened that the left PT is also recruited during word reading (Buchsbaum et al. 2005;Hesling et al. 2019). These changes occur regardless of when literacy was acquired (adulthood vs childhood) (Dehaene et al. 2010). The inferior frontal gyrus, harbouring 'Broca's area', has been associated with various linguistic processes, including phonological segmentation, syntactic processing, and unification (Burton et al. 2000;Friederici 2002;Hagoort 2005). Common denominators to these processes are the aspects of segmenting and linking different types of linguistic information. Reading single words, in our study, highly automated words (i.e. days of the week) is unlikely to engage complex semantic and syntactic processing. Still, single-word reading requires linking phonemic sequences with motor gestures (Flinker et al. 2015) and verbal working memory, which may also recruit the inferior frontal gyrus (Jobard et al. 2003).
Literacy, therefore, requires a network of language-processing, auditory, and visual regions to interact. These relevant regions are distant from each other and, as such, rely on a network of connections. Classical language and auditory areas in the frontal and temporal lobes are connected by the arcuate fasciculus, the most prominent pathway for language functions (for review, see Forkel et al. 2021). The fronto-temporal branch, also known as the long segment of the arcuate fasciculus, has been shown to be relevant for learning new words via the phonological route (López-Barroso et al. 2013). The long segment is also associated with reading skills in children (Yeatman et al. 2011), while the parietal-temporal branch, also known as the posterior segment, has been associated with reading skills in adults . Another network that connects occipital-temporal-frontal cortices includes the inferior longitudinal fasciculus (ILF, occipito-temporal connection) (Catani et al. 2003), inferior fronto-occipital fasciculus (IFOF, occipital-frontal connection) , and the vertical occipital fasciculus (e.g. VOF, intraoccipital connection) (Forkel et al. 2015;Sihvonen et al. 2021;Vergani et al. 2014;Yeatman et al. 2014). Longitudinal studies in children acquiring literacy revealed the left IFOF (Vanderauwera et al. 2018) and ILF ) as mediators of orthographic (visual) but not phonological (auditory) processes. The VOF has been suggested as the link between the dorsal and ventral visual streams and to contribute to reading through the integration of eye movement with word recognition (anterior branch) and communication with visual processes (posterior branch) (Yeatman et al. 2014). In addition to the intrahemispheric communication between dorsal and ventral networks, information integration between both hemispheres (i.e. interhemispheric interaction) is critical to literacy. Reading relies on visual inputs from both eyes that must be synthesised in the brain. Accordingly, the posterior corpus callosum connects the parietal and occipital lobes to mediate this integration (Carreiras et al. 2009). While all these studies sporadically highlight pieces of brain networks that orchestrate the language, visual, and auditory processes necessary for reading, their relationship with the networks whose damage will contribute critically to reading impairment is not unequivocal.
Clinical neuroanatomical studies in patients who temporarily or chronically lost the ability to read after brain lesions only partially overlap with the networks related to literacy acquisition.
For instance, phonological dyslexia, an acquired dyslexia characterised primarily by difficulties in reading pronounceable nonwords, is usually caused by sizeable left hemisphere perisylvian lesions in the middle cerebral artery territory (Luzzatti et al. 2001). Surface dyslexia, the difficulty with whole word recognition and spelling, has been associated with lesions to the anterior lateral temporal lobe (Patterson and Hodges 1992). While letter-by-letter dyslexia, characterised by prolonged reading and a linear word length effect, presents after lesions to the occipital and inferior temporal regions Fiset et al. 2006). Pure alexia, defined as pure word blindness by Dejerine (1892), has been initially characterised by lesions to the angular gyrus. However, this view has been challenged, and current models suggest an occlusion of the left posterior cerebral artery subserving parietal and occipital-temporal regions to cause the symptom (Cohen et al. 2000;Dehaene and Cohen 2011;Philipose et al. 2007;Starrfelt and Shallice 2014). Hence the relationship between the areas functionally activated when reading, the network required for literacy acquisition, and the critical network for reading disorders is missing.
A potential way to fill this gap is to take advantage of recent developments in functional-structural connectivity coupling that provides methods to decode functional MRI activations based on the pattern of brain disconnection derived from stroke-the disconnectome (Thiebaut de Schotten et al. 2020). This method reveals the white matter circuitry connecting activated areas and typical lesions that should interrupt this circuitry, thereby identifying the critical white matter circuitry for a specific cognitive function-in this case, reading. The disconnectome has never been applied in full to the reading circuitry. Therefore, in the present study, we reanalysed data  in light of the disconnectome to reveal the reading circuitry and typical lesions that should impair its functioning.

Methods
A methodological overview is provided in Fig. 1. Briefly, we used previously published data (grey boxes) for a functional reading task ) and the disconnectome 1 3 maps (Thiebaut de Schotten et al. 2020) to calculate the group-level statistical map for reading and use the disconnectome to decode the critical white matter network for reading.

Functional MRI
The present study included a sample of 144 healthy righthanded volunteers (72 women, 27 ± 6 years, age range 19-53 years) from the BIL&GIN database (www. gin. cnrs. fr/ en/ curre nt-resea rch/ axis2/ bilgin-en/), which is a multimodal database harbouring imaging, psychometric, and genetic data to study the neural correlates of brain lateralisation (Mazoyer et al. 2016). All participants presented without a history of developmental, neurological, or psychiatric diagnoses. They were not on medication and showed no abnormality on their structural brain MRI. All participants were French native speakers (Mazoyer et al. 2016). On average, participants spent 4 years at a university. All participants (n = 144) were right-handed, and the group was balanced by sex (72 women).
We used the functional MRI data (fMRI) previously published in Hesling et al. 2019 to generate t-maps of the significant voxels of functional activations during the word-list reading task (Fig. 2). This word reading task is part of a larger study whose runs alternated this word reading task with a sentence reading task .
The stimuli were lists of over-learnt words, making it possible to decrease the load of lexico-semantic and syntactic processing. The lists included the days of the week, the months of the year, and the seasons. This task consisted of 13 trials lasting 14 s. Response times were recorded during the fMRI experiment.
The fMRI experiment included further tasks not analysed in this study. For each run, the participants were shown a scrambled drawing for 1 s, immediately followed by a central fixation crosshair. Then the participants performed a reading Word-List task and had to click on the response pad after task completion. This was followed by a low-level control task in which participants fixated on the central cross and pressed the response pad when the cross changed to a square (both stimuli covered a 0.8° × 0.8° visual area). This part of the trial lasted at least half the total trial duration. It was aimed at refocusing the participant's attention to a nonverbal stimulus and controlling for manual motor response activation and lasted 6-10 s.
For computational purposes, t-maps of functional activations during the reading task were further characterised by measuring the average level of activations in subcortical

T-maps functional MRI activations
The acquisition parameters and the pre-processing steps have been fully described by Hesling et al. 2019. Briefly, the data were acquired on a 3 T Philips Intera Achieva scanner. After normalisation to the BIL&GIN template in the MNI space, the T2*-weighted echo-planar images (EPI) (TR = 2 s; TE = 35 ms; flip angle = 80°; 31 axial slices; field of view, 240 × 240 mm 2 ; isotropic voxel size, 3.75 × 3.75 × 3.75 mm 3 ) were initially smoothed using a 6 mm full-width at halfmaximum Gaussian filter to each functional volume (194 T2*-weighted volumes). The data were processed using global linear modelling using SPM12 (https:// www. fil. ion. ucl. ac. uk/ spm/). For each participant, the BOLD variations of the Word-List reading task were contrasted against the cross-change detection task.
The average brain activation map of BOLD signal variations (at the voxel level) during the reading task was calculated. This activation map corresponded to the t values extracted from the Student's t test for each voxel across the 144 participants. Finally, the resulting average brain map was thresholded to conserve only the significant and positive t values (Bonferroni corrected: a voxel is significantly activated if its p value is < 0.05/n, with n = 188,561 being the number of grey matter voxels for which there is a signal across all participants).  Fig. 2 Schematic description of the functional MRI paradigm for word-list reading. An event was initiated by presenting a scrambled picture for one second, followed by the presentation of a list of words (days, months, seasons) that participants were instructed to read. Upon task completion, participants responded with a button press. Subsequently, they had to indicate when the central cross changed into a square by clicking the button. Institute space (2 × 2 × 2 mm, MNI; http:// www. bic. mni. mcgill. ca) to allow direct individual comparisons. Notably, a second set of disconnectome components exist and can be used to replicate subsequent analyses. The disconnectome components exist as a string of values corresponding to the average level of disconnection in subcortical areas (Catani and Thiebaut de Schotten 2012) as well as in areas derived from a multimodal atlas of the brain surface (Glasser et al. 2016), and as neuroanatomical maps in the MNI152 space. Two sets of component maps-original and replicationhave been provided for replication convenience of any results derived from the disconnectome (see Data availability).

Lesion identification
1333 ischemic lesions were provided by the authors (PN) and were manually delineated on T1-weighted images. Lesion masks are derived by an unsupervised patch-based lesion-segmentation algorithm followed by manual curation as described in Xu and colleagues. The process yields binary masks non-linearly registered to SPM's tissue probability map (see Xu et al. 2018 for details).

Decoding the functional activations with the disconnectome component maps
Linear regression was applied to identify the statistical contribution of each component map (i.e. independent variables) to real group-level fMRI maps for reading (i.e. the dependent variable) (see Fig. 1) using SPSS 22 (https:// www. ibm. com/ produ cts/ spss-stati stics). The group-level fMRI maps (Fig. 3A) were derived from a previously collected reading task in 144 healthy right-handed participants ). The 46 disconnectome components (Thiebaut de Schotten et al. 2020) explained 50% of the variance of reading activation. The linear component provided a loading (Fig. 3B) for each component that was subsequently used as an independent variable in a permuted linear regression (https:// fsl. fmrib. ox. ac. uk/ fsl/ fslwi ki/ Rando mise). Using FSL Randomise produced a map of the white matter that significantly contributes to the functional group-level reading maps.
Since we have two sets of component maps, we replicated the results mapped onto the white matter. The replicability of this analysis was assessed using a Pearson correlation between the results and indicated excellent replicability for fMRI reading mapping onto the white matter with a score of 83%.

Anatomical labelling
The white matter and grey matter were labelled manually by expert anatomists (SJF/MTdS) according to the Atlas of Human Brain Connections (Catani and Thiebaut de Schotten 2012).

Results
The reading task activated a bilateral network of regions in the frontal, temporal, occipital lobe, and subcortical areas (Fig. 3A). In the frontal lobe, activations are centred on Exner's area (ExA)-located just between the inferior frontal gyrus' language area and the frontal eye field (saccadic eye movements) at the intersection of the middle frontal gyrus with the precentral gyrus. Previous studies have shown that words in the visual field activate a widespread neuronal network congruent with classical language areas (Price et al. 1996). Our results replicate this observation by showing activations in the inferior frontal and superior temporal gyrus (Fig. 3A).
Functional activations are also evident in the medial superior frontal gyri, the supplementary motor area (SMA). The SMA has been shown to mediate motor speech control, integrate processing between the two hemispheres and interact with subcortical structures (Hertrich et al. 2016). As a reminder, participants were asked to covertly read each word at a quiet pace to avoid them from quickly extracting the global structure of the list without actually reading it. The striatum is a subcortical structure and part of the basal ganglia that has been shown to play a role in sequencing. The frontostriatal circuit may, therefore, contribute to linguistic sequencing in a domain-specific manner (Chan et al. 2013), and its dysregulation has been linked to reading disorders (Hancock et al. 2017). It is thus likely that their reading included subvocal articulation. The superior, middle, and inferior occipital cortices have also been activated bilaterally.
Mapping the profiles of disconnections identified reading clusters in the occipital lobe centered around the VOF. Another white matter cluster was located is the short U-shaped fibres connecting the precentral to the middle frontal gyrus (Fig. 3B).
Finally, a typical lesion was identified with a disconnection profile that could explain 33% of the variance in the reading task (Fig. 4). Accordingly, this lesion will likely disconnect the white matter pathways 'feeding' into the activation pattern for reading. The typical lesion Fig. 4 Distribution of reading activation correlations with disconnection pattern maps in 1333 patients. A According to our experimental framework, a "pure lesion" to the reading network would only occur in 0.75% of cases. By pure lesion, we mean lesions whose disconnection profile correlates with the reading pattern of activation with a large effect size (r > 0.5). B An example of a large effect size is the correlation between the pattern of functional reading task activations and the pattern of brain disconnection of a single patient. Each dot corresponds to a brain region. This disconnection profile explains 33% of the variance functional reading task activations. C Display of the lesion pattern of the "pure lesion" to the reading network that causes a reading deficit is in the left hemisphere and affects Exner's Area and the precentral gyrus in the frontal lobe; the angular gyrus in the inferior parietal lobe; the superior and middle temporal gyri, and the superior, middle and inferior occipital cortices (Fig. 4C). According to the distribution of disconnection pattern maps in 1333 patients' correlations with reading activation, such a lesion should be considered as rare-occurring in less than 0.75% of the stroke population-as indicated in Fig. 4. This rare occurrence of "pure lesion" to the reading network may explain the rare occurrence of "pure alexia" patients in stroke clinics (Dejerine 1895;Epelbaum et al. 2008;Gaillard et al. 2006;Rupareliya et al. 2017).

Discussion
This study combined the latest research in functional MRI of reading with the disconnectome framework to revisit the neurobiology of reading. This approach couples functional MRI and structural connectivity to disentangle the circuitry, supporting a reading network in healthy controls. We propose a new, more comprehensive model of the reading network based on the convergence of functional, structural, and lesion data. First, we confirmed the functional areas activated by reading tasks during fMRI tasks. Second, we identified the critical white matter network supporting reading. Third, we revealed the typical brain areas essential for reading deficits if damaged by a lesion that disconnects a reading circuitry, which would involve phylogenetically recent white matter tracts.
The cortical networks for reading have been studied for over 30 years using functional MRI (Carreiras et al. 2009;Hesling et al. 2019;Labache et al. 2019;McCrory et al. 2005;Price 2000Price , 2012Price et al. 1999;Price and Mechelli 2005). Multiple cortical and subcortical areas have been identified as relevant for reading processes. For example, the pre-SMA and SMA-proper were associated with sequencing abstract motor plans and motor output control, respectively (Alario et al. 2006;Tremblay and Gracco 2006). The subcortical structures, for example, the thalamus, may be involved in lexical information manipulation processes (Llano 2013). The striatum is a subcortical structure and part of the basal ganglia associated with phonological processing and implicit sequence learning relevant to early language reading acquisition. Hyperactivation in this network has been linked to reading disorders (Hancock et al. 2017). As identified by our analysis, the network of cortical areas involved in reading matches the reading literature and is in line with the dual route theory of reading (Jobard et al. 2003). The areas of highest activation were centred in the posterior frontal lobe, anterior insula, and the occipital lobes in both hemispheres (Fig. 3A).
Reading familiar words, such as the days of the week and months of the year, would not necessitate sublexical assembly as the process is facilitated by lexical knowledge of the whole word. It does, however, rely on the conversion of orthographic to articulatory codes (i.e. spelling to sound conversion). As such, activations in the occipital cortex (visual), the superior temporal gyrus (auditory), and the frontal network of the IFG, MFG, preC, and SMA (articulation) are very plausible (Ekert et al. 2021;Kaestner et al. 2021). fMRI does highlight all areas involved in a task by providing correlational evidence for structural-functional relationships (Damoiseaux and Greicius 2009). Recent non-invasive and invasive brain stimulation methods allow disentangling the critical cortical regions by drawing causal inferences. Using such an approach, the precentral gyrus, which is a strongly activated hub in our fMRI analysis, is relevant for silent reading (Kaestner et al. 2021). A recent review consolidated the literature and summarised the contributions of the individual areas into a dual stream model where pseudoword reading was propagated along a dorsal parietal-frontal route while reading familiar words was processed along a ventral occipito-temporal-frontal route (Turker and Hartwigsen 2021). In this model, the frontal areas are relevant for sound and sequencing, parietal regions for grapheme-phoneme conversion and word meaning, and the temporal areas for phonological analyses match the fMRI result we presented. A recent study identified a reading subnetwork using network-based lesionsymptom mapping that included a ventral stream with hubs in occipital, temporal, parietal, and frontal areas (Li et al. 2021). The difference in the functional tasks may partially explain the different anatomy identified in this study. Li and colleagues used an overt reading task of semantically loaded content (e.g. famous people, animals). By contrast, in the task from Hesling and colleagues, no semantic content was probed, and as such, the fMRI activations lack the classical cortical semantic hubs and circuits (Binder et al. 2009).
Within the large cortical network for reading, another area that stands out is Exner's Area. First postulated as a neuroanatomical 'writing centre' by Exner (Exner 1881), the posterior part of the middle frontal gyrus (MFG), also known as 'Exner's Area', is one of the regions involved in the writing process that is highlighted in the literature (Anderson et al. 1990;Exner 1881;Lubrano et al. 2004;Planton et al. 2013;Roux et al. 2009). Despite significant advancements in neuroimaging techniques, few studies describe patients with isolated lesions in 'Exner's Area' or look at the underlying white matter connections involved in reading and writing processes. The small number of patients in the original description, in tandem with the absence of pure agraphia symptoms plus the interindividual variability in the lesion anatomy, is partially why Exner's hypothesis is still debated (Roux et al. 2010). There are isolated reports of Exner's area involved in reading and writing (Roux et al. 2009) and word learning in the blind (Mizuochi-Endo et al. 2021). In our study, Exner's area has been identified as a contributing cortical region in the fMRI analysis, its disconnection from the inferior frontal gyrus and the precentral gyrus has been shown as a critical white matter connection in the disconnectome analysis, and the typical lesion that would explain most the variance in reading has also been identified in the triangle between the IFG/MFG and preC gyrus (Figs. 3A,  B, 4C).
Accordingly, the disconnectome analysis revealed connections to Exner's area and Wernicke's VOF as essential white matter circuits for reading. The frontal U-shaped fibres of the dorsal precentral gyrus connect the motor cortex and the middle frontal gyrus near Exner's area (Catani et al. 2017) and will be referred to as 'Exner's U-shaped' fibres. The VOF connects the dorsolateral and ventrolateral visual cortex (Takemura et al. 2017;Vergani et al. 2014;Yeatman et al. 2014), including the Visual Word Form Area (VWFA), a part of the ventral occipito-temporal cortex specialised in processing visual formation from words and reading (Dehaene et al. 2002;Wandell et al. 2012). A ventral temporal-occipital disconnection by lesions or surgical interventions has been linked to pure alexia (Epelbaum et al. 2008;Ng et al. 2021;Starrfelt et al. 2010). The anatomy of the VOF has been verified using tractography (Briggs et al. 2018;Keser et al. 2016;Panesar et al. 2019;Schurr et al. 2019;Yeatman et al. 2014), post-mortem dissections (Güngör et al. 2017;Palejwala et al. 2020;Vergani et al. 2014), and comparative anatomy (Takemura et al. 2017). These descriptions were scrutinised against historical postmortem descriptions from Wernicke (Yeatman et al., 2014), Sachs (Forkel et al. 2015;Vergani et al. 2014) and the Dejerines (Bugain et al. 2021). Anatomically, the vertical occipital fibre system is distinct from the parietal-temporal branch of the arcuate fasciculus (e.g. posterior segment) in its trajectories and cortical terminations. The intralobar vertical occipital fasciculus projects between the ventral and dorsal occipital cortices while the interlobar posterior segment of the arcuate fasciculus connects posterior temporal and inferior parietal regions (i.e. supramarginal and angular gyrus) . A recent review identified an association between the VOF and visual tasks, including object recognition and integrating the dorsal and ventral stream during hand actions (Budisavljevic et al. 2018;Forkel et al. 2021). The VOF terminations have also been shown to predict reading-related responses in the human ventral temporal cortex (Grotheer et al. 2021) and early literacy skills in children (Broce et al. 2019). In the clinical setting, the VOF has recently also been associated with visual scores in assessments of patients with multiple sclerosis (Abdolalizadeh et al. 2022). Hence, the disconnectome identification of U-shaped fibres and the VOF as critical connections supporting literacy is a plausible explanation, although other U-shaped fibres, yet unidentified, may also overlap with this result.
These U-shaped fibres might have evolved later on in the phylogenetic tree. The emergence of literacy necessitated accessing our knowledge of spoken language through a novel modality, one that was never anticipated by evolution: vision. While we as a species have been using representations of visual objects and spoken sounds using visual signs (e.g. hieroglyphs) for a long time, our brains mostly evolved for millions of years in a world where literacy, or any form thereof, did not exist (Hodgson, 2019). It is argued that reading developed some 5000 years ago, a time frame too short for significant evolutionary changes. As a consequence, some say that reading would rely on brain regions that initially evolved to process vision (e.g. VOF) and spoken language (e.g. Exner's U-shaped fibres) (Dehaene et al. 2010) and the connections at their intersection (Yeatman and White 2021). However, the specificity of the VWFA has been challenged Devlin 2004, 2003). Some might argue that the reading brain does not exist specifically but emerged through brain plasticity in the grey and the white matter (Carreiras et al. 2009;Dehaene et al. 2010;Thiebaut de Schotten et al. 2014).
This study identified the critical role of fronto-parietal U-shaped fibres underneath Exner's area. The superficial white matter, which includes short-ranged intergyral U-shaped fibres, is more variable between people and a more recent addition to brain evolution (Croxson et al. 2018). If we consider the emergence of literacy as a late appearance in the evolution of the brain, it would be no surprise that reading relies on a more adaptable and flexible white matter network, i.e. U-shaped fibres. This is because there is a gradient of variability where the more superficial white matter is more variable while the deeper white matter is more conserved across individuals. The U-shaped fibres are the most superficial connections in the human brain. They might consequently adjust to new demands as the most variable areas appears to be associated with higher cognitive functions in humans (Croxson et al. 2018). The VOF is also a short-ranging intracortical connection and, as such, is located superficially. The critical involvement of U-shaped fibres and short-range connections in reading may reflect a recent adaptation to the advancement of human skills.
We have demonstrated in this study that the most likely lesion to disrupt the reading circuitry was located in the inferior parietal lobe, the posterior, superior, and middle temporal cortex, the superior, middle, inferior occipital cortex, and Exner's area of the left hemisphere (Fig. 4). The typical lesion affected the right hemisphere's white matter behind the inferior temporal gyrus. This lesion would impact the language areas identified from the fMRI task, the U-shaped fibres, and the vertical occipital fasciculus, unifying the results determined by functional imaging and disconnectome mapping. This result is in line with recent neuroimaging studies suggesting that learned orthographic knowledge is associated with the left inferior temporal-occipital cortex and fusiform gyrus . As well as, the literature suggests that lesions to that area should lead to reading impairments. For example, a voxel-based lesion-symptom study in 111 stroke patients identified the critical lesion areas for reading (left occipito-temporal cortex, and inferior parietal cortex), writing (left inferior parietal lobe, superior temporal cortex including Heschl's gyrus, and sensorimotor cortex), single-word reading (left temporo-occipital cortex, posterior lateral temporal, and inferior parietal cortex), and sentence reading (left temporo-occipital cortex, lateral temporal cortex, superior portion of the temporal pole, inferior and posterior insula) (Baldo et al. 2018). The typical lesion identified from our cohort of 1333 patients matches these results for single-word reading with the addition of Exner's area and the cerebellum. The cerebellum has been discussed in the clinical dyslexia literature as the cerebellar deficit hypothesis, which suggests that cerebellar anomalies could give rise to developmental dyslexia by impairing articulatory/phonological monitoring (Nicolson et al. 2001). The discussion of the cerebellum's role in language functions has recently gained momentum (Mariën et al. 2014;Mariën and Borgatti 2018;Murdoch 2010). An fMRI study in healthy adults with differing reading proficiency has also identified cerebellar activations associated with motor/articulatory processing (Cullum et al. 2019).
Decoding task-based functional MRI imaging using a disconnectome framework is a new and promising method to link structural connectivity with function and behaviour using brain disconnections. This method has already been successfully applied to various cognitive functions and behaviours (Dulyan et al. 2021;Thiebaut de Schotten et al. 2020) and offers an indirect measure of the structural connectivity supporting functional activations. To directly predict reading impairments and estimate measures of reading impairments, we provide the maps associated with this study on Neurovault for future research.
Nevertheless, caution should be applied to these maps as patient results might vary according to demographics, reading abilities or literacy skills. Further, we employed the reading task based on automatic stimuli and tapped into basic language processes. Accordingly, this task is not loaded linguistically and likely bypasses complex language processes. However, complex linguistic fMRI tasks would likely activate the whole brain rendering our decoding impossible. Finally, the low requirement of our reading task makes it suitable for the clinical settings opening up new opportunities for future clinical investigations.
Of all human inventions, the ingenuity of the written language has been momentous and far-reaching in human evolution. Literacy has enabled a considerable increase in the rate and efficacy of cultural evolution. The cultural evolution of society and individuals relies on knowledge transmission within and across generations, which was greatly facilitated by literacy. Neuroimaging offers the tools to uncover the brain circuitry that enables literacy, and this paper is a novel contribution to deciphering the neurobiology of reading. Overall, we demonstrated that the reading circuitry can be decoded in healthy participants' functional activations using a disconnectome approach and identified Exner's U-shaped fibres and the vertical occipital fasciculus as critical white matter connections for reading.

Conflict of interest The authors declare no conflict of interest.
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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.