Skip to main content

ORIGINAL RESEARCH article

Front. Psychiatry, 17 June 2022
Sec. Neuroimaging
This article is part of the Research Topic Temporal Lobe Dysfunction in Neuropsychiatric Disorder View all 22 articles

Abnormal Ventral Somatomotor Network Homogeneity in Patients With Temporal Lobe Epilepsy

\r\nDongbin Li,Dongbin Li1,2Ruoshi LiuRuoshi Liu3Lili Meng,Lili Meng4,5Pingan XiongPingan Xiong6Hongwei RenHongwei Ren7Liming Zhang*Liming Zhang1*Yujun Gao*Yujun Gao8*
  • 1Department of Neurology, The First Affiliated Hospital of Harbin Medical University, Harbin, China
  • 2First Department of Neurology and Neuroscience Center, Heilongjiang Provincial Hospital, Harbin, China
  • 3Department of Neurology, The Fourth Affiliated Hospital of Harbin Medical University, Harbin, China
  • 4Department of Psychiatry, Wuhan Mental Health Center, Wuhan, China
  • 5Department of Sleep, Wuhan Hospital for Psychotherapy, Wuhan, China
  • 6Department of Taihe Hospital Reproductive Medicine Center Affiliated To Hubei University of Medicine, Shiyan, China
  • 7Department of Medical Imaging, Tianyou Hospital Affiliated To Wuhan University of Science and Technology, Wuhan, China
  • 8Department of Psychiatry, Renmin Hospital of Wuhan University, Wuhan, China

Background: Abnormalities of functional connectivity in the somatomotor network have been thought to play an essential role in the pathophysiology of epilepsy. However, there has been no network homogeneity (NH) study about the ventral somatomotor network (VSN) in patients with temporal lobe epilepsy (TLE). Therefore, we explored the NH of the VSN in TLE patients in this study.

Methods: The sample included 52 patients with left temporal lobe epilepsy, 83 patients with right temporal lobe epilepsy, and 68 healthy controls. The NH method was utilized to analyze the resting-state functional magnetic resonance imaging data.

Results: Compared to the controls, rTLE patients had significantly higher NH in the bilateral postcentral gyrus, and significantly lower NH in the bilateral Rolandic operculum and the right superior temporal gyrus (STG). The NH values of the left postcentral gyrus were significantly higher in lTLE patients than in the healthy controls, and lTLE patients had lower NH in the right Rolandic operculum. The altered NH in the postcentral gyrus was negatively correlated with the illness duration, and the decreased NH in the left Rolandic operculum was negatively correlated with the executive control reaction time (ECRT).

Conclusion: Our findings suggest that altered NH of the postcentral gyrus, Rolandic operculum and STG might be associated with the pathophysiology of TLE, and thus, highlight the contribution of the VSN to the pathophysiology of TLE.

Introduction

Epilepsy is one of the most common neurological disorders, affecting over 70 million people worldwide and putting a considerable strain on health-care infrastructure and the economy (1). Temporal lobe epilepsy (TLE) is the most common type of partial epilepsy referred for surgery, accounting for more than 40% of surgical cases, because of the failure of antiepileptic drug therapy (2). TLE is characterized by cognitive dysfunction, which includes memory, executive functioning, and general intellectual functioning disorders, thus contributing to a poor quality of life (3). Although temporal lobectomy results in seizure-free status in 70% of TLE patients (4), surgical resection is still invasive and has undesirable side effects. Recent work has linked TLE with network-level disruption (5), and resting-state functional magnetic resonance imaging (rs-fMRI) studies have revealed its neural connections. MRI is thought to be the most direct method for demonstrating functional connectivity among separated regions of the brain at rest and is the primary method used for studying brain networks.

The somatomotor network is a resting-state network composed of bilateral pre- and postcentral gyri. Based on large amounts of resting-state fMRI data and a data-driven clustering approach (6), Biswal et al. (7) found that the somatomotor network belongs to the conventional group of seven cortical neuronal networks. Several studies have confirmed that the somatomotor network has strong positive connectivity with various brain areas, including the ventral attention networks, frontoparietal networks, and default mode networks (8). In addition, the somatomotor network is closely correlated with the occurrence and progression of diseases such as autism spectrum disorder, schizophrenia, and major depressive disorder (810). However, there have been a few reports on the somatomotor network in patients with epilepsy. Furthermore, the symptoms of epilepsy, such as myotonia, myoclonia, and atonic seizures, are closely related to the motor system. In this study, we focused our work on the potential functional mechanism of the ventral somatomotor network (VSN), which is a subdivision of the somatomotor network, in patients with TLE.

With the development of imaging technology in recent years, many brain structural and functional differences have been reported in patients with TLE and healthy controls. For example, Gao et al. (11) found a decrease in functional connectivity and structural deficits in the alerting network of patients with right-sided TLE (rTLE) by using the seed-based functional connectivity method. Mankinen et al. (12) found that interictal epileptiform activity may lead to the reorganization of resting-state brain networks by using independent component analysis (ICA). The seed-based region-of-interest (ROI) functional connectivity method and ICA are two approaches that are often used to assess brain networks. However, these analytical methods have several disadvantages. For example, ROI seed-based methods are critical for selecting an a priori ROI within a network. Unlike the ROI seed-based method, ICA does not require an a priori definition of seed regions; however, the results are highly dependent on the number of components the algorithm is asked to produce (13). To address these issues, an unbiased method for assessing imaging data is critical.

Network homogeneity (NH) (14) is a survey method that has been widely used in investigating many diseases, such as attention-deficit/hyperactivity disorder and depression, and has also been used to assess connectivity among brain networks (14, 15). The NH method is a novel method that combines the advantages of ICA and ROI seed-based functional connectivity. Thus, it provides an unbiased, hypothesis-driven measure for evaluating a specific brain network associated with a pathophysiologic process or a disorder. The NH method investigates a given network without requiring prior knowledge of the location of network abnormalities. As a voxel-wise survey, NH can be utilized to assess the connectivity of a voxel with the other brain voxels in a predefined network. Homogeneity is defined as the average connectivity of a given voxel. However, VSN homogeneity in patients with TLE has not been reported. In the present study, we analyzed rs-fMRI data from patients with TLE to investigate abnormalities in the NH of the VSN and explore the potential mechanism of impaired somatomotor function in TLE.

Materials and Approaches

Subjects

Fifty-two patients with left temporal lobe epilepsy (lTLE) and 83 patients with right temporal lobe epilepsy (rTLE) were recruited from the Department of Neurology, Tianyou Hospital Affiliated to Wuhan University of Science and Technology. Sixty-eight healthy people were recruited from those who underwent a standard physical examination at the medical examination center of the Tianyou Hospital Affiliated to Wuhan University of Science and Technology. Patients were diagnosed with TLE based on the diagnostic manual from the International League Against Epilepsy (16). The inclusion criteria for TLE were as follows (epilepsy patients who met any two of the following symptoms): (1) the clinical onset of symptoms suggested epileptic foci in the temporal lobe, such as psychiatric symptoms, abnormal emotional experiences, automatisms, epigastric rising, or dystonic posturing of the limb; (2) the imaging results showed atrophy or sclerosis in the right/left temporal lobe; and (3) the interictal electroencephalographic traces suggested an abnormality in the right/left temporal lobe. The exclusion criteria for all subjects were as follows: (1) left-handed; (2) any lifetime psychiatric disorder; (3) history of serious medical diseases or other neurological illness; and (4) Mini-Mental State Examination scores (MMSE) < 24. All subjects gave written, informed consent before participating in the study. All participants were right-handed, and the groups were matched by age, education level, and sex ratio. Our study was performed according to the Declaration of Helsinki and approved by the Medical Ethics Committee of the Tianyou Hospital Affiliated to Wuhan University of Science and Technology.

Behavioral Paradigm

The executive function was assessed by the attentional network test (ANT) (17). The stimulus signals of ANT visually appear on a screen, and the subjects were required to correctly and quickly identify the orientation in which a central target arrow pointed. The reaction time (RT) of all the subjects were recorded, and the executive control reaction time (ECRT) was calculated by subtracting the consistent arrow direction RT from the inconsistent arrow direction RT (17). A longer ECRT indicated inferior executive control performance.

Scan Acquisition

Images were acquired with an Achieva 3T MRI scanner (Philips, Amsterdam, the Netherlands) for resting-state functional magnetic resonance imaging (rsfMRI). “Participants were instructed to lie down with their eyes closed and remain awake. A prototype quadrature birdcage head coil filled with foam padding was used to limit the head motion. The scanning parameters were as follows: ratio of repetition time to echo time (TR/TE) (2,000/30 ms), slice thickness (5 mm), pitch (1 mm), field of view (240 × 240 mm), and flip angle (90°). On the structural scan (T1-weighted), the following settings were used: spin-echo sequence, repetition time (TR) = 20 ms, echo time (TE) = 3.5 ms, slice thickness = 1 mm, and field of view (FOV) = 24 × 24 cm.” [excerpted from our previous study (15)].

Data Preprocessing

Imaging data of rs-fMRI were preprocessed using the DPARSF software (18) in MATLAB (Mathworks). “After signal stabilization, head motion and slice-timing correction were conducted (19, 20). The subjects had a maximal translation ≤ 2 mm in the x, y, or z direction and an angular rotation ≤ 2° on each axis. The functional images were normalized to the standard template in Montreal Neurological Institute (MNI) template and spatially resampled to a voxel size of 3 mm × 3 mm × 3 mm. The head motion parameters obtained by rigid body correction, the white matter signal, and the cerebrospinal fluid signal were removed from the images by linear regression. The signal was bandpass filtered (0.01–0.1 Hz) and linearly detrended to reduce high-frequency physiological noise and low-frequency drift. The global signal removal may introduce artifacts into the data and distort resting-state connectivity patterns. Furthermore, the regression of the global signal may significantly distort results when studying clinical populations. Therefore, the global signal was preserved (21, 22)” [excerpted from our previous study (15)].

Ventral Somatomotor Network Identification

The toolbox GIFT1 was used to pick out VSN as a mask from all participants through the group ICA method. “Three steps from the GIFT toolbox were used as following: data reduction, separation of independent components, and back rebuilding. On the consideration of every component, the voxel-wise one-sample t-test set a statistical map and a threshold. Based on Gaussian random field (GRF) theory, p < 0.01 represents a significant statistical modification of multiple comparisons. Voxel significance: p < 0.01, and cluster significance: p < 0.01). Masks were created for the VSN components. Finally, the masks were combined to generate a VSN mask utilized in the following NH analysis.” [excerpted from our previous study (15)].

Network Homogeneity Analysis

MATLAB was used for NH analysis (14). “For each patient, the correlation coefficients were obtained in a given voxel with all other voxels within the VSN mask. The mean correlation coefficient was defined as the homogeneity of the given voxel, and subsequently changed into z-value by using z-transformation to improve the normal distribution as described. The resultant values generated the NH map that finally underwent z-transformation for group comparison” [excerpted from our previous study (15)].

Statistical Analyses

Demographic information, including sex, age, education degree, and imaging data, were calculated between the patient and control groups. Non-parametric Kruskal-Wallis test was used to compare the distributions between multiple groups because not all samples were in compliance with normal distribution. Categorical data were analyzed with a chi-square test using the IBMSPSS Statistics 22.0 software. Analyses of covariance were executed to compare differences across the three groups on voxel-based VSN maps. Then, post-hoc t-tests were executed to compare VSN differences between every two groups. Sex, age, years of education, and head motion were applied as covariates in group comparisons to limit the possible effects of these components. The significance level was set at the corrected p < 0.01 for multiple comparisons using the Gaussian Random Field (GRF) theory (GRF-corrected, voxel significance: p < 0.001, cluster significance: p < 0.01). Correlations between clinical variables were analyzed using partial correlations with head motion as a covariate. Bonferroni correction for multiple comparisons was used in the correlation analysis.

Results

Demographics and Clinical Characteristics of Subjects

In this study, 135 TLE patients (52 lTLE patients and 83 rTLE patients) and 68 age- and sex-matched healthy controls were recruited for the study. The demographic and clinical characteristics of the study subjects are provided in Table 1. No significant differences were observed among the three groups regarding age, sex, and years of education. The patient group showed a longer executive control reaction time.

TABLE 1
www.frontiersin.org

Table 1. Characteristics of the participants.

Ventral Somatomotor Network Maps Determined by Group Independent Component Analysis

By employing ICA, the VSN masks were chosen from the control group. The parts involved in the VSN were the bilateral postcentral and Rolandic operculum (Figure 1). The VSN was used as a mask in the following NH analysis.

FIGURE 1
www.frontiersin.org

Figure 1. Ventral somatomotor network (left: axial, coronal, and sagittal images of VSN; right: the detailed axial images of VSN. Based on group ICA with a threshold at z ≥ 5).

Group Differences in Ventral Somatomotor Network Regarding Network Homogeneity

Significant group differences in NH values between the patients (rTLE/lTLE) and the controls within the VSN mask were observed via two-sample t-tests. Compared to the controls, rTLE patients had significantly higher NH in the bilateral postcentral gyrus, and significantly lower NH in the bilateral Rolandic operculum and right superior temporal gyrus (STG) (Figure 2). The NH values of the left postcentral gyrus were significantly higher in lTLE patients than in the healthy controls, and lTLE patients had lower NH in the right Rolandic operculum (Figure 3). The NH values of the right Rolandic operculum and left postcentral gyrus were significantly higher in rTLE patients than in the lTLE patients (Figure 4).

FIGURE 2
www.frontiersin.org

Figure 2. Statistical maps showing NH differences in the bilateral postcentral gyrus, the bilateral Rolandic operculum and the right STG between the rTLE group and control group. (Blue indicates lower NH and the color bar indicates the t-values from the two-sample t-test).

FIGURE 3
www.frontiersin.org

Figure 3. Statistical maps showing NH differences in the left postcentral gyrus and the right Rolandic operculum between the lTLE group and the control group. (Blue indicates lower NH and the color bar indicates the t-values from the two-sample t-test).

FIGURE 4
www.frontiersin.org

Figure 4. Statistical maps showing NH differences in the right Rolandic operculum and the left postcentral gyrus between the rTLE group and the lTLE group. (Blue indicates lower NH and the color bar indicates the t-values from the two-sample t-test).

Correlation of Network Homogeneity With Clinical Variables

The correlations between abnormal NH and clinical variables in the patients were examined. The increased NH of the left postcentral gyrus was negatively correlated with the illness duration in the lTLE group (r = –0.393, p = 0.004) (Figure 5A). The increased NH of the right postcentral gyrus was negatively correlated with the illness duration in the rTLE group (r = –0.345, p = 0.001) (Figure 5C). And the decreased NH in the left Rolandic operculum was negatively correlated with the ECRT in the rTLE group (r = –0.326, p = 0.003) (Figure 5B). No other correlations were observed in the participants.

FIGURE 5
www.frontiersin.org

Figure 5. Correlations between abnormal NH and clinical variables. (A) Negative correlation between the NH values in the left postcentral gyrus and illness duration in the lTLE group. (B) Negative correlation between the NH values in the left Rolandic operculum and ECRT in the rTLE group. (C) Negative correlation between the NH values in the right postcentral gyrus and illness duration in the rTLE group.

Discussion

Temporal lobe epilepsy is the most common drug-resistant form of epilepsy in adults (23) and is the most common indication for surgical intervention. However, surgical treatment is invasive and has several undesirable and serious side effects, and not every patient is a candidate for surgical intervention. Given the limited treatment options available for those with drug-resistant TLE, it is crucial to understand the pathophysiologic mechanism of TLE. Because TLE patients usually have motor and sensory impairment (24), we aimed to investigate the abnormal ventral somatomotor network in patients with rTLE/lTLE. In this study, we used an NH method to investigate the VSN region, which is associated with motor and somatosensory function, in patients with TLE (rTLE/lTLE). Compared to healthy controls, rTLE patients had significantly higher NH in the bilateral postcentral gyrus and significantly lower NH in the bilateral Rolandic operculum. Compared to healthy controls, lTLE patients had significantly higher NH in the left postcentral gyrus and lower NH in the right Rolandic operculum. These patients had a longer ECRT. In addition, the altered NH in the postcentral gyrus was negatively correlated with the illness duration, and the decreased NH in the left Rolandic operculum was negatively correlated with the ECRT. We speculated that abnormal NH in the postcentral gyrus and Rolandic operculum play a critical role in the pathophysiology of TLE.

The postcentral gyrus, which is located on the lateral surface of the parietal lobe between the central sulcus and postcentral sulcus and caudal to the central sulcus, corresponds to Brodmann areas 3b, 1, and 2 (25), which contain the primary somatosensory network and are thereby crucial for proprioception and motor control (26). Furthermore, the postcentral gyrus also includes the secondary somatosensory network, which is involved in integrating memories with somatosensory stimuli (27). Because of its unique characteristics, the postcentral gyrus is often affected in psychiatric illness. Li et al. (28) and Kilts et al. (29) found that increased activity of the postcentral gyrus is associated with decreased social anxiety symptoms. Zhuo et al. (30) reported decreased density in the resting-state global functional connectivity of the postcentral gyrus in patients with major depressive disorder. Larabi et al. (31) observed lowered activation of the postcentral gyrus during suppression, which resulted in poorer self-reflectiveness in schizophrenia patients. Jalbrzikowski et al. (32) confirmed that increased postcentral surface area is associated with less severe negative symptoms and better executive cognition in bipolar disorder. With a whole-brain voxel-based unbiased resting-state functional connectivity method, Cheng et al. (33) demonstrated reduced connectivity in the bilateral postcentral gyrus of patients with autism. Song et al. (34) demonstrated that destructive lesions in the postcentral gyrus, as well as changes in its outflow pathways to neighboring motor areas, could lead to epileptic negative myoclonus. These findings could be attributed to the location, structure and function of the postcentral gyrus. As reported by DiGuiseppi et al. (35), axons from the ventral posterolateral nucleus travel from the thalamus through the posterior limb of the internal capsule and terminate in the appropriate region of the postcentral gyrus. The postcentral gyrus has numerous connections with other brain areas, including the insula (36), amygdala (37), limbic system (38), cerebellum (39), and parietal lobe (40). Because of these numerous connections, the postcentral gyrus is able to perform a variety of functions, including those involved in somatic perceptual processes. The postcentral gyrus can also integrate sensory-motor connections associated with poor attentional set shifting and participate in suppression and cognitive insight because emotion-specific representations convey emotions and emotion processing. Our results showed significantly higher NH in the bilateral postcentral gyrus in patients with rTLE and higher NH in the left postcentral gyrus in patients with lTLE, which is consistent with previous studies (4143). Furthermore, the altered NH values of the left postcentral gyrus were negatively correlated with the duration of the illness in the lTLE group. And similar finding was observed in the rTLE group that altered NH values of the right postcentral gyrus were negatively correlated with the illness duration. We speculate that as rTLE/lTLE progressed, the correlation of the postcentral region with all other regions in the VSN decreased, resulting in functional loss in the postcentral region. We hypothesize that altered NH values in the postcentral gyrus may be responsible for cognitive dysfunctions in TLE patients, including executive cognition (32), attention alertness (44), emotion processing (29), and primary information processing.

The Rolandic operculum is defined as the confluence of the most caudal parts of the pre- and postcentral gyri, is contiguous to the oropharyngeal muscle control area of M1S1 and adjacent to the insula, and corresponds to Brodmann areas 6, 4, and 43 (45, 46). The most well-known function of the Rolandic operculum is its contribution to articulation and tongue movement during speech production, as it includes the ventral portion of the somatotopic tongue and lip representations (47). Furthermore, many other functions of the Rolandic operculum have been discovered in recent years. Fink et al. (48) found that an activated Rolandic operculum supports multimodal input processing from various motor, sensory, and perceptual sources. Ventre-Dominey (49) reported that the Rolandic operculum overlaps with a cortical network that is thought to be associated with self-referential processes that involve self-location in space. Blefari et al. (50) confirmed that the Rolandic operculum is crucial for interoceptive awareness and bodily self-consciousness by processing integrated exteroceptive-interoceptive signals. The findings of Wu et al. (51) suggest that functional impairments of the left Rolandic operculum in patients with schizophrenia are related to delusional thoughts, and patients with schizophrenia have been found to have increased mean connectivity in the left Rolandic operculum (52). Shan et al. (53) discovered that an activated Rolandic operculum is associated with cognitive disabilities in Alzheimer’s disease and mild cognitive impairment. Wang et al. (54) reported that the Rolandic operculum is associated with the neural mechanisms of tic generation. Decades ago, it was discovered that critical electrical discharges in temporal lobe epilepsy always affect extratemporal structures, such as the Rolandic operculum (55). In the present study, lower NH was found in the bilateral Rolandic operculum of patients with rTLE and in the right Rolandic operculum of patients with lTLE. Furthermore, we found that decreased NH in the left Rolandic operculum was negatively correlated with the ECRT between the control group and rTLE group. A longer ECRT indicates inferior executive control performance, and thus, we speculated that abnormal NH in the left Rolandic operculum played a critical role in the executive function of patients with rTLE.

In the present study, we observed that rTLE and lTLE patients had altered NH values in the postcentral gyrus and Rolandic operculum. Because the postcentral gyrus and Rolandic operculum in play critical roles in cognitive function, we speculate that network homogeneity abnormalities caused by temporal lobe epilepsy may account for cognitive dysfunctions in TLE patients. Additionally, this may be the cause of abnormal electrical discharge in patients with TLE. However, our results are consistent with those of some previous studies (56, 57). Li et al. (56) found that increased functional connectivity of the Rolandic operculum was lateralized in patients with rTLE, and Zhou et al. (57) reported that decreased functional connectivity of the postcentral gyrus was lateralized in patients with rTLE. However, we found that the altered NH of the Rolandic operculum and postcentral gyrus were bilateral in patients with rTLE. We speculate that the following factors may account for this phenomenon. (1) There might be a compensatory mechanism that occurs in patients with lTLE/rTLE at rest. As part of the VSN, the postcentral gyrus and the Rolandic operculum share many similarities and overlaps in function; for example, they are both closely related to depression, schizophrenia, epilepsy, and emotion and sensation processing. Bettus et al. (58) reported decreased basal functional connectivity within epileptogenic networks with concurrent contralaterally increased connectivity, which confirmed our suspicions. Campo et al. (59) discovered decreased effective connectivity in the medial temporal lobe of the lesional hemisphere and increased effective connectivity in the medial temporal lobe/inferior frontal cortex of the contralesional hemisphere. Similar to these findings regarding functional connectivity in patients with epilepsy, our results demonstrated increased NH values in the ipsilesional postcentral gyrus and decreased NH values in the contralesional Rolandic operculum. (2) Because TLE is a chronic disorder, epileptic discharges propagating into the contralateral gyrus and recurrent epileptic activities might be responsible for the contralateral structural and functional decline as TLE progressed. Jokeit et al. (60) confirmed the TLE patients had altered glucose metabolism in the contralateral hippocampus, which was related to seizure frequency. Additionally, the duration of TLE may have a greater impact on the contralateral measure of hippocampal metabolism than the ipsilateral measure. Following a multiple regression analysis, there was a subsequent decrease in the hippocampal volume contralateral to the primary temporal seizure focus (61, 62). We believe that similar processing occurs in the postcentral gyrus and Rolandic operculum, which is strongly supported by the results of our study.

The STG plays a key role in language function (63), and study showed that the activity of STG was strongly associated with the listeners’ subjective experience (64). A fMRI study enrolled 43 TLE patients with language-impaired and 42 TLE patients with non-language-impaired confirmed that activations within the STG was the most predictive of language impairment (65). A Chinese tasks based fMRI study showed that the main effect region of auditory naming and picture naming tasks was in right STG, and the main effect region of semantic fluency task was in left STG (66). Meanwhile, to protect the language function in the postoperative patients, the task-state fMRI and intraoperative electroencephalography were suggested to be used to develop a personalized surgical plan for epilepsy treatment (66). Moreover, the dysfunction in the STG was related to the cognitive dysfunction in epilepsy patients. A rs-fMRI study found that decreased functional connectivity of the left precentral gyrus with the bilateral STG has a significant effect on intelligence in children and adolescents with focal epilepsy (67). The right STG was found blunted neural response to emotional faces in pediatric epilepsy, which was regarded as a marker of risk for social cognitive deficits (68). In our study, lower NH in the right STG was found in the rTLE, and no significant correlation was found between the NH and clinical variables including the ECRT. Based on these prior studies, we speculate that the decreased NH was associated with the language function in TLE, though relevant language task hadn’t perform in our study.

There were several limitations in the study: (1) Physiological noise cannot be completely removed. (2) We focused only on abnormalities of the VSN in patients with TLE, and although illuminating the pathophysiological contribution of the VSN is critical, other significant brain networks might have been neglected. (3) Antiepileptic drugs may also have effects on functional networks in epilepsy, and we did not account for this confounding factor.

Conclusion

In conclusion, we used the NH method to analyze resting-state fMRI data in patients with TLE, and abnormal NH values in the VSN were confirmed. The importance of NH alterations in the VSN implies the importance of the postcentral gyrus, STG and Rolandic operculum in the progression of TLE and provides new insights into the pathophysiological mechanism of TLE.

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 the Medical Ethics Committee of the Tianyou Hospital Affiliated to Wuhan University of Science and Technology. 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

DL and LM: writing the article. PX, HR, and RL: collection and assembly of data. LZ and YG: research concept and design, data analysis, and interpretation. All authors contributed to the article and approved the submitted version.

Funding

The investigation was supported by the grant from the Health Commission of Hubei Province Scientific Research Project (Grant No. 2020CFB512) and the Epilepsy Research Fund of China Association Against Epilepsy (CW-2022-003).

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.

Footnotes

  1. ^ http://mialab.mrn.org/software/#gica

References

1. Thijs RD, Surges R, O’Brien TJ, Sander JW. Epilepsy in adults. Lancet. (2019) 393:689–701. doi: 10.1016/S0140-6736(18)32596-0

CrossRef Full Text | Google Scholar

2. Sheilabi MA, Takeshita LY, Sims EJ, Falciani F, Princivalle AP. The sodium channel B4-subunits are dysregulated in temporal lobe epilepsy drug-resistant patients. Int J Mol Sci. (2020) 21:2955. doi: 10.3390/ijms21082955

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Chen YY, Huang S, Wu WY, Liu CR, Yang XY, Zhao HT, et al. Associated and predictive factors of quality of life in patients with temporal lobe epilepsy. Epilepsy Behav. (2018) 86:85–90. doi: 10.1016/j.yebeh.2018.06.025

PubMed Abstract | CrossRef Full Text | Google Scholar

4. McIntosh AM, Wilson SJ, Berkovic SF. Seizure outcome after temporal lobectomy: current research practice and findings. Epilepsia. (2001) 42:1288–307. doi: 10.1046/j.1528-1157.2001.02001.x

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Bernhardt BC, Bonilha L, Gross DW. Network analysis for a network disorder: the emerging role of graph theory in the study of epilepsy. Epilepsy Behav. (2015) 50:162–70. doi: 10.1016/j.yebeh.2015.06.005

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Yeo BT, Krienen FM, Sepulcre J, Sabuncu MR, Lashkari D, Hollinshead M, et al. The organization of the human cerebral cortex estimated by intrinsic functional connectivity. J Neurophysiol. (2011) 106:1125–65. doi: 10.1152/jn.00338.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Biswal B, Yetkin FZ, Haughton VM, Hyde JS. Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn Reson Med. (1995) 34:537–41. doi: 10.1002/mrm.1910340409

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Korgaonkar MS, Goldstein-Piekarski AN, Fornito A, Williams LM. Intrinsic connectomes are a predictive biomarker of remission in major depressive disorder. Mol Psychiatry. (2020) 25:1537–49. doi: 10.1038/s41380-019-0574-2

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Lee M, Sehatpour P, Hoptman MJ, Lakatos P, Dias EC, Kantrowitz JT, et al. Neural mechanisms of mismatch negativity dysfunction in schizophrenia. Mol Psychiatry. (2017) 22:1585–93. doi: 10.1038/mp.2017.3

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Stoodley CJ. Distinct regions of the cerebellum show gray matter decreases in autism. ADHD, and developmental dyslexia. Front Syst Neurosci. (2014) 8:92. doi: 10.3389/fnsys.2014.00092

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Gao Y, Zheng J, Li Y, Guo D, Wang M, Cui X, et al. Decreased functional connectivity and structural deficit in alertness network with right-sided temporal lobe epilepsy. Medicine (Baltimore). (2018) 97:e0134. doi: 10.1097/MD.0000000000010134

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Mankinen K, Jalovaara P, Paakki JJ, Harila M, Rytky S, Tervonen O, et al. Connectivity disruptions in resting-state functional brain networks in children with temporal lobe epilepsy. Epilepsy Res. (2012) 100:168–78. doi: 10.1016/j.eplepsyres.2012.02.010

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Fox MD, Raichle ME. Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging. Nat Rev Neurosci. (2007) 8:700–11. doi: 10.1038/nrn2201

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Uddin LQ, Kelly AM, Biswal BB, Margulies DS, Shehzad Z, Shaw D, et al. Network homogeneity reveals decreased integrity of default-mode network in ADHD. J Neurosci Methods. (2008) 169:249–54. doi: 10.1016/j.jneumeth.2007.11.031

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Gao Y, Wang M, Yu R, Li Y, Yang Y, Cui X, et al. Abnormal default mode network homogeneity in treatment-naive patients with first-episode depression. Front Psychiatry. (2018) 9:697. doi: 10.3389/fpsyt.2018.00697

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Manford M, Fish DR, Shorvon SD. An analysis of clinical seizure patterns and their localizing value in frontal and temporal lobe epilepsies. Brain. (1996) 119:17–40. doi: 10.1093/brain/119.1.17

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Fan J, McCandliss BD, Sommer T, Raz A, Posner MI. Testing the efficiency and independence of attentional networks. J Cogn Neurosci. (2002) 14:340–7. doi: 10.1162/089892902317361886

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Yan CG, Zang YF. DPARSF: a MATLAB toolbox for “Pipeline” data analysis of resting-state fMRI. Front Syst Neurosci. (2010) 4:13. doi: 10.3389/fnsys.2010.00013

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Nichols TE, Das S, Eickhoff SB, Evans AC, Glatard T, Hanke M, et al. Best practices in data analysis and sharing in neuroimaging using MRI. Nat Neurosci. (2017) 20:299–303. doi: 10.1038/nn.4500

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Winkler AM, Ridgway GR, Douaud G, Nichols TE, Smith SM. Faster permutation inference in brain imaging. Neuroimage. (2016) 141:502–16. doi: 10.1016/j.neuroimage.2016.05.068

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Saad ZS, Gotts SJ, Murphy K, Chen G, Jo HJ, Martin A, et al. Trouble at rest: how correlation patterns and group differences become distorted after global signal regression. Brain Connect. (2012) 2:25–32. doi: 10.1089/brain.2012.0080

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Hahamy A, Calhoun V, Pearlson G, Harel M, Stern N, Attar F, et al. Save the global: global signal connectivity as a tool for studying clinical populations with functional magnetic resonance imaging. Brain Connect. (2014) 4:395–403. doi: 10.1089/brain.2014.0244

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Blümcke I, Thom M, Aronica E, Armstrong DD, Bartolomei F, Bernasconi A, et al. International consensus classification of hippocampal sclerosis in temporal lobe epilepsy: a task force report from the ILAE commission on diagnostic methods. Epilepsia. (2013) 54:1315–29. doi: 10.1111/epi.12220

PubMed Abstract | CrossRef Full Text | Google Scholar

24. McLachlan RS. A brief review of the anatomy and physiology of the limbic system. Can J Neurol Sci. (2009) 36:S84–7.

PubMed Abstract | Google Scholar

25. Akselrod M, Martuzzi R, Serino A, van der Zwaag W, Gassert R, Blanke O. Anatomical and functional properties of the foot and leg representation in areas 3b, 1 and 2 of primary somatosensory cortex in humans: a 7T fMRI study. Neuroimage. (2017) 159:473–87. doi: 10.1016/j.neuroimage.2017.06.021

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Borich MR, Brodie SM, Gray WA, Ionta S, Boyd LA. Understanding the role of the primary somatosensory cortex: opportunities for rehabilitation. Neuropsychologia. (2015) 79:246–55. doi: 10.1016/j.neuropsychologia.2015.07.007

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Chen TL, Babiloni C, Ferretti A, Perrucci MG, Romani GL, Rossini PM, et al. Human secondary somatosensory cortex is involved in the processing of somatosensory rare stimuli: an fMRI study. Neuroimage. (2008) 40:1765–71. doi: 10.1016/j.neuroimage.2008.01.020

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Li Y, Meng Y, Yuan M, Zhang Y, Ren Z, Zhang Y, et al. Therapy for adult social anxiety disorder: a meta-analysis of functional neuroimaging studies. J Clin Psychiatry. (2016) 77:e1429–38. doi: 10.4088/JCP.15r10226

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Kilts CD, Kelsey JE, Knight B, Ely TD, Bowman FD, Gross RE, et al. The neural correlates of social anxiety disorder and response to pharmacotherapy. Neuropsychopharmacology. (2006) 31:2243–53. doi: 10.1038/sj.npp.1301053

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Zhuo C, Zhu J, Wang C, Qu H, Ma X, Qin W. Different spatial patterns of brain atrophy and global functional connectivity impairments in major depressive disorder. Brain Imaging Behav. (2017) 11:1678–89. doi: 10.1007/s11682-016-9645-z

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Larabi DI, van der Meer L, Pijnenborg GHM, Ćurčić-Blake B, Aleman A. Insight and emotion regulation in schizophrenia: a brain activation and functional connectivity study. Neuroimage Clin. (2018) 20:762–71. doi: 10.1016/j.nicl.2018.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Jalbrzikowski M, Freedman D, Hegarty CE, Mennigen E, Karlsgodt KH, Olde Loohuis LM, et al. Structural brain alterations in youth with psychosis and bipolar spectrum symptoms. J Am Acad Child Adolesc Psychiatry. (2019) 58:1079–91. doi: 10.1016/j.jaac.2018.11.012

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Cheng W, Rolls ET, Gu H, Zhang J, Feng J. Autism: reduced connectivity between cortical areas involved in face expression, theory of mind, and the sense of self. Brain. (2015) 138:1382–93. doi: 10.1093/brain/awv051

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Song IU, Lee DG, Kim JS, An JY, Lee SB, Kim YI, et al. Unilateral epileptic negative myoclonus following focal lesion of the postcentral cerebral cortex due to acute middle cerebral infarction. J Clin Neurol. (2006) 2:272–5. doi: 10.3988/jcn.2006.2.4.272

PubMed Abstract | CrossRef Full Text | Google Scholar

35. DiGuiseppi J, Tadi P. Neuroanatomy, postcentral gyrus. StatPearls. (Treasure Island, FL: StatPearls Publishing) (2021).

PubMed Abstract | Google Scholar

36. Addicott MA, Luber B, Nguyen D, Palmer H, Lisanby SH, Appelbaum LG. Low– and high-frequency repetitive transcranial magnetic stimulation effects on resting-state functional connectivity between the postcentral gyrus and the insula. Brain Connect. (2019) 9:322–8. doi: 10.1089/brain.2018.0652

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Höistad M, Barbas H. Sequence of information processing for emotions through pathways linking temporal and insular cortices with the amygdala. Neuroimage. (2008) 40:1016–33. doi: 10.1016/j.neuroimage.2007.12.043

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Purves D, Augustine GJ, Fitzpatrick D, Hall WC, LaMantia AS, White LE. Neuroscience. 5th ed. Sunderland: Sinauer Associates (2012).

Google Scholar

39. Ramos TC, Balardin JB, Sato JR, Fujita A. Abnormal cortico-cerebellar functional connectivity in autism spectrum disorder. Front Syst Neurosci. (2019) 12:74. doi: 10.3389/fnsys.2018.00074

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Zhang L, Qiu F, Zhu H, Xiang M, Zhou L. Neural efficiency and acquired motor skills: an fmri study of expert athletes. Front Psychol. (2019) 10:2752. doi: 10.3389/fpsyg.2019.02752

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Jiang Y, Song L, Li X, Zhang Y, Chen Y, Jiang S, et al. Dysfunctional white-matter networks in medicated and unmedicated benign epilepsy with centrotemporal spikes. Hum Brain Mapp. (2019) 40:3113–24. doi: 10.1002/hbm.24584

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Guo L, Bai G, Zhang H, Lu D, Zheng J, Xu G. Cognitive functioning in temporal lobe epilepsy: a BOLD-fMRI study. Mol Neurobiol. (2017) 54:8361–9. doi: 10.1007/s12035-016-0298-0

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Zhu H, Zhu J, Bao FS, Liu H, Zhu X, Wu T, et al. Statistical parametric mapping for analyzing interictal magnetoencephalography in patients with left frontal lobe epilepsy. Seizure. (2016) 34:38–43. doi: 10.1016/j.seizure.2015.11.003

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Liu J, Zhou X, Zhang Z, Qin L, Ye W, Zheng J. Disrupted functional network in patients with temporal lobe epilepsy with impaired alertness. Epilepsy Behav. (2019) 101:106573. doi: 10.1016/j.yebeh.2019.106573

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Tzourio-Mazoyer N, Landeau B, Papathanassiou D, Crivello F, Etard O, Delcroix N, et al. Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage. (2002) 15:273–89. doi: 10.1006/nimg.2001.0978

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Brown S, Ngan E, Liotti M. A larynx area in the human motor cortex. Cereb Cortex. (2008) 18:837–45. doi: 10.1093/cercor/bhm131

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Brown S, Laird AR, Pfordresher PQ, Thelen SM, Turkeltaub P, Liotti M. The somatotopy of speech: phonation and articulation in the human motor cortex. Brain Cogn. (2009) 70:31–41. doi: 10.1016/j.bandc.2008.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Fink A, Bay JU, Koschutnig K, Prettenthaler K, Rominger C, Benedek M, et al. Brain and soccer: functional patterns of brain activity during the generation of creative moves in real soccer decision-making situations. Hum Brain Mapp. (2019) 40:755–64. doi: 10.1002/hbm.24408

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Ventre-Dominey J. Vestibular function in the temporal and parietal cortex: distinct velocity and inertial processing pathways. Front Integr Neurosci. (2014) 8:53. doi: 10.3389/fnint.2014.00053

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Blefari ML, Martuzzi R, Salomon R, Bello-Ruiz J, Herbelin B, Serino A, et al. Bilateral rolandic operculum processing underlying heartbeat awareness reflects changes in bodily self-consciousness. Eur J Neurosci. (2017) 45:1300–12. doi: 10.1111/ejn.13567

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Wu C, Zheng Y, Li J, She S, Peng H, Li L. Cortical gray matter loss, augmented vulnerability to speech-on-speech masking, and delusion in people with schizophrenia. Front Psychiatry. (2018) 9:287. doi: 10.3389/fpsyt.2018.00287

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Andreou C, Nolte G, Leicht G, Polomac N, Hanganu-Opatz IL, Lambert M, et al. Increased resting-state gamma-band connectivity in first-episode schizophrenia. Schizophr Bull. (2015) 41:930–9. doi: 10.1093/schbul/sbu121

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Shan Y, Wang JJ, Wang ZQ, Zhao ZL, Zhang M, Xu JY, et al. Neuronal specificity of acupuncture in Alzheimer’s disease and mild cognitive impairment patients: a functional MRI study. Evid Based Complement Alternat Med. (2018) 2018:7619197. doi: 10.1155/2018/7619197

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Wang Z, Maia TV, Marsh R, Colibazzi T, Gerber A, Peterson BS. The neural circuits that generate tics in Tourette’s syndrome. Am J Psychiatry. (2011) 168:1326–37. doi: 10.1176/appi.ajp.2011.09111692

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Bossi L, Munari C, Stoffels C, Bonis A, Bacia T, Talairach J, et al. Manifestations motrices et posturales dans les crises d’épilepsie d’origine temporale [Motor and postural manifestations of temporal lobe epilepsy seizure]. Rev Electroencephalogr Neurophysiol Clin. (1982) 12:101–6. doi: 10.1016/s0370-4475(82)80032-4

CrossRef Full Text | Google Scholar

56. Li J, Chen X, Ye W, Jiang W, Liu H, Zheng J. Alteration of the alertness-related network in patients with right temporal lobe epilepsy: a resting state fMRI study. Epilepsy Res. (2016) 127:252–9. doi: 10.1016/j.eplepsyres.2016.09.013

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Zhou X, Zhang Z, Liu J, Qin L, Pang X, Zheng J. Disruption and lateralization of cerebellar-cerebral functional networks in right temporal lobe epilepsy: a resting-state fMRI study. Epilepsy Behav. (2019) 96:80–6. doi: 10.1016/j.yebeh.2019.03.020

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Bettus G, Guedj E, Joyeux F, Confort-Gouny S, Soulier E, Laguitton V, et al. Decreased basal fMRI functional connectivity in epileptogenic networks and contralateral compensatory mechanisms. Hum Brain Mapp. (2009) 30:1580–91. doi: 10.1002/hbm.20625

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Campo P, Garrido MI, Moran RJ, García-Morales I, Poch C, Toledano R, et al. Network reconfiguration and working memory impairment in mesial temporal lobe epilepsy. Neuroimage. (2013) 72:48–54. doi: 10.1016/j.neuroimage.2013.01.036

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Jokeit H, Ebner A, Arnold S, Schüller M, Antke C, Huang Y, et al. Bilateral reductions of hippocampal volume, glucose metabolism, and wada hemispheric memory performance are related to the duration of mesial temporal lobe epilepsy. J Neurol. (1999) 246:926–33. doi: 10.1007/s004150050484

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Barr WB, Ashtari M, Schaul N. Bilateral reductions in hippocampal volume in adults with epilepsy and a history of febrile seizures. J Neurol Neurosurg Psychiatry. (1997) 63:461–7. doi: 10.1136/jnnp.63.4.461

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Salmenperä T, Kälviäinen R, Partanen K, Pitkänen A. Hippocampal damage caused by seizures in temporal lobe epilepsy. Lancet. (1998) 351:35. doi: 10.1016/S0140-6736(05)78092-2

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Bigler ED, Mortensen S, Neeley ES, Ozonoff S, Krasny L, Johnson M, et al. Superior temporal gyrus, language function, and autism. Dev Neuropsychol. (2007) 31:217–38. doi: 10.1080/87565640701190841

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Yi HG, Leonard MK, Chang EF. The encoding of speech sounds in the superior temporal gyrus. Neuron. (2019) 102:1096–110. doi: 10.1016/j.neuron.2019.04.023

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Kaestner E, Reyes A, Macari AC, Chang YH, Paul BM, Hermann BP, et al. Identifying the neural basis of a language-impaired phenotype of temporal lobe epilepsy. Epilepsia. (2019) 60:1627–38. doi: 10.1111/epi.16283

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Wang P, Du F, Li J, Yu H, Tang C, Jiang R. Functional magnetic resonance imaging based on Chinese tasks to protect language function in epileptics. Brain Behav. (2021) 11:e01979. doi: 10.1002/brb3.1979

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Songjiang L, Tijiang Z, Heng L, Wenjing Z, Bo T, Ganjun S, et al. Impact of brain functional network properties on intelligence in children and adolescents with focal epilepsy: a resting-state MRI study. Acad Radiol. (2021) 28:225–32. doi: 10.1016/j.acra.2020.01.004

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Morningstar M, Hung A, Grannis C, French RC, Mattson WI, Ostendorf AP, et al. Blunted neural response to emotional faces in the fusiform and superior temporal gyrus may be marker of emotion recognition deficits in pediatric epilepsy. Epilepsy Behav. (2020) 112:107432. doi: 10.1016/j.yebeh.2020.107432

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: temporal lobe epilepsy, ventral somatomotor network, network homogeneity, resting-state functional magnetic resonance imaging, executive function

Citation: Li D, Liu R, Meng L, Xiong P, Ren H, Zhang L and Gao Y (2022) Abnormal Ventral Somatomotor Network Homogeneity in Patients With Temporal Lobe Epilepsy. Front. Psychiatry 13:877956. doi: 10.3389/fpsyt.2022.877956

Received: 17 February 2022; Accepted: 16 May 2022;
Published: 17 June 2022.

Edited by:

Linda J. Larson-Prior, University of Arkansas for Medical Sciences, United States

Reviewed by:

Jinping Liu, Guilin Medical University, China
Jinou Zheng, Guangxi Medical University, China

Copyright © 2022 Li, Liu, Meng, Xiong, Ren, Zhang and Gao. 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: Liming Zhang, zfx001@yahoo.com; Yujun Gao, gaoyujun19820214@163.com

Disclaimer: 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.