Abstract

Although evidence from studies on blind adults indicates that visual deprivation early in life leads to structural and functional disruption and reorganization of the brain, whether young blind people show similar patterns remains unknown. Therefore, this study is aimed at exploring the structural and functional alterations of the brain of early-blind adolescents (EBAs) compared to normal-sighted controls (NSCs) and investigating the effects of residual light perception on brain microstructure and function in EBAs. We obtained magnetic resonance imaging (MRI) data from 23 EBAs (8 with residual light perception (LPs), 15 without light perception (NLPs)) and 21 NSCs (age range 11-19 years old). Whole-brain voxel-based analyses of diffusion tensor imaging metrics and region-of-interest analyses of resting-state functional connectivity (RSFC) were performed to compare patterns of brain microstructure and the corresponding RSFC between the groups. The results showed that structural disruptions of LPs and NLPs were mainly located in the occipital visual pathway. Compared with NLPs, LPs showed increased fractional anisotropy (FA) in the superior frontal gyrus and reduced diffusivity in the caudate nucleus. Moreover, the correlations between FA of the occipital cortices or mean diffusivity of the lingual gyrus and age were consistent with the development trajectory of the brain in NSCs, but inconsistent or even opposite in EBAs. Additionally, we found functional, but not structural, reorganization in NLPs compared with NSCs, suggesting that functional neuroplasticity occurs earlier than structural neuroplasticity in EBAs. Altogether, these findings provided new insights into the mechanisms underlying the neural reorganization of the brain in adolescents with early visual deprivation.

1. Introduction

It is well established that early-blind adults, i.e., those individuals who lose sight at birth or within a short period after birth, show alterations in the neural structure and function due to the absence of visual inputs. Structurally, the morphological changes of the grey matter (GM) in early-blind adults have been studied extensively and include decreased GM volume in the primary visual area [15]; increased GM volume in the occipital, frontal, and entorhinal cortices [4] and the globus pallidus [6]; and increased cortical thickness of the early visual areas [2, 7, 8]. In addition to GM alterations, white matter (WM) alterations were also identified based on the diffusion tensor imaging (DTI) datasets of early-blind adults. In detail, microstructural alterations were mainly found in the optic radiation and the lateral geniculate nucleus with decreased WM volume [2, 9, 10] or density [9, 11]. Other microstructural disruptions due to early visual deprivation were also found in the inferior longitudinal fasciculus [12] and the occipital/temporal/parietal thalamocortical projections [4, 13]. In addition, Leporé et al. [3], using tensor-based morphometry, reported interesting hypertrophy in the prefrontal and parietal WM and in a section of the splenium of the corpus callosum in early-blind adults. Interestingly, some studies [11, 14] also identified increased fractional anisotropy in the corticospinal tract in early-blind adults by tract-based quantitative analysis, suggesting structural neuroplasticity. To some extent, such morphological changes reflect deprivation-induced processes, both neurodegenerative and neurodevelopmental, following the loss of visual input and/or altered visual experience [3, 5].

In addition to structural alterations, functional changes were also identified in early-blind adults in many studies. Task-based research using functional magnetic resonance imaging (fMRI) demonstrated that the occipital cortex could process nonvisual information, such as tactile [15, 16] and auditory input [1620], and some complex cognitive- and verbal-related input [17, 2123]. These results collectively suggested that the visual cortex of blind individuals exhibits cross-modal reorganization properties. Moreover, activity in the primary, secondary, and tertiary visual cortices of blind adults was increased, spatially and in amplitude, after sensory substitution training, even in short-term (10 min) [24]. Using resting-state fMRI (rs-fMRI), both decreased functional connectivity [7, 22, 2528] and increased functional connectivity of the occipital cortex [25, 29, 30] were identified in early-blind adults compared with sighted controls, providing extra evidence of functional neuroplasticity in early-blind adults.

However, these results about structural and functional disruptions and neuroplasticity were derived from early-blind adults. Recently, we investigated resting-state interhemispheric functional connectivity in early-blind adolescents (EBAs) and found significantly lower voxel-mirrored homotopic connectivity (VMHC) in the primary visual cortex, visual association cortex, and somatosensory association cortex [31]. Although the study focused on EBAs, it is still unclear if this group would show structural and functional alteration patterns similar to those of early-blind adults. Moreover, most of the previous studies on early blindness included subjects with residual light perception, thus failing to exclude possible effects of residual light perception on the structure and function of the blind brain. Therefore, we aimed to (1) investigate whether and how the brain’s structure and function are altered in EBAs and (2) investigate the effects of residual light perception on brain structure and function. Given the evidence just discussed, we hypothesized that structural and functional alterations in EBAs would mostly be identified in the brain regions associated with vision. To verify this hypothesis, we recruited a cohort of 23 EBAs (8 with residual light perception (LPs), 15 without light perception (NLPs)) and 21 NSCs (age range: 11–19 years old) to explore the structural and functional brain changes in EBAs using the whole-brain voxel-based analysis (VBA) and the resting-state functional connectivity (RSFC).

2. Materials and Methods

2.1. Subjects

Twenty-three EBAs with onset year were enrolled from the Guangzhou City Blind School. Among them, 8 EBAs with residual light perception were included in the LP group (6 males and 2 females) and 15 EBAs without light perception were included in the NLP group (8 males and 7 females). Twenty-one age- and sex-matched normal-sighted volunteers were recruited as the NSC group (10 males and 11 females). The demographics of the three groups are summarized in Table 1. All subjects met the following inclusion criteria: (1) right-handed and (2) age ranging from 11 to 19 years. The exclusion criteria were (1) any history of psychiatric or neurologic diseases, (2) symptomatic or atypical neuralgia, and (3) identifiable MRI abnormalities, such as demyelination, vascular malformations, or tumors. All individuals and their guardians signed a written informed consent form prior to the MRI examinations. All research procedures were approved by the Ethics Committee of the First Affiliated Hospital of Jinan University.

2.2. Data Acquisition

A 3.0T MRI scanner (Discovery MR750 System; General Electric, Milwaukee, WI, USA) with an 8-channel head coil was used. The subjects were scanned in the supine position and were placed head first into the scanner. The position of the head was fixed using several foam cushions. For each subject, the DTI and rs-fMRI data were acquired. For DTI, we acquired 75 diffusion-weighted images (/mm2) and 5 nondiffusion-weighted images (/mm2) using a spin-echo echo-planar imaging sequence with the following parameters: number of , repetition time , echo time , acquisition , and voxel . The array spatial sensitivity encoding technique was used with an acceleration factor of 2 to reduce acquisition time and anamorphosis. The rs-fMRI data were acquired using an echo-planar imaging pulse sequence, sensitive to blood-oxygen-level-dependent contrast (, , flip , matrix , field of , slice , , 240 time points, sequence , and 41 slices in the axial plane). During rs-fMRI, subjects were instructed to stay awake, not to move, to think of nothing in particular, to relax, and to keep their eyes closed.

2.3. VBA and Statistics Analysis

We performed DTI preprocessing using the FMRIB Software Library (FSL5.0; http://www.fmrib.ox.ac.uk/fsl), including eddy-current and motion correction, brain mask creation for the DTI data, exclusion of voxels outside the braincase, and diffusion tensor reconstruction with the “DTIFIT” toolbox [10]. All diffusion-related parameters including fractional anisotropy (FA), mean diffusivity (MD), radial diffusivity (RD), and axial diffusivity (AD) were calculated. Following this, individual FA maps were analyzed using the statistical parametric mapping software (SPM8; Wellcome Department of Cognitive Neurology, London, UK; http://www.fil.ion.ucl.ac.uk/spm). The FA maps were normalized to the standard Montreal Neurological Institute (MNI) template and smoothed with a 6 mm full-width at half-maximum (FWHM) Gaussian kernel to reduce the effects of misregistration in spatial normalization [32]. Two-sample -tests were conducted to assess intergroup differences in FA between the groups with age and sex as covariates. In order to account for multiple comparisons, the AlphaSim method was used with a voxel-level threshold of and a cluster-level threshold of . VBA of the MD, RD, and AD maps was performed similarly.

The clusters showing significant differences in the DTI parameters between the groups were identified according to the standard brain atlas, and the DTI values from these significant clusters of the corresponding tensor parameter maps (FA, MD, RD, and AD) in the subjects were extracted. Subsequently, correlation analysis followed by an independent two-sample -test was performed using the Statistical Package for the Social Sciences (SPSS) software (version 23.0; IBM Corporation, NY, USA). For correlations between DTI metric values and age, the Pearson correlation coefficient was computed separately for each group. Two-sample -tests were used to detect differences in the DTI values between the groups.

2.4. Resting-State fMRI Data Preprocessing

The rs-fMRI data were preprocessed using the toolbox for Data Processing & Analysis for Brain Imaging (DPABI, http://rfmri.org/dpabi) [33]. The preprocessing steps were (1) discarding the first 10 time points to avoid transient signal changes that occurred before magnetization reached the steady state and subjects’ adaptation to the scanning noise; (2) slice timing correction; (3) head motion realignment; (4) normalization to the standard MNI template in the DPABI package; (5) smoothing with a 6 mm Gaussian kernel; (6) linear regression to remove confounding factors, including six head motion parameters, as well as the ventricle, WM, and global signals; (7) scrubbing with interpolation to remove volumes with head motion spatial shift larger than 1.5 mm and rotation greater than 1.5° in any direction; and (8) band-pass filtering from 0.01 to 0.1 Hz. The resulting images were manually inspected by two experienced researchers for image quality.

2.5. RSFC and Statistics Analysis

To investigate the corresponding functional alterations, the clusters showing significant group difference in VBA were selected as the seed ROIs for RSFC analysis. Firstly, the peak coordinates of each significant cluster obtained from VBA were determined and used to create a ROI with a radius of 6 mm. Secondly, the RSFC was examined with a seed-voxel correlation approach, in which the correlations of the mean time-course signal in a seed region with all other voxels in the whole brain were calculated. Before statistical analysis, the correlation coefficients were transformed into values using the Fisher -to- transformation to improve the normality. Two-sample -tests between the groups were then performed using the DPABI toolbox, with age, sex, and head motion parameters as covariates. The statistical results were corrected for multiple comparisons using the AlphaSim method, with a voxel-level threshold of and a cluster-level threshold of . Finally, the results were projected onto a smoothed International Consortium for Brain Mapping (ICBM152) surface template using the BrainNet Viewer (http://www.nitrc.org/projects/bnv/) [34].

3. Results

3.1. Demographic Characteristics

Six male and two female subjects were included in the LP group and eight male and seven female subjects in the NLP group; we also recruited 21 age- and sex-matched NSCs. Analysis of variance (ANOVA) followed by two-sample -tests and chi-square tests were performed to detect age and sex differences, respectively, between the various groups (Table 2). No significant differences in the demographic characteristics were identified between the three groups (age: , comparing the LP and NLP groups; , , LP vs. NSC; , , NLP vs. NSC; sex: , LP vs. NLP; , LP vs. NSC; and , NLP vs. NSC).

3.2. VBA Results

The two-sample t-test revealed a reduced FA in the left occipital lobe/subgyral, right parahippocampal gyrus, and right superior occipital gyrus in NLPs compared to NSCs. We also found higher MD in the bilateral lingual gyrus, as well as higher RD in the right lingual gyrus and right parahippocampal gyrus in NLPs, compared with NSCs (, AlphaSim corrected; Table 3 and Figure 1). No significant AD alterations were observed between these two groups. Compared with NSCs, LPs showed significant FA reduction in the bilateral middle occipital gyrus, left lingual gyrus, left cuneus, left parahippocampal gyrus, and right parietal lobe/subgyral (, AlphaSim corrected; Table 3 and Figure 2). No significant differences in MD, RD, and AD were found when comparing the LP and NSC groups. No elevated FA or reduced diffusivity parameters were found for any of the clusters in NLPs or LPs compared with NSCs, even with a loose threshold (, AlphaSim corrected). Unfortunately, the comparison of all diffusion tensor maps between the LP and NLP groups showed no results at the voxel level, at AlphaSim-corrected . However, some voxels in the brain of the LP group exhibited trends (at a voxel-level AlphaSim-corrected ) towards increased/decreased DTI parameters compared with NLPs, including reduced FA in the right fusiform gyrus, right middle temporal gyrus, and superior occipital gyrus; increased FA in the right superior frontal gyrus; increased MD, RD, and AD in the right fusiform gyrus and middle temporal gyrus; and reduced MD, RD, and AD in the bilateral caudate (Table SI and Figure SI in the Supplementary Materials).

3.3. Correlation Analysis Results

Correlation analysis showed positive and significant correlations between the MD and RD values of the right lingual gyrus and age in NLPs (MD: , ; RD: , ), but negative and significant correlations in NSCs (MD: , ; RD: , ). In NSCs, the FA values of the left occipital lobe/subgyral (, ) and the right superior occipital gyrus (, ) were positively correlated with age, while the MD values of the left lingual gyrus were negatively correlated with age (, ) (Figure 3).

3.4. Head Motion Effects

Head motion time courses were computed by estimating the shift in each direction and the rotations about each axis for each of the 230 consecutive volumes. To assess the head motion confounders, we calculated the mean framewise displacement among the three groups (LP: ; NLP: ; and NSC: ). Moreover, ANOVA was used to identify potential differences in head motion between the 3 groups, but none was found. (ANOVA: , ).

3.5. RSFC Results

The whole-brain RSFC of the brain regions with significant microstructural alterations was analyzed and compared between the NLP and NSC groups. The left occipital lobe/subgyral of NLPs showed enhanced RSFC with the left inferior frontal triangular areas, left middle frontal gyrus, and right inferior frontal opercular areas; the right parahippocampal gyrus of NLPs showed enhanced RSFC with the left superior frontal gyrus; the right superior occipital gyrus of NLPs showed enhanced RSFC with the right supramarginal gyrus/left superior frontal gyrus and decreased RSFC with the left fusiform gyrus (, AlphaSim corrected; Table 4 and Figure 4). When comparing LPs with NSCs, only the left middle occipital gyrus showed a significantly decreased RSFC in LPs with the right middle occipital gyrus, when setting the bilateral middle occipital gyrus, left parahippocampal gyrus, left lingual gyrus, left cuneus, and right parietal lobe/subgyral as ROI seeds (, AlphaSim corrected; Table 4).

To investigate the effects of residual light perception on the brain function in blind brains, the differences in the whole-brain RSFC of the brain regions showing microstructural alterations between the LP and NLP groups were analyzed and compared. A trend of RSFC change in LPs relative to NLPs could also be observed and specifically enhanced RSFC between the right fusiform gyrus and the right superior frontal orbital areas/right middle temporal gyrus and between the right superior frontal gyrus and the left inferior temporal gyrus. On the contrary, the RSFC between the right superior frontal gyrus and the right postcentral gyrus/superior frontal gyrus was decreased in the LP compared to the NLP group (, AlphaSim-corrected; Table SII and Figure SII in the Supplementary Materials).

4. Discussion

We investigated the alterations of brain microstructure and their corresponding functional connectivity in EBAs with and without light perception using the DTI and rs-fMRI datasets. Compared to NSCs, both EBA groups showed significant brain microstructural disruptions mainly in the posterior visual pathway of the occipital lobe and the extrastriate visual cortex. Moreover, the correlations between bilateral occipital cortex FA or bilateral lingual gyrus MD and age were consistent with the developmental trajectory of the brain in NSCs, but inconsistent or even opposite in EBAs. Additionally, using brain regions with significant DTI parameter changes in our study as the ROIs, we identified enhanced RSFC in the occipital cortex with the bilateral frontal cortex and the homolateral parietal cortex and reduced intrahemispheric RSFC in the right visual cortex in NLPs, as well as reduced interhemispheric RSFC within the middle occipital cortex in LPs compared with NSCs. Our findings provide additional evidence that early visual deprivation may lead to functional neuroplasticity earlier than structural neuroplasticity in EBAs.

4.1. Disruptions in the Brain Microstructure of EBAs

The FA reduction in the bilateral visual pathway of the occipital cortex in both EBA groups, elevated MD of the bilateral lingual gyrus, and elevated RD of the right lingual gyrus in NLPs vs. NSCs were consistent with the findings of previous research on young or middle-aged blind adults [2, 9, 12, 35, 36]. These structural disruptions were also supported by other evidence of reduced anatomical connectivity associated with the visual cortex [12, 37] and decreased GM and WM volume in the occipital cortex of blind individuals [6, 38]. Our study also validated the GM and WM changes in the visual cortex of EBAs at the microstructural level, which can be explained by transneuronal degeneration and/or immaturity due to early visual deprivation [9, 39]. Moreover, we found diminished FA in the left parahippocampal gyrus of LPs and the right parahippocampal gyrus of NLPs compared to NSCs. The right parahippocampal gyrus is important for learning and visuospatial configuration of objects [40], while the left parahippocampal gyrus is involved in cumulative verbal memory [41]. Considering the asymmetry function of the bilateral parahippocampal gyri during learning, their inconsistent alteration in the two groups may provide evidence of different learning patterns between the LP and NLP groups.

Furthermore, to investigate the effects of residual light perception on brain microstructure, we compared the DTI parameters between LPs and NLPs. The results showed increased FA in the right superior frontal gyrus and reduced diffusivity parameters in the bilateral caudate nucleus in LPs compared with NLPs. Since the superior frontal gyrus is anatomically connected with the cognitive control network and functionally involved in complex cognitive processing [42, 43], the increased FA of the right superior frontal gyrus in LPs vs. NLPs may suggest that residual light perception could play a role in complex cognitive function or reduce the damage to cognitive function induced in EBAs by the lack of visual information. As the caudate nucleus is known to integrate spatial information with motor behavior and is associated with some learning processing [4446], its decreased diffusivity parameters in LPs vs. NLPs demonstrate that the function of guiding motor performance and learning in the caudate nucleus may be partly mediated by light perception.

We also calculated the correlations between age and DTI metrics of the altered brain regions in the three groups. The FA of the bilateral occipital cortex was positively correlated with age in NSCs but not in NLPs, and the MD and RD of the right lingual gyrus were negatively correlated with age in NSCs, but positively in NLPs. The trends of DTI parameters alteration with age in NSCs were consistent with the developmental trajectory of WM in the young [4749], while these trends were not observed or were reversed in NLPs. This inconsistency suggests that microstructural disruptions caused by early visual deprivation in the corresponding brain regions may progressively accumulate in adolescents. Given the continuous development of WM in adolescence, a probably more reasonable explanation for the trend shown by the NLPs is the coexistence of blindness-induced damage and normal cerebral development.

4.2. Alterations in Functional Connectivity

Changes in brain microstructure preferentially point to abnormalities in brain functional connectivity [50, 51]. Resting-state fMRI can reflect the integrated features of intrinsic functional connectivity in the brains of visually impaired persons [27]. Many previous studies reported changes in RSFC between the occipital cortex and other sensory [7, 2527, 52], somatosensory [25, 26, 28, 52, 53], and motor cortices [25, 26, 28] in early-blind adults. In the current study, the occipital gyrus in NLPs showed enhanced RSFC with some frontal areas, including the left inferior frontal triangle, left middle frontal gyrus, left superior frontal gyrus, and right inferior frontal opercular areas. Such findings are consistent with those of previous studies [54, 55]. In healthy humans, these frontal regions are classically considered as high-level areas related to language and cognition and involved in verbal fluency [56, 57], single word and motor speech production [58, 59], and error detection and imitation [60, 61]. Beyond these classical functions, they also participate in verbal working memory [62, 63]. In blind individuals, the occipital visual cortex and the frontal language and cognitive areas were activated when performing Braille reading tasks [64, 65]. These findings indicated that these frontal areas might play an important role, related to the language network, in the blind brain. Additionally, compared with sighted controls, early-blind participants had more extensive activity in the occipital visual areas and frontal language areas when performing a vibrotactile memory task [64, 65]. All evidence may explain why the RSFC between the occipital areas and frontal areas was increased in the NLPs. From another perspective, the elevated RSFC between the occipital visual areas and the frontal language regions can provide compelling evidence for the existence of functional neuroplasticity in NLPs.

We also observed reduced interhemispheric RSFC between the bilateral middle occipital gyrus in LPs and between the right superior occipital gyrus and the left fusiform gyrus in NLPs compared with NSCs, extending the findings of many previous studies [22, 25, 31, 66]. According to some published reports, such reduced interhemispheric communication may be interpreted as extensive abnormalities in WM integrity, particularly in the corpus callosum [67, 68]. Many existing studies found WM abnormalities in the splenium of the corpus callosum, thought to play a crucial role in the communication between the two cerebral hemispheres, in early-blind young, middle-aged, or elderly adults [3, 10, 35, 69, 70]. However, no WM structural abnormalities were found in the splenium of the corpus callosum in the present study. One possible reason may be that the posterior part of the corpus callosum continues to develop during adolescence and increases in density [47, 48, 71, 72], likely counteracting the neurodegenerative effects of visual deprivation [73]. Based on prior and present findings, we speculated that the RSFC changes in EBAs might result not only from a combination mechanism of general loss and complementary plasticity [25] but also from an antagonistic mechanism of neural development. Thus, our findings provide new information regarding the functional profiles that emerge in the early life of blind adolescents.

Furthermore, we evaluated the RSFC differences between the LP and NLP groups to investigate the differences in functional connectivity patterns depending on residual light perception. Compared with NLPs, LPs showed enhanced RSFC in the right fusiform with the homolateral middle temporal gyrus and reduced RSFC in the right superior frontal gyrus with the right postcentral gyrus. In sighted subjects, the fusiform and middle temporal gyri have been reported to be involved in the ventral and dorsal visual streams, respectively, and to be robustly activated during object shape/location detection and visual-motor tasks [74, 75]. There is no doubt that light perception is crucial for such processing. Therefore, the residual perception of light in LPs may account for the enhanced RSFC between these regions. Given that the postcentral gyrus is well-known as the primary somatosensory area and is associated with fine touch sense, the decreased RSFC of this area in LPs may suggest that the touch sense of early-blind individuals with residual light perception is not as strong as that of those without light perception.

4.3. Limitation

Our findings of the effects of residual light perception on brain microstructure and function in EBAs were derived at a loose level of statistical significance, presumably due to the small sample size. Therefore, these results should be interpreted with caution, and further studies with larger sample size will be needed to provide solid evidence of these phenomena.

5. Conclusion

The present study demonstrated significant microstructural and functional alterations in EBAs with and without residual light perception when compared to NSCs. Our findings provide additional evidence that early visual deprivation may lead to functional neuroplasticity earlier than structural neuroplasticity in EBAs. The structural alterations in EBAs contain complex mechanisms including general loss, complementary plasticity, and neural development. Collectively, these results provided new insights into the mechanisms underlying the reorganization of the brain in adolescents with early visual deprivation.

Data Availability

The datasets generated and/or analyzed during the current study are available from the corresponding author Dr. Hengguo Li ([email protected]) on reasonable request.

Conflicts of Interest

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.

Authors’ Contributions

ZZ, XL, FH, JZ, and JL acquired and analyzed data. JX, XL, and HL conceived this study and designed the experiments. ZZ wrote the article with the help of JX, LS, and QH. All authors were involved in data interpretation and critically revising the manuscript.

Acknowledgments

This work was supported by the Shenzhen Sanming Project (SZSM201512038), the Science and Technology Innovation Program of Shenzhen (JCYJ20160531184531506), and the Natural Science Foundation of Guangdong Province (Nos. 2017A030313744 and 2016A030313180).

Supplementary Materials

Table SI: brain regions with significant changes in diffusion tensor imaging metrics between LPs and NLPs. Table SII: brain regions with significant changes of RSFC between the LP and NLP groups. Figure SI: group differences of fractional anisotropy (FA), mean diffusivity (MD), radial diffusivity (RD), and axial diffusivity (AD) using voxel-based analyses in LPs compared to NLPs (AlphaSim corrected with a voxel-level threshold of and a cluster-level threshold of ). Blue and red regions denote decreased and increased DTI parameters, respectively. The numbers at the top indicate the value of MNI coordinates. Figure SII: brain regions that showed altered functional connectivity with FG.R (A) and SFG.R (B), respectively, in LPs compared to NLPs. Two sample t-tests were performed to explore the intergroup differences. The results were corrected by the AlphaSim method at a voxel level of . All abbreviations of the brain regions are shown in Tables SI and SII. (Supplementary Materials)