# White Matter Microstructure is Associated with Auditory and Tactile Processing in Children with and without Sensory Processing Disorder

Yi-Shin Chang, Mathilde Gratiot, Julia P. Owen, Anne Brandes-Aitken, Shivani S. Desai, Susanna S. Hill, Anne B. Arnett, Julia Harris, Elysa J. Marco, Pratik Mukherjee

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

## To cite this version:

Yi-Shin Chang, Mathilde Gratiot, Julia P. Owen, Anne Brandes-Aitken, Shivani S. Desai, et al. White Matter Microstructure is Associated with Auditory and Tactile Processing in Children with and without Sensory Processing Disorder. Frontiers in Neuroanatomy, 2016, 9, pp.169. 10.3389/fnana.2015.00169. [hal-01274069](https://hal-fr.sorbonne-universite.fr/hal-01274069)

HAL Id: hal-01274069

https://hal.sorbonne-universite.fr/hal-01274069

Submitted on 17 Feb 2016

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

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

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

# White Matter Microstructure is Associated with Auditory and Tactile Processing in Children with and without Sensory Processing Disorder

Yi-Shin Chang, Mathilde Gratiot, Julia P. Owen, Anne Brandes-Aitken, Shivani S. Desai, Susanna S. Hill, Anne B. Arnett, Julia Harris, Elysa J. Marco, Pratik Mukherjee

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

## To cite this version:

Yi-Shin Chang, Mathilde Gratiot, Julia P. Owen, Anne Brandes-Aitken, Shivani S. Desai, et al. White Matter Microstructure is Associated with Auditory and Tactile Processing in Children with and without Sensory Processing Disorder. Frontiers in Neuroanatomy, Frontiers, 2016, 9, pp.169. <10.3389/fnana.2015.00169>. <hal-01274069>

HAL Id: hal-01274069

http://hal.upmc.fr/hal-01274069

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_6.png" /></figure>

CrossMark

# White Matter Microstructure is Associated with Auditory and Tactile Processing in Children with and without Sensory Processing Disorder

Yi-Shin Chang<fnref n="1" />, Mathilde Gratiot<fnref n="2" />, Julia P. Owen<fnref n="1" />, Anne Brandes-Aitken³,
Shivani S. Desai³, Susanna S. Hill³, Anne B. Arnett³, Julia Harris³, Elysa J. Marco³,⁴,⁵* and
Pratik Mukherjee<fnref n="1" />

Sensory processing disorders (SPDs) affect up to 16% of school-aged children, and contribute to cognitive and behavioral deficits impacting affected individuals and their families. While sensory processing differences are now widely recognized in children with autism, children with sensory-based dysfunction who do not meet autism criteria based on social communication deficits remain virtually unstudied. In a previous pilot diffusion tensor imaging (DTI) study, we demonstrated that boys with SPD have altered white matter microstructure primarily affecting the posterior cerebral tracts, which subserve sensory processing and integration. This disrupted microstructural integrity, measured as reduced white matter fractional anisotropy (FA), correlated with parent report measures of atypical sensory behavior. In this present study, we investigate white matter microstructure as it relates to tactile and auditory function in depth with a larger, mixed-gender cohort of children 8–12 years of age. We continue to find robust alterations of posterior white matter microstructure in children with SPD relative to typically developing children (TDC), along with more spatially distributed alterations. We find strong correlations of FA with both parent report and direct measures of tactile and auditory processing across children, with the direct assessment measures of tactile and auditory processing showing a stronger and more continuous mapping to the underlying white matter integrity than the corresponding parent report measures. Based on these findings of microstructure as a neural correlate of sensory processing ability, diffusion MRI merits further investigation as a tool to find biomarkers for diagnosis, prognosis and treatment response in children with SPD. To our knowledge, this work is the first to demonstrate associations of directly measured tactile and non-linguistic auditory function with white matter microstructural integrity – not just in children with SPD, but also in TDC.

**Keywords:** diffusion tensor imaging, sensory processing disorders, auditory processing, tactile processing, white matter

OPEN ACCESS

**Edited by:**
Yun-Qing Li,
The Fourth Military Medical University,
China

**Reviewed by:**
Richard S. Nowakowski,
Florida State University, USA
Margaret Lang Bauman,
Boston University School of Medicine,
USA

*Correspondence:
Elysa J. Marco
elysa.marco@ucsf.edu

Received: 25 September 2015

Accepted: 31 December 2015

Published: 26 January 2016

**Citation:**
Chang Y-S, Gratiot M, Owen JP, Brandes-Aitken A, Desai SS, Hill SS, Arnett AB, Harris J, Marco EJ and Mukherjee P (2016) White Matter Microstructure is Associated with Auditory and Tactile Processing in Children with and without Sensory Processing Disorder. *Front. Neuroanat.* 9:169.
doi: 10.3389/fnana.2015.00169

<footnote n="1">¹ Department of Radiology and Biomedical Imaging, University of California, San Francisco, San Francisco, CA, USA,</footnote>

<footnote n="2">² Department of Cognitive Neuroscience, Université Pierre et Marie Curie, Paris, France, ³ Department of Neurology, School</footnote>

<footnote>of Medicine, University of California, San Francisco, San Francisco, CA, USA, ⁴ Department of Psychiatry, University of</footnote>

<footnote>California, San Francisco, San Francisco, CA, USA, ⁵ Department of Pediatrics, University of California, San Francisco, San</footnote>

<footnote>Francisco, CA, USA</footnote>
# INTRODUCTION

Hypo-and/or hyper responsiveness to sensory stimulation is estimated to occur in 5-16% of children within the general population, and 40-80% of children with neurodevelopmental disorders (Ahn et al., 2004). Such sensory dysfunction can hinder a child's ability to accomplish practical, daily activities and age-appropriate learning tasks, thus resulting in long-term impairment of intellectual and social abilities. Sensory processing disorder (SPD) is reported to be highly co-incident with attention deficit hyperactivity disorder (ADHD) and autism spectrum disorders (ASDs). However, it is also clear that children can have sensory processing dysfunction without the degree of attention, language, or social challenges that would meet criteria for ADHD or ASD. This has been referred to in the literature as isolated SPD. The Diagnostic Classification of Mental Health and Developmental Disorders in Infancy and Early Childhood includes a diagnostic label for Regulation Disorders of Sensory Processing (Zero to Three, 2005), but the Diagnostic and Statistical Manual 5 (DSM-V) does not include SPDs as a standalone category. They do now, however, include hyper- or hyporeactivity to sensory input or unusual interest in sensory aspects of the environment in their revised ASD criteria. These sensory processing differences are being increasingly investigated in the field of autism research and recognized as a core and critical clinical feature—however, children with SPD who do not also have social communication deficits that meet ASD criteria can provide insight into the neural underpinnings of sensory processing in particular.

The present literature on SPD primarily utilizes parent/caregiver report measures that describe sensory-related behaviors and physiological measures that provide information about arousal and sensory reactivity. Recently, our group has published two studies using diffusion tensor imaging (DTI) to better define the neural correlates of these sensory processing deficits. Our first study took a whole-brain, data-driven approach to demonstrate decreased fractional anisotropy (FA) and increased mean diffusivity (MD) and radial diffusivity (RD), reflecting reduced microstructural integrity, in the posterior white matter tracts of 16 boys with SPD compared to 24 neurotypically developing boys (Owen et al., 2013). In addition, we found that FA in affected brain regions correlated with atypical auditory, multisensory, and attention-related behaviors as reported by parents on the Sensory Profile. In a subsequent report, we established that children with SPD and those with ASD both demonstrate decreased FA in parietal-occipital tracts whereas only children with ASD show differences in temporal tracts subserving social-emotional processing (Chang et al., 2014). We recently investigated whether direct assessment measures of tactile and auditory processing might also inform our understanding and evaluation of children with atypical sensory related behaviors both with and without an ASD diagnosis. We found that children with both ASD and SPD show impairment in tactile processing (right hand graphesthesia) whereas only the ASD group showed significant impairment in a measure of cortical auditory processing (Demopoulos et al., 2015). However, correlations including all children (controls,

SPD, and ASD) showed a significant association between a direct measure of auditory processing impairment and a parent report measure of real-world communication ability.

As our previous imaging findings were limited to a small cohort of affected boys, we seek to investigate these results in a larger mixed-gender cohort sample. We hypothesize that boys and girls with SPD will show impaired white matter microstructural integrity, with a posterior predominance, relative to typically developing children (TDC). We further hypothesize that this microstructural integrity will correlate with parent report as well as with direct measurements of sensory processing, but that the direct measurements will show stronger correlation with the underlying microstructure.

# MATERIALS AND METHODS

## Demographic, Sensory, Cognitive and Behavioral Data

Children ages 8-12 years were enrolled under an institutional review board approved protocol. SPD subjects were recruited from the UCSF Sensory Neurodevelopment and Autism Program (SNAP) and from local online parent board listings. TDC were recruited from online parent group listings as well as referrals from affiliated sensory neurodevelopment and autism research groups. Informed consent was obtained from the parents or legal guardians, with the assent of all participants. Exclusion criteria were brain malformation or injury, movement disorder, bipolar disorder, psychotic disorder, hearing impairment, full-scale IQ (FSIQ) score <70 on the Wechsler Intelligence Scale for Children-Fourth Edition (WISC-IV; Wechsler, 2003), or meeting criteria for ASD. Subjects were also excluded for any anomalies or artifacts on MR imaging or DTI using criteria explained below in Section "DTI Analysis." A total of 40 right-handed children with SPD (32 male, 8 female) and 41 right-handed TDC (28 male, 13 female) met all inclusion and exclusion criteria for the study.

Autism spectrum disorders was initially screened for using the Social Communication Questionnaire (SCQ; Rutter et al., 2003), a parent report screening measure for autism symptoms, with a score of 15 or above being suggestive of ASD. There were four boys in the SPD cohort who exceeded the SCQ screening cut off score of 15, which prompted a standardized structured play session following the Autism Diagnostic Observation Schedule (ADOS; Lord et al., 2000) and clinical review. All four boys scored less than 3 for the social and communication total score with a score greater than 7 meeting concern for ASD and a score greater than 10 meeting concern for a full autism disorder diagnosis based on DSM IV-TR criteria. Clinical interview and review by Dr. Marco, a pediatric cognitive and behavioral child neurologist was also non-consistent with a clinical ASD diagnosis. None of the TDC cohort and none of the SPD girls met screening criteria for ASD based on the SCQ.

All TDC and SPD subjects were assessed with the Sensory Profile (Dunn and Westman, 1997), a parent report questionnaire which measures atypical sensory related behaviors. A definite difference in each sensory domain is defined as a score that
is greater than two standard deviations from the mean: $\le 25$ out of 40 for auditory processing, $\le 64$ out of 90 for tactile processing, $\le 26$ out of 45 for visual processing, and $\le 23$ out of 35 for multisensory integration. There was one child without an auditory processing score, one without tactile, one without multisensory, two without visual scores, and five without sensory profile total scores given incomplete data from the parent survey. A child was included in the SPD cohort if they carried an outside diagnosis of SPD from a community occupational therapist and they had a score in the definite difference range on at least one of the Sensory Profile subscores listed above. Six of the TDC group scored in the probable difference range for one of the subscales, while the remainder scored in the normative range on all subscales. Finally, two SPD kids and one control were born early (gestational ages 36, 35, 33 weeks, respectively).

The Acoustic Index of the Differential Screening Test for Processing (DSTP) was used to assess auditory processing (Richard and Ferre, 2006) for all subjects. This test was administered by the lead study coordinator in accordance with manualized instructions and with training and supervision by a cognitive and behavioral neurologist. The acoustic index is calculated from totaling correct items in: (1) dichotic listening, an index of interhemispheric auditory processing, in which the participant hears and repeats different numbers simultaneously presented through headphones to each ear; (2) temporal patterning, in which the participant reports the order of high and low tones presented in a sequence; and (3) auditory discrimination, in which a participant repeats nonsense syllables presented in background noise. Thirty four out of the 41 controls and 35 out of the 40 SPD subjects received the DSTP. The graphesthesia subtest of The Sensory Integration Praxis Tests (Ayres, 1989) was used to assess tactile proprioception by asking participants to recreate seven designs (neither numbers nor letters) drawn on the dorsum of each hand with closed eyes. Each drawing is then scored for accuracy. This test was administered by the lead study coordinator in accordance with manualized instructions and with training and supervision by a cognitive and behavioral neurologist. 33 out of the 41 controls and 34 out of the 40 SPD subjects received the graphesthesia assessment.

Summary demographic information is included in Table 1. No significant differences were found in any of the demographic or sensory variables between boys and girls, for either TDC or SPD kids.

## Image Acquisition

MR imaging was performed on a 3T Tim Trio scanner (Siemens, Erlangen, Germany) using a 12 channel head coil. Structural MR imaging of the brain was performed with an axial 3D magnetization prepared rapid acquisition gradient-echo T1-weighted sequence (TE = 2.98 ms, TR = 2300 ms, TI = 900 ms, flip angle of 90°) with in-plane resolution of 1 × 1 mm on a 256 × 256 matrix and 160 1.0 mm contiguous partitions. Whole-brain diffusion imaging was performed with a multislice 2D single-shot twice-refocused spin echo echo-planar sequence with 64 diffusion-encoding directions, diffusion-weighting strength of $b = 2000 \text{ s/mm}^2$, iPAT reduction factor of 2, TE/TR = 109/8000 ms, NEX = 1, interleaved 2.2 mm-thick axial

slices with no gap, and in-plane resolution of 2.2 mm × 2.2 mm on a 100 × 100 matrix. An additional image volume was acquired with no diffusion weighting ($b = 0 \text{ s/mm}^2$). The total diffusion acquisition time was 8.7 min. Structural MRI for all children was reviewed by Dr. Pratik Mukherjee, a pediatric neuroradiologist, who was blind to cohort assignment. No structural anomalies or other clinically significant findings were reported.

## DTI Analysis

### Pre-processing

The diffusion-weighted images were corrected for motion and eddy currents using FMRIB's Linear Image Registration Tool (FLIRT<fnref n="1" />) with 12-parameter linear image registration (Jenkinson et al., 2002). All diffusion-weighted volumes were registered to the reference $b = 0 \text{ s/mm}^2$ volume. To evaluate subject movement, we calculated a scalar parameter quantifying the transformation of each diffusion volume to the reference. Sixteen children (6 TDC, 10 SPD) were excluded for DTI artifacts and/or median relative displacement between volumes greater than 2 mm, where a volume represents a single diffusion directional measurement of the entire brain. This left a total of 81 children (40 SPD, 41 TDC) with DTI datasets meeting quality control criteria. A heteroscedastic two-sample Student's t-test verified that there were no significant differences between these SPD and TDC groups in movement during the DTI scan ($p > 0.05$). The non-brain tissue was removed using the Brain Extraction Tool (BET<fnref n="2" />). FA, mean diffusivity (MD), and radial diffusivity (RD) were calculated using FSL's DTIFIT at every voxel, yielding FA, MD, and RD maps for each subject.

## Group Differences

Tract-Based Spatial Statistics (TBSS) in FSL (Smith et al., 2006) was used to skeletonize and register the diffusion maps for each subject in order to perform voxel-wise comparisons along the white matter skeleton. First, each subject's FA map was non-linearly registered to each other subject's FA map to identify the most representative FA map as a registration target. The registered maps were then averaged and skeletonized to the center of the white matter. Next, each subject's FA data was projected onto this mean skeleton to obtain skeletonized FA maps per subject. MD and RD maps were then registered and projected onto the white matter skeleton. Finally, voxelwise statistics were performed on the skeletonized maps to assess for group differences with non-parametric permutation testing using the *randomize* function from FSL. Based on our prior work (Owen et al., 2013), we tested the contrasts of: TDC > SPD for FA, SPD > TDC for MD, and SPD > TDC for RD. The resulting group difference maps for each comparison were corrected for multiple comparisons over the 3D image volume with threshold-free cluster enhancement (TFCE) (Smith and Nichols, 2009) using a significance threshold of $p < 0.05$. It is important to note that, in TFCE, the cluster, and not the individual voxel, is the ultimate object of statistical inference and therefore every voxel in the cluster has exactly the same level of statistical significance

<footnote n="1">¹www.fmrib.ox.ac.uk/fsl/flirt</footnote>

<footnote n="2">²http://www.fmrib.ox.ac.uk/analysis/research/bet</footnote>
<table><caption>TABLE 1 | Demographic information and sensory scores.</caption><thead><tr><th></th><th>#TDC/#SPD</th><th>TDC (mean ± standard deviation)</th><th>SPD (mean ± standard deviation)</th><th>p</th></tr></thead><tbody><tr><td>Age (years)</td><td>41/40</td><td>10.1 ± 1.1</td><td>9.6 ± 1.2</td><td>0.066</td></tr><tr><td>FSIQ</td><td>41/40</td><td>116 ± 10</td><td>112 ± 13</td><td>0.077</td></tr><tr><td>SP – Auditory</td><td>41/39</td><td>33.8 ± 3.5</td><td>23.2 ± 5.00</td><td><strong>1.2E-16</strong></td></tr><tr><td>SP – Tactile</td><td>41/39</td><td>83.5 ± 5.8</td><td>62.4 ± 12.2</td><td><strong>1.5E-13</strong></td></tr><tr><td>DSTP</td><td>34/35</td><td>36.1 ± 3.53</td><td>32.2 ± 5.6</td><td><strong>0.0017</strong></td></tr><tr><td>Graphesthesia</td><td>33/34</td><td>21.5 ± 3.91</td><td>19.1 ± 4.8</td><td><strong>0.0025</strong></td></tr></tbody></table>

Bolded p values indicate statistically significant group differences using an unpaired student t-test (p < 0.05).

in the final results. The Johns Hopkins University (JHU) ICBM-DTI-81 White-Matter Labeled Atlas (Wakana et al., 2004) was used to determine the anatomic locations of white matter regions.

## Effects of Age and Gender

For each subject, mean FA values were obtained within the significant voxels of each of the four white matter regions implicated in the group difference analyses – the left and right posterior thalamic radiations (PTR), splenium of the corpus callosum (SCC), and retrolenticular limb of the right internal capsule (RLIC), as defined by the JHU white matter atlas. Then, general linear models (GLMs) of FA as a function of group (TDC or SPD), age, and gender were created for the each of these ROIs:

$$FA_{roi} = \beta_{0} + \beta_{1}*group + \beta_{2}*age + \beta_{3}*gender$$

where group is 1 for TDC and 0 for SPD, and gender is 1 for females and 0 for males. The coefficient values ($\beta_1 - \beta_3$) and their significance levels were assessed for each ROI.

<figure><img src="image_7.png" /><figcaption>FIGURE 1 | Tract-Based Spatial Statistics (TBSS) results of group differences in FA, MD, and RD between TDC and SPD. Blue regions indicate voxels of significant decreases in SPD relative to TDC, while yellow regions indicate voxels of significant increases in SPD relative to TDC.</figcaption></figure>
<table><caption>TABLE 2 | Coefficient estimates and *p*-values for the general linear model of FA in a few significant regions as a function of group, age, and gender.</caption><thead><tr><th rowspan="2"></th><th colspan="2">Group (TDC = 1, SPD = 0)</th><th colspan="2">Age</th><th colspan="2">Gender (F = 1, M = 0)</th></tr><tr><th>b1</th><th>p</th><th>b2</th><th>p</th><th>b3</th><th>p</th></tr></thead><tbody><tr><td>PTR-L</td><td>0.028</td><td><strong>0.00034</strong></td><td>0.0066</td><td><strong>0.038</strong></td><td>0.012</td><td>0.14</td></tr><tr><td>PTR-R</td><td>0.028</td><td><strong>0.00053</strong></td><td>0.0080</td><td><strong>0.017</strong></td><td>0.016</td><td>0.063</td></tr><tr><td>RLIC-R</td><td>0.014</td><td><strong>0.011</strong></td><td>0.0033</td><td>0.16</td><td>0.010</td><td>0.12</td></tr><tr><td>SCC</td><td>0.019</td><td><strong>0.00033</strong></td><td>0.0049</td><td><strong>0.022</strong></td><td>0.0039</td><td>0.49</td></tr></tbody></table>

PTR, posterior thalamic radiation; RLIC, retrolenticular limb of the internal capsule; SCC, splenium of the corpus callosum. Bolded *p* values indicate statistically significant coefficient estimates (*p* < 0.05).

<figure><img src="image_8.png" /><figcaption>FIGURE 2 | Tract-Based Spatial Statistics results of correlations of the Sensory Profile tactile score and Graphesthesia with FA, including regression of motion.</figcaption></figure>

## Sensory Correlations with DTI

In order to limit the number of statistical comparisons, FA alone was examined for correlations with our sensory variables. FA was tested for correlations with the Sensory Profile auditory score and DSTP separately to assess the relationship between white matter microstructure and auditory processing. FA was tested for correlations with the Sensory Profile tactile score and Graphesthesia separately to assess the relationship between white matter microstructure and tactile processing. The correlation analyses were performed on a voxel-wise basis along the white matter skeleton. *Randomize* was used to assess for significant positive correlations of FA with each parent report or direct assessment, with regression of the motion parameters, and the resultant statistical maps were corrected for multiple comparisons using TFCE with *p* < 0.05.

As a *post hoc* analysis to determine contributions of group, age, and gender to the correlational results, for each sensory metric, FA was averaged across the voxels of significant clusters separately for each subject, and separately for several JHU white matter regions. Then, GLMs were constructed with the cognitive metric

as the response variable, and FA, group (TDC vs. SPD), age, and gender as predictor variables. The significance of each predictor variable for each model was determined.

To further investigate the contributions of different types of auditory processing to the correlational results found with DSTP, an additional *post hoc* analysis of correlations between the three subscores of the DSTP acoustic subtest and FA was conducted using the same approach described above.

## RESULTS

### Group Differences

Significantly lower FA, and higher MD and RD values, were found in the SPD cohort relative to TDC (Figure 1). Visual reference with the JHU white matter atlas reveals that the differences in FA are primarily localized to the bilateral PTRs, the SCC, and the right RLIC. MD and RD showed extensive elevations throughout the white matter. Consistent with the FA results, MD and RD are elevated only in the posterior corpus callosum. Also consistent
<table><caption>TABLE 3 | Number of significantly correlated voxels in several ROIs, along with results of the GLMs of Sensory Profile tactile score and Graphesthesia as functions of group, FA, age, and gender.</caption><thead><tr><th rowspan="2"></th><th colspan="3">SP Tactile</th><th colspan="3">Graphesthesia</th></tr><tr><th># sig vox</th><th>p_FA</th><th>p_TDCvSPD</th><th># sig vox</th><th>p_FA</th><th>p_TDCvSPD</th></tr></thead><tbody><tr><td>ACR-L</td><td>659</td><td>0.074</td><td><strong>4.1E-13</strong></td><td>361</td><td>0.11</td><td><strong>0.024</strong></td></tr><tr><td>ACR-R</td><td>372</td><td>0.14</td><td><strong>2.0E-12</strong></td><td>938</td><td><strong>0.039</strong></td><td><strong>0.020</strong></td></tr><tr><td>SCR-L</td><td>179</td><td><strong>0.0037</strong></td><td><strong>5.6E-15</strong></td><td>181</td><td><strong>0.034</strong></td><td><strong>0.012</strong></td></tr><tr><td>SCR-R</td><td>24</td><td>0.85</td><td><strong>9.0E-14</strong></td><td>170</td><td><strong>0.0072</strong></td><td><strong>0.0080</strong></td></tr><tr><td>PCR-L</td><td>221</td><td><strong>0.031</strong></td><td><strong>8.4E-14</strong></td><td>232</td><td><strong>0.023</strong></td><td><strong>0.012</strong></td></tr><tr><td>PCR-R</td><td>229</td><td>0.11</td><td><strong>1.8E-13</strong></td><td>398</td><td>0.056</td><td><strong>0.031</strong></td></tr><tr><td>ALIC-L</td><td>267</td><td>0.29</td><td><strong>9.0E-13</strong></td><td>-</td><td>-</td><td>-</td></tr><tr><td>ALIC-R</td><td>308</td><td>0.10</td><td><strong>2.3E-12</strong></td><td>234</td><td><strong>0.014</strong></td><td><strong>0.034</strong></td></tr><tr><td>PLIC-L</td><td>590</td><td><strong>0.014</strong></td><td><strong>5.9E-14</strong></td><td>49</td><td>0.20</td><td><strong>0.0077</strong></td></tr><tr><td>PLIC-R</td><td>506</td><td><strong>0.012</strong></td><td><strong>2.1E-13</strong></td><td>322</td><td><strong>0.010</strong></td><td><strong>0.016</strong></td></tr><tr><td>RLIC-L</td><td>393</td><td><strong>0.033</strong></td><td><strong>3.9E-13</strong></td><td>375</td><td>0.721</td><td><strong>0.011</strong></td></tr><tr><td>RLIC-R</td><td>333</td><td>0.091</td><td><strong>3.4E-13</strong></td><td>386</td><td>0.301</td><td><strong>0.021</strong></td></tr><tr><td>PTR-L</td><td>720</td><td><strong>0.0041</strong></td><td><strong>3.4E-12</strong></td><td>910</td><td><strong>0.026</strong></td><td>0.056</td></tr><tr><td>PTR-R</td><td>644</td><td><strong>0.032</strong></td><td><strong>1.5E-12</strong></td><td>425</td><td>0.057</td><td>0.062</td></tr><tr><td>GCC</td><td>19</td><td>0.069</td><td><strong>1.8E-14</strong></td><td>1039</td><td><strong>0.031</strong></td><td><strong>0.022</strong></td></tr><tr><td>BCC</td><td>398</td><td>0.21</td><td><strong>1.6E-13</strong></td><td>1803</td><td><strong>0.046</strong></td><td><strong>0.012</strong></td></tr><tr><td>SCC</td><td>1010</td><td><strong>0.027</strong></td><td><strong>4.8E-12</strong></td><td>1654</td><td><strong>0.0012</strong></td><td>0.055</td></tr><tr><td>CGC-L</td><td>172</td><td>0.052</td><td><strong>2.4E-13</strong></td><td>73</td><td><strong>0.011</strong></td><td><strong>0.038</strong></td></tr><tr><td>CGC-R</td><td>2</td><td>0.41</td><td><strong>9.0E-14</strong></td><td>-</td><td>-</td><td>-</td></tr><tr><td>EC-L</td><td>602</td><td><strong>0.0058</strong></td><td><strong>2.5E-14</strong></td><td>60</td><td><strong>0.021</strong></td><td><strong>0.014</strong></td></tr><tr><td>EC-R</td><td>278</td><td><strong>0.033</strong></td><td><strong>5.3E-14</strong></td><td>370</td><td><strong>0.033</strong></td><td><strong>0.019</strong></td></tr><tr><td>SLF-L</td><td>98</td><td><strong>0.0019</strong></td><td><strong>5.9E-15</strong></td><td>549</td><td><strong>0.0059</strong></td><td><strong>0.012</strong></td></tr><tr><td>SLF-R</td><td>102</td><td>0.44</td><td><strong>4.9E-13</strong></td><td>161</td><td><strong>0.0077</strong></td><td><strong>0.010</strong></td></tr><tr><td>SS-L</td><td>156</td><td><strong>0.043</strong></td><td><strong>1.6E-13</strong></td><td>186</td><td>0.099</td><td><strong>0.013</strong></td></tr><tr><td>SS-R</td><td>187</td><td>0.17</td><td><strong>5.3E-13</strong></td><td>169</td><td><strong>0.021</strong></td><td><strong>0.043</strong></td></tr></tbody></table>

The significance of group and FA to the GLMs are displayed. The significance of age and gender are described in the main text. ACR, anterior corona radiata; SCR, superior corona radiata; PCR, posterior corona radiata; ALIC, anterior limb of the internal capsule; PLIC, posterior limb of the internal capsule; RLIC, retrolenticular limb of the internal capsule; PTR, posterior thalamic radiation; GCC, genu of the corpus callosum; BCC, body of the corpus callosum; SCC, splenium of the corpus callosum; CGC, cingulate gyrus part of the cingulum; EC, external capsule; SLF, superior longitudinal fasciculus; SS, sagittal stratum. Bolded p values indicate statistically significant effects (p < 0.05).

with the FA results, the most significant elevations of both MD and RD occur in the bilateral PTR and right RLIC.

## Effects of Age and Gender

The coefficient estimates and *p*-values for the GLM of FA as a function of group, age, and gender, are displayed in Table 2. The effect of group remained strongly significant when accounting for age and gender. As expected, increasing age contributed to higher FA in the bilateral PTR and SCC. There were no significant gender effects.

## Sensory Correlations with DTI
### Tactile Correlations with DTI

There were widespread, significant positive associations of FA across groups with the Sensory Profile tactile score and with graphesthesia, after regression of motion parameters (Figure 2). The number of significantly correlated voxels in several ROIs is included in Table 3.

The significance of the group and FA predictor variables in the GLMs for prediction of the Sensory Profile tactile score and Graphesthesia are included in Table 3. The effects of group

(TDC vs. SPD) for the Sensory Profile tactile score model were much more strongly significant than the group effects in the Graphesthesia model. The FA effect was significant for many ROIs for both the Sensory Profile tactile score and for Graphesthesia, though the effect was significant for more ROIs with Graphesthesia than with the Sensory Profile tactile score (12 ROIs for Sensory Profile tactile score, 16 ROIs for Graphesthesia). For the Sensory Profile tactile response variable, neither age nor gender was significant for any of the models. For graphesthesia, gender was significant for every model, with females achieving higher scores than males.

Representative plots of the Sensory Profile tactile score model and Graphesthesia vs. FA in significant voxels of the bilateral PTR and SCC are displayed in Figure 3. The regression lines for TDC and SPD (of sensory score vs. FA) are more concordant with one another for Graphesthesia than for the Sensory Profile tactile score.

### Auditory Correlations with DTI

There were widespread, significant positive associations of FA across groups with the Sensory Profile auditory score and with
<figure><img src="image_9.png" /><figcaption>FIGURE 3 | Representative scatterplots and regression lines of the Sensory Profile tactile score and Graphesthesia (mean-centered) versus FA in significant voxels (differing between the two sensory variables) of the posterior thalamic radiations (PTRs) (both sides combined) and splenium of the corpus callosum (SCC).</figcaption></figure>

DSTP, after regression of motion parameters (Figure 4). The number of significantly correlated voxels in several ROIs is included in Table 4.

The significance of the group and FA predictor variables in the GLMs for prediction of the Sensory Profile auditory score and DSTP are included in Table 4. Similarly to the tactile model, the group effects (TDC vs. SPD) in the Sensory Profile auditory score model were much more strongly significant than group effects in the DSTP model, with the DSTP group effects only reaching weak significance in a few ROIs. The FA effect was significant for many ROIs for DSTP, but only in a few ROIs for the Sensory Profile auditory score. For the Sensory Profile auditory response variable, neither age nor gender was significant for any of the models. For DSTP, age was significant for every model, with increasing age giving rise to higher scores.

Representative plots of the Sensory Profile tactile score model and Graphesthesia versus FA in significant voxels of the bilateral

PTR and SCC are displayed in Figure 5. The regression lines for TDC and SPD (of sensory score vs. FA) are more concordant with one another for Graphesthesia than for the Sensory Profile tactile score.

## DSTP Subscore Correlations with DTI

Given the extensive and robust correlations of FA with the DSTP acoustic subscore, additional post-hoc tests for correlations were performed between FA and the three subscores of the DSTP acoustic test – dichotic digits, temporal patterning, and auditory discrimination. The number of significantly correlated voxels in several ROIs is included in Table 5. Overall, the dichotic digits task was significantly associated with the largest number of voxels, as compared with temporal patterning and auditory discrimination.

The significance of the group and FA predictor variables in the GLMs for prediction of each DSTP subscore are included in
<figure><img src="image_10.png" /><figcaption>FIGURE 4 | Tract-Based Spatial Statistics results of correlations of the Sensory Profile auditory score and DSTP with FA, including regression of motion.</figcaption></figure>

**Table 5.** The group effect (TDC vs. SPD) was only significant for the temporal patterning task, while the FA effect was significant for all tasks, but most strongly for dichotic digits. For the dichotic digits test and temporal patterning test, age was a significant contributor for most models, with increasing age giving rise to higher scores, but gender did not contribute significantly to the model. For auditory discrimination, neither age nor gender was significant.

Representative plots of each DSTP subscore versus FA in significant voxels of the bilateral PTR and SCC are displayed in **Figure 6**. The regression lines for TDC and SPD (of sensory score vs. FA) are concordant with one another for dichotic digits and auditory discrimination, but not temporal patterning. A ceiling effect of the dichotic digits and temporal patterning tasks can be observed from these plots, limiting the range of sensory function that can be differentiated between subjects, particularly in the TDC cohort.

DISCUSSION

These findings corroborate and generalize our previous work demonstrating the role of disrupted posterior white matter microstructure in SPD. Furthermore, the larger, mixed-gender cohort unmasks a more extensive distribution of white matter differences which includes anterior white matter. More importantly, to our knowledge, this work is the first to demonstrate a relationship between direct measurements of tactile function and non-linguistic auditory function with white

matter microstructural integrity not just in SPD, but also in TDC.

Group Differences

As in our prior work (Owen et al., 2013), our results show strong decreases of white matter microstructural integrity in the posterior projection and commissural tracts of the bilateral PTR and the SCC, which respectively contain all primary sensory projection pathways excluding olfaction, and connect homologous sensory cortical regions (Huang et al., 2005; Hofer and Frahm, 2006). Our findings additionally reveal marked reductions of microstructural integrity in the RLIC, which contains visual and auditory projection fibers. Overall, microstructural integrity is disrupted extensively, but with a posterior bias, throughout the white matter.

Sensory Correlations with DTI

We find that correlations of white matter microstructure with direct measurements of tactile and auditory processing are stronger than the correlations of microstructure with parent report measures, likely due to the more objective nature of the direct measurements. Furthermore, the stronger concordance between the TDC and SPD regression lines of FA with the direct sensory measurements of Graphesthesia and DSTP, as compared with the sensory profile tactile and auditory scores, suggest that these direct measurements map more closely to the underlying biology. The offset between the TDC and SPD regression lines of the sensory profile metrics with FA may reflect biased parent reporting as a function of the presence or absence of an SPD
<table><caption>TABLE 4 | Number of significantly correlated voxels in several ROIs, along with results of the GLMs of Sensory Profile auditory score and DSTP as functions of group, FA, age, and gender.</caption><thead><tr><th rowspan="2"></th><th colspan="3">SP - Auditory</th><th colspan="3">DSTP</th></tr><tr><th>Num vox</th><th>p_FA</th><th>p_TDCvSPD</th><th>Num vox</th><th>p_FA</th><th>p_TDCvSPD</th></tr></thead><tbody><tr><td>ACR-L</td><td>562</td><td>0.13</td><td><strong>9.9E-15</strong></td><td>1155</td><td><strong>0.014</strong></td><td>0.077</td></tr><tr><td>ACR-R</td><td>573</td><td>0.26</td><td><strong>9.7E-15</strong></td><td>1371</td><td><strong>0.00036</strong></td><td>0.12</td></tr><tr><td>SCR-L</td><td>79</td><td>0.17</td><td><strong>1.7E-15</strong></td><td>148</td><td><strong>0.067</strong></td><td><strong>0.038</strong></td></tr><tr><td>SCR-R</td><td>61</td><td>0.77</td><td><strong>2.5E-15</strong></td><td>464</td><td><strong>0.029</strong></td><td><strong>0.043</strong></td></tr><tr><td>PCR-L</td><td>48</td><td>0.26</td><td><strong>2.5E-15</strong></td><td>236</td><td><strong>0.081</strong></td><td><strong>0.050</strong></td></tr><tr><td>PCR-R</td><td>298</td><td>0.11</td><td><strong>4.8E-15</strong></td><td>518</td><td><strong>0.016</strong></td><td>0.068</td></tr><tr><td>ALIC-L</td><td>276</td><td>0.06</td><td><strong>1.1E-14</strong></td><td>398</td><td><strong>0.0074</strong></td><td>0.18</td></tr><tr><td>ALIC-R</td><td>374</td><td>0.20</td><td><strong>1.8E-14</strong></td><td>576</td><td><strong>0.010</strong></td><td>0.19</td></tr><tr><td>PLIC-L</td><td>339</td><td>0.28</td><td><strong>7.8E-16</strong></td><td>491</td><td><strong>0.016</strong></td><td>0.058</td></tr><tr><td>PLIC-R</td><td>173</td><td>0.64</td><td><strong>2.8E-15</strong></td><td>553</td><td><strong>0.011</strong></td><td>0.082</td></tr><tr><td>RLIC-L</td><td>-</td><td>-</td><td>-</td><td>539</td><td><strong>0.016</strong></td><td>0.075</td></tr><tr><td>RLIC-R</td><td>5</td><td>0.69</td><td><strong>6.4E-16</strong></td><td>475</td><td><strong>0.039</strong></td><td>0.11</td></tr><tr><td>PTR-L</td><td>396</td><td><strong>0.026</strong></td><td><strong>6.8E-14</strong></td><td>647</td><td><strong>0.029</strong></td><td>0.23</td></tr><tr><td>PTR-R</td><td>429</td><td><strong>0.029</strong></td><td><strong>6.5E-14</strong></td><td>768</td><td><strong>0.013</strong></td><td>0.29</td></tr><tr><td>GCC</td><td>551</td><td>0.25</td><td><strong>8.1E-15</strong></td><td>1234</td><td><strong>0.020</strong></td><td>0.086</td></tr><tr><td>BCC</td><td>1718</td><td>0.055</td><td><strong>9.2E-16</strong></td><td>2173</td><td><strong>0.048</strong></td><td><strong>0.046</strong></td></tr><tr><td>SCC</td><td>1293</td><td>0.096</td><td><strong>5.3E-14</strong></td><td>1032</td><td><strong>0.027</strong></td><td>0.18</td></tr><tr><td>CGC-L</td><td>102</td><td><strong>0.038</strong></td><td><strong>1.7E-15</strong></td><td>255</td><td><strong>0.022</strong></td><td>0.067</td></tr><tr><td>CGC-R</td><td>7</td><td><strong>0.021</strong></td><td><strong>3.9E-15</strong></td><td>10</td><td>0.26</td><td><strong>0.041</strong></td></tr><tr><td>EC-L</td><td>399</td><td><strong>0.00038</strong></td><td><strong>2.7E-17</strong></td><td>383</td><td>0.091</td><td><strong>0.048</strong></td></tr><tr><td>EC-R</td><td>289</td><td><strong>0.010</strong></td><td><strong>1.7E-16</strong></td><td>1104</td><td><strong>0.002</strong></td><td>0.063</td></tr><tr><td>SLF-L</td><td>-</td><td>-</td><td>-</td><td>560</td><td><strong>0.018</strong></td><td>0.068</td></tr><tr><td>SLF-R</td><td>141</td><td>0.062</td><td><strong>5.1E-15</strong></td><td>852</td><td><strong>0.00040</strong></td><td>0.060</td></tr><tr><td>SS-L</td><td>-</td><td>-</td><td>-</td><td>246</td><td><strong>0.0055</strong></td><td><strong>0.047</strong></td></tr><tr><td>SS-R</td><td>-</td><td>-</td><td>-</td><td>308</td><td><strong>0.00016</strong></td><td>0.36</td></tr></tbody></table>

The significance of group and FA to the GLMs are displayed. The significance of age and gender are described in the main text. ACR, anterior corona radiata; SCR, superior corona radiata; PCR, posterior corona radiata; ALIC, anterior limb of the internal capsule; PLIC, posterior limb of the internal capsule; RLIC, retrolenticular limb of the internal capsule; PTR, posterior thalamic radiation; GCC, genu of the corpus callosum; BCC, body of the corpus callosum; SCC, splenium of the corpus callosum; CGC, cingulate gyrus part of the cingulum; EC, external capsule; SLF, superior longitudinal fasciculus; SS, sagittal stratum. Bolded p values indicate statistically significant effects (p < 0.05).

diagnosis. More explicitly, if a child without an SPD label and a child with an SPD label exhibit the same level of function for a given sensory processing domain, parents of the child who has not been attributed an SPD label may be less likely to report abnormalities than the parents of the child who has been clinically assigned an SPD label.

The relative lack of anatomic specificity of these correlations of DTI with sensory processing measures may be due, in part, to the high degree of microstructural covariance between different white matter tracts (Wahl et al., 2010), with a spatially global factor accounting for almost half of the variance in FA (Penke et al., 2010). It is also likely reflective of the multiple functional brain networks required for each of these tests, as reviewed below.

## Tactile Processing

Diffusion tensor imaging has previously been used to link degree of periventricular white matter injury in the PTR, as assessed by size reduction of white matter tracts on visual inspection, with contralateral touch threshold, proprioception, and motor severity in children with cerebral palsy (Hoon et al., 2009). However, this prior work used a subjective rater system of white matter

injury. One study has reported a correlation between tactile defensiveness and FA in the inferior longitudinal fasciculus of children with ASD (Pryweller et al., 2014). However, the metric of tactile defensiveness is an assessment of sensory-related behavior, rather than early-stage sensory function. While we are not aware of other studies reporting a relationship between white matter microstructure and primary or secondary tactile processing, there are a multitude of functional imaging studies have demonstrated primarily fronto-parietal network activations in tasks involving either active or passive touch discrimination of macrogeometric object features such as shape or length (Bodegård et al., 2001; Van de Winckel et al., 2005). These circuits include the ventral and dorsal premotor cortex, secondary somatosensory area, superior parietal lobe, anterior part of the intraparietal sulcus (AIP), and supramarginal gyrus (Kawashima et al., 1994; O'Sullivan et al., 1994; Hadjikhani and Roland, 1998; Roland and Zilles, 1998; Binkofski et al., 1999; Bodegård et al., 2000, 2001; Servos et al., 2001; Stoeckel et al., 2003, 2004; Stoesz et al., 2003).

Our results are the first to demonstrate associations of white matter microstructure with tactile processing, both among children with SPD as well as among TDC. Both the Sensory
<figure><img src="image_11.png" /></figure>

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

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

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

FIGURE 5 | Representative scatterplots and regression lines of the Sensory Profile auditory score and DSTP (mean-centered) versus FA in significant voxels (differing between the two sensory variables) of the posterior thalamic radiations (PTR) (both sides combined) and splenium of the corpus callosum (SCC).

Profile tactile score and Graphesthesia are associated with FA in regions subserving primary sensory processing, including the posterior projection tracts of the superior and posterior corona radiata, the posterior limbs of the internal capsule, and the PTR. Both are also associated with the SCC, which connects homologous sensory areas. However, they are also associated with microstructure in associational tracts such as the external capsules, superior longitudinal fasciculus, sagittal stratum. Unlike the Sensory Profile tactile score, Graphesthesia demonstrates associations with the frontal regions of the right anterior corona radiata, anterior limb of the RLIC, and the genu of the corpus callosum. The widespread nature of these correlations may in part reflect the non-specificity of the assessments used, in addition to the previously mentioned contribution of high microstructural covariance of white matter (from the beginning of the Discussion). For example, in addition to primary tactile processing, Graphesthesia engages secondary modalities requiring synthesis and interpretation of the primary

tactile inputs, including the spatio-temporal analysis of these inputs to form a visual representation of what was drawn on the back of the hand, and the motor planning and coordination to re-create this image.

The tactile sense develops the earliest among all sensory systems, with somatosensory responses detected in utero as early as 8 weeks of gestation (Montagu, 1986). Tactile feedback from mechanoreceptors in the skin and joints critically guide the development of motor skills during childhood, both in gross motor (Metcalfe et al., 2005) and fine motor (Soechting and Flanders, 2008) functions. Tactile processing has been further suggested to play a central role in the development of social and communicative behaviors, with tactile-centered therapies shown to effectively modulate arousal, attention, and sensory defensiveness (Cascio, 2010). It is thus critical to understand the neural correlates of tactile processing in order to better understand the downstream effects of its abnormalities as well as to better design and evaluate tactile-centered therapies.
## Auditory Processing

One prior study reported associations between FA and performance on auditory processing tasks in TDC (Schmithorst et al., 2011). This prior study reports positive correlations of FA around the prefrontal regions with performance in a speech-in-noise task. However, all measures of auditory processing in this previous study involved manipulations of language stimuli, unlike the acoustic test of the DSTP which involves manipulations of non-linguistic stimuli, therefore being more likely to probe basic auditory processing function instead of the integrity of language networks.

To our knowledge, our study is the first to demonstrate associations of white matter microstructure with non-linguistic auditory processing, both among children with SPD as well as among TDC. Both the Sensory Profile auditory score and DSTP are associated with FA in the PTR, which contains the primary auditory projection pathway. However, they are also associated with microstructure in associational tracts such as the external capsules and the cingulate portion of the cingulum bundles. Unlike the Sensory Profile auditory score,

DSTP demonstrates widespread associations with both frontal and posterior projection and commissural pathways, along with the associational tracts of the superior longitudinal fasciculus and sagittal stratum. Similar to tactile processing, one contributor to the extensive regions of significant correlations with the DSTP task may lie in the test's additional recruitment of attentional processes. For example, dichotic listening tasks have long been used to test different neural models of attention (Broadbent, 1954; Treisman, 1964), and have been used to demonstrate enhanced right ear advantage in individuals with mild Alzheimer's disease (Duchek and Balota, 2005), decreased right ear advantage in children with dyslexia (Helland and Asbjornsen, 2001), and challenges with attention shifting in sleep deprivation (Johnson et al., 2002). The dichotic listening portion of the DSTP demonstrates the most robust correlations with FA of the three DSTP subtests, showing strong FA effects across projection, commissural, and association tracts (Table 5). Temporal patterning shows strong effects in certain projection and association tracts. It has previously been suggested that children with specific language impairment perform worse on

<table><caption>TABLE 5 | Number of significantly correlated voxels in several ROIs, along with results of the three DSTP subscores as functions of group, FA, age, and gender.</caption><thead><tr><th rowspan="2"></th><th colspan="3">DSTPdd</th><th colspan="3">DSTPtP</th><th colspan="3">DSTPad</th></tr><tr><th># sig vox</th><th>p_FA</th><th>p_TDCvSPD</th><th># sig vox</th><th>p_FA</th><th>p_TDCvSPD</th><th># sig vox</th><th>p_FA</th><th>p_TDCvSPD</th></tr></thead><tbody><tr><td>ACR-L</td><td>734</td><td><strong>0.016</strong></td><td>0.32</td><td>582</td><td>0.10</td><td><strong>0.010</strong></td><td>-</td><td>-</td><td>-</td></tr><tr><td>ACR-R</td><td>1190</td><td><strong>0.0025</strong></td><td>0.34</td><td>918</td><td><strong>0.028</strong></td><td><strong>0.011</strong></td><td>1176</td><td><strong>0.00057</strong></td><td>0.52</td></tr><tr><td>SCR-L</td><td>334</td><td><strong>0.037</strong></td><td>0.15</td><td>7</td><td>0.23</td><td><strong>0.0022</strong></td><td>-</td><td>-</td><td>-</td></tr><tr><td>SCR-R</td><td>460</td><td><strong>0.022</strong></td><td>0.18</td><td>326</td><td><strong>0.037</strong></td><td><strong>0.0033</strong></td><td>83</td><td><strong>0.029</strong></td><td>0.19</td></tr><tr><td>PCR-L</td><td>255</td><td><strong>0.013</strong></td><td>0.21</td><td>52</td><td>0.27</td><td><strong>0.0061</strong></td><td>-</td><td>-</td><td>-</td></tr><tr><td>PCR-R</td><td>499</td><td><strong>0.0034</strong></td><td>0.27</td><td>341</td><td><strong>0.046</strong></td><td><strong>0.0086</strong></td><td>40</td><td><strong>0.012</strong></td><td>0.25</td></tr><tr><td>ALIC-L</td><td>268</td><td><strong>0.038</strong></td><td>0.36</td><td>316</td><td><strong>0.0074</strong></td><td><strong>0.021</strong></td><td>-</td><td>-</td><td>-</td></tr><tr><td>ALIC-R</td><td>521</td><td><strong>0.013</strong></td><td>0.50</td><td>475</td><td><strong>0.037</strong></td><td><strong>0.028</strong></td><td>326</td><td><strong>0.0054</strong></td><td>0.64</td></tr><tr><td>PLIC-L</td><td>584</td><td><strong>0.0050</strong></td><td>0.21</td><td>306</td><td>0.077</td><td><strong>0.0048</strong></td><td>-</td><td>-</td><td>-</td></tr><tr><td>PLIC-R</td><td>602</td><td><strong>0.0025</strong></td><td>0.30</td><td>338</td><td>0.29</td><td><strong>0.0070</strong></td><td>18</td><td><strong>0.0063</strong></td><td>0.55</td></tr><tr><td>RLIC-L</td><td>455</td><td><strong>0.0091</strong></td><td>0.28</td><td>440</td><td><strong>0.043</strong></td><td><strong>0.0087</strong></td><td>-</td><td>-</td><td>-</td></tr><tr><td>RLIC-R</td><td>453</td><td><strong>0.0083</strong></td><td>0.43</td><td>363</td><td>0.14</td><td><strong>0.012</strong></td><td>111</td><td><strong>0.026</strong></td><td>0.62</td></tr><tr><td>PTR-L</td><td>596</td><td><strong>0.0075</strong></td><td>0.67</td><td>463</td><td>0.32</td><td><strong>0.019</strong></td><td>-</td><td>-</td><td>-</td></tr><tr><td>PTR-R</td><td>894</td><td><strong>0.00043</strong></td><td>0.99</td><td>583</td><td>0.12</td><td><strong>0.031</strong></td><td>305</td><td><strong>0.0072</strong></td><td>0.81</td></tr><tr><td>GCC</td><td>1116</td><td><strong>0.014</strong></td><td>0.33</td><td>928</td><td>0.095</td><td><strong>0.0088</strong></td><td>426</td><td><strong>0.043</strong></td><td>0.44</td></tr><tr><td>BCC</td><td>2048</td><td><strong>0.030</strong></td><td>0.18</td><td>1788</td><td>0.085</td><td><strong>0.0044</strong></td><td>529</td><td><strong>0.041</strong></td><td>0.23</td></tr><tr><td>SCC</td><td>939</td><td><strong>0.028</strong></td><td>0.43</td><td>871</td><td>0.074</td><td><strong>0.021</strong></td><td>129</td><td><strong>0.025</strong></td><td>0.45</td></tr><tr><td>CGC-L</td><td>268</td><td><strong>0.0064</strong></td><td>0.28</td><td>153</td><td><strong>0.025</strong></td><td><strong>0.0091</strong></td><td>-</td><td>-</td><td>-</td></tr><tr><td>CGC-R</td><td>8</td><td><strong>0.11</strong></td><td>0.20</td><td>69</td><td>0.080</td><td><strong>0.0080</strong></td><td>-</td><td>-</td><td>-</td></tr><tr><td>EC-L</td><td>223</td><td><strong>0.13</strong></td><td>0.18</td><td>439</td><td><strong>0.049</strong></td><td><strong>0.0041</strong></td><td>-</td><td>-</td><td>-</td></tr><tr><td>EC-R</td><td>911</td><td><strong>0.0054</strong></td><td>0.23</td><td>598</td><td><strong>0.013</strong></td><td><strong>0.0045</strong></td><td>752</td><td><strong>0.00054</strong></td><td>0.44</td></tr><tr><td>SLF-L</td><td>702</td><td><strong>0.0014</strong></td><td>0.24</td><td>9</td><td>0.14</td><td><strong>0.0024</strong></td><td>-</td><td>-</td><td>-</td></tr><tr><td>SLF-R</td><td>981</td><td><strong>0.00046</strong></td><td>0.18</td><td>715</td><td><strong>0.014</strong></td><td><strong>0.0067</strong></td><td>206</td><td><strong>0.0013</strong></td><td>0.45</td></tr><tr><td>SS-L</td><td>177</td><td><strong>0.023</strong></td><td>0.20</td><td>235</td><td><strong>0.0068</strong></td><td><strong>0.0051</strong></td><td>5</td><td><strong>0.0026</strong></td><td>0.18</td></tr><tr><td>SS-R</td><td>259</td><td><strong>0.0054</strong></td><td>0.56</td><td>284</td><td><strong>0.016</strong></td><td><strong>0.032</strong></td><td>212</td><td><strong>0.00044</strong></td><td>0.98</td></tr></tbody></table>

The significance of group and FA to the GLMs are displayed. The significance of age and gender are described in the main text. DSTPdd, dichotic digits; DSTPtP, temporal patterning; DSTPad, auditory discrimination. ACR, anterior corona radiata; SCR, superior corona radiata; PCR, posterior corona radiata; ALIC, anterior limb of the internal capsule; PLIC, posterior limb of the internal capsule; RLIC, retrolenticular limb of the internal capsule; PTR, posterior thalamic radiation; GCC, genu of the corpus callosum; BCC, body of the corpus callosum; SCC, splenium of the corpus callosum; CGC, cingulate gyrus part of the cingulum; EC, external capsule; SLF, superior longitudinal fasciculus; SS, sagittal stratum. Bolded p values indicate statistically significant effects (p < 0.05).
<figure><img src="image_15.png" /></figure>

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

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

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

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

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

FIGURE 6 | Representative scatterplots and regression lines of the DSTP subscores (mean-centered) – dichotic digits, temporal patterning, and auditory discrimination – vs. FA in significant voxels (differing between the three sensory variables) of the PTR (both sides combined) and splenium of the corpus callosum (SCC).
temporal patterning tasks, but only when the interstimulus interval between tones is short (Tallal and Piercy, 1973). Interestingly, the left RLIC, with which temporal patterning demonstrates an FA effect in our present study (but not with the right RLIC), has been previously associated with reading scores in individuals with dyslexia by a DTI study (Klingberg et al., 2000). Auditory discrimination, which involves a strong element of attentional control, demonstrates the strongest FA effects with the right anterior corona radiata, which is concordant with Schmithorst et al.'s (2011) finding of correlations of language discrimination in noise with FA in white matter of prefrontal regions.

Auditory processing is of primary importance for language acquisition, with speech perception requiring the ability to determine spectral shape, to discriminate modulation of amplitude and spectral frequencies, and to do this at varying temporal resolutions (Bailey and Snowling, 2002). These auditory capabilities are likely sufficient by early infancy to support the discrimination of phonetic elements in language, though attentional capabilities may take longer to develop (Bailey and Snowling, 2002). Given the wealth of evidence that developmental language-based learning disorders can often be traced to non-verbal auditory processing deficits (Tallal et al., 1997), it is highly relevant to characterize the biological underpinnings of these deficits.

## Study Limitations and Future Directions

Despite the larger number of subjects in this study, we are still limited by sample size in our ability to harness DTI as a clinically utilizable tool for the diagnosis, prognosis, and treatment of SPD. Going forward, larger sample sizes and multimodal imaging biomarkers from DTI, fMRI, and MEG may aid in better definition and diagnosis of SPD. This could be of particular use if these biomarkers can identify individuals at risk for SPD

at early ages before clinical symptoms are apparent, allowing for early intervention and recruitment of support services. In addition to diagnosis, larger scale longitudinal studies will allow us to evaluate the utility of quantitative imaging biomarkers, as compared with clinical metrics, neuropsychological testing, or direct sensory testing, for the prognostication of behavioral and cognitive development of individuals with SPD. Finally, interventional studies will allow us to evaluate the utility of quantitative imaging biomarkers for the monitoring of behavioral and psychopharmacological interventions, as well as for the prediction of interventional outcome. These biomarkers may further aid in the design of interventions if they can be used to stratify the SPD population into subgroups that will better respond to particular interventions. Overall, future studies will aim to shift SPD from a clinical diagnosis to a “biomarker diagnosis,” with imaging, and in particular DTI, metrics among the most promising of these biomarkers.

## AUTHOR CONTRIBUTIONS

Conception and design: JO, EM, PM. Acquisition: AB-A, SD, SH, ABA, JH. Analysis: Y-SC, MG, JO, EM, PM. Interpretation: YSC, MG, JPO, EJM, PM.

## ACKNOWLEDGMENTS

This work was funded by grants from the Wallace Research Foundation to EJM and to PM and a gift from Toby Mickelson and Donald Brody to EJM. EJM has also received neuroimaging support that contributed to this work from NIH K23 MH083890. We also received generous support from the SPD community of family and friends through gifts large and small to our UCSF Sensory Neurodevelopment and Autism Program (SNAP).

## REFERENCES

Ahn, R. R., Miller, L. J., Milberger, S., and McIntosh, D. N. (2004). Prevalence of parents' perceptions of sensory processing disorders among kindergarten children. *Am. J. Occup. Ther.* 58, 287-293. doi: 10.5014/ajot.58.3.287

Ayres, A. J. (1989). *Sensory Integration and Praxis Tests: SIPT Manual*. Los Angeles, CA: Western Psychological Services.

Bailey, P. J., and Snowling, M. J. (2002). Auditory processing and the development of language and literacy. *Br. Med. Bull.* 63, 135-146. doi: 10.1093/bmb/63.1.135

Binkofski, F., Buccino, G., Posse, S., Seitz, R. J., Rizzolatti, G., and Freund, H. (1999). A fronto-parietal circuit for object manipulation in man: evidence from an fMRI-study. *Eur. J. Neurosci.* 11, 3276-3286. doi: 10.1046/j.1460-9568.1999.00753.x

Bodegård, A., Geyer, S., Grefkes, C., Zilles, K., and Roland, P. E. (2001). Hierarchical processing of tactile shape in the human brain. *Neuron* 31, 317-328. doi: 10.1016/S0896-6273(01)00362-2

Bodegård, A., Ledberg, A., Geyer, S., Naito, E., Zilles, K., and Roland, P. E. (2000). Object shape differences reflected by somatosensory cortical activation. *J. Neurosci.* 20, RC51.

Broadbent, D. E. (1954). The role of auditory localization in attention and memory span. *J. Exp. Psychol.* 44, 51-55. doi: 10.1037/h0056491

Cascio, C. J. (2010). Somatosensory processing in neurodevelopmental disorders. *J. Neurodev. Disord.* 2, 62-69. doi: 10.1007/s11689-010-9046-3

Chang, Y. S., Owen, J. P., Desai, S. S., Hill, S. S., Arnett, A. B., Harris, J., et al. (2014). Autism and sensory processing disorders: shared white matter disruption in sensory pathways but divergent connectivity in social-emotional pathways. *PLoS ONE* 9:e103038. doi: 10.1371/journal.pone.0103038

Demopoulos, C., Brandes-Aitken, A. N., Desai, S. S., Hill, S. S., Antovich, A. D., Harris, J., et al. (2015). Shared and Divergent Auditory and Tactile Processing in Children with Autism and Children with Sensory Processing Dysfunction Relative to Typically Developing Peers. *J. Int. Neuropsychol. Soc.* 6, 1-11. doi: 10.1017/s1355617715000387

Duchek, J. M., and Balota, D. A. (2005). Failure to control prepotent pathways in early stage dementia of the Alzheimer's type: Evidence from dichotic listening. *Neuropsychology* 19, 687-695. doi: 10.1037/0894-4105.19.5.687

Dunn, W., and Westman, K. (1997). The sensory profile: the performance of a national sample of children without disabilities. *Am. J. Occup. Ther.* 51, 25-34. doi: 10.5014/ajot.51.1.25

Hadjikhani, N., and Roland, P. E. (1998). Cross-modal transfer of information between the tactile and the visual representations in the human brain: a positron emission tomographic study. *J. Neurosci.* 18, 1072-1084.

Helland, T., and Asbjornsen, A. (2001). Brain asymmetry for language in dyslexic children. *Laterality* 6, 289-301. doi: 10.1080/713754422

Hofer, S., and Frahm, J. (2006). Topography of the human corpus callosum revisited – Comprehensive fiber tractography using diffusion tensor magnetic resonance imaging. *Neuroimage* 32, 989-994. doi: 10.1016/j.neuroimage.2006.05.044
Hoon, A. H. Jr., Stashinko, E. E., Nagae, L. M., Lin, D. D., Keller, J., Bastian, A., et al. (2009). Sensory and motor deficits in children with cerebral palsy born preterm correlate with diffusion tensor imaging abnormalities in thalamocortical pathways. *Dev. Med. Child Neurol.* 51, 697–704. doi: 10.1111/j.1469-8749.2009.03306.x

Huang, H., Zhang, J., Jiang, H., Wakana, S., Poetscher, L., Miller, M. I., et al. (2005). DTI tractography based parcellation of white matter: application to the mid-sagittal morphology of corpus callosum. *Neuroimage* 26, 195–205. doi: 10.1016/j.neuroimage.2005.01.019

Jenkinson, M., Beckmann, C. F., Behrens, T. E. J., Woolrich, M. W., and Smith, S. M. (2002). Fsl. *Neuroimage* 62, 782–790. doi: 10.1016/j.neuroimage.2011.09.015

Johnson, B. H., Laberg, J. C., Eid, J., and Hugdahl, K. (2002). Dichotic listening and sleep deprivation. *Scand. J. Psychol.* 43, 413–417. doi: 10.1111/1467-9450.00309

Kawashima, R., Roland, P. E., and O'Sullivan, B. T. (1994). Activity in the human primary motor cortex related to ipsilateral hand movements. *Brain Res.* 663, 251–256. doi: 10.1016/0006-8993(94)91270-X

Klingberg, T., Hedehus, M., Temple, E., Salz, T., Gabrieli, J. D., Moseley, M. E., et al. (2000). Microstructure of temporo-parietal white matter as a basis for reading ability: evidence from diffusion tensor magnetic resonance imaging. *Neuron* 25, 493–500. doi: 10.1016/S0896-6273(00)80911-3

Lord, C., Risi, S., Lambrecht, L., Cook, E. H. Jr., Leventhal, B. L., DiLavore, P. C., et al. (2000). The autism diagnostic observation schedule-generic: a standard measure of social and communication deficits associated with the spectrum of autism. *J. Autism. Dev. Disord.* 30, 205–223. doi: 10.1023/A:1005592 401947

Metcalfe, J. S., McDowell, K., Chang, T. Y., Chen, L. C., Jeka, J. J., and Clark, J. E. (2005). Development of somatosensory-motor integration: and event-related analysis of infant posture in the first year of independent walking. *Dev. Psychobiol.* 46, 19–35. doi: 10.1002/dev.20037

Montagu, A. (1986). *Touching: the Human Significance of the Skin*. (New York: Perennial Library).

O'Sullivan, B. T., Roland, P. E., and Kawashima, R. (1994). A PET study of somatosensory discrimination in Man. *Microgeometry versus macrogeometry*. *Eur. J. Neurosci.* 6, 137–148. doi: 10.1111/j.1460-9568.1994.tb00255.x

Owen, J. P., Marco, E. J., Desai, S., Fourie, E., Harris, J., Hill, S. S., et al. (2013). Abnormal white matter microstructure in children with sensory processing disorders. *Neuroimage Clin.* 2, 844–853. doi: 10.1016/j.nicl. 2013.06.009

Penke, L., Muñoz Maniega, S., Murray, C., Gow, A. J., Hernández, M. C., Clayden, J. D., et al. (2010). A general factor of brain white matter integrity predicts information processing speed in healthy older people. *J. Neurosci.* 30, 7569– 7574. doi: 10.1523/JNEUROSCI.1553-10.2010

Pryweller, J. R., Schauder, K. B., Anderson, A. W., Heacock, J. L., Foss-Feig, J. H., Newsom, C. R., et al. (2014). White matter correlates of sensory processing in autism spectrum disorders. *Neuroimage Clin.* 6, 379–387. doi: 10.1016/j.nicl.2014.09.018

Richard, G., and Ferre, J. (2006). *Differential Screening Test for Processing*. East Moline, IL: LinguiSystems, Inc.

Roland, P. E., and Zilles, K. (1998). Structural divisions and functional fields in the human cerebral cortex. *Brain Res. Rev.* 26, 87–105. doi: 10.1016/S0165- 0173(97)00058-1

Rutter, M., Bailey, A., and Lord, C. (2003). *Social Communication Questionnaire (SCQ)*. Los Angeles, CA: Western Psychological Services.

Schmithorst, V. J., Holland, S. K., and Plante, E. (2011). Diffusion tensor imaging reveals white matter microstructure correlations with auditory processing ability. *Ear Hear.* 32, 156–167. doi: 10.1097/AUD.0b013e3181f7a481

Servos, P., Lederman, S., Wilson, D., and Gati, J. (2001). fMRI-derived cortical maps for haptic shape, texture, and hardness. *Cogn. Brain Res.* 12, 307–313. doi: 10.1016/S0926-6410(01)00041-6

Smith, S. M., Jenkinson, M., Johansen-Berg, H., Rueckert, D., Nichols, T. E., Mackay, C. E., et al. (2006). Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion data. *Neuroimage* 31, 1487–1505. doi: 10.1016/j.neuroimage.2006.02.024

Smith, S. M., and Nichols, T. E. (2009). Threshold-free cluster enhancement: addressing problems of smoothing, threshold dependence and localisation in cluster inference. *Neuroimage* 44, 83–98. doi: 10.1016/j.neuroimage.2008.03.061

Soechting, J. F., and Flanders, M. (2008). Sensorimotor control of contact force. *Curr. Opin. Neurobiol.* 18, 565–572. doi: 10.1016/j.conb.2008.11.006

Stoeckel, M. C., Weder, B., Binkofski, F., Buccino, G., Jon Shah, N., and Seitz, R. J. (2003). A fronto-parietal circuit for tactile object discrimination: an event-related fMRI study. *Neuroimage* 19, 1103–1114. doi: 10.1016/S1053- 8119(03)00182-4

Stoeckel, M. C., Weder, B., Binkofski, F., Choi, H. J., Amunts, K., Pieperhoff, P., et al. (2004). Left and right superior parietal lobule in tactile object discrimination. *Eur. J. Neurosci.* 19, 1067–1072. doi: 10.1111/j.0953-816X.2004.03185.x

Stoesz, M. R., Zhang, M., Weisser, V. D., Prather, S. C., Mao, H., and Sathian, K. (2003). Neural networks active during tactile form perception: common and differential activity during macrospatial and microspatial tasks. *Int. J. Psychophysiol.* 50, 41–49. doi: 10.1016/S0167-8760(03)00123-5

Tallal, P., Miller, S. L., Jenkins, W. M., and Merzenich, M. M. (1997). “The role of temporal processing in developmental language-based learning disorders: research and clinical implications,” in *Foundations of Reading Acquisition and Dyslexia: Implications for Early Intervention*, ed. B. A. Blachman (Mahwah, NJ: Lawrence Erlbaum), 49–66.

Tallal, P., and Piercy, M. (1973). Developmental aphasia: impaired rate of non-verbal processing as a function of sensory modality. *Neuropsychologia* 11, 389–398. doi: 10.1016/0028-3932(73)90025-0

Treisman, A. (1964). Verbal cues, language, and meaning in selective attention. *Am. J. Psychol.* 77, 206–209. doi: 10.2307/1420127

Van de Winckel, A., Sunaert, S., Wenderoth, N., Peeeters, R., Van Hecke, P., Feys, H., et al. (2005). Passive somatosensory discrimination tasks in healthy volunteers: differential networks involved in familiar versus unfamiliar shape and length discrimination. *Neuroimage* 26, 441–453. doi: 10.1016/j.neuroimage.2005.01.058

Wahl, M., Li, Y. O., Ng, J., Lahue, S. C., Cooper, S. R., Sherr, E. H., et al. (2010). Microstructural correlations of white matter tracts in the human brain. *Neuroimage* 51, 531–541. doi: 10.1016/j.neuroimage.2010.02.072

Wakana, S., Jiang, H., Nagae-Poetscher, L. M., Van Zijl, P. C. M., and Mori, S. (2004). Fiber Tract-based Atlas of Human White Matter Anatomy. *Radiology* 230, 77–87. doi: 10.1148/radiol.2301021640

Wechsler, D. (2003). *Wechsler Intelligence Scale for Children*, 4th Edn. San Antonio, TX: Harcourt Assessment.

Zero to Three (2005). *Diagnostic Classification of Mental Health and Developmental Disorders of Infancy and Early Childhood: Revised Edition*. Washington, DC: ZERO TO THREE Press.

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2016 Chang, Gratiot, Owen, Brandes-Aitken, Desai, Hill, Arnett, Harris, Marco and Mukherjee. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.