prelingual deafness is the hearing loss that occurs at birth or before the onset of speech . due to the deprivation of auditory inputs , cross - modal plasticity  has been suggested in deaf subjects , which is represented by the phenomenon that the auditory cortex can be activated when deaf subjects perform various tasks , such as speech [ 1 , 2 ] and visual tasks [ 3 , 4 ] . many morphological studies did not find structural changes in the primary auditory cortex [ 58 ] , indicating the atrophy of the auditory cortex due to hearing loss could be compensated by the use of this cortex for other stimuli . however , brain regions involved in visual and speech processing have been found to change in deaf subjects [ 7 , 8 ] , indicating that the sensory systems for vision and speech might participate in the tasks which are supposed to activate the auditory system in healthy controls . therefore , we speculated that the cooperative manner for different brain regions would be altered in prelingually deaf subjects when dealing with a complex task . the human brain is a highly complex system with synchronized neural activity from different brain regions . the concept of  connectome  was first proposed by sporns et al . in 2005 , which represents the human brain as an interconnected network . the network architecture identified by the analysis of functional connectivity could be an effective pattern to present the cooperative manner for brain regions , and it reflects the potential anatomical connections between brain regions as well . to further quantitatively measure the brain network , network topological properties could be evaluated at both global and regional levels using graph theory , which becomes a promising tool for analyzing brain networks in recent years [ 1021 ] . small - worldness , a concept that originated from social network , quantifies the effectiveness of information transferring within networks and has been successfully used to characterize brain networks [ 1416 ] . besides , nodal topological parameters based on graph theory analysis are used to illustrate the properties for nodes which are defined as brain regions , identifying the role of specific regions in transferring information . it has been demonstrated that these measures of brain network are sensitive to aging [ 17 , 18 ] as well as various neuropsychiatric diseases , such as schizophrenia [ 16 , 19 ] and alzheimer 's disease [ 20 , 21 ] . therefore , investigation of brain network properties provides a new insight into brain reorganization and is critical to understand the working mechanism of brains with hearing loss . recently , a few studies have emerged to investigate the structural or functional connectivity between brain regions in deaf subjects . kim et al .   examined morphological brain network in deaf adults using tissue density on mri and analyzed the network properties using graph theory and network filtration . they found altered morphological network in prelingually deaf adults compared to normal controls but not in postlingually deaf adults , concluding that auditory experience could affect the morphology of brain networks in deaf adults . employed resting - state fmri to investigate the effect of deafness on the intra- and interregional synchronization of different parts of superior temporal sulcus and revealed that the intrinsic function of these different parts are distinctly impacted by deafness . however , these previous studies focused on alterations in grown adult brains of deaf subjects but not developing brain . adolescence is an important transitional stage for brain development , in which brain 's structure and function are in developmental changes to accommodate to external environment . therefore , adolescents with prelingual deafness are supposed to go through a different pattern of brain reorganization to adapt to the silent world relative to normal adolescents , which is supported by one of our previous studies   using structural mri to investigate grey matter connectivity within and between auditory , language , and visual systems in deaf adolescents . to the best of our knowledge , the whole - brain functional network and its properties have not been investigated in prelingually deaf adolescents yet . we hypothesized the architectures or properties of brain network would change in prelingually deaf adolescents . in the present study , we used resting - state fmri and graph theory analysis to investigate alterations of whole - brain functional networks in prelingually deaf adolescents . functional brain network was constructed based on mean time series extracted from 90 brain regions . functional connectivity patterns were studied in deaf and control groups , respectively , and compared between two groups . the global network property of small - worldness was evaluated in each group and group differences were investigated as well . besides , nodal topological parameters were calculated and brain hubs were identified for deaf and control groups , respectively . nodal properties were finally compared in all the hubs between prelingually deaf adolescents and normal controls . in this study , we recruited 16 prelingually deaf adolescents and 16 normal controls with matched age ( p = 0.815 ) and gender ( chi - squared p = 1.0 ) , which were the same as the data used in our previous studies [ 7 , 8 , 24 ] . all prelingually deaf adolescents suffered from severe sensorineural hearing loss after the birth and had a mean pure tone audiometry average air conduction threshold greater than 90  db of hearing loss for the better ear and no single frequency better than 45  db of hearing loss between 500 and 4000  hz . furthermore , all subjects had neither history of central nervous system disease nor significant head trauma . the data from one prelingually deaf adolescent and one normal control were discarded due to excessive head movement during scanning . therefore , 15 deaf subjects ( age 14.32  2.24 years , 10~18 years ; 8 males ) and 15 normal controls ( age 14.81  2.07 years , 10~18 years ; 8 males ) remained , and these two groups were still matched in age ( p = 0.769 ) and gender ( chi - squared p = 1.0 ) . detailed biographical information of all subjects is shown in table 1 . this project was approved by the committee at beijing tongren hospital , and all subjects and their parents gave written informed consent before inclusion . resting - state fmri data were collected from all subjects on a 3-tesla mr imaging scanner ( ge medical system , milwaukee , wi , usa ) with an eight - channel phased - array head coil . to acquire resting - state fmri data , all subjects lay on their backs and were instructed explicitly to keep their eyes closed and not to think of anything in particular . the scanning parameters were as follows : repetition time ( tr ) = 2000  ms , echo time ( te ) = 30  ms , field of view ( fov ) = 24  24  cm , matrix size = 64  64 , slice thickness = 5  mm , and flip angle = 90 , yielding 28 axial slices with in - plane resolution of 3.75  3.75  cm . the fmri scanning lasted 400  s and 200 volumes were obtained for each subject . furthermore , high - resolution t1-weighted structural brain images were collected for structural reference using a 3d spgr sequence ( tr = 9  ms , te = 3.5  ms , inversion time ( ti ) = 450  ms , fov = 24  24  cm , matrix size = 256  256 , slice thickness = 1  mm , and flip angle = 13 , yielding 196 sagittal slices with in - plane resolution of 0.9375  0.9375  cm ) . the preprocessing steps were performed using the statistical parametric mapping ( spm8 , http://www.fil.ion.ucl.ac.uk/spm , wellcome department of cognitive neurology , london , uk , 2008 ) on a matlab 7.9 platform ( mathworks , natick , ma , usa ) . the first 10 volumes of the resting - state fmri data were discarded because of the initial transient effects and the adaption of the subjects to the environment . the remaining 190 time points were corrected for within - scan acquisition time differences between slices and realigned to the first image for head movement correction , resulting in 3 translational and 3 rotational parameters . the subjects with excessive head movement ( translation > 2  mm and rotation degree > 3 ) were thus excluded . subsequently , the functional scans were spatially normalized to a standard space ( montreal neurological institute , mni ) using the normalization parameters from t1 image with high resolution to structural t1 template in mni space and resampled to 3  3  3  mm . the images were then smoothed with a 6 mm full width at half maximum ( fwhm ) isotropic gaussian kernel . after that , the preprocessed images were linearly drifted and passed through a band - pass filter ( 0.010.08  hz ) to remove the effects of low frequency drift and physiological high frequency respiratory and cardiac noise . we regressed out six parameters of head movement as well as mean signals of the whole brain , white matter , and cerebral spinal fluid . the functional images were parcellated into 90 brain regions ( 45 regions for each hemisphere ) based on automated anatomical labelling ( aal ) template ( see table 2 ) . to construct the whole - brain network , time courses were first extracted from all voxels within each roi and averaged . then , pearson 's correlation coefficient was calculated between pairs of regions throughout the whole brain , resulting in a 90  90 correlation matrix . larger correlation coefficients indicate more synchronized time courses between the corresponding pairs of regions , implying stronger functional connectivity between these two regions . brain network properties were analyzed using graph theory based on the correlation matrix , in which 90 rois and their connections were considered as nodes and edges , respectively . the correlation matrix was binarized by a given threshold , resulting in a sparse matrix . using graph theory analysis , we characterized the global topological properties of brain functional networks using the parameters of small - worldness (  ) , which was defined as ( 1)==cpreal / cprandlpreal / lprand , where  and  are normalized clustering coefficient and normalized path length , respectively . the normalized clustering coefficient is the ratio of clustering coefficient of a real network cp to that of a random network cp , and the normalized path length is defined as the ratio of characteristic path length of a real network lp to that of a random network lp . here , clustering coefficient cp is a measure of the degree to which nodes in a network cluster together and is defined as(2)cp=1ni=1n2eidnodidnodi1,where dnod(i ) is the degree of the node i , which is defined in the following part of this section , ei is the total number of edges connecting the node i with the nodes which are the nearest neighbors to the node i , and n is the number of nodes in the network . characteristic path length lp is the average shortest path length between two nodes over all pairs of nodes and is defined as(3)lp=ijlijnn1,where n is the number of nodes in the network and lij is the shortest path length between nodes i and j. a network is said to have the property of small - worldness if it satisfies   1 and   1 or  = /  1 . it has been demonstrated that the sparsity of 12% is close to an optimal threshold at which most informative network edges are retained and disconnections between nodes are rare [ 26 , 27 ] . therefore , we chose 12% as the sparsity threshold for analysis of nodal topological properties , including betweenness centrality , nodal degree , and nodal efficiency , which are defined as follows . betweenness centrality cb is a measure of the number of times that a node is along the shortest path between other two nodes , reflecting the importance of the individual node in the overall network structure . it is computed as follows:(4)cb=sitnstist , where n is the set of all nodes in the network , st is the total number of shortest paths from node s to node t , and st(i ) is the number of the paths that pass through the node i. 
 ( 2 ) nodal degree . nodal degree of the node i ( dnod(i ) ) in a network is the number of connections linking the node with all other nodes and is defined as(5)dnodi=jnaij , where n is the set of all nodes in the network and aij is the connection status between node i and node j ( aij = 1 when the connection exists , and aij = 0 otherwise ) . nodal efficiency of the node i ( enod(i ) ) is the mean of the inverse shortest path length from node i to all other nodes and is defined as(6)enodi=1nji1lij , where n is the number of nodes in the network and lij is the shortest path length between nodes i and j. when constructing whole - brain network using resting - state fmri , the correlation matrix characterized the brain functional connectivity pattern for each subject . to identify significant connectivity pattern in each group , one - sample t - test was performed to find significantly positive correlation in deaf group and normal controls , respectively , indicating significant connections between pairs of regions throughout the whole brain . to compare differences of connectivity patterns between groups , fisher 's r - to - z transformation was first performed and then two - sample t - test was conducted on the correlation matrices between prelingually deaf adolescents and normal controls . these statistical results were corrected for multiple comparisons using partial bonferroni correction which considers the correlation between the dependent variables . this correction method was implemented using the simple interactive statistical analysis bonferroni tool ( http://www.quantitativeskills.com/sisa/calculations/bonfer.htm ) , which optimally balanced type - i and type - ii errors . when performing one - sample t - test , the mean correlation coefficient between all variables in the matrix was 0.2236 in control group and 0.2274 in deaf group , leading to an equivalent corrected p of 7.94e  5 for control group and 8.63e  5 for deaf group ( number of tests = 90  89/2 = 4005 ) , respectively . for analysis of group differences , the mean correlation coefficient between pairs of variables for all subjects was 0.1624 , resulting in an equivalent corrected p of 4.83e  5 . to present the properties of brain functional network for each group , the property of small - worldness was analyzed at different sparsities ( from 10% to 45% ) for prelingually deaf adolescents and normal controls , respectively . then , two - sample t - test was performed between two groups to figure out the group differences . as described above , we chose 12% as the sparsity threshold for analysis of nodal topological properties . first , the mean betweenness centrality , mean nodal degree , and mean nodal efficiency were calculated for each node in deaf and normal groups , respectively . then , we ranked them among all the nodes with the highest score ( 90 ) for the highest betweenness centrality , nodal degree , and nodal efficiency , respectively , and the scores for each measure were summed up for each node . the top 20% nodes ( 18 out of 90 nodes ) with high scores were identified as hub nodes in each group . the three nodal topological parameters were then compared between two groups for the hubs using two - sample t - test and the results were corrected for multiple comparisons using partial bonferroni correction as well . using one - sample t - test , significant functional connectivity was detected for both prelingually deaf adolescents and normal controls after partial bonferroni correction ( as shown in figure 1 ) . the 3d view in figure 1 was visualized with the brainnet viewer ( http://www.nitrc.org/projects/bnv ) . the figure shows functional connectivity was widely distributed throughout the whole brain in both groups . to identify group differences of functional connectivity throughout the whole brain , connectivity matrices , we can see that prelingually deaf adolescents show weaker connectivity between the regions within the temporal lobe , but stronger connectivity between regions in the occipital lobe and regions in the limbic system and temporal lobe . after partial bonferroni correction , prelingually deaf adolescents showed significantly increased functional connectivity between the right superior gyrus ( spg ) and right insula ( ins ) and between the left middle temporal gyrus ( mtg ) and right posterior cingulate gyrus ( pcg ) when compared with normal controls ( see figure 3 ) . besides , significantly decreased connectivity was found in deaf group between the right spg and left middle frontal gyrus ( orbital part , mfgorb ) and between the right postcentral gyrus ( pocg ) and left inferior frontal gyrus ( opercular part , ifgoper ) . small - worldness  of brain functional network was measured at different sparsity ( from 10% to 45% with an interval of 1% ) in deaf group and control group , respectively . as shown in figure 4 , we found   1 across all the sparsities in both groups , indicating that the network has the property of small - worldness . we did not find significant differences of the small world index between prelingually deaf adolescents and normal controls . besides the global network property , the regional topological properties , including betweenness centrality , nodal degree , and nodal efficiency , figure 5 shows the mean values for these regional properties in control and deaf groups . the hub in the network is supposed to have relatively high betweenness centrality , nodal degree , and nodal efficiency . based on the method described in section 2 , we identified 18 regions out of 90 rois as the hubs in each group at the sparsity of 12% since a network with this sparsity tends to retain the edges with most information and ensure disconnection is rare . as shown in figure 6 and table 3 , ten brain regions , including the right rolandic operculum ( rol ) , left gyrus rectus ( reg ) , bilateral calcarine fissure ( cal ) , right lingual gyrus ( ling ) , left putamen ( put ) , right pallidum ( pal ) , bilateral heschl gyrus ( he s ) , and right superior temporal gyrus ( stg ) , were the hubs in both deaf and control groups . besides , the middle and superior frontal gyri and cuneus ( cun ) were hubs in deaf group but not in control group , while the right precentral gyrus ( precg ) , bilateral hippocampus ( hipp ) , and supramarginal gyrus ( smg ) were hubs only in control group . there were a total of 26 brain regions identified as hubs either in deaf or in control groups . the nodal topological parameters , including betweenness centrality , nodal degree , and nodal efficiency , were compared for these nodes between prelingually deaf adolescents and normal controls . significant differences did not survive after correction for multiple comparisons . however , we found trends for group differences , as listed in table 3 . betweenness centrality is higher in the left superior frontal gyrus ( orbital part , sfgorb ) ( p = 0.0204 , not corrected ) and is lower in the left rol ( p = 0.0480 ) , smg ( p = 0.0251 ) , and pal ( p = 0.0221 ) in deaf group relative to control group . besides , prelingually deaf adolescents have higher nodal degree in the left middle frontal gyrus ( orbital part , mfgorb ) ( p = 0.0484 ) and left superior frontal gyrus ( medial part , sfgmed ) ( p = 0.0302 ) and lower nodal efficiency in the left rol ( p = 0.0367 ) , left hippocampus ( p = 0.0442 ) , and left pallidum ( p = 0.0015 ) . the present study employed resting - state fmri and graph theory analysis to investigate changes of brain functional network in prelingually deaf adolescents . functional connectivity was widely distributed in both deaf and control groups , and significant differences were found for connections between the right spg and right insula , the left mtg and right pcg , the right spg and left mfgorb , and the right pocg and left ifgoper . the whole - brain network has the property of small - worldness in both prelingually deaf adolescents and normal controls . ten regions were identified as the hubs in both groups . while the middle and superior frontal gyri and cun were the hubs that emerged only in deaf group , the right precg , bilateral hipp , and smg were hubs only in control group . these findings provide new evidence for brain reorganization in prelingually deaf adolescents and could help to understand the working mechanism in the period of brain development after hearing loss . in prelingually deaf adolescents , functional connectivity was widely distributed as observed in normal controls but differed in the connectivity pattern . significantly increased connectivity was detected between the right spg and right insula and between the left mtg and right pcg . the insula and pcg are in limbic system , which are thought to play an important role in recruiting relevant brain regions for sensory information processing . the spg is supposed to receive considerable visual input and is involved in visual - spatial relations . besides , the mtg has been reported to be one of the critical nodes in the brain 's language network , which can access word meaning while reading . therefore , our findings of increased connectivity between these regions suggest that the limbic system in prelingually deaf adolescents could reinforce the processing of the visual and verbal information due to the deprivation of auditory sensory processing . additionally , we found decreased functional connectivity between the right spg and left mfgorb and between the right pocg and left ifgoper . evidence suggests that the left mfg is the specific brain region for chinese reading , and the left ifg is extremely important for language comprehension and production . the decreased connectivity might be caused by the inferior skills of reading or speaking observed in most of the deaf subjects [ 34 , 35 ] . although it is not significant , we found a trend of significance for group differences in functional connectivity pattern . connectivity between the regions within the temporal lobe was found to be weaker in prelingually deaf adolescents , indicating that the auditory cortex are less activated to auditory stimuli . in contrast , stronger connectivity was observed between regions in the occipital lobe and regions in the temporal lobe , which were identified as primary visual and auditory cortices , respectively . this result is consistent with one of our findings in the previous study , which detected significantly stronger structural connections between visual and auditory systems in prelingually deaf adolescents . using graph theory analysis , whole - brain network holds the property of small - worldness in both prelingually deaf adolescents and normal controls . the small world index was almost the same across all the sparsities from 10% to 45% ( with an interval of 1% ) between the two groups ( see figure 3 ) . this result suggests that the effectiveness of information transferring of the brain network was not affected even without auditory inputs in prelingually deaf adolescents . in the brain network , hubs are regions playing an important role in facilitating communication among parallel , distributed brain networks . the hubs included the right precg , bilateral rol , left reg , bilateral hipp , bilateral cal , right ling , bilateral put , bilateral pal , bilateral he s , and right stg . compared with healthy adolescents , middle and superior frontal gyri , left ling , left cun , left stg , and right reg were hubs only observed in prelingually deaf adolescents . evidences show that the stg is the primary auditory cortex and is involved in phonological processing for speech perception and production as well [ 37 , 38 ] . although not activated by auditory stimuli , auditory regions could be activated by other stimuli [ 14 ] as stated in section 1 . therefore , auditory regions might become a brain hub connecting with other brain regions . in the hubs of deaf subjects , the ling and cun are associated with visual processing [ 39 , 40 ] . to adapt to deficient auditory input , prelingually our findings of the hubs in ling and cun suggest the improvement of visual performance in deaf adolescents , supporting cross - modal changes in visual processing . of note , it has been demonstrated that the middle and superior frontal gyri can be activated during observing and imitating csl . we speculated that the hubs of mfg and sfg only observed in prelingually deaf subjects might be caused by the use of csl . we recruited 16 prelingually deaf adolescents and 16 normal controls , and 15 subjects remained for each group after removing subjects with excessive head movement . second , the effects of clinical parameters , such as the information for the use of sign language and hearing aids , were not analyzed in the present study , which will be taken into consideration in future . the present study employed resting - state fmri and graph theory to investigate alterations of brain functional network in prelingually deaf adolescents compared to normal controls . increased connectivity was significantly found between the limbic system and visual and language - related regions , suggesting reinforcement of the processing for the visual and verbal information in deaf adolescents . decreased connectivity was detected between the visual regions and language regions possibly due to inferior reading or speaking skills in deaf subjects graph theory analysis revealed that brain functional network holds small - worldness in both prelingually deaf adolescents and normal controls . compared with healthy adolescents , middle and superior frontal gyri , left ling , left cun , left stg , and right reg were hubs only observed in prelingually deaf adolescents . these hubs were involved in visual , language , and auditory processing , reflecting brain reorganization to adapt to deficient auditory input .
previous neuroimaging studies suggested structural or functional brain reorganizations occurred in prelingually deaf subjects . 
 however , little is known about the reorganizations of brain network architectures in prelingually deaf adolescents . 
 the present study aims to investigate alterations of whole - brain functional network using resting - state fmri and graph theory analysis . 
 we recruited 16 prelingually deaf adolescents ( 10~18 years ) and 16 normal controls matched in age and gender . 
 brain networks were constructed from mean time courses of 90 regions . 
 widely distributed network was observed in deaf subjects , with increased connectivity between the limbic system and regions involved in visual and language processing , suggesting reinforcement of the processing for the visual and verbal information in deaf adolescents . decreased connectivity was detected between the visual regions and language regions possibly due to inferior reading or speaking skills in deaf subjects . using graph theory analysis 
 , we demonstrated small - worldness property did not change in prelingually deaf adolescents relative to normal controls . 
 however , compared with healthy adolescents , eight regions involved in visual , language , and auditory processing were identified as hubs only present in prelingually deaf adolescents . 
 these findings revealed reorganization of brain functional networks occurred in prelingually deaf adolescents to adapt to deficient auditory input .