Exploring memory function in earthquake trauma survivors with resting-state fMRI and machine learning
BMC Psychiatry volume 20, Article number: 43 (2020)
Traumatized earthquake survivors may develop poor memory function. Resting-state functional magnetic resonance imaging (rs-fMRI) and machine learning techniques may one day aid the clinical assessment of individual psychiatric patients. This study aims to use machine learning with Rs-fMRI from the perspectives of neurophysiology and neuroimaging to explore the association between it and the individual memory function of trauma survivors.
Rs-fMRI data was acquired for eighty-nine survivors (male (33%), average age (SD):45.18(6.31) years) of Wenchuan earthquakes in 2008 each of whom was screened by experienced psychiatrists based on the clinician-administered post-traumatic stress disorder (PTSD) scale (CAPS), and their memory function scores were determined by the Wechsler Memory Scale-IV (WMS-IV). We explored which memory function scores were significantly associated with CAPS scores. Using simple multiple kernel learning (MKL), Rs-fMRI was used to predict the memory function scores that were associated with CAPS scores. A support vector machine (SVM) was also used to make classifications in trauma survivors with or without PTSD.
Spatial addition (SA), which is defined by spatial working memory function, was negatively correlated with the total CAPS score (r = − 0.22, P = 0.04). The use of simple MKL allowed quantitative association of SA scores with statistically significant accuracy (correlation = 0.28, P = 0.03; mean squared error = 8.36; P = 0.04). The left middle frontal gyrus and the left precuneus contributed the largest proportion to the simple MKL association frame. The SVM could not make a quantitative classification of diagnosis with statistically significant accuracy.
The use of the cross-sectional study design after exposure to an earthquake and the leave-one-out cross-validation (LOOCV) increases the risk of overfitting.
Spontaneous brain activity of the left middle frontal gyrus and the left precuneus acquired by rs-fMRI may be a brain mechanism of visual working memory that is related to PTSD symptoms. Machine learning may be a useful tool in the identification of brain mechanisms of memory impairment in trauma survivors.
Trauma survivors are at high risk of developing mental disorders, including a series of psychiatric symptoms (e.g., re-experiencing , depression [1, 2], and anxiety [1, 3, 4]) as well as cognitive decline (e.g., memory impairment [5, 6]). Nevertheless, these symptoms are difficult to measure, especially cognitive decline, which may lead to loss of well-being in later life .
Trauma can change brain function and structure, which may, in turn, lead to memory impairment. For example, memory impairment in mice may persist for extended periods beyond the presentation of predatory stress . Besides, some studies have suggested that patients with PTSD develop impairments in working memory, delayed memory, instantaneous memory, attention, figural memory and visuospatial ability [9,10,11]. Furthermore, trauma survivors who express higher levels of early PTSD symptoms have impaired immediate figural memory and delayed figural memory, which indicates that they are more likely to develop full-blown PTSD in the future . Previous memory models indicate that dissociation and inadequate memory encoding and processing may play a causal role in the development of PTSD [13,14,15,16,17]. Trauma changes the behavior and underlying neurobiology of the individual regardless of whether he or she develops PTSD . Research has demonstrated that trauma survivors have more substance- and alcohol-abuse problems [19,20,21,22]. Regarding changes in underlying neurobiology, previous studies have indicated trauma-induced changes in brain function, i.e., abnormalities in the function and structure of the brain, which can lead to specific cognitive abnormalities, including memory damage and emotional processing abnormalities. For example, Lui et al. demonstrated that in participants who have survived earthquakes, the activity of the frontolimbic and striatal areas increases, while the connectivity among limbic and striatal networks decreases . Earthquake survivors have been found to have a decreased gray matter volume in the bilateral insula, hippocampus, left caudate and putamen and increased gray matter volume in the bilateral orbitofrontal cortex and the parietal lobes 13–25 days after the earthquake . Besides, Golub et al. found that due to the shrinkage of axonal protrusions, a traumatic experience in mice causes a reduction in hippocampal and central amygdala volumes, which are important areas in memory function . Therefore, dividing trauma survivors into two groups (with and without PTSD) based on diagnostic categories can cause brain, cognition, and behavioral changes in trauma survivors without critical PTSD to go unnoticed.
To provide effective precision therapeutics for trauma survivors, it may be helpful to study the effects of trauma on the brain and on an individual’s behavior to further investigate the neuromechanistic and behavioral indicators regardless of the diagnosis, also known as the research domain criteria (RDoC) . To make a contribution to the RDoC in PTSD, Gong et al. quantitatively predicted individual clinical scores of PTSD in earthquake trauma survivors using the resting-state mean amplitude of spontaneous low-frequency fluctuation (mALFF) combined with multivariate machine learning techniques . The spontaneous low-frequency fluctuation (ALFF) is a neuroimaging method used to measure regional spontaneous fluctuations in the BOLD-fMRI signal intensity in fMRI, and it is the averaged square root of activity in spontaneous low-frequency (0.01–0.08 Hz) fluctuations . The mALFF value reflects the degree of its raw ALFF value divided by the average ALFF value of the whole brain, which can avoid the influence of individual differences on brain activity levels . Their study suggested that the prefrontal, parietal and occipital regions make a significant contribution to the association of 17-item PTSD checklist (PCL-17) scores, which challenges the traditional hypothesis that the frontolimbic network is the most important contributor to PTSD symptoms. Gong’s study was the first to suggest that RDoC can be used at the individual level and to explore the essential neuroimaging mechanism of trauma in the brain. They found that a number of regions outside this frontolimbic network, which is traditionally associated with PTSD, contribute to the association with clinical PTSD scores. However, this study only predicted PTSD symptoms.
Noticeably, memory deficits can be objective behavioral indicators related to the severity of trauma-related disorders. Elzinga et al. have indicated that both psychological and neurobiological data create a model for a trauma-related disorder as a disorder of memory [22, 30]. According to general memory models [3, 16, 30,31,32], inadequate memory encoding and processing can affect visual sensory information in the visual-spatial template caused by trauma-related stimuli and result in flashbacks  and other symptoms . Therefore, it would be more efficient to explore the neuroimaging mechanisms, such as mALFF, of specific abnormal cognitive functions than to explore the relationship between symptoms and neuroimaging mechanisms. In addition, cognitive function is more measurable than symptoms estimated by a physician’s evaluation, and it could be the internal neural cognitive basis of post-traumatic mental symptoms , which could be used as supplementary objective criteria for PTSD clinical diagnosis in the future. Although the neuromechanisms underlying trauma-related memory defects remain unclear, it is necessary to further examine whether neuroimaging could be applied to investigate memory functions in trauma-exposed individuals.
The use of multivariate machine learning techniques  has been a recent way to try to increase the translational applicability of functional neuroimaging. Compared with standard mass-univariate analytical methods correlated with cognitive performance measures, these techniques can provide results at the individual level for cognitive function without clinical cognitive function measures as well as an ideal framework for investigating cognitive impairment, which could include a distributed network of regions. Importantly, we used simple multiple kernel learning (simple MKL)  in the current study, which is one of the MKL methods. MKL is a set of machine learning methods using a predefined set of kernels. Kernels are specified by researchers. In the current study, we defined regions of interest (ROIs) as kernels. The ROIs were defined by the AAL template (Additional file 1: Part-1). Kernel methods in MKL are class of algorithms for pattern analysis which can be operated in a high-dimensional, and learn a combination of kernels as part of the algorithm [37,38,39]. Simple MKL is a method that can use support vector machines (SVMs), which is supervised learning models with associated learning algorithms that analyze data used for classification and regression analysis, to simultaneously learn kernels and the associated predictor in supervised learning settings and is based on mixed-norm regularization .
In the present study, we used the Wechsler Memory Scale-IV (WMS-IV)  to investigate the memory function correlated with CAPS-IV  scores in trauma survivors. Consequently, we evaluated the potential of rs-fMRI in making accurate associations about memory function [41,42,43,44,45,46]. Then, we tried to predict the memory function related to PTSD symptoms at an individual level and to explore the mechanism underlying memory deficits. In addition, to compare the accuracy of rs-fMRI in classifying by PTSD diagnosis with the accuracy of rs-fMRI in associating PTSD symptoms with memory deficits, we applied SVM and simple multiple kernel learning (MKL) (Additional file 1: Part-2). We hypothesized that memory function, especially visual memory and working memory, is related to the CAPS-IV score. Besides, we calculated the quantitative associations of individual visual memory functions in trauma survivors by machine learning using rs-fMRI mALFF. Furthermore, we tried to explore the spontaneous activity of brain regions that could be part of the fundamental neuromechanism underlying memory deficits in trauma survivors.
The study was approved by the Medical Research Ethics Committee of West China Hospital, Sichuan University, and written informed consent was obtained from all participants before the study. We acquired whole-brain resting-state fMRI for 89 individuals with a history of trauma (aged between 21 and 60 years, with an average of 45.18 years, all right-handed). According to the RDoC approach, we included subclinical PTSD to investigate the whole-brain mALFF as well as memory functions regardless of the diagnosis. The current study include (a) Individuals with CAPS-IV scores from 0 to 19 (asymptomatic/few symptoms (n = 41, mean = 5.60, SD = 5.97), (b) Individuals with CAPS-IV scores from 20 to 39 (subthreshold PTSD, n = 8, mean = 27.57, SD = 5.19), (c) Individuals with CAPS-IV scores above 39 (moderate to extreme PTSD, n = 40, mean = 85.33, SD =27.20)  (Additional file 1: Part-3). All participants were recruited seven years after the Wenchuan earthquake hit Sichuan, China; participants were recruited through the Mental Health Center of the Western China Hospital, Chengdu, China (Table 1). In addition, all participants had similar socioeconomic and cultural backgrounds. They were all residents of Qingchuan Village (the epicenter of the earthquake), and they were all present at Qingchuan Village when the earthquake occurred. In addition, all participants were interviewed to confirm the following: no history of psychiatric illness among their first-degree relatives; no history of head injury or loss of consciousness (> 1 h) or neurologic disorders; no present or past Axis-I psychiatric disorders other than PTSD; no history of psychotherapy; no learning or developmental disorders; and no history of drug or alcohol abuse in the six months preceding the scan. All subjects were assessed by DSM-IV structured clinical interview (SCID)  and CAPS-IV  by a consensus between two attending psychiatrists and a trained interviewer.
The WMS-IV  is a battery of tests designed to evaluate immediate and delayed recall, working memory, learning, and recognition of information that is presented in visual or verbal modalities. WMS-IV consists of five subtests: logical memory subtest (LM) to assess the narrative memory, vocabulary paired association (VPA) to evaluate speech memory, design (DE) to assess the sensory visual stimulation, visual reproduction (VR) to measure the memory of nonverbal visual stimuli, and spatial addition (SA) to evaluate visual working memory. The primary subtest yields four WMS-IV index scores (visual working memory, auditory memory, visual memory, and delayed memory) and an overall full-scale memory quotient (Additional file 1: Part-4).
Experiments were performed on a 3.0 T magnetic resonance scanner (Siemens 3.0 T Trio Tim, Germany) with an eight-channel phased-array head coil. Functional images were acquired using a single shot, gradient-recalled echo-planar imaging sequence (repetition time (TR)/echo time (TE) = 2000/30 ms; flip angle = 90°). Five dummy scans were collected prior to the actual MRI scans, and the first 5 volumes of the MRI time series were discarded for magnetization stabilization. The slice thickness was 5 mm with no gap, 64 × 64 matrix size and a field of view of 240 × 240 mm2, resulting in a voxel size of 3.75 × 3.75 × 5 mm3. Each brain volume comprised 30 axial slices, and each functional run contained 205 image volumes. Subjects were instructed to relax, keep their eyes closed, and let their minds wander without falling asleep during the 6.8 min scan .
The data were preprocessed using Statistical Parametric Mapping (SPM8, http://www.fil.ion.ucl.ac.uk/spm). The first 5 volumes were discarded for scanner calibration and participant acclimation to the scanning conditions. All data were corrected for slice timing. Besides, head motion correction of the functional scans was performed for the remaining 200 consecutive volumes. The three translational and three rotational motion parameters were first computed during the realignment step. Then, we generated the mean framewise displacement (FD), which reflected the volume-to-volume changes in the head position . Data from 9 trauma survivors in 98 subjects were discarded when the mean FD exceeded 0.25 mm or when translational or rotational parameters exceeded ±1.5 mm or ± 1.5°. Therefore, 89 subjects were included in the analysis of the current study. Demographic data and clinical symptom scores of 98 subjects were shown in Additional file 1: Part-5. We coregistered high-resolution structural images to the mean functional image and segmented them into white matter, gray matter and cerebrospinal fluid in MNI space using “new-segment and DARTEL” in the data processing assistant for resting-state fMRI (DPARSF) after slice timing and realignment. Nuisance covariates, including the head motion (Friston 24-parameter model) , cerebrospinal fluid signals and white matter signals, were regressed out. Next, we removed the linear trend after spatial normalization. Finally, the images were normalized to the MNI space (voxel size: 3 mm3) with the DARTEL tool, and images were smoothed using a 6-mm full width half maximum (FWHM) isotropic Gaussian filter.
The result data were further temporally bandpass filtered (0.01–0.08 Hz) to reduce the effects of low-frequency drift and high-frequency physiological noise. The time series were transformed into the frequency domain using the Fourier transform (FFT) to obtain the power spectrum after linear-trend removal. To calculate the ALFF, the square-rooted power spectrum was obtained. Then, the ALFF of each voxel was divided by the global mean of ALFF values for standardization to obtain the mALFF, which was calculated using the data processing and analysis for brain imaging (DPABI)  (http://rfmri.org/dpabi) software.
The score and subscale scores of the CAPS_IV and cognitive test scores were analyzed by partial correlation analyses, controlling for the potential influences of age, gender, and educational background, achieving multiple comparison corrections by Bootstrap analysis. P < 0.05 was considered to be a statistically significant difference. All tests were performed using SPSS software, version 17.0 for Windows.
Univariate SPM analysis
A standard, univariate approach was used in SPM8 software to investigate the relationship between the resting-state mALFF and the memory subscale scores that are significantly correlated with CAPS-IV scores.
Simple multiple kernel learning analysis
The relationship between the cognitive test scores and the mALFF at rest were examined using simple MKL  as implemented in PRoNTo v2.1 (http://www.mlnl.cs.ucl.ac.uk/pronto/) running under MATLAB (Mathworks, 2010 release). The simple MKL algorithm is based on the objective value of gradient descent on the support vector machine (SVM). Then, gradient descent wrapping is used to iteratively determine the combination of kernels . In the current study, a linear combination of 116 basic kernels was considered for the final synthetic nuclear space. We also regress out the mean FD in the analysis, a leave-one-out cross-validation was applied across participants to obtain estimates for each participant. The accuracy of simple MKL association was calculated, defined as the Pearson’s correlation coefficient and the mean squared error (MSE) between the actual and predicted values of the cognitive scores (Additional file 1: Part-6).
A randomization test, defined as a permutation test, was used to estimate the distribution of correlation and MSE values under the null hypothesis of no association between mALFF and cognitive ability by randomly pairing the input-target data and the simple MKL rerun 1000 times, which provided an estimated P-value for both the correlation coefficient and the observed MSE.
To facilitate visualization, the current study established a table to sort the contribution of different ROIs to the size of the synthesis of the nuclear space. The expected ROI of different ROIs sorted by leave-one-out cross-validation is calculated by the synthetic kernel space. The smaller the expected sort value, the larger the average weight of the ROI is and the higher the ranking is.
Support vector machine (SVM) classification analysis
The mALFF classification of trauma survivors with (CAPS_IV score ≥ 40) or without PTSD (CAPS_IV score < 20) was examined using simple SVM  as implemented in PRoNTo running under MATLAB. SVM was used to investigate the potential of whole-brain mALFF for discriminating among trauma survivors with or without PTSD. Besides, the multi-kernel combination strategy  was used to effectively combine different feature vectors. The weights of different kernels in the SVM were learned based on the training samples . We also regress out the mean FD in the analysis.
A leave-one-out cross-validation was applied to validate the performance of our proposed approach. Thus, when each sample was designated as a test sample (all other samples were left out for the test), the remaining samples were used to train the classification function. In this manner, it could derive an approximately unbiased estimator of the model. To quantify the performance of compared methods, balanced accuracy, sensitivity and specificity are reported, which are defined based on the association results of LOOCV. Balanced accuracy takes the number of samples in each class into account, and it gives equal weight to the accuracies obtained on the test samples of each class.
In the discrimination maps, each voxel carried a certain weightage, indicating its contribution to the classification function. In this way, a discrimination map could be generated.
The details on the kernels and the information on the tuning parameters are in the Additional file 1: Part 1.
Psychological and behavioral data
The study included 89 subjects. Subject demographic data and clinical symptom scores are shown in Table 1.
Correlation between cognitive dimensions and CAPS_IV scale
Table 2 shows the correlations between cognitive performance and the symptoms assessed by CAPS_IV for the trauma survivors. Most of the cognitive functions evaluated in the WMS-IV did not correlate with the total scores of CAPS_IV. The spatial addition (SA) subtest is the only subtest that displayed a negative correlation with the total CAPS_IV scores. The visual memory index, visual working memory index, immediate memory index, visual reproduction, and SA displayed a negative correlation with the symptoms of re-experiencing as measured by CAPS_IV. Moreover, SA displayed a negative correlation with the symptoms of avoidance as measured by CAPS_IV.
Univariate SPM analysis
There are no regions that showed a significant association with SA scores at P < 0.05 (corrected for multiple comparisons using family-wise error (FWE)) (Additional file 1: Part 7).
Simple MKL analysis
As mentioned above, the SA subtest is the only subtest that displayed a negative correlation with the total CAPS_IV scores. The application of simple MKL allowed quantitative association of SA scores by whole-brain resting-state mALFF in statistically significant accuracy (correlation = 0.28, P-value = 0.026; mean squared error = 8.36, P-value = 0.035) (corrected for multiple comparisons using the permutation test, both P < 0.05 was the significant level) (Fig. 1) in trauma survivors. The association was based on functional alterations across the whole brain, particularly in the left frontal middle gyrus and left precuneus, in addition to parietal lobes and occipital regions (Table 3, Fig. 2, and Additional file 1: Part 8). Table 3 expresses the neuroanatomical regions with a contribution to the association frame above 2% across all regions for the rs-fMRI-based MKL used to accurately predict SA. However, the simple MKL to the whole-brain resting-state mALFF data failed to make a statistically significant accurate quantitative association of SA scores (correlation = − 0.04, P-value = 0.497; mean squared error = 9.25, P-value = 0.712) (corrected for multiple comparisons using a permutation test, both P < 0.05 was the significant level) in PTSD patients (CAPS_IV score ≥ 40).
In discriminating among trauma survivors with moderate to extreme PTSD (n = 40) or with asymptomatic/few symptoms (n = 41), the SVM analysis of the whole-brain resting-state mALFF data failed to classify trauma survivors with a statistically significant accuracy (Total accuracy = 55.56%, balanced accuracy = 55.46%, P value = 0.142; sensitivity =47.50%, P value = 0.42; specificity =63.41%, P value = 0.103; positive predictive value = 55.88%; negative predictive value = 55.32%; kappa = 0.158) (corrected for multiple comparisons using permutation test, both P < 0.05 was the significant level).
The current study indicated that the SA score is the only subtest score in WSM-IV that is associated with PTSD symptoms. The current study explored the quantitative association of visual working memory by whole-brain resting-state mALFF at the individual level with statistically significant accuracy. It examined the pathophysiological mechanisms of the visual working memory in trauma survivors. Importantly, the current work used a model that predicts visual working memory and further studied the neuromechanisms of visual working memory in the distributed network of brain regions of trauma survivors, laying the foundation for RDoC .
This study found that visual working memory was negatively correlated with the severity of PTSD symptoms (see Additional file 1: Part 9). These results are consistent with previous studies . Stein et al. suggested that visual-spatial ability is significantly reduced in patients with PTSD who experienced rape compared with that in healthy controls . Besides, among war asylum seekers, the visual-spatial ability in patients with PTSD is significantly lower than that in war survivors without PTSD . Moreover, a negative relationship has also been found between PTSD symptom severity and working memory performance after controlling for age, gender and educational background . We can speculate that the visual working memory impairment may have a detrimental role in the failure of visual information processing [3, 16, 30,31,32] and that the visual working memory impairment is responsible for disrupted declarative memory of visual information [59, 60], leading to PTSD symptoms . Consequently, our data indicated an important role for visual working memory in PTSD.
In the present study, we extended these early results to show that the application of simple MKL in the rs-fMRI can be used for the quantitative association of cognitive scores by whole-brain resting-state mALFF at the individual level in statistically significant accuracy. In recent work on trauma survivors, the combination of multivariate machine learning techniques and rs-fMRI showed the potential for predicting and investigating neuromechanisms at the individual level . The failure in predicting SA scores at the individual level in PTSD patients may be due to the small sample size. Moreover, as the mean score of SA was much lower in PTSD patients, a floor effect may have also resulted in the failure to predict these scores, but the influence might not be large (see Additional file 1: Part 10). This further indicates the importance of the current study, which explored a way to measure memory by neuroimaging, which can be fast, accurate and appropriate for most people. Interestingly, the low accuracy classification of PTSD in trauma survivors is consistent with previous studies. Qiongmin Zhang et al.  classified trauma survivors with or without PTSD with low accuracy (the best accuracy was 64.86%) when using ALFF. Although they applied a combined multimodal feature approach (combining Reho, GMV, and ALFF), the accuracy improved to only 67.57% with a sensitivity of 52.94% , but this accuracy is still insufficient. In contrast, in the identification of PTSD versus healthy controls (HCs), 89.19% of individuals were correctly assigned to the appropriate diagnostic category, which is much better than the identification of trauma survivors with or without PTSD . In addition, when Gong et al.  compared survivors with or without PTSD, they found that gray matter allowed discrimination with an accuracy of 67% (p < 0.001); however, the two groups could not be distinguished based on white matter. This accuracy is much lower than that obtained when comparing PTSD and HCs 76% (p < 0.001) and 85% (p < 0.001) based on gray and white matter, respectively) . Moreover, the application of MKL to the whole-brain resting-state mALFF data did not allow a quantitative association of CAPS_IV scores (see Additional file 1: Part 10) in the current study. This evidence indicates that it may be better to predict or classify based on RDoC than on the diagnosis, and it is more efficient to explore the whole-brain resting-state mALFF of specific abnormal cognitive functions than to explore the relationship between symptoms and whole-brain resting-state mALFF.
The current study found that the left medial frontal gyrus and bilateral precuneus, but mainly the left precuneus (the left precuneus contributed 20% to the association frame versus the 8% of the right precuneus), both of which are in the default mode network (DMN), contributed to a substantial proportion of the association model of SA score and the whole-brain resting-state mALFF in the current study. The DMN is associated with autobiographical memory, stimulus independence and internally focused thought . Specifically, the precuneus is involved in memory processing and spatial location encoding [64, 65]. Moreover, previous studies have observed that greater activation in the left medial frontal gyrus is involved in the processing of working memory [66, 67]. Besides, our results indicated that the frontal lobe and temporal lobe contribute substantially to the association frame, which is consistent with our previous study . Geuze et al. found that, compared with healthy controls, during the coding phase, patients with PTSD showed underactivation of the frontal cortex and overactivation of the temporal cortex . In the retrieval process, decreased activation of the right frontal cortex, the bilateral middle temporal gyri, and the left posterior hippocampus/parahippocampal gyrus in veterans with PTSD were observed . This difference suggests the possible involvement of the frontal lobe - temporal lobe network in the cognitive deficits seen in PTSD , which is in line with the current findings.
Our previous study found that the gray matter volume (GMV) difference of the bilateral middle occipital gyrus, left superior frontal gyrus and bilateral middle frontal gyrus are the most discriminative regions for distinguishing PTSD from HC . Gong et al. have suggested that PTSD and HC could be discriminated based on gray and white matter in several prefrontal, temporal, parietal and occipital regions . In the current study, the functions of these regions are also important to the association frame, which indicated that these regions are vulnerable brain regions related to the dysfunction of trauma survivors, not only for GMV prediction but also for ALFF estimation. This finding indicated that it may be more sensitive in finding the trauma-related brain regions by predicting the SA score. However, previous studies have been limited to the symptoms, for example, predicting whether the individuals have PTSD or not .
In another study, Gong et al. employed multivariate machine learning techniques for the quantitative association of clinical scores (PCL-17) by whole-brain resting-state mALFF in trauma survivors between 10 and 15 months after the event with statistically significant accuracy (correlation = 0.32, P-value = 0.006; the mean sum of squares = 176.88, P-value = 0.001) . The accuracy in this study was higher compared with the current work, which may be due to the smaller sample size used in the current study compared to Gong’s study (188 trauma survivors) and the current study investigated trauma survivors 7 years after the event, the long time may also plays a role in reduce the accuracy. Nevertheless, the present study aimed to investigate trauma survivors 7 years after the event; all subjects were assessed by SCID , CAPS_IV  and the neuropsychological assessment performed by psychiatrists aimed to predict long-term prognostic indicators of trauma survivors (such as CAPS_IV score and memory performance) by machine learning. Noticeably, whole-brain ALFF predicted SA more accurately than it predicted the scores of the CAPS_IV scale (Additional file 1: Part-11). We can speculate that the relationship between neuroimaging and cognitive function is more significant than that between neuroimaging and symptoms. Additionally, the results indicated that during a long time, the neuroimaging of ALFF was more relevant to the SA function than the PTSD symptoms in trauma survivors. Gong et al. demonstrated that ALFF functional activation in several prefrontal, parietal, and occipital regions is the basis of accurate prediction, which is in line with the current study . In the present study, the left frontal middle gyrus and parietal lobes contributed above 32% to the frame. The results indicated the importance of the frontal-parietal network in predicting visual working memory. Saar-Ashkenazy et al. have reported that the brain activation patterns of visual working memory in survivors with PTSD are different from those of non-PTSD trauma survivors . The frontal cortex is structurally and functionally associated with the parietal lobe , regulating spatial memory and visual-spatial processing [71, 72]. Furthermore, Clark et al. have suggested that the encoding of visual working memory in trauma survivors with PTSD is more dependent on the spatial coding ability of the parietal lobe and less dependent on the executive control function of the frontal lobe . Therefore, the dysfunctional network between the frontal lobe and parietal lobe may be the pathophysiological mechanism behind the disorder of visual working memory and post-traumatic mental disorder, especially PTSD; this current finding is in line with previous findings. The relationship between the prefrontal and parietal lobes and whether the frontal-parietal network is a common mechanism behind PTSD visual working memory impairment and PTSD symptoms needs to be further investigated.
Noticeably, in addition to the brain regions found in the above-reported studies, the current work found that the supplementary motor area (SMA) had a certain proportion of the association of SA, especially the right SMA that had an association weight of over 5% in the association model. MacNamara et al. found that patients with PTSD increased the activation of the motor area in the process of emotional regulation, unlike trauma survivors without PTSD . SMA is believed to have a key role in the network of neural regions mediating top-down control of negative affect , and it may be involved in implementing dorsolateral prefrontal cortex-initiated emotion regulatory effects . Cunnington et al. found that SMA is essential in the early component of premovement activity, which is strongly influenced by higher cognitive factors . In addition, Whalley et al. found that individuals with PTSD exhibited flashback-specific increases in the SMA .
In contrast, the univariate analysis of the resting-state mALFF data did not reveal any regions that were significantly associated with SA scores, which is in line with a previous study  that found that univariate analysis was not significantly associated with clinical scores while multivariate methods were. This may be due to possible sample heterogeneity of mALFF at the group level, for example, the heterogeneity of brain regions and the individual; however, the algorithm is applied to separate individuals, so it is not a problem in machine learning analysis . In addition, the standard univariate approach explored linear correlations; however, some correlations could be found in nonlinear relationships . For example, in the current study, we applied nonlinear transformations by using a kernel function. The results suggest that multivariate methods are more sensitive to the stable diffuse alterations observed in psychiatric disorders. Thus, compared with standard mass-univariate techniques, multivariate methods are more suitable for development as a real-world clinical diagnostic tool .
Importantly, Lianne et al. found that multivariable machine learning techniques, which allow for individual associations based on high-dimensional data, are more sensitive to spatially distributed effects  and changes in brain regions exposed to trauma compared with standard quality - single variable techniques  used in previous studies . Therefore, these multivariable machine learning techniques might be more suitable for clinical application. Moreover, the simple MKL showed the results related to the average contribution of each brain region to the model, rather than just showing these regions by setting the threshold to 30% of the maximum weight value .
This study has several limitations as well. One limitation is the use of cross-sectional study design with participants who were exposed to an earthquake. It is not possible to determine whether the observed cognitive function and the variability in brain function reflect potential pre-existing plastic changes in individual psychological vulnerability or if they occurred after the earthquake. Another limitation is that, because gender is one of the risk factors of PTSD, the potential neuropsychological mechanism of PTSD may differ, and future research needs to separate these mechanisms based on gender. Moreover, both the SA score and the resting-state mALFF may be influenced by the symptoms and age; however, as the covariates that are correlated with the targets could not be regressed out in the model—as this could lead to biases (positive or negative) in the obtained results , we could not exclude the influence of age and PTSD symptoms on the model. We are also aware that the effect size in our study is not large; however, if we exclude individuals with a very low SA score (SA = 1), the effect size of the frame is much better (correlation = 0.36, P-value = 0.010; mean squared error = 5.09, P-value = 0.010; corrected for multiple comparisons using the permutation test, both P < 0.05 was the significant level) (Additional file 1: Part-12), but we do not have the evidence to exclude them. Thus, larger sample size is needed in the future, and the influence of outliers will then be reduced. Besides, the multimodal prediction is a promising field. For example, we will explore whether both rs-fMRI and SA could inform PTSD symptoms in the future. Moreover, although calculating permutation confidence intervals was a computationally difficult problem, a method for performing this calculation in the two-sample problem was presented . However, as far as we know, there is no software package in neuroimaging machine learning that can do this. Besides, we only explored the association frame based on the AAL brain atlas. No consensus atlas has been proven to be far superior to others; AAL remains the most-used atlas in machine learning analysis, although it has many drawbacks. Although AAL is widely accepted in neuroimaging studies, new brain atlases have frequently been used in neuroimaging and machine learning studies, such as the Power 264-region atlas [82, 83] and the Dosenbach’s 160 functional atlas . Future studies should verify our results using these atlases in the analyses of different brain networks. Finally, although the LOOCV was helpful in developing the association model and discovering the critical features derived from rs-fMRI, it increases the risk of overfitting. However, in the current study, LOOCV was adopted to prevent the training set from deviating too much from the overall population for the relatively smaller sample (89 subjects in total) of subjects included in our study. Out-of-sample validation could be the best method to use in this study, but we did not include it in the study design. In the future, we will increase the sample size so that out-of-sample validation could be used in the study to avoid the risk of overfitting.
In conclusion, this study investigated the association between memory and PTSD symptoms, and it indicated that visual working memory impairment was related to PTSD symptoms. We also explored individual associations between mALFF and visual working memory—which may be suitable for development as a real-world clinical memory assistive assessment tool—and the neuromechanism of the spatial overlay dimension of traumatic exposure using the multicore learning method. In addition, our findings indicated that, from the perspective of the whole-brain pattern, brain mechanisms (found in the frontal cortex and parietal cortex) largely contribute to the impairment of visual working memory of trauma survivors, which, to our knowledge, has not been reported before this study. These brain areas are related to memory processing and spatial location encoding. The failure to process visual memory may be due to the dysfunction of these brain areas, which in turn may be related to the symptoms of PTSD [13,14,15,16,17]. This indicates that RDoC studies are more effective in discovering mechanisms than studies that pay more attention to individuals with PTSD diagnosis. Working memory-related training may benefit the functioning of these brain areas  and may be helpful to trauma-exposed people.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Clinician-Administered post-traumatic stress disorder Scale
Statistical Manual of Mental Disorders
Family Wise Error
The Hamilton Anxiety Rating Scale
Hamilton Depression Rating Scale
logical memory subtest
Leave One Out Cross-Validation
The mean amplitude of spontaneous low-frequency fluctuation
Multiple kernel learning
17-item PTSD Checklist
Post-traumatic stress disorder
Research Domain Criteria
Regions of interest
Resting-state functional magnetic resonance imaging
DSM-IV Structured Clinical Interview
Support vector machine
Vocabulary paired association
Wechsler Memory Scale-IV
Farooqui M, Quadri SA, Suriya SS, et al. Posttraumatic stress disorder: a serious post-earthquake complication. Trends Psychiatry Psychother. 2017;39(2):135–43.
Golier JA, Harvey PD, Legge J, Yehuda R. Memory performance in older trauma survivors: implications for the longitudinal course of PTSD. Ann N Y Acad Sci. 2006;1071:54–66.
Breslau N, Kessler RC, Chilcoat HD, Schultz LR, Davis GC, Andreski P. Trauma and posttraumatic stress disorder in the community: the 1996 Detroit area survey of trauma. Arch Gen Psychiatry. 1998;55(7):626–32.
Yehuda RFJD. Differentiating biological correlates of risk, PTSD, and resilience following trauma exposure. J Trauma Stress. 2007;20:435–47.
Buckley TC, Blanchard EB, Neill WT. Information processing and PTSD: a review of the empirical literature. Clin Psychol Rev. 2000;20(8):1041–65.
Tempesta D, Mazza M, Iaria G, De Gennaro L, Ferrara M. A specific deficit in spatial memory acquisition in post-traumatic stress disorder and the role of sleep in its consolidation. Hippocampus. 2012;22(5):1154–63.
Wilson RS, Boyle PA, Segawa E, et al. The influence of cognitive decline on well-being in old age. Psychol Aging. 2013;28(2):304–13.
El HW, Griebel G, Belzung C. Long-term impaired memory following predatory stress in mice. Physiol Behav. 2006;87(1):45–50.
Terada S, Sato S, Nagao S, et al. Trail making test B and brain perfusion imaging in mild cognitive impairment and mild Alzheimer's disease. Psychiatry Res. 2013;213(3):249–55.
MacKinnon MMA, Siegel LS, Perry NE, Weinberg J. Reactivity to stress and the cognitive components of math disability in grade 1 children. J Learn Disabil. 2014;47(4):349–65.
Butollo W, Karl R, König J, Rosner R. A randomized controlled clinical trial of dialogical exposure therapy versus cognitive processing therapy for adult outpatients suffering from PTSD after type I trauma in adulthood. Psychother Psychosom. 2016;85(1):16–26.
Dalia Brandesa GB, Ea AG. PTSD symptoms and cognitive performance in recent trauma survivors. Psychiatry Res. 2002;110:231–8.
Brewin CR. A cognitive neuroscience account of posttraumatic stress disorder and its treatment. Behav Res Ther. 2001. (39): 373–393. [PubMed: 11280338].
Foa EB, Riggs DS. Posttraumatic stress disorder in rape victims. In: Oldham J, Riba MB, Tasman A, editors. , vol. 12. Washington DC: American Psychiatric Press;American Psychiatric Press Review of Psychiatry; 1993. p. 273–303.
Krystal J, Southwick SM, Charney DS. Postrauamtic stress sdisorder: psychobiological mechanisms of traumatic rememberance. Memory distortions: how minds, brains, and societies reconstruct the past. Cambridge: Harvard University Press; 1995.
Terr LC. Childhood trauma: an outline and overview. Am J Psychiatr. 1991;148:10–20.
van der Kolk, Rita Fisler BA. Dissociation and the fragmentary nature of traumatic memories: overview and exploratory study. J Trauma Stress. 1995;8(4):505–25.
Rinne-Albers MA, van der Wee NJ, Lamers-Winkelman F, Vermeiren RR. Neuroimaging in children, adolescents and young adults with psychological trauma. Eur Child Adolesc Psychiatry. 2013;22(12):745–55.
Mills KLTM, Ross JPL. Trauma, PTSD, and substance use disorders: findings from the Australian National Survey of Mental Health and Well-Being. Am J Psychiatry. 2006;163, 4:652–8.
Delker BC. FJJ. From betrayal to the bottle: investigating possible pathways from trauma to problematic substance use. J Trauma Stress. 2014;27(5):576–84.
Mills KL, Teesson M, Ross J, Peters L. Trauma, PTSD, and substance use disorders: findings from the Australian National Survey of mental health and well-being. Am J Psychiatry. 2006;163(4):652–8.
Elzinga BM, Bremner JD. Are the neural substrates of memory the final common pathway in posttraumatic stress disorder (PTSD). J Affect Disord. 2002;70(1):1–17.
Lui S, Huang X, Chen L, et al. High-field MRI reveals an acute impact on brain function in survivors of the magnitude 8.0 earthquake in China. Proc Natl Acad Sci U S A. 2009;106(36):15412–7.
Lui S, Chen L, Yao L, et al. Brain structural plasticity in survivors of a major earthquake. J Psychiatry Neurosci. 2013;38(6):381–7.
Golub Y, Kaltwasser SF, Mauch CP, et al. Reduced hippocampus volume in the mouse model of posttraumatic stress disorder. J Psychiatr Res. 2011;45(5):650–9.
Insle TCB, Garvey M. Ea. research domain criteria (RDoC): toward a new classification framework for research on mental disorders. Am J Psychiatry. 2010;167(7):748–51.
Gong Q, Li L, Du M, et al. Quantitative prediction of individual psychopathology in trauma survivors using resting-state FMRI. Neuropsychopharmacology. 2014;39(3):681–7.
Qi-Hong Zou CZ, Yihong Yang XZ, Xiang-Yu Long QC, Yu-FengWang YZ. An improved approach to detection of amplitude of low-frequency fluctuation (ALFF) for resting-state fMRI: fractional ALFF. J Neurosci Methods. 2008;172:137–41.
Peiyu Huang MX, Quanquan Gu XY, Xiaojun Xu WL, Zhang M. Abnormal amygdala function in Parkinson’s disease patients and its relationship to depression. J Affect Disord. 2015;183:263–8.
van der Kolk BAFR. Dissociation and the fragmentary nature of traumatic memories: overview and exploratory study. J Trauma Stress. 1995;8(4):505–25.
Foa EB, Riggs DS. Post-traumatic stress disorder in rape victims. In: Oldham J, Riba MB, Tasman A, editors. Review of psychiatry, vol. 12. Washington,DC: American Psychiatric Press; 1993. p. 273–303.
Krystal J, Southwick SM, Charney DS. Postrauamtic stress sdisorder: psychobiological mechanisms of traumatic rememberance. In: Schacter DL, editor. Memory distortions: how minds, brains, and societies reconstruct the past. Harvard University press. Cambridge: MA; 1995. p. 150–72.
Brewin CR. Episodic memory, perceptual memory, and their interaction: foundations for a theory of posttraumatic stress disorder. Psychol Bull. 2014;140(1):69–97.
MBaLA Z. Dissociation and memory fragmentation in posttraumatic stress disorder: an evaluation of the dissociative encoding hypothesis. Memory. 2012;20(3):277.
Brammer M. The role of neuroimaging in diagnosis and personalized medicine–current position and likely future directions. Dialogues Clin Neurosci. 2009;11:389–96.
Rakotomamonjy A, Bach F, Canu S, Grandvalet Y. SimpleMKL. J Machine Learn Res. 2008;9:2491–521.
Lanckriet T, Bie D, Cristianini N, Jordan M, Noble W. A statistical framework for genomic data fusion. Bioinformatics. 2004a;20:2626–35.
Lanckriet GRG, Christianini N, Bartlett PL, Ghaoui LE, Jordan MI. Learning the Kernel Matrix with Semi-Definite Programming 2002 : 323–330.
Gert RG, Lanckriet NC, LEG PB, Jordan MI. Learning the kernel matrix with semi-definite programming. J Mach Learn Res. 2004;5:27–72.
Blake D WF, Nagy L ea. Clinician-administered PTSD scale (CAPS). Boston (Mass).. 1990. 7.
Buckley TC, Blanchard EB, Neill W, Trammell NW. Information processing and PTSD: a review of the empirical literature. Clin Psychol Rev. 2000;20:1041–65.
Klonoff H, McDougall G, Clark C, et al. The neuropsychological, psychiatric, and physical effects of prolonged and severe stress: 30 years later. J Nervous Mental Dis. 1976;(63):246–52.
Goldfield AE, Mollica RF, Pesavento BH, et a. The physical and psychological effects of torture: symptomatology and diagnosis 1988. 259: 2725–2729.
Uddo M, Vasterling JJ, Brailey K, Sutker PB. Memory and attention in combat-related posttraumatic stress disorder. PTSD J Psychopathol Behav Assess. 1993;15:43–52.
Wolfe J, Schlesinger LK. Performance of PTSD patients on standard tests of memory. Implications for trauma. In: Yehuda, R., McFarlane, A.C. (Eds.), Psychobiology of Posttraumatic Stress Disorder. New York Academy of Sciences. 1997. 821: 208–218.
Weber MKWD, Rosso IMBJC, WMRea SZJ. Voxel-based morphometric gray matter correlates of posttraumatic stress disorder. J Anxiety Disord. 2013;27:413–9.
Devon E. Hinton /C, Vuth Pich SAS, Stefan G. Hofmann aMHP. A randomized controlled trial of cognitive-behavior therapy for Cambodian refugees with treatment-resistant PTSD and panic attacks: A cross-over design 2005. 18(6): 617–629.
Koo HJ, Han DH, Park SY, Kwon JH. The structured clinical interview for DSM-5 internet gaming disorder: development and validation for diagnosing IGD in adolescents. Psychiatry Investig. 2017;14(1):21–9.
Wechsler D. Wechsler Memory Scale–Fourth Edition. San Antonio, TX: Pearson; 2009.
Zhu HZJ, Zhan WQC, Wu R, MY a e. Altered spontaneous neuronal activity of visual cortex and medial anterior cingulate cortex in treatment-naïve posttraumatic stress disorder. Comprehensive Psychiatry. 2014;55(2014):1688–95.
JD BKA P, AZ ea S. Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion [J]. Neuroimage. 2012;59(3):2142.
Friston K J WS, Howard R ea. Movement-related effects in fMRI time-series. Magn Reson Med. 1996. 35(3): 346.
Yan CG, Wang XD, Zuo XN, Zang YF. DPABI: Data Processing & Analysis for (resting-state) brain imaging. Neuroinformatics. 2016;14(3):339–51.
Alain Rakotomamonjy FRB, Canu S. SimpleMKL. J Mach Learn Res. 2008;9:2491–521.
Liu F, Wee CY, Chen H, Shen D. Inter-modality relationship constrained multi-modalitymulti-task feature selection for Alzheimer’s disease and mild cognitive impairment identification. Neuroimage. 2014;84:466–75. https://doi.org/10.1016/j.neuroimage.2013.09.015.
Ainamani HE, Elbert T, Olema DK, Hecker T. PTSD symptom severity relates to cognitive and psycho-social dysfunctioning - a study with Congolese refugees in Uganda. Eur J Psychotraumatol. 2017;8(1):1283086.
Stein MB, Kennedy CM, Twamley EW. Neuropsychological function in female victims of intimate partner violence with and without posttraumatic stress disorder. Biol Psychiatry. 2002;52(11):1079–88.
Kivling-Bodén G, Sundbom E. Cognitive abilities related to post-traumatic symptoms among refugees from the former Yugoslavia in psychiatric treatment. Nord J Psychiatry. 2003;57(3):191–8.
Roediger HL 3rd. Implicit memory: retention without remembering. Am Psychol. 1990;45(9):1043.
Schacter DL. Searching for memory: the brain, the mind, and the past. BasicBooks. 1996.
Maxfield L, Hyer L. The relationship between efficacy and methodology in studies investigating EMDR treatment of PTSD. J Clin Psychol. 2002;58(1):23–41.
Zhang Q, Wu Q, Zhu H, et al. Multimodal MRI-based classification of trauma survivors with and without post-traumatic stress disorder. Front Neurosci. 2016;10:292.
Gong Q, Li L, Tognin S, et al. Using structural neuroanatomy to identify trauma survivors with and without post-traumatic stress disorder at the individual level. Psychol Med. 2014;44(1):195–203.
Qin P, Northoff G. How is our self related to midline regions and the default-mode network. Neuroimage. 2011;57(3):1221–33.
Ahmed F, Spottiswoode BS, Carey PD, Stein DJ, Seedat S. Relationship between neurocognition and regional brain volumes in traumatized adolescents with and without posttraumatic stress disorder. Neuropsychobiology. 2012;66(3):174–84.
Ragland JD, Turetsky BI, Gur RC, et al. Working memory for complex figures: an fMRI comparison of letter and fractal n-back tasks. Neuropsychology. 2002;16(3):370–9.
Owen AM, Evans AC, Petrides M. Evidence for a two-stage model of spatial working memory processing within the lateral frontal cortex: a positron emission tomography study. Cereb Cortex. 1996;6(1):31–8.
Geuze E, Vermetten E, Ruf M, de Kloet CS, Westenberg HG. Neural correlates of associative learning and memory in veterans with posttraumatic stress disorder. J Psychiatr Res. 2008;42(8):659–69.
Saar-Ashkenazy R, Veksler R, Guez J, et al. Breakdown of inter-hemispheric connectivity is associated with posttraumatic symptomatology and memory impairment. PLoS One. 2016;11(2):e0144766.
Goldman-Rakic PS. Regional and cellular fractionation of working memory. Proc Natl Acad Sci U S A. 1996;93(24):13473–80.
Jonides J, Smith EE, Koeppe RA, Awh E, Minoshima S, Mintun MA. Spatial working memory in humans as revealed by PET. Nature. 1993;363(6430):623–5.
Gentili ERB. Neural correlates of spatial working memory in humans: a functional magnetic resonance imaging study comparing visual and tactile processes. Neuroscience. 2006;139(1):339–49.
Clark CR, Mcfarlane AC, Morris P, et al. Cerebral function in posttraumatic stress disorder during verbal working memory updating: a positron emission tomography study. Biol Psychiatry. 2003;53(6):474–81.
MacNamara A, Rabinak CA, Kennedy AE, et al. Emotion regulatory brain function and SSRI treatment in PTSD: neural correlates and predictors of change. Neuropsychopharmacology. 2016;41(2):611–8.
Kohn NESB, Scheller MLAR, Fox PTHU. Neural network of cognitive emotion regulation- an ALE metaanalysis and MACM analysis. Neuroimage. 2014;87:345–55.
Cunnington RWC, Moser E. Premovement activity of the pre-supplementary motor area and the readiness for action: studies of timeresolved event-related functional MRI. Hum Mov Sci. 2005;(24):644–56.
Whalley MG, Kroes MC, Huntley Z, Rugg MD, Davis SW, Brewin CR. An fMRI investigation of posttraumatic flashbacks. Brain Cogn. 2013;81(1):151–9.
Schmaal L, Marquand AF, Rhebergen D, et al. Predicting the naturalistic course of major depressive disorder using clinical and multimodal neuroimaging information: a multivariate pattern recognition study. Biol Psychiatry. 2015;78(4):278–86.
Mourao-Miranda J, Reinders A, Rocha-Rego V, et al. Individualized prediction of illness courseat the first psychoticepisode:as upport vector machine MRI study 2012. 42: 1037–1047.
Miller GA, Chapman JP. Misunderstanding analysis of covariance. J Abnorm Psychol. 2001;110(1):40–8.
Abd-Elfattah EF. Saddlepoint p-values and confidence intervals for a class of two sample permutation tests for current status and panel count data. Lifetime Data Anal. 2011;17(3):461–72.
Zhang D, Wang J, Liu X, Chen J&L, B Aberrant Brain Network Efficiency in Parkinson’s Disease Patients with Tremor: A Multi-Modality Study Front Aging Neurosci 2015. 7: 169.
Power JD, Cohen AL, Nelson SM, et al. Functional network organization of the human brain. Neuron. 2011;72(4):665–78.
Dosenbach NU, Nardos B, Cohen AL, et al. Prediction of individual brain maturity using fMRI. Science. 2010;329(5997):1358–61.
Stevens MCGA, Bessette KL. Ea. a preliminary study of the effects of working memory training on brain function. Brain Imaging Behav. 2015;10(2):1–21.
We thank the participants and their families for their engagement with this study.
This work was supported by National Key Research & Development Program of China [grant number 2016YFC1307200]; the National Natural Science Foundation of China [grant number 81701328, 81371484, and 81621003]; the Project on Natural Chronic Noninfectious Diseases [grant number 2016YFC1307201]; and the China Postdoctoral Science Foundation [grant number 2017 M612972]; Support Plan of Sichuan [grant number 2014GZ0143].
Ethics approval and consent to participate
All procedures performed in studies involving human participants were in accordance with the ethical standards of the Institutional Review Board of West China Hospital and Sichuan University and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. Informed consent was obtained from all individual participants included in the study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Li, Y., Zhu, H., Ren, Z. et al. Exploring memory function in earthquake trauma survivors with resting-state fMRI and machine learning. BMC Psychiatry 20, 43 (2020). https://doi.org/10.1186/s12888-020-2452-5
- Trauma survivor
- Machine learning