Longitudinal Changes in Diffusion Tensor Imaging Following Mild Traumatic Brain Injury and Correlation With Outcome
Bo Yin1,2 
Dan-Dong Li2 
Huan Huang3 Cheng-Hui Gu2 
Guang-Hui Bai3 Liu-Xun Hu2 Jin-Fei Zhuang4 Ming Zhang1*- 1Department of Medical Imaging, The First Affiliated Hospital of Xi’an Jiaotong University, Xi’an, China
- 2Department of Neurosurgery, The Second Affiliated Hospital and Yuying Children’s Hospital of Wenzhou Medical University, Wenzhou, China
- 3Department of Radiology, The Second Affiliated Hospital and Yuying Children’s Hospital of Wenzhou Medical University, Wenzhou, China
- 4Department of Rehabilitation Medicine, The Second Affiliated Hospital and Yuying Children’s Hospital of Wenzhou Medical University, Wenzhou, China
The chronic consequences of traumatic brain injury (TBI) may contribute to the increased risk for early cognitive decline and dementia, primarily due to diffusion axonal injury. Previous studies in mild TBI (mTBI) have been controversial in describing the white matter tract integrity changes occurring at acute and subacute post-injury. In this prospective longitudinal study, we aim to investigate the longitudinal changes of white matter (WM) using diffusion tensor imaging (DTI) and their correlations with neuropsychological tests. Thirty-three patients with subacute mTBI and 31 matched healthy controls were studied with an extensive imaging and clinical battery. Neuroimaging was obtained within 7 days post-injury for acute scans and repeated at 1 and 3 months post-injury. Using a region-of-interest-based approach, tract-based spatial statistics was used to conduct voxel-wise analysis on diffusion changes in mTBI and was compared to those of healthy matched controls, scanned during the same time period and rescanned with an interval similar to that of patients. We found decreased fractional anisotropy (FA) values in the left anterior limb of internal capsule (ALIC) and right inferior fronto-occipital fasciculus (IFOF) during the 7 days post-injury, which showed longitudinal evidence of recovery following 1 month post-injury. Increased FA values in these two tracts at 1 month post-injury were positively associated with better performance on cognitive information processing speed at initial assessment. By contrast, there were also some tracts (right anterior corona radiata, forceps major, and body of corpus callosum) exhibiting the continuing loss of integrity sustaining even beyond 3 months, which can predict the persisting post-concussion syndromes. Continuing loss of structural integrity in some tracts may contribute to the persistent post-concussion syndromes in mTBI patients, suggesting certain tracts providing an objective biomarker for tracking the pathological recovery process following mTBI.
Introduction
The United States Centers for Disease Control and Prevention estimates that mild traumatic brain injury (TBI) is experienced in 70–90% of TBI-related emergency department (ED) visits (Taylor et al., 2017). The diagnosis of mild TBI (mTBI) still lacks an interdisciplinary consensus regarding what constitutes an mTBI, and its determination has been largely epidemiological in nature. Since the majority of patients always obtained negative results from both CT and conventional MRI, there is still limited information about the neurobiological changes that are associated with clinical recovery.
Despite the negative findings in conventional imaging assessments, distinctive white matter tracts (WMTs), also known as diffuse axonal injury (DAI), after mTBI have been confirmed in both human biopsies (Bigler, 2010) and animal models (Mac Donald et al., 2007; Spain et al., 2010). Diffusion tensor imaging (DTI) has gained wide acceptance in clinical-research-based detection of abnormalities in mTBI patients (Basser and Pierpaoli, 1998; Arfanakis et al., 2002; Bigler, 2010). Water molecules in brain tissue are influenced by cellular structures and macromolecules by promoting anisotropy diffusion due to its various disparities in different orientations. Detection of water molecule diffusion fractional anisotropy (FA) can indicate the integrity of a cellular structure (Beaulieu, 2002). Although consistent observations indicate decreased anisotropic diffusion at more chronic periods post-injury, considerable debate remains on the direction (i.e., increased or decreased) of diffusion metrics on the acute to semi-acute phase (usually less than 2 months after injury) of mTBI (Arfanakis et al., 2002; Bazarian et al., 2007; Lipton et al., 2009; Chu et al., 2010; Mayer et al., 2010; Henry et al., 2011; Messé et al., 2012). Mayer et al. (2010) investigated 22 adult patients with mTBI within weeks post-injury and found increased FA noted in the genu, left superior coronal radiation (CR), and left uncinate fasciculus (UF), while decreased radial diffusivity (RD) was noted only in the genu, left CR, and left UF. Wilde et al. (2008) indicated the increased FA in the entire corpus callosum (CC) in a sample of 10 pediatric mTBI patients scanned between 1 and 6 days post-injury. Another study presented decreased FA in the splenium and posterior CC in 20 adult mTBI patients between 1 and 10 days post-injury (Inglese et al., 2005). Some studies have reported null findings during the semi-acute stages of mTBI. In a study, 12 adult mTBI patients were scanned within 3 months of injury (Rutgers et al., 2008). There were no significant differences between mTBI patients and controls in the CC.
The existing data on diffusion measurements in mTBI patients have key limitations regarding experimental design and analysis methods, motivating further investigation of injury patterns. Acute effects of injury, particularly edema, can confound such analysis. Longitudinal studies documenting serial changes in DTI metrics following mild TBI have yet to be performed (Meier et al., 2016; Churchill et al., 2017). One of the studies demonstrated that concussed athletes had significantly increased FA in the superior longitudinal fasciculus (SLF) than healthy athletes both at acute (within 1 week) and chronic periods (1 month). Further analyses indicated that this increase of FA in the SLF was driven by the decrease in RD rather than an increase in axial diffusivity (RAD; Meier et al., 2016). Another study found regions showing reduced average FA for athletes at acute concussion and at return-to-play relative to controls, including right corona radiata and bilaterally in posterior limbs of the internal capsule (Churchill et al., 2017). In addition, some researchers have enrolled small numbers of patients (e.g., <20) or had highly variable between-scan intervals, making the interpretation of between-subject variability difficult. Furthermore, most studies have lacked appropriate control groups for a longitudinal investigation (Meier et al., 2015). This is essential for modeling normal time-related changes in the brain measurements and neuropsychological ratings.
In the present study, DTI was collected to measure regional (i.e., voxel-wise) changes in FA values in the initial 7 days post-injury and the 1- and 3-month follow-up phase in a sample of mild TBI patients. To our knowledge, no studies have prospectively examined serial DTI changes during the normal course of recovery that typifies patients with mild TBI from these time course phase. To control the between-scan intervals, we confined this at a narrow interval. We mainly focused on the main time effect of the structural integrity of white matter (WM) fibers and its relationship with the initial neuropsychological measures. We hypothesized that the structural integrity loss in some tracts would normalize in patients with mild TBI as they transitioned from the acute to subacute injury phase and predicted that these deficits would relate to outcome measures, while other tracts will have persistent injury even beyond the first 3 months.
Materials and Methods
Participants
All consecutive patients with non-contrast head CT due to acute head trauma enrolled from the local ED formed the initial population. Inclusion criteria for all mild TBI patients were based on the World Health Organization’s Collaborating Centre for Neurotrauma Task Force (Borg et al., 2004): (i) Glasgow Coma Scale (GCS) score of 13–15 on presentation to the ED; (ii) one or more/any of the following: loss of consciousness (LOC) for less than 30 min, posttraumatic amnesia (PTA) for 24 h or less, and/or other transient neurological abnormalities such as focal signs, seizure, and intracranial lesion not requiring surgery; (iii) within 1 week after onset of a mild TBI (concussion); and (iv) age 16 years or older. Mild TBI patients were excluded for the following: (1) a history of a previous brain injury, neurological disease, long-standing psychiatric condition, or concurrent substance or alcohol abuse; (2) a structural abnormality on neuroimaging (CT and MRI); (3) intubation and/or presence of a skull fracture and administration of sedatives; (4) the manifestation of mild TBI due to medications by other injuries (e.g., systemic injuries, facial injuries, or spinal cord injury); (5) other problems (e.g., psychological trauma, language barrier, or coexisting medical conditions); and (6) TBI caused by penetrating craniocerebral injury. All of the patients were also screened for litigation to avoid bias in the assessment of neuropsychological tests (NPTs).
Healthy controls (HCs) with no history of neurological or psychiatric disorder were also recruited via the local imaging research facilities. Thirty-three patients with mTBI and 31 sex-, age-, and education-matched HC participated in the study. All participants were right-handed according to the Edinburgh Handedness Inventory (Oldfield, 1971).
Standard Protocol
All the subjects gave written informed consent in person approved by a local institutional review board; the research procedures were approved by the Ethical Committee of the Second Affiliated Hospital of Wenzhou Medical University and conducted in accordance with the Declaration of Helsinki. MRI and neuropsychological assessments were performed on patients at both initial visit within 7 days post-injury (T1; median, 2 days; range, 0–5 days) and follow-up at both 1 month (T2; median, 37 days; range, 27–35 days) and 3 months (T3; median, 104 days; range 85–105 days). NPTs were performed within 48 h of MRI. The HCs received the first scanning within the same time range as the patients and follow-up 1 month (median, 35 days; range, 28–43 days) and 3 months (median, 108 days; range, 92–111 days) later. DTI and neuropsychological assessments were acquired in the full dataset of 33 patients and 31 HC at all three time points.
Image Acquisition
A non-contrast CT scan was performed on all consecutive patients following acute head injury with a 64-row CT scanner (GE, Lightspeed VCT). The MRI scans were acquired with the use of a 3.0-T MRI scanner (GE 750). A custom-built head holder was used to prevent head movements. The MRI protocol involved the following: the high-resolution T1-weighted 3D MPRAGE sequence [echo time (TE) = 3.17 ms, repetition time (TR) = 8.15 ms, flip angle = 9°, slice thickness = 1 mm, field of view (FOV) = 256 mm × 256 mm, matrix size = 256 × 256], axial FLAIR (TR = 9,000 ms, TE = 95 ms, flip angle = 150°, thickness = 5 mm, slices = 20, FOV = 240 mm × 240 mm, matrix size = 173 × 256), axial susceptibility weighted imaging (TR = 37.8 ms, TE = 25 ms, flip angle = 15°, thickness = 2 mm, slices = 70, FOV = 230 mm × 230 mm, matrix size = 512 × 512), and diffusion-weighted imaging (TR = 7,300 ms, TE = 99 ms, flip angle = 90°, thickness = 3 mm, slices = 50, FOV = 256 mm × 256 mm, matrix size = 128 × 128, two averages, voxel size = 2 mm × 2 mm × 3 mm). DTI scan (b = 1,000 s/mm2) were acquired with 30 diffusion gradient orientations and the b = 0 repeated two times. The presence of focal lesions and cerebral microbleeds was independently determined by experienced clinical neuroradiologists (with 9 and 10 years’ experience) who assessed multiple modalities of neuroimaging data acquired at baseline [T1-flair, T2-flair, T2, susceptibility weighted imaging (SWI)]. Any disagreement between these two observers was resolved by consensus. None of the patients were with visible contusion lesions using conventional neuroimaging techniques or exhibited cerebral microbleeds on SWI.
DTI Data Analysis
Preprocessing of the raw DTI data was performed using FSL software1. Motion and eddy current corrections were carried out by means of affine registration to the reference volume, and the corrected data were brain-extracted using FSL’s Brain Extraction Tool (BET). FA images were then created by fitting a tensor model to the raw diffusion data using the FMRIB Diffusion Toolbox (FDT). Voxel-wise statistical analysis of the FA data was performed using tract-based spatial statistics (TBSS), a tool included in the FSL software to examine WM diffusion in a whole-brain voxel-based manner. All the subjects’ diffusion metrics were subsequently aligned into a stereotactic coordinate system with the MNI152 template, using the FMRIB non-linear registration tool (FNIRT) with B-spline representation of the registration warp field. Next, all FA maps were averaged to produce a group mean FA image. The resultant mean FA map was then thinned to create a mean FA skeleton, representing the centers of all white-matter tracts in both study groups with a threshold FA value of 0.2. Statistical analysis was performed using the “randomize” command in FSL. The number of permutations was set to 10,000, and correction for multiple comparisons was achieved using threshold-free cluster enhancement (TFCE) with a family-wise error (FWE) rate of P < 0.05.
Effect of White Matter Hyperintensities on Diffusion Tensor Imaging Analysis
Automated lesion segmentation was performed using the T1-weighted and T2-FLAIR image modalities on LST (Lesion Segmentation Tool) software package implemented in SPM12. LST determines gray matter (GM), WM, and cerebrospinal fluid (CSF) segmentations from T1-weighted images and computes the FLAIR intensity distributions of these tissue classes. The amount of “hyperintensity” of each voxel in terms of distance from the mean intensity of the WM, GM, and CSF distributions in the FLAIR image is crucial for defining a conservative lesion belief map (obtained by thresholding the GM belief map, initial threshold kappa was set to κ = 0.15) and a liberal lesion belief map (consisting of the sum of the three lesion belief maps). Lesion growing is then performed iteratively between the conservative and the liberal belief maps, until no more voxels are added to the lesions (Valverde et al., 2015). Finally, the subject-specific WM mask, generated by FAST in FSL (Zhang et al., 2001), was used to calculate the volume of WM hyperintensity based on the lesion belief map. Then, we tested whether the between-group differences were influenced by the individual WM hyperintensities.
Neuropsychological Tests
Comprehensive NPTs were assessed: (i) Trail-making test part A (TMT-A) (Arnett and Labovitz, 1995) and WAIS-III Digit Symbol Coding score (DSC) (Wechsler et al., 1997) to examine cognitive information processing speed; (ii) forward digit span (FDS) and backward digit span (BDS) of the Wechsler Adult Intelligence Scale WAIS-III (Harman-Smith et al., 2013) to assess working memory; and (iii) verbal fluency test (VFT; Troyer et al., 1997) to assess verbal fluency including language ability, semantic memory, and executive function. Self-reported symptomatology included the following: the Rivermead post-concussion symptoms questionnaire (RPCS; King et al., 1995), Posttraumatic stress disorder Checklist – Civilian version (PCL-C) (Bliese et al., 2008), Beck Depression Inventory 2nd edition BDI-II (beck; Beck et al., 1996), Fatigue Severity Scale (FSS; Krupp et al., 1989), and Insomnia Severity Index (ISI; Bastien et al., 2001). We have conducted the “reliable change indices” analysis for all of NPTs in the present study. This method provides an estimate of the probability that a given difference score would not be obtained by chance; that is, the score would not be due to measurement error (Iverson, 2001; Baxendale and Thompson, 2005). The standard error of difference (Sdiff) provides the clinician with an estimate of possible measurement error relating to test–retest scores.
Statistical Analyses
SPSS 19.0 software was used (IBM Corp., Armonk, NY, United States). The Shapiro–Wilk W-test was used to test for normality distribution of all continuous variables. The independent two-sample t-test and the Mann–Whitney test were used to compare group differences based on data normality, respectively. Chi-square analyses were applied to assess categorical variables. Effect sizes (Cohen’s d) were computed to demonstrate the magnitude of observed differences. Longitudinal analysis, via repeated-measures ANOVA, was conducted to examine changes in both neuropsychological assessments and imaging FA metric as a function of recovery. This analysis only included the group and time as terms with no further control variables. Spearman’s correlation coefficient was used to examine the association between FA metric of WMTs and neuropsychology test scores. The significance level was adjusted by using the Bonferroni correction with P < 0.05. None of the patients were with visible contusion lesions using neuroimaging techniques. Three patients exhibited cerebral microbleeds on SWI. All of the following analyses were also conducted in patients without microbleeds to investigate the influence of the presence of microbleeds on final results. We have also conducted the between-group difference in the TBSS analysis and used the individual WM hyperintensities observed in the T2-FLAIR sequence as regressors.
Results
Demographic and Initial Neuropsychological Data
There were no significant differences in the age (Cohen’s d = 0.01, P = 0.956), years of education (Cohen’s d = -0.49, P = 0.083), or gender (Cohen’s d = 0.16, P = 0.36) between the mTBI and control groups. The mTBI patients’ NPT differed significantly from controls on a number of tests at initial admission (RPCS, PCL-C, beck, and ISI, all for Cohen’s d > 0.9 and P < 0.005); however, some tests had no significant difference (TMT-A, DSC, FDS, BDS, VFT, and FSS all for P > 0.05). All demographic and clinical characteristics for patients with mild TBI and HC were presented in Table 1.
Table 1. Demographic and neuropsychological data for mTBI and control participants.
Neuropsychological Performance and Changes Over Time
Our results show that the majority of NPTs were sorted to the “most reliable” group, with the HCs having the most reliable change indices, given the relatively small practice effects (Supplementary Figure S1). However, only the Trail-making Test A presented the practice effects across the three time points. Thus, we were more cautious in explaining the corresponding changes in the Trail-making Test A in patients. For the tests with reliable change indices, the percentages of patients who demonstrated a significant improvement or decline were illustrated in Supplemetary Figure S2. For the DCS, about one-third of patients deteriorated and another one-third of patients improved between initial acute phase (T1) and follow-up 1 month post-injury (T2). For the FDS and BDS, one-third of patients presented improvement between T1 and T2. For the VFT, 20% of patients deteriorated and 15% of patients improved between T1 and T2.
Longitudinal analyses were first conducted to examine changes in the NPT as a function of recovery. A main effect of time was only observed in the DSC [F(2,60) = 4.14, P = 0.018]. For the DSC, scores were significantly improved at T2 and T3 relative to T1 (P = 0.003 and P < 0.001) and differences became non-significant at T2 and T3 following multiple comparison correction (P > 0.1). The main effect of group was significant for the RPCS, PCL-C, and beck (all for P < 0.05). Scores on the RPCS, PCL-C, and beck were significantly higher than HCs at T1 (P < 0.001), T2 (P < 0.005), and T3 (P < 0.002) following multiple comparison correction.
Diffusion Metrics Differences and Recovery
There was no significant interaction effect (group × time) for the WMTs. A significant main effect of time was located in the left anterior limb of internal capsule (ALIC) [F(1,61) = 9.49; P = 0.0023], right inferior fronto-occipital fasciculus (IFOF) [F(1,61) = 9.33; P = 0.0025], and body of CC [F(1,61) = 8.94; P = 0.0031] (Figure 1). For these three tracts, post hoc analyses with Bonferroni correction revealed that FA values were significantly lower at T1 (P < 0.001) and T2 (P < 0.005) relative to T3 and presented salient difference between T1 and T2 (P < 0.01). Mild TBI data from the left ALIC, right IFOF, and body of CC were next compared with those of HCs using post hoc analyses with Bonferroni correction. We found significant differences in the left ALIC and right IFOF FA values in patients relative to controls only at T1 (P < 0.005) but presented no significant difference at both T2 (P > 0.1 for all) and T3 (P > 0.08 for all) compared with HCs. The body of CC exhibited persistent loss of integrity at T1, T2, and T3 (P < 0.001) compared to HCs.
Figure 1. Main effect of time and post hoc results in the fractional anisotropy (FA). (A) FA showed significant main effect of time in the body of CC, left anterior limb of internal capsule (ALIC), and right inferior fronto-occipital fasciculus (IFOF). (B) For these three tracts, post hoc analyses with Bonferroni correction revealed that FA values were significantly lower at T1 (P < 0.001) and T2 (P < 0.005) relative to the T3, and presented salient difference between T1 and T2 (P < 0.005). (C) Significant differences in the left ALIC and right IFOF FA values in patients relative to controls only at T1 (P < 0.005) but returned to the normal level at both T2 (P > 0.1 for all) and T3 (P > 0.08 for all). The body of CC exhibited the continuing loss of integrity at T1, T2, and T3 (P < 0.001) compared to healthy controls (HCs). Error bars represent 95% confidence intervals. ∗Significant Bonferroni corrected for multiple comparisons (∗P < 0.005; ∗∗P < 0.001).
Significant effect of group was presented in the right anterior corona radiata (ACR), left posterior thalamic radiation (PTR), forceps major, body of CC, and splenium of CC (Figure 2). Post hoc analyses with Bonferroni correction further revealed that FA values in the left PTR, body of CC, and splenium of CC were reduced at T1, T2, and T3 relative to the HCs (P < 0.005 for all). The right ACR presented loss of integrity only at T3 (P < 0.001) and not at T1 and T2. For the forceps major, the reduced FA only presented at T2 (P = 0.0012) and T3 (P < 0.001) and not at T1. The combination of multiple parameters (MD, RD, and AD) from DTI was very helpful to understand the possible physiopathology mechanism following mild TBI. We also supplemented other diffusion metrics that presented reduced FA in mild TBI patients, compared with controls (Supplementary Table S1). Our findings showed that most between-group differences were not influenced by the individual WM hyperintensities. Only the between-group difference in the left PTR at 1 month post-injury following mild TBI disappeared after regressing for the WM hyperintensities (Supplementary Table S2).
Figure 2. Main effect of group and post hoc results in recovery of white matter tracts (WMTs). (A) Significant effect of group was presented in the right ACR, left posterior thalamic radiation (PTR), forceps major, body of corpus callosum (CC), and splenium of CC. (B) Post hoc analyses with Bonferroni correction further revealed that FA values in the left PTR, body of CC, and splenium of CC were reduced at T1, T2, and T3 relative to the HCs (P < 0.005 for all). The right ACR presents only loss of integrity at the later phase T3 rather than T1 and T2 (P < 0.001). For the forceps major, the reduced FA only presented at later phase at T2 (P = 0.0012) and T3 (P < 0.001) but not at T1. ∗Significant Bonferroni corrected for multiple comparisons (∗P < 0.005; ∗∗P < 0.001).
Relationship Between Outcome and Other Variables
The purpose of correlation analysis focused on the assessment of the relationship between the altered integrity of WMTs and neuropsychological performances within each time point and across time points. We were mainly concerned about both neuropsychological testing and WMTs showing the main effect of time. For the NPTs, a main effect of time was only observed in the DSC [F(2,60) = 4.14, P = 0.018]. For the fibers, a significant main effect of time was located only in the left ALIC [F(1,61) = 9.49; P = 0.0023], right IFOF [F(1,61) = 9.33; P = 0.0025], and body of CC [F(1,61) = 8.94; P = 0.0031]. The following analysis was conducted. Firstly, a correlation analysis between the clinical measure (DSC) and these three tracts for each time point was carried out. The FA value of these three tracts presented no significant associations with the DSC at any specific time points (all for P > 0.2). Additional analysis found that there was no significant correlation between altered DSC scores and FA values between T1 and T2, between T1 and T3, and between T2 and T3 in these three fibers (all for P > 0.2). Finally, considering that the left ALIC and right IFOF showed significant between-group differences at initial assessment but showed no significant difference at T2 and T3, we wanted to explore whether there was any association of the baseline DSC performance with the WM recovery in these two tracts. Exploratory analyses (Bonferroni correction at P < 0.05) were thus performed to assess the relationship between T1 clinical measure (DSC) with the left ALIC and right IFOF tracts at both T2 and T3, respectively. Results indicated that initial DSC performance can predict better recovery in the structural integrity of the left ALIC (r = 0.48, P = 0.005) and right IFOF (r = 0.54, P = 0.0012) at T2 after Bonferroni correction (P < 0.05/2, Figure 3). These findings may infer that an individual patient with better performance on the DSC at the acute phase can obtain better recovery in these two WMTs 1 month post-injury following mild TBI.
Figure 3. Correlation between the FA of WMTs with clinical outcomes. Patient’s cognitive information processing speed performance (rated by digit symbol coding, DSC) at T1 was positively related with FA value of the left anterior limb of internal capsule (ALIC) at T2 (r = 0.48, P < 0.005). Correlations of the DSC at T1 and FA value in the right IFOF at T2 (r = 0.5, P = 0.003) also showed significant positive values after Bonferroni correction.
We were also concerned with both neuropsychological testing and WMTs showing the main effect of group. For the NPTs, a main effect of time was only observed in the RPCS, PCL-C, and Beck depression (all for P < 0.05). For the fibers, a significant main effect of time was only located in the right ACR, left PTR, forceps major, body of CC, and splenium of CC. For each time point, correlation analyses showed that the FA value of these five tracts presented no significant correlations with these syndromes at any specific time point (P > 0.05). We also want to explore whether the baseline assessments of these NPTs were associated with the later integrity of WMTs at T2 and T3 follow-up phases. Results showed that the significance only existed at T3 (Figure 4). Results indicated a negative relationship of the right ACR with both baseline post-concussion severity (rated by RPCS, r = -0.53, P = 0.0017) and posttraumatic stress (PCL-C) (r = -0.45, P = 0.0095) but did not survive correction for multiple comparisons (P < 0.05/3 × 5). Correlations of the forceps major with RPCS and Beck depression (r = -0.5, P = 0.003; r = -0.57, P < 0.001) also showed significant negative values after Bonferroni correction. The loss of structural integrity in the forceps major also presented more complaints in the baseline PCL-C but did not survive correction for multiple comparisons (r = -0.46, P = 0.007). These findings inferred that patients with more complaints on the post-concussion syndromes at baseline have more loss of structural integrity in the right ACR and forceps major. Other correlations between baseline syndromes and later WM injury did not obtain significance.
Figure 4. Correlation between the FA of WMTs with clinical outcomes. At T3, the results indicated a negative relationship of the right anterior corona radiata (ACR) with both the post-concussion severity (rated by RPCS, r = -0.48, P < 0.005) and posttraumatic stress (PCL-C; r = -0.49, P = 0.0037). Correlations with the forceps major and RPCS, PCL-C, and Beck depression (r = -0.5, P = 0.003; r = -0.49, P = 0.0037; r = -0.57, P < 0.001) also showed significant negative values after Bonferroni correction.
Discussion
To our knowledge, this is the first study to investigate longitudinal changes in WMTs from the acute through subacute phases to 3 months post-injury following mTBI. In a homogeneous sample of mTBI, we found decreased FA values in the left ALIC and right IFOF during the 7 days post-injury, which showed longitudinal evidence of recovery following 1 month post-injury. Increased FA values in these two tracts at 1 month post-injury were positively associated with better performance on cognitive information processing speed at initial assessment. By contrast, there were also some tracts (right ACR, forceps major, and body of CC) exhibiting the persistent loss of integrity sustaining even beyond 3 months, which can also predict the persisting post-concussion syndromes, posttraumatic stress, and depression even beyond 3 months post-injury. Crucially, the persistent loss of structural integrity in these tracts may contribute to the persistent post-concussion syndromes in mild TBI patients.
Diffusion tensor imaging is widely used in vivo studies on WM injury of TBI (Shenton et al., 2012). Among severe and chronic TBI patients, DTI scans can determine the altered FA values in fiber tracts, wherein the mechanism had been illustrated in previous studies with animal models or humans and was considered as the result of diffuse degenerative changes such as Wallerian degeneration, axonal collapse, and myelin degeneration (Albensi et al., 2000; Arfanakis et al., 2002; Inglese et al., 2005; Li et al., 2011). However, due to the preliminary nature of these studies, there is still controversy regarding the DTI metric changes in mild TBI patients. Indeed, there are several differences regarding the pathological changes between chronic and acute and subacute TBI patients (Povlishock and Katz, 2005; Newcombe et al., 2007). The cellular injury in the acute phase is volatile and may, in cases, stay the same, ameliorate gradually, or deteriorate; thus, it is thought to be a very dynamic process (Albensi et al., 2000; Povlishock and Katz, 2005). In addition to this characteristic of axon injury during the acute and subacute phases, the injuries of mTBI patients also differ from those of moderate to severe TBI patients. Until now, it is still unclear how these different neuropathological changes (individually or combined) affect the DTI test results. Indeed, many studies have reported differences in axial diffusivity and have been inconsistent, with some groups finding increases (Sidaros et al., 2008; Tasker et al., 2010), decreases (Henry et al., 2011), and no change (Mac Donald et al., 2007), and some regions showed axial diffusivity increases while others showed decreases (Farbota et al., 2012). Therefore, it is premature to suggest a physiopathology mechanism based on imaging data alone. This is a limitation for our study as well as most other current studies. These pathological changes occur at different points of time and have different effects on axons and neurons, which may explain the variable results of the FA values. To overcome these variabilities, all of the participants in the present study were evaluated using both neurological and imaging assessments at a narrow time interval.
The acute phase is thought to be a very dynamic process, and previous studies present bidirectional changes in FA of mTBI patients (Sidaros et al., 2008; Tasker et al., 2010; Henry et al., 2011). Histopathological studies with animal and human tissues show the acute stage of mTBI, and the pathological change of axon microstructure is a dynamic process with multiple variable effects. These inconsistencies may be partly due to the highly variable post-injury days for acute phase, such as 1–6 days (Chu et al., 2010; Yallampalli et al., 2013), within 72 h (Bazarian et al., 2007; Toth et al., 2013), and within 24 h (Arfanakis et al., 2002), or limited samples sizes, such as 5 patients in Arfanakis’s study (Arfanakis et al., 2002), 6 in Bazarian’s study (Bazarian et al., 2007), 10 in Chu’s study (Chu et al., 2010), and 11 in Yallampalli’s study (Yallampalli et al., 2013), making the interpretation of between-subject variability difficult. In the present study, decreased FA values were observed in the left ALIC and right IFOF within acute 7 days post-injury and no bidirectional changes were detected. The result presented in this report highlights the need of strictly standardized image acquisition time points for mTBI patients.
For mTBI injuries, partial pathological changes are usually transient, but may persist in some cases. The pathological changes at 1 month post-injury are, however, not the research focus in previous studies. We found that FA values of the right left ALIC and right IFOF in patients were firstly reduced at the initial acute phase but presented no significant difference after 1 month post-injury. Previously, Toth et al. (2013) found that the mTBI patients showed decreased FA values in nearly all the WMTs 72 h after injury; however, after 1 month of injury, the FA values remained decreased in some tracts. Similarly, in a homogeneous sample of collegiate athletes with sports-related mTBI, Meier et al. (2015) observed decreased CBF in the right dorsal middle insula cortex and superior temporal sulcus during the first week post-injury, which showed recovery at 1 month post-injury. Additionally, in athletes with concussion, Vagnozzi et al. (2010) demonstrated that a reduction in the ratios of N-acetylaspartate to creatine and choline is pronounced at 3 days post-concussion but fully recovered by 1 month. The recovery of FA could reflect reorganization within the WM, due to axonal recovery or even regrowth (Sidaros et al., 2008). We found some support for this proposal by specifically examining the left ALIC and right IFOF, where the increased FA values at 1 month post-injury were positively associated with better performance on cognitive information processing speed at initial assessment. The ALIC connects the thalamus and prefrontal cortex, which conveys cognitive and corticothalamic fibers, and the right IFOF connects occipital, temporal, and frontal lobes, which has an important role in reading, attention, visual perception, processing and memory, and language; thus, such connections likely influence higher-order cognition (Martino et al., 2010). This result therefore provides support for the performance that 1 month post-injury is an important phase when the normalization of axial diffusivity in some WM occurs. Although statistical unification of neuroimaging data obtained acutely and subacutely (e.g., 1 month) after injury does not seem justifiable, these findings suggest that the MRI-detectable abnormality developing in the acute phase resolves dynamically in the first month after mild injury. Therefore, to determine whether the injuries obtained by complex and variable pathogenic mechanisms during the early stage might persist throughout the later stages of disease, 1 month post-injury was rather important.
By contrast, there were also some tracts (right ACR, forceps major, and body of CC) exhibiting the persistent loss of integrity sustaining beyond 3 months. Crucially, the correlations between FA and improved cognitive deficits at 3 months were positive: higher FA in right ACR and forceps major was associated with better outcomes. ACR connects to the internal capsule and forceps major connects to prefrontal and fronto-orbital regions, which may be related to perceptual and cognitive functions. Symptoms such as dizziness, fatigue, poor concentration, and memory problems experienced by chronic mTBI patients are thought to result from these lesions (Messé et al., 2011; Khong et al., 2016). Hence, the persistent loss of structural integrity in these tracts may contribute to the persistent post-concussion syndromes in mTBI patients, which can also predict the persisting post-concussion syndromes, posttraumatic stress, and depression even beyond 3 months post-injury. The identification of such predictive tracts (right ACR and forceps major) may help to better stratify patients early and to refine the concept of mTBI severity beyond traditional symptoms of alteration.
There were some limitations in the current study. Voxel-based analyses (i.e., TBSS) included the alignment of WMTs. In fact, such voxel-wise analysis is based on the assumption that clinically heterogeneous patients have a homogeneous (i.e., high degree of spatial overlap) spatial overlap pattern of WM abnormalities limited to detect the relative subtle diffuse “signals” influenced by the heterogeneity in injury location of mild TBI. Some studies adopted a new approach for quantifying diffusion abnormalities through a metric similar to lesion load (White et al., 2009; Ling et al., 2012). Specifically, clusters of abnormally high or low anisotropic diffusion were determined on a voxel-wise basis and then summed to represent total burden of distributed pathology. Further studies are needed to capture the individual-based “fine” diffuse patterns of WM abnormalities.
Conclusion
In conclusion, this study characterizes the dynamic variation rules of WM fibers within the first 3 months following mTBI. The results indicated that the pathological changes of fiber tracts at the acute phase were unstable and subject to a dynamic changing process. Therefore, we suggest a strict homogeneity of samples, and strictly standardized image acquisition time points for mTBI patients are required in subsequent studies. We found that FA values of the right left ALIC and right IFOF in patients were firstly reduced at the initial acute phase but presented no significant difference after 1 month post-injury compared with HCs. Hence, we suggest that 1 month post-injury is an important phase to determine whether the injuries obtained during the acute stage might persist throughout the chronic stages of mTBI. Finally, we also suggest that decreased FA values in right ACR and forceps major may help distinguish patients with posttraumatic stress disorder (PTSD) from mTBI patients at the acute and subacute stage.
Data Availability
All datasets generated for this study are included in the manuscript and/or the Supplementary Files.
Ethics Statement
All the subjects gave written, informed consent in person approved by a local institutional review board; the research procedures were approved by the Ethical Committee of The Second Affiliated Hospital of Wenzhou Medical University and conducted in accordance with the Declaration of Helsinki.
Author Contributions
BY and MZ contributed to the conception and design of the study. D-DL and HH organized the database. L-XH performed the statistical analysis. D-DL and C-HG wrote the first draft of the manuscript. J-FZ and G-HB wrote the sections of the manuscript. All authors contributed to the manuscript revision, read, and approved the submitted version.
Funding
This research was supported by the National Natural Science Foundation of China under Grant No. 81571640, the Zhejiang Natural Science Foundation (Grant No. LY15H090016), and Wenzhou Municipal Sci-Tech Bureau (Y20140577).
Conflict of Interest Statement
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncir.2019.00028/full#supplementary-material
Footnote
References
Albensi, B. C., Knoblach, S. M., Chew, B. G., O’Reilly, M. P., Faden, A. I., and Pekar, J. J. (2000). Diffusion and high resolution MRI of traumatic brain injury in rats: time course and correlation with histology. Exp. Neurol. 162, 61–72. doi: 10.1006/exnr.2000.7256
Arfanakis, K., Haughton, V. M., Carew, J. D., Rogers, B. P., Dempsey, R. J., and Meyerand, M. E. (2002). Diffusion tensor MR imaging in diffuse axonal injury. AJNR Am. J. Neuroradiol. 23, 794–802.
Arnett, J. A., and Labovitz, S. S. (1995). Effect of physical layout in performance of the trail making test. Psychol. Assess. 7:220. doi: 10.1037/1040-3590.7.2.220
Basser, P. J., and Pierpaoli, C. (1998). A simplified method to measure the diffusion tensor from seven MR images. Magn. Reson. Med. 39, 928–934. doi: 10.1002/mrm.1910390610
Bastien, C. H., Vallières, A., and Morin, C. M. (2001). Validation of the Insomnia Severity Index as an outcome measure for insomnia research. Sleep Med. 2, 297–307. doi: 10.1016/S1389-9457(00)00065-4
Baxendale, S., and Thompson, P. (2005). Defining meaningful postoperative change in epilepsy surgery patients: measuring the unmeasurable? Epilepsy Behav. 6, 207–211. doi: 10.1016/j.yebeh.2004.12.009
Bazarian, J. J., Zhong, J., Blyth, B., Zhu, T., Kavcic, V., and Peterson, D. (2007). Diffusion tensor imaging detects clinically important axonal damage after mild traumatic brain injury: a pilot study. J. Neurotrauma 24, 1447–1459. doi: 10.1089/neu.2007.0241
Beaulieu, C. (2002). The basis of anisotropic water diffusion in the nervous system—a technical review. NMR Biomed. 15, 435–455. doi: 10.1002/nbm.782
Beck, A. T., Steer, R. A., Ball, R., and Ranieri, W. F. (1996). Comparison of beck depression inventories-IA and -II in psychiatric outpatients. J. Pers. Assess. 67, 588–597. doi: 10.1207/s15327752jpa6703_13
Bigler, E. D. (2010). Neuroimaging in mild traumatic brain injury. Psychol. Inj. Law 3, 36–49. doi: 10.1007/s12207-010-9064-1
Bliese, P. D., Wright, K. M., Adler, A. B., Cabrera, O., Castro, C. A., and Hoge, C. W. (2008). Validating the primary care posttraumatic stress disorder screen and the posttraumatic stress disorder checklist with soldiers returning from combat. J. Consult. Clin. Psychol. 76, 272–281. doi: 10.1037/0022-006X.76.2.272
Borg, J., Holm, L., Cassidy, J. D., Peloso, P. M., Carroll, L. J., von Holst, H., et al. (2004). Diagnostic procedures in mild traumatic brain injury: results of the WHO collaborating centre task force on mild traumatic brain injury. J. Rehabil. Med. 36, 61–75. doi: 10.1080/16501960410023822
Chu, Z., Wilde, E. A., Hunter, J. V., McCauley, S. R., Bigler, E. D., Troyanskaya, M., et al. (2010). Voxel-based analysis of diffusion tensor imaging in mild traumatic brain injury in adolescents. AJNR Am. J. Neuroradiol. 31, 340–346. doi: 10.3174/ajnr.A1806
Churchill, N. W., Hutchison, M. G., Richards, D., Leung, G., Graham, S. J., and Schweizer, T. A. (2017). Neuroimaging of sport concussion: persistent alterations in brain structure and function at medical clearance. Sci. Rep. 7:8297. doi: 10.1038/s41598-017-07742-3
Farbota, K. D., Bendlin, B. B., Alexander, A. L., Rowley, H. A., Dempsey, R. J., and Johnson, S. C. (2012). Longitudinal diffusion tensor imaging and neuropsychological correlates in traumatic brain injury patients. Front. Hum. Neurosci. 6:160. doi: 10.3389/fnhum.2012.00160
Harman-Smith, Y. E., Mathias, J. L., Bowden, S. C., Rosenfeld, J. V., and Bigler, E. D. (2013). Wechsler adult intelligence scale—third edition profiles and their relationship to self-reported outcome following traumatic brain injury. J. Clin. Exp. Neuropsychol. 35, 785–798. doi: 10.1080/13803395.2013.824554
Henry, L. C., Tremblay, J., Tremblay, S., Lee, A., Brun, C., Lepore, N., et al. (2011). Acute and chronic changes in diffusivity measures after sports concussion. J. Neurotrauma 28, 2049–2059. doi: 10.1089/neu.2011.1836
Inglese, M., Makani, S., Johnson, G., Cohen, B. A., Silver, J. A., Gonen, O., et al. (2005). Diffuse axonal injury in mild traumatic brain injury: a diffusion tensor imaging study. J. Neurosurg. 103, 298–303. doi: 10.3171/jns.2005.103.2.0298
Iverson, G. L. (2001). Interpreting change on the WAIS-III/WMS-III in clinical samples. Arch. Clin. Neuropsychol. 16, 183–191. doi: 10.1093/arclin/16.2.183
Khong, E., Odenwald, N., Hashim, E., and Cusimano, M. D. (2016). Diffusion tensor imaging findings in post-concussion syndrome patients after mild traumatic brain injury: a systematic review. Front. Neurol. 7:156. doi: 10.3389/fneur.2016.00156
King, N. S., Crawford, S., Wenden, F. J., Moss, N. E., and Wade, D. T. (1995). The rivermead post concussion symptoms questionnaire: a measure of symptoms commonly experienced after head injury and its reliability. J. Neurol. 242, 587–592. doi: 10.1007/BF00868811
Krupp, L. B., LaRocca, N. G., Muir-Nash, J., and Steinberg, A. D. (1989). The fatigue severity scale: application to patients with multiple sclerosis and systemic lupus erythematosus. Arch. Neurol. 46, 1121–1123. doi: 10.1001/archneur.1989.00520460115022
Li, J., Li, X. Y., Feng, D. F., and Gu, L. (2011). Quantitative evaluation of microscopic injury with diffusion tensor imaging in a rat model of diffuse axonal injury. Eur. J. Neurosci. 33, 933–945. doi: 10.1111/j.1460-9568.2010.07573.x
Ling, J. M., Pena, A., Yeo, R. A., Merideth, F. L., Klimaj, S., Gasparovic, C., et al. (2012). Biomarkers of increased diffusion anisotropy in semi-acute mild traumatic brain injury: a longitudinal perspective. Brain 135, 1281–1292. doi: 10.1093/brain/aws073
Lipton, M. L., Gulko, E., Zimmerman, M. E., Friedman, B. W., Kim, M., Gellella, E., et al. (2009). Diffusion-tensor imaging implicates prefrontal axonal injury in executive function impairment following very mild traumatic brain injury. Radiology 252, 816–824. doi: 10.1148/radiol.2523081584
Mac Donald, C. L., Dikranian, K., Bayly, P., Holtzman, D., and Brody, D. (2007). Diffusion tensor imaging reliably detects experimental traumatic axonal injury and indicates approximate time of injury. J. Neurosci. 27, 11869–11876. doi: 10.1523/JNEUROSCI.3647-07.2007
Martino, J., Brogna, C., Robles, S. G., Vergani, F., and Duffau, H. (2010). Anatomic dissection of the inferior fronto-occipital fasciculus revisited in the lights of brain stimulation data. Cortex 46, 691–699. doi: 10.1016/j.cortex.2009.07.015
Mayer, A. R., Ling, J., Mannell, M. V., Gasparovic, C., Phillips, J. P., Doezema, D., et al. (2010). A prospective diffusion tensor imaging study in mild traumatic brain injury. Neurology 74, 643–650. doi: 10.1212/WNL.0b013e3181d0ccdd
Meier, T. B., Bellgowan, P. S., Singh, R., Kuplicki, R., Polanski, D. W., and Mayer, A. R. (2015). Recovery of cerebral blood flow following sports-related concussion. JAMA Neurol. 72, 530–538. doi: 10.1001/jamaneurol.2014.4778
Meier, T. B., Bergamino, M., Bellgowan, P. S., Teague, T. K., Ling, J. M., Jeromin, A., et al. (2016). Longitudinal assessment of white matter abnormalities following sports-related concussion. Hum. Brain Mapp. 37, 833–845. doi: 10.1002/hbm.23072
Messé, A., Caplain, S., Paradot, G., Garrigue, D., Mineo, J. F., Soto Ares, G., et al. (2011). Diffusion tensor imaging and white matter lesions at the subacute stage in mild traumatic brain injury with persistent neurobehavioral impairment. Hum. Brain Mapp. 32, 999–1011. doi: 10.1002/hbm.21092
Messé, A., Caplain, S., Pélégrini-Issac, M., Blancho, S., Montreuil, M., Lévy, R., et al. (2012). Structural integrity and postconcussion syndrome in mild traumatic brain injury patients. Brain Imaging Behav. 6, 283–292. doi: 10.1007/s11682-012-9159-2
Newcombe, V. F., Williams, G. B., Nortje, J., Bradley, P. G., Harding, S. G., Smielewski, P., et al. (2007). Analysis of acute traumatic axonal injury using diffusion tensor imaging. Br. J. Neurosurg. 21, 340–348. doi: 10.1080/02688690701400882
Oldfield, R. C. (1971). The assessment and analysis of handedness: the Edinburgh inventory. Neuropsychologia 9, 97–113. doi: 10.1016/0028-3932(71)90067-4
Povlishock, J. T., and Katz, D. I. (2005). Update of neuropathology and neurological recovery after traumatic brain injury. J. Head. Trauma. Rehabil. 20, 76–94. doi: 10.1097/00001199-200501000-00008
Rutgers, D. R., Fillard, P., Paradot, G., Tadié, M., Lasjaunias, P., and Ducreux, D. (2008). Diffusion tensor imaging characteristics of the corpus callosum in mild, moderate, and severe traumatic brain injury. AJNR Am. J. Neuroradiol. 29, 1730–1735. doi: 10.3174/ajnr.A1213
Shenton, M. E., Hamoda, H. M., Schneiderman, J. S., Bouix, S., Pasternak, O., Rathi, Y., et al. (2012). A review of magnetic resonance imaging and diffusion tensor imaging findings in mild traumatic brain injury. Brain Imaging Behav. 6, 137–192. doi: 10.1007/s11682-012-9156-5
Sidaros, A., Engberg, A. W., Sidaros, K., Liptrot, M. G., Herning, M., Petersen, P., et al. (2008). Diffusion tensor imaging during recovery from severe traumatic brain injury and relation to clinical outcome: a longitudinal study. Brain 131, 559–572. doi: 10.1093/brain/awm294
Spain, A., Daumas, S., Lifshitz, J., Rhodes, J., Andrews, P. J., Horsburgh, K., et al. (2010). Mild fluid percussion injury in mice produces evolving selective axonal pathology and cognitive deficits relevant to human brain injury. J. Neurotrauma 27, 1429–1438. doi: 10.1089/neu.2010.1288
Tasker, R. C., Westland, A. G., White, D. K., and Williams, G. B. (2010). Corpus callosum and inferior forebrain white matter microstructure are related to functional outcome from raised intracranial pressure in child traumatic brain injury. Dev. Neurosci. 32, 374–384. doi: 10.1159/000316806
Taylor, C. A., Bell, J. M., Breiding, M. J., and Xu, L. (2017). Traumatic brain injury–related emergency department visits, hospitalizations, and deaths—United States, 2007 and 2013. MMWR Surveill. Summ. 66, 1–16. doi: 10.15585/mmwr.ss6609a1
Toth, A., Kovacs, N., Perlaki, G., Orsi, G., Aradi, M., Komaromy, H., et al. (2013). Multi-modal magnetic resonance imaging in the acute and sub-acute phase of mild traumatic brain injury: can we see the difference? J. Neurotrauma 30, 2–10. doi: 10.1089/neu.2012.2486
Troyer, A. K., Moscovitch, M., and Winocur, G. (1997). Clustering and switching as two components of verbal fluency: evidence from younger and older healthy adults. Neuropsychology 11, 138–146. doi: 10.1037/0894-4105.11.1.138
Vagnozzi, R., Signoretti, S., Cristofori, L., Alessandrini, F., Floris, R., Isgrò, E., et al. (2010). Assessment of metabolic brain damage and recovery following mild traumatic brain injury: a multicentre, proton magnetic resonance spectroscopic study in concussed patients. Brain 133, 3232–3242. doi: 10.1093/brain/awq200
Valverde, S., Oliver, A., Roura, E., Pareto, D., Vilanova, J. C., Ramió-Torrentà, L., et al. (2015). Quantifying brain tissue volume in multiple sclerosis with automated lesion segmentation and filling. Neuroimage Clin. 9, 640–647. doi: 10.1016/j.nicl.2015.10.012
Wechsler, D., Coalson, D. L., and Raiford, S. E. (1997). WAIS-III: Wechsler Adult Intelligence Scale. San Antonio, TX: NCS Pearson.
White, T., Schmidt, M., and Karatekin, C. (2009). White matter ‘potholes’ in early-onset schizophrenia: a new approach to evaluate white matter microstructure using diffusion tensor imaging. Psychiatry Res. 174, 110–115. doi: 10.1016/j.pscychresns.2009.04.014
Wilde, E. A., McCauley, S. R., Hunter, J. V., Bigler, E. D., Chu, Z., Wang, Z. J., et al. (2008). Diffusion tensor imaging of acute mild traumatic brain injury in adolescents. Neurology 70, 948–955. doi: 10.1212/01.wnl.0000305961.68029.54
Yallampalli, R., Wilde, E. A., Bigler, E. D., McCauley, S. R., Hanten, G., Troyanskaya, M., et al. (2013). Acute white matter differences in the fornix following mild traumatic brain injury using diffusion tensor imaging. J. Neuroimaging 23, 224–227. doi: 10.1111/j.1552-6569.2010.00537.x
Keywords: mild traumatic brain injury, diffusion tensor imaging, longitudinal changes, fractional anisotropy, neuropsychological test
Citation: Yin B, Li D-D, Huang H, Gu C-H, Bai G-H, Hu L-X, Zhuang J-F and Zhang M (2019) Longitudinal Changes in Diffusion Tensor Imaging Following Mild Traumatic Brain Injury and Correlation With Outcome. Front. Neural Circuits 13:28. doi: 10.3389/fncir.2019.00028
Received: 05 March 2019; Accepted: 01 April 2019;
Published: 07 May 2019.
Edited by:
Tuo Zhang, Northwestern Polytechnical University, ChinaReviewed by:
Yibin Xi, Fourth Military Medical University, ChinaJun Liu, Central South University, China
Copyright © 2019 Yin, Li, Huang, Gu, Bai, Hu, Zhuang and Zhang. 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: Ming Zhang, zhangming01@mail.xjtu.edu.cn