# Machine Learning Classification of Mild Traumatic Brain Injury Using Whole-Brain Functional Activity: A Radiomics Analysis

Authors: Xiaoping Luo; Dezhao Lin; Shengwei Xia; Dongyu Wang; Xinmang Weng; Wenming Huang; Hongda Ye
Journal: Disease Markers
 (2021)

---

## Abstract

Objectives. To investigate the classification performance of support vector machine in mild traumatic brain injury (mTBI) from normal controls. Methods. Twenty-four mTBI patients (15 males and 9 females; mean age, 38.88±13.33 years) and 24 age and sex-matched normal controls (13 males and 11 females; mean age, 40.46±11.4 years) underwent resting-state functional MRI examination. Seven imaging parameters, including amplitude of low-frequency fluctuation (ALFF), fractional amplitude of low-frequency fluctuation (fALFF), regional homogeneity (ReHo), degree centrality (DC), voxel-mirrored homotopic connectivity (VMHC), long-range functional connectivity density (FCD), and short-range FCD, were entered into the classification model to distinguish the mTBI from normal controls. Results. The ability for any single imaging parameters to distinguish the two groups is lower than multiparameter combinations. The combination of ALFF, fALFF, DC, VMHC, and short-range FCD showed the best classification performance for distinguishing the two groups with optimal AUC value of 0.778, accuracy rate of 81.11%, sensitivity of 88%, and specificity of 75%. The brain regions with the highest contributions to this classification mainly include bilateral cerebellum, left orbitofrontal cortex, left cuneus, left temporal pole, right inferior occipital cortex, bilateral parietal lobe, and left supplementary motor area. Conclusions. Multiparameter combinations could improve the classification performance of mTBI from normal controls by using the brain regions associated with emotion and cognition.


---

## Body

## 1. Introduction

Traumatic brain injury (TBI), a major public health problem and a leading cause of disability, affects half the world’s population [1]. Approximately 70%-90% of TBI patients are mild TBI (mTBI), and 30-40% of whom cannot fully recover even at 6 months postinjury [1, 2]. Patients with mild head injury often manifest as dizziness, headache, and memory and attention deficit, which was considered to be associated with abnormal changes of brain networks [3]. Recently, functional and structural neuroimaging methods have been widely used to address the functional and morphological changes of mTBI [4–11]. Zhou et al. found abnormal functional connectivity within the default mode network in mTBI patients, which was associated with cognitive neurological dysfunction and posttraumatic symptoms (i.e., depression, anxiety, fatigue, and postconcussion syndrome) [12]. Nakamura et al. found that mTBI was associated with changes in the “small world” networks [13]. Zhan et al. found decreased ReHo value in the left insula, left pre-/postcentral gyrus, and left supramarginal gyrus in mTBI patients [14]. However, the potential neurobiological mechanism of the mTBI left unclear.Most current studies focus attentions on investigating group differences between two different labels (knowing the classes of all subject before statistics); however, group-based methods cannot classify different types for individual classification and are not sensitive for feature selection [15]. Support vector machine (SVM) classifier is an efficient and sensitive neuroimaging biological indicator for feature selection and classification. There is a growing application of the SVM algorithm into several diseases, such as insomnia [16, 17], epilepsy [15], and autistic spectrum disorder [18]. However, the mTBI has not been studied. Differences in brain regions in mTBI were not the same when we analyzed the between-group differences by different neuroimaging methods, which may be associated with the sensitivity of different methods in searching features (brain areas). Therefore, we hypothesized that the combination of different neuroimaging methods may improve the sensitivity for feature selection. To address these hypotheses, the present study is the first to apply the SVM algorithm to perform the classification for mTBI.


## 2. Materials and Methods

### 2.1. Subjects

This case-control study comprised 170 subjects from our hospital between May 2014 and May 2021, among whom a total of 146 subjects were excluded, including 139 subjects unmatched diagnosis with mTBI, 4 mTBI with more than 1.5 mm maximum translation inx, y, or z directions and/or 1.5 degree of motion rotation, and 3 mTBI with missing data. Finally, 24 patients with acute mTBI (15 males and 9 females; mean age, 38.88±13.33 years; mean years of education, 8.88±3.58 years; and mean time of postinjury, 3.58±3.28 days) and 24 age and sex-matched (13 males and 11 females; mean age, 40.46±11.4 years; and mean years of education, 8.54±3.41 years) healthy controls were included. All subjects were asked to complete the following questionnaires, including the Glasgow Coma Scale (GCS), Disability Rating Scale (DRS), Motor Assessment Scale (MAS), Agitated Behavior Scale (ABS), Hamilton Anxiety Scale (HAMA), Clinical Dementia Rating (CDR), Mini Mental State Examination (MMSE), Activates of Daily Living (ADL), and Beck Depression Inventory (BDI).Inclusion criteria for patients with acute mTBI were as follows: (a) have a diagnosis of mTBI within two weeks, (b) age between 18 and 65 years, (c) time of lack of consciousness less than 30 min, and (d) time of posttraumatic amnesia less than 24 hours. Exclusion criteria for patients with acute mTBI were as follows: (a) involvement in litigation, (b) a history of psychiatric disorders, (c) a history of addiction, and (d) a history of traumatic brain injury. This study was approved by the Human Research Ethics Committee in accordance with the Declaration of Helsinki, and written informed consent was obtained.


### 2.2. MRI Parameters

MRI data were acquired with a clinical 3-Tesla MRI scanner (Trio Tim, SIEMENS, Erlangen, Germany), including T1WI, T2WI, T2-FLAIR, high-resolution T1WI, functional MRI, and SWI. A total of 176 three-dimensional high-resolution anatomical T1-weighted volumes were acquired in a sagittal orientation (rapid-gradient-echo sequence,repetitiontime=1900 ms, echotime=2.26 ms, thickness=1.0 mm, matrix=256×256, and fieldofview=240mm×240mm). For functional images, a total of 250 volumes (Echo-Planar Imaging pulse sequence, 30 transverse slices, repetitiontime=2000 ms, echotime=40 ms, thickness=4.0 mm, matrix=64×64, fieldofview=240mm×240mm, and flip angle 90°) were acquired.


### 2.3. Data Processing

All functional MRI data preprocessing were performed with DPABI (version 2.1,http://rfmri.org/DPABI) toolbox. First, the first ten volumes were deleted, and the remaining volumes were converted their data format. The following steps of slice timing, head motion correction, spatial normalization, smooth (Gaussian kernel of 8×8×8 mm3), linear regression of possible spurious covariates, linearly detrended, and temporally band-pass filtered (0.01-0.1 Hz) were performed for data preprocessing. After the step of head motion correction, a “head motion scrubbing regressors” procedure was implemented, and the subjects who had more than 1.5 degree of motion rotation and/or 1.5 mm maximum translation in x, y, or z directions were excluded. Furthermore, the head motion effect was regressed out with Friston 24 head motion parameter model. During the step of spatial normalization, all data were spatially normalized to Montreal Neurological Institute (MNI) space and resampled at a resolution of 3×3×3 mm3.


### 2.4. Feature Selection and Binary Classification

We calculated seven MRI parameters, including ALFF, fALFF, ReHo, degree centrality, long-term FCD, short-term FCD, and VMHC. The maps of MRI parameters were segmented into 116 regions of interest (ROIs) using the automated anatomical labeling (AAL) atlas. The total of 812 features was extracted in the following classification with multivariate pattern analysis (MVPA).We used a LIBSVM toolbox (http://www.csie.ntu.edu.tw/~cjlin/libsvm/) to perform the classification, and a 5-fold cross-validation was used to validate the classification performance of the classifier. Permutation test was used to evaluate the probability of the classification performance for 5000 times randomly. The clusters of brain regions with higher than 70% of classification accuracy were considered as accuracies. The area under curve (AUC), sensitivity, and specificity of the classifier were quantified.


### 2.5. Statistical Analyses

Comparisons of demographic factors were performed using two-samplet-tests. Chi-square (χ2) test was used for categorical data. Statistical analysis was performed using IBM SPSS 21.0 version. Data are presented as mean±standarddeviation. All the quoted results are two-tailed values, and p<0.05 was considered as statistically significant.


## 2.1. Subjects

This case-control study comprised 170 subjects from our hospital between May 2014 and May 2021, among whom a total of 146 subjects were excluded, including 139 subjects unmatched diagnosis with mTBI, 4 mTBI with more than 1.5 mm maximum translation inx, y, or z directions and/or 1.5 degree of motion rotation, and 3 mTBI with missing data. Finally, 24 patients with acute mTBI (15 males and 9 females; mean age, 38.88±13.33 years; mean years of education, 8.88±3.58 years; and mean time of postinjury, 3.58±3.28 days) and 24 age and sex-matched (13 males and 11 females; mean age, 40.46±11.4 years; and mean years of education, 8.54±3.41 years) healthy controls were included. All subjects were asked to complete the following questionnaires, including the Glasgow Coma Scale (GCS), Disability Rating Scale (DRS), Motor Assessment Scale (MAS), Agitated Behavior Scale (ABS), Hamilton Anxiety Scale (HAMA), Clinical Dementia Rating (CDR), Mini Mental State Examination (MMSE), Activates of Daily Living (ADL), and Beck Depression Inventory (BDI).Inclusion criteria for patients with acute mTBI were as follows: (a) have a diagnosis of mTBI within two weeks, (b) age between 18 and 65 years, (c) time of lack of consciousness less than 30 min, and (d) time of posttraumatic amnesia less than 24 hours. Exclusion criteria for patients with acute mTBI were as follows: (a) involvement in litigation, (b) a history of psychiatric disorders, (c) a history of addiction, and (d) a history of traumatic brain injury. This study was approved by the Human Research Ethics Committee in accordance with the Declaration of Helsinki, and written informed consent was obtained.


## 2.2. MRI Parameters

MRI data were acquired with a clinical 3-Tesla MRI scanner (Trio Tim, SIEMENS, Erlangen, Germany), including T1WI, T2WI, T2-FLAIR, high-resolution T1WI, functional MRI, and SWI. A total of 176 three-dimensional high-resolution anatomical T1-weighted volumes were acquired in a sagittal orientation (rapid-gradient-echo sequence,repetitiontime=1900 ms, echotime=2.26 ms, thickness=1.0 mm, matrix=256×256, and fieldofview=240mm×240mm). For functional images, a total of 250 volumes (Echo-Planar Imaging pulse sequence, 30 transverse slices, repetitiontime=2000 ms, echotime=40 ms, thickness=4.0 mm, matrix=64×64, fieldofview=240mm×240mm, and flip angle 90°) were acquired.


## 2.3. Data Processing

All functional MRI data preprocessing were performed with DPABI (version 2.1,http://rfmri.org/DPABI) toolbox. First, the first ten volumes were deleted, and the remaining volumes were converted their data format. The following steps of slice timing, head motion correction, spatial normalization, smooth (Gaussian kernel of 8×8×8 mm3), linear regression of possible spurious covariates, linearly detrended, and temporally band-pass filtered (0.01-0.1 Hz) were performed for data preprocessing. After the step of head motion correction, a “head motion scrubbing regressors” procedure was implemented, and the subjects who had more than 1.5 degree of motion rotation and/or 1.5 mm maximum translation in x, y, or z directions were excluded. Furthermore, the head motion effect was regressed out with Friston 24 head motion parameter model. During the step of spatial normalization, all data were spatially normalized to Montreal Neurological Institute (MNI) space and resampled at a resolution of 3×3×3 mm3.


## 2.4. Feature Selection and Binary Classification

We calculated seven MRI parameters, including ALFF, fALFF, ReHo, degree centrality, long-term FCD, short-term FCD, and VMHC. The maps of MRI parameters were segmented into 116 regions of interest (ROIs) using the automated anatomical labeling (AAL) atlas. The total of 812 features was extracted in the following classification with multivariate pattern analysis (MVPA).We used a LIBSVM toolbox (http://www.csie.ntu.edu.tw/~cjlin/libsvm/) to perform the classification, and a 5-fold cross-validation was used to validate the classification performance of the classifier. Permutation test was used to evaluate the probability of the classification performance for 5000 times randomly. The clusters of brain regions with higher than 70% of classification accuracy were considered as accuracies. The area under curve (AUC), sensitivity, and specificity of the classifier were quantified.


## 2.5. Statistical Analyses

Comparisons of demographic factors were performed using two-samplet-tests. Chi-square (χ2) test was used for categorical data. Statistical analysis was performed using IBM SPSS 21.0 version. Data are presented as mean±standarddeviation. All the quoted results are two-tailed values, and p<0.05 was considered as statistically significant.


## 3. Results

### 3.1. Sample Characteristics

There were no significant differences in mean age (t=−0.442, p=0.66), sex (χ2=0.343, p=0.558), and educational level (t=0.33, p=0.743) between the healthy controls and patients with mTBI. Compared with healthy controls, patients with mTBI had higher HAMA score (t=5.077, p<0.001), ADL score (t=4.654, p<0.001), and BDI score (t=3.808, p=0.001), and a lower MMSE score (t=−2.284, p=0.03). The mean time between injury and MRI examination of patients with mTBI was 3.58±3.28 days. The mean GCS score, DRS score, MAS score, and ABS score in patients with mTBI were 14.42±0.88, 2.58±2.36, 44.38±5.86, and 14.42±0.78, respectively. The details are shown in Table 1.Table 1
Demographic and clinical features of patients with acute mTBI and healthy controls.
mTBIHealthy controlst valuep valueAge, years38.88±13.3340.46±11.40-0.4420.66Sex (male, female)24 (15, 9)24 (13, 11)0.3430.558Education, years8.88±3.588.54±3.410.3300.743Postinjury, days3.58±3.28N/AN/AN/AGCS14.42±0.88N/AN/AN/ADRS2.58±2.36N/AN/AN/AMAS44.38±5.86N/AN/AN/AABS14.42±0.78N/AN/AN/AHAMA3.83±3.610.08±0.285.077<0.001MMSE29.04±1.6329.83±0.20-2.2840.03ADL21.71±8.0714.04±6.064.654<0.001BDI1.58±1.910.08±0.283.8080.001


### 3.2. Classification Performance

First, we compared the classification performances of the seven MRI parameters and found they could not differentiate well between healthy controls and patients with mTBI (AUC:0.66±0.03, range, 0.61~0.69; accuracy rate: 66.4%±3.4%, range, 60.2%~70.9%; sensitivity: 64.1%±7.9%, range, 49.0%~75.0%; and specificity: 68.4%±5.6%, range, 61.0%~75.0%). Second, we combined these MRI parameters and found the features with the highest contributions to the classification to discriminate between mTBI and healthy controls. We found that the combination with ALFF, fALFF, DC, VMHC, and short-term FCD significantly reached up the classification accuracy, sensitivity, and specificity and received the highest classification performances among all combination with classification accuracy of 81.1% (p<0.001), sensitivity of 88.0% (p<0.001), and specificity of 75.0% (p<0.001) (Figure 1).Figure 1
Schematic diagram overview of machine learning classification framework. Note: this figure shows the classification of the combination of ALFF, fALFF, DC, VMHC, and short-term FCD in distinguishing the mTBI from the normal controls. The classification received the highest (a) AUC value, (b) classification accuracy, (c) sensitivity, and (d) specificity among all combination.
(a)(b)(c)(d)


### 3.3. Consensus Features and Region Weight

In this study, all consensus features were mapped to AAL116 template (116 brain regions), and each of the 116 brain regions was given a weight value which indicates the contribution to classification model. For the combination with ALFF, fALFF, DC, VMHC, and short-term FCD, Table2 shows the weight ranking of the 116 brain regions from highest to lowest.Table 2
Weight ranking of the 116 brain regions to the classification of the combination with ALFF, fALFF, DC, VMHC, and short-term FCD.
ROI weightVoxel sizeVermis_101.94834Cerebellum_9_R1.470156Cerebellum_10_L1.46540Cerebellum_9_L1.455158Frontal_Mid_Orb_L1.417224Frontal_Sup_Orb_L1.398280Cuneus_L1.353472Cerebellum_Crus2_R1.329539Cerebellum_7b_L1.26598Cerebellum_Crus2_L1.257543Frontal_Mid_Orb_L1.256273Temporal_pole_Mid_L1.222177Occipital_Inf_R1.220316Parietal_Sup_L1.136575Paracentral_lobule_R1.134221Frontal_Sup_R1.1181120Cuneus_R1.117416Cerebellum_Crus1_R1.103723Occipital_Inf_L1.103263Vermis_71.10254Calcarine_L1.076649Occipital_Sup_R1.074407Rectus_L1.070258Postcentral_R1.0551050Paracentral_lobule_L1.053340Precentral_R1.009941Parietal_Inf_R1.008397Occipital_Sup_L0.998373Cerebellum_10_R0.99737Cerebellum_7b_R0.97578Cerebellum_8_L0.961303Cerebellum_6_L0.960524Vermis_90.95950Temporal_Inf_R0.9591076Occipital_Mid_L0.954947Cerebellum_Crus1_L0.950725Lingual_L0.945662Supp_motor_area_L0.937630Frontal_Mid_R0.9311448Calcarine_R0.925528Temporal_Mid_L0.9221437Parietal_Sup_R0.921569Cerebellum_4_5_L0.919352Frontal_Sup_medial_L0.918847Lingual_R0.916683Angular_R0.916511Temporal_pole_Sup_R0.918325Cerebellum_6_R0.908532Precuneus_R0.908927Temporal_Sup_R0.907942Frontal_Sup_L0.890987Angular_L0.860341Precuneus_L0.8511008Cingulum_post_L0.850111Frontal_Inf_Tri_L0.848675Frontal_Mid_L0.8471323Temporal_pole_Sup_L0.828329Temporal_Sup_L0.825694Temporal_pole_Mid_R0.822264Cerebellum_8_R0.810298Cerebellum_3_R0.80165Occipital_Mid_R0.796578Supp_motor_area_R0.790695Vermis_4_50.788176Frontal_Sup_medial_R0.787589Frontal_Inf_Tri_R0.783560Supramarginal_R0.779562Precentral_L0.764931Heschl_R0.76360Frontal_Mid_Orb_R0.762296Frontal_Sup_Orb_R0.752296Frontal_Inf_Orb_L0.752504Cerebellum_3_L0.75042Supramarginal_L0.750357Fusiform_L0.747665Temporal_Inf_L0.745948Vermis_1_20.7389Rectus_R0.728208Parietal_Inf_L0.723687Cerebellum_4_5_R0.719239Frontal_Inf_Oper_R0.715396Caudate_L0.703270Postcentral_L0.7021069Fusiform_R0.688759Pallidum_L0.68576Vermis_60.67587Amygdala_L0.67263Putamen_L0.660280Frontal_Inf_Orb_R0.660498Frontal_Mid_Orb_R0.657271Vermis_80.65260Insula_R0.646497Rolandic_Oper_L0.645301Cingulum_Mid_L0.635579Olfactory_L0.63280Thalamus_R0.631296Frontal_Inf_Oper_L0.628309Parahippocampal_L0.624298Pallidum_R0.62167Cingulum_Ant_R0.621385Temporal_Mid_R0.6191311Cingulum_Ant_L0.598425Thalamus_L0.596280Cingulum_Mid_R0.574612Vermis_30.54962Rolandic_Oper_R0.548404Heschl_L0.54872Olfactory_R0.53888Cingulum_post_R0.52269Caudate_R0.511287Parahippocampal_R0.506318Hippocampus_L0.493279Amygdala_R0.47473Insula_L0.465545Putamen_R0.436309Hippocampus_R0.411282Among the 116 brain regions, a total of 51 brain regions showed higher contributions to the classification than the average weight value (contribution), including the bilateral cerebellum, left orbitofrontal cortex, left cuneus, left temporal pole, right inferior occipital gyrus, bilateral parietal lobe, and left supplementary motor area (Table2).


## 3.1. Sample Characteristics

There were no significant differences in mean age (t=−0.442, p=0.66), sex (χ2=0.343, p=0.558), and educational level (t=0.33, p=0.743) between the healthy controls and patients with mTBI. Compared with healthy controls, patients with mTBI had higher HAMA score (t=5.077, p<0.001), ADL score (t=4.654, p<0.001), and BDI score (t=3.808, p=0.001), and a lower MMSE score (t=−2.284, p=0.03). The mean time between injury and MRI examination of patients with mTBI was 3.58±3.28 days. The mean GCS score, DRS score, MAS score, and ABS score in patients with mTBI were 14.42±0.88, 2.58±2.36, 44.38±5.86, and 14.42±0.78, respectively. The details are shown in Table 1.Table 1
Demographic and clinical features of patients with acute mTBI and healthy controls.
mTBIHealthy controlst valuep valueAge, years38.88±13.3340.46±11.40-0.4420.66Sex (male, female)24 (15, 9)24 (13, 11)0.3430.558Education, years8.88±3.588.54±3.410.3300.743Postinjury, days3.58±3.28N/AN/AN/AGCS14.42±0.88N/AN/AN/ADRS2.58±2.36N/AN/AN/AMAS44.38±5.86N/AN/AN/AABS14.42±0.78N/AN/AN/AHAMA3.83±3.610.08±0.285.077<0.001MMSE29.04±1.6329.83±0.20-2.2840.03ADL21.71±8.0714.04±6.064.654<0.001BDI1.58±1.910.08±0.283.8080.001


## 3.2. Classification Performance

First, we compared the classification performances of the seven MRI parameters and found they could not differentiate well between healthy controls and patients with mTBI (AUC:0.66±0.03, range, 0.61~0.69; accuracy rate: 66.4%±3.4%, range, 60.2%~70.9%; sensitivity: 64.1%±7.9%, range, 49.0%~75.0%; and specificity: 68.4%±5.6%, range, 61.0%~75.0%). Second, we combined these MRI parameters and found the features with the highest contributions to the classification to discriminate between mTBI and healthy controls. We found that the combination with ALFF, fALFF, DC, VMHC, and short-term FCD significantly reached up the classification accuracy, sensitivity, and specificity and received the highest classification performances among all combination with classification accuracy of 81.1% (p<0.001), sensitivity of 88.0% (p<0.001), and specificity of 75.0% (p<0.001) (Figure 1).Figure 1
Schematic diagram overview of machine learning classification framework. Note: this figure shows the classification of the combination of ALFF, fALFF, DC, VMHC, and short-term FCD in distinguishing the mTBI from the normal controls. The classification received the highest (a) AUC value, (b) classification accuracy, (c) sensitivity, and (d) specificity among all combination.
(a)(b)(c)(d)


## 3.3. Consensus Features and Region Weight

In this study, all consensus features were mapped to AAL116 template (116 brain regions), and each of the 116 brain regions was given a weight value which indicates the contribution to classification model. For the combination with ALFF, fALFF, DC, VMHC, and short-term FCD, Table2 shows the weight ranking of the 116 brain regions from highest to lowest.Table 2
Weight ranking of the 116 brain regions to the classification of the combination with ALFF, fALFF, DC, VMHC, and short-term FCD.
ROI weightVoxel sizeVermis_101.94834Cerebellum_9_R1.470156Cerebellum_10_L1.46540Cerebellum_9_L1.455158Frontal_Mid_Orb_L1.417224Frontal_Sup_Orb_L1.398280Cuneus_L1.353472Cerebellum_Crus2_R1.329539Cerebellum_7b_L1.26598Cerebellum_Crus2_L1.257543Frontal_Mid_Orb_L1.256273Temporal_pole_Mid_L1.222177Occipital_Inf_R1.220316Parietal_Sup_L1.136575Paracentral_lobule_R1.134221Frontal_Sup_R1.1181120Cuneus_R1.117416Cerebellum_Crus1_R1.103723Occipital_Inf_L1.103263Vermis_71.10254Calcarine_L1.076649Occipital_Sup_R1.074407Rectus_L1.070258Postcentral_R1.0551050Paracentral_lobule_L1.053340Precentral_R1.009941Parietal_Inf_R1.008397Occipital_Sup_L0.998373Cerebellum_10_R0.99737Cerebellum_7b_R0.97578Cerebellum_8_L0.961303Cerebellum_6_L0.960524Vermis_90.95950Temporal_Inf_R0.9591076Occipital_Mid_L0.954947Cerebellum_Crus1_L0.950725Lingual_L0.945662Supp_motor_area_L0.937630Frontal_Mid_R0.9311448Calcarine_R0.925528Temporal_Mid_L0.9221437Parietal_Sup_R0.921569Cerebellum_4_5_L0.919352Frontal_Sup_medial_L0.918847Lingual_R0.916683Angular_R0.916511Temporal_pole_Sup_R0.918325Cerebellum_6_R0.908532Precuneus_R0.908927Temporal_Sup_R0.907942Frontal_Sup_L0.890987Angular_L0.860341Precuneus_L0.8511008Cingulum_post_L0.850111Frontal_Inf_Tri_L0.848675Frontal_Mid_L0.8471323Temporal_pole_Sup_L0.828329Temporal_Sup_L0.825694Temporal_pole_Mid_R0.822264Cerebellum_8_R0.810298Cerebellum_3_R0.80165Occipital_Mid_R0.796578Supp_motor_area_R0.790695Vermis_4_50.788176Frontal_Sup_medial_R0.787589Frontal_Inf_Tri_R0.783560Supramarginal_R0.779562Precentral_L0.764931Heschl_R0.76360Frontal_Mid_Orb_R0.762296Frontal_Sup_Orb_R0.752296Frontal_Inf_Orb_L0.752504Cerebellum_3_L0.75042Supramarginal_L0.750357Fusiform_L0.747665Temporal_Inf_L0.745948Vermis_1_20.7389Rectus_R0.728208Parietal_Inf_L0.723687Cerebellum_4_5_R0.719239Frontal_Inf_Oper_R0.715396Caudate_L0.703270Postcentral_L0.7021069Fusiform_R0.688759Pallidum_L0.68576Vermis_60.67587Amygdala_L0.67263Putamen_L0.660280Frontal_Inf_Orb_R0.660498Frontal_Mid_Orb_R0.657271Vermis_80.65260Insula_R0.646497Rolandic_Oper_L0.645301Cingulum_Mid_L0.635579Olfactory_L0.63280Thalamus_R0.631296Frontal_Inf_Oper_L0.628309Parahippocampal_L0.624298Pallidum_R0.62167Cingulum_Ant_R0.621385Temporal_Mid_R0.6191311Cingulum_Ant_L0.598425Thalamus_L0.596280Cingulum_Mid_R0.574612Vermis_30.54962Rolandic_Oper_R0.548404Heschl_L0.54872Olfactory_R0.53888Cingulum_post_R0.52269Caudate_R0.511287Parahippocampal_R0.506318Hippocampus_L0.493279Amygdala_R0.47473Insula_L0.465545Putamen_R0.436309Hippocampus_R0.411282Among the 116 brain regions, a total of 51 brain regions showed higher contributions to the classification than the average weight value (contribution), including the bilateral cerebellum, left orbitofrontal cortex, left cuneus, left temporal pole, right inferior occipital gyrus, bilateral parietal lobe, and left supplementary motor area (Table2).


## 4. Discussion

In this case-control study, we documented two novel findings. First, we developed an SVM classifier that was a useful neuroimaging biomarker for mTBI classification. We found that the combination with ALFF, fALFF, DC, VMHC, and short-term FCD received the highest classification performances among all combination (accuracy=81.1%, sensitivity=88.0%, and specificity=75.0%). Second, the consensus brain regions with the highest contributions to classification were located in the bilateral cerebellum, left orbitofrontal cortex, left cuneus, left temporal pole, right inferior occipital gyrus, bilateral parietal lobe, and left supplementary motor areas (contribution above the average value among 116 brain regions).Our study is the first to apply the SVM classifier to find a promising model for mTBI classification. Although several previous studies have offered insights into brain functional and structural abnormalities of mTBI using traditional group-level statistical differences based on one single imaging method, they could not be translated into predictive or diagnostic neurobiological biomarkers for mTBI. The emergence of radiomics has broadened the scope of routine medical imaging, which carried multimodality medical information to reflect the development and progression of diseases [19, 20]. Machine learning classification based on the radiomics strategy allows detecting subtle, nonstrictly localized effects that may remain invisible to the conventional analysis with univariate statistics [21, 22], which are being increasingly used in functional MRI data [15, 16]. These findings could explain the high classification performance of the SVM classifier.Cerebellum is associated with emotion, motor, and advanced cognitive function [23]. The cerebellum anterior lobe is associated with sensorimotor function, and the cerebellum posterior lobe is associated with the regulation of coordinating movement, balance and sleep, and emotional changes [24–28]. Brain volume atrophy and reduction of metabolism functional activity can be found in subjects after TBI [29–31]. Peskind et al. found that soldiers with mTBI showed reduction of glucose metabolism in the cerebellar vermis, cerebellar hemisphere, and pons and functional deficits in attention, language, and working memory [31]. In addition, cerebellar activation was also significantly reduced during auditory-related task stimulation [30]. These studies suggest that the cerebellum plays an important role in the neuropathological basis of mTBI, which supports our findings of high contributions of the cerebellum to the SVM classifier.The prefrontal lobe is one of the brain areas that are most vulnerable to the mTBI. Even minor brain damage can easily cause a damage of the frontal lobe. Studies have found that abnormal functional changes in the frontal lobe are one of the neural mechanisms of emotional numbness, attention, planning, high alertness, and psychological avoidance in patients with posttraumatic injury [32–34]. Keightley et al. found that adolescents with mTBI showed weaker working memory and language function and reduced brain activity in supplementary motor areas, dorsolateral prefrontal lobe, and superior parietal lobe than that of healthy adolescents [35]. Pardini et al. and Jantzen et al. found that parietal lobe and orbitofrontal cortex are associated with severity of mTBI and postconcussion symptoms [36, 37]. Our findings support these studies. Therefore, the abnormal functional changes in the frontal-parietal lobe may be associated with the posttraumatic injury severity and symptoms, which contribute to the high contributions to the SVM classifier.Abnormal functional connectivity between temporal pole and parietal lobe and decreased glucose metabolism in these two areas were found in mTBI patients relative to normal controls [31, 38, 39]. The temporal pole is closely related to the functions such as social interaction, face recognition, semantic memory, mental speculation, and emotion and is responsible for the synthesis of complex and finely processed perceptual input of internal emotions [40]. The abnormal function of the temporal pole in mTBI patients will help us understand the biological mechanism of daily life disorders of mTBI.


## 5. Conclusions

In this study, we developed an SVM classifier that can be severed as a promising sensitive neuroimaging biomarker for mTBI classification based on a combination of multiple imaging indicators. Our analysis using the model showed that the bilateral cerebellum, left orbitofrontal cortex, left cuneus, left temporal pole, right inferior occipital gyrus, bilateral parietal lobe, and left supplementary motor areas exhibited the highest contributions to the classification model. These findings may expand our understanding of the neurobiological mechanism of mTBI. However, there are several limitations that should be addressed. First, the sample size of our study was relatively small. A larger number of sample sizes and multiple center studies are necessary to corroborate our findings. Second, the data of subacute mTBI and follow-up were scarce. Third, this study only used SVM to perform the classification, and other classification methods should be introduced to compare their performances. Fourth, location and size of the lesion, disease of severity, and subtype of mild traumatic brain injury were not considered in the classification.



---

*Source: 3015238-2021-11-18.xml*