Skip to main content


Front. Neurol., 09 June 2020
Sec. Applied Neuroimaging

Investigating Brain Network Changes and Their Association With Cognitive Recovery After Traumatic Brain Injury: A Longitudinal Analysis

  • 1CNNP Lab, Interdisciplinary Complex Systems Group, School of Computing, Newcastle University, Newcastle upon Tyne, United Kingdom
  • 2Faculty of Medical Sciences, Institute of Neuroscience, Newcastle University, Newcastle upon Tyne, United Kingdom
  • 3Department of Neurosurgery, Newcastle upon Tyne Hospitals NHS Foundation Trust, Newcastle upon Tyne, United Kingdom
  • 4Institute of Cellular Medicine, Newcastle MR Centre, Newcastle University, Newcastle upon Tyne, United Kingdom

Traumatic brain injury (TBI) can result in acute cognitive deficits and diffuse axonal injury reflected in white matter brain network alterations, which may, or may not, later recover. Our objective is to first characterize the ways in which brain networks change after TBI and, second, investigate if those changes are associated with recovery of cognitive deficits. We aim to make initial progress in discerning the relationships between brain network changes, and their (dys)functional correlates. We analyze longitudinally acquired MRI from 23 TBI patients (two time points: 6 days, 12 months post-injury) and cross-sectional data from 28 controls to construct white matter brain networks. Cognitive assessment was also performed. Graph theory and regression analysis were applied to identify changed brain network metrics after injury that are associated with subsequent improvements in cognitive function. Sixteen brain network metrics were found to be discriminative of different post-injury phases. Eleven of those explain 90% (adjusted R2) of the variability observed in cognitive recovery following TBI. Brain network metrics that had a high contribution to the explained variance were found in frontal and temporal cortex, additional to the anterior cingulate cortex. Our preliminary study suggests that network reorganization may be related to recovery of impaired cognitive function in the first year after a TBI.


Traumatic brain injury (TBI), which frequently involves white matter connectivity damage, is the leading cause of morbidity, death among children, and individuals under the age of 45 (1, 2). Every year in England and Wales, around 1.4 million patients attend hospital after sustaining a recent head injury, which represents 10% of all emergency admissions (3). Although, mild TBI (mTBI) patients usually make good recoveries, a significant proportion experience persistent cognitive deficits (4).

Diffuse axonal injury (DAI), one of the most common pathologies in TBI (5, 6), is triggered by mechanical disruption of axons, resulting in complex, and diverse effects on brain function

(7). Diffusion tensor imaging (DTI) is particularly suited to the study of DAI and has been used to investigate white matter brain connectivity changes after TBI (8, 9). To understand how widespread DAI lesions affect brain function, it may be necessary to analyze the global impact of these lesions on the whole-brain network (9).

Interpreting widespread network changes is challenging. Network neuroscience techniques (e.g., graph theory) allow the description and analysis of human brain network properties and have successfully been applied in cognitive neuroscience (911). Those techniques used anatomically defined gray matter regions as nodes and a measure of association (e.g., number of tracts, mean diffusion along tracts, etc.) between pairs of nodes based on the white matter tracts or edges connecting them (12). Those tracts can be estimated from tractography techniques or extracted from white matter atlases. All pairwise associations between regions are compiled in a connectivity matrix. The approach offers alternative analysis of how brain network changes may lead to cognitive variability over time.

Despite the importance of understanding temporal variability, most TBI studies are cross-sectional in design, which has limited value in understanding longitudinal recovery (13). Since TBI is a heterogeneous disorder with a dynamic behavior, longitudinal studies are vital to capture brain changes over time, and establish longitudinal relationships (13, 14).

In this study, we identified brain network properties that change over the first 12 months after a TBI using longitudinally acquired data and investigated if these properties relate to changes in cognitive functioning.

Materials and Methods

Data Acquisition and Subjects

Twenty-three TBI patients were scanned twice on a 3T MRI scanner (Phillips Achieva MRI) initially a mean of 6 days after injury (early phase) and again 1 year later (late phase). The study protocol consisted of structural T1-weighted sequences [magnetization-prepared rapid-acquisition gradient echo, repetition time (TR) = 8.1 ms, echo time (TE) = 4.6 ms, matrix size 240 × 216 × 180, isotropic 1 mm resolution] and diffusion-weighted images (TR/TE = 2,524/71 ms; 24 slices; b = 0; 1000 s mm–2; 16 diffusion directions; 2 × 2 × 6 mm3 resolution), and it was consistent across all sessions. Acquisition methods have been described previously, as has the effect of injury on the basic DTI-derived metrics in these same participants (15). Patients were classified as mild or moderate TBI based on Glasgow Coma Scale (GCS). Five patients had moderate TBI [Glasgow Comma Sale (GCS): median = 12, IQR = 10–12] and 18 had a mild TBI (GCS: median = 14, IQR = 14–15) (Table 1). Twenty-eight healthy controls matched for age, gender, level of education, and National Adult Reading Test (NART, a proxy for pre-injury educational status) were scanned once. Detailed description of imaging protocols and subject's clinical data (e.g., age, gender, level of education, NART, injury mechanism, Glasgow coma scale, loss of consciousness, and post-traumatic amnesia) can be found in Croall et al. (15) and Extended Data (Tables S1, S2, S3).


Table 1. Clinical data of controls and patients with mild or moderate TBI.

DTI Data Processing

Fiber orientation reconstruction was carried out using FMRIB Diffusion Toolbox to create fractional anisotropic (FA) images. FA images were then normalized to standard space by a combination of linear and non-linear transformations (Figure 1) (16). All registrations were visually inspected. In standard space, tracts from the recently published structural connectome atlas space were projected to each individuals' FA image using DSI Studio ( (17). This atlas contains 550,000 white-matter tracts verified manually by experienced neuroanatomists (17). Mean FA along the paths between each of the 90 regions of the AAL atlas were computed and saved in a connectivity matrix for each individual. This approach is inspired by a previous TBI study (18).


Figure 1. Processing pipeline. (A) Input data comprise FA and T1-weighted images and neuropsychological tests. (B) Connectome and AAL atlases were mapped to each individual FA image in ICBM-152 space. Mean FA along the tracts of the connectome atlas between the AAL regions was computed to create a connectivity matrix for each individual. (C) Network metrics were calculated from those matrices. (D) Statistical tests and genetic algorithm based random forest method (GARF) were used to classify patients in different TBI post-injury phases and, thus, identify the network metrics assessing changes following TBI. The change over time of those metrics, ΔX, were further investigated to explain changes in cognitive function, ΔY, using regression models. The change in cognitive function, ΔY, represents the change of the principal components obtained after applying PCA to the neuropsychological tests.

Network metrics were estimated from the connectivity matrices to quantify network integration, segregation, centrality, and resilience to perturbations (19). A detailed description of all network metrics is available in Table S4 in the Extended Data. Network integration and segregation measures describe the ability of a network to integrate information from distributed regions and the ability for local specialized information processing, respectively. Examples of those metrics are strength (the sum of all connection strengths to a given region), local clustering coefficient (measure the connectedness of neighbors of a region), local efficiency (the inverse of the average shortest path length of a node to their nodes), local assortativity (measure the contribution of a region to maintain the integration of information in the network after disruption), eigenvector and page rank centrality (measure the influence of a region has on a network), betweenness centrality (measure how often a region is traversed by the shortest paths in the network), and closeness centrality (measure the average shortest path of a region to all other regions in the network). More detailed descriptions are provided in references (19).

Lesions Probability Map

Lesion masks were manually drawn for 19 patients by a neurosurgeon on the baseline scans, as described in Aribisala et al. (20). Lesion masks are defined as a binary mask where the lesioned voxels are labeled as 1 and the remaining voxels are labeled as 0. The masks include any visible contusion, hematoma, or edema. Visible lesions were not present for the remaining four patients. Lesion masks were then co-registered to the MNI-152 space using the warp fields of T1-MRI registration to standard space, which was conducted with linear and non-linear transformations in FSL (16).

To create a probability map, lesion masks were summed, and normalized to 0–100% range with 100% (0%) indicating the presence of a lesion in all (no) patients. Due to the heterogeneity of TBI, we applied a threshold to keep the most prevalent lesions and remove more sporadic ones. Therefore, only the lesions observed in at least three patients were considered.

We additionally calculated the Sørensen–Dice similarity coefficient to measure the overlap/similarity ratio of the Lesions Probability Map as an ROI with each AAL ROIs. This allowed us to inspect the impact of injuries to the observed changes in network metrics at a region level in our cohort.

Components of Cognitive Function

Patients and controls underwent a full battery of standardized neuropsychological tests sensitive to cognitive impairments in mild TBI at each time of scanning. The tests included assessments of attention, memory, executive functions, and semantic knowledge. Due to missing data for some follow-up assessments, we restricted our analysis to the following tests: Speed of Information Processing (SOIP), Paced-Auditory-Serial-Addition Test (PASAT), D-KEFS Color-Word-Interference Test (CWIT), and List Learning (2123). Early and late post-injury neuropsychological scores were compared with controls by inferential tests with false discovery rate (FDR) correction.

Following Irimia et al. (24) and Kuceyeski et al. (25), principal component analysis (PCA) was used to combine multiple cognitive scores in a set of principal components that maximize the variance of the neuropsychological tests and represent main cognitive functions. Since PCA is sensitive to skewness (26), cognitive scores with z-skewness higher than 1.96 were prior transformed using Box-Cox transformation (27).

Network Metrics Selection and Association With Cognitive Function

We conducted a two-step feature selection to identify the most discriminative network metrics between different TBI phases using R ( This allowed us to identify a set of network metrics that changed over time and we interpreted as being causally related to the injury. Since the two-step feature selection was a cross-sectional analysis, a sample size of 46 was used (i.e., 23 samples for each TBI phase/time point or class).

Finally, using a regression model, we investigate how changes in those network metrics may be predictive of cognitive function over time.

Connectivity Metrics of MTBI

Initially, we used inferential tests to rank standardized network metrics that differed significantly between mTBI patients in the early and late post-injury phases. Following Thatcher et al. (28), Barik et al. (29), and Haury et al. (30), correction for multiple comparison was considered not relevant as the goal of this step was to reduce the number of features by separating the most significant metrics from the less significant and not draw inferential conclusion.

In the second step, a cross-sectional and multivariate study was conducted to select combinations of the previously selected metrics that best discriminate between patients in early or late post-injury time phases. To that end, a 5-fold cross-validation with genetic algorithm random forest (GARF) approach was used to classify patients in the different TBI phases. Although classification algorithms are often used to generate predictive models, our main intention was to identify the most discriminative network metrics between TBI phases. We opted to use GARF as it boosts the Random Forests performance while reducing considerably the number of metrics needed for classification (31).

Longitudinal Analysis

After identifying network metrics discriminating between different TBI phases, we investigated their association with cognitive score. This was achieved by a linear regression model, as described in the following equation:

ΔY=β0+β1ΔX1+β2ΔX2++βNΔXN    (1)

where Δ represents the difference between late and early post-injury data for the N standardized network metrics (X) and the new neuropsychological score obtained by PCA (Y). The fitting of the linear regression model was achieved using stats package in R, which uses the least-squares method to minimize the sum of the squares of residuals. Initially, we included all selected network metrics in the linear regression model. However, to overcome multicollinearity and remove unnecessary metrics, we discarded, at each round, the metric with the highest p-value, until the regression model and its coefficients were statistically significant.

Experimental Design and Statistical Analysis

Inferential statistical tests were used to compare neuropsychological scores, PCA components, and network metrics for both groups. Depending on whether parametric assumptions were satisfied, paired t-test or Wilcoxon signed-rank and independent t or Mann–Whitney U were used at a significance level of 5%. In multiple comparisons, FDR correction was applied at a significance level of 5%.

The regression model was evaluated by normalized root mean square error (nRMSE), adjusted R2, and residuals inspection. To verify whether some of the underlying assumptions of regression had been violated, variance inflation factor (VIF) was used to identify multicollinearity among regression network metrics, and the inspection of residuals included the analysis of the scatter plot of residuals vs. predicted scores, lag, and normal probability plots of residuals. F test was used to test the null hypothesis that the fitting either using intercept-only model or our proposed regression model is the same. The rejection of the null hypothesis suggests a better fitting by the proposed model and its predictors.

The final regression model was investigated to examine the presence of bias, overfitting, and the significance of the network metrics selected by GARF. The former was achieved by randomly selecting 20 patients for training and three for testing. We ran this analysis 1,000 times and tested different split sizes (20–3, 18–5, 19–4, etc.). The significance of the network metrics selected by GARF was assessed by random selection of metrics followed by fitting of a linear regression model to predict changes in new neuropsychological score. The network metrics were selected from the set obtained in the first step of feature selection. We ran this analysis 1,000 times and, for each run, computed the nRMSE, and adjusted R2.

Data Accessibility

Data and code will be made available after acceptance of the manuscript.


Network-Based Metrics of MTBI

The GARF model selected a total of 16 network metrics involving 12 brain regions and showed a 5-fold cross-validation accuracy of 83.3%, confirming robust discrimination between early and late post-injury time phases. Figure 2 shows the brain regions and properties that were selected by the model, with a majority located in the fronto-temporal cortex.


Figure 2. Network metrics that showed strongest power in discriminating between early and late post-injury time phases: strength in left posterior and right cingulum, Eigenvector centrality in left middle temporal; Page Ranks centrality in left inferior temporal and left orbital medial frontal; Local Assortativity in right pallidum, left rectus, left orbital superior frontal, left olfactory, and left medial superior frontal; Clustering Coefficient in right orbital superior frontal; Betweenness Centrality in right thalamus and left medial superior frontal; Local Efficiency in right anterior cingulum and left medial superior frontal and Closeness Centrality in left medial superior frontal. Colors illustrate the network metrics identified in each region. L, Left Side; R, Right Side.

Cognitive Function Assessment

Inferential statistics were used to compare cognitive functioning between patients in early or late post-injury phase along with controls (Figure 3). Cognitive functioning was assessed by PCA to reduce Type 1 errors commonly found in multiple statistical comparisons. This approach resulted in one main component, explaining 62.3% of the variance. PASAT and CWIT showed the highest weightings/contributions to the first PCA component, supporting the interpretation that the first PCA component predominantly represents executive functions. For ease of interpretation, we will refer to the first PCA component as cognitive function component (CFC). In our recent study, CFC showed to be correlated with a multivariate measure for intrinsic injury severity (32). An increase over time in CFC suggests an overall improvement of cognitive function. A total of 18 patients showed an increase in cognitive functioning at 12 months, four showed a low decay, while one did not show change in its cognitive functioning over 1 year.


Figure 3. Cognitive function component (CFC) for all subjects. The light gray lines display the change of CFC between early and late post-injury phases display for each patient. Significant differences were found between early post-injury CFC and both controls (p-value = 0.001) and late post-injury CFC (p-value < 0.001). As expected, no significant difference was observed between controls and late post-injury CFC (p-value = 0.268), due to the recovery of most patients.

Patients did not show significant differences in the CWIT inhibition trial when compared with controls. In contrast, they showed significant disability in reading and naming colored patches, as quantified by the baseline trials of CWIT. Since reading and naming also involve cognitive pathways (33, 34), we include those measures in the PCA.

Longitudinal Relationships Between Network-Based Metrics of MTBI and Cognitive Functions

After identifying the network metrics associated with mTBI, we performed a longitudinal analysis to investigate associated changes in cognitive performance, namely, in CFC. A significant multiple linear regression equation (F test: p < 0.0001) was achieved by a total of 11 significant network metrics between nine brain regions. All metrics had an acceptable (35) VIF lower than 10. Figure 4A shows the 11 network metrics with the greatest power in explaining changes in cognitive function and their correlation sign with recovery (Figure 4B). The importance of these 11 features to explain changes in cognitive function is ranked based on the standardized coefficients of the regression model, which are available in Table S5 in Extended Data. From the 11 network metrics identified by the model, 45% measured local assortativity, and the most important network change to explain changes in cognitive function was local assortativity in left medial superior frontal. For most of the network metrics, late post-injury patients in our study showed higher or lower median values than controls. However, an increase on the network metrics was observed early after the injury in anterior cingulate cortex (ACC) and gyrus rectus (GR), followed by a normalization to controls values in the late post-injury phase. The most common of the network metrics assesses local assortativity. Global assortativity quantifies the tendency of regions to be connected to regions with similar strength. A network with positive assortativity is likely to maintain its ability to integrate information after disruption (i.e., it is more resilient to insults), as it may contain interconnected high-strength regions (19). Local assortativity measures the contribution of each node to network global assortativity (36).


Figure 4. Changed network metrics associated with cognitive improvement. Colors in (A) illustrate the network metrics identified in each region and the circles' radius indicates the importance for prediction based on the standardized coefficients of regression model. The list of network changes associated with changes in cognitive function were strength in right cingulum; eigenvector centrality in left middle temporal; page ranks centrality in left orbital medial frontal; local assortativity in right pallidum, left rectus, left orbital superior frontal, left olfactory, and left medial superior frontal; clustering coefficient in right orbital superior frontal; betweenness centrality in left medial superior frontal; and local efficiency in left medial superior frontal cortices. Red circles in (B) represent regions in which an increase in their network metrics is negatively correlated with recovery. Possible justifications to these changes are secondary brain injuries, inability to recover after the trauma, or brain maladaptation. Green circles in (B) illustrate the regions in which an increase in their network metrics over 1 year is positively correlated with improvement, therefore suggesting brain adaptation or recovery. L, Left Side; R, Right Side.

The linear regression equation proposed in this study showed an adjusted R2 of 90% and nRMSE of 0.06. These results indicate that the model explains 90% of the variability found in the CFC over time in patients. Figure 5A shows the bootstrap line close to the ideal line, when the observed and expected responses are equal. This result implies consistency of our model performance and no presence of overfitting. Similar results were achieved with different split sizes (Figure 5B). Figure 5C shows the distribution of nRMSE and adjusted R2 by taking random selections of different sets of 11 network metrics after the first feature selection step. The vertical lines correspond to the results of the significant regression equation. The low p-values suggest that network metrics in Figure 4A are associated with changes following TBI and not random variability.


Figure 5. Regression model of changed network metrics explains changes in cognitive function. (A) Bootstrap line (blue) closely approximates to the ideal line in red (adjusted R2 = 76%), when observed and expected responses are equal. The vertical lines indicate the prediction standard deviation for each patient across 1,000 runs. (B) Bootstrap analysis: for each different split size, we saved the mean predictions after 1,000 runs and repeated this process 10 times. The vertical lines indicate the standard deviation of the adjusted R2 and nRMSE after repeating the bootstrap analysis 10 times. (C) Distribution of adjusted R2 and nRMSE by taking random selection of different sets of discriminative network metrics between early and late post-injury stages 1,000 times.

Figure 6A shows the Dice overlap of the AAL regions assessed by the 11 network metrics with the Probability Map of Lesions (Table S5; Extended Data). Although most of the areas showed overlap, frontal brain areas were particularly affected (Figure 6B).


Figure 6. Lesions overlap with (A) regions associated with cognitive function component (CFC) change and (B) across patients. (A) Dice overlap between Lesions Probability Map in (B) and AAL regions found to be predictive of changes in cognitive function. Node radius indicates the amount of Dice overlap for the predictive AAL regions, which can be found in extended data (Figure 4). (B) Lesions Probability Map for different Z voxel coordinates of MNI space, showing the visible overlapping amount across patients. The lesions are found mainly in the frontal lobe, temporal lobe, and slightly in the anterior cingulate cortex. The colormap ranges between an overlap of 16 and 32% patients.


In this preliminary study of 23 patients, we investigated the existence of longitudinal relationships between changes in structural connectivity and changes in cognitive function following TBI. Using graph theory and a regression model, we identified 11 changes in structural connectivity explaining 90% of the observed variability in cognitive function change over 1 year. The observed brain reorganization suggests possible involvement of frontal regions and the anterior cingulate cortex in cognitive function recovery. The regression model also combines information from other brain regions, emphasizing the value of graph-theoretical network-wide perspectives in cognitive neuroscience.

Limitations of our work include the relatively small size of our dataset, which reflects the challenges of recruiting large longitudinal cohorts with repeated clinical and imaging data. Furthermore, the patient imaging and psychological data were acquired at only two time points, within days of injury and a year later, so we could not establish the starting and stopping point for each different reorganization pattern or its causality. Another limitation in this study is the resolution of DTI images, which could negatively affect fiber tracking. We overcome it by a two-step co-registration followed by quantification of FA values along a connectome atlas. The pathways comprised in this atlas were manually examined and labeled by experienced neuroanatomists (17). A similar approach has also been used in previous severe TBI studies (18). Finally, our dataset comprises patients with both mild and moderate TBI, which could potentially prejudice our understanding about recovery in mild TBI. Using the same dataset, Croall et al. (15) did not observe any significant differences between patients with mild or moderate TBI in any diffusion metrics. Since only one moderate patient had a GCS lower than 10 and the majority had 12, the moderate group is close to the milder end of injuries. As in Croall et al. (15), we assume that these patients have microstructural changes similar to mild injuries.

Areas Associated With Cognitive Function

The regions identified by our regression model are consistent with previous cognitive function localization literature (3742). For instance, Zhou et al. (43) found WM volume loss in the left and right rostral anterior cingulum to be correlated with changes in memory and attention. Furthermore, left medial superior frontal cortex (MSFC), gyrus rectus (GR), and middle temporal cortex (MTC) are known to be associated with specific tasks involved in the neuropsychological tests assessed in this study, such as auditory verbal attention, scores on sequencing, semantic retrieval, semantic memory, and semantic control (4447). Changes in the left olfactory cortex (OC) also contributed to the prediction model: to our knowledge, no cognitive functions are directly supported by it. Although two patients lost sense of smell, the OC has connections to both orbitofrontal cortex and basal ganglia and is implicated in executive function pathways.

Network Topology Changes Over Time After MTBI

The network metric changes observed in MTC, pallidum, and OC showed a positive correlation with ΔCFC (Figure 4B). Those changes suggest strengthening to regions that are highly connected, which may be a brain response to improve or establish a “rich-club organization” for information integration (9), either due to local lesions or disruptions in cognitive subnetworks. Since the pallidum is highly connected to other high-connected regions in the brain (39), the increase in local assortativity suggests that it connects to other similar regions, preferably high-strength regions (Figure 7A). Similarly, a positive correlation between an increase of local assortativity in OC and ΔCFC may indicate a strengthening of its neighbors' connections involved in executive function circuits, such as orbitofrontal and basal ganglia (Figures 7B,C) (40). Some possible explanations for network metric changes, which are negatively correlated with ΔCFC, may be continuous deterioration, incapacity to recover, maladaptive plasticity, or compensatory pathways recruitment, which are no longer required after brain recovery. Increased local assortativity in MSFC was observed to be negatively correlated with ΔCFC (Figure 4B). Since the MSFC is highly susceptible to injury (41), internal disconnections may contribute to the reduction in the average strength differences, which leads to increased local assortativity (Figures 7A,C).


Figure 7. Effect on local disassortativity by changes in regions' strength. (A,B) indicate some possible mechanisms of local disassortativity change. (A) Local disassortativity of the red region decreases when the strength differences between red regions and its neighboring regions is reduced. This reduction may occur by strengthening or weakening of internal connections. When disassortativity decreases, the red region is more assortative. (B) Similarly, local disassortativity decreases with peripheral weakening or peripheral strengthening. When the red region is connected to other regions with lower strength, disassortativity may decrease by strengthening of its low-strength neighbors. However, when connected to regions with higher strength, strengthening of those neighbors will increase disassortativity. (C) Examples on how to calculate local disassortativity for regions colored in red.

We observed an increase in early post-injury phase in anterior cingulate cortex strength and gyrus rectus local assortativity, followed by a decrease in 1 year. This is in agreement with a longitudinal study by five, who also observed increased structural connectivity in a subnetwork within 7 days post-injury in mTBI patients, including pathways from or to anterior cingulate cortex, and gyrus rectus.


Our findings demonstrate a longitudinal relationship between changes in structural brain connectivity and changes in cognitive functions following TBI. The detailed graph theoretical analysis suggests that a combination of different network metrics in distinct brain regions captures most of the longitudinal variance in cognitive performance.

Future longitudinal investigations should assess patients' cognitive outcome and brain networks at more time points to enable personalized predictions of optimal rehabilitation strategies based on network metrics changes at higher temporal resolution.

Data Availability Statement

Data and code to reproduce main figures can be found at:

Ethics Statement

The studies involving human participants were reviewed and approved by Newcastle University Ethics Committee. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

NM: study design and concept, analysis and interpretation of data, statistical analysis, and drafting the manuscript. CC: data acquisition and lesions masks drawing. AB: data acquisition and revised the manuscript. RF: study design and concept, interpreted the data, and revised the manuscript. PT: study design and concept, interpreted the data, imaging processing, and revised the manuscript.


Patient recruitment and acquisition of all MRI scan and experimental data were funded by a Sir Jules Thorn Biomedical Research Award (AMB), SJT/05. PT was funded by Wellcome Trust (105617/Z/14/Z).

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.


The authors would like to thank Tom Kelly and Yujiang Wang for their advice and assistance.

Supplementary Material

The Supplementary Material for this article can be found online at:


1. Forsyth R, Kirkham F. Predicting outcome after childhood brain injury. Cmaj. (2012) 184:1257–64. doi: 10.1503/cmaj.111045

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Werner C, Engelhard K. Pathophysiology of traumatic brain injury. Br J Anaesth. (2007) 99:4–9. doi: 10.1093/bja/aem131

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Lawrence T, Helmy A, Bouamra O, Woodford M, Lecky F, Hutchinson PJ. Traumatic brain injury in england and wales: prospective audit of epidemiology, complications and standardised mortality. BMJ Open. (2016) 6:e012197. doi: 10.1136/bmjopen-2016-012197

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Carroll L, Cassidy JD, Peloso P, Borg J, Von Holst H, Holm L, et al. Prognosis for mild traumatic brain injury: results of the WHO collaborating centre task force on mild traumatic brain injury. J Rehabil Med. (2004) 36:84–105. doi: 10.1080/16501960410023859

CrossRef Full Text | Google Scholar

5. Smith DH, Meaney DF, Shull WH. Diffuse axonal injury in head trauma. J Head Trauma Rehabil. (2003) 18:307–16. doi: 10.1097/00001199-200307000-00003

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Newcombe VF, Correia MM, Ledig C, Abate MG, Outtrim JG, Chatfield D, et al. Dynamic changes in white matter abnormalities correlate with late improvement and deterioration following TBI: a diffusion tensor imaging study. Neurorehabil Neural Repair. (2016) 30:49–62. doi: 10.1177/1545968315584004

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Sharp DJ, Scott G, Leech R. Network dysfunction after traumatic brain injury. Nat Rev Neurol. (2014) 10:156. doi: 10.1038/nrneurol.2014.15

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Dall'Acqua P, Johannes S, Mica L, Simmen HP, Glaab R, Fandino J, et al. Functional and structural network recovery after mild traumatic brain injury: a 1-year longitudinal study. Front Hum Neurosci. (2017) 11:280. doi: 10.3389/fnhum.2017.00280

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Fagerholm ED, Hellyer PJ, Scott G, Leech R, Sharp DJ. Disconnection of network hubs and cognitive impairment after traumatic brain injury. Brain. (2015) 138:1696–709. doi: 10.1093/brain/awv075

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Aerts H, Fias W, Caeyenberghs K, Marinazzo D. Brain networks under attack:robustness properties and the impact of lesions. Brain. (2016) 139:3063–83. doi: 10.1093/brain/aww194

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Bullmore ET, Bassett DS. Brain graphs: graphical models of the human brain connectome. Annu Rev Clin Psychol. (2011) 7:113–40. doi: 10.1146/annurev-clinpsy-040510-143934

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Bullmore E, Sporns O. Complex brain networks: graph theoretical analysis of structural and functional systems. Nat Rev Neurosci. (2009) 10:186–98. doi: 10.1038/nrn2575

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Solem RC. Limitation of a cross-sectional study. Am J Orthod Dentofacial Orthop. (2015) 148:205. doi: 10.1016/j.ajodo.2015.05.006

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Cole JH, Jolly A, de Simoni S, Bourke N, Patel MC, Scott G, et al. Spatial patterns of progressive brain volume loss after moderate-severe traumatic brain injury. Brain. (2018) 141:822–36. doi: 10.1093/brain/awx354

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Croall ID, Cowie CJ, He J, Peel A, Wood J, Aribisala BS, et al. White matter correlates of cognitive dysfunction after mild traumatic brain injury. Neurology. (2014) 83:494–501. doi: 10.1212/WNL.0000000000000666

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Jenkinson M, Beckmann CF, Behrens TE, Woolrich MW, Smith SM. FSL. Neuroimage. (2012) 62:782–90. doi: 10.1016/j.neuroimage.2011.09.015

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Yeh F-C, Panesar S, Fernandes D, Meola A, Yoshino M, Fernandez-Miranda JC, et al. Population-averaged atlas of the macroscale human structural connectome and its network topology. Neuroimage. (2018) 178:57–68. doi: 10.1016/j.neuroimage.2018.05.027

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Squarcina L, Bertoldo A, Ham TE, Heckemann R, Sharp DJ. A robust method for investigating thalamic white matter tracts after traumatic brain injury. Neuroimage. (2012) 63:779–88. doi: 10.1016/j.neuroimage.2012.07.016

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Rubinov M, Sporns O. Complex network measures of brain connectivity: uses and interpretations. Neuroimage. (2010) 52:1059–69. doi: 10.1016/j.neuroimage.2009.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Aribisala BS, Cowie CJ, He J, Wood J, Mendelow DA, Mitchell P, et al. A histogram-based similarity measure for quantitative magnetic resonance imaging: application in acute mild traumatic brain injury. J Comput Assist Tomogr. (2014) 38:915–23. doi: 10.1097/RCT.0000000000000143

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Delis DC, Kaplan E, Kramer JH. Delis-kaplan executive function system. Examiner's Man. (2001). doi: 10.1037/t15082-000

CrossRef Full Text | Google Scholar

22. Gronwall DMA. Paced auditory serial-addition task: a measure of recovery from concussion. Percept Mot Skills. (1977) 44:367–73. doi: 10.2466/pms.1977.44.2.367

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Coughhlan AK, Oddy M, Crawford JR. Birt Memory and Information Processing Battery London. Wakefield: Brain Injury Rehabilitation Trust (2007).

Google Scholar

24. Irimia A, Wang B, Aylward SR, Prastawa MW, Pace DF, Gerig G, et al. Neuroimaging of structural pathology and connectomics in traumatic brain injury: toward personalized outcome prediction. NeuroImage Clin. (2012) 1:1–17. doi: 10.1016/j.nicl.2012.08.002

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Kuceyeski A, Maruta J, Niogi SN, Ghajar J, Raj A. The generation and validation of white matter connectivity importance maps. Neuroimage. (2011) 58:109–21. doi: 10.1016/j.neuroimage.2011.05.087

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Hubert M, Rousseeuw P, Verdonck T. Robust PCA for skewed data and its outlier map. Comput Stat Data Anal. (2009) 53:2264–74. doi: 10.1016/j.csda.2008.05.027

CrossRef Full Text | Google Scholar

27. Rigon J, Burro R, Guariglia C, Maini M, Marin D, Ciurli P, et al. Self-awareness rehabilitation after traumatic brain injury: a pilot study to compare two group therapies. Restor Neurol Neurosci. (2017) 35:115–27. doi: 10.3233/RNN-150538

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Thatcher RW, North D, Biver C. EEG and intelligence: relations between EEG coherence, EEG phase delay, and power. Clin Neurophysiol. (2005) 116:2129–41. doi: 10.1016/j.clinph.2005.04.026

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Barik K, Daimi SN, Jones R, Bhattacharya J, Saha G. A machine learning approach to predict perceptual decisions: an insight into face pareidolia. Brain Informatics. (2019) 6:2. doi: 10.1186/s40708-019-0094-5

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Haury A-C, Gestraud P, Vert J-P. The influence of feature selection methods on accuracy, stability and interpretability of molecular signatures. PLoS ONE. (2011) 6:e28210. doi: 10.1371/journal.pone.0028210

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Bader-El-Den M, Gaber M. GARF: towards self-optimised random forests. In: International Conference on Neural Information Processing. Doha: Springer (2012). p. 506–15. doi: 10.1007/978-3-642-34481-7_62

CrossRef Full Text | Google Scholar

32. Taylor PN, da Silva NM, Blamire A, Wang Y, Forsyth R. Early deviation from normal structural connectivity: a novel intrinsic severity score for mild TBI. Neurology. (2020) 94:e1021–6. doi: 10.1212/WNL.0000000000008902

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Kim V. Executive functions and word reading fluency: a brief intervention with economically disadvantaged secondary students. PCOM Psychology Dissertations (2017). p. 401. Available online at:

Google Scholar

34. Messer D, Henry LA, Nash G. The relation between executive functioning, reaction time, naming speed, and single word reading in children with typical development and language impairments. Br J Educ Psychol. (2016) 86:412–28. doi: 10.1111/bjep.12115

PubMed Abstract | CrossRef Full Text | Google Scholar

35. O'brien RM. A caution regarding rules of thumb for variance inflation factors. Qual Quant. (2007) 41:673–90. doi: 10.1007/s11135-006-9018-6

CrossRef Full Text | Google Scholar

36. Thedchanamoorthy G, Piraveenan M, Kasthuriratna D, Senanayake U. Node assortativity in complex networks: an alternative approach. Procedia Comput Sci. (2014) 29:2449–61. doi: 10.1016/j.procs.2014.05.229

CrossRef Full Text | Google Scholar

37. Barbey AK, Koenigs M, Grafman J. Orbitofrontal contributions to human working memory. Cerebral Cortex. (2010) 21:789–95. doi: 10.1093/cercor/bhq153

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Bubb EJ, Metzler-Baddeley C, Aggleton JP. The cingulum bundle: anatomy, function, and dysfunction. Neurosci Biobehav Rev. (2018). 92:104–27. doi: 10.1016/j.neubiorev.2018.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Leisman G, Braun-Benjamin O, Melillo R. Cognitive-motor interactions of the basal ganglia in development. Front Syst Neurosci. (2014) 8:16. doi: 10.3389/fnsys.2014.00016

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Mai JK, Paxinos G. The Human Nervous System. (2011). Academic Press.

Google Scholar

41. McAllister TW. Neurobiological consequences of traumatic brain injury. Dialogues Clin Neurosci. (2011) 13:287. doi: 10.1176/appi.books.9781585624201.js15

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Peng K, Steele SC, Becerra L, Borsook D. Brodmann area 10: collating, integrating and high level processing of nociception and pain. Prog Neurobiol. (2018) 161:1–22. doi: 10.1016/j.pneurobio.2017.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Zhou Y, Kierans A, Kenul D, Ge Y, Rath J, Reaume J, et al. Mild traumatic brain injury: longitudinal regional brain volume changes. Radiology. (2013) 267:880–90. doi: 10.1148/radiol.13122542

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Elderkin-Thompson V, Hellemann G, Pham D, Kumar A. Prefrontal brain morphology and executive function in healthy and depressed elderly. Int J Geriatr Psychiatry. (2009) 24:459–68. doi: 10.1002/gps.2137

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Nakai T, Kato C, Matsuo K. An fMRI study to investigate auditory attention: a model of the cocktail party phenomenon. Magn Reson Med Sci. (2005) 4:75–82. doi: 10.2463/mrms.4.75

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Shallice T, Stuss DT, Alexander MP, Picton TW, Derkzen D. The multiple dimensions of sustained attention. Cortex. (2008) 44:794–805. doi: 10.1016/j.cortex.2007.04.002

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Xu J, Wang J, Fan L, Li H, Zhang W, Hu Q, et al. Tractography-based parcellation of the human middle temporal gyrus. Sci Rep. (2015) 5:18883. doi: 10.1038/srep18883

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: TBI, longitudinal, graph theoretical analysis, network, structural connectivity

Citation: Moreira da Silva N, Cowie CJA, Blamire AM, Forsyth R and Taylor PN (2020) Investigating Brain Network Changes and Their Association With Cognitive Recovery After Traumatic Brain Injury: A Longitudinal Analysis. Front. Neurol. 11:369. doi: 10.3389/fneur.2020.00369

Received: 22 October 2019; Accepted: 14 April 2020;
Published: 09 June 2020.

Edited by:

Christian Gaser, Friedrich Schiller University Jena, Germany

Reviewed by:

Yongxia Zhou, University of Southern California, United States
Yunliang Cai, Worcester Polytechnic Institute, United States
Catherine D. Chong, Mayo Clinic Arizona, United States

Copyright © 2020 Moreira da Silva, Cowie, Blamire, Forsyth and Taylor. 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: Peter Neal Taylor,

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.