Preliminary Evidence of Sex Differences in Cortical Thickness Following Acute Mild Traumatic Brain Injury.

The main objective of this study was to evaluate sex differences in cortical thickness after acute mild traumatic brain injury (mTBI) and its associations with clinical outcomes. Thirty-two patients with mTBI at acute phase (2.4 ± 1.3 days post-injury) and 30 healthy controls were enrolled. All the participants underwent comprehensive neurocognitive assessments and MRI to assess cortical thickness. Significant sex differences were determined by using variance analysis of factorial design. Relations between the cortical thickness and clinical assessments were measured with the Spearman Correlation. Results revealed that patients with mTBI had significantly reduced cortical thickness in the left entorhinal cortex while increased cortical thickness in the left precuneus cortex and right lateral occipital cortex, compared with healthy controls. The interaction effect of the group × sex on cortical thickness was significant. Female patients had significant thicker cortical thickness in the left caudal anterior cingulate cortex (ACC) than male patients and had higher scores on Posttraumatic stress disorder Checklist-Civilian Version (PCL-C). Spearman correlational analysis showed a significantly positive correlations between the cortical thickness of the left caudal ACC and PCL-C ratings in female patients. Sex differences in cortical thickness support its potential as a neuroimaging phenotype for investigating the differences in clinical profiles of mild TBI between women and men.

Introduction

Traumatic brain injury (TBI) is an important global health issue, of which 75–90% are classified as mild TBI (mTBI) (1). Although most patients with mTBI become asymptomatic within days to weeks, some develop a series of persistent symptoms that have been called as “persistent post-concussive syndrome.” There are several factors associated with those various outcomes, but one of the most controversial and interesting factors is sex (2). Sex differences in outcomes after mTBI have been addressed in lots of studies, of which some found females have a poorer outcome than males (3–8). Another study have found that females have a higher risk of developing epilepsy, suicide, and use more health care services and males have a higher risk for schizophrenia after mTBI (9). While other studies have shown no substantial difference in outcome with regard to sex (10, 11). The effect of sex on outcome after mTBI is less clear. Observation studies may be confounded by many factors, including different symptoms reporting among women and men (12, 7). Therefore, a more objective measurement, such as imaging findings, is crucial in avoiding such bias. Understanding sex differences of brain injury mechanism after mild TBI may change the future diagnostic work-up in patients with mTBI and lead to separate management strategies for patients of different sex.

Sex differences in brain activities after TBI have been reported in recent years. An fMRI study in veterans with TBI reveals that males showed increased functional connectivity between the left orbitofrontal cortex (OFC) and the right mid frontal cortex compared with females. In the meantime, a significant negative association is found on the overall score on the Buss Perry Aggression Questionnaire with functional connectivity between the left OFC and left angular region in male veterans (13). Diffusion tensor imaging (DTI) shows that male patients with mTBI have decreased fractional anisotropy (FA) in the uncinate fasciculus (UF) compared with female patients and was negatively correlated with time to symptom resolution (14). However, there is lack of findings on the cortical thickness alternations following mild TBI. Cortical thickness is an intrinsic biological parameter and should be independent of external factors such as the MRI scanner type, imaging sequence, spatial resolution and/or field strength (15, 16). Sex differences in cortical thickness have been well-documented in healthy adults (17). Cortical thickness alternations are also reported in the chronic phase of various TBI (from vary mild, moderate to more serve TBI) (18–20), there is no study focused on the sex differences in cortical thickness following acute mTBI. Moreover, gender difference in outcomes, favoring females as endogenous neuroprotectants, has been documented in TBI. However, a consistent finding in the research literature on general traumatic experience is that women exhibit twice the rate of the disorder as men, in spite of men experiencing greater lifetime exposure to traumatic events overall (21–23). We hypothesized that there was a significant interaction effect of gender and diagnosis in the traumatic complaints and associated cortical thickness alternations following mTBI.

Materials and methods

Participants

All consecutive patients with non-contrast head CT due to acute head trauma enrolling from the local emergency department (ED) formed the initial population. Inclusion criteria for all mTBI patients were based on the World Health Organization's Collaborating Center for Neurotrauma Task Force (24): (i) Glasgow Coma Scale (GCS) score of 13–15, (ii) one or more/any of the following: loss of consciousness (LOC) for < 30 min, posttraumatic amnesia (PTA) for 24 or less hours, and/or other transient neurological abnormalities such as focal signs, seizure, and intracranial lesions not requiring surgery, (iii) within 1 week after onset of a mTBI, (iv) were aged 16 years or older. Mild TBI patients were excluded for: (1) a history of a previous brain injury, neurological disease, long-standing psychiatric condition, 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 mTBI 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), (6) caused by penetrating craniocerebral injury.

Thirty-two patients with mTBI (18 males) and 30 sex-, age-, and education-matched healthy controls (14 males) without neurologic impairment or psychiatric disorders participated in the study. Participants were all right-handed according to the Edinburgh Handedness Inventory (25). All the subjects gave written, informed consent in person approved by a local institutional review board and conducted in accordance with the Declaration of Helsinki.

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). All the patients with mTBI went through the MRI scans the day they were recruited in the group. The MRI scans were acquired with the use of 3T MRI scanner (GE 750). A custom-built head holder was used to prevent head movements. The MRI protocol involved the high-resolution T1-weighted 3D BRAVO sequence (echo time = 3.4 ms, repetition time = 7.7 ms, flip angle = 9°, slice thickness = 1 mm, field of view = 256 × 256 mm, matrix size = 256 × 256). 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). Any disagreement between these two observers was resolved by consensus. None of patients were with visible contusion lesions using conventional neuroimaging techniques or exhibited cerebral micro-bleeds on SWI.

Clinical assessments

Clinical assessments were performed within 48 h of MR imaging for all the participants. Based on these previous publications, a limited set of neuropsychological tests were analyzed in the current study, to reduce multiple testing issues. This selection was based on our previous work demonstrating sensitivity to TBI-related alterations to brain structure (26, 27). The following tests were assessed: (i) Trail-Making Test Part A and Digit Symbol coding score from the Wechsler Adult Intelligence Scale III (WAIS-III) to examine cognitive information processing speed (28); (ii) Forward Digit Span and Backward Digit Span from the WAIS-III to assess immediate auditory span, working memory, and executive function (29); (iii) Verbal Fluency Test to assess verbal fluency including language ability, semantic memory and executive function (30); (iv) Posttraumatic stress disorder Checklist—Civilian Version (31). In addition, post concussive symptoms (PCS) were measured with the Rivermead Post-Concussion Symptom Questionnaire (RPCS) (32).

Cortical thickness analyses

We used FreeSurfer version 5.3.0 (33) (https://surfer.nmr.mgh.harvard.edu/fswiki) to extract surface-based features from the high-resolution T1-weighted images. The reliability of obtaining cortical thickness measurements from MRI scans has been well validated (15, 34, 35). The high- resolution T1-weighted MR volume for each participant was bias corrected, skull stripped, and segmented into white matter, gray matter, and cerebrospinal fluid before surface-based morphometry (15). Then, we conducted to tessellate the gray-white boundary, perform automated topology correction, and perform surface deformation to locate the gray-white and gray-pial boundaries (36, 37). Cortical thickness was calculated as the closet distance between the gray-white matter boundary and the pial mesh at each vertex on the tessellated surface (37). The surface was then anatomically parcellated by using the Desikan-Killiany atlas into 66 structures (33 structures for each hemisphere) (38). Accuracy for automated processing was inspected by an expert (with 4 years of experience in editing data from more than 200 examinations and trained in this field) and manual corrections were applied if necessary.

We mapped these structures onto a spherical space to achieve point-to-point correspondence for each subject (39). The final segmentation of surface-based labeling was based on both a subject-independent probabilistic atlas and on subject-specific measured values. Combining the cortical-thickness map and surface-based labels, we computed the average cortical thickness for each region.

Statistical analysis

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. P < 0.05 were considered to indicate a significant difference. Effect sizes (Cohen's d) were computed to demonstrate the magnitude of observed differences. Two-sample t-tests were used to explore cortical thickness differences between groups from each native-surface region of interest (ROIs), and results were assessed for significance after controlling the false-discovery rate (FDR) at < 0.05 to correct for multiple comparisons. The 2 × 2 (Group × Sex) mixed measures ANOVAs were performed to test the interactions and group effects, respectively. P < 0.05 were considered to indicate a significant difference. Simple effect was restricted to the ROIs showing significant interaction effect of group and sex. All regional results were Bonferroni-corrected by a factor of number of ROIs (N) showing significant interaction effect, corresponding to a corrected α of P < 0.05/N after controlling for age and education level. Spearman's correlations were conducted between clinical assessments and the region-of-interest variables, since the data were not normally distributed.

Results

Demographic and clinical characteristics

Thirty-seven patients with mTBI participated in this study. Of which, data from five patients were excluded because of poor MR imaging quality (n = 3), and excessive head motions (n = 2). Finally, 32 patients (18 male) were included. Thirty matched healthy controls (14 male) were also recruited. No significant difference was showed between patients with mTBI and healthy controls regarding age and education level. The average age was 33.4 years (range 14–54 years) in healthy controls, and was 31.0 years (range 13–59 years) in patients with mTBI [F(1, 60) = 0.62, P = 0.43, Cohen' s d = −0.20]. The average education level was 11.8 years (range 1–18 years) in healthy controls, and was 9.6 years (range 1–16 years) in patients with mTBI [F(1, 60) = 2.88, P = 0.1, Cohen' s d = −0.43]. No significant difference was in sex (χ2 = 0.57, P = 0.45). No significant difference was found in age and education level among the four groups (i.e., male and female controls and male and female patients with mTBI) [for age: F(3, 58) = 1.98, P = 0.13; for education level: F(3, 58) = 1.32, P = 0.27]. A detailed demographic data and clinical characteristics were summarized in the Tables 1, 2. The major mechanism of trauma was a motor vehicle accident (11 of 18 male patients [61%], 10 of 14 female patients [71%]), followed by an assault (two of 18 male patients [11%], two of 14 female patients [14%]) (P = 0.4).

Table 1

Demographic and Clinical assessments in patients of mTBI and healthy controls.

	Patients of mTBI				Healthy control subjects
Demographic characterstics	Female (n = 14)	Male (n = 18)	P	Cohen's d	Female (n = 16)	Male (n = 14)	P	Cohen's d
Age	33.1 ± 14.3 (13~59)	29.3 ± 10.3 (16~48)	0.41	0.3	37.5 ± 12.3 (20~54)	28.6 ± 8.4 (14~53)	0.19	0.8
Education level	9.5 ± 4.5 (1~16)	9.7 ± 3.2 (5~15)	0.87	0.05	10.7 ± 4.7 (3~18)	12.4 ± 4.8 (1~18)	0.34	0.4
MECHANISM OF INJURY
Motor vehicle accident	10 (71.4%)	11 (61.1%)
Assault	2 (14.2%)	2 (11.1%)
Fall	1 (7%)	2 (11.1%)
Other	1 (7%)	3 (16.7%)
CLINICAL ASSESSMENTS
Trail making test A	55.1 ± 34.9	47.6 ± 22.5	0.712	0.3	49.4 ± 24.1	32.8 ± 22.4	0.093	0.7
RPCS	13.9 ± 8.7	6.9 ± 6.0	0.001	0.9	2.8 ± 2.8	1.1 ± 1.9	0.391	0.7
PCL-C	26.4 ± 8.6	21.0 ± 2.9	0.001	0.8	17.0 ± 0.0	17.0 ±0.0	1.000	Non
DSC	37.4 ± 17.0	35.2 ± 14.3	0.685	0.1	42.8 ± 17.0	54.7 ± 11.8	0.035	−0.8
Forward DS	8.2 ± 1.9	8.0 ± 1.5	0.705	0.1	7.7 ± 1.7	9.1 ± 1.1	0.020	−0.9
Backward DS	3.6 ± 1.7	4.3 ± 1.8	0.247	−0.39	3.8 ± 0.9	5.7 ± 2.0	0.002	−1.3
VF	17.2 ± 5.0	17.1 ± 6.0	0.931	0.01	17.2 ± 5.7	21.3 ± 5.7	0.052	−0.7

RPCS, Rivermead Post-Concussion Symptom Questionnaire; PCL-C, Posttraumatic stress disorder Checklist—Civilian Version; DSC, Digit Symbol coding; DS, Digit Span; VF, Verbal Fluency. Effect sizes reported are Cohen's d-values.

Table 2

Demographic and Clinical assessments in female participants and male participants.

	Female participants				Male participants
Characteristic	MTBI patients	Healthy control	P	Cohen's d	MTBI patients	Healthy control	P	Cohen's d
Age	33.1 ± 14.3 (13~59)	37.5 ± 12.3 (20~54)	0.18	−0.3	29.3 ± 10.3 (16~48)	28.6 ± 8.4 (14~53)	0.93	0.07
Education level	9.5 ± 4.5 (1~16)	10.7 ± 4.7 (3~18)	0.49	−0.26	9.7 ± 3.2 (5~15)	12.4 ± 4.8 (1~18)	0.07	−0.7
CLINICAL ASSESSMENTS
TMT-A	55.1 ± 34.9	49.4 ± 24.1	0.862	0.2	47.6 ± 22.5	32.8 ± 22.4	0.124	0.7
RPCS	13.9 ± 8.7	2.8 ± 2.8	0.000	1.7	6.9 ± 6.0	1.1 ± 1.9	0.004	1.3
PCL-C	26.4 ± 8.6	17.0 ± 0.0	0.000	1.5	21.0 ± 2.9	17.0 ± 0.0	0.012	1.9
DSC	37.4 ± 17.0	42.8 ± 17.0	0.342	−0.3	35.2 ± 14.3	54.7 ± 11.8	0.001	−1.5
Forward DS	8.2 ± 1.9	7.7 ± 1.7	0.366	0.3'	8.0 ± 1.5	9.1 ± 1.1	0.062	−0.8
Backward DS	3.6 ± 1.7	3.8 ± 0.9	0.860	−0.1	4.3 ± 1.8	5.7 ± 2.0	0.023	−0.7
VF	17.3 ± 5.0	17.2 ± 5.7	0.962	0.01	17.1 ± 6.0	21.3 ± 5.7	0.042	−0.7

TMT-A, Trail making test A; RPCS, Rivermead Post-Concussion Symptom Questionnaire; PCL-C, Posttraumatic stress disorder Checklist—Civilian Version; DSC, Digit Symbol coding; DS, Digit Span; VF, Verbal Fluency. Effect sizes reported are Cohen's d-values.

ANCOVAs (Group × Sex) analysis on the clinical assessments were conducted. The interaction effect of the group × sex was significant for the PCL-C [F(1, 58) = 5.99, P = 0.017] with simple effect testing suggested that females presented more complaints in the PCL-C compared with male counterparts only in the patient group but not in the control group (P = 0.001 after Bonferroni-correction) (Table 2). For patients, more complaints on the PCL-C was presented in both female (P < 0.001) and male (P = 0.012) compared with their corresponding controls.

Cortical thickness results

Patients with mTBI presented prominently reduced cortical thickness than healthy controls in the left entorhinal cortex, while significantly increased cortical thickness in the left precuneus cortex and right lateral occipital cortex (P < 0.05, FDR corrected, Figure 1). The interaction effect of the group × sex on cortical thickness was significant in the left caudal anterior cingulate cortex (caudal ACC) [F(1, 57) = 5.976, P = 0.018], fusiform cortex [F(1, 57) = 10.13, P = 0.002], insula cortex [F(1, 57) = 7.35, P = 0.009] and superior frontal cortex (SFC) [F(1, 57) = 5.131, P = 0.027] (Figure 2). Simple effect testing indicated that female patients had significant increased cortical thickness than male patients in the left caudal ACC (P = 0.004). Increased cortical thickness in the left SFC and fusiform cortex (all for P = 0.005) were also presented in female controls compared with male controls. Female patients had non-significant tendency of increased cortical thickness than male patients in the insula cortex (P = 0.036). Other comparison did not obtain the significance.

Figure 1

Regions of significant cortical thinning (blue) and thickening (red) in mild traumatic brain injury patients, compared with healthy controls. Regions in the left entorhinal cortex, left precuneus cortex, and right lateral occipital cortex were significant at P < 0.05 after FDR correction.

Figure 2

Regions of significant increased cortical thickness (blue) in female participants, compared with male participants. Left column represents in patients with mTBI and right column represents in healthy controls. Regions in the caudal anterior cingulate cortex (caudal ACC), fusiform cortex and superior frontal cortex (SFC) in the left hemisphere were significant at P < 0.0125 after Bonferroni-correction. Region in the left insula cortex was not significant at P = 0.036 after Bonferroni-correction.

Correlation analysis results

Spearman correlation analysis was only restricted into clinical assessment (PCL-C) and regional cortical variable showing a significant interaction effect. There was a significantly positive correlations between the cortical thickness of the left caudal ACC and PCL-C ratings only in female patients (r = 0.594, P = 0.011). No other correlation was presented in either female, male patients or whole patient group.

Discussion

To our knowledge, this is the first study to examine sex differences in cortical thickness during the very early acute post-injury of mTBI. The present study indicated a salient modulatory effect of sex on both self-reported symptomatology (PCL-C) and regional cortical thickness following acute mTBI. Female patients had significant increased cortical thickness than male patients in the left caudal ACC. The increased cortical thickness in the left caudal ACC was positively related with more complaints in the PCL-C ratings only in female patients. These findings may provide the clue to the management strategies for mTBI patients of different sexes.

Previous studies about sex differences in outcomes of mTBI revealed confounded results (4–6, 11, 13). A prior research suggested that it may be easier to admit their concussion symptoms after TBI for female athletes than male (12). Similarly, male concussion athletes presented more difficult to have a willing to report their symptoms than female athletes because of social norms (40) and the pressure to return to sports (41). Outcome differences in patients with mTBI could be masked by the subjective assessments. Our study avoids it by evaluating an objective measurement of the underlying neuroimaging detected injury using brain cortical thickness. To date, sex differences in cortical thickness have been reported in healthy individuals in several researches though no such findings were reported in mTBI (41). We found that the female healthy controls had significant increased cortical thickness than male controls in the left fusiform cortex, while this difference disappears modulated by the mild TBI (42). In addition, interaction effect of group and sex in the cortical thickness were primarily located in the four regions, including the left caudal ACC, fusiform cortex, insula cortex, and SFC. Based on these findings, sex difference in cortical thickness may be modulated by the injury.

We observed that the cortical thickness of female patients in the left caudal anterior cingulate cortex was positively related with the PCL-C scores. Previous work has shown that caudal anterior cingulate cortex is involved in motor control (43). For instance, it was reported that post-traumatic stress disorder (PTSD) is an anxiety disorder associated with the anterior cingulate cortex (ACC) (44). Results exhibited that decreased functional connectivity was observed between the caudal ACC and the precentral gyrus in veterans with PTSD compared to healthy controls. So we supposed that the thickened cortical thickness in left caudal anterior cingulate cortex may affect functional connectivity in female survivors. This possibility needs to be tested by future study to focus on dynamic structural and functional changes after mTBI. We have also found non-significant trend of thicker thickness in the left insula in the female patients, compared with male patients. Simple effect testing suggested that female patients had significant increased cortical thickness in this region than female controls that didn't found in the male groups (P = 0.016) after Bonferroni-correction. Such finding may be limited by the relatively small sample sizes. So we still speculated that it may due to the significant increased cortical thickness in female patients than male patients, compared with respective controls.

We did not explore the mechanisms that may contribute to cortical thickening in the current study, but several possibilities in this context can be considered. In animal studies, regional micro-edema has been found in thickened cortical regions within a day after cortical impact (45). Gray matter changes were in the form of increased, not decreased cortical thickness, which may have resulted from neuroinflammatory or other trophic process related to endocrine changes or functional compensation (46). Acute cerebral inflammatory reactions have been found to recover within months after injury and animal studies revealed thickened cortical regions became thinning over days with reduction of micro-edema (45). Unfortunately, we still need to conform such changes in cortical thickness after mTBI in a more chronical follow-up.

There are several limitations to our study. Post-injury time may be not long enough to observe cortical changing at acute phase and longitudinal analysis needs to be involved using following-up data. Furthermore, we did not evaluate the heterogeneity of injury, future studies should use additional outcome measures, including diffusion tensor imaging for structure integrity, resting state functional connectivity study for dynamic changes in functional networks, perfusion of cerebral blood flow (CBF) for brain metabolism using Arterial Spin Labeling Technology (ASL), which may be helpful in understanding the underlying pathophysiology and causes of sex differences in mTBI. Considering the selection of control may influence the detected injury pattern following mTBI, further study needs to enroll both orthopedically-injured patients and healthy subjects as different control group for comparison.

Conclusions

In conclusion, the study presented the abnormal cortical thickness changes related to sex in patients with mTBI, which correlated with the more possibility to develop PTSD and impairments in the information processing speed. Thus, our results indicated a role for cortical thickness as a metric for evaluating the sex difference of mTBI injuries and may predicting subsequent clinical outcome.

Author contributions

MS, LB, ZY, BY, and GB contributed to the conception of the study. MS, JC, and LB contributed significantly to analysis and manuscript preparation. MS, JC, LB, and SW performed the data analyses and wrote the manuscript. WH, SW, CS, SG, LY, DZ, CG, and LH helped perform the analysis with constructive discussions.

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.