Research Paper Volume 11, Issue 22 pp 10684—10696

Disrupted brain functional hub and causal connectivity in acute mild traumatic brain injury

Fengfang Li1, *, , Liyan Lu1, *, , Huiyou Chen1, , Peng Wang1, , Yu-Chen Chen1, , Hong Zhang2, , Xindao Yin1, ,

  • 1 Department of Radiology, Nanjing First Hospital, Nanjing Medical University, Nanjing, China
  • 2 Department of Radiology, The Affiliated Jiangning Hospital of Nanjing Medical University, Nanjing, China
* Equal contribution

Received: August 21, 2019       Accepted: November 8, 2019       Published: November 20, 2019      

https://doi.org/10.18632/aging.102484
How to Cite

Copyright © 2019 Li et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

There have been an increasing number of functional magnetic resonance imaging (fMRI) reports on brain abnormalities in mild traumatic brain injury (mTBI) at different phases. However, the neural bases and cognitive impairment after acute mTBI are unclear. This study aimed to identify brain functional hubs and connectivity abnormalities in acute mTBI patients and their correlations with deficits in cognitive performance. Within seven days after brain injury, mTBI patients (n=55) and age-, sex-, and educational -matched healthy controls (HCs) (n=41) underwent resting-state fMRI scans and cognitive assessments. We derived functional connectivity (FC) strength of the whole-brain network using degree centrality (DC) and performed Granger causality analysis (GCA) to analyze causal connectivity patterns in acute mTBI. Compared with HCs, acute mTBI patients had significantly decreased network centrality in the left middle frontal gyrus (MFG). Additionally, acute mTBI showed decreased inflows from the left MFG to bilateral middle temporal gyrus (MTG), left medial superior frontal gyrus (mSFG), and left anterior cingulate cortex (ACC). Correlation analyses revealed that changes in network centrality and causal connectivity were associated with deficits in cognitive performance in mTBI. Our findings may help to provide a new perspective for understanding the neuropathophysiological mechanism of acute cognitive impairment after mTBI.

Introduction

Traumatic brain injury (TBI) is a traumatic brain structural injury or an injury in which other changes in brain function caused by external forces occur; TBI is a major reason for sustained disability and morbidity, both in military populations and civilians [1]. Mild TBI (mTBI) is the most common type of traumatic brain injury, representing 80%–90% of TBI cases [2]. Many mTBI patients may develop a wide range of cognitive deficits, such as deficits in executive function, attention, and working memory [3, 4]. Although most mTBI patients recover quickly, approximately 10%-15% of patients have disabling problems that persist for a long time [5]. However, to date, no obvious structural damage has been identified in mTBI patients when using traditional brain imaging techniques (i.e., computed tomography (CT) and magnetic resonance imaging (MRI)) [6]. Therefore, the cognitive problems of many mTBI patients may be underrated and thus have a substantial impact on the patients’ lives and social interactions. Nonetheless, the pathophysiologic mechanisms of mTBI that cause the disorder to be persistent remain poorly understood.

Previous neuroimaging studies using positron emission tomography (PET), magnetoencephalography (MEG), and functional magnetic resonance imaging (fMRI) have investigated the neuropathophysiological mechanisms implicated in subacute and chronic mTBI patients [79]. Resting-state fMRI (rs-fMRI), which reflects brain activity while an individual is resting and not performing a specific task, may reveal damage to cognitive function caused by mTBI [10]. Using rs-fMRI, multiple resting-state brain intrinsic networks relevant to the neural mechanisms of mTBI have been demonstrated, such as the default mode network (DMN) [11, 12], executive control network (ECN) [13], and motor network [14]. However, these results were variable due to the application of different approaches to rs-fMRI data analysis and the selection of pre- and post-processing procedures and heterogeneous samples. Mayer et al. implied that there was decreased functional connectivity (FC) within the DMN in chronic mTBI patients using seed-based rs-fMRI analysis [15]. Through independent component analysis (ICA), Palacios et al. found that mTBI patients in the subacute stage showed reduced connectivity in the ECN, the frontal nodes of the DMN and the parietal areas of the dorsal attentional network, as well as an increase in connectivity in the visual network [10]. Through the fractional amplitude of low-frequency fluctuations (fALFF), Madhavan et al. observed decreased fALFF in the motor, language, visual and salience networks for patients with higher symptom severity scores [16]. Heterogeneity among mTBI populations may lead to changes in network status. Based on these latest theories and observations, we speculate that the pathophysiological basis of acute mTBI may be due to disrupted brain function and network hubs.

The latest rs-fMRI studies identified abnormalities in the prefrontal cortex (PFC) associated with mTBI-related cognitive deficits, such as deficits in executive function, attention and working memory [1720]. The PFC is a collection of interconnected new cortical regions that sends and receives projections from cortical motor systems, sensory systems, and many subcortical structures. The PFC helps to improve cognitive function and has been shown to be an important brain structure for performing executive functions [21]. Previous neuroimaging studies have confirmed the involvement of the PFC in the regulation and adaption of cognitive impairment after mTBI [22, 23]. Additionally considering the theory that the PFC is vulnerable to neural contusion, the PFC is considered a key region involved in mTBI. Therefore, exploring the role of the PFC in the brain functional network hubs of mTBI patients will provide valuable insights into the neural mechanisms of mTBI.

The ICA and seed-based FC approaches have been proved extremely useful in exploring FC patterns for specific components of interest. However, few studies have investigated the acute mTBI-related changes of whole-brain FC pattern or large-scale brain network. Degree centrality (DC) is a voxel-wise data-driven method for quantifying the importance of each node in a brain network [24]. This graph theory based on network analysis can evaluate the centrality of a network without selecting interested nodes or networks in advance. In addition, in order to examine the directional connectivity network involved in acute mTBI, Granger causality analysis (GCA) will be applied to evaluate the directional effects between simultaneously recorded time series. GCA allows the directional connections between different brain regions to be studied and the changes in brain networks to be dynamically observed, and GCA has been used to reveal causal relationships among brain regions in individuals with various diseases, such as autism, schizophrenia, Alzheimer’s disease (AD), and presbycusis [2528]. However, it remains to be determined whether such DC and causal connection abnormalities can be observed at the early stages of mTBI.

Thus, in the present study, in order to unravel the details of the network hub and causal connectivity of brain function in mTBI patients in the acute stage, we used DC to quantify the importance of each node in the brain network in patients with acute mTBI, based on the strength of its connection to all the other voxels and then used GCA to analyze the causal connections and understand the directional aspects of these changes. We further aimed to determine whether DC and causal connection abnormalities are associated with cognitive dysfunction performance in acute mTBI patients. We hypothesized that (1) the intrinsic disconnection pattern of PFC plays a pivotal role in the brain functional network of mTBI patients at the acute stage. (2) DC and causal connection abnormalities would be associated with cognitive dysfunction performance in mTBI patients at the acute stage.

Results

Demographics and neuropsychological characteristics

During this study, 95 patients with well-documented mTBI were identified. Among these patients, 70 patients were willing to participate in this study. Then, 6 patients were excluded because of excessive head motion artifacts, 5 patients were excluded because of a previous head injury, and 4 patients were excluded because of MRI contraindications. The remaining 55 mTBI patients obtained an MRI scan at an average of 3.47 days (range, 0–7 days) after head injury. A total of 41 healthy subjects matched by age, sex, and years of education were also assessed.

The basic demographic and cognitive performance characteristics of the mTBI group and HCs group are given in Table 1. There were no significant differences in age, sex or educational attainment between the mTBI group and the HCs group. Compared with the HCs, the mTBI patients at the acute stage performed worse on the MoCA scale (p=0.001). Among the subcategories of the MoCA test, our findings indicated that mTBI patients only showed significantly worse performance in visuospatial/executive (p=0.035) and attention (p=0.008) subcategories than did the HCs. There were no other significant differences in performance in naming, language, memory, abstraction or orientation between the mTBI group and HCs in this study.

Table 1. Demographic Characteristics and cognitive performance in patients with mTBI and HCs.

CharacteristicmTBI (n=55)Controls(n=41)p-value
Age(years)40.65±11.4443.24±11.720.280
Education(years)12.15±3.2713.22±3.450.124
Sex(Female/ Male)28/2722/190.838
GCS Score15--
Time since injury(d)3.47±1.89--
MoCA scores24.25±2.2625.88±2.400.001*
Visuospatial/executive3.69±0.844.10±1.020.035*
Naming2.85±0.412.73±0.450.163
Attention5.36±0.835.76±0.490.008*
Language2.35±0.652.49±0.590.272
Abstraction1.89±0.321.98±0.160.117
Memory2.31±1.172.68±1.650.197
Orientation5.80±0.405.90±0.300.175
Data are the mean ± standard deviation; *p<0.05. mTBI, mild traumatic brain injury; HCs, healthy controls; GCS, Glasgow Coma Scale; MoCA, Montreal Cognitive Assessment.

Degree centrality analysis

In both mTBI patients (Figure 1A) and HCs (Figure 1B), the spatial distribution of the weighted DC was mostly localized in the PFC and temporal lobe areas. Specifically, after two-sample t-test analysis, significantly decreased DC within the left middle frontal gyrus (MFG) was found in mTBI patients at the acute stage compared to that of the HCs (Figure 1C and Table 2).

(A) Spatial distribution of the degree centrality (DC) in mTBI patients at the acute stage. (B) Spatial distribution of the DC in the healthy controls (HCs). (C) Significantly decreased DC within the left middle frontal gyrus (MFG) between acute mTBI patients and HCs. mTBI, mild traumatic brain injury; HCs, healthy controls.

Figure 1. (A) Spatial distribution of the degree centrality (DC) in mTBI patients at the acute stage. (B) Spatial distribution of the DC in the healthy controls (HCs). (C) Significantly decreased DC within the left middle frontal gyrus (MFG) between acute mTBI patients and HCs. mTBI, mild traumatic brain injury; HCs, healthy controls.

Table 2. Regions showing significant differences in degree centrality and causal connectivity between mTBI patients and HCs.

Brain regionsBAPeak MNI coordinates x, y, z (mm)Peak T valueVoxels
Decreased DC
Left MFG10-33, 45, 18-4.126108
Decreased FC (left MFG → whole brain)
Left MTG20-57, -12, -24-3.457755
Left mSFG10-15, 63, 18-4.0393111
Left ACC32-6, 24, 30-4.6747113
Right MTG3954, -57, 3-3.6263157
Increased FC (whole brain → left MFG)
Left ACC32-12, 24, 304.131681
A corrected threshold of p<0.05 was determined by Monte Carlo simulation. DC, degree centrality; FC, functional connectivity; BA, Brodmann’s area; MNI, Montreal Neurological Institute; MFG, middle frontal gyrus; IFG, inferior frontal gyrus; MTG, middle temporal gyrus; mSFG, medial superior frontal gyrus; ACC, anterior cingulate cortex.

Causal connectivity analysis

The left MFG-to-whole brain GCA showed that patients with acute mTBI, compared with the HCs, demonstrated significantly decreased causal connectivity from the left MFG to the left middle temporal gyrus (MTG), right MTG, left medial superior frontal gyrus (mSFG), and left ACC. Next, whole brain-to-left MFG GCA showed that the acute mTBI group showed increased connectivity from the left ACC to the left MFG compared with the HC group, and we did not find a decreased feedback effect to the left MFG in the acute mTBI patients (Figure 2 and Table 2).

Aberrant causal connectivity from the bilateral SFG in mTBI patients at the acute stage. (A) Decreased causal connectivity from the left MFG to the L-MTG, R-MTG, L-ACC, and L-mSFG. (B) Increased causal connectivity from the left MFG to the left ACC. L, left; R, right; MFG, middle frontal gyrus; MTG, middle temporal gyrus; ACC, anterior cingulate cortex; mSFG, medial superior frontal gyrus.

Figure 2. Aberrant causal connectivity from the bilateral SFG in mTBI patients at the acute stage. (A) Decreased causal connectivity from the left MFG to the L-MTG, R-MTG, L-ACC, and L-mSFG. (B) Increased causal connectivity from the left MFG to the left ACC. L, left; R, right; MFG, middle frontal gyrus; MTG, middle temporal gyrus; ACC, anterior cingulate cortex; mSFG, medial superior frontal gyrus.

Correlation results

The decreased DC of the left MFG was positively correlated with the total MoCA score (r=0.426, p=0.002) in patients with mTBI in the acute stage, and another positive correlation between deficits in attention (r=0.346, p=0.012) and the DC of the left MFG was also revealed (Figure 3A). In addition, the visuospatial/executive score was negatively correlated with decreases in causal connectivity from the left MFG to the left ACC (r=-0.360, p=0.006) (Figure 3B). The other regions with abnormal causal connectivity revealed no significant correlations with the MoCA scores. However, no significance was observed after Bonferroni multiple comparisons correction due to the relatively strict calculation.

Correlations between the changes in degree centrality (DC) and causal connectivity and deficits in cognitive performance in mTBI patients at the acute stage. (A) The decreased DC in left MFG was positively correlated with the total MoCA scores (r = 0.426, p = 0.002) and the attention scores (r = 0.346, p = 0.012). (B) The visuospatial/executive score was negatively correlated with the decreased causal connectivity from the left MFG to left ACC (r = 0.360, p = 0.006). The full line represents the decreased causal connectivity and the dotted line represents the increased causal connectivity. L, left; R, right; MFG, middle frontal gyrus; MTG, middle temporal gyrus; ACC, anterior cingulate cortex; mSFG, medial superior frontal gyrus; MoCA, Montreal Cognitive Assessment.

Figure 3. Correlations between the changes in degree centrality (DC) and causal connectivity and deficits in cognitive performance in mTBI patients at the acute stage. (A) The decreased DC in left MFG was positively correlated with the total MoCA scores (r = 0.426, p = 0.002) and the attention scores (r = 0.346, p = 0.012). (B) The visuospatial/executive score was negatively correlated with the decreased causal connectivity from the left MFG to left ACC (r = 0.360, p = 0.006). The full line represents the decreased causal connectivity and the dotted line represents the increased causal connectivity. L, left; R, right; MFG, middle frontal gyrus; MTG, middle temporal gyrus; ACC, anterior cingulate cortex; mSFG, medial superior frontal gyrus; MoCA, Montreal Cognitive Assessment.

Discussion

This is the first study to use both DC and GCA approaches to explore intrinsic brain functional hubs and causal connectivity in mTBI patients at the acute stage. Using DC analysis, we found significantly decreased network centrality within the left MFG region in acute mTBI patients than in HCs. Using the GCA algorithm, we employed the left MFG as a seed to examine its causal effect with the whole brain. Specifically, the left MFG showed decreased causal connectivity with several regions, including the prefrontal, temporal and cingulate cortex. In addition, the left ACC showed significantly increased causal connectivity to the left MFG. Of note, the DC abnormalities in the left MFG and the disrupted causal connectivity between left ACC and left MFG were significantly correlated with cognitive function performance in mTBI patients at the acute stage. The current study could provide deeper insights into this issue and further confirm the concept that the acute mTBI is a pathophysiological process related to aberrant brain connectivity.

DC is the most direct measure to characterize node centrality in a network analysis. The higher the degree centrality a network node has, the more important the node is in the network. Numerous studies over the past century have confirmed the idea that the human brain network consists of neurons connected in complex patterns that support different cognition functions [29]. Our findings showed significantly decreased DC in the left MFG region of mTBI patients than in that of HCs, suggesting that the left MFG is the main brain functional hub affected by acute mTBI, which is consistent with the hypothesis that the dysconnectivity pattern of the PFC is involved in acute mTBI. However, only unilateral DC abnormalities were found in our study. A possible reason for this result is that the locations of the forces that caused head trauma in the patients included in this study were not completely consistent and need to be further explored. PFC is considered a key area by many TBI studies, which shows that the PFC is associated with many aspects of posttraumatic symptoms and cognitive impairment in individuals with mTBI [30, 22]. Harm et al. demonstrated significantly lower activation levels within the medial PFC in patients with mTBI in the subacute stage postinjury compared to HCs by using GLM analysis [17]. Moreover, Sullivan and colleagues found that chronic blast-related mTBI patients showed higher deactivation levels within the left dorsomedial PFC [13].

In the current study, we observed that decreased DC in the left MFG region was positively associated with the total MoCA score and attention score among patients with mTBI in the acute stage. These findings were consistent with the role of the PFC in the development of persistent cognitive and emotional impairment after a concussion. Consistent with this notion, previous studies have also reported functional alterations within the PFC region, both in moderate/severe TBI and in mild TBI, and such changes have been linked to cognitive performance [3133]. Moreover, there is also evidence of structural and neurochemical abnormalities in the PFG and a relationship between the structural and neurochemical changes and cognitive performance [21]. Dall et al. demonstrated increased thickness within the PFC regions and weak correlations between the amount of PFC thickening and cognitive performance in chronic mTBI patients [31]. Popesccu et al. indicated that mTBI patients showed reduced alpha-band power in the bilateral MFG by using resting-state MEG, and the alpha-band power in the MFG was negatively correlated with scores reflecting symptoms of emotional numbing [34]. Therefore, all related studies have suggested that the PFC, specifically the MFG, might be an important part of cognitive and emotional regulation circuits; changes in this region may be associated with characteristic symptoms of cognitive impairment. Meanwhile, consistent with these valuable findings, our results are instructive for understanding the specific role of MFG abnormalities in patients with acute mTBI and suggest that MFC may provide valuable insights into the neural mechanism of acute mTBI.

In the current study, we used the GCA method to study the directed functional connection of the brain network in mTBI patients at the acute stage. The left MFG-to-whole brain GCA showed that left MFG activity predicted subsequent decreases in activity in the frontal cortex (ACC and mSFG) and temporal cortex (MTG). Previous neuroimaging studies have found that abnormalities in the frontal cortex may be the cause of cognitive impairment after mTBI. First, it has been revealed that the ACC regulates activity in the executive areas of the brain that direct attention and produce responses [35]. Changes in the ACC have been found in the MEG signals between combat-related mTBI patients and HCs by using MEG and N-back working memory tasks [36]. Furthermore, we found a negative correlation between the visuospatial/executive score and decreases in connectivity from the MFG to the ACC, which also emphasized the pivotal role of the ACC in acute mTBI. However, the exact mechanism for this correlation result is unclear. We speculated that this negative correlation may also be a compensatory effect during the first few days following mTBI drawing on the previously study [37], but this compensatory mechanism needs to be confirmed in future studies. Our current results also indicated that acute mTBI patients showed decreased causal connectivity from the left MFG to the left mSFG. This finding was in line with that of another study with rs-fMRI, which showed significantly decreased ALFFs in the SFG of mTBI patients compared with those in the SFG of HCs [38]. Popescu and colleagues’ analysis of MEG data showed that the power of the alpha band in the SFG was generally lower in mTBI patients than in HCs, and the changes in alpha-band power were negatively associated with scores reflecting symptoms of emotional numbing [34]. All these studies supported the idea that the left MFG is an important part of information processing and can coordinate working memory (e.g., short-term memory) and other sources of information that are particularly important for attention and emotional maintenance [39].

Previous fMRI studies have suggested that MTG located in the posterior temporal lobe plays a vital role in auditory and motor processing functions [40]. However, recent neuropsychological and imaging literature suggested that MTG may be important for various aspects of working memory, visual and emotional cognition perception [41]. In the present study, we also found that the mTBI group showed reduced causal connectivity from the left MFG to both sides of the MTG. One previous study provided evidence that the mTBI group had significantly higher levels of fractional anisotropy and axial diffusivity in the MTG white matter region, and the diffusion alterations were correlated with acute symptom burden in mTBI patients [42]. By using brain SPECT imaging, Bonne and colleagues found that mTBI patients in the chronic stage had reduced relative blood perfusion in the MTG [43]. All these studies suggested that the structural and functional abnormalities in the MTG may be one of the underlying mechanisms of cognitive function impairment, and given the interactions between multiple functional networks that are known to exist in the brain, this change seems reasonable. Therefore, we speculated that the MTG may be a modal semantic center that binds information from different sensory modes together. Nevertheless, the association of the MTG with mTBI has not been substantially elucidated, and further research is needed to clarify the role of the MTG in mTBI.

Additionally, the whole brain-to-left MFG GCA also revealed that acute mTBI patients had significantly increased causal connectivity from the left ACC to the left MFG, suggesting that activity in the left ACC may predict subsequent increases in activity in the left MFG. Hyper-connectivity in TBI patients has previously been reported in task-related brain function and resting-state brain function. Chen et al. found that mTBI patients showed increased activity in several areas of the same region, including the pre-exercise area and the motor cortex, assuming that the activity in these areas represented a compensation mechanism because TBI patients performed as well on memory tasks as healthy controls [37]. Sullivan et al. showed enhanced negative coupling FC between the dorsomedial PFC and the dorsal ACC in blast-related mTBI patients, and the authors postulated that the enhanced negative coupling in mTBI patients serves as a compensatory mechanism for successful task performance [13]. In a recent study of cognitive control in civilians with mTBI, the enhanced activation in the mTBI group was also interpreted to reflect a compensatory mechanism, but this activation was observed in the parietal subcortex, which is not part of the functional network discussed in the current study [44]. The use of compensatory mechanisms that substitute for cognitive resources is also thought to explain functional imaging results in studies in moderate to severe TBI patients [45]. Consistent with the above notion, our results provide further evidence that compensatory activity may allow normal levels of some cognitive function performance to be maintained in acute mTBI patients.

There are several limitations in our study. First, our experimental design is cross-sectional and limited in sample size, so it is difficult to make direct causal inferences about the relationship between brain functional network structures and cognitive impairment characteristics in patients with mTBI. Additional longitudinal studies involving more subjects are required in the future. Second, our present study is limited to patients with mTBI in different populations, who showed different injury mechanisms and different brain injury sites. Heterogeneity of the injury mechanisms may minimize the generalization of observed results to any particular injury mechanism. Distinguishing the damage models and mechanisms of mTBI and maintaining homogeneity in the study population will allow better inferences to be made for different types of injuries. Third, the meaning of causal connectivity by GCA in the static state is not completely clear, and there are still some deficiencies. It is unclear whether low-pass filtering is the best method for resting fMRI data of GCA, and further studies are needed to clarify this meaningful view. Functional magnetic resonance imaging and electrophysiology require further research to elucidate the link between causal connections and neuronal activity. Finally, cognitive impairment after mTBI is very subtle, and our MoCA cognitive assessment method may not be sufficiently sensitive to identify these differences comprehensively. Therefore, in the future, it will be worthwhile to conduct additional studies using neuropsychological assessments and more sensitive cognitive performance indicators.

Conclusions

To summarize, we have observed crucial functional network hub and causal connectivity aberrations within several regions revealed by DC and GCA in patients with mTBI at the acute stage. Compared to HCs, acute mTBI patients showed decreased FC strength in the left MFG region and abnormal directionality of influence both from and to the left MFG. Furthermore, decreased DC values in the MFG and altered causal connectivity from the left MFG to the left ACC were correlated with deficits in cognitive performance. The abnormalities in brain functional network structures are helpful in better understanding the neuropathophysiological mechanism of cognitive impairment after mTBI at the acute stage, which may help diagnose and evaluate the degree of cognitive impairment and recovery of patients to a certain extent.

Materials and Methods

Subjects and clinical data

The current study was approved by the Institutional Review Board of Nanjing Medical University. Written informed consent was obtained from each subject before their participation in the study protocol.

We recruited 55 patients with mTBI (27 males and 28 females, 40.65 ± 11.44 years) at the acute stage from the emergency department of our hospital between January 2018 and May 2019. The inclusion criteria for patients were as follows: (1) age of 20 or older; (2) right-handedness before injury; (3) the presence of impact to the head; (4) an initial Glasgow Coma Score [GCS] of 13–15 in the emergency department; and (5) mTBI that was evaluated initially in the emergency room (posttraumatic amnesia <24 h, loss of consciousness<30 minutes, or recorded alteration of mental status (i.e., confused, dazed, or disoriented)). In addition, all patients were required to have a CT scan as part of their clinical evaluation. The exclusion criteria were as follows: (1) a documented history of a previous brain injury, the use of neuropsychological or psychoactive medications, a neurological disorder, or concurrent substance abuse; (2) history of drug or alcohol abuse; (3) history of the use of a sedative before a hospital or emergency department visit; and (4) MRI contraindications. Additionally, 41 age-, sex-, and educational attainment-matched healthy subjects (all right-handed) (19 males and 22 females; 43.24 ± 11.72 years) were recruited as controls. The control participants met the same exclusion criteria applied to the patient group. All control subjects underwent the same neuroimaging protocol as the patients. Demographic characteristics and outcome measurements of the mTBI and healthy control groups were summarized in Table 1.

Cognitive function assessment

Given the emergency care setting, a complete neuropsychological assessment was not feasible. Therefore, we used a short instrument called the Montreal Cognitive Assessment (MoCA) to assess all mTBI patients’ neurocognitive status. The MoCA was developed by Nasreddine et al. based on clinical experience and Mini-Mental State Examination cognitive ratings and items [46]. It is a brief yet comprehensive cognitive instrument used to assess the level of impairment in neurological populations. The test lasts approximately 10 minutes and has a maximum score of 30 points. A score greater than or equal to 26 indicates a cognitively normal state, and a lower score indicates poorer cognitive performance [46]. The MoCA evaluates several categories of functionality, including visuospatial/execution, language, attention, naming, abstraction, memory (short-term immediate and delayed recall), and location of time and place [47]. All subjects underwent the same cognitive function assessment within seven days after mTBI.

MRI acquisition

MRI data were acquired using a 3.0 T MRI scanner (Ingenia, Philips Medical Systems, Netherlands) at our hospital. The type of the receiver coil was digital head coil 3.0T and the parallel imaging was employed. Scanner noise and head motion were reduced using earplugs and foam padding. All subjects were asked to close their eyes, lie quietly without falling asleep, not think about anything in particular, and avoid any head movements during the scan. Structural images in this study were acquired using a three-dimensional turbo fast echo (3D-TFE) T1 WI sequence with high resolution parameters as follows: repetition time (TR) = 8.1 ms; echo time (TE) = 3.7 ms; thickness = 1 mm; slices = 170; gap = 0 mm; flip angle (FA)= 8°; FOV = 256 mm × 256 mm; and acquisition matrix = 256 mm × 256 mm. The structural sequence scan took 5 min and 28 s. Axial functional images were obtained using gradient echo-planer imaging (EPI) sequences, and the scanning parameters were as follows: slices = 36; thickness = 4 mm; TR/TE = 2000/30 ms; gap = 0 mm; acquisition matrix = 64 × 64; field of view (FOV) = 240 mm × 240 mm; and FA = 90°. The fMRI sequence scan took 8 min and 6 s.

Functional data preprocessing

Functional data analyses were conducted with the Data Processing & Analysis for Resting-State Brain Imaging (DPABI_V2.3_170105) [48] and statistical parametric mapping (SPM8, http://www.fil.ion.ucl.ac.uk/spm) packages. A total of 240 volumes were scanned, and the first 10 volumes were discarded to balance the initial MR signal and subject's adaptation to the scanner. Afterwards, the remaining 230 consecutive volumes were used for data analysis. Then, the following procedures were performed: slice-timing adjustment, head motion correction adjustment, Montreal Neurological Institute (MNI) template space normalization (resampling voxel size = 3 × 3 × 3 mm3), isotropic Gaussian kernel smoothing (full width at halt maximum (FWHM)= 6 mm), and detrending and filtering (0.01-0.08 Hz). Subjects with head motion > 2.0° of rotation or 2.0 mm of translation in any direction were excluded. Finally, several sources of spurious variances were removed by linear regression, which included six head motion parameters, and average signals from white matter, cerebrospinal fluid, and whole brain.

Degree centrality analysis

We limited voxel-wise centrality analyses to a predefined gray matter (GM) mask that included tissues with a GM probability >20%, as previously described [24]. In the mask, a separate network centrality map was generated in voxel-wise mode. First, whole brain correlation analysis based on the voxels was carried out for functional operation after preprocessing. The time course of each individual voxel of each participant was correlated with the time course of another individual voxel, resulting in a correlation matrix. Then, the undirected adjacency matrix was obtained by thresholding each association at r > 0.25 [49]. We chose a high threshold to prevent voxels with low temporal correlations due to signal noise from being counted. Then, we calculated the binary DC value of the whole brain network. A graph of the DC values for each gray matter voxel of each participant was obtained. Finally, these voxel-wise DC maps were converted into a z-score map for group comparisons. As previously reported, negative correlations were not included in DC calculations because their interpretation is not clear and has an adverse effect on the reliability of test-retest [50].

First, we estimated the spatial distribution of the mean DC in the mTBI group and the healthy group, respectively. The z-value of each individual was entered into the SPM8 software, and a random effects one-sample t-test was performed for each voxel to generate an average DC map for each group. The multiple comparison corrections were performed using false discovery rate (FDR) criterion and set at a threshold of p<0.01. In addition, to find the damaged brain center region, a two-tailed two-sample t-test was performed to study the differences in DC maps between mTBI patients and healthy controls (HCs). Between-group comparisons of the DC maps were performed in SPM8 software using general linear model (GLM) analysis, and age, sex, and educational attainment were included as unwanted covariates. Multiple comparison corrections were performed using the AlphaSim program determined by Monte Carlo simulation. Statistical maps of the two-sample t-test were created using a combined threshold of p<0.01 and a minimum cluster size of 40 voxels, yielding a corrected threshold of p<0.05. For between-group analyses, a mask was created by combining the significant clusters in the two groups, which were obtained from the one-sample t-test results.

Causal connectivity analysis

To further investigate the effects of FC directionality, we applied GCA to assess alterations in causal connectivity. Based on the results of the DC analysis, we selected a seed region that showed a significant difference between mTBI patients and HCs (left MFG: MNI coordinates (x, y, z) -33, 45, 18). The causal connectivity was analyzed by using REST-GCA in the REST toolbox [51]. In this study, the time series of the left MFG was defined as the seed time series x, and the time series y represented the time series of all voxels in the brain. The linear direct effect of x on y (Fx→y) and the linear direct effect of y on x (Fy→x) were voxels calculated by voxels in the brain. Therefore, two Granger causality maps were generated based on the impact metric for each participant. The residual-based F was normalized (F’) and normalized to the Z score for each voxel (Zy→x and Zx→y), subtracting the global mean F’ values, divided by standard deviation.

For group analyses of causal connectivity, the mean values of the Zy→x and Zx→y maps were calculated for each group. All four Granger causality maps, including two for each direction and two for each group, were obtained (the left MFG included both Z y→x and Z x→y for both the mTBI patients and HCs). These Granger causality maps were then entered into the SPM8 software for group comparisons. A random effects two-sample t-test was performed for each voxel to determine the difference in the causal connectivity of MFG between mTBI patients and HCs, with sex, age, and educational attainment included as unwanted covariates. Multiple comparison corrections were performed by Monte Carlo simulation using the AlphaSim program, with a corrected threshold of p < 0.01 and a minimum cluster size of 40 voxels, yielding a corrected threshold of p<0.05.

Correlation analysis

To explore the relationship between clinical cognitive deficit characteristics of mTBI patients and DC and causal connectivity measures, the areas showing significant differences in causal connections or DC between groups were extracted. We calculated the correlation between the mean z-values within these clusters and the cognitive deficit characteristics of each mTBI patient using Pearson correlation analysis in SPSS software (version 18.0; SPSS, Chicago, IL, USA). P < 0.05 was considered statistically significant and was corrected for age, sex, and educational attainment. Bonferroni correction for multiple comparisons was applied to the correlation analysis.

Author Contributions

Conceptualization, F.L. and L.L.; data curation, F.L. and L.L.; formal analysis, F.L. and L.L.; funding acquisition, Y.-C.C. and X.Y.; investigation, F.L., L.L., H.C., P.W. and H.Z.; methodology, F.L. and L.L.; project administration, Y.-C.C. and H.Z.; resources, Y.-C.C. and X.Y.; software, F.L. and L.L; supervision, Y.-C.C. and H.Z.; validation, F.L. and L.L.; visualization, H.C. and P.W.; writing—original draft, F.L. and L.L.; writing—review and editing, Y.-C.C. and H.Z.

Conflicts of Interest

There are no conflicts of interest to declare.

Funding

This work was funded by the Nanjing Medical University grant (No.2017NJMU124), Jiangsu Provincial Special Program of Medical Science (BE2017614), Youth Medical Talents of Jiangsu Province (No. QNRC2016062), 14th “Six Talent Peaks” Project of Jiangsu Province (No. YY-079), Nanjing Outstanding Youth Fund (No. JQX17006), and 333 High-level Talents Training Project of Jiangsu Province.

References

  • 1. Blennow K, Brody DL, Kochanek PM, Levin H, McKee A, Ribbers GM, Yaffe K, Zetterberg H. Traumatic brain injuries. Nat Rev Dis Primers. 2016; 2:16084. https://doi.org/10.1038/nrdp.2016.84 [PubMed]
  • 2. Eme R. Neurobehavioral Outcomes of Mild Traumatic Brain Injury: A Mini Review. Brain Sci. 2017; 7:7. https://doi.org/10.3390/brainsci7050046 [PubMed]
  • 3. de Freitas Cardoso MG, Faleiro RM, de Paula JJ, Kummer A, Caramelli P, Teixeira AL, de Souza LC, Miranda AS. Cognitive Impairment Following Acute Mild Traumatic Brain Injury. Front Neurol. 2019; 10:198. https://doi.org/10.3389/fneur.2019.00198 [PubMed]
  • 4. Karr JE, Areshenkoff CN, Garcia-Barrera MA. The neuropsychological outcomes of concussion: a systematic review of meta-analyses on the cognitive sequelae of mild traumatic brain injury. Neuropsychology. 2014; 28:321–36. https://doi.org/10.1037/neu0000037 [PubMed]
  • 5. Emery CA, Barlow KM, Brooks BL, Max JE, Villavicencio-Requis A, Gnanakumar V, Robertson HL, Schneider K, Yeates KO. A Systematic Review of Psychiatric, Psychological, and Behavioural Outcomes following Mild Traumatic Brain Injury in Children and Adolescents. Can J Psychiatry. 2016; 61:259–69. https://doi.org/10.1177/0706743716643741 [PubMed]
  • 6. Iraji A, Benson RR, Welch RD, O’Neil BJ, Woodard JL, Ayaz SI, Kulek A, Mika V, Medado P, Soltanian-Zadeh H, Liu T, Haacke EM, Kou Z. Resting State Functional Connectivity in Mild Traumatic Brain Injury at the Acute Stage: Independent Component and Seed-Based Analyses. J Neurotrauma. 2015; 32:1031–45. https://doi.org/10.1089/neu.2014.3610 [PubMed]
  • 7. Byrnes KR, Wilson CM, Brabazon F, von Leden R, Jurgens JS, Oakes TR, Selwyn RG. FDG-PET imaging in mild traumatic brain injury: a critical review. Front Neuroenergetics. 2014; 5:13. https://doi.org/10.3389/fnene.2013.00013 [PubMed]
  • 8. Li L, Arakaki X, Harrington M, Zouridakis G. Source Connectivity Analysis Can Assess Recovery of Acute Mild Traumatic Brain Injury Patients. Conf Proc IEEE Eng Med Biol Soc. 2018; 2018:3165–3168. https://doi.org/10.1109/EMBC.2018.8513045 [PubMed]
  • 9. Yamakawa GR, Lengkeek C, Salberg S, Spanswick SC, Mychasiuk R. Behavioral and pathophysiological outcomes associated with caffeine consumption and repetitive mild traumatic brain injury (RmTBI) in adolescent rats. PLoS One. 2017; 12:e0187218. https://doi.org/10.1371/journal.pone.0187218 [PubMed]
  • 10. Palacios EM, Yuh EL, Chang YS, Yue JK, Schnyer DM, Okonkwo DO, Valadka AB, Gordon WA, Maas AI, Vassar M, Manley GT, Mukherjee P. Resting-State Functional Connectivity Alterations Associated with Six-Month Outcomes in Mild Traumatic Brain Injury. J Neurotrauma. 2017; 34:1546–57. https://doi.org/10.1089/neu.2016.4752 [PubMed]
  • 11. Dretsch MN, Rangaprakash D, Katz JS, Daniel TA, Goodman AM, Denney TS, Deshpande G. Strength and Temporal Variance of the Default Mode Network to Investigate Chronic Mild Traumatic Brain Injury in Service Members with Psychological Trauma. J Exp Neurosci. 2019; 13:1179069519833966. https://doi.org/10.1177/1179069519833966 [PubMed]
  • 12. Zhu DC, Covassin T, Nogle S, Doyle S, Russell D, Pearson RL, Monroe J, Liszewski CM, DeMarco JK, Kaufman DI. A potential biomarker in sports-related concussion: brain functional connectivity alteration of the default-mode network measured with longitudinal resting-state fMRI over thirty days. J Neurotrauma. 2015; 32:327–41. https://doi.org/10.1089/neu.2014.3413 [PubMed]
  • 13. Sullivan DR, Hayes JP, Lafleche G, Salat DH, Verfaellie M. Functional Brain Alterations Associated With Cognitive Control in Blast-Related Mild Traumatic Brain Injury. J Int Neuropsychol Soc. 2018; 24:662–72. https://doi.org/10.1017/S1355617718000279 [PubMed]
  • 14. Wang S, Hu L, Cao J, Huang W, Sun C, Zheng D, Wang Z, Gan S, Niu X, Gu C, Bai G, Ye L, Zhang D, et al. Sex Differences in Abnormal Intrinsic Functional Connectivity After Acute Mild Traumatic Brain Injury. Front Neural Circuits. 2018; 12:107. https://doi.org/10.3389/fncir.2018.00107 [PubMed]
  • 15. Mayer AR, Mannell MV, Ling J, Gasparovic C, Yeo RA. Functional connectivity in mild traumatic brain injury. Hum Brain Mapp. 2011; 32:1825–35. https://doi.org/10.1002/hbm.21151 [PubMed]
  • 16. Madhavan R, Joel SE, Mullick R, Cogsil T, Niogi SN, Tsiouris AJ, Mukherjee P, Masdeu JC, Marinelli L, Shetty T. Longitudinal Resting State Functional Connectivity Predicts Clinical Outcome in Mild Traumatic Brain Injury. J Neurotrauma. 2019; 36:650–60. https://doi.org/10.1089/neu.2018.5739 [PubMed]
  • 17. Kaltiainen H, Liljeström M, Helle L, Salo A, Hietanen M, Renvall H, Forss N. Mild Traumatic Brain Injury Affects Cognitive Processing and Modifies Oscillatory Brain Activity during Attentional Tasks. J Neurotrauma. 2019; 36:2222–32. https://doi.org/10.1089/neu.2018.6306 [PubMed]
  • 18. Mayer AR, Yang Z, Yeo RA, Pena A, Ling JM, Mannell MV, Stippler M, Mojtahed K. A functional MRI study of multimodal selective attention following mild traumatic brain injury. Brain Imaging Behav. 2012; 6:343–54. https://doi.org/10.1007/s11682-012-9178-z [PubMed]
  • 19. Rockswold SB, Burton PC, Chang A, McNally N, Grant A, Rockswold GL, Low WC, Eberly LE, Yacoub E, Lenglet C. Functional Magnetic Resonance Imaging and Oculomotor Dysfunction in Mild Traumatic Brain Injury. J Neurotrauma. 2019; 36:1099–105. https://doi.org/10.1089/neu.2018.5796 [PubMed]
  • 20. Santhanam P, Wilson SH, Oakes TR, Weaver LK. Effects of mild traumatic brain injury and post-traumatic stress disorder on resting-state default mode network connectivity. Brain Res. 2019; 1711:77–82. https://doi.org/10.1016/j.brainres.2019.01.015 [PubMed]
  • 21. Funahashi S, Andreau JM. Prefrontal cortex and neural mechanisms of executive function. J Physiol Paris. 2013; 107:471–82. https://doi.org/10.1016/j.jphysparis.2013.05.001 [PubMed]
  • 22. van der Horn HJ, Liemburg EJ, Aleman A, Spikman JM, van der Naalt J. Brain Networks Subserving Emotion Regulation and Adaptation after Mild Traumatic Brain Injury. J Neurotrauma. 2016; 33:1–9. https://doi.org/10.1089/neu.2015.3905 [PubMed]
  • 23. Witt ST, Lovejoy DW, Pearlson GD, Stevens MC. Decreased prefrontal cortex activity in mild traumatic brain injury during performance of an auditory oddball task. Brain Imaging Behav. 2010; 4:232–47. https://doi.org/10.1007/s11682-010-9102-3 [PubMed]
  • 24. Hoptman MJ, Zuo XN, D'Angelo D, Mauro CJ, Butler PD, Milham MP, Javitt DC. Decreased interhemispheric coordination in schizophrenia: a resting state fMRI study. Schizophr Res. 2012; 141:1–7. https://doi.org/10.1016/j.schres.2012.07.027 [PubMed]
  • 25. Di Martino A, Zuo XN, Kelly C, Grzadzinski R, Mennes M, Schvarcz A, Rodman J, Lord C, Castellanos FX, Milham MP. Shared and distinct intrinsic functional network centrality in autism and attention-deficit/hyperactivity disorder. Biol Psychiatry. 2013; 74:623–32. https://doi.org/10.1016/j.biopsych.2013.02.011 [PubMed]
  • 26. Iwabuchi SJ, Palaniyappan L. Abnormalities in the effective connectivity of visuothalamic circuitry in schizophrenia. Psychol Med. 2017; 47:1300–10. https://doi.org/10.1017/S0033291716003469 [PubMed]
  • 27. Chen YC, Yong W, Xing C, Feng Y, Haidari NA, Xu JJ, Gu JP, Yin X, Wu Y. Directed functional connectivity of the hippocampus in patients with presbycusis. Brain Imaging Behav. 2019. [Epub ahead of print]. https://doi.org/10.1007/s11682-019-00162-z [PubMed]
  • 28. Liu Z, Zhang Y, Bai L, Yan H, Dai R, Zhong C, Wang H, Wei W, Xue T, Feng Y, You Y, Tian J. Investigation of the effective connectivity of resting state networks in Alzheimer’s disease: a functional MRI study combining independent components analysis and multivariate Granger causality analysis. NMR Biomed. 2012; 25:1311–20. https://doi.org/10.1002/nbm.2803 [PubMed]
  • 29. Medaglia JD. Functional Neuroimaging in Traumatic Brain Injury: From Nodes to Networks. Front Neurol. 2017; 8:407. https://doi.org/10.3389/fneur.2017.00407 [PubMed]
  • 30. Shin LM, Rauch SL, Pitman RK. Amygdala, medial prefrontal cortex, and hippocampal function in PTSD. Ann N Y Acad Sci. 2006; 1071:67–79. https://doi.org/10.1196/annals.1364.007 [PubMed]
  • 31. Dall’Acqua P, Johannes S, Mica L, Simmen HP, Glaab R, Fandino J, Schwendinger M, Meier C, Ulbrich EJ, Müller A, Jäncke L, Hänggi J. Prefrontal Cortical Thickening after Mild Traumatic Brain Injury: A One-Year Magnetic Resonance Imaging Study. J Neurotrauma. 2017; 34:3270–79. https://doi.org/10.1089/neu.2017.5124 [PubMed]
  • 32. Messé A, Caplain S, Pélégrini-Issac M, Blancho S, Lévy R, Aghakhani N, Montreuil M, Benali H, Lehéricy S. Specific and evolving resting-state network alterations in post-concussion syndrome following mild traumatic brain injury. PLoS One. 2013; 8:e65470. https://doi.org/10.1371/journal.pone.0065470 [PubMed]
  • 33. Zhou Y, Milham MP, Lui YW, Miles L, Reaume J, Sodickson DK, Grossman RI, Ge Y. Default-mode network disruption in mild traumatic brain injury. Radiology. 2012; 265:882–92. https://doi.org/10.1148/radiol.12120748 [PubMed]
  • 34. Popescu M, Hughes JD, Popescu EA, Riedy G, DeGraba TJ. Reduced prefrontal MEG alpha-band power in mild traumatic brain injury with associated posttraumatic stress disorder symptoms. Clin Neurophysiol. 2016; 127:3075–85. https://doi.org/10.1016/j.clinph.2016.06.004 [PubMed]
  • 35. Markela-Lerenc J, Ille N, Kaiser S, Fiedler P, Mundt C, Weisbrod M. Prefrontal-cingulate activation during executive control: which comes first? Brain Res Cogn Brain Res. 2004; 18:278–87. https://doi.org/10.1016/j.cogbrainres.2003.10.013 [PubMed]
  • 36. Huang MX, Nichols S, Robb-Swan A, Angeles-Quinto A, Harrington DL, Drake A, Huang CW, Song T, Diwakar M, Risbrough VB, Matthews S, Clifford R, Cheng CK, et al. MEG Working Memory N-Back Task Reveals Functional Deficits in Combat-Related Mild Traumatic Brain Injury. Cereb Cortex. 2019; 29:1953–1968. https://doi.org/10.1093/cercor/bhy075 [PubMed]
  • 37. Chen JK, Johnston KM, Petrides M, Ptito A. Neural substrates of symptoms of depression following concussion in male athletes with persisting postconcussion symptoms. Arch Gen Psychiatry. 2008; 65:81–89. https://doi.org/10.1001/archgenpsychiatry.2007.8 [PubMed]
  • 38. Xiong KL, Zhang JN, Zhang YL, Zhang Y, Chen H, Qiu MG. Brain functional connectivity and cognition in mild traumatic brain injury. Neuroradiology. 2016; 58:733–39. https://doi.org/10.1007/s00234-016-1675-0 [PubMed]
  • 39. Niogi SN, Mukherjee P, Ghajar J, Johnson CE, Kolster R, Lee H, Suh M, Zimmerman RD, Manley GT, McCandliss BD. Structural dissociation of attentional control and memory in adults with and without mild traumatic brain injury. Brain. 2008; 131:3209–21. https://doi.org/10.1093/brain/awn247 [PubMed]
  • 40. Clark RE. Current Topics Regarding the Function of the Medial Temporal Lobe Memory System. Curr Top Behav Neurosci. 2018; 37:13–42. https://doi.org/10.1007/7854_2017_481 [PubMed]
  • 41. Jeneson A, Squire LR. Working memory, long-term memory, and medial temporal lobe function. Learn Mem. 2011; 19:15–25. https://doi.org/10.1101/lm.024018.111 [PubMed]
  • 42. Babcock L, Yuan W, Leach J, Nash T, Wade S. White matter alterations in youth with acute mild traumatic brain injury. J Pediatr Rehabil Med. 2015; 8:285–96. https://doi.org/10.3233/PRM-150347 [PubMed]
  • 43. Bonne O, Gilboa A, Louzoun Y, Kempf-Sherf O, Katz M, Fishman Y, Ben-Nahum Z, Krausz Y, Bocher M, Lester H, Chisin R, Lerer B. Cerebral blood flow in chronic symptomatic mild traumatic brain injury. Psychiatry Res. 2003; 124:141–52. https://doi.org/10.1016/S0925-4927(03)00109-4 [PubMed]
  • 44. Mayer AR, Hanlon FM, Dodd AB, Ling JM, Klimaj SD, Meier TB. A functional magnetic resonance imaging study of cognitive control and neurosensory deficits in mild traumatic brain injury. Hum Brain Mapp. 2015; 36:4394–406. https://doi.org/10.1002/hbm.22930 [PubMed]
  • 45. Ricker JH, Hillary FG, DeLuca J. Functionally activated brain imaging (O-15 PET and fMRI) in the study of learning and memory after traumatic brain injury. J Head Trauma Rehabil. 2001; 16:191–205. https://doi.org/10.1097/00001199-200104000-00007 [PubMed]
  • 46. Nasreddine ZS, Phillips NA, Bédirian V, Charbonneau S, Whitehead V, Collin I, Cummings JL, Chertkow H. The Montreal Cognitive Assessment, MoCA: a brief screening tool for mild cognitive impairment. J Am Geriatr Soc. 2005; 53:695–99. https://doi.org/10.1111/j.1532-5415.2005.53221.x [PubMed]
  • 47. de Guise E, Alturki AY, LeBlanc J, Champoux MC, Couturier C, Lamoureux J, Desjardins M, Marcoux J, Maleki M, Feyz M. The Montreal Cognitive Assessment in persons with traumatic brain injury. Appl Neuropsychol Adult. 2014; 21:128–35. https://doi.org/10.1080/09084282.2013.778260 [PubMed]
  • 48. Yan CG, Wang XD, Zuo XN, Zang YF. DPABI: Data Processing & Analysis for (Resting-State) Brain Imaging. Neuroinformatics. 2016; 14:339–51. https://doi.org/10.1007/s12021-016-9299-4 [PubMed]
  • 49. Yan CG, Craddock RC, Zuo XN, Zang YF, Milham MP. Standardizing the intrinsic brain: towards robust measurement of inter-individual variation in 1000 functional connectomes. Neuroimage. 2013; 80:246–62. https://doi.org/10.1016/j.neuroimage.2013.04.081 [PubMed]
  • 50. Chen HJ, Jiang LF, Sun T, Liu J, Chen QF, Shi HB. Resting-state functional connectivity abnormalities correlate with psychometric hepatic encephalopathy score in cirrhosis. Eur J Radiol. 2015; 84:2287–95. https://doi.org/10.1016/j.ejrad.2015.08.005 [PubMed]
  • 51. Zang ZX, Yan CG, Dong ZY, Huang J, Zang YF. Granger causality analysis implementation on MATLAB: a graphic user interface toolkit for fMRI data processing. J Neurosci Methods. 2012; 203:418–26. https://doi.org/10.1016/j.jneumeth.2011.10.006 [PubMed]