Skip to main content

ORIGINAL RESEARCH article

Front. Hum. Neurosci., 27 November 2018
Sec. Brain Imaging and Stimulation
Volume 12 - 2018 | https://doi.org/10.3389/fnhum.2018.00474

Functional and Structural Brain Plasticity in Adult Onset Single-Sided Deafness

  • 1Department of Otolaryngology-Head and Neck Surgery, University of California, San Francisco, San Francisco, CA, United States
  • 2Department of Otorhinolaryngology, Peking Union Medical College Hospital, Beijing, China
  • 3Department of Radiology and Biomedical Imaging, University of California, San Francisco, San Francisco, CA, United States
  • 4Department of Psychiatry, University of California, San Francisco, San Francisco, CA, United States

Single-sided deafness (SSD) or profound unilateral hearing loss obligates the only serviceable ear to capture all acoustic information. This loss of binaural function taxes cognitive resources for accurate listening performance, especially under adverse environments or challenging tasks. We hypothesized that adults with SSD would manifest both functional and structural brain plasticity compared to controls with normal binaural hearing. We evaluated functional alterations using magnetoencephalographic imaging (MEGI) of brain activation during performance of a moderately difficult auditory syllable sequence reproduction task and assessed structural integrity using diffusion tensor imaging (DTI). MEGI showed the SSD cohort to have increased induced oscillations in the theta band over the left superior temporal cortex and decreased induced gamma band oscillations over the frontal and parietal cortices between 175 and 475 ms following stimulus onset. DTI showed the SSD cohort to have extensive fractional anisotropy (FA) reduction in both auditory and non-auditory tracts and regions. Overlaying functional and structural changes revealed by the two imaging techniques demonstrated close registration of cortical areas and white matter tracts that expressed brain plasticity. Hence, complete loss of input from one ear in adulthood triggers both functional and structural alterations to dorsal temporal and frontal-parietal areas.

Introduction

Bilateral deafness or anacusis enforces undesirable auditory isolation and triggers widespread central network plasticity. Those with deafness in both ears compared to controls with normal hearing (NH) have increased activation of auditory association cortex when attending to visual stimuli and performing working memory tasks (Finney et al., 2001, 2003; Lomber et al., 2010; Meredith et al., 2011; Ding et al., 2015) and can develop superior visual perception skills (Dye et al., 2009; Ding et al., 2015). In subjects who acquire bilateral deafness early, the extent of plasticity in auditory cortex appears to be correlated with level of performance enhancement on visuospatialtasks (Ding et al., 2015). Here, we examine whether widespread brain plasticity observed in bilateral deafness may be similarly expressed in unilateral deafness, despite continued access to sound in the only hearing ear.

Unilateral deafness or profound hearing loss is commonly referred to as single-sided deafness (SSD). It is prevalent in adults and children alike and may be congenital or acquired. SSD differs from bilateral deafness in that the only hearing ear can capture sound information for further processing by the central auditory system and related networks. However, an important functional constraint of monaural hearing is loss of panoramic sound access due to the head shadow effect. Some behavioral adaptations to single sensor hearing are deployment of cognitive and attentional resources to monitor the acoustic scene, readiness to orient to sounds of interest by making head or body adjustments, and use of visually guided lip-reading strategies to bolster accuracy. Despite those steps to mitigate challenges imposed by SSD, troublesome clinical consequences include greater difficulty with sound localization, higher signal-to-noise ratio requirement for speech discrimination, and reduced ease of listening (Rigby et al., 1997; Noble and Gatehouse, 2004; Tufarelli et al., 2006; Douglas et al., 2007). Those and other chronic listening difficulties in SSD may stress cognitive networks that subserve attention, working memory, and executive function. Children with SSD exhibit delayed acquisition of oral language and verbal intelligence skills, decreased achievement in school, and impaired executive control when challenged by irrelevant verbal information (Bess and Tharpe, 1984; Brookhouser et al., 1991; Borg et al., 2002; Lieu et al., 2012; Vila and Lieu, 2015). Current treatments for SSD include routing acoustic information to the better ear (Peters et al., 2015), transferring sounds to the contralateral cochlea by bone conduction (Kim et al., 2017), and stimulating the deaf ear by cochlear implantation (Vlastarakos et al., 2014). While the latter intervention reintroduces auditory information through the deaf ear, with promise for restoration of binaural hearing, treatment outcomes remain highly variable (van Zon et al., 2015; Cabral Junior et al., 2016).

Given the wide ranging clinical deficits observed in SSD and lack of reliable predictors to assist with patient selection for treatment by cochlear implantation, there is a need for greater understanding of the neurophysiological and neuroanatomical consequences of unilateral deafness, and more broadly, asymmetric hearing loss. In non-human primates with noise-induced asymmetric hearing loss, microelectrode-derived primary auditory cortical maps show plastic change in canonical microcircuits for the sound processing. That is, auditory cortex contralateral to the poorer ear undergoes time-dependent realignment of interaural frequency maps in association with elevated cortical thresholds to sound input from the better ear (Cheung et al., 2009). A follow-up study that contrast functional organization of the contralateral and ipsilateral hemispheres relative to the poorer ear demonstrates loss of normal mirror-image relationship across the two hemispheres, suggesting that whole brain alterations triggered by asymmetric hearing loss may account for impaired performance in tasks that require binaural input alignment or interhemispheric processing (Cheung et al., 2017). Parallel human brain imaging studies that aim to identify brain plasticity in SSD asymmetric hearing loss are limited. Magnetoencephalographic imaging (MEGI) studies that have used simple sound stimuli and examined low frequency phase-locked evoked activity in auditory cortices of adults with SSD. Interhemispheric differences in the M100 response provide evidence for temporal and spatial plasticity (Pross et al., 2015; Chang et al., 2016). Beyond auditory cortex, hearing restricted to one ear may have wider impact on functional networks that support attention and working memory demands during difficult listening tasks of complex acoustic stimuli. Cortical plasticity of such non-phase locked neural oscillatory activity across multiple frequency bands have to date not yet been examined in SSD.

In this study, we evaluate stimulus-induced neural oscillatory activity differences between adult onset SSD subjects with one NH ear against controls with two normal ears using a moderately difficult auditory syllable sequence reproduction task. We use MEGI to measure power fluctuations of non-phase locked activity across the whole brain for the entire spectrum of neural oscillations. Moreover, we use diffusion tensor imaging (DTI) to assess whole brain white matter structural integrity. We hypothesize that restriction of auditory input to one ear in adults with SSD results in both functional plasticity during accurate listening and structural plasticity of attention and working memory networks.

Materials and Methods

Subjects

This study protocol was approved by the Committee on Human Research at University of California, San Francisco (UCSF). Consent obtained from all subjects to participate in the study was both informed and written. Twenty-six adult onset SSD subjects [all values: mean ± standard deviation; age (years) = 48.7 ± 14.6; 12 female; one left handed] with unilateral hearing loss longer than 1 year duration were recruited from UCSF Otolaryngology and Audiology clinics. Of these 26 subjects, 15 were deaf in the left ear and 11 were deaf in the right ear. The average duration of deafness was 8.6 ± 8.7 years, ranging from 1 to 31 years. The average audiometric threshold of pure tone frequencies at 0.5, 1, 2, and 4 kHz for the deaf ear was 96.8 ± 18.4 dB HL, and the only hearing ear was 6.1 ± 5.8 dB HL. The etiology of sensorineural SSD was idiopathic in 18 subjects and acoustic neuroma in eight subjects that was managed by observation or gamma knife radiosurgery (no open craniotomy). Fifteen NH control subjects were recruited [age (years) = 46.6 ± 12.4; seven female; one left handed] from the UCSF community, matching age, gender, and handedness to the cohort of SSD patients (Table 1). The average audiometric threshold of pure tone frequencies at 0.5, 1, 2, and 4 kHz in control subjects was ≤25 dB HL in both ears. Exclusion criteria included implanted bone conduction devices, significant head trauma, seizure disorder, cerebrovascular accident, magnetic artifact, and contraindication to MRI. After data quality assessments and exclusions were applied, MEGI data was analyzable in 24 SSD patients and 14 NH controls and DTI data was analyzable in 20 SSD patients and 12 NH controls.

TABLE 1
www.frontiersin.org

TABLE 1. Demographic data of SSD patients and normal hearing controls.

MEGI Acquisition and Task

Magnetoencephalographic imaging data were collected using a whole-head 275-channel CTF Omega 2000 biomagnetometer with third-order gradient correction (CTF Inc., Port Coquitlam, Canada), at 1200 Hz sampling rate. Radio-emitting coils were placed at the nasion and 1 cm anterior to the left and right preauricular points for head position tracking and co-registration with MRI (MPRAGE) data. MEGI sessions in which head movement exceeded 2 mm were discarded and repeated to maximize high data quality.

During MEGI, subjects performed an auditory syllable reproduction task (Herman et al., 2013). Subjects were instructed to listen to and vocally reproduce two- or four-syllable utterances presented binaurally during the imaging session. Stimuli were pre-recorded from a single female speaker and consisted of permutations of the syllables /ba/, /da/, and /pa/ (e.g., /ba da/ or /ba da pa ba/). At the beginning of each trial, subjects listened to the target presentation, waited for a visual “go” cue to respond (jittered between 2050 and 2150 ms after stimulus onset), and then vocally repeated back the stimulus pattern they heard. Syllables lasted 470 ms on average and were separated by 150 ms within a trial. After a go cue, subjects had up to 3 s to complete their response, after which the experiment proceeded to the next trial. Within two-syllable trials, no syllable was repeated, and within four-syllable trials, no syllable pair was repeated. Subjects had no previous knowledge or clues about the content of the upcoming trial. Eighty two- and four-syllable target trials were presented in pseudorandomized order. Correct and incorrect syllable repetition was recorded. Trials were labeled correct if the target syllables were repeated in proper order within 3 s of the go cue. Behavioral data (reaction time, accuracy) were analyzed using 2 × 2 ANOVA, by cohort (SSD, NH) and condition (two-syllable, four-syllable).

MRI Acquisition

Structural MRI and DTI data were acquired from each subject on a 3-T MRI scanner (Discovery MR750, GE Medical system, Waukesha, WI, United States). Structural T1-weighted FSPGR BRAVO scans were obtained (120 axial slices, FOV = 512 mm × 512 mm, TR = 7232 ms, TE = 2.78 ms, in-plane voxel dimensions 0.5 mm × 0.5 mm, slice thickness = 1.5 mm). For the DTI protocol, a 12-min DTI (HARDI) scan was collected in 55 diffusion weighted directions (b = 3000 s/mm2, TR = 9313 ms, TE = 76.6 ms, FOV = 256 mm × 256 mm, in-plane voxel dimensions = 2 mm × 2 mm, axial slice thickness = 2 mm).

MRI-DTI Analysis

Preprocessing of DTI data was performed with FSL (FMRIB Software Library, version 4.1.9, Oxford, United Kingdom) (Smith et al., 2004; Jenkinson et al., 2012), including correction for eddy current induced distortion and head motion. Diffusion-weighted images were registered to the non-diffusion weighted volume (b = 0) with six degrees of transformation (Smith, 2002). Fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD) maps were analyzed using tract-based spatial statistics (TBSS) in FSL (Smith et al., 2006). The diffusion tensor was estimated on a voxel-by-voxel basis using FSL’s DTIfit toolbox within the FMRIB Diffusion Toolbox of FSL (FSLv4.1.71). The FA maps were then transformed in Montreal Neurological Institute (MNI) space to an anatomical template (FMRIB58_FA) using non-linear registration in FSL’s FNIRT. These parameters were then applied to the MD, AD, and RD maps for each subject.

A study-specific “skeleton,” representing the center of all fiber bundles common to all subjects was generated using an optimized FA threshold of 0.2 and confirmed by visual inspection. Each subject’s aligned FA, MD, AD, and RD maps were then projected onto this skeleton and entered into voxel-wise cross-subject statistics. Comparisons (independent t-tests) between groups (SSD and NH) were performed using permutation testing (5000 random permutations) applied to the general linear model implemented in FSL. Statistical maps were thresholded at p < 0.01 corrected for multiple comparisons at a cluster level using the TFCE (threshold-free cluster enhancement) approach. This threshold was relaxed in order to maximize conjunction with MEGI source reconstructions (see below), although significant group effects were still present at the p < 0.05 TFCE corrected level. Anatomical labels for white matter tracts were identified using the Johns Hopkins University (JHU) Atlas. To identify cortical structures impacted by white matter tract changes, we also used the Harvard-Oxford cortical atlas.

MEGI Analysis

Magnetoencephalographic imaging data were epoched into -1 to 7 s trials relative to onset of the first syllable. Channels and trials with high frequency activity consistently >1.5 pT or in which the subject spoke during this interval were discarded. A forward model lead field describing the magnetic field strength at each sensor arising from a dipole source at each voxel was computed using a multiple local-spheres spherical volume conductor model (Dalal et al., 2008). All data were denoised using a novel algorithm, dual signal subspace projection (DSSP) for removing large interference and movement artifacts arising from speech production (Sekihara et al., 2016). Default choices were made for the spatial subspace dimension (N = 30) and temporal subspace dimension (N = 5). To avoid mislocalizations attributable to temporally correlated sources between the two hemispheres, data from sensors covering each hemisphere were analyzed separately (Dalal et al., 2008, 2011). The focus of our analysis was on induced neural oscillations during auditory stimulus encoding across the neocortical mantle. We evaluated all induced oscillations across all frequencies. After notch filtering ∼60 Hz, we segmented the data into four bands [4–8 Hz (θ), 8–13 Hz (α), 13–30 Hz (β), and 30–50 Hz (γ)] with a 60 and 120 Hz 1.5-Hz notch filter. Induced, phase-independent activity in each band was localized to the subject’s spatially normalized MRI using the NUTMEG time–frequency beam-forming spatially adaptive filter algorithm, which has been previously described in detail (Dalal et al., 2011). A time–frequency optimized beam-forming inverse solution for the dipole moment dependent on the lead field and sensor covariance was then computed for each voxel for each frequency band at every time window, averaged across overlapping time windows. Localizations were computed using the shared computing cluster at the California Institute for Quantitative Biomedical Research.2 For activations, noise-corrected pseudo-F ratios were computed between active windows (i.e., stimulus) and a prestimulus baseline. Window sizes were frequency-band optimized (4–8 Hz: 400 ms; 8–13 Hz: 300 ms; 13–30 Hz: 200 ms; 30–50 Hz: 150 ms) with an overlap of 50 ms. The time windows centered from 75 to 475 ms after the beginning of stimulus were analyzed. Activations were computed from averaged single-trial data covariance for each time window and frequency band. Since we were primarily interested in induced oscillations during stimulus encoding, we grouped together two- and four-syllable trials to increase power of the response patterns.

For group analyses, tomographic volumes of potential dipolar source locations (voxels) were first generated in the subject’s native MRI and subsequently normalized to the template space. To enable neural-source localization, high-resolution anatomical MR images for each subject were spatially normalized to a standard MNI brain template using SPM23 software. SnPM (statistical non-parametric mapping) unpaired t-tests were performed to compare the difference between SSD and NH cohorts with noise-corrected pseudo-F ratios (Dalal et al., 2008). In brief, time–frequency beam-formed images for each subject were first spatially normalized to the MNI template. The three-dimensional average and variance maps across subjects were calculated for each time–frequency window, and variance maps were smoothed with a 20 mm × 20 mm × 20 mm Gaussian kernel. From this image, a pseudo-t statistic was obtained at each voxel, time window, and frequency band. Non-parametric null distributions were created by permuting voxel labels (2N permutation, where N is the number of subjects) to derive p-values for the true image that were then corrected for multiple comparisons across all voxels, frequency bands, and time points using the cluster correction procedure. The threshold for significant changes was set to k > 20 whole brain cluster level with a cluster building threshold of p < 0.001 uncorrected at voxel level.

Results

Behavioral Performance

In both cohorts, subject accuracy was higher (p < 0.05) and reaction times were faster (p < 0.05) for two-syllable trials compared to four-syllable trials, consistent with previous reports on this task (Herman et al., 2013). Comparing SSD and NH cohorts, there were no significant differences in performance on the two- or four-syllable repetition trials (Figure 1, all p > 0.05). Based on comparable psychoacoustic performance in the two cohorts, we interpret functional and structural imaging differences between SSD and NH as whole brain changes attributable to SSD not confounded by task performance differences.

FIGURE 1
www.frontiersin.org

FIGURE 1. Comparison of reaction times and accuracy rates to two- and four-syllable phoneme tasks between single-sided deafness (SSD) and normal hearing control cohorts. (A,B) Reaction times and accuracy rates, respectively, are indistinguishable. ms, milliseconds; error bars, standard error of the mean.

Neural Activation During Auditory Encoding

We first reconstructed induced oscillatory neural activity during the stimulus period, computing averages with respect to a pre-stimulus baseline for SSD patients and NH controls (Figure 2). To illustrate that the two cohorts share the same activation pattern, we reconstructed time–frequency domain neural activity patterns for each cohort, 14 SSD patients and 14 NH controls, pairwise matched for age, gender, and handedness. This within group analysis revealed that the two cohorts displayed qualitatively similar activation patterns across all frequency bands. During stimulus encoding, induced power changes were positive in the theta and alpha bands from 50 to 200 ms post-stimulus, but were negative in the beta and gamma bands from 150 to 500 ms post-stimulus. However, the across group contrast of 24 SSD patients and 14 NH controls revealed significant quantitative differences in induced oscillatory activity in the theta (4–8 Hz) and gamma (30–55 Hz) bands, mainly over the left hemisphere (Figure 2 and Table 2). Specifically, between 175 and 475 ms following stimulus onset, the SSD cohort showed increased induced theta band oscillations in the left superior temporal cortex and decreased induced gamma band oscillations in the left frontal, superior parietal, and occipital cortices (Figure 2 and Table 2).

FIGURE 2
www.frontiersin.org

FIGURE 2. Induced theta band (4–8 Hz) and gamma band (30–50 Hz) neural oscillations during syllable encoding in normal hearing control and single-sided deafness (SSD) cohorts, and differences between the two cohorts. (A) Theta and (B) gamma band oscillations show similar recruitment of cortical regions in the two cohorts (left and middle columns). For within group contrasts, warm colors indicate increased magnitude relative to a pre-stimulus baseline in the theta band and cool colors indicate increased magnitude relative to a pre-stimulus baseline in the gamma band. (A) SSD cohort has increased recruitment (right column) of higher order auditory cortex in the theta band. (B) SSD cohort has decreased recruitment attention and working memory networks (frontal, parietal, and occipital cortical regions) in the gamma band. In SSD and control cohort contrasts, warm colors in the theta band indicate regions where activity is increased in SSD and cool colors in the gamma band indicate regions where activity is decreased in SSD. ms, milliseconds.

TABLE 2
www.frontiersin.org

TABLE 2. Tabulation of all brain regions showing significant differences between SSD and control cohorts that include a cluster size >20 voxels and uncorrected p-value <0.001.

White Matter Structural Integrity

Tract-based spatial statistics analysis showed the SSD cohort to have significantly reduced FA for multiple white matter skeleton clusters (Figure 3) and increased MD and RD for more extensive tracts that overlapped with the region of reduced FA (Figures 4, 5). By contrast, there were only modest increases in AD that were restricted to a very limited area (data not shown). When mapped to the JHU white matter atlas, FA reduction was observed in the genu, body and splenium of the corpus callosum, fornix, bilateral anterior, superior and posterior corona radiata, anterior and posterior limb of internal capsule, posterior thalamic radiation (including optic radiation), fornix (cres)/stria terminalis, superior longitudinal fasciculus, cingulum (cingulate gyrus), and right superior fronto-occipital fasciculus and external capsule (Table 3). When mapped to the Harvard-Oxford cortical atlas, FA reduction was identified in the frontal pole, insular cortex, superior, middle and inferior frontal gyruses, precentral gyrus, posterior division of superior temporal gyrus, temporo-occipital part of middle and inferior temporal gyrus, postcentral gyrus, superior parietal lobule, supramarginal gyrus, angular gyrus, lateral occipital cortex, intracalcarine cortex, frontal medial cortex, subcallosal cortex, paracingulate gyrus, cingulate gyrus, precuneus cortex, cuneal cortex, frontal orbital cortex, lingual gyrus, occipital fusiform gyrus, frontal operculum cortex, central operculum cortex, planum polare, Heschl’s gyrus, planum temporale, supracalcarine cortex, and occipital pole (Table 3). The SSD cohort showed reduction in FA across many regions, indicating deafness in one ear broadly impacted white matter structural integrity throughout the brain.

FIGURE 3
www.frontiersin.org

FIGURE 3. Tract-based spatial statistics for fractional anisotropy (FA). Widespread tracts with decreased FA in SSD patients compared to control subjects (blue) are superimposed on the skeleton (green) of tracts common to both cohorts. The z coordinate in the Montreal Neurological Institute (MNI) space for each axial section is noted at bottom. L, left; R, right.

FIGURE 4
www.frontiersin.org

FIGURE 4. Tract-based spatial statistics for medial diffusivity (MD). Widespread tracts with increased MD in SSD patients compared to control subjects (red) are superimposed on the skeleton (green) of tracts common to both cohorts. The z coordinate in the MNI space for each axial section is noted at bottom.

FIGURE 5
www.frontiersin.org

FIGURE 5. Tract-based spatial statistics for radial diffusivity (RD). Widespread tracts with increased RD in SSD patients compared to control subjects (red) are superimposed on the skeleton (green) of tracts common to both cohorts. The z coordinate in the MNI space for each axial section is noted at bottom.

TABLE 3
www.frontiersin.org

TABLE 3. White matter regions with decreased FA in SSD patients compared to the normal hearing controls.

Overlay of MEGI and DTI Maps

To examine functional changes in the context of white matter structural integrity in SSD, we combined gamma band findings on MEGI with white matter results from DTI by generating an overlay of these two statistical maps (Figure 6). In the SSD cohort, there were two regions where reductions in gamma (30–55 Hz) power (Figure 6, green) overlapped with reductions in white matter integrity (Figure 6, blue). First, the cluster of reduced gamma power over the left frontal pole was adjacent to reduced FA in the superior longitudinal fasciculus. Second, the region of reduced gamma power over the superior division of parietal cortex on MEGI was also adjacent to the superior longitudinal fasciculus. The SSD cohort had reduced gamma band oscillatory power over the frontal and parietal cortices that were accompanied by impaired integrity of related white matter tracts. Long-term obligatory monaural hearing impacted cortical induced activity and subcortical white matter microstructure of the dorsal attention network (frontal and parietal regions), providing evidence for co-occurrence of functional and structural plasticity in SSD.

FIGURE 6
www.frontiersin.org

FIGURE 6. Relationship between functional and structural differences in single-sided deafness (SSD). In the SSD cohort, tracts with decreased FA are subjacent to cortical regions with decreased activation during performance of the auditory working memory task (gamma band, time window 375 ms). The z and x coordinates in the MNI space are shown at the bottom of axial and sagittal images. (A) Axial. (B) Sagittal.

Discussion

In a first of its kind adult SSD multimodal structure-function imaging study, we measured cortical activity during an auditory stimulus encoding task on millisecond time and sub-centimeter spatial resolution scales and assessed whole brain white matter structural integrity. We identified functional plasticity in SSD, manifested by enhanced activation of higher order auditory cortices in the left hemisphere, suggesting hyperactivity to support comparable performance accuracy relative to NH in a difficult listening task. We also observed corresponding reductions in task-related induced gamma band oscillatory activity in the left frontal and parietal cortices, brain regions implicated in supporting attention and working memory (Rauschecker and Scott, 2009; Rauschecker, 2012, 2018), findings consistent with disengagement of these cortical networks in SSD. Documented functional changes were associated with structural changes in white matter tracts neighboring these regions.

Individuals afflicted by SSD often face considerable auditory processing demands posed by rapid sequential or overlapping stimuli in real-life communicative interactions. We assessed the neural correlates of monaural hearing under cognitive duress by choosing a demanding auditory syllable encoding task that challenged working memory. We established comparable reaction time and accuracy performance in the two cohorts to remove any confounds introduced by unequal hearing competencies. Therefore, our study design directly queried changes in oscillatory network activity due to unilateral deafness.

In earlier adult SSD studies, we found evidence for functional plasticity of early sensory cortices using MEGI. Auditory cortical responses to simple tones revealed spatial plasticity, expressed as increased activation spread distance in the hemisphere contralateral to the only hearing ear and decreased distance in the ipsilateral hemisphere (Chang et al., 2016) and temporal plasticity, expressed as loss of normal interhemispheric latency difference of the M100 signal (Pross et al., 2015). In the current study, the SSD cohort expressed enhanced theta-band oscillations of higher order auditory cortex during encoding of multisyllabic sequences, showing increased recruitment for accurate listening performance and another facet of within-modality sensory cortical plasticity.

Our MEGI observation of task-related gamma band activity in frontal, parietal, and occipital regions is consistent with published accounts of electrocorticography recordings during a working memory task that show recruitment of these regions (Jerbi et al., 2010; Ossandon et al., 2011, 2012; Ramot et al., 2012). Attention and working memory demands of the syllable sequence reproduction task activated these brain networks in both cohorts. However, the SSD cohort expressed reduced gamma band activity in those networks, a broader manifestation of cortical plasticity revealed by the phoneme replication task. This finding suggests that chronically overtaxed cognitive subsystems in SSD can no longer allocate comparable resources relative to NH subjects to perform the challenging auditory task. Taken together, decreased recruitment of overburdened attention and working memory network along with increased recruitment of higher order auditory cortices may be neuroimaging biomarkers for impairment of ease of communication in SSD (Rigby et al., 1997; Gatehouse and Noble, 2004).

Another form of cortical plasticity, namely, cross-modality audiovisual, has been found in few task-based and resting-state fMRI studies. In children with SSD, lower activation of secondary visual processing regions on an audiovisual task has been reported (Schmithorst et al., 2014). In adults with SSD, resting-state fMRI reveals decreased functional connectivity between primary auditory cortex and visual cortex (Liu et al., 2015) and enhanced synchronous output entropy connectivity between left primary auditory cortex with certain regions of visual networks (Zhang et al., 2016). The current study did not explicitly explore audiovisual integration, and our data cannot exclude the possibility of such plasticity expressed as change in induced neural oscillation.

We also evaluated whether the observed functional modifications may be accompanied by structural changes. There were widespread reductions of white matter integrity in multiple tracts, including those that connect regions of altered activity as revealed by MEGI. First, FA reduction was extensive and observed not only in auditory tracts but also in multiple non-auditory tracts. Second, tracts with FA reduction were mainly located in the dorsal cortical regions. It has been well accepted that there are ventral and dorsal processing streams in auditory cortex, responsible for the perception of “what” and “where” aspects of sound processing (Rauschecker and Tian, 2000; Hickok and Poeppel, 2004; Rauschecker and Scott, 2009). We surmise that the dorsal auditory pathways are more likely impacted in SSD. Furthermore, with combined functional and structural alterations revealed by a challenging auditory working memory task, immediate follow up questions on the causal relationships and sequence of events between structural and functional brain plasticity in SSD will require additional studies.

Whereas several studies have used DTI to investigate structural changes in SSD, they focused mainly on auditory pathways; much less attention was devoted to non-auditory pathways (Wu et al., 2009; Rachakonda et al., 2014; Vos et al., 2015). We deployed whole brain TBSS analyses without predefining voxels or tracts of interest (Smith et al., 2006). To the best of our knowledge, no SSD study has used this approach for DTI analyses. Our analyses revealed the SSD cohort had significant reductions of FA, driven by increased MD and RD but not AD in several non-auditory tracts, implicating axonal demyelination as opposed to fiber loss as the main mechanism of structural change. Other SSD studies in the literature also report decreased FA, with inconsistent effects on other diffusion parameters (Lin et al., 2008; Wu et al., 2009; Miao et al., 2013; Hribar et al., 2014). Differences between the present study and prior studies may be due to variations in study populations, demographic features, onset of deafness, and congenital versus acquired deafness.

We recognize important limitations to this study. The number of subjects in both SSD and NH cohorts were relatively modest, resulting in reduced sensitivity to make group level inferences that may be more significant for diffusion imaging data. And, inadequate number of subjects to address right versus left deafness in the SSD cohort necessitated analysis of the entire SSD cohort without differentiation of deafness laterality to achieve suitable statistical power. The literature indicates that aspects brain plasticity in SSD may be dependent on side of deafness (Schmithorst et al., 2014; Zhang et al., 2016). Furthermore, covariates of duration of deafness, age of deafness onset, labyrinthine dysfunction, namely tinnitus and vertigo, were not explicitly studied. Those clinical features may also alter functional and structural manifestations of brain plasticity (Schmidt et al., 2013; Benson et al., 2014; Seydell-Greenwald et al., 2014; Hinkley et al., 2015; Ryu et al., 2015). Future investigations should include much larger cohorts with more detailed clinical assessments to address both issues.

Author Contributions

SN, PM, and SC contributed to the conception and design of the study. CG and DM organized the database. YS, LH, CC, KS, Y-SC, and JO performed the statistical analysis. YS wrote the first draft of the manuscript. SN and SC wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.

Funding

This work was supported by Ministry of Education of China (YS), Steven M. Bauer Research Fund (SC), Hearing Research, Inc. (SC), and Coleman Memorial Fund (SC).

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.

Footnotes

  1. ^ https://www.fmrib.ox.ac.uk/fsl
  2. ^ www.qb3.org
  3. ^ http://www.fil.ion.ucl.ac.uk/spm/software/spm2/

References

Benson, R. R., Gattu, R., and Cacace, A. T. (2014). Left hemisphere fractional anisotropy increase in noise-induced tinnitus: a diffusion tensor imaging (DTI) study of white matter tracts in the brain. Hear. Res. 309, 8–16. doi: 10.1016/j.heares.2013.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Bess, F. H., and Tharpe, A. M. (1984). Unilateral hearing impairment in children. Pediatrics 74, 206–216.

Google Scholar

Borg, E., Risberg, A., McAllister, B., Undemar, B. M., Edquist, G., Reinholdson, A. C., et al. (2002). Language development in hearing-impaired children. Establishment of a reference material for a ‘Language test for hearing-impaired children’, LATHIC. Int. J. Pediatr. Otorhinolaryngol. 65, 15–26. doi: 10.1016/S0165-5876(02)00120-9

CrossRef Full Text | Google Scholar

Brookhouser, P. E., Worthington, D. W., and Kelly, W. J. (1991). Unilateral hearing loss in children. Laryngoscope 101(12 Pt 1), 1264–1272. doi: 10.1002/lary.5541011202

CrossRef Full Text | Google Scholar

Cabral Junior, F., Pinna, M. H., Alves, R. D., Malerbi, A. F., and Bento, R. F. (2016). Cochlear implantation and single-sided deafness: a systematic review of the literature. Int. Arch. Otorhinolaryngol. 20, 69–75. doi: 10.1055/s-0035-1559586

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, J. L., Pross, S. E., Findlay, A. M., Mizuiri, D., Henderson-Sabes, J., Garrett, C., et al. (2016). Spatial plasticity of the auditory cortex in single-sided deafness. Laryngoscope 126, 2785–2791. doi: 10.1002/lary.25961

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheung, S. W., Atencio, C. A., Levy, E. R. J., Froemke, R. C., and Schreiner, C. E. (2017). Anisomorphic cortical reorganization in asymmetric sensorineural hearing loss. J. Neurophysiol. 118, 932–948. doi: 10.1152/jn.00119.2017

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheung, S. W., Bonham, B. H., Schreiner, C. E., Godey, B., and Copenhaver, D. A. (2009). Realignment of interaural cortical maps in asymmetric hearing loss. J. Neurosci. 29, 7065–7078. doi: 10.1523/JNEUROSCI.6072-08.2009

CrossRef Full Text | Google Scholar

Dalal, S. S., Guggisberg, A. G., Edwards, E., Sekihara, K., Findlay, A. M., Canolty, R. T., et al. (2008). Five-dimensional neuroimaging: localization of the time-frequency dynamics of cortical activity. Neuroimage 40, 1686–1700. doi: 10.1016/j.neuroimage.2008.01.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Dalal, S. S., Zumer, J. M., Guggisberg, A. G., Trumpis, M., Wong, D. D., Sekihara, K., et al. (2011). MEG/EEG source reconstruction, statistical evaluation, and visualization with NUTMEG. Comput. Intell. Neurosci. 2011:758973. doi: 10.1155/2011/758973

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, H., Qin, W., Liang, M., Ming, D., Wan, B., Li, Q., et al. (2015). Cross-modal activation of auditory regions during visuo-spatial working memory in early deafness. Brain 138(Pt 9), 2750–2765. doi: 10.1093/brain/awv165

PubMed Abstract | CrossRef Full Text | Google Scholar

Douglas, S. A., Yeung, P., Daudia, A., Gatehouse, S., and O’Donoghue, G. M. (2007). Spatial hearing disability after acoustic neuroma removal. Laryngoscope 117, 1648–1651. doi: 10.1097/MLG.0b013e3180caa162

PubMed Abstract | CrossRef Full Text | Google Scholar

Dye, M. W., Hauser, P. C., and Bavelier, D. (2009). Is visual selective attention in deaf individuals enhanced or deficient? The case of the useful field of view. PLoS One 4:e5640. doi: 10.1371/journal.pone.0005640

PubMed Abstract | CrossRef Full Text | Google Scholar

Finney, E. M., Clementz, B. A., Hickok, G., and Dobkins, K. R. (2003). Visual stimuli activate auditory cortex in deaf subjects: evidence from MEG. Neuroreport 14, 1425–1427. doi: 10.1097/01.wnr.0000079894.11980.6a

PubMed Abstract | CrossRef Full Text | Google Scholar

Finney, E. M., Fine, I., and Dobkins, K. R. (2001). Visual stimuli activate auditory cortex in the deaf. Nat. Neurosci. 4, 1171–1173. doi: 10.1038/nn763

PubMed Abstract | CrossRef Full Text | Google Scholar

Gatehouse, S., and Noble, W. (2004). The speech, spatial and qualities of hearing scale (SSQ). Int. J. Audiol. 43, 85–99. doi: 10.1080/14992020400050014

CrossRef Full Text | Google Scholar

Herman, A. B., Houde, J. F., Vinogradov, S., and Nagarajan, S. S. (2013). Parsing the phonological loop: activation timing in the dorsal speech stream determines accuracy in speech reproduction. J. Neurosci. 33, 5439–5453. doi: 10.1523/JNEUROSCI.1472-12.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Hickok, G., and Poeppel, D. (2004). Dorsal and ventral streams: a framework for understanding aspects of the functional anatomy of language. Cognition 92, 67–99. doi: 10.1016/j.cognition.2003.10.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Hinkley, L. B., Mizuiri, D., Hong, O., Nagarajan, S. S., and Cheung, S. W. (2015). Increased striatal functional connectivity with auditory cortex in tinnitus. Front. Hum. Neurosci. 9:568. doi: 10.3389/fnhum.2015.00568

PubMed Abstract | CrossRef Full Text | Google Scholar

Hribar, M., Suput, D., Carvalho, A. A., Battelino, S., and Vovk, A. (2014). Structural alterations of brain grey and white matter in early deaf adults. Hear. Res. 318, 1–10. doi: 10.1016/j.heares.2014.09.008

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Jerbi, K., Vidal, J. R., Ossandon, T., Dalal, S. S., Jung, J., Hoffmann, D., et al. (2010). Exploring the electrophysiological correlates of the default-mode network with intracerebral EEG. Front. Syst. Neurosci. 4:27. doi: 10.3389/fnsys.2010.00027

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, G., Ju, H. M., Lee, S. H., Kim, H. S., Kwon, J. A., and Seo, Y. J. (2017). Efficacy of bone-anchored hearing aids in single-sided deafness: a systematic review. Otol. Neurotol. 38, 473–483. doi: 10.1097/MAO.0000000000001359

PubMed Abstract | CrossRef Full Text | Google Scholar

Lieu, J. E., Tye-Murray, N., and Fu, Q. (2012). Longitudinal study of children with unilateral hearing loss. Laryngoscope 122, 2088–2095. doi: 10.1002/lary.23454

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, Y., Wang, J., Wu, C., Wai, Y., Yu, J., and Ng, S. (2008). Diffusion tensor imaging of the auditory pathway in sensorineural hearing loss: changes in radial diffusivity and diffusion anisotropy. J. Magn. Reson. Imaging 28, 598–603. doi: 10.1002/jmri.21464

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, B., Feng, Y., Yang, M., Chen, J. Y., Li, J., Huang, Z. C., et al. (2015). Functional connectivity in patients with sensorineural hearing loss using resting-state MRI. Am. J. Audiol. 24, 145–152. doi: 10.1044/2015_AJA-13-0068

PubMed Abstract | CrossRef Full Text | Google Scholar

Lomber, S. G., Meredith, M. A., and Kral, A. (2010). Cross-modal plasticity in specific auditory cortices underlies visual compensations in the deaf. Nat. Neurosci. 13, 1421–1427. doi: 10.1038/nn.2653

PubMed Abstract | CrossRef Full Text | Google Scholar

Meredith, M. A., Kryklywy, J., McMillan, A. J., Malhotra, S., Lum-Tai, R., and Lomber, S. G. (2011). Crossmodal reorganization in the early deaf switches sensory, but not behavioral roles of auditory cortex. Proc. Natl. Acad. Sci. U.S.A. 108, 8856–8861. doi: 10.1073/pnas.1018519108

PubMed Abstract | CrossRef Full Text | Google Scholar

Miao, W., Li, J., Tang, M., Xian, J., Li, W., Liu, Z., et al. (2013). Altered white matter integrity in adolescents with prelingual deafness: a high-resolution tract-based spatial statistics imaging study. AJNR Am. J. Neuroradiol. 34, 1264–1270. doi: 10.3174/ajnr.A3370

PubMed Abstract | CrossRef Full Text | Google Scholar

Noble, W., and Gatehouse, S. (2004). Interaural asymmetry of hearing loss, speech, spatial and qualities of hearing scale (SSQ) disabilities, and handicap. Int. J. Audiol. 43, 100–114. doi: 10.1080/14992020400050015

PubMed Abstract | CrossRef Full Text | Google Scholar

Ossandon, T., Jerbi, K., Vidal, J. R., Bayle, D. J., Henaff, M. A., Jung, J., et al. (2011). Transient suppression of broadband gamma power in the default-mode network is correlated with task complexity and subject performance. J. Neurosci. 31, 14521–14530. doi: 10.1523/JNEUROSCI.2483-11.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

Ossandon, T., Vidal, J. R., Ciumas, C., Jerbi, K., Hamame, C. M., Dalal, S. S., et al. (2012). Efficient ”pop-out” visual search elicits sustained broadband gamma activity in the dorsal attention network. J. Neurosci. 32, 3414–3421. doi: 10.1523/JNEUROSCI.6048-11.2012

CrossRef Full Text | Google Scholar

Peters, J. P., Smit, A. L., Stegeman, I., and Grolman, W. (2015). Review: bone conduction devices and contralateral routing of sound systems in single-sided deafness. Laryngoscope 125, 218–226. doi: 10.1002/lary.24865

PubMed Abstract | CrossRef Full Text | Google Scholar

Pross, S. E., Chang, J. L., Mizuiri, D., Findlay, A. M., Nagarajan, S. S., and Cheung, S. W. (2015). Temporal cortical plasticity in single-sided deafness: a functional imaging study. Otol. Neurotol. 36, 1443–1449. doi: 10.1097/MAO.0000000000000821

PubMed Abstract | CrossRef Full Text | Google Scholar

Rachakonda, T., Shimony, J. S., Coalson, R. S., and Lieu, J. E. (2014). Diffusion tensor imaging in children with unilateral hearing loss: a pilot study. Front. Syst. Neurosci. 8:87. doi: 10.3389/fnsys.2014.00087

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramot, M., Fisch, L., Harel, M., Kipervasser, S., Andelman, F., Neufeld, M. Y., et al. (2012). A widely distributed spectral signature of task-negative electrocorticography responses revealed during a visuomotor task in the human cortex. J. Neurosci. 32, 10458–10469. doi: 10.1523/JNEUROSCI.0877-12.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Rauschecker, J. P. (2012). Ventral and dorsal streams in the evolution of speech and language. Front. Evol. Neurosci. 4:7. doi: 10.3389/fnevo.2012.00007

PubMed Abstract | CrossRef Full Text | Google Scholar

Rauschecker, J. P. (2018). Where, when, and how: are they all sensorimotor? towards a unified view of the dorsal pathway in vision and audition. Cortex 98, 262–268. doi: 10.1016/j.cortex.2017.10.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Rauschecker, J. P., and Scott, S. K. (2009). Maps and streams in the auditory cortex: nonhuman primates illuminate human speech processing. Nat. Neurosci. 12, 718–724. doi: 10.1038/nn.2331

PubMed Abstract | CrossRef Full Text | Google Scholar

Rauschecker, J. P., and Tian, B. (2000). Mechanisms and streams for processing of ”what” and ”where” in auditory cortex. Proc. Natl. Acad. Sci. U.S.A. 97, 11800–11806. doi: 10.1073/pnas.97.22.11800

PubMed Abstract | CrossRef Full Text | Google Scholar

Rigby, P. L., Shah, S. B., Jackler, R. K., Chung, J. H., and Cooke, D. D. (1997). Acoustic neuroma surgery: outcome analysis of patient-perceived disability. Am. J. Otol. 18, 427–435.

PubMed Abstract | Google Scholar

Ryu, C. W., Park, M. S., Byun, J. Y., Jahng, G. H., and Park, S. (2015). White matter integrity associated with clinical symptoms in tinnitus patients: a tract-based spatial statistics study. Eur. Radiol. 26, 2223–2232. doi: 10.1007/s00330-015-4034-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt, S. A., Akrofi, K., Carpenter-Thompson, J. R., and Husain, F. T. (2013). Default mode, dorsal attention and auditory resting state networks exhibit differential functional connectivity in tinnitus and hearing loss. PLoS One 8:e76488. doi: 10.1371/journal.pone.0076488

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmithorst, V. J., Plante, E., and Holland, S. (2014). Unilateral deafness in children affects development of multi-modal modulation and default mode networks. Front. Hum. Neurosci. 8:164. doi: 10.3389/fnhum.2014.00164

PubMed Abstract | CrossRef Full Text | Google Scholar

Sekihara, K., Kawabata, Y., Ushio, S., Sumiya, S., Kawabata, S., Adachi, Y., et al. (2016). Dual signal subspace projection (DSSP): a novel algorithm for removing large interference in biomagnetic measurements. J. Neural. Eng. 13:036007. doi: 10.1088/1741-2560/13/3/036007

PubMed Abstract | CrossRef Full Text | Google Scholar

Seydell-Greenwald, A., Raven, E. P., Leaver, A. M., Turesky, T. K., and Rauschecker, J. P. (2014). Diffusion imaging of auditory and auditory-limbic connectivity in tinnitus: preliminary evidence and methodological challenges. Neural. Plast. 2014:145943. doi: 10.1155/2014/145943

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, S. M. (2002). Fast robust automated brain extraction. Hum. Brain Mapp. 17, 143–155. doi: 10.1002/hbm.10062

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, S. M., Jenkinson, M., Woolrich, M. W., Beckmann, C. F., Behrens, T. E., Johansen-Berg, H., et al. (2004). Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage 23(Suppl. 1), S208–S219. doi: 10.1016/j.neuroimage.2004.07.051

PubMed Abstract | CrossRef Full Text | Google Scholar

Tufarelli, D., Meli, A., Alesii, A., De Angelis, E., Badaracco, C., Falcioni, M., et al. (2006). Quality of life after acoustic neuroma surgery. Otol. Neurotol. 27, 403–409. doi: 10.1097/00129492-200604000-00018

CrossRef Full Text | Google Scholar

van Zon, A., Peters, J. P., Stegeman, I., Smit, A. L., and Grolman, W. (2015). Cochlear implantation for patients with single-sided deafness or asymmetrical hearing loss: a systematic review of the evidence. Otol. Neurotol. 36, 209–219. doi: 10.1097/MAO.0000000000000681

PubMed Abstract | CrossRef Full Text | Google Scholar

Vila, P. M., and Lieu, J. E. (2015). Asymmetric and unilateral hearing loss in children. Cell Tissue Res. 361, 271–278. doi: 10.1007/s00441-015-2208-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Vlastarakos, P. V., Nazos, K., Tavoulari, E. F., and Nikolopoulos, T. P. (2014). Cochlear implantation for single-sided deafness: the outcomes. An evidence-based approach. Eur. Arch. Otorhinolaryngol. 271, 2119–2126. doi: 10.1007/s00405-013-2746-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Vos, S. B., Haakma, W., Versnel, H., Froeling, M., Speleman, L., Dik, P., et al. (2015). Diffusion tensor imaging of the auditory nerve in patients with long-term single-sided deafness. Hear. Res. 323, 1–8. doi: 10.1016/j.heares.2015.01.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, C. M., Ng, S. H., and Liu, T. C. (2009). Diffusion tensor imaging of the subcortical auditory tract in subjects with long-term unilateral sensorineural hearing loss. Audiol. Neurootol. 14, 248–253. doi: 10.1159/000191282

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, G. Y., Yang, M., Liu, B., Huang, Z. C., Li, J., Chen, J. Y., et al. (2016). Changes of the directional brain networks related with brain plasticity in patients with long-term unilateral sensorineural hearing loss. Neuroscience 313, 149–161. doi: 10.1016/j.neuroscience.2015.11.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: auditory memory, diffusion tensor imaging, magnetoencephalography, plasticity, single-sided deafness

Citation: Shang Y, Hinkley LB, Cai C, Subramaniam K, Chang Y-S, Owen JP, Garrett C, Mizuiri D, Mukherjee P, Nagarajan SS and Cheung SW (2018) Functional and Structural Brain Plasticity in Adult Onset Single-Sided Deafness. Front. Hum. Neurosci. 12:474. doi: 10.3389/fnhum.2018.00474

Received: 26 June 2018; Accepted: 12 November 2018;
Published: 27 November 2018.

Edited by:

Hidehiko Okamoto, International University of Health and Welfare, Japan

Reviewed by:

Shinichi Furuya, Sony Computer Science Laboratories, Japan
Dona M. P. Jayakody, Ear Science Institute Australia, Australia

Copyright © 2018 Shang, Hinkley, Cai, Subramaniam, Chang, Owen, Garrett, Mizuiri, Mukherjee, Nagarajan and Cheung. 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) and the copyright owner(s) 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.

*Correspondence: Steven W. Cheung, steven.cheung@ucsf.edu

Download