- 1Department of Neurosurgery, The Affiliated Nanchong Central Hospital of North Sichuan Medical College, Nanchong, China
- 2Department of Neurosurgery, Southwest Hospital, Army Medical University, Chongqing, China
- 3Department of Neurosurgery, West China Hospital of Sichuan University, Chengdu, China
- 4Department of Geriatrics, The Affiliated Nanchong Central Hospital of North Sichuan Medical College, Nanchong, China
- 5Department of Radiology, West China Hospital of Sichuan University, Chengdu, China
- 6Department of Radiology, The Affiliated Nanchong Central Hospital of North Sichuan Medical College, Nanchong, China
Objective: Acoustic neuroma (AN) is a common benign tumor. Little is known of neuropsychological studies in patients with acoustic neuroma, especially cognitive neuropsychology, and the neuropsychological abnormalities of patients affect their life quality. The purpose of this study was to explore the changes in the cognitive function of patients with acoustic neuroma, and the possible mechanism of these changes by structural magnetic resonance imaging.
Materials and methods: We used a neuropsychological assessment battery to assess cognitive function in 69 patients with acoustic neuroma and 70 healthy controls. Then, we used diffusion tensor imaging data to construct the structural brain network and calculate topological properties based on graph theory, and we studied the relation between the structural brain network and cognitive function. Moreover, three different subnetworks (short-range subnetwork, middle-range subnetwork, and long-range subnetwork) were constructed by the length of nerve fibers obtained from deterministic tracking. We studied the global and local efficiency of various subnetworks and analyzed the correlation between network metrics and cognitive function. Furthermore, connectome edge analysis directly assessed whether there were differences in the number of fibers in the different brain regions. We analyzed the relation between the differences and cognitive function.
Results: Compared with the healthy controls, the general cognitive function, memory, executive function, attention, visual space executive ability, visual perception ability, movement speed, and information processing speed decreased significantly in patients with acoustic neuroma. A unilateral hearing loss due to a left acoustic neuroma had a greater impact on cognitive function. The results showed that changes in the global and local metrics, the efficiency of subnetworks, and cognitively-related fiber connections were associated with cognitive impairments in patients with acoustic neuroma.
Conclusion: Patients exhibit cognitive impairments caused by the decline of the structure and function in some brain regions, and they also develop partial compensation after cognitive decline. Cognitive problems are frequent in patients with acoustic neuroma. Including neuropsychological aspects in the routine clinical evaluation and appropriate treatments may enhance the clinical management and improve their life quality.
Introduction
Acoustic neuroma (AN), also known as vestibular schwannoma, is a benign tumor with a high incidence of 85% in the cerebellopontine angle, which mostly originates from the 8th cranial nerve sheath (Fortnum et al., 2009). About 90% of the patients present gradual progressive unilateral hearing loss (UHL)(Suzuki et al., 2010). AN is a common cause in UHL patients. UHL is defined as an asymmetric loss of hearing in which one side has a hearing loss and the other has a normal hearing (Vincent et al., 2015). The incidence of UHL is 7.2% in the United States (Golub et al., 2018). The difficulty of understanding and communicating in a noisy environment and the decreased ability of the sound source localization in patients with UHL, which affect the ability to receive important information from both sides simultaneously (Nelson et al., 2019).
Previous studies on patients with bilateral deafness have confirmed that the severity of hearing loss was positively correlated with cognitive impairment. Cognitive impairment in presbycusis with moderate and severe hearing loss is 1.4 and 1.6 times higher than in healthy people, respectively (Davies et al., 2017). Presbycusis is related to the increased risk of cognitive impairment, including decreased executive function (Gurgel et al., 2014), memory (Lin et al., 2011), and psychomotor processing disorders (Chen et al., 2018). Early intervention of the hearing can improve the cognitive function and behavior state in bilateral presbycusis patients (Ma et al., 2016; Adrait et al., 2017). Studies have shown that a partial or unilateral decrease in hearing in patients with UHL, which can affect the central auditory cortex and remodel vision and other sensory systems (Propst et al., 2010), even cognitively-related brain regions (Zhang et al., 2018), resulting in the decline in advanced cognitive functions such as memory, language, and learning. Children with UHL have poorer language development and cognition than normal children, and they have problems such as language retardation and inattention (Carew et al., 2018; van Wieringen et al., 2019), among whom the proportion of children with behavioral problems (25%) is higher than the normal children (Lieu et al., 2012). The connections in the brain network of executive function, language comprehension, and cognition in children with UHL are different from those of normal children (Jung et al., 2017).
To date, there are few reports on whether the hearing loss in UHL, especially AN patients, will cause cognitive dysfunction; the conclusions from previous studies were different; the number of reported cases was low; and little was known about the mechanism. In this study, we analyzed the cognitive function between AN patients and healthy controls (HCs) using various neuropsychological tests. Furthermore, we first systematically explored the possible underlying mechanism of cognitive impairment in AN patients using diffusion MRI. We constructed the structural brain network based on diffusion tensor imaging (DTI) data, and analyzed the network based on graph theory and other related analyses of brain networks, in order to study the brain connectomics of cognitive dysfunction in AN patients.
Materials and methods
Participants
Sixty-nine right-handed AN patients (44 females and 25 males, age range: 19–76 years) were recruited from the outpatient and ward of Neurosurgery of West China Hospital of Sichuan University, from October 2019 to July 2020. Seventy right-handed hearing controls (48 females and 22 males, age range: 26–74 years) were also enrolled in our study. The demographic information for these subjects is shown in Table 1; Supplementary Tables 2, 3. The experiment was approved by the hospital ethics committee, and all participants signed the informed consent.
The inclusion criteria of AN were as follows: (1) unoperated unilateral acoustic neuroma; (2) no history of mental or neurological diseases; (3) right-handedness; (4) no contraindications of MRI (such as spatial claustrophobia and metal implants). Healthy control (HC) group: (1) normal hearing; (2) right-handedness; (3) no history of craniocerebral trauma; (4) 18–75 years old.
Exclusion criteria were as follows: (1) previous diseases such as chronic otitis media that affect hearing threshold; (2) excluding the previous history of ear surgery, hearing loss caused by ototoxic drugs, or wearing hearing aids; (3) bilateral acoustic neuroma, such as neurofibromatosis type 2; (4) conductive hearing loss; (5) craniocerebral trauma; (6) tumors with other intracranial sites and nature; (7) excluding the history of craniocerebral surgery; (8) those who could not cooperate to complete the neuropsychological test (including daltonism and color weakness).
The average air conduction thresholds at four frequencies (0.5, 1, 2, and 4 kHz) were calculated as the pure tone average (PTA), representing the hearing levels of the subjects. According to the World Health Organization, hearing loss was classified as mild (PTA 26-40 dB HL), moderate (PTA 41–60 dB HL), severe (PTA 61–80 dB HL), profound (PTA > 81 dB HL).
Patients with acoustic neuroma are often complicated with tinnitus symptoms. In this study, AN patients accompanying tinnitus were assessed using the tinnitus handicap inventory (THI) scale (Newman et al., 1996). Higher scores indicate greater severity and greater impact on daily life.
Cognitive functions of all subjects were assessed, including the following measures: Montreal cognitive assessment (MoCA), Rey auditory verbal learning test (RAVLT) immediate memory and delayed memory, Stroop color-word test A, B, and C (Stroop A, B, and C), symbol digit modalities test (SDMT), Trail-Making Test A and B (TMT A and B), Hamilton depression scale (HAMD), and Hamilton anxiety scale (HAMA).
MRI data acquisition
The data acquisition parameters were as follows: DTI and 3D high-resolution T1WI were acquired on GE 750 W 3.0Tesla magnetic resonance equipment (General Electric Medical System, Milwaukee, WI, United States) with a 32-channel head coil. Acquisition parameters used for DTI sequence were as follows: TR = 7,000 ms, TE = 72 ms, FOV = 24 cm × 24 cm,acquisition matrix = 256 × 256,flip angle = 90°,slice thickness / slice spacing = 4.0 mm/0.0 mm (no intervals). A total of 34 slices were scanned. Diffusion-sensitive gradients were applied in 51 directions (b = 1,000 s/mm2) and one b0 (b = 0 s/mm2) images. T1 scanning parameters were as follows: slice thickness = 1 mm, scanning matrix = 512 × 512, voxel size = 0.5 × 0.5 × 1.0 mm3.
Data preprocessing
The MRI data preprocessing and network construction were performed using PANDA (www.nitrc.org/projects/panda; Cui et al., 2013). The main steps were as follows: data quality check, data format conversion, and head eddy-current effect correction. The whole brain fiber bundle was tracked using deterministic fiber tracking (Mori et al., 1999). The main results tracked by PANDA software were fiber number (FN) and fiber length (FL). The next steps were to explore the possible mechanism of cognitive function changes in patients with AN from different perspectives based on these results. Firstly, the structural brain network was constructed using the fiber number, and we analyzed the network properties and the correlation between the network properties and cognition function. Secondly, we directly compared the FN in any two brain regions between AN patients and healthy controls. Thirdly, the FL was used to construct short-range, middle-range, and long-range subnetworks, and the relationships between subnetwork topological properties and cognition decline were discussed.
Network construction and graph theory metrics calculation
The FN was calculated by the deterministic fiber tracking and was registered to the individual Anatomical Automatic Labeling 90 regions (AAL90). Due to tumors compressing the cerebellum and brainstem in some patients, it may affect the study of the cerebellum, so we constructed networks based on the 90 cerebral regions of the AAL atlas, excluding the cerebellar regions. The FN value of every two brain regions constructed a weighted matrix. The connections were considered effectively structurally connected if at least three fibers in two brain regions in 80% of the subjects (Shu et al., 2015), transforming the weighting matrix into a binary matrix. Network topological properties were calculated using the GRETNA toolbox (http://www.nitrc.org/projects/gretna/), including global network metrics (Shu et al., 2015; Wang et al., 2015a): 1. small world property: shortest path length (Lp), clustering coefficient (Cp), and small-world parameters (λ, γ, and σ). 2. global efficiency (Eg), local efficiency (Eloc); local metrics: node local efficiency, node efficiency, node shortest path, node clustering coefficient, and degree centrality (see Supplementary Table 1).
Meanwhile, for the weighted FN matrix, we compared the FN between 90 × 90 brain regions using network-based statistic (NBS) multiple comparison correction (Zalesky et al., 2010). The NBS procedure was performed as follows: (1) admit weighted connections with statistical analysis surpassing p < 0.001 between groups; (2) seek clusters that suprathreshold connections; and (3) via permutation testing (5,000 permutations), compute p values (FWE-corrected) for each cluster. Significant differences between groups were performed using p < 0.05 to control the FWE after NBS correction.
The FL was classified according to the following criteria: 1. short-range subnetwork: fiber length less than 40 mm 2. middle-range subnetwork: fiber length 40-80 mm; 3. long-range subnetwork: fiber length was greater than 80 mm (Wang et al., 2013). The global and local efficiency of these subnetworks were calculated, respectively. And Spearman correlation was assessed between the network metrics and the cognitive scale.
Statistical analysis
SPSS (version 23.0, IBM, United States) was used for statistical analysis. Data are presented as mean ± standard deviation for normally distributed variables or median (interquartile range) for not normally distributed variables. Categorical variables were expressed as percentage. When comparing two groups of independent data, t-test (two-tailed) was used if the data were satisfied normal distribution with homogeneous variance, otherwise, Mann–Whitney U nonparametric test was performed. According to the comparison among the two subgroups of AN (LAN&RAN) and the HC groups, one-way ANOVA or Kruskal-Wallis H test was used to test whether the data were satisfied normal distribution and homogeneous variance. The qualitative data were compared by χ 2 test. Spearman correlation analysis was used to explore the relationship between clinical data and cognitive function. The statistically significant p value was less than 0.05.
The global metrics of the brain network were statistically analyzed by using the SPSS23.0 and the GRETNA software (http://www.nitrc.org/projects/gretna/; Wang et al., 2015b) was used for statistical analysis for local metrics. The comparison of local metrics is p < 0.05, using a False Discovery Rate (FDR) for multiple comparison correction.
Results
Demographic characteristics
A total of 139 subjects were included in this study, including 69 patients with AN (LAN: RAN = 44: 25) and 70 HC patients. No significant differences in gender, age, and years of education between the two groups or among the left acoustic neuroma (LAN), right acoustic neuroma (RAN), and HC groups were attested (p > 0.05), as reported in Table 1, Supplementary Table 2. There were no significant differences in the course of the disease, PTA on the affected side, and THI scores between the LAN and RAN groups, as displayed in Table 1.
Comparison of cognitive function and correlation analysis between AN and HC groups
Compared with the HC group, the patients with AN performed worse in MoCA, RAVLT, Stroop, SDMT, and TMT (p < 0.05; See Supplementary Tables 2, 3), and both groups (LAN and RAN) also performed worse in MoCA, RAVLT, Stroop B, C, SDMT, and TMT (p < 0.05; see Table 1; Supplementary Table 4). Spearman correlation analysis showed that the left-sided PTA of patients with LAN was negatively correlated with MoCA subscores on visuospatial executive and delayed recall, and SDMT, and positively correlated with Stroop A and B. In RAN patients, the right-sided PTA was negatively correlated with MoCA subscores on orientation, RAVLT, and SDMT, while positively correlated with Stroop B and TMT-A. See Figure 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/970159/fnagi-14-970159-HTML/image_m/fnagi-14-970159-g001.jpg)
Figure 1. Relation between pure tone average (PTA) and cognition. (A) Scatter plots representing the distribution of PTA in the left acoustic neuroma (LAN) and right acoustic neuroma (RAN) groups. The numbers in the box represent the average PTA of the affected side in LAN and RAN patients, respectively. (B) Scatterplots showing the associations between right PTA and cognitive scale in the RAN patients. (C) Scatterplots showing the correlations between left PTA and cognitive scale in the LAN patients. The correlation analysis between PTA and cognition shows that hearing loss in AN patients affects cognitive functions such as visuospatial executive, memory, attention, motor speed, visual perception, and executive control. Error bars in (A) represent mean ± standard deviation. The gray regions in (B,C) show a 95% confidence interval. *, p < 0.05;**, p < 0.01.
Comparison of cognitive function among patients with different grades of AN and HC groups
According to Koos et al., (1993), all AN patients were graded into four groups: grade 1 (tumor diameter < 1 cm) in 1 cases; grade 2 (tumor diameter 1-2 cm) in 16 cases; grade 3 (tumor diameter 2–3 cm) in 17 cases; grade 4 (tumor diameter > 3 cm) in 35 cases. Because the number of patients with grade 1 was too small to analyze statistically, the cognitive functions of patients of other grades were compared with those of the HC group. The results showed that compared with the HC group, the cognitive function of patients with grades 2–4 decreased, as listed in Supplementary Tables 5, 6. To investigate whether tumor size affects cognitive function, tumor size was defined on the basis of the longest diameter, and Spearman correlation analysis was performed between tumor diameter and cognitive function, and no significant correlation was found between tumor size and cognitive function.
Comparison of cognitive function among patients with different degrees of hearing loss and HCs
According to WHO grade (1997), the hearing of the affected side (tumor side) of the patients with AN were as follows: normal hearing (PTA < 25 dB HL) in 5 cases, mild loss (PTA 26–40 dB HL) in 12 cases, moderate loss (PTA 41–60 dB HL) in 9 cases, severe loss (PTA 61–80 dB HL) in 13 cases, profound loss (PTA > 81 dB HL) in 10 cases. Compared with the HC group, the cognitive function of AN patients with mild to profound hearing loss decreased to various degrees, see Supplementary Tables 7, 8.
Results of global metrics
Small world property
The small world property with LAN was tested by one-sample t-test, the same for RAN. The results were p < 0.001, showing that both left and right AN patients have obvious “small world property.” No significant differences were found among the LAN, RAN, and HC groups. The small-world networks can highly effectively integrate, segregate, and transmit information with low cost. The above results indicated that patients with auditory neuroma still have the characteristics of high efficacy and less energy loss. Compared with the HC group, the Cp and λ were significantly lower in RAN patients (p < 0.001, p = 0.002, respectively), and the Lp was significantly higher in LAN patients (p = 0.003). See Figure 2.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/970159/fnagi-14-970159-HTML/image_m/fnagi-14-970159-g002.jpg)
Figure 2. The comparisons of global metrics in three groups. (A) Scatter plots representing the distribution of global efficiency (Eg) and local efficiency (Eloc) in LAN, RAN, and HC groups. (B) Scatter plots representing the distribution of small world properties (σ) in the LAN, RAN, and HC groups. (C) Eg and Eloc of structural brain network for three groups. Eg and Eloc in AN patients decreased compared to the HC group. (D) Shortest path length (Lp), clustering coefficient (Cp), and small-world parameters λ (λ = Lpreal / Lprandom, Lpreal: the Lp of the real network; Lprandom: the Lp of random network) of structural brain network for three groups. The numbers in the box represent mean values. Error bars represent mean ± standard deviation.
Global efficiency (Eg) and local efficiency (Eloc).
Compared with the HC group, the Eg of the LAN and RAN groups decreased (p = 0.009, 0.003, respectively); the Eloc of the RAN group decreased significantly (p < 0.001). The Eloc of the RAN group was lower than that of the LAN group (p = 0.036). See Figure 2.
Results of local metrics
The results of local metrics were as follows:
Compared with the node efficiency of the HC group, both left and right AN patients showed a more extensive decrease, mainly in the frontal lobe, occipital lobe, parietal lobe, limbic system, basal ganglia, thalamus, and so on. Only in the LAN group, the node efficiency of the left middle and inferior temporal gyrus increase, which was associated with compensation for auditory deprivation. The differential brain regions of node efficiency are shown in Figure 3B; Supplementary Tables 8, 9. Furthermore, compared with the HC group, all the other local metrics decreased, such as node local efficiency, node shortest path, node clustering coefficient, and degree centrality, see Supplementary information (p < 0.05, FDR corrected).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/970159/fnagi-14-970159-HTML/image_m/fnagi-14-970159-g003.jpg)
Figure 3. Correlations between global metrics and cognition. (A) Correlations between global metrics (Eg and Eloc) and cognitive scales (MOCA and RAVLT delay recall); Statistically significant Spearman correlation coefficients (n = 139, p < 0.05) are shown (MOCA: Montreal cognitive assessment; RAVLT: Rey auditory verbal learning test). (B) Schematic diagram of the differential node efficiency among the LAN, RAN, and HC groups. AN patients showed a more extensive decrease of node efficiency, mainly in the frontal lobe, occipital lobe, parietal lobe, limbic system, basal ganglia, thalamus.
Correlations between graph theory metrics and cognitive performance
The results showed that the Eg and Eloc were positively correlated with the MoCA (r = 0.180, p = 0.034; r = 0.207, p = 0.014, respectively), and the Eloc was positively correlated with RAVLT delayed recall (r = 0.209, p = 0.014), as shown in Figure 3A.
The node efficiency of 90 brain regions can reflect the inherent property of each brain region to some extent, so Spearman correlation analysis was performed between the node efficiency of differential brain regions and cognitive performance. Widespread decreased node efficiency affected general cognitive function mainly in the frontal lobe, parietal lobe, insular, and limbic system, involving the default mode network (DMN), frontoparietal network, and salience network, which related to advanced cognition. The results are shown in Supplementary Tables 9, 10.
Results of connectome edge analysis
Using connectome edge analysis, fiber number in the brain region between AN and HC groups was directly compared. We defined positive or negative connections according to the increase or decrease in the number of fibers between brain regions. We found that compared with the HC group, the LAN and RAN patients had both positive and negative connections. See Figure 4 for details.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/970159/fnagi-14-970159-HTML/image_m/fnagi-14-970159-g004.jpg)
Figure 4. The results of connectome edge analysis. Ninety brain regions constitute the circle part of the circos diagram. Red lines represent the increase in the number of white matter fibers compared to the HC group, which are defined as positive connections, and blue lines represent the decrease in the number of fibers, which are defined as negative connections. The red and blue lines have the same meaning in (A–D), (A,B) Significant differences in the number of white matter fibers in 90 brain regions between LAN and HC groups are shown. The brain regions of negative connections mainly involve the default mode network, sensorimotor network, dorsolateral prefrontal loop, basal ganglia and thalamus, and medial temporal lobe memory (MTL) system, and are related to cognition declines such as memory, language, abstraction, and visuospatial ability. The positive connections mainly involve brain regions between the precuneus and occipital lobe fibers and memory-related fiber connections, suggesting that the body compensates for hearing loss and cognitive decline. (C,D) Significant differences between RAN and HC groups are shown. Fibers between the left lenticular putamen and the left hippocampus were the only negative connection in RAN patients. The fibers belonging to the MTL system are associated with memory. Positive connection regions mainly involve the occipital lobe, the frontostriatal circuits, the salience network, and the MTL system, and are associated with the visual system, attention, cognition switching, and memory, revealing that the body presents anatomical fiber connection enhancement to compensate for cognitive declines such as attention and memory.
Eg and Eloc of the three subnetworks
Compared with the HC group, the Eg of the LAN group decreased significantly in the short-range subnetwork (p = 0.002), middle-range subnetwork (p < 0.001), and long-range subnetwork (p = 0.001), and the Eloc of LAN patients in the short-range subnetwork decreased significantly (p = 0.002). There was no significant change in RAN patients. See Figure 5.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/970159/fnagi-14-970159-HTML/image_m/fnagi-14-970159-g005.jpg)
Figure 5. Correlations between the global metrics (Eg, Eloc) of short-range subnetwork and cognitive scale. (A) Correlations between the Eg of short-range subnetwork and cognitive scale. The Eg of short-range subnetwork was negatively correlated with cognition function (n = 139, p < 0.05). (B) Correlations between the Eloc of short-range subnetwork and TMT. MOCA: Montreal cognitive assessment; RAVLT: Rey auditory verbal learning test; SDMT: symbol digit modalities test; Stroop: Stroop color-word test; TMT: trail-making test.
Correlation analysis between metrics of subnetworks and cognitive scale
The Eg of short-range and middle-range subnetwork was positively correlated with MoCA, RAVLT, and SDMT, and negatively correlated with Stroop and TMT. The Eg of the long-range subnetwork and the Eloc of the short-range subnetwork were negatively correlated with Stroop C and TMT, respectively. See Figure 5; Supplementary Figure 3.
Discussion
The cognitive performance in AN patients
The LAN and RAN groups performed worse on MoCA, RAVLT, Stroop, SDMT, and TMT tests than those of the HC group, revealing that both LAN and RAN patients developed cognitive dysfunction such as memory, attention, executive function, movement speed, and information processing speed. Previous studies have established that UHL children are more likely to develop cognition dysfunction such as hypoprosexia, language development retardation, and poor learning. Goebel and Mehdorn (2018) enrolled 27 patients with AN and 18 patients with posterior fossa meningioma, and the authors found that most of the patients (69%) develop cognitive impairment, the most common performance is attention (such as alertness) and visual movement speed, which is consistent with our findings. Fan et al., (2020) investigated 25 LAN and 15 RAN patients and demonstrated that the general cognitive function is normal in the AN group, but the attention, information processing efficiency, executive function, and memory decreased. Our conclusion is inconsistent with theirs. We found that general cognitive function, attention, executive function, memory, visuospatial and perception, movement speed, and information processing speed in AN patients were significantly lower than those of normal people. In our study, the AN patients performed worse in more neuropsychological tests, maybe because our sample was larger and the statistical efficiency was relatively efficient. The course of the disease was negatively correlated with cognition, showing that with the development of the course of the disease, the ability of visuospatial executive, attention, memory, motor speed, and information processing speed decrease gradually, so early intervention may delay the effect of disease on cognitive function (Supplementary Figures 1, 2). We found that patients with different grades of AN had different degrees of cognitive decline, however, there was no correlation between tumor size and cognitive function, demonstrating that tumor size was not the main factor for cognitive decline in AN patients. Meanwhile, this study also found that cognitive function decreased in patients with different degrees of hearing loss, and there was a significant correlation between PTA and cognitive function. It is suggested that hearing loss may have a greater effect on the cognitive function of patients with AN. Better hearing protection or reconstruction, may be better cognitive function.
The relation between global metrics and cognitive function
We found that AN patients had small-world properties, however, the average shortest path was higher than that of normal people. The Eg was lower than that of normal people, and the Eloc decreased in RAN patients. The small-world networks can highly effectively integrate, segregate, and transmit information with low cost (Bullmore and Sporns 2012). Generally, the fundament of brain organization efficiently is two basic principles: functional separation and integration (Tononi et al., 1994; Sporns 2011). The clustering coefficient and Eloc reflect the functional separation; the average shortest path and Eg reflect the functional integration. In our study, the clustering coefficient, the average shortest path, Eg, and Eloc changed, suggesting that functional separation and integration decreased in AN patients. The correlation analysis between global and local efficiency and cognitive scale in AN patients showed that global and local efficiency (Eg, Eloc) were positively correlated with MoCA, and local efficiency (Eloc) was also positively correlated with RAVLT delayed recall. The decrease in the global and local efficiency of patients led to general cognitive impairment, and the decrease in delayed recall ability was related to the decrease in local efficiency of the brain network. The global efficiency reflects the ability of the brain network to transmit information and the degree of network polymerization, and the local efficiency is used to measure the efficiency of distributed information processing in the whole brain structure network. The decrease in global and local efficiency in AN patients led to the decline of whole brain information transmission ability and information processing efficiency, which may be the reason for cognitive impairment.
The relation between local metrics and cognitive function
The brain regions with decreased node efficiency affected MoCA (representing general cognitive function) mainly in the frontal lobe, parietal lobe, insular, and limbic system, involving the default model network (DMN), frontoparietal network, salience network which are related to advanced cognition. We first found that the decline of immediate and delayed recall in RAVLT was mainly related to the decrease of node efficiency in the right supramarginal gyrus, dorsolateral superior frontal gyrus, opercular part of inferior frontal gyrus, rolandic operculum, insula, precentral gyrus, postcentral gyrus, superior parietal gyrus, cuneus, angular gyrus, and precuneus (Supplementary Tables 9, 10). The process of memory includes encoding and retrieval, and the memory retrieval brain region includes the right frontal lobe, anterior cingulate cortex, mid-parietal regions, and thalamus (Nyberg et al., 1996). Tulving et al., (1994) found that the left frontal lobe is more active in memory encoding and the right frontal lobe in memory retrieval, which is called HERA model by authors. We also found the brain regions related to RAVLT immediate and delayed recall were mainly located in the frontal lobe, insular, and parietal lobe. Smith et al., (1991) demonstrated that the visuospatial components of working memory are significantly activated in the right hemisphere, including the prefrontal lobe, premotor area, parietal lobe, and occipital lobe. In this study, the brain regions related to immediate memory were also located in the prefrontal lobe, premotor area, and parietal lobe, which was consistent with the literature. The decline in memory in AN patients was related to the decrease of node efficiency in brain regions involved in memory retrieval and working memory.
Connectome edge analyses and cognitive function
Previous studies have supported that the right parahippocampal gyrus, temporal pole, and fusiform gyrus are involved in the composition of memory (Aminoff et al., 2013). We had confirmed that decreased memory occurred in LAN patients in our study. The fiber connections between the right precuneus and the occipital lobe were strengthened, indicating that the compensation of the visual pathway increased after left auditory deprivation in LAN patients. The negative connections in LAN patients were as follows: (1) Right precuneus and left calcarine cortex, left posterior cingulate gyrus, postcentral gyrus. The precuneus and posterior cingulate gyrus are the key brain regions of the DMN, and their functions include memorizing, making social inferences, and looking forward to the future (Buckner and DiNicola 2019). The internal connections of DMN in LAN patients were weakened, which may affect the corresponding cognitive function. (2) Right precentral gyrus and supplementary motor area, dorsolateral superior frontal gyrus; right paracentral lobule and precentral gyrus and postcentral gyrus. These brain regions belong to the somatosensory motor system. Studies have affirmed that the sensorimotor system initiates and regulates the sensation and movement of the body, and plays an important role in speech processing such as vocabulary, phonetics, sentences, and chapter processing (Fischer and Zwaan 2006). During phonological processing, hearing syllables can activate the motor cortex corresponding to the mouth and tongue related to pronunciation (Pulvermüller et al., 2006). When understanding certain movements, the premotor area or motor cortex can be specifically activated (Shtyrov et al., 2014). Activations in the sensorimotor area increase during reading action-related texts (Speer et al., 2005). Therefore, we speculate that the decrease in sensorimotor network connections may lead to the decline of language and abstract ability in LAN patients. (3) Right thalamus and pallidum; left thalamus and pallidum, putamen. The fibers which involve projection fibers from the dorsolateral prefrontal lobe loop originate from the frontal lobe, pass through the caudate nucleus, pallidum, thalamus, and final project to the dorsolateral area of the frontal lobe. The function is to mediate executive function, and executive dysfunction is observed after injury (Grahn et al., 2009). The decline of executive function in AN patients may be associated with the negative connections in this loop. The subthalamic globus pallidus and globus pallidus subthalamic fibers are the afferent and efferent fibers of the globus pallidus, respectively, belonging to the basal ganglion nucleus. Numerous studies have confirmed that basal ganglia lesions may lead to cognitive function, including decreased memory (Grahn et al., 2009), visuospatial impairment (Su et al., 2007), and so on. (4) Left putamen and left hippocampus. The fibers belong to the medial temporal lobe system and are involved in memory (Sekeres et al., 2018).
The positive connections in RAN patients were as follows: (1) Right precuneus and right superior and middle occipital gyrus. The enhanced connections may indicate that the activation of visual fibers in RAN patients compensates for auditory deprivation. (2) Left orbital part of inferior frontal gyrus and anterior cingulate and paracingulate gyri, gyrus rectus. The fibers involve the frontostriatal circuits, including the orbitofrontal gyrus, anterior cingulate gyrus, striatum, ventromedial prefrontal lobe, and other nerve loops, which are mainly related to attention function (Fan et al., 2005). The enhancement of the connections between the left orbitofrontal and anterior cingulate gyrus may be linked to compensation of function after attention dysfunction in RAN patients. (3) Left anterior cingulate and paracingulate gyri and insular. Both regions are key nodes of the salience network (Seeley et al., 2007), which mainly plays a key role in the switching between cognitive function-related networks (Sidlauskaite et al., 2016). In RAN patients, the connections between the left anterior cingulate gyrus and insular were enhanced, suggesting that patients may need additional resources to maintain attention to the outside world after auditory deprivation. Fibers between the left lenticular putamen and the left hippocampus were the only negative connection in RAN patients. The fibers belonging to the medial temporal lobe system are associated with memory (Lisman et al., 2017). The negative connection may play a role in memory dysfunction in RAN patients.
The study of subnetworks and cognitive function
This study found that the Eg of the three subnetworks and the Eloc of the short-range subnetwork decreased in LAN patients, showing a decline in the ability of information integration and transmission in LAN patients, while no significant changes in RAN patients (Figure 6). The short-range subnetwork is usually located in the same brain region and reflects changes in local brain regions, and the middle-range subnetwork is distributed in the same functional region, while the long-range subnetwork connects different brain regions and is mainly responsible for information transmission (Su et al., 2015). These subnetworks ensure the separation and integration of brain functions. The subnetworks had changed in AN patients, affecting information transmission in the brain. In LAN patients, the Eg of the three subnetworks and the Eloc of the short-range subnetwork decreased, resulting in a decline in the ability of brain information transmission and integration.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/970159/fnagi-14-970159-HTML/image_m/fnagi-14-970159-g006.jpg)
Figure 6. The results of subnetworks. (A) Short-range subnetwork of one subject (B) middle-range subnetwork of the same subject. (C) long-range subnetwork of the same subject. (D) Comparison of global metrics (Eg and Eloc) in the three groups. The global metrics of the three subnetworks in AN patients significantly decreased. The numbers above the error bars represent median values. Error bars represent median ± 95%CI. * p < 0.05, ** p ≤ 0.01, *** p ≤ 0.001; LAN: left acoustic neuroma; RAN: right acoustic neuroma; HC: health control.
In this study, we found that the effects of LAN and RAN on the brain may be different. The LAN tends to decrease in structure and function, while the RAN tends to be compensated in structure and function, which is more resistant to hearing loss and more stable. The LAN patients may be more sensitive to auditory deprivation, and the RAN patients may tend to be functionally compensated. The RAN mainly affected the left dominant hemisphere and may have a stronger compensatory mechanism, so the RAN had less impact on patients and the performance was more stable. This was consistent with the literature (Wang et al., 2014; Zhang et al., 2015).
The shortcomings of this study are as follows: 1. This study is a cross-sectional study, which failed to follow up on the changes in cognitive function with the development of disease course and hearing level. However, the correlation analysis part provides a gradual decline in cognitive function with the progress of the disease course. 2. This study failed to follow up on whether the cognitive function of patients with acoustic neuroma recovered after hearing improvement.
Conclusion
In summary, in this study, we found that the cognitive function decreased in AN patients, and first systematically studied its possible mechanism. The results showed that widespread changes occurred in the cognition-related structural and functional regions of AN patients. The global efficiency and local efficiency of the structural brain network decrease in AN patients, which is closely related to cognitive dysfunction. In this study, using the edge analysis of brain structural networks in AN patients for the first time, we found that AN affected the fiber connections between cognitive-related brain regions. The global and local efficiency of subnetworks based on fiber length decreased, which affected patients’ general cognition, memory, execution, attention, processing speed, and visuospatial ability. The mechanisms of the left and right acoustic neuromas affecting the brain network are different. The LAN tends to decrease in structure and function, while the RAN tends to be compensated in structure and function, which is more resistant to hearing loss and more stable. Cognitive problems are frequent in AN patients. Including neuropsychological aspects in the routine clinical evaluation and appropriate treatment may enhance clinical management and improve the quality of life of patients.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by Ethics Committee of Nanchong Central Hospital. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author contributions
XD, XH, and LL designed the study. XD, JL, and HY collected the data. XD, SH, JS, and ZG analyzed the data. XD wrote the paper. WL, LL, and XH drafted the paper. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by Bureau of Science and Technology Nanchong city (grant no. 19SXHZ0273 and 20YFZJ0115), Nanchong Social Science Federation (grant no.NC21B188), Sichuan Province Medical Youth Innovative Research Project Program (grant no. Q21029), Primary Health Development Research Center of Sichuan Province (grant no. SWFZ20-C-069), and Special Funding for Postdoctoral Research Projects of Chongqing (grant no. 2021XM3012). Open access publication fees received from Nanchong Central Hospital.
Conflict 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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2022.970159/full#supplementary-material
Abbreviations
AN, Acoustic neuroma; DMN, Default ode network; DTI, Diffusion tensor imaging; Eg, Global efficiency; Eloc, Local efficiency; FN, Fiber number; FL, Fiber length; HC, Healthy control; LAN, Left acoustic neuroma; MoCA, Montreal cognitive assessment; PTA, Pure tone average; RAN, Right acoustic neuroma; RAVLT, Rey Auditory Verbal Learning Test; SDMT, Symbol-digital modalities test; TMT, Trail-making test; UHL, Unilateral hearing loss.
References
Adrait, A., Perrot, X., Nguyen, M. F., Gueugnon, M., Petitot, C., Collet, L., et al. (2017). Do hearing aids influence behavioral and psychological symptoms of dementia and quality of life in hearing impaired Alzheimer's disease patients and their caregivers? J. Alzheimer's Dis. 58, 109–121. doi: 10.3233/JAD-160792
Aminoff, E., Kveraga, K., and Bar, M. (2013). The role of the parahippocampal cortex in cognition. Trends Cogn. Sci. 17, 379–390. doi: 10.1016/j.tics.2013.06.009
Buckner, R., and DiNicola, L. (2019). The brain's default network: updated anatomy, physiology and evolving insights. Nat. Rev. Neurosci. 20, 593–608. doi: 10.1038/s41583-019-0212-7
Bullmore, E., and Sporns, O. (2012). The economy of brain network organization. Nat. Rev. Neurosci. 13, 336–349. doi: 10.1038/nrn3214
Carew, P., Mensah, F., Rance, G., Flynn, T., Poulakis, Z., and Wake, M. (2018). Mild-moderate congenital hearing loss: secular trends in outcomes across four systems of detection. Child Care Health Dev. 44, 71–82. doi: 10.1111/cch.12477
Chen, Y.-C., Chen, H., Jiang, L., Bo, F., Xu, J. J., Mao, C. N., et al. (2018). Presbycusis disrupts spontaneous activity revealed by resting-state functional MRI. Front. Behav. Neurosci. 12:44. doi: 10.3389/fnbeh.2018.00044
Cui, Z., Zhong, S., Xu, P., He, Y., and Gong, G. (2013). PANDA: a pipeline toolbox for analyzing brain diffusion images. Front. Hum. Neurosci. 7:42. doi: 10.3389/fnhum.2013.00042
Davies, H., Cadar, D., Herbert, A., Orrell, M., and Steptoe, A. (2017). Hearing impairment and incident dementia: findings from the English longitudinal study of ageing. J. Am. Geriatr. Soc. 65, 2074–2081. doi: 10.1111/jgs.14986
Fan, Z., Fan, Z., Li, S., Shi, Y., Liu, Y., and Ding, D. (2020). The effects of acoustic neuroma on cognitive function. Chin. J. Clin. Neurosci. 1, 52–58
Fan, J., McCandliss, B., Fossella, J., Flombaum, J., and Posner, M. (2005). The activation of attentional networks. NeuroImage 26, 471–479. doi: 10.1016/j.neuroimage.2005.02.004
Fischer, M., and Zwaan, R. (2006). Embodied language: a review of the role of the motor system in language comprehension. Q. J. Exp. Psychol. 2008, 825–850. doi: 10.1080/17470210701623605
Fortnum, H., O'Neill, C., Taylor, R., Lenthall, R., Nikolopoulos, T., Lightfoot, G., et al. (2009). The role of magnetic resonance imaging in the identification of suspected acoustic neuroma: a systematic review of clinical and cost effectiveness and natural history. Health Technol. Assess. 13, 1–154. doi: 10.3310/hta13180
Goebel, S., and Mehdorn, H. (2018). A missing piece? Neuropsychiatric functioning in untreated patients with tumors within the cerebellopontine angle. J. Neuro-Oncol. 140, 145–153. doi: 10.1007/s11060-018-2944-z
Golub, J., Lin, F., Lustig, L., and Lalwani, A. (2018). Prevalence of adult unilateral hearing loss and hearing aid use in the United States. Laryngoscope 128, 1681–1686. doi: 10.1002/lary.27017
Grahn, J., Parkinson, J., and Owen, A. (2009). The role of the basal ganglia in learning and memory: neuropsychological studies. Behav. Brain Res. 199, 53–60. doi: 10.1016/j.bbr.2008.11.020
Gurgel, R., Ward, P., Schwartz, S., Norton, M., Foster, N., and Tschanz, J. (2014). Relationship of hearing loss and dementia: a prospective, population-based study. Otology & neurotology : official publication of the American Otological Society, American Neurotology Society [and] European academy of. Otol. Neurotol. 35, 775–781. doi: 10.1097/MAO.0000000000000313
Jung, M., Colletta, M., Coalson, R., Schlaggar, B., and Lieu, J. (2017). Differences in interregional brain connectivity in children with unilateral hearing loss. Laryngoscope 127, 2636–2645. doi: 10.1002/lary.26587
Koos, W, Spetzler, R, and Lang, J. Color atlas of microneurosurgery : Stuttgart:Thieme Germany, 1993.
Lieu, J., Tye-Murray, N., and Fu, Q. (2012). Longitudinal study of children with unilateral hearing loss. Laryngoscope 122, 2088–2095. doi: 10.1002/lary.23454
Lin, F., Ferrucci, L., Metter, E., An, Y., Zonderman, A., and Resnick, S. (2011). Hearing loss and cognition in the Baltimore longitudinal study of aging. Neuropsychology 25, 763–770. doi: 10.1037/a0024238
Lisman, J., Buzsáki, G., Eichenbaum, H., Nadel, L., Ranganath, C., and Redish, A. (2017). Viewpoints: how the hippocampus contributes to memory, navigation and cognition. Nat. Neurosci. 20, 1434–1447. doi: 10.1038/nn.4661
Ma, W., Li, M., Gao, F., Zhang, X., Shi, L., Yu, L., et al. (2016). DTI analysis of Presbycusis using voxel-based analysis. Am. J. Neuroradiol. 37, 2110–2114. doi: 10.3174/ajnr.A4870
Mori, S., Crain, B., Chacko, V., and van Zijl, P. (1999). Three-dimensional tracking of axonal projections in the brain by magnetic resonance imaging. Ann. Neurol. 45, 265–269. doi: 10.1002/1531-8249(199902)45:2<265::AID-ANA21>3.0.CO;2-3
Nelson, E., Reeder, R., Holden, L., and Firszt, J. (2019). Front- and rear-facing horizontal sound localization results in adults with unilateral hearing loss and normal hearing. Hear. Res. 372, 3–9. doi: 10.1016/j.heares.2018.03.016
Newman, C., Jacobson, G., and Spitzer, J. (1996). Development of the tinnitus handicap inventory. Arch. Otolaryngol. Head Neck Surg. 122, 143–148.
Nyberg, L., Cabeza, R., and Tulving, E. (1996). PET studies of encoding and retrieval: the HERA model. Psychon. Bull. Rev. 3, 135–148. doi: 10.3758/BF03212412
Propst, E., Greinwald, J., and Schmithorst, V. (2010). Neuroanatomic differences in children with unilateral sensorineural hearing loss detected using functional magnetic resonance imaging. Arch. Otolaryngol. Head Neck Surg. 136, 22–26. doi: 10.1001/archoto.2009.208
Pulvermüller, F., Huss, M., Kherif, F., Moscoso del Prado Martin, F., Hauk, O., and Shtyrov, Y. (2006). Motor cortex maps articulatory features of speech sounds. Proc. Natl. Acad. Sci. U. S. A. 103, 7865–7870. doi: 10.1073/pnas.0509989103
Seeley, W., Menon, V., Schatzberg, A., Keller, J., Glover, G. H., Kenna, H., et al. (2007). Dissociable intrinsic connectivity networks for salience processing and executive control. J. Neurosci. 27, 2349–2356. doi: 10.1523/JNEUROSCI.5587-06.2007
Sekeres, M., Winocur, G., and Moscovitch, M. (2018). The hippocampus and related neocortical structures in memory transformation. Neurosci. Lett. 680, 39–53. doi: 10.1016/j.neulet.2018.05.006
Shtyrov, Y., Butorina, A., Nikolaeva, A., and Stroganova, T. (2014). Automatic ultrarapid activation and inhibition of cortical motor systems in spoken word comprehension. Proc. Natl. Acad. Sci. U. S. A. 111, E1918–E1923. doi: 10.1073/pnas.1323158111
Shu, N., Li, X., Ma, C., Zhang, J., Chen, K., Liang, Y., et al. (2015). Effects of APOE promoter polymorphism on the topological organization of brain structural connectome in nondemented elderly. Hum. Brain Mapp. 36, 4847–4858. doi: 10.1002/hbm.22954
Sidlauskaite, J., Sonuga-Barke, E., Roeyers, H., and Wiersema, J. (2016). Altered intrinsic organisation of brain networks implicated in attentional processes in adult attention-deficit/hyperactivity disorder: a resting-state study of attention, default mode and salience network connectivity. Eur. Arch. Psychiatry Clin. Neurosci. 266, 349–357. doi: 10.1007/s00406-015-0630-0
Smith, E, Jonides, J, and Koeppe, R. Dissociating verbal and spatial working memory using PET. Cerebral cortex (New York, NY: 1991) 1996:11–20.
Speer, N, Zacks, J, and Reynolds, J. (2005). “Neural activity during reading reflects changes in the situation described by the text.” in Proceedings of the Society for Neuroscience Annual Meeting. Washington, DC.
Sporns, O. (2011). The human connectome: a complex network. Ann. N. Y. Acad. Sci. 1224, 109–125. doi: 10.1111/j.1749-6632.2010.05888.x
Su, C., Chen, H., Kwan, A., Lin, Y., and Guo, N. (2007). Neuropsychological impairment after hemorrhagic stroke in basal ganglia. Arch. Clin. Neuropsychol. 22, 465–474. doi: 10.1016/j.acn.2007.01.025
Su, T., Hsu, T., Lin, Y., and Lin, C. (2015). Schizophrenia symptoms and brain network efficiency: a resting-state fMRI study. Psychiatry Res. 234, 208–218. doi: 10.1016/j.pscychresns.2015.09.013
Suzuki, M., Hashimoto, S., Kano, S., and Okitsu, T. (2010). Prevalence of acoustic neuroma associated with each configuration of pure tone audiogram in patients with asymmetric sensorineural hearing loss. Ann. Otol. Rhinol. Laryngol. 119, 615–618. doi: 10.1177/000348941011900908
Tononi, G., Sporns, O., and Edelman, G. (1994). A measure for brain complexity: relating functional segregation and integration in the nervous system. Proc. Natl. Acad. Sci. U. S. A. 91, 5033–5037. doi: 10.1073/pnas.91.11.5033
Tulving, E., Kapur, S., Craik, F., Moscovitch, M., and Houle, S. (1994). Hemispheric encoding/retrieval asymmetry in episodic memory: positron emission tomography findings. Proc. Natl. Acad. Sci. U. S. A. 91, 2016–2020. doi: 10.1073/pnas.91.6.2016
van Wieringen, A., Boudewyns, A., Sangen, A., Wouters, J., and Desloovere, C. (2019). Unilateral congenital hearing loss in children: challenges and potentials. Hear. Res. 372, 29–41. doi: 10.1016/j.heares.2018.01.010
Vincent, C., Arndt, S., Firszt, J., Fraysse, B., Kitterick, P. T., Papsin, B. C., et al. (2015). Identification and evaluation of cochlear implant candidates with asymmetrical hearing loss. Audiol. Neurootol. 20, 87–89. doi: 10.1159/000380754
Wang, X., Fan, Y., Zhao, F., Wang, Z., Ge, J., Zhang, K., et al. (2014). Altered regional and circuit resting-state activity associated with unilateral hearing loss. PLoS One 9:e96126. doi: 10.1371/journal.pone.0096126
Wang, J., Wang, X., He, Y., Yu, X., Wang, H., and He, Y. (2015a). Apolipoprotein E ε4 modulates functional brain connectome in Alzheimer's disease. Hum. Brain Mapp. 36, 1828–1846. doi: 10.1002/hbm.22740
Wang, J., Wang, X., Xia, M., Liao, X., Evans, A., and He, Y. (2015b). GRETNA: a graph theoretical network analysis toolbox for imaging connectomics. Front. Hum. Neurosci. 9:386. doi: 10.3389/fnhum.2015.00386
Wang, J., Zuo, X., Dai, Z., Xia, M., Zhao, Z., Zhao, X., et al. (2013). Disrupted functional brain connectome in individuals at risk for Alzheimer's disease. Biol. Psychiatry 73, 472–481. doi: 10.1016/j.biopsych.2012.03.026
Zalesky, A., Fornito, A., and Bullmore, E. (2010). Network-based statistic: identifying differences in brain networks. NeuroImage 53, 1197–1207. doi: 10.1016/j.neuroimage.2010.06.041
Zhang, Y., Mao, Z., Feng, S., Wang, W., Zhang, J., and Yu, X. (2018). Convergent and divergent functional connectivity patterns in patients with long-term left-sided and right-sided deafness. Neurosci. Lett. 665, 74–79. doi: 10.1016/j.neulet.2017.11.050
Keywords: cognition, acoustic neuroma, vestibular schwannoma, diffusion tensor imaging, brain network
Citation: Deng X, Liu L, Li J, Yao H, He S, Guo Z, Sun J, Liu W and Hui X (2022) Brain structural network to investigate the mechanism of cognitive impairment in patients with acoustic neuroma. Front. Aging Neurosci. 14:970159. doi: 10.3389/fnagi.2022.970159
Edited by:
Kimberly Diggle Mueller, University of Wisconsin-Madison, United StatesReviewed by:
Jun Yang, Shanghai Jiaotong University School of Medicine, ChinaRandall Porter, Barrow Neurological Institute, United States
Copyright © 2022 Deng, Liu, Li, Yao, He, Guo, Sun, Liu and Hui. 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: Xuhui Hui, huixuhuiwch@126.com
†These authors have contributed equally to this work and share first authorship