Associations of reading skills and properties of cerebral white matter pathways in 8-year-old children born preterm
===================================================================================================================

* Sarah E Dubner
* Michal Ben-Shachar
* Aviv Mezer
* Heidi M Feldman
* Katherine E Travis

## ABSTRACT

**AIM** Children born preterm (PT) experience perinatal white matter injury and later reading deficits at school age. We used two complementary neuroimaging modalities to determine if reading skills would be associated with contemporaneous white matter properties in school-aged PT children.

**METHOD** In 8-year-old PT children (N=29), we measured diffusivity (fractional anisotropy, FA), from diffusion MRI, and myelin content (relaxation rate, R1) from quantitative relaxometry. We assessed reading (Gray’s Oral Reading Test, Fifth Edition) in each child. Whole-brain deterministic tractography coupled with automatic segmentation and quantification were applied to extract FA and R1 along four tracts and assess their statistical association with reading scores.

**RESULTS** Reading-FA correlations were not significant along the four analyzed tracts. Reading-R1 correlations were significantly positive in portions of the left superior longitudinal fasciculus, right uncinate fasciculus, and left inferior longitudinal fasciculus. FA positively correlated with R1 in limited areas of reading-R1 associations, but did not contribute to the variance in reading scores.

**INTERPRETATION** Combining complementary neuroimaging approaches identified relations between reading and white matter properties not found using a single MRI measure. Associations of reading skills and white matter properties may vary across white matter tracts and metrics in PT children.

**What this paper adds**

*   ▪ Preterm children’s reading was associated with white matter myelin content.

*   ▪ Preterm children’s reading was not associated with white matter diffusivity.

Children born very or extremely preterm (PT), at < 32 weeks’ gestation, are at risk for reading problems (1). Deficits have been attributed to cerebral white matter injury and dysmaturity, part of a complex of encephalopathic changes incurred at birth (2). White matter differences consistent with encephalopathy of prematurity can be identified on neonatal diffusion magnetic resonance imaging scans (dMRI) as decreased directional diffusivity (fractional anisotropy, FA) or increased mean diffusivity (3). Differences in dMRI white matter metrics persist from childhood (4) through adulthood (5). White matter changes in PT children are thought to impact neurocognitive functions, such as reading, by altering structural and functional connectivity (6). However, studies of the associations of reading and white matter metrics in PT children show inconsistent findings (7). Identifying a neurobiological explanation for variability in findings is challenging because FA is influenced by multiple tissue properties (8). Complementary imaging methods have been used to interpret dMRI metric associations with neurocognitive outcomes in children born full term (FT) (9). Here we apply this novel approach to understanding the neurobiological basis of reading skills in PT children.

FT children’s reading skills are associated with white matter microstructural metrics derived from dMRI (10,11). Associations in dorsal pathways, including the arcuate fasciculus (Arc) and superior longitudinal fasciculus (SLF) (7,10–12), are thought to relate to transmission of phonological information for converting visual word forms to spoken words (13). Associations in ventral pathways, including the inferior longitudinal fasciculus (ILF) and the uncinate fasciculus (UF) are thought to relate to transmission of semantic and orthographic information (7,12,14,15). Individual differences in reading among PT children have also been associated with white matter microstructural metrics (16–18). However, different patterns of association of reading and white matter microstructural metrics have been found in PT versus FT children (7). The findings suggest that a different balance of white matter properties may lead to individual differences in reading abilities in PT versus FT children.

Quantitative relaxometry (qR) is an alternative MRI technique that can be used to assess white matter microstructure. qR measures the longitudinal relaxation rate of water protons in a magnetic field (R1, seconds-1). R1 is highly correlated with myelin content in white matter areas (19). Developmental changes identified by R1 are presumed to reflect maturational increases in myelin and related molecular factors (20– 22). Structural differences measured by R1 in PT and FT children are presumed to reflect primarily differences in the amount of myelin (23). qR metrics have been related to reading in typically developing FT children (9,24). R1 has not been used to examine individual variations in neurocognitive abilities in PT children. Further, studies have yet to combine both dMRI and qR to assess the contributions of these metrics in explaining variation in reading skills in PT children.

The present study assessed a sample of PT children who were part of a longitudinal study of the neural bases of reading development (15). When these children turned age 8 years, qR became available at our imaging center. To interrogate tissue properties associated with reading in this PT sample, we collected both dMRI and qR in these children. The aims of this study were: (1) To assess the associations between reading skills and FA of reading-related dorsal and ventral white matter pathways. (2) To assess the associations between reading skills and R1 in the same pathways. We hypothesized that we would not observe associations between reading skills and FA (12). We hypothesized that we would observe positive associations between reading skills and R1, given group differences in R1 in a PT and FT sample (23) and the assumption that these the group differences impact function. We assessed the degree of associations between FA and R1 in regions of significant association of reading and one or both white matter metrics and the relative contributions of each measure to reading scores.

## METHOD

### Participants

8-year-old PT children (*n*=43) were recruited between 2012 to 2015 as part of a longitudinal study (15). The present study included participants with useable dMRI and qR data and complete cognitive and reading measures. 14 participants were excluded because they underwent a different MRI protocol (n=2), did not undergo qR imaging (n=9), or showed excessive motion (n=3). The final sample consisted of 29 participants. We did not include a FT comparison group because we have previously documented differences in reading-FA associations (12).

PT birth was defined as gestational age (GA) *≤* 32 weeks. Below this GA, neonates are at greatest risk for white matter injury (2). English-speaking children were recruited from the Lucile Packard Children’s Hospital Stanford High-Risk Infant Follow-Up Clinic, local parent groups, and the surrounding San Francisco Bay Area. Exclusion criteria included neurological factors unrelated to prematurity that would account for white matter differences or impact learning to read English, including congenital anomalies, active seizure disorder, hydrocephalus, sensorineural hearing loss, visual impairment, intelligence quotient <80. The Stanford University Institutional Review Board approved the experimental protocol (#IRB-22233). A parent or legal guardian provided informed written consent. Children provided assent at age 8. Participants were compensated for participation.

Demographic characteristics included socio-economic status, measured using a modified 4-Factor Hollingshead Index (25). Handedness information was measured by the Edinburgh Handedness Inventory (26). Detailed medical information was available for 26 of 29 participants. A neuroradiologist assessed T1w scans collected as part of the current study for 5 features associated with white matter injury (15).

### Neurocognitive Measures

Reading ability measures were obtained as part of a battery of norm-referenced behavioral tests (12,15,18). We analyzed the Oral Reading Index (ORI) generated from the Gray’s Oral Reading Tests, 5th Edition (27). The ORI, a composite score composed of fluency (rate and accuracy) and comprehension, is reported as a standard score (mean=100, standard deviation=15). General intelligence was assessed using the Wechsler Abbreviated Scale of Intelligence (WASI-II) (28). Study assessments were completed in English.

### MRI Acquisition and Analysis

MRI data were acquired on a 3T Discovery MR750 scanner (General Electric Healthcare, Milwaukee, WI, USA) equipped with a 32-channel head coil (Nova Medical, Wilmington, MA, USA) at the Center for Cognitive and Neurobiological Imaging at Stanford University ([www.cni.stanford.edu](http://www.cni.stanford.edu)). Participants were scanned without sedation. High-resolution T1-weighted (T1w) anatomical images were collected using an inversion recovery (IR)-prep 3D fast-spoiled gradient-echo (FSPGR) sequence collected in the sagittal plane (0.9 x 0.9 x 0. 9 mm3 voxel size). The T1w image was used as an anatomical reference for the alignment of the diffusion tensor image (DTI) maps.

dMRI data were acquired with a diffusion-weighted, dual-spin echo, echo-planar imaging sequence with full brain coverage. Diffusion weighted gradients were applied at 30 non-collinear directions with a b-value=1,000 s/mm2. Three volumes were acquired at *b*=0 at the beginning of the scan. We collected 70 axial slices in each participant (TR=8300 ms; TE=83.1 ms; FOV=220 mm; matrix size of 256 x 256, 0.8594 x 0.8594 x 2 mm3 voxel size).

qR relaxometry was acquired with SPGR echo images acquired at four flip angles (*α* = 4°, 10°, 20°, and 30°; TR=14 ms; TE=2 ms) (29). SPGR scans were acquired sequentially, without breaks to ensure that scan parameters were similar across flip angles. A 0.9375 x 0.9375 x 1.5 mm3 voxel size was chosen to acquire the four SPGR scans in a reasonably short time frame, appropriate for a pediatric clinical sample. To obtain an unbiased R1 map, we acquired a spin-echo inversion-recovery (SEIR) sequence with echo-planar imaging (EPI) read-out (30) measured at lower resolution (1.875 x 1.875 x 4 mm3 voxel size). Four inversion times were measured (2400, 1200, 400 and 50 ms; TR = 3 s; TE = 47 ms). Unbiased R1 maps obtained from SEIR-EPI scans were used to correct for RF transmit bias in the higher-resolution SPGR scans as described and validated by (31).

Image preprocessing, registration, deterministic tractography, fiber tract identification, segmentation and quantification are described in the **Supplemental Information**. We computed fractional anisotropy (FA) from dMRI scans and derived R1 from qR scans. The validity of the current methods have been demonstrated in developmental studies (20). Additionally, evidence shows reduced R1 in patients with clinical conditions in which myelin content is decreased (23,31). We characterized white matter metrics from two dorsal tracts (Arc-L and bilateral SLF) and two bilateral ventral white matter tracts (ILF and UF). We selected these pathways *a priori* based on evidence implicating their involvement in phonological and semantic processes (13,32).

**Figure 1** shows the left hemisphere tracts in a representative participant. We generated tract profiles that measured FA or R1 at 30 equidistant locations along the central portion of each fiber tract, bounded by the same ROIs used for tract segmentation. We averaged FA or R1 values across all 30 locations to generate a mean metric for each tract.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/14/2020.12.11.20247965/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2020/12/14/2020.12.11.20247965/F1)

Figure 1. Tractography of selected white matter tracts.
Left hemisphere tract renderings are displayed on a mid-sagittal T1 image from a representative participant. Right hemisphere tract renderings not shown. A) The arcuate fasciculus is shown in blue, B) the superior longitudinal fasciculus is shown in green, C) the uncinate fasciculus is shown in yellow, D) the inferior longitudinal fasciculus is shown in orange. Dashed lines represent the location of the regions of interest (ROIs) used to segment each pathway from the whole-brain tractogram.

### Statistical Approach

Statistical analyses were conducted using SPSS (version 26.0, IBM Corp., 2019). Statistical significance was set at *p*<0.05. Normality of clinical and neurobiological data was assessed using the Shapiro–Wilk test. We used parametric tests for all associations. A False Discovery Rate (FDR) of 5% was calculated for zero-order associations to account for multiple comparisons (33).

### Associations between reading and white matter microstructure metrics

We computed Pearson correlations to assess associations between ORI and mean FA or mean R1. To identify relations obscured by mean measures, we computed correlations between ORI and FA or R1 for 30 equidistant nodes along each tract. The tract profile analysis utilized a nonparametric permutation-based method to control for 30 along-tract comparisons (34), producing a family-wise error cluster size and a critical r-value for each tract. Along-tract correlations were considered significant if (1) correlations at an uncorrected level of p<0.05 occurred in adjacent nodes *≥* the critical cluster size (range 6-9 nodes) or (2) correlations *≥* the critical r-value occurred within a cluster size *≥* 3 nodes. For FA and R1, we computed the mean r-value from the longest tract segment meeting these criteria and applied a false discovery rate (FDR) at p<0.05, to control for multiple comparisons across 7 tracts.

In order to understand whether myelin was the dominant tissue property represented by FA in regions of reading-white matter metric correlations, we computed Pearson correlations of mean FA and mean R1 for the segments with significant reading-white matter correlations. We also conducted a series of multiple linear regression models to assess the unique contributions of mean R1 and mean FA to ORI scores in the segments along each tract with significant reading-white matter correlations.

## RESULTS

### Participant Characteristics

Participant characteristics are shown in **Table 1**. Four of 29 participants had mildly abnormal T1w scans. Mean reading and IQ scores fell within the normal range.

View this table:
[Table 1.](http://medrxiv.org/content/early/2020/12/14/2020.12.11.20247965/T1)

Table 1. Participant Characteristics

### Reading and FA

Mean FA values are shown in **Supplemental Table S1**. No significant correlations were identified between reading and mean FA in any tract (**Table 2)**. Reading and along-tract FA correlations confirmed that mean measures did not obscure significant associations in any tract.

View this table:
[Table 2.](http://medrxiv.org/content/early/2020/12/14/2020.12.11.20247965/T2)

Table 2. Correlations among mean-tract FA, mean-tract R1, and Reading

### Reading and R1

Mean R1 values are shown in **Supplemental Table S1**. Moderate positive correlations were found between reading and mean R1 in the SLF-L, UF-L, UF-R, and ILF-L (**Table 2)**, but those did not reach statistical significance. Along-tract analysis revealed significant positive correlations between reading and R1 in the SLF-L (nodes 1-5; rmean=0.42), UF-R (nodes 13-19; rmean=0.45), and ILF-L (nodes 20-23 rmean=0.53) (p<0.05, FWE corrected for 30 nodes within each tract), as shown in **Figure 2**. These associations remained significant (p<0.05) after controlling for the FDR across 7 tracts. *FA and R1*

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/14/2020.12.11.20247965/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2020/12/14/2020.12.11.20247965/F2)

Figure 2. Associations between along-tract R1 and ORI at age 8 years in PT children.
Each panel includes a rendering of a specific tract with the strength of the Pearson correlation (r) between R1 and ORI at 30 equidistant nodes along the tract represented as a heat-map cylinder surrounding the tract. Each panel also includes a scatter plot representing the strength of the association between mean R1 of the adjacent nodes in the longest region of the tract with significant correlations between R1 and ORI. Panels show the results for the SLF-L, UF-R, and ILF-L. R1 is significantly positively associated with reading within regions of these tracts. Brackets (**[**) indicate approximate regions of the tract with significant positive correlations between FA and R1.

R1 = Relaxation rate (s-1); ORI = Oral Reading Index; SLF-L = Left Superior Longitudinal Fasciculus; UF-R = Right Uncinate Fasciculus; ILF-L = Inferior Longitudinal Fasciculus.

We calculated the strength of association between mean FA and mean R1 of the nodes in the segments of the SLF-L, UF-R, and ILF-L where reading was significantly associated with R1 (**Supplemental Figure S1)**. In the SLF-L, the FA-R1 associations were significant in the region of R1-reading association (r = 0.52). Correlations between mean FA and mean R1 in the UF-R and ILF-L did not reach statistical significance. In separate multiple linear regression models we confirmed that R1 made a significant contribution to ORI all three tracts - SLF-L R2 15% (p = 0.022), UF-R R2 17% (p = 0.017), ILF-L R2 25% (p = 0.003) (**Supplementary Table 2**). FA did not contribute unique variance in any of the models.

## DISCUSSION

In this investigation of 8-year-old PT children, we interrogated the associations of reading skills and concurrent white matter microstructure metrics from two complementary imaging methods, dMRI and qR. We did not find any associations between reading and FA, consistent with previous findings. A novel finding was that we found significant associations between PT children’s reading abilities and R1, a measure from qR imaging that is highly sensitive to myelin content within white matter tracts. The location of the positive associations included a dorsal tract (SLF-L) and two ventral tracts (UF-R, and ILF-L) that have been implicated in reading in many studies of FT children. We found that significant positive associations between FA and R1 localized to the region of reading-R1 associations within the SLF-L, but not within the UF-R or ILF-L. However, in the SLF-L, FA did not make a unique contribution to the variance in reading scores. These findings contribute substantially to understanding the reading outcomes of preterm birth at school age.

### Correlations between reading and white matter metrics

We did not observe associations between reading skills and FA within any of the tracts. This finding is consistent with our previous work that found FA at age 6 was related to reading skills at both age 6 and 8 years in FT children but not in PT children (12,15). The findings suggested that individual differences in reading among PT children are not associated with directional diffusivity at this age. However, positive associations between reading skills and FA have been found at older ages in PT children (7).

In this sample of PT children, R1 proved to be more sensitive to individual variation in reading than did FA. R1 within white matter tracts is considered a proxy for myelin content of the tract. Long-term neurodevelopmental outcomes in PT children are currently thought to result from injury and/or aberrant regeneration and repair of oligodendrocyte precursor cell injury that are part of the encephalopathic changes of PT birth. R1 may be more directly associated with these early processes than is FA. Generally, myelin is thought to contribute to directional diffusion by restricting water diffusion in directions perpendicular to the main direction of the axon. However, other factors also contribute to FA, including axonal diameter, fiber coherence, and proportion of crossing fibers. It is possible that these other factors have a greater impact on FA or neutralize the influence of myelin content on FA in PT children, thereby reducing the strength of reading-FA associations. Alternatively, part of the R1 effect could be related to iron or lipid content that effect R1 and not FA (21,22).

### Tracts in which associations were found

The location of reading-R1 associations is consistent with studies of white matter microstructural metrics derived from diffusion studies in FT children. Several reading skills have been associated with SLF white matter microstructural metrics in samples of FT children (7,12). Reading comprehension has been associated with UF white matter microstructural metrics (14). In 6-year-old FT children, reading measures have been associated with FA of both the UF-R and SLF-L, but associations were not found in these tracts in the present PT cohort at age 6 (15,18). ILF maturation relates to children’s reading abilities (12) and ILF-L microstructural metrics have been associated with reading comprehension in FT children with weaker decoding skills (14). Here, associations between reading and R1 of these tracts in PT children, suggest that these tracts are also important for reading in PT children. However, the white matter microstructural characteristics related to reading in PT children may be different from those in FT children.

### FA-R1 associations

Previous studies of FT children or adults have shown FA to be predominately affected by fiber crossings (35) with minimal correlations between FA and qR metrics (36). In PT children the sequelae of prematurity may disrupt the normal progression of myelination and axonal development, potentially increasing variability in the tissue composition of myelin and other axonal properties in PT children (2). Here, FA and R1 were not consistently correlated in areas of reading-R1 correlation, highlighting the complexity of reading-white matter associations. Likely, not all white matter tracts operate in the same manner in the service of reading. Additionally, white matter properties may not be equally represented by FA and R1 across tracts. Other properties, such as axon diameter, may contribute to FA. Longitudinal studies are needed to clarify how white matter maturation and prematurity affect the balance of tissue properties underlying variations in dMRI metrics throughout childhood. After accounting for the contribution of myelin to variance in reading, FA did not explain significant additional variance in reading in this sample. Future studies combining analytic approaches for measuring fiber dispersion and packing from dMRI (37) with structural MRI techniques for assessing axonal diameter (38) or g-ratio (22) will aid in understanding how myelin and other axonal properties relate to variability in FA and reading within and across clinical populations.

Limitations include a modest sample size and reading scores within the normal range. Future studies with larger samples with greater diversity in reading skills could confirm the findings described here and explore whether associations exist in other white matter pathways. dMRI with increased number of directions would allow use of probabilistic tractography.

## Conclusions

In summary, by combining complementary structural neuroimaging approaches we identified relations between reading and concurrent white matter properties that were not found using a single measure from dMRI in isolation. Variation in relations between FA and R1 in regions of significant reading-R1 correlation suggested that the white matter properties represented by FA vary by fiber group and that myelin plays an important role in reading abilities in PT children compared to other white matter properties. Further exploration of white matter pathways using multiple imaging techniques will contribute to understanding reading-brain associations in the aftermath of preterm birth.

## Supporting information

Supplemental Information [[supplements/247965_file05.pdf]](pending:yes)

Supplemental Figure S1 [[supplements/247965_file06.pdf]](pending:yes)

Supplemental Table S1 [[supplements/247965_file07.pdf]](pending:yes)

Supplemental Table S1 [[supplements/247965_file08.pdf]](pending:yes)

## Data Availability

Anonymized data will be made available upon request to qualified investigators for purposes of replicating procedures and results from the corresponding author.

## Acknowledgments

The authors received funding support from the National Institute of Child Health and Human Development (Grant # R01HD069162 [Feldman], 5K99HD084749 [Travis]), a Young Investigator Award to Dr. Travis, from the Society of Developmental and Behavioral Pediatrics (2014), and a Tashia and John Morgridge endowed fellowship through the Stanford Maternal Child Health Research Institute Clinical Trainee Grant to Dr. Dubner. Funders were not involved in study design, data collection, data analysis, manuscript preparation or publication decisions.

## ABBREVIATIONS

AFQ
:   Automated Fiber Quantification
Arc-L
:   left arcuate fasciculus
Arc-R
:   right arcuate fasciculus
dMRI
:   diffusion magnetic resonance imaging
DTI
:   Diffusion Tensor Imaging
FSPGR
:   fast-spoiled gradient-echo
FA
:   fractional anisotropy
FT
:   full term
GA
:   gestational age
ILF-L
:   left inferior longitudinal fasciculus
ILF-R
:   right inferior longitudinal fasciculus
MRI
:   magnetic resonance imaging
ORI
:   Oral Reading Index
PT
:   preterm
qR
:   Quantitative Relaxometry
R1
:   relaxation rate (seconds-1)
ROI
:   region of Interest
SES
:   socioeconomic status
SLF-L
:   left superior longitudinal fasciculus
SLF-R
:   right superior longitudinal fasciculus
T1w
:   T1-weighted
UF-L
:   left uncinate fasciculus
UF-R
:   Right uncinate fasciculus
WASI-II
:   Wechsler Abbreviated Scale of Intelligence, 2nd Edition.

*   Received December 11, 2020.
*   Revision received December 11, 2020.
*   Accepted December 14, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.Kovachy VN, Adams JN, Tamaresis JS, Feldman HM. Reading abilities in school-aged preterm children: a review and meta-analysis. Dev Med Child Neurol. 2015 May;57(5):410–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/dmcn.12652&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25516105&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 

2.  2.Back SA. White matter injury in the preterm infant: pathology and mechanisms. Acta Neuropathol. 2017 Sep;134(3):331–49.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00401-017-1718-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28534077&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 

3.  3.Ball G, Boardman JP, Aljabar P, Pandit A, Arichi T, Merchant N, et al. The influence of preterm birth on the developing thalamocortical connectome. Cortex. 2013 Jun;49(6):1711–21.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cortex.2012.07.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22959979&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000321169200021&link_type=ISI) 

4.  4.Thompson DK, Lee KJ, Egan GF, Warfield SK, Doyle LW, Anderson PJ, et al. Regional white matter microstructure in very preterm infants: predictors and 7 year outcomes. Cortex. 2014 Mar;52:60–74.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cortex.2013.11.010&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24405815&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 

5.  5.Eikenes L, Løhaugen GC, Brubakk A-M, Skranes J, Håberg AK. Young adults born preterm with very low birth weight demonstrate widespread white matter alterations on brain DTI. NeuroImage. 2011;54(3):1774–85.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2010.10.037&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20965255&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000286302000002&link_type=ISI) 

6.  6.Lubsen J, Vohr B, Myers E, Hampson M, Lacadie C, Schneider KC, et al. Microstructural and Functional Connectivity in the Developing Preterm Brain. Seminars in Perinatology. 2011 Feb 1;35(1):34–43.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1053/j.semperi.2010.10.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21255705&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 

7.  7.Travis KE, Ben-Shachar M, Myall NJ, Feldman HM. Variations in the neurobiology of reading in children and adolescents born full term and preterm. Neuroimage Clin. 2016;11:555–65.
    
    

8.  8.Basser PJ, Pierpaoli C. Microstructural and physiological features of tissues elucidated by quantitative-diffusion-tensor MRI. Journal of Magnetic Resonance - Series B. 1996;111(3):209–19.
    
    

9.  9.Geeraert BL, Chamberland M, Lebel RM, Lebel C. Multimodal principal component analysis to identify major features of white matter structure and links to reading. PLoS One. 2020;15(8):e0233244.
    
    

10. 10.Wang Y, Mauer MV, Raney T, Peysakhovich B, Becker BLC, Sliva DD, et al. Development of Tract-Specific White Matter Pathways During Early Reading Development in At-Risk Children and Typical Controls. Cereb Cortex. 2017 01;27(4):2469–85.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cercor/bhw095&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27114172&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 

11. 11.Yeatman JD, Dougherty RF, Ben-Shachar M, Wandell BA. Development of white matter and reading skills. Proceedings of the National Academy of Sciences of the United States of America. 2012 Oct 30;109(44):E3045–53.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTA5LzQ0L0UzMDQ1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTIvMTQvMjAyMC4xMi4xMS4yMDI0Nzk2NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

12. 12.Bruckert L, Borchers LR, Dodson CK, Marchman VA, Travis KE, Ben-Shachar M, et al. White Matter Plasticity in Reading-Related Pathways Differs in Children Born Preterm and at Term: A Longitudinal Analysis. Front Hum Neurosci. 2019;13:139.
    
    

13. 13.Saur D, Kreher BW, Schnell S, Kümmerer D, Kellmeyer P, Vry M-S, et al. Ventral and dorsal pathways for language. Proc Natl Acad Sci USA. 2008 Nov 18;105(46):18035–40.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTA1LzQ2LzE4MDM1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTIvMTQvMjAyMC4xMi4xMS4yMDI0Nzk2NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

14. 14.Arrington CN, Kulesz PA, Juranek J, Cirino PT, Fletcher JM. White matter microstructure integrity in relation to reading proficiency☆. Brain Lang. 2017;174:103–11.
    
    

15. 15.Dodson CK, Travis KE, Borchers LR, Marchman VA, Ben-Shachar M, Feldman HM. White matter properties associated with pre-reading skills in 6-year-old children born preterm and at term. Developmental Medicine and Child Neurology. 2018;60(7):695–702.
    
    

16. 16.Andrews JS, Ben-Shachar M, Yeatman JD, Flom LL, Luna B, Feldman HM. Reading performance correlates with white-matter properties in preterm and term children. Developmental medicine and child neurology. 2010 Jun;52(6):e94–100.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1469-8749.2009.03456.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19747208&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 

17. 17.Odegard TN, Farris EA, Ring J, McColl R, Black J. Brain connectivity in non-reading impaired children and children diagnosed with developmental dyslexia. Neuropsychologia. 2009 Jul;47(8–9):1972–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuropsychologia.2009.03.009&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19428430&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000266701600022&link_type=ISI) 

18. 18.Travis KE, Adams JN, Kovachy VN, Ben-Shachar M, Feldman HM. White matter properties differ in 6-year old Readers and Pre-readers. Brain Structure and Function. 2017;222(4).
    
    

19. 19.Stüber C, Morawski M, Schäfer A, Labadie C, Wähnert M, Leuze C, et al. Myelin and iron concentration in the human brain: a quantitative study of MRI contrast. Neuroimage. 2014 Jun;93 Pt 1:95–106.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2014.02.026&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24607447&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 

20. 20.Yeatman JD, Wandell BA, Mezer AA. Lifespan maturation and degeneration of human brain white matter. Nat Commun. 2014 Sep 17;5:4932.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms5932&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25230200&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 

21. 21.Filo S, Shtangel O, Salamon N, Kol A, Weisinger B, Shifman S, et al. Disentangling molecular alterations from water-content changes in the aging human brain using quantitative MRI. Nat Commun. 2019 30;10(1):3403.
    
    

22. 22.Weiskopf N, Mohammadi S, Lutti A, Callaghan MF. Advances in MRI-based computational neuroanatomy: from morphometry to in-vivo histology. Current Opinion in Neurology. 2015 Aug;28(4):313–322.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/WCO.0000000000000222&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26132532&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 

23. 23.Travis KE, Castro MRH, Berman S, Dodson CK, Mezer AA, Ben-Shachar M, et al. More than myelin: Probing white matter differences in prematurity with quantitative T1 and diffusion MRI. NeuroImage: Clinical. 2019;
    
    

24. 24.Bruckert L, Travis KE, Mezer AA, Ben-Shachar M, Feldman HM. Associations of Reading Efficiency with White Matter Properties of the Cerebellar Peduncles in Children. Cerebellum. 2020 Dec;19(6):771–7.
    
    

25. 25.Hollingshead A. Four Factor Index of Social Status. New Haven, CT: Yale University Department of Psychology; 1975.
    
    

26. 26.Oldfield RC. The assessment and analysis of handedness: the Edinburgh inventory. Neuropsychologia. 1971 Mar;9(1):97–113.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0028-3932(71)90067-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=5146491&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1971J199600013&link_type=ISI) 

27. 27.Wiederholt JL, Bryant BR. Gray oral reading tests—fifth edition (GORT-5). Austin, TX: Pro-Ed. 2012;
    
    

28. 28.Wechsler D, Hsiao-pin C. WASI-II: Wechsler abbreviated scale of intelligence-(WASI-II). San Antonio, TX: NCS Pearson. 2011;
    
    

29. 29.Fram EK, Herfkens RJ, Johnson GA, Glover GH, Karis JP, Shimakawa A, et al. Rapid calculation of T1 using variable flip angle gradient refocused imaging. Magnetic Resonance Imaging. 1987;
    
    

30. 30.Stikov N, Perry LM, Mezer A, Rykhlevskaia E, Wandell BA, Pauly JM, et al. Bound pool fractions complement diffusion measures to describe white matter micro and macrostructure. NeuroImage. 2011;
    
    

31. 31.Mezer A, Yeatman JD, Stikov N, Kay KN, Cho NJ, Dougherty RF, et al. Quantifying the local tissue volume and composition in individual brains with magnetic resonance imaging. Nature Medicine. 2013;19(12):1667–72.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nm.3390&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24185694&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 

32. 32.Vandermosten M, Boets B, Wouters J, Ghesquière P. A qualitative and quantitative review of diffusion tensor imaging studies in reading and dyslexia. Neurosci Biobehav Rev. 2012 Jul;36(6):1532–52.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neubiorev.2012.04.002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22516793&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 

33. 33.Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B (Methodological). 1995;57(1):289–300.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/2346101&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:A1995QE4&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995QE45300017&link_type=ISI) 

34. 34.Nichols TE, Holmes AP. Nonparametric permutation tests for functional neuroimaging: A primer with examples. Human Brain Mapping. 2002;15(1):1–25.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hbm.1058&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11747097&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000172887500001&link_type=ISI) 

35. 35.Jeurissen B, Leemans A, Tournier J-D, Jones DK, Sijbers J. Investigating the prevalence of complex fiber configurations in white matter tissue with diffusion magnetic resonance imaging. Hum Brain Mapp. 2013 Nov;34(11):2747–66.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hbm.22099&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22611035&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 

36. 36.De Santis S, Drakesmith M, Bells S, Assaf Y, Jones DK. Why diffusion tensor MRI does well only some of the time: Variance and covariance of white matter tissue microstructure attributes in the living human brain. Neuroimage. 2014 Apr 1;89:35–44.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2013.12.003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24342225&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 

37. 37.Zhang H, Schneider T, Wheeler-Kingshott CA, Alexander DC. NODDI: practical in vivo neurite orientation dispersion and density imaging of the human brain. Neuroimage. 2012 Jul 16;61(4):1000–16.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2012.03.072&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22484410&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000305920600028&link_type=ISI) 

38. 38.Assaf Y, Blumenfeld-Katzir T, Yovel Y, Basser PJ. AxCaliber: a method for measuring axon diameter distribution from diffusion MRI. Magn Reson Med. 2008 Jun;59(6):1347–54.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mrm.21577&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18506799&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F14%2F2020.12.11.20247965.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000256266400017&link_type=ISI)