Personalised structural connectomics for moderate-to-severe traumatic brain injury
==================================================================================

* Phoebe Imms
* Adam Clemente
* Evelyn Deutscher
* Ahmed M. Radwan
* Hamed Akhlaghi
* Paul Beech
* Peter H Wilson
* Andrei Irimia
* Govinda Poudel
* Juan F Domínguez D
* Karen Caeyenberghs

## Abstract

Graph theoretical analysis of the structural connectome has been employed successfully to characterise brain network alterations in patients with traumatic brain injury (TBI). However, heterogeneity in neuropathology is a well-known issue in the TBI population, such that group comparisons of patients against controls are confounded by within-group variability. Recently, novel single-subject profiling approaches have been developed to capture inter-patient heterogeneity. We present a personalised connectomics approach that examines structural brain alterations in six chronic patients with moderate-to-severe TBI who underwent anatomical and diffusion magnetic resonance imaging (MRI). We generated individualised profiles of lesion characteristics and network measures (including personalised graph metric ‘GraphMe’ plots, and nodal and edge-based brain network alterations) and compared them against healthy reference cases (N=12) to assess brain damage qualitatively and quantitatively at the individual level. Our findings revealed clinically significant alterations of brain networks with high variability between patients. Our profiling can be used by clinicians to formulate a neuroscience-guided integrative rehabilitation program for TBI patients, and for designing personalised rehabilitation protocols based on their unique lesion load and connectome.

Keywords
*   traumatic brain injury
*   structural connectomics
*   graph theory
*   personalised medicine
*   personalised connectomics
*   lesion filling

Moderate-to-severe traumatic brain injury (TBI) can result in diverse focal lesions and white matter pathology. The locations of these lesions greatly contribute to functional outcomes following TBI, whereby cognitive functions that rely on broadly distributed circuits in the brain are affected due to disruptions to axonal pathways and cortical structures1,2,3. In TBI patients, diffusion weighted MRI (dMRI) studies have shown altered topological properties of structural brain networks, as indexed by graph metrics at the group level4-7. In our recent meta-analysis8, we found that only two of 14 graph metrics (characteristic path length and normalised clustering coefficient) showed significant differences in TBI patients compared with controls, reflecting the heterogeneous nature of TBI patients. This heterogeneity, including complex structural profiles, variation in lesion location, severity, response to treatment, as well as varied secondary injury trajectories, poses a challenge for the prediction of functional and cognitive symptoms of TBI patients. As a result, there is growing impetus for subject-tailored approaches that enable injury characterization and treatment planning9-11.

Recent studies have addressed heterogeneity in clinical cohorts by performing individualised analyses of dMRI-derived fractional anisotropy (FA), T1-derived cortical thickness, and streamline counts11-13 at the level of white matter (WM) tracts or grey matter (GM) regions, respectively. For example, Lv and colleagues13 found no group consensus in anatomic locations of lower FA and reduced cortical thickness in schizophrenia patients, and as such group-level FA and cortical thickness maps were not representative of individuals. To date, however, few studies have analysed brain networks at the level of individual patients, an approach known as *personalised connectomics*10.

Pioneered by Irimia and colleagues10, personalised connectomics enables the use of an individual’s brain network as a ‘fingerprint’ of brain network topology14,15. Personalised connectomics allows the visualization of individual white matter atrophy profiles (as indexed by dMRI-inferred streamline counts) using circular plots and considering patients scores relative to a healthy cohort. These individualized graphs can be used by clinicians to develop personalised rehabilitation programs, by detailing network level abnormalities9 that may indicate specific cognitive deficits following injury. No study to date has examined TBI patients’ network alterations using graph metrics, whereby a literature-driven selection of graph metrics that summarise segregation, integration, and centrality are represented for individual patients16. Since graph metrics were recently shown to have prognostic potential7,17, this type of approach could provide valuable information to clinicians, leading to neuroimaging-guided strategies to improve functional outcomes of TBI patients. However, personalised connectomics in moderate to severe TBI cohorts with diverse brain injuries pose a serious technical challenge, as the available tools for MRI processing to generate connectomes fail in such18.

The present study introduces personalized measurement and analysis of individual connectomic profiles in five chronic moderate-to-severe TBI patients with varying lesion loads, mechanisms of injury, age at injury, and burden of neural/cognitive symptoms. Our implementation extends current methods by addressing the long-standing and prominent challenge of analysing TBI structural profiles when automatic sub/cortical segmentation or parcellation of MRIs fail in the presence of lesions18. Significantly, this problem is addressed here by synergizing connectomic analysis with *virtual brain repair*, where the lesion is replaced by healthy looking tissue in the T1-weighted images (lesion inpainting). The capabilities of our implementation of personalised connectomics in TBI include: (*i*) *lesion masking* undertaken in a semi-automated manner from anatomical T1 MRI scans to identify the affected brain regions in individual patients; (*ii*) the use of the recently developed *Virtual Brain Grafting* (VBG) toolbox to overcome the challenges of segmentation and parcellation of focal lesions using lesion inpainting19; (*iii*) graphical representation of the structural connectome using innovative tools for graph metric profiling (GraphMe plots) to delineate subject-specific changes in brain network integration, segregation, and centrality; and i*v*) regional assessment of network hub regions and edge alterations in individual TBI cases. Together, these innovative solutions overcome major, longstanding methodological impediments in the field of macroscale TBI profiling. Our implementation is the first to allow the comprehensive generation of lesion-aware connectomic profiles that may be used to inform clinical decision making, thus moving closer to the crucial aim of translating research findings into clinical practice.

## Materials and Methods

### Participants

Patients with chronic moderate-to-severe TBI were recruited from St Vincent’s Hospital in Melbourne (SVHM). The definition of moderate-to-severe TBI was based on (i) a Glasgow Coma Scale score at the time of hospital admission between 3-1220; (ii) loss of consciousness of at least 30 minutes; (iii) post-traumatic amnesia of at least 24 hours21: and (iv) positive findings of gross injury on MRIs as per evaluation by a neuroradiologist (PB). Patients who met the following inclusion criteria were contacted to take part in the study: (a) between 18 and 65 years of age; (b) no history of head injury prior to the TBI for which they were included in this study; (c) fluency in English, (d) no history of psychiatric illness prior to the TBI, and (e) no contra-indications for MRI. Ten moderate to severe TBI patients who had sustained closed head injuries due to sports or motor-vehicle accidents more than 6 months prior to the study were recruited. Informed written consent was obtained from each subject in accordance with the Declaration of Helsinki. Due to time constraints during scanning, dMRI were not acquired from four TBI patients, who were subsequently removed from further analysis (see Table 1). One further participant was removed from personalised connectome construction due to excess movement in the scanner during dMRI, which caused a severe motion artifact (see Supplementary Material 2). For the reference group, 12 healthy controls were recruited from the general population using flyers and the snowball method (see Table 1). Ethical approval was granted by the SVHM ethics committee for human research (project #250/17).

View this table:
[Table 1.](http://medrxiv.org/content/early/2022/03/04/2022.03.02.22271654/T1)

Table 1. 
Participant Demographics and Injury Characteristics

### Data acquisition

MRI scans were acquired at the Royal Children’s Hospital using a 3T Siemens PRISMA with a 64-channel head coil. dMRI data were acquired using a single-shot echo planar imaging sequence (twice-reinforced spin echo, multi-band acceleration factor of 2, 70 contiguous sagittal slices) and a high angular resolution diffusion imaging (HARDI) gradient scheme with 66 non-collinear gradient directions (total acquisition time (TA) = 6:17 mins, *b* = 3000 s/mm2, field of view (FOV) = 260 mm2, voxel size = 2.3 mm isotropic, repetition time (TR) = 3500 ms, echo time (TE) = 67 ms, seven volumes with *b* = 0, two reverse phase-encoded volumes with *b* = 0, *b* being the constant of diffusion weighting). T1-weighted MRIs were also acquired using a magnetisation-prepared rapid acquisition gradient-echo (TA = 5:48 mins, 208 contiguous slices, FOV = 256 mm2, voxel size = 0.8 mm isotropic, TR = 2100 ms, TE = 2.22 ms, flip angle = 8°).

### Lesion masking

Manual lesion delineation for computation of lesion load and for improvement of anatomical segmentation was performed by an assessor (ED), who was trained in lesion identification by neuroradiologist (PB). Lesions were drawn in the T1 native space using Fsleyes version 0.27.3 in FSL version 6.0.1 ([https://fsl.fmrib.ox.ac.uk/fsl/fslwiki](https://fsl.fmrib.ox.ac.uk/fsl/fslwiki)). An in-house systematic search method and lesion identification protocol was implemented by JD, KC, ED, and PB. Abnormalities resulting in tissue loss, such as regions of encephalomalacia and damage from surgical drainage tracts were included in binarised lesion masks. Enlarged ventricles and hyperintensities often occurring in proximity to the skull (e.g., from surgical craniotomies) were not included in the lesion masks. Lesion load was computed (in cm3) as the total volume of the binary lesion masks in FSL.

### Personalised connectome construction

Our connectome processing pipeline is showcased in Figure 1, Supplementary Material 1, and in our previous publication23. Our personalised connectomics implementation performs state-of-the-art, single-subject analyses of structural MRI scans. Briefly, raw dMRI data were processed using MRtrix3Tissue (v5.2.8; [https://3tissue.github.io](https://3tissue.github.io)), a fork of MRtrix324. White matter fibre orientation distributions were estimated using single-shell 3-tissue constrained spherical deconvolution (SS3T-CSD)25,26. Whole-brain, anatomically-constrained tractography was performed27 and twenty-two million streamlines were generated per subject28. The SIFT2 algorithm was applied to match the fibre density of the reconstructed streamlines to that of the underlying white matter structures28-30. Thus, edges encode filtered streamlines count.

![Fig. 1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/04/2022.03.02.22271654/F1.medium.gif)

[Fig. 1](http://medrxiv.org/content/early/2022/03/04/2022.03.02.22271654/F1)

Fig. 1 
Overview of the processing pipeline for connectome mapping. (A) In the raw diffusion images, noise33,34, Gibbs ringing artefacts35, as well as distortions induced by motion, eddy current artefacts and EPI/susceptibility distortions were detected and corrected36,37. (B) Concurrently, T1 volumes were registered to diffusion volumes. The advanced normalisation tools package (ANTS)38 was used to remove non-brain structures from the T1 weighted images for white matter extraction39. FSL FLIRT40,41 was used to perform the boundary-based registration between brain-extracted anatomical and diffusion images. Registered images are provided to (i) 5ttgen (brain extracted), to create priors for anatomically constrained tractography (ACT), and (ii) FreeSurfer (non-brain extracted), to parcellate the nodes for the connectome analysis. All subcortical grey-matter structures were segmented42; image intensity normalised43; pial surfaces and the grey-white matter boundaries estimated44; and the entire brain “inflated” to smooth the gyri and sulci45. (C) Lesion maps of subjects who failed the quality assessment after FreeSurfer parcellation were provided along with the T1 image to VBG. (D) Average response functions for white matter, grey matter, and cerebrospinal fluid were estimated from the dMRI data using an automated unsupervised approach46,47. Pre-processed data were upsampled to a voxel size of 1.3 mm3 to assume higher spatial resolution for image registration before binary masks were created. Fibre orientation distributions (FODs) were estimated from the group average response functions on upsampled images, and corrected for intensity inhomogeneity and global intensity level differences48. (E) Anatomically constrained tractography (ACT) was performed using the FODs from (D) and the 5ttgen images from (B(i)). The FOD cut-off threshold, step size, and angle were determined to attain a reasonable trade-off between false negatives and false positives (seed points=dynamic; maximum length=250mm; minimum length=5 mm; step size=1.25; angle=45°; FOD amplitude = 0.08). Spherically informed filtering of tractograms (SIFT2) is applied to make the weight of the streamlines proportional to the underlying fibre orientation distribution. (F) The connectome is created using the FreeSurfer parcellation and the sifted tractogram.

T1 anatomical MRIs were parcellated into 84 regions of the Desikan-Killiney atlas31 using FreeSurfer’s *recon-all* function (v6.0; [http://surfer.nmr.mgh.harvard.edu/](http://surfer.nmr.mgh.harvard.edu/))32. Two patients (TBI3 and TBI4) had significant segmentation failures due to gross pathology, and were therefore processed utilising VBG v0.3719. Rather than lesion masking and manual editing which are subjective and time-consuming, VBG automatically fills uni- and bilateral brain lesions using synthetic healthy donor tissue to permit or to improve segmentation. To illustrate the performance of VBG in TBI, we included a report on VBG outcome for patient TBI2, who was excluded from personalised connectomics due to movement during HARDI acquisition but otherwise had a quality control compliant T1-weighted volume (Supplementary Material 3). Given VBG artificially reconstructs lesioned nodes, part of our quality control also included ensuring streamlines were not aberrantly assigned to these nodes. Connectivity matrices were generated using edge weights from SIFT2 and nodes defined as brain regions from FreeSurfer and VBG. Global network properties were quantified in terms of *strength, global efficiency, characteristic path length, navigation efficiency, average local efficiency, clustering coefficient, normalised clustering coefficient*, and *average betweenness centrality* (Table 2) using the Brain Connectivity Toolbox16. These graph metrics were chosen from all available metrics as the most clinically informative/intuitive according to previous findings from our meta-analysis8. Grading of diffuse axonal injury (DAI) was performed by expert raters PB and ED (Table 1)22.

View this table:
[Table 2.](http://medrxiv.org/content/early/2022/03/04/2022.03.02.22271654/T2)

Table 2. 
Graph Metric Descriptions and Interpretations

### Brain network profiles

Graph metrics spiderplots (‘GraphMe’ plots) show results for each TBI patient in a concise and intuitive manner relative to mean scores from the healthy controls with 95% confidence intervals. Selected graph metrics (characteristic path length, normalised clustering coefficient, and betweenness centrality – Table 2) were inverted (1/*x*) to facilitate interpretation (so that *higher* scores on any graph metric denote *better* brain network structure). Graph metrics of individual patients were categorised as follows: *normal* (if the scores/metrics fell within the 95% confidence interval); *supra-normal* (higher than the 95% confidence interval); or *infra-normal* (lower than the 95% confidence interval)13.

### Regional brain network analyses

A key component of personalised connectomics is to localise network alterations in the brain relative to a healthy cohort. Nodal hubs and weakest edges were also examined for each individual patient based on comparison to the healthy controls. *Betweenness centrality* was used to identify brain regions essential for communication within the brain network16,59. The top 10% (*n*=8) highest-scoring nodes were identified as hubs; for the healthy control group these are shown in Figure 2, and for the TBI patients these are shown in Figures 3(a) to (d).

![Fig. 2](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/04/2022.03.02.22271654/F2.medium.gif)

[Fig. 2](http://medrxiv.org/content/early/2022/03/04/2022.03.02.22271654/F2)

Fig. 2 
Healthy control hubs (top 10% of nodes with highest betweenness centrality), in teal. Larger nodes represent higher values. Hubs (bilaterally) were the superior frontal gyrus (BCleft=1493; BCright=1533), superior parietal gyrus (BCleft=610; BCright=665), precentral gyrus (BCleft=588; BCright=616), and thalamus (BCleft=336; BCright=346). The strongest edges (0.5th percentile) are coloured by strength (yellow=weaker; red=stronger). Visualisation in NeuroMArVL ([https://immersive.erc.monash.edu/neuromarvl/](https://immersive.erc.monash.edu/neuromarvl/)).

![Fig. 3](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/04/2022.03.02.22271654/F3.medium.gif)

[Fig. 3](http://medrxiv.org/content/early/2022/03/04/2022.03.02.22271654/F3)

Fig. 3 
Personalised connectome profile for TBI1 including (A) lesion profile; (B) quality assessment; (C) radar plot showing the patient’s personalised connectome profile (red indicates patient’s scores, dark blue indicates healthy control average and the 95% Cis are represented by the blue shade); and (D) regional analysis (blue: edges lower than the healthy control average; red: edges stronger than the healthy control average; thicker edges: more standard deviations away from the healthy mean).

### Z-score Matrix for Edge Analysis

An edge analysis scrutinised the white matter connections that drive overall differences in the network properties in greater detail60. A *z*-score matrix *Z**i,j* was derived, which describes the distance from the healthy control mean, divided by the healthy control standard deviation, between each subject’s connectivity matrix *T**i,j* and the controls *H**i,j* according to equations from a previous edgewise analysis60: ![Formula][1]</img>  Positive scores represent stronger edges in the TBI patient compared to controls, while negative scores represent weaker edges. To visualise the *z*-score matrix, only edges with a score >4 (i.e., edges with weights that are 4 standard deviations away from the healthy control mean, corresponding to a highly conservative *p*<.0001) are displayed on a glass brain, while all other edges are discarded.

## Results

### TBI1

TBI1 (Figure 3) had a relatively small lesion load (0.75 cm3) spanning the posterior segment of the right superior frontal gyrus and right precentral gyrus, and a DAI grade of 0. Registration between structural and diffusion images was unaffected by this lesion. There were no failures in the FreeSurfer pipeline and there was no need for VBG. FODs were generated at the site of the lesion (see red arrow) but did not meet streamline criteria for ACT. The GraphMe plot indicated that TBI1 has slightly weaker integration than healthy controls, including infra-normal lower navigation efficiency, as well as lower strength and clustering coefficient. Regional brain network analyses revealed, interestingly, that the right superior frontal gyrus was not a hub node in TBI1. In addition, four other alterations in the hub arrangement for TBI1 were observed, whereby the nucleus accumbens (BCleft=1570; BCright=1546), palladium (BCleft=1382; BCright=978) and right putamen were bilateral hubs (BCright=1578), and the bilateral precentral gyri, thalami, and right superior frontal gyrus did not meet the hub threshold. Weaker edges (*z* < 4) (*n*=43 out of a total of 7056 edges) were observed projecting across frontal, parietal, temporal, and subcortical areas, in particular the edges between (A) the left posterior cingulate cortex and the right frontal pole (*z*=-8.32); (B) the left superior temporal gyrus and the left frontal pole (*z*=-6.66); (C) the left lateral orbitofrontal gyrus and the left temporal pole (*z*=-7.98); and (D) the left medial frontal and left temporal pole (*z*=-6.90). Some stronger (*z* > 4) edges (*n*=4) were also observed, including the connection between the right superior temporal gyrus and the right temporal pole (*z*=5.92).

### TBI3

TBI3 (Figure 4) had a relatively large lesion load (15.46cm3) involving primarily frontal regions (predominantly on the left), white matter hyperintensities in the medial right parietal lobe and the corpus callosum, and a DAI grade of 2. Prior to VBG, 10 nodes failed the quality assessment: VBG repaired 9 nodes for parcellation. Registration showed that streamlines were not assigned to lesioned nodes. The GraphMe plot demonstrated an infra-normal graph metric profile in all domains except centrality, which was normal. Two hub alterations were observed, whereby the bilateral putamen (BCleft=871; BCright=932) were hubs, and the bilateral precentral gyri were not. Weaker edges (*n*=62) projected across the whole brain, especially the frontal regions, including between the left frontal pole and the left middle temporal gyrus (*z*=8.45), right superior frontal gyrus (*z*=8.12), right lateral orbitofrontal gyrus (*z*=8.17), and right putamen (*z*=8.41); between left medial orbitofrontal gyrus and the left amygdala (*z*=8.67); right lateral orbitofrontal gyrus and right nucleus accumbens (*z*=8.03) and right caudate nucleus (*z*=8.33); and right pars orbitalis and right lingual gyrus (*z*=8.10). No stronger edges were observed.

![Fig. 4](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/04/2022.03.02.22271654/F4.medium.gif)

[Fig. 4](http://medrxiv.org/content/early/2022/03/04/2022.03.02.22271654/F4)

Fig. 4 
Personalised connectome profile for TBI3 including (A) lesion profile, (B) quality assessment, (C) radar plot showing the patient’s personalised connectome profile, and (D) regional analysis (blue: edges weaker than healthy control average; thicker edges: more standard deviations away from the healthy mean).

### TBI4

TBI4 (Figure 5) had a relatively large lesion load (17.59 cm3) involving bilateral frontal lesions and right temporal lesions, white matter hyperintensities in the medial right parietal lobe and the corpus callosum. However, the DAI grade was low (0/1). Prior to VBG, 9 nodes failed the quality assessment. All lesions overlapping with these nodes were repaired by VBG. Alignment between VBG-repaired nodes and streamlines indicated that any aberrant streamlines generated in areas with oedema/haemorrhage were not assigned to lesioned nodes. This patient exhibited supranormal graph metrics in all domains except the normalised clustering coefficient of the network (which was normal), and centrality (which was infra-normal). Four alterations in the hub arrangement were observed, whereby the bilateral putamen (BCleft=2246; BCright=1550), left palladium (BCleft=1210) and left inferior parietal (BCright=902) were hubs, and the bilateral precentral gyri and thalamic regions were not hubs. Weaker edges (*n*=26) projected across the left hemisphere, including between the entorhinal and lingual gyrus (*z*=10.87), pericalcarine (*z*=9.55), superior parietal (*z*=9.46), and lateral occipital regions (*z*=9.33); the temporal pole and the insula (*z*=9.37); and the nucleus accumbens and the posterior cingulate cortex (*z*=6.71), insula (*z*=6.30), and rostral anterior cingulate cortex (*z*=5.98). Some stronger edges (*n*=4) were also observed in the right hemisphere, including between the pars triangularis and postcentral gyrus (*z*=5.92), the putamen and lateral orbitofrontal (*z*=4.86), the pallidum and the thalamus (*z*=4.24); and the left pallidum and the amygdala (*z*=4.55).

![Fig. 5](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/04/2022.03.02.22271654/F5.medium.gif)

[Fig. 5](http://medrxiv.org/content/early/2022/03/04/2022.03.02.22271654/F5)

Fig. 5 
Personalised connectome profile for TBI5 including (A) lesion profile, (B) quality assessment, (C) radar plot showing the patient’s personalised connectome profile, and (D) regional analysis (blue: edges weaker than healthy control average; red: edges stronger than healthy control average; thicker edges: more standard deviations away from the healthy mean).

### TBI5

TBI5 (Figure 6) had no MRI-discernible lesion load and a DAI grade of 0. Many weaker edges were observed relative to healthy controls that connected the parietal, temporal, and subcortical lobes. There were no failures in the FreeSurfer pipeline, and no manual edits were necessary. FODs were generated correctly and registration between segmentation and tractography was free of error. The GraphMe plot revealed infra-normal strength, navigation, and normalised clustering. Two alterations in hub arrangement were observed, whereby the bilateral putamina were hubs (BCleft=1182; BCright=1110), whereas the bilateral thalami were not. Weaker edges (*n*=33) projected inter-hemispherically across parietal, temporal, and subcortical areas. The weakest edges (compared to healthy controls) were in the left hemisphere, between the amygdala and the temporal pole (*z*=-7.19); the amygdala and the inferior temporal gyrus (*z*=-7.80); the inferior temporal gyrus and the hippocampus (*z*=-6.11); and the inferior temporal gyrus and the thalamus (*z*=-6.40). Edges stronger than in healthy controls (*n*=3) were also observed, including the set of connections between the left postcentral gyrus and the left lateral occipital gyrus (*z*=5.87).

![Fig. 6](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/04/2022.03.02.22271654/F6.medium.gif)

[Fig. 6](http://medrxiv.org/content/early/2022/03/04/2022.03.02.22271654/F6)

Fig. 6 
Personalised connectome profile for TBI5 including (A) lesion profile, (B) quality assessment, (C) radar plot showing the patient’s personalised connectome profile, and (D) regional analysis (blue=edges lower than healthy control average; red = edges stronger than healthy control average; and thicker edges=larger number of standard deviations away from healthy mean).

### TBI6

TBI6 (Figure 7) had a small lesion in the splenium of the corpus callosum (0.5 cm3), and a DAI grade of 2. There were no failures in the FreeSurfer pipeline, and no manual edits were necessary. FODs were generated at the site of the lesion but did not meet streamline criteria for ACT. The GraphMe plot showed that the global connectivity properties of TBI6 were normal. Three hub alterations were observed, whereby the right caudate nucleus (BCright=722), right hippocampus (BCright=606) and right inferior parietal gyrus (BCright=680) were hubs, and the bilateral precentral and right superior parietal regions were not hubs. No edges met the stringent threshold of being at least 4 standard deviations away from the healthy control mean.

![Fig. 7](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/04/2022.03.02.22271654/F7.medium.gif)

[Fig. 7](http://medrxiv.org/content/early/2022/03/04/2022.03.02.22271654/F7)

Fig. 7 
Personalised connectome profile for TBI6 including (A) lesion profile, (B) quality assessment, (C) radar plot showing the patient’s personalised connectome profile, and (D) regional analysis (blue=edges lower than healthy control average; red = edges stronger than healthy control average; and thicker edges=larger number of standard deviations away from healthy mean).

## Discussion

For the first time, we showcase an implementation of personalised connectomics in chronic moderate-to-severe TBI patients. In the following sections we discuss the defining characteristics of our single-subject profiles and explore ways in which our approach can contribute to improving existing methods of personalised structural connectome analyses in TBI patients.

### Single-subject network profiling observations

Our observations highlight a major caveat to approaches that attempt to identify a single graph metric that can be used as an adequate and parsimonious descriptor of structural network alterations in TBI patients8. In accordance with Lv and colleagues13, we observed that each TBI patient showed a unique *pattern* of graph metric alterations, regardless of lesion load. For example, while both TBI1 and TBI5 patients had small lesion loads, patient TBI1 had lower brain network integration and segregation measures, compared with non-significant deviation from the normal range for patient TBI5. By comparison, TBI3 and TBI4 both had much larger lesions, but patient TBI3’s brain network profile showed infra-normal integration and segregation measures, while patient TBI4’s brain network was supra-normal. Our results highlight the benefit of using a multivariate *profile* of graph metrics that reveal individual differences in brain network topology otherwise obscured by group-level analyses. Importantly, with the incorporation of individual edge and hub comparisons, the location of the lesion can be compared to edge deterioration in single patients.

### Improving methods for personalised connectomics

The present study introduces the structural connectivity aspect of a novel framework called STREAMLINES (Single-subjecT pRofiling of cognitivE impAirMents and structuraL braIN NEtwork metricS), which operates on state-of-the-art single-subject analyses of structural MRI scans that may provide clinicians with a novel user-friendly framework for leveraging graph metrics to benefit the individual patient. Advancing individual brain network profiling has the potential to inform neuroimaging-guided personalised rehabilitation programs by providing network-based summary statistics with prognostic capabilities61-63. More precisely, our approach can help to assess network alternations in TBI patients in the following three ways. First, regional connectome maps can be used as profiles of patients’ brain network topographies, thereby providing clinicians with time-efficient visual summaries of network disruption, asymmetry, hub alterations, and overall reductions in strength. Second, by comparing an individual patient to a healthy control reference group, we can observe portions of brain networks that are topologically altered but correspond to brain regions beyond the site of initial injuries. Finally, the GraphMe plots can be used longitudinally to map how the brain undergoes progressive secondary damage, recovery, and/or functional reorganisation over time64,65.

The current ‘best practice’ methods (lesion masking and manual editing66) are time-consuming and have low inter-rater reliability67. By contrast, use of the semi-automated lesion inpainting program VBG reduces the burden imposed by having to manually delineate lesions and avoids the exclusion of cases with large focal lesions that fail segmentation (e.g., from FreeSurfer)19. Furthermore, we observed that the SS3T-CSD model46,68 was suitable for constructing connectomes in the presence of lesions in all our TBI patients. SS3T-CSD removes the contributions from GM and cerebrospinal fluid (CSF) components to increase the specificity of FODs to the WM, while avoiding over-estimation into GM and CSF signal from the lesioned area26. Combined with anatomically constrained tractography27, streamlines are not generated in lesioned areas (e.g., see TBI1, Figure 3a panel C), and therefore anatomically disconnected regions do not have to be removed from connectivity matrices. This allowed us to calculate graph metrics from connectivity matrices of the same dimensions as those extracted from the healthy controls.

## Limitations

The implementation of personalised connectomics requires extensive validation and assessment of test-retest reliability. However, our study provides initial validation of this approach using six TBI patients and a small healthy control reference group (*N*=12)12. Normative analysis with a large (*N*>100), stratified sample of healthy individuals will allow stronger inferences to be made using techniques such as quartile regression13, 11. Personalised connectomics should also include a patient group as an additional reference cohort (*N*>100), to help clinicians understand how a patient is evolving with reference not only to healthy controls but also to patients with the same condition. Furthermore, our study only utilised T1 images for lesion identification; in the future, other structural imaging modalities such as fluid attenuated inversion recovery (FLAIR) and susceptibility weighted imaging (SWI) should also be used in accordance with best practice guidelines for lesion identification69. Despite multiple expert raters and use of an established procedure22, DAI grading remains subjective and requires independent confirmation of reliability. Hubs were defined according to previous work citing the use of betweenness centrality – however, other work has shown that finding ‘consensus’ between multiple centrality metrics may be more stable for hub definition70.Finally, cognitive outcomes associated with graph measures are still being evaluated – this progress will be essential for providing clinically informative personalised connectomes23 and is included as part of the STREAMLINES platform.

## Conclusions

Our results emphasize the translational potential for single-subject network analyses in the study of brain injury. Profiling individual patients based on their unique presentation provides insights into brain network disruption that are otherwise obscured by group-level approaches. The GraphMe profiling can provide clinicians with a novel user-friendly framework for leveraging graph metrics to benefit the individual patient by characterising network-level brain alterations with potential prognostic relevance. This study therefore enables us to progress towards a personalised medicine approach which, alongside group-based comparisons of patients against controls, is essential for translating connectomics to evidence-based clinical practice.

## Supporting information

Supplementary Materials [[supplements/271654_file03.pdf]](pending:yes)

## Data Availability

All data produced in the present study are available upon reasonable request to the authors

## Authors’ contributions

*Imms, P:* Conceptualisation; Methodology; Software; Formal Analysis; Investigation; Data Curation; Writing – Original Draft; Writing – Review & Editing; Visualisation. *Clemente, A:* Investigation; Resources; Data Curation; Writing – Review & Editing. *Deutscher, E:* Formal Analysis; Writing – Review & Editing. *Ahmed Radwan:* Software; Resources. *Hamed Akhlaghi:* Resources; Data curation. *Beech, P:* Formal Analysis; Writing – Review & Editing. *Wilson, PH:* Conceptualisation; Resources; Writing – Review & Editing; Supervision; Funding acquisition. *Irimia, A:* Writing – Review & Editing. *Poudel, G:* Conceptualisation; Methodology; Software; Validation; Formal Analysis; Resources; Data Curation; Writing – Review & Editing; Visualisation; Supervision. *Domínguez D, JF:* Conceptualisation; Methodology; Software; Validation; Formal Analysis; Resources; Data Curation; Writing – Review & Editing; Supervision. *Caeyenberghs, K:* Conceptualisation; Methodology; Validation; Formal Analysis; Investigation; Resources; Writing – Original Draft; Writing – Review & Editing; Supervision; Project Administration; Funding Acquisition.

## Conflicts of Interest

*Imms, P (corresponding):* None. *Clemente, A:* None. *Ahmed Radwan:* None. *Beech, P:* None. *Wilson, PH:* None. *Irimia, A:* None. *Poudel, G:* None. *Domínguez D, JF (senior):* None. *Caeyenberghs, K (senior):* None.

## Acknowledgements

We would like to thank Dr. Alex Burmester for his help with the Hierarchical Drift Diffusion Modelling; and Michael Kean and the radiographers at the Royal Children’s Hospital.

## Footnotes

*   * Joint senior authors

*   *Email*: Adam.Clemente{at}acu.edu.au

*   *Email*: edeutscher{at}deakin.edu.au;

*   *Email*: ahmed.radwan{at}kuleuven.be;

*   *Email*: Hamed.AKHLAGHI{at}svha.org.au;

*   *Email*: p.beech{at}alfred.org.au;

*   *Email*: PeterH.Wilson{at}acu.edu.au;

*   *Email*: irimia{at}usc.edu;

*   *Email*: Govinda.Poudel{at}acu.edu.au;

*   *Email*: juan.dominguezduque{at}deakin.edu;

*   *Email*: k.caeyenberghs{at}deakin.edu.au;

*   Received March 2, 2022.
*   Revision received March 2, 2022.
*   Accepted March 4, 2022.


*   © 2022, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  1.Bressler, S.L. and Menon, V. (2010). Large-scale brain networks in cognition: Emerging methods and principles. Trends in Cognitive Sciences 14, 277–290.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tics.2010.04.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20493761&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000279089300011&link_type=ISI) 

2.  2.Catani, M. and Ffytche, D.H. (2005). The rises and falls of disconnection syndromes. Brain 128, 2224–2239.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awh622&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16141282&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000232278600003&link_type=ISI) 

3.  3.Hampshire, A., Hellyer, P.J., Parkin, B., Hiebert, N., MacDonald, P., Owen, A.M., Leech, R. and Rowe, J. (2016). Network mechanisms of intentional learning. NeuroImage 127, 123–134.
    
    

4.  4.Caeyenberghs, K., Leemans, A., Leunissen, I., Gooijers, J., Michiels, K., Sunaert, S. and Swinnen, S.P. (2014). Altered structural networks and executive deficits in traumatic brain injury patients. Brain structure & function 219, 193–209.
    
    

5.  5.Kim, J., Parker, D., Whyte, J., Hart, T., Pluta, J., Ingalhalikar, M., Coslett, H.B. and Verma, R. (2014). Disrupted structural connectome is associated with both psychometric and real-world neuropsychological impairment in diffuse traumatic brain injury. Journal of the International Neuropsychological Society 20, 887–896.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S1355617714000812&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25287217&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 

6.  6.Raizman, R., Tavor, I., Biegon, A., Harnof, S., Hoffmann, C., Tsarfaty, G., Fruchter, E., Tatsa-Laur, L., Weiser, M. and Livny, A. (2020). Traumatic brain injury severity in a network perspective: A diffusion MRI based connectome study. Scientific Reports 10, 9121.
    
    

7.  7.Van Der Horn, H.J., Kok, J.G., De Koning, M.E., Scheenen, M.E., Leemans, A., Spikman, J.M. and Van Der Naalt, J. (2017). Altered Wiring of the Human Structural Connectome in Adults with Mild Traumatic Brain Injury. Journal of Neurotrauma 34, 1035–1044.
    
    

8.  8.Imms, P., Clemente, A., Cook, M., D’Souza, W., Wilson, P.H., Jones, D.K. and Caeyenberghs, K. (2019). The structural connectome in traumatic brain injury: A meta-analysis of graph metrics. Neuroscience & Biobehavioral Reviews 99, 128–137.
    
    

9.  9.Irimia, A., Chambers, M.C., Torgerson, C.M., Filippou, M., Hovda, D.A., Alger, J.R., Gerig, G., Toga, A.W., Vespa, P.M., Kikinis, R. and Van Horn, J.D. (2012). Patient-tailored connectomics visualization for the assessment of white matter atrophy in traumatic brain injury. Frontiers in Neurology 3, 10.
    
    

10. 10.Irimia, A., Wang, B., Aylward, S.R., Prastawa, M.W., Pace, D.F., Gerig, G., Hovda, D.A., Kikinis, R., Vespa, P.M. and Van Horn, J.D. (2012). Neuroimaging of structural pathology and connectomics in traumatic brain injury: Toward personalized outcome prediction. NeuroImage: Clinical 1, 1–17.
    
    

11. 11.Jolly, A.E., Balaet, M., Azor, A., Friedland, D., Sandrone, S., Graham, N.S.N., Zimmerman, K. and Sharp, D.J. (2021). Detecting axonal injury in individual patients after traumatic brain injury. Brain 144, 92–113.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awaa372&link_type=DOI) 

12. 12.Attyé, A., Renard, F., Baciu, M., Roger, E., Lamalle, L., Dehail, P., Cassoudesalle, H. and Calamante, F. (2020). TractLearn: A geodesic learning framework for quantitative analysis of brain bundles. medRxiv.
    
    

13. 13.Lv, J., Di Biase, M., Cash, R.F.H., Cocchi, L., Cropley, V.L., Klauser, P., Tian, Y., Bayer, J., Schmaal, L., Cetin-Karayumak, S., Rathi, Y., Pasternak, O., Bousman, C., Pantelis, C., Calamante, F. and Zalesky, A. (2020). Individual deviations from normative models of brain structure in a large cross-sectional schizophrenia cohort. Molecular Psychiatry.
    
    

14. 14.Sanz Leon, P., Knock, S., Woodman, M., Domide, L., Mersmann, J., McIntosh, A. and Jirsa, V. (2013). The Virtual Brain: A simulator of primate brain network dynamics. Frontiers in Neuroinformatics 7.
    
    

15. 15.Schirner, M., Rothmeier, S., Jirsa, V.K., McIntosh, A.R. and Ritter, P. (2015). An automated pipeline for constructing personalized virtual brains from multimodal neuroimaging data. NeuroImage 117, 343–357.
    
    

16. 16.Rubinov, M. and Sporns, O. (2010). Complex network measures of brain connectivity: Uses and interpretations. NeuroImage 52, 1059–1069.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2009.10.003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19819337&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000280181800027&link_type=ISI) 

17. 17.Roine T, Mohammadian M, Hirvonen J, Kurki T, Posti JP, Takala RSK, Newcombe VF, Tallus J, Katila AJ, Maanpää HR, Frantzen J, Menon D and O., T. (2022). Structural Brain Connectivity Correlates with Outcome in Mild Traumatic Brain Injury. Journal Of Neurotrauma.
    
    

18. 18.King, D.J., Novak, J., Shephard, A.J., Beare, R., Anderson, V.A. and Wood, A.G. (2020). Lesion Induced Error on Automated Measures of Brain Volume: Data From a Pediatric Traumatic Brain Injury Cohort. Frontiers in Neuroscience 14.
    
    

19. 19.Radwan, A.M., Emsell, L., Blommaert, J., Zhylka, A., Kovacs, S., Theys, T., Sollmann, N., Dupont, P. and Sunaert, S. (2021). Virtual brain grafting: Enabling whole brain parcellation in the presence of large lesions. NeuroImage 229, 117731.
    
    

20. 20.Teasdale, G. and Jennett, B. (1974). Assessment of coma and impaired conciousness: A practical scale. The Lancet 304, 81–84.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s0140-6736(74)91639-0&link_type=DOI) 

21. 21.Rabinowitz, A.R. and Levin, H.S. (2014). Cognitive sequelae of traumatic brain injury. Psychiatric Clinics of North America 37, 1–11.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.psc.2013.11.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24529420&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 

22. 22.Adams, J.H., Graham, D.I., Murray, L.S. and Scott, G. (1982). Diffuse axonal injury due to nonmissile head injury in humans: An analysis of 45 cases. Annuals of Neurology 12, 557–563.
    
    

23. 23.Imms, P., Dominguez, D.J.F., Burmester, A., Seguin, C., Clemente, A., Dhollander, T., Wilson, P.H., Poudel, G. and Caeyenberghs, K. (2021). Navigating the link between processing speed and network communication in the human brain. Brain structure & function 226, 1281–1302.
    
    

24. 24.Tournier, J.D., Smith, R., Raffelt, D., Tabbara, R., Dhollander, T., Pietsch, M., Christiaens, D., Jeurissen, B., Yeh, C.-H. and Connelly, A. (2019). MRtrix3: A fast, flexible and open software framework for medical image processing and visualisation. NeuroImage 202, 116137.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.NEUROIMAGE.2019.116137&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 

25. 25.Dhollander, T. and Connelly, A. (2016). A novel iterative approach to reap the benefits of multi-tissue CSD from just single-shell (+b=0) diffusion MRI data. In: International Society for Magnetic Resonance in Medicine: Singapore, pps. 3010.
    
    

26. 26.Khan, W., Egorova, N., Khlif, M.S., Mito, R., Dhollander, T. and Brodtmann, A. (2020). Three-tissue compositional analysis reveals in-vivo microstructural heterogeneity of white matter hyperintensities following stroke. NeuroImage 218, 116869.
    
    

27. 27.Smith, R.E., Tournier, J.-D., Calamante, F. and Connelly, A. (2012). Anatomically-constrained tractography: Improved diffusion MRI streamlines tractography through effective use of anatomical information. NeuroImage 62, 1924–1938.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2012.06.005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22705374&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 

28. 28.Yeh, C., Smith, R., Liang, X., Calamante, F. and A, C. (2018). Investigating the streamline count required for reproducible structural connectome construction across a range of brain parcellation resolutions. In: International Society for Magnetic Resonance in Medicine: Paris, pps. 1558.
    
    

29. 29.Smith, R.E., Tournier, J., Calamante, F. and Connelly, A. (2015). SIFT2: Enabling dense quantitative assessment of brain white matter connectivity using streamlines tractography. NeuroImage 119, 338–351.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 

30. 30.Smith, R.E., Tournier, J.-D., Calamante, F. and Connelly, A. (2015). The effects of SIFT on the reproducibility and biological accuracy of the structural connectome. NeuroImage 104, 253–265.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2014.10.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25312774&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 

31. 31.Desikan, R.S., Ségonne, F., Fischl, B., Quinn, B.T., Dickerson, B.C., Blacker, D., Buckner, R.L., Dale, A.M., Maguire, R.P., Hyman, B.T., Albert, M.S. and Killiany, R.J. (2006). An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. NeuroImage 31, 968–980.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2006.01.021&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16530430&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000238358500003&link_type=ISI) 

32. 32.Fischl, B. and Dale, A. (2000). Measuring the thickness of the human cerebral cortex from magnetic resonance images. Proceedings of the National Academy of Sciences 97, 11050–11055.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiOTcvMjAvMTEwNTAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMy8wNC8yMDIyLjAzLjAyLjIyMjcxNjU0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

33. 33.Cordero-Grande, L., Christiaens, D., Hutter, J., Price, A.N. and Hajnal, J.V. (2019). Complex diffusion-weighted image estimation via matrix recovery under general noise models. NeuroImage 200, 391–404.
    
    

34. 34.Veraart, J., Novikov, D.S., Christiaens, D., Ades-aron, B., Sijbers, J. and Fieremans, E. (2016). Denoising of diffusion MRI using random matrix theory. NeuroImage 142, 394–406.
    
    

35. 35.Kellner, E., Dhital, B., Kiselev, V.G. and Reisert, M. (2016). Gibbs-ringing artifact removal based on local subvoxel-shifts. Magnetic Resonance in Medicine 76, 1574–1581.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mrm.26054&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26745823&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 

36. 36.Andersson, J.L.R., Skare, S. and Ashburner, J. (2003). How to correct susceptibility distortions in spin-echo echo-planar images: Application to diffusion tensor imaging. NeuroImage 20, 870–888.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1053-8119(03)00336-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14568458&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000186260700020&link_type=ISI) 

37. 37.Andersson, J.L.R. and Sotiropoulos, S.N. (2016). An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging. NeuroImage 125, 1063–1078.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2015.10.019&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26481672&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 

38. 38.Avants, B.B., Tustison, N. and Song, G. (2009). Advanced normalization tools (ANTS). Insight Journal 2, 1–35.
    
    

39. 39.Zhang, H., Liu, J., Zhu, Z. and Li, H. (2011). An automated and simple method for brain MR image extraction. BioMedical Engineering OnLine 10, 81.
    
    

40. 40.Jenkinson, M. and Smith, S. (2001). A global optimisation method for robust affine registration of brain images. Medical Image Analysis 5, 143–156.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1361-8415(01)00036-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11516708&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000169672400005&link_type=ISI) 

41. 41.Jenkinson, M., Bannister, P., Brady, M. and Smith, S. (2002). Improved optimization for the robust and accurate linear registration and motion correction of brain images. NeuroImage 17, 825–841.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1053-8119(02)91132-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12377157&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000178642000027&link_type=ISI) 

42. 42.Fischl, B., Salat, D.H., Busa, E., Albert, M., Dieterich, M., Haselgrove, C., van der Kouwe, A., Killiany, R., Kennedy, D., Klaveness, S., Montillo, A., Makris, N., Rosen, B. and Dale, A.M. (2002). Whole brain segmentation: Automated labeling of neuroanatomical structures in the human brain. Neuron 33, 341–355.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0896-6273(02)00569-X&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11832223&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000173643200006&link_type=ISI) 

43. 43.Sled, J.G., Zijdenbos, A.P. and Evans, A.C. (1998). A nonparametric method for automatic correction of intensity nonuniformity in MRI data. Institute of Electrical & Electronics Engineers Transactions on Medical Imaging 17, 87–97.
    
    

44. 44.Dale, A.M., Fischl, B. and Sereno, M.I. (1999). Cortical Surface-Based Analysis: I. Segmentation and Surface Reconstruction. NeuroImage 9, 179–194.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/nimg.1998.0395&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9931268&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000078608900001&link_type=ISI) 

45. 45.Fischl, B., Sereno, M.I. and Dale, A.M. (1999). Cortical surface-based analysis. II: Inflation, flattening, and a surface-based coordinate system. NeuroImage 9, 195–207.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/nimg.1998.0396&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9931269&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000078608900002&link_type=ISI) 

46. 46.Dhollander, T., Mito, R., Raffelt, D. and Connelly, A. (2019). Improved white matter response function estimation for 3-tissue constrained spherical deconvolution. In: International Society of Magnetic Resonance in Medicine: Montréal, Québec, Canada, pps. 555.
    
    

47. 47.Dhollander, T., Raffelt, D. and Connelly, A. (2016). Unsupervised 3-tissue response function estimation from single-shell or multi-shell diffusion MR data without a co-registered T1 image In: ISMRM Workshop on Breaking the Barriers of Diffusion MRI Lisbon, Portugal, pps. 5.
    
    

48. 48.Raffelt, D., Dhollander, T., Tournier, J.-D., Tabbara, R., Smith, R.E., Pierre, E. and Connelly, A. (2017). Bias field correction and intensity normalisation for quantitative analysis of apparent fibre density. In: International Society of Magnetic Resonance in Medicine: Honolulu, pps. 3541.
    
    

49. 49.Watts, D.J. and Strogatz, S.H. (1998). Collective dynamics of ‘small-world’ networks. Nature 393, 440–442.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/30918&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9623998&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000074020000035&link_type=ISI) 

50. 50.Hellyer, P.J., Scott, G., Shanahan, M., Sharp, D.J. and Leech, R. (2015). Cognitive flexibility through metastable neural dynamics is disrupted by damage to the structural connectome. J. Neurosci. 35, 9050–9063.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjEwOiIzNS8yNC85MDUwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDMvMDQvMjAyMi4wMy4wMi4yMjI3MTY1NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

51. 51.Latora, V. and Marchiori, M. (2001). Efficient behavior of small-world networks. Physical Review Letters 87, 198701.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1103/PhysRevLett.87.198701&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11690461&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 

52. 52.Boguna, M., Krioukov, D. and Claffy, K.C. (2009). Navigability of complex networks. Nature Physics 5, 74–80.
    
    

53. 53.Seguin, C., van den Heuvel, M.P. and Zalesky, A. (2018). Navigation of brain networks. Proceedings of the National Academy of Sciences 115, 6297–6302.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTE1LzI0LzYyOTciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMy8wNC8yMDIyLjAzLjAyLjIyMjcxNjU0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

54. 54.Wang, X., Seguin, C., Zalesky, A., Wong, W.-w., Chu, W.C.-w. and Tong, R.K.-y. (2019). Synchronization lag in post stroke: Relation to motor function and structural connectivity. Network Neuroscience 3, 1121–1140.
    
    

55. 55.Caeyenberghs, K., Leemans, A., De Decker, C., Heitger, M., Drijkoningen, D., Linden, C.V., Sunaert, S. and Swinnen, S.P. (2012). Brain connectivity and postural control in young traumatic brain injury patients: A diffusion MRI based network analysis. Neuroimage: Clinical 1, 106–115.
    
    

56. 56.Verhelst, H., Vander Linden, C., De Pauw, T., Vingerhoets, G. and Caeyenberghs, K. (2018). Impaired rich club and increased local connectivity in children with traumatic brain injury: Local support for the rich? Human Brain Mapping 39, 2800–2811.
    
    

57. 57.Jolly, A.E., Scott, G.T., Sharp, D.J. and Hampshire, A.H. (2020). Distinct patterns of structural damage underlie working memory and reasoning deficits after traumatic brain injury. Brain 143, 1158–1176.
    
    

58. 58.Caeyenberghs, K., Leemans, A., De Decker, C., Heitger, M., Drijkoningen, D., Vander Linden, C., Sunaert, S. and Swinnen, S.P. (2012). Brain connectivity and postural control in young traumatic brain injury patients: A diffusion MRI based network analysis. NeuroImage: Clinical 1, 106–115.
    
    

59. 59.Freeman, L.C. (1978). Centrality in social networks conceptual clarification. Social Networks 1, 215–239.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0378-8733(78)90021-7&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1979GL19400002&link_type=ISI) 

60. 60.Wills, P. and Meyer, F.G. (2020). Metrics for graph comparison: A practitioner’s guide. PloS One 15, e0228728–e0228728.
    
    

61. 61.Dichter, G.S., Sikich, L., Song, A., Voyvodic, J. and Bodfish, J.W. (2012). Functional neuroimaging of treatment effects in psychiatry: Methodological challenges and recommendations. International Journal of Neuroscience 122, 483–493.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22471393&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 

62. 62.Stoeckel, L.E., Garrison, K.A., Ghosh, S.S., Wighton, P., Hanlon, C.A., Gilman, J.M., Greer, S., Turk-Browne, N.B., deBettencourt, M.T. and Scheinost, D. (2014). Optimizing real time fMRI neurofeedback for therapeutic discovery and development. NeuroImage: Clinical 5, 245–255.
    
    

63. 63.Wing, B.H., Tucker, B.J., Fong, A.K. and Allen, M.D. (2017). Developing the standard of care for post-concussion treatment: Neuroimaging-guided rehabilitation of neurovascular coupling. The Open Neuroimaging Journal 11, 58.
    
    

64. 64.Meningher, I., Bernstein-Eliav, M., Rubovitch, V., Pick, C.G. and Tavor, I. (2020). Alterations in network connectivity after traumatic brain injury in mice. Journal of Neurotrauma 37.
    
    

65. 65.Osmanlioglu, Y., Alappatt, J.A., Parker, D., Kim, J. and Verma, R. (2019). A graph based similarity measure for assessing altered connectivity in traumatic brain injury. In: International MICCAI Brainlesion Workshop: Granada, Spain, pps. 189–198.
    
    

66. 66.Siegel, J.S., Shulman, G.L. and Corbetta, M. (2017). Measuring functional connectivity in stroke: Approaches and considerations. Journal of Cerebral Blood Flow & Metabolism 37, 2665–2678.
    
    

67. 67.Beelen, C., Phan, T.V., Wouters, J., Ghesquière, P. and Vandermosten, M. (2020). Investigating the added value of FreeSurfer’s manual editing procedure for the study of the reading network in a pediatric population. Frontiers in human neuroscience 14, 143–143.
    
    

68. 68.Dhollander, T., Clemente, A., Singh, M., Boonstra, F., Civier, O., Egorova, N., Enticott, P.G., Fuelscher, I., Gajamange, S., Genc, S., Gottlieb, E., Hyde, C., Imms, P., Kelly, C., Kirkovski, M., Kolbe, S., Liang, X., Malhotra, A., Poudel, G.R., Silk, T.J., Vaughan, D.N., Zanin, J., Raffelt, D. and Caeyenberghs, K. (2020). Fixel-based analysis of diffusion MRI: Methods, applications, challenges and opportunities. OSF Preprints.
    
    

69. 69.Olsen, A., Babikian, T., Bigler, E.D., Caeyenberghs, K., Conde, V., Dams-O’Connor, K., Dobryakova, E., Genova, H., Grafman, J., Håberg, A.K., Heggland, I., Hellstrøm, T., Hodges, C.B., Irimia, A., Jha, R.M., Johnson, P.K., Koliatsos, V.E., Levin, H., Li, L.M., Lindsey, H.M., Livny, A., Løvstad, M., Medaglia, J., Menon, D.K., Mondello, S., Monti, M.M., Newcombe, V.F.J., Petroni, A., Ponsford, J., Sharp, D., Spitz, G., Westlye, L.T., Thompson, P.M., Dennis, E.L., Tate, D.F., Wilde, E.A. and Hillary, F.G. (2020). Toward a global and reproducible science for brain imaging in neurotrauma: The ENIGMA adult moderate/severe traumatic brain injury working group. Brain Imaging & Behavior.
    
    

70. 70.van den Heuvel, M.P., Mandl, R.C., Stam, C.J., Kahn, R.S. and Hulshoff Pol, H.E. (2010). Aberrant frontal and temporal complex network structure in schizophrenia: A graph theoretical analysis. The Journal of Neuroscience 30, 15915–15926.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjExOiIzMC80Ny8xNTkxNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAzLzA0LzIwMjIuMDMuMDIuMjIyNzE2NTQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

71. 71.Hagmann, P., Cammoun, L., Gigandet, X., Meuli, R., Honey, C.J., Wedeen, V.J. and Sporns, O. (2008). Mapping the structural core of human cerebral cortex. PLoS Biology 6, e159.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pbio.0060159&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18597554&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom) 

72. 72.Andersson, J.L.R., Graham, M.S., Zsoldos, E. and Sotiropoulos, S.N. (2016). Incorporating outlier detection and replacement into a non-parametric framework for movement and distortion correction of diffusion MR images. NeuroImage 141, 556–572.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2016.06.058&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F04%2F2022.03.02.22271654.atom)

 [1]: /embed/graphic-7.gif